TrilinosSolver.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2018-2020 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2020 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2018-2020 Total, S.A
 * Copyright (c) 2019-     GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */

#ifndef GEOSX_LINEARALGEBRA_INTERFACES_TRILINOSSOLVER_HPP_
#define GEOSX_LINEARALGEBRA_INTERFACES_TRILINOSSOLVER_HPP_

#include "linearAlgebra/utilities/LinearSolverParameters.hpp"
#include "linearAlgebra/utilities/LinearSolverResult.hpp"

namespace geosx
{

class DofManager;
class EpetraVector;
class EpetraMatrix;

class TrilinosSolver
{
public:

  TrilinosSolver( LinearSolverParameters parameters );

  ~TrilinosSolver();

  void solve( EpetraMatrix & mat,
              EpetraVector & sol,
              EpetraVector & rhs,
              DofManager const * const dofManager = nullptr );

  LinearSolverResult const & result()
  {
    return m_result;
  }

private:

  LinearSolverParameters m_parameters;
  LinearSolverResult m_result;

  void solve_direct( EpetraMatrix & mat,
                     EpetraVector & sol,
                     EpetraVector & rhs );

  void solve_krylov( EpetraMatrix & mat,
                     EpetraVector & sol,
                     EpetraVector & rhs );

};

} // end geosx namespace

#endif /* TRILINOSSOLVER_HPP_ */