| 57 | } |
| 58 | |
| 59 | void test_matrix_norm() |
| 60 | { |
| 61 | la::MatrixCSR A0 = create_operator<double>(MPI_COMM_SELF); |
| 62 | la::MatrixCSR A1 = create_operator<double>(MPI_COMM_WORLD); |
| 63 | CHECK(A1.squared_norm() == Catch::Approx(A0.squared_norm()).epsilon(1e-8)); |
| 64 | } |
| 65 | |
| 66 | void test_matrix_apply() |
| 67 | { |