Skip to content
Snippets Groups Projects
Commit 879658df authored by wenqing's avatar wenqing
Browse files

Enabled PETSc 3.5

parent 3e7f467d
No related branches found
No related tags found
No related merge requests found
...@@ -117,7 +117,11 @@ IF(OGS_USE_LIS) ...@@ -117,7 +117,11 @@ IF(OGS_USE_LIS)
ENDIF() ENDIF()
IF(OGS_USE_PETSC) IF(OGS_USE_PETSC)
MESSAGE(STATUS "PETSc version is ${PETSC_VERSION}")
ADD_DEFINITIONS(-DUSE_PETSC) ADD_DEFINITIONS(-DUSE_PETSC)
IF(PETSC_VERSION VERSION_GREATER 3.4)
ADD_DEFINITIONS(-DPETSC_V35)
ENDIF(PETSC_VERSION VERSION_GREATER 3.4)
ENDIF() ENDIF()
IF(OGS_USE_EIGEN) IF(OGS_USE_EIGEN)
......
...@@ -41,7 +41,11 @@ bool PETScLinearSolver::solve(const PETScVector &b, PETScVector &x) ...@@ -41,7 +41,11 @@ bool PETScLinearSolver::solve(const PETScVector &b, PETScVector &x)
PetscMemoryGetCurrentUsage(&mem1); PetscMemoryGetCurrentUsage(&mem1);
#endif #endif
#ifdef PETSC_V35
KSPSetOperators(_solver, _A.getRawMatrix(), _A.getRawMatrix());
#else
KSPSetOperators(_solver, _A.getRawMatrix(), _A.getRawMatrix(), DIFFERENT_NONZERO_PATTERN); KSPSetOperators(_solver, _A.getRawMatrix(), _A.getRawMatrix(), DIFFERENT_NONZERO_PATTERN);
#endif
KSPSolve(_solver, b.getData(), x.getData()); KSPSolve(_solver, b.getData(), x.getData());
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment