From 0026ef4bdd8e40eaa12a3dbc5b441718cfbd6da6 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 3 Jul 2013 13:19:33 +0000
Subject: [PATCH] Whitespace
[[Imported from SVN: r11715]]
---
dune/solvers/solvers/cgsolver.cc | 24 ++++++++++++------------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/dune/solvers/solvers/cgsolver.cc b/dune/solvers/solvers/cgsolver.cc
index b715cf60..35ead789 100644
--- a/dune/solvers/solvers/cgsolver.cc
+++ b/dune/solvers/solvers/cgsolver.cc
@@ -81,25 +81,25 @@ void CGSolver<MatrixType, VectorType>::solve()
// Perform one iteration step
// minimize in given search direction p
- matrix_->mv(p,q); // q=Ap
- alpha = p*q; // scalar product
- lambda = rholast/alpha; // minimization
- x_->axpy(lambda,p); // update solution
- rhs_->axpy(-lambda,q); // update defect
+ matrix_->mv(p,q); // q=Ap
+ alpha = p*q; // scalar product
+ lambda = rholast/alpha; // minimization
+ x_->axpy(lambda,p); // update solution
+ rhs_->axpy(-lambda,q); // update defect
// determine new search direction
- q = 0; // clear correction
+ q = 0; // clear correction
- // apply preconditioner
+ // apply preconditioner
preconditioner_->setProblem(*matrix_,q,b);
preconditioner_->iterate();
q = preconditioner_->getSol();
- rho = q*b; // orthogonalization
- beta = rho/rholast; // scaling factor
- p *= beta; // scale old search direction
- p += q; // orthogonalization with correction
- rholast = rho; // remember rho for recurrence
+ rho = q*b; // orthogonalization
+ beta = rho/rholast; // scaling factor
+ p *= beta; // scale old search direction
+ p += q; // orthogonalization with correction
+ rholast = rho; // remember rho for recurrence
// write iteration to file, if requested
if (this->historyBuffer_!="")
--
GitLab