ab81795ef9a23a0e8c49f94bf5b44b4fb153f869,dipy/reconst/tests/test_shm.py,,test_real_sph_harm_fibernav,#,131

Before Change

    cos = np.cos
    assert_array_almost_equal(rsh(0, 0, 0, 0),
                              0.5 / sqrt(pi))
    assert_array_almost_equal(rsh(2, 2, pi / 3, pi / 5),
                              0.25 * sqrt(15. / (2. * pi)) *
                              (sin(pi / 5.)) ** 2. *
                              sin(0 + 2. * pi / 3) * sqrt(2))
    assert_array_almost_equal(rsh(-2, 2, pi / 3, pi / 5),
                              0.25 * sqrt(15. / (2. * pi)) *
                              (sin(pi / 5.)) ** 2. *
                              cos(2. * pi / 3) * sqrt(2))

After Change

    // This test should do for now
    // The mrtrix basis should be the same as re-ordering and re-scaling the
    // fibernav basis
    new_order = [0, 5, 4, 3, 2, 1, 14, 13, 12, 11, 10, 9, 8, 7, 6]
    sphere = hemi_icosahedron.subdivide(2)
    basis, m, n = real_sph_harm_mrtrix(4, sphere.theta, sphere.phi)
    expected = basis[:, new_order]
    expected *= np.where(m == 0, 1., np.sqrt(2))

    fibernav_basis, m, n = real_sph_harm_fibernav(4, sphere.theta, sphere.phi)
    assert_array_almost_equal(fibernav_basis, expected)

Instances