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)