/* Copyright 2012 SINTEF ICT, Applied Mathematics. This file is part of the Open Porous Media project (OPM). OPM is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. OPM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with OPM. If not, see . */ #if HAVE_CONFIG_H #include "config.h" #endif #include #if HAVE_SUITESPARSE_UMFPACK_H #include #endif #if HAVE_DUNE_ISTL #include #endif #if HAVE_AGMG #include #endif #include #include #include namespace Opm { LinearSolverFactory::LinearSolverFactory() { #if HAVE_SUITESPARSE_UMFPACK_H solver_.reset(new LinearSolverUmfpack); #elif HAVE_DUNE_ISTL solver_.reset(new LinearSolverIstl); #else THROW("No linear solver available, you must have UMFPACK or dune-istl installed to use LinearSolverFactory."); #endif } LinearSolverFactory::LinearSolverFactory(const parameter::ParameterGroup& param) { const std::string ls = param.getDefault("linsolver", "umfpack"); if (ls == "umfpack") { #if HAVE_SUITESPARSE_UMFPACK_H solver_.reset(new LinearSolverUmfpack); #endif } else if (ls == "istl") { #if HAVE_DUNE_ISTL solver_.reset(new LinearSolverIstl(param)); #endif } else if (ls == "agmg") { #if HAVE_AGMG solver_.reset(new LinearSolverAGMG(param)); #endif } else { THROW("Linear solver " << ls << " is unknown."); } if (! solver_) { THROW("Linear solver " << ls << " is not enabled in " "this configuration."); } } LinearSolverFactory::~LinearSolverFactory() { } LinearSolverInterface::LinearSolverReport LinearSolverFactory::solve(const int size, const int nonzeros, const int* ia, const int* ja, const double* sa, const double* rhs, double* solution) const { return solver_->solve(size, nonzeros, ia, ja, sa, rhs, solution); } void LinearSolverFactory::setTolerance(const double tol) { solver_->setTolerance(tol); } double LinearSolverFactory::getTolerance() const { return solver_->getTolerance(); } } // namespace Opm