slight change.

AMDiS/src/UmfPackSolver.h

@@ -90,8 +90,8 @@ namespace AMDiS {
 	
 	// std::cout << "UmfPackSolver: A = \n" << A << "\n";
 	int code= (*solver)(x, b);
-	mtl::dense_vector<value_type> r= b - A * x;
-	std::cout << "UmfPackSolver: ||b-Ax|| = " << two_norm(r) << "\n";
+	mtl::dense_vector<value_type> r(b); r-= A * x; residual= two_norm(r);
+	std::cout << "UmfPackSolver: ||b-Ax|| = " << residual << "\n";
 	return code;
     }