OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed: avoid segfaulting if petsc is requested and support is not available

Browse Source 
git-svn-id: http://svn.sintef.no/trondheim/IFEM/trunk@1392 e10b68d5-8a6e-419e-a041-bce267b0401d
This commit is contained in:
akva 2012-01-13 12:33:54 +00:00 committed by Knut Morten Okstad
parent 7b723f518b
commit 06f3017b56
1 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/LinAlg/SystemMatrix.C
View File

@ -49,9 +49,16 @@ void SystemVector::copy(const SystemVector& x)
SystemMatrix* SystemMatrix::create (Type matrixType, const LinSolParams& spar)
{
if (matrixType == PETSC)
#ifdef HAS_PETSC
if (matrixType == PETSC) return new PETScMatrix(spar);
return new PETScMatrix(spar);
#else
{
std::cerr << "PETSc support not compiled in, bailing" << std::endl;
exit(1);
}
#endif
return SystemMatrix::create(matrixType);
}
@ -64,11 +71,14 @@ SystemMatrix* SystemMatrix::create (Type matrixType, int num_thread_SLU)
case SPR : return new SPRMatrix();
case SPARSE: return new SparseMatrix(SparseMatrix::SUPERLU,num_thread_SLU);
case SAMG : return new SparseMatrix(SparseMatrix::S_A_M_G);
#ifdef HAS_PETSC
case PETSC :
#ifdef HAS_PETSC
// Use default PETSc settings when no parameters are provided by user
static LinSolParams defaultPar;
return new PETScMatrix(defaultPar);
#else
std::cerr << "PETSc support not compiled in, bailing" << std::endl;
exit(1);
#endif
default:
std::cerr <<"SystemMatrix::create: Unsupported matrix type "