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)