def test_calling_cartesian_laplacian_with_precomputed_matrices(
radial_order=4, time_order=2, ut=2e-3, us=np.r_[1e-3, 2e-3, 3e-3]):
ind_mat = qtdmri.qtdmri_index_matrix(radial_order, time_order)
part4_reg_mat_tau = qtdmri.part4_reg_matrix_tau(ind_mat, 1.)
part23_reg_mat_tau = qtdmri.part23_reg_matrix_tau(ind_mat, 1.)
part1_reg_mat_tau = qtdmri.part1_reg_matrix_tau(ind_mat, 1.)
S_mat, T_mat, U_mat = mapmri.mapmri_STU_reg_matrices(radial_order)
laplacian_matrix_precomputed = qtdmri.qtdmri_laplacian_reg_matrix(
ind_mat, us, ut, S_mat, T_mat, U_mat, part1_reg_mat_tau,
part23_reg_mat_tau, part4_reg_mat_tau)
laplacian_matrix_regular = qtdmri.qtdmri_laplacian_reg_matrix(
ind_mat, us, ut)
assert_array_almost_equal(laplacian_matrix_precomputed,
laplacian_matrix_regular)
def test_calling_cartesian_laplacian_with_precomputed_matrices(
radial_order=4, time_order=2, ut=2e-3, us=np.r_[1e-3, 2e-3, 3e-3]):
ind_mat = qtdmri.qtdmri_index_matrix(radial_order, time_order)
part4_reg_mat_tau = qtdmri.part4_reg_matrix_tau(ind_mat, 1.)
part23_reg_mat_tau = qtdmri.part23_reg_matrix_tau(ind_mat, 1.)
part1_reg_mat_tau = qtdmri.part1_reg_matrix_tau(ind_mat, 1.)
S_mat, T_mat, U_mat = mapmri.mapmri_STU_reg_matrices(radial_order)
laplacian_matrix_precomputed = qtdmri.qtdmri_laplacian_reg_matrix(
ind_mat, us, ut, S_mat, T_mat, U_mat,
part1_reg_mat_tau, part23_reg_mat_tau, part4_reg_mat_tau
)
laplacian_matrix_regular = qtdmri.qtdmri_laplacian_reg_matrix(
ind_mat, us, ut)
assert_array_almost_equal(laplacian_matrix_precomputed,
laplacian_matrix_regular)
def test_calling_spherical_laplacian_with_precomputed_matrices(
radial_order=4, time_order=2, ut=2e-3, us=np.r_[2e-3, 2e-3, 2e-3]):
ind_mat = qtdmri.qtdmri_isotropic_index_matrix(radial_order, time_order)
part4_reg_mat_tau = qtdmri.part4_reg_matrix_tau(ind_mat, 1.)
part23_reg_mat_tau = qtdmri.part23_reg_matrix_tau(ind_mat, 1.)
part1_reg_mat_tau = qtdmri.part1_reg_matrix_tau(ind_mat, 1.)
part1_uq_iso_precomp = (
mapmri.mapmri_isotropic_laplacian_reg_matrix_from_index_matrix(
ind_mat[:, :3], 1.))
laplacian_matrix_precomp = qtdmri.qtdmri_isotropic_laplacian_reg_matrix(
ind_mat,
us,
ut,
part1_uq_iso_precomp=part1_uq_iso_precomp,
part1_ut_precomp=part1_reg_mat_tau,
part23_ut_precomp=part23_reg_mat_tau,
part4_ut_precomp=part4_reg_mat_tau)
laplacian_matrix_regular = qtdmri.qtdmri_isotropic_laplacian_reg_matrix(
ind_mat, us, ut)
assert_array_almost_equal(laplacian_matrix_precomp,
laplacian_matrix_regular)
def test_calling_spherical_laplacian_with_precomputed_matrices(
radial_order=4, time_order=2, ut=2e-3, us=np.r_[2e-3, 2e-3, 2e-3]):
ind_mat = qtdmri.qtdmri_isotropic_index_matrix(radial_order, time_order)
part4_reg_mat_tau = qtdmri.part4_reg_matrix_tau(ind_mat, 1.)
part23_reg_mat_tau = qtdmri.part23_reg_matrix_tau(ind_mat, 1.)
part1_reg_mat_tau = qtdmri.part1_reg_matrix_tau(ind_mat, 1.)
part1_uq_iso_precomp = (
mapmri.mapmri_isotropic_laplacian_reg_matrix_from_index_matrix(
ind_mat[:, :3], 1.
)
)
laplacian_matrix_precomp = qtdmri.qtdmri_isotropic_laplacian_reg_matrix(
ind_mat, us, ut,
part1_uq_iso_precomp=part1_uq_iso_precomp,
part1_ut_precomp=part1_reg_mat_tau,
part23_ut_precomp=part23_reg_mat_tau,
part4_ut_precomp=part4_reg_mat_tau)
laplacian_matrix_regular = qtdmri.qtdmri_isotropic_laplacian_reg_matrix(
ind_mat, us, ut)
assert_array_almost_equal(laplacian_matrix_precomp,
laplacian_matrix_regular)
