Commit c4e453f9 authored by Thomas Witkowski's avatar Thomas Witkowski
Browse files

Added UMFPACK 5.4

parent 8c3aebe9
...@@ -88,9 +88,11 @@ namespace AMDiS { ...@@ -88,9 +88,11 @@ namespace AMDiS {
else else
solver->update(); solver->update();
// std::cout << "UmfPackSolver: A = \n" << A << "\n";
int code= (*solver)(x, b); int code= (*solver)(x, b);
print(x);
std::cout << std::endl;
mtl::dense_vector<value_type> r(b); r-= A * x; residual= two_norm(r); mtl::dense_vector<value_type> r(b); r-= A * x; residual= two_norm(r);
// std::cout.precision(5);
std::cout << "UmfPackSolver: ||b-Ax|| = " << residual << "\n"; std::cout << "UmfPackSolver: ||b-Ax|| = " << residual << "\n";
return code; return code;
} }
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment