From 137d710d496405013fc8884e34b4fad169f2bec2 Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Wed, 14 Dec 2011 15:55:57 +0100
Subject: [PATCH] moved #pragma omp parall some lines below
---
MathLib/LinAlg/Solvers/CGParallel.cpp | 18 +++++++++---------
1 file changed, 9 insertions(+), 9 deletions(-)
diff --git a/MathLib/LinAlg/Solvers/CGParallel.cpp b/MathLib/LinAlg/Solvers/CGParallel.cpp
index 8c74cb401b4..98014dbad6d 100644
--- a/MathLib/LinAlg/Solvers/CGParallel.cpp
+++ b/MathLib/LinAlg/Solvers/CGParallel.cpp
@@ -71,17 +71,17 @@ unsigned CGParallel(CRSMatrix<double,unsigned> const * mat, double const * const
std::cout << "Step " << l << ", resid=" << resid / nrmb << std::endl;
#endif
+ // r^ = C r
+ // rhat = r
+// blas::copy(N, r, rhat);
+ #pragma omp parallel for
+ for (k = 0; k < N; k++) {
+ rhat[k] = r[k];
+ }
+ mat->precondApply(rhat);
+
#pragma omp parallel
{
- // r^ = C r
- // rhat = r
-// blas::copy(N, r, rhat);
- #pragma omp parallel for
- for (k = 0; k < N; k++) {
- rhat[k] = r[k];
- }
- mat->precondApply(rhat);
-
// rho = r * r^;
rho = scpr(r, rhat, N);
--
GitLab