deallog << "relative norm distance = "
<< x.l2_norm() / ub.l2_norm()
<< std::endl;
+ Assert (x.l2_norm() / ub.l2_norm() < 1e-12,
+ ExcInternalError());
deallog << "Block Sparse Factorization" << std::endl;
umfpackb.factorize(Bb);
deallog << "relative norm distance = "
<< b.l2_norm() / ubb.l2_norm()
<< std::endl;
+ Assert (b.l2_norm() / ubb.l2_norm() < 1e-12,
+ ExcInternalError());
}
DEAL:cg::Starting value 3.12141
DEAL:cg::Convergence step 24 value 0
DEAL:cg::Starting value 3.12141
-DEAL:cg::Convergence step 24 value 0
+DEAL:cg::Convergence step 25 value 0
DEAL::Sparse Factorization
DEAL::absolute norms = 13.7931 13.7931
DEAL::relative norm distance = 0
#include <lac/sparse_direct.h>
-
-template <int dim>
void test ()
{
std::vector<unsigned int> size (2,2U);
ubb(i) = i;
umfpackb.solve(ubb);
+
+ for (unsigned int i=0; i<4; ++i)
+ Assert (std::fabs (ubb(i) - i) < 1e-12,
+ ExcInternalError());
+
+ deallog << "OK" << std::endl;
}
deallog.depth_console(0);
deallog.threshold_double(1.e-10);
- test<1> ();
- test<2> ();
- test<3> ();
+ test ();
}
-DEAL::1d
-DEAL::Number of dofs = 97
-DEAL:cg::Starting value 0.215517
-DEAL:cg::Convergence step 111 value 0
-DEAL:cg::Starting value 0.215517
-DEAL:cg::Convergence step 111 value 0
-DEAL::Sparse Factorization
-DEAL::absolute norms = 6.45748 6.45748
-DEAL::relative norm distance = 0
-DEAL::Block Sparse Factorization
-DEAL::absolute norms = 6.45748 6.45748
-DEAL::relative norm distance = 0
-DEAL::2d
-DEAL::Number of dofs = 994
-DEAL:cg::Starting value 0.595296
-DEAL:cg::Convergence step 58 value 0
-DEAL:cg::Starting value 0.595296
-DEAL:cg::Convergence step 58 value 0
-DEAL::Sparse Factorization
-DEAL::absolute norms = 10.0596 10.0596
-DEAL::relative norm distance = 0
-DEAL::Block Sparse Factorization
-DEAL::absolute norms = 10.0596 10.0596
-DEAL::relative norm distance = 0
-DEAL::3d
-DEAL::Number of dofs = 669
-DEAL:cg::Starting value 3.12141
-DEAL:cg::Convergence step 24 value 0
-DEAL:cg::Starting value 3.12141
-DEAL:cg::Convergence step 24 value 0
-DEAL::Sparse Factorization
-DEAL::absolute norms = 13.7931 13.7931
-DEAL::relative norm distance = 0
-DEAL::Block Sparse Factorization
-DEAL::absolute norms = 13.7931 13.7931
-DEAL::relative norm distance = 0
+DEAL::OK