From 56cf274688c5d6baf85659bdf46ae7759ccf5444 Mon Sep 17 00:00:00 2001
From: Peter Gottschling <peter.gottschling@simunova.com>
Date: Wed, 3 Jun 2009 14:03:14 +0000
Subject: [PATCH] Setzt AMDiS den x-Wert vor dem Loesen???
---
AMDiS/src/ITL_OEMSolver.h | 10 ++++++++++
1 file changed, 10 insertions(+)
diff --git a/AMDiS/src/ITL_OEMSolver.h b/AMDiS/src/ITL_OEMSolver.h
index 3c4c0a1f..13257068 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);
}
--
GitLab