In line 55 in spherecluster/util.py, in the _sample_orthonormal_to function, it reads:
proj_mu_v = mu * np.dot(mu, v) / np.linalg.norm(mu)
Shouldn't it instead be:
proj_mu_v = mu * np.dot(mu, v) / np.linalg.norm(mu)**2
If norm(mu)=1 then it doesn't make a difference, but otherwise they are quite different.
In line 55 in spherecluster/util.py, in the _sample_orthonormal_to function, it reads:
proj_mu_v = mu * np.dot(mu, v) / np.linalg.norm(mu)
Shouldn't it instead be:
proj_mu_v = mu * np.dot(mu, v) / np.linalg.norm(mu)**2
If norm(mu)=1 then it doesn't make a difference, but otherwise they are quite different.