diff --git a/AMDiS/src/ITL_OEMSolver.h b/AMDiS/src/ITL_OEMSolver.h index 3c4c0a1fe200d57476dc86484a50ade47ecf57e9..13257068e73f6c4f8f1cf07134e611376cb2b800 100644 --- a/AMDiS/src/ITL_OEMSolver.h +++ b/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); }