linear_algebra_types.h

Go to the documentation of this file.

// ------------------------------------------------------
//
// SPDX-License-Identifier: GPL-2.0-or-later
//
// SPDX-FileCopyrightText: Copyright (C) 2011-2026 by the ASPECT authors.
//
// This file is part of ASPECT.
//
// Detailed license information governing the source code and
// contributions can be found in the folder LICENSES and in
// LICENSE and CONTRIBUTING.md at the top level directory.
//
// ------------------------------------------------------


#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