test_matrix.add(-1, reference_matrix);
+ double frobenius_norm = 0;
for (unsigned int i=0; i<2; ++i)
for (unsigned int j=0; j<2; ++j)
- {
- const double frobenius_norm = test_matrix.block(i,j).frobenius_norm();
+ frobenius_norm += test_matrix.block(i,j).frobenius_norm();
- // the data should add up exactly (unfortunately, there is some roundoff due
- // to the cell similarity detection, but there should not be any similarity
- // for the hypershell geometry)
- AssertThrow(frobenius_norm == 0., ExcInternalError());
- }
+ deallog << "log error in matrix norm: " << std::log(frobenius_norm) << std::endl;
test_rhs.add(-1., reference_rhs);
- AssertThrow(test_rhs.l2_norm() == 0., ExcInternalError());
-
- deallog << "OK" << std::endl;
+ deallog << "log error in vector norm: " << std::log(test_rhs.l2_norm()) << std::endl;
}
-DEAL:2d::OK
-DEAL:2d::OK
-DEAL:2d::OK
-DEAL:3d::OK
-DEAL:3d::OK
-DEAL:3d::OK
+DEAL:2d::log error in matrix norm: -inf
+DEAL:2d::log error in vector norm: -inf
+DEAL:2d::log error in matrix norm: -inf
+DEAL:2d::log error in vector norm: -inf
+DEAL:2d::log error in matrix norm: -inf
+DEAL:2d::log error in vector norm: -inf
+DEAL:3d::log error in matrix norm: -inf
+DEAL:3d::log error in vector norm: -inf
+DEAL:3d::log error in matrix norm: -inf
+DEAL:3d::log error in vector norm: -inf
+DEAL:3d::log error in matrix norm: -inf
+DEAL:3d::log error in vector norm: -inf
// the data should add up exactly (unfortunately, there is some roundoff due
// to the cell similarity detection, but there should not be any similarity
// for the hypershell geometry)
- AssertThrow(frobenius_norm == 0., ExcInternalError());
+ deallog << "log error in matrix norm: " << std::log(frobenius_norm) << std::endl;
test_rhs.add(-1., reference_rhs);
- AssertThrow(test_rhs.l2_norm() == 0., ExcInternalError());
-
- deallog << "OK" << std::endl;
+ deallog << "log error in vector norm: " << std::log(test_rhs.l2_norm()) << std::endl;
}
-DEAL:2d::OK
-DEAL:2d::OK
-DEAL:2d::OK
-DEAL:3d::OK
-DEAL:3d::OK
-DEAL:3d::OK
+DEAL:2d::log error in matrix norm: -inf
+DEAL:2d::log error in vector norm: -inf
+DEAL:2d::log error in matrix norm: -inf
+DEAL:2d::log error in vector norm: -inf
+DEAL:2d::log error in matrix norm: -inf
+DEAL:2d::log error in vector norm: -inf
+DEAL:3d::log error in matrix norm: -inf
+DEAL:3d::log error in vector norm: -inf
+DEAL:3d::log error in matrix norm: -inf
+DEAL:3d::log error in vector norm: -inf
+DEAL:3d::log error in matrix norm: -inf
+DEAL:3d::log error in vector norm: -inf