Watch this Author's files
Thank you very much for the very nice and useful OMP code.
I believe there might be a mistake in the orthogonalization loop of "atom_new", when working in fast mode. I think the equation:
atom_new = atom_new - (atom_new'*A_T(:,j))*A_T(:,j);
should be modified to:
atom_new = atom_new - (A_T(:,j)'*atom_new)*A_T(:,j);
This turns out significant when working with complex numbers. Otherwise, one may notice different behavior in slow and fast modes.