From 72dd56044573d3ae178c4e834dec4be82b18c130 Mon Sep 17 00:00:00 2001
From: Jonathan Lambrechts <jonathan.lambrechts@uclouvain.be>
Date: Fri, 25 Mar 2011 12:11:41 +0000
Subject: [PATCH] dg : parallel block jacobi matrix free
---
Solver/linearSystemPETSc.hpp | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/Solver/linearSystemPETSc.hpp b/Solver/linearSystemPETSc.hpp
index 73af243fa2..dc66e63041 100644
--- a/Solver/linearSystemPETSc.hpp
+++ b/Solver/linearSystemPETSc.hpp
@@ -103,7 +103,9 @@ void linearSystemPETSc<scalar>::allocate(int nbRows)
_try(MatAppendOptionsPrefix(_a, this->_parameters["petscPrefix"].c_str()));
_try(MatSetFromOptions(_a));
_try(MatGetOwnershipRange(_a, &_localRowStart, &_localRowEnd));
- _try(MatGetSize(_a, &_globalSize, &_localSize));
+ int nbColumns;
+ _localSize = _localRowEnd - _localRowStart;
+ _try(MatGetSize(_a, &_globalSize, &nbColumns));
// preallocation option must be set after other options
_try(VecCreate(PETSC_COMM_WORLD, &_x));
_try(VecSetSizes(_x, nbRows, PETSC_DETERMINE));
--
GitLab