Setzt AMDiS den x-Wert vor dem Loesen???

AMDiS/src/ITL_OEMSolver.h

@@ -23,6 +23,7 @@
 #ifndef AMDIS_ITL_OEM_SOLVER_H
 #define AMDIS_ITL_OEM_SOLVER_H
 
+#include <iostream>
 #include "OEMSolver.h"
 #include "ITL_Preconditioner.h"
 #include "SolverMatrix.h"
@@ -63,6 +64,15 @@ namespace AMDiS {
     {
       itl::cyclic_iteration<value_type> iter(b, this->max_iter, this->relative, 
 					     this->tolerance, this->print_cycle);
+#if 0
+      std::cout << "A is " << num_rows(A) << " x " << num_cols(A)
+		<< ", one_norm(A) = " << one_norm(A)
+		<< ", b is " << num_rows(b) 
+		<< ", one_norm(b) = " << one_norm(b)
+		<< ", x is " << num_rows(x)
+		<< ", one_norm(x) = " << one_norm(x) << "\n";
+      x= 0.0;
+#endif 
       return ITLSolver()(A, x, b, *this->leftPrecon, *this->rightPrecon, iter);
     }