linear_algebra_types.h

Go to the documentation of this file.

/*
  Copyright (C) 2011 - 2024 by the authors of the ASPECT code.
 
  This file is part of ASPECT.
 
  ASPECT 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 2, or (at your option)
  any later version.
 
  ASPECT 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 ASPECT; see the file LICENSE.  If not see
  <http://www.gnu.org/licenses/>.
*/
 
 
#ifndef _aspect_linear_algebra_types_h
#define _aspect_linear_algebra_types_h
 
#include <aspect/global.h>
 
DEAL_II_DISABLE_EXTRA_DIAGNOSTICS
 
#ifdef ASPECT_USE_TPETRA
#include <deal.II/base/memory_space.h>
#include <deal.II/lac/trilinos_tpetra_precondition.h>
#include <deal.II/lac/trilinos_tpetra_sparsity_pattern.h>
#include <deal.II/lac/block_sparsity_pattern.h>
#else
#include <deal.II/lac/generic_linear_algebra.h>
#include <deal.II/lac/trilinos_precondition.h>
#endif
 
DEAL_II_ENABLE_EXTRA_DIAGNOSTICS
 
namespace aspect
{
  namespace LinearAlgebra
  {
#ifdef ASPECT_USE_TPETRA
    using PreconditionBase = ::LinearAlgebra::TpetraWrappers::PreconditionBase<double>;
 
    using PreconditionAMG = ::LinearAlgebra::TpetraWrappers::PreconditionAMGMueLu<double>;
 
    using PreconditionILU = ::LinearAlgebra::TpetraWrappers::PreconditionILU<double>;
 
    using PreconditionJacobi = ::LinearAlgebra::TpetraWrappers::PreconditionJacobi<double>;
#else
    using PreconditionBase = ::TrilinosWrappers::PreconditionBase;
 
    using PreconditionAMG = ::TrilinosWrappers::PreconditionAMG;
 
    using PreconditionILU = ::TrilinosWrappers::PreconditionILU;
 
    using PreconditionJacobi = ::TrilinosWrappers::PreconditionJacobi;
#endif
 
    using BlockDynamicSparsityPattern = ::BlockDynamicSparsityPattern;
 
    using DynamicSparsityPattern = ::DynamicSparsityPattern;
  }
}
 
 
 
#endif