def _sample_orthonormal_to(mu):
v = np.random.randn(mu.shape[<NUM_LIT:0>])<EOL>proj_mu_v = mu * np.dot(mu, v) / np.linalg.norm(mu)<EOL>orthto = v - proj_mu_v<EOL>return orthto / np.linalg.norm(orthto)<EOL>
Sample point on sphere orthogonal to mu.