generic_linear_algebra.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 2012 - 2026 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#ifndef dealii_generic_linear_algebra_h
#define dealii_generic_linear_algebra_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/lac/block_sparse_matrix.h>
#include <deal.II/lac/block_sparsity_pattern.h>
#include <deal.II/lac/block_vector.h>
#include <deal.II/lac/precondition.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/solver_gmres.h>
#include <deal.II/lac/sparse_matrix.h>
#include <deal.II/lac/vector.h>
 
#ifdef DEAL_II_WITH_PETSC
#  include <deal.II/lac/petsc_block_sparse_matrix.h>
#  include <deal.II/lac/petsc_precondition.h>
#  include <deal.II/lac/petsc_solver.h>
#  include <deal.II/lac/petsc_sparse_matrix.h>
#endif
 
#ifdef DEAL_II_WITH_TRILINOS
#  include <deal.II/lac/trilinos_block_sparse_matrix.h>
#  include <deal.II/lac/trilinos_precondition.h>
#  include <deal.II/lac/trilinos_solver.h>
#  include <deal.II/lac/trilinos_sparse_matrix.h>
#endif
 
 
DEAL_II_NAMESPACE_OPEN
 
namespace LinearAlgebraDealII
{
  using Vector = Vector<double>;
 
  using BlockVector = BlockVector<double>;
 
  using SparseMatrix = SparseMatrix<double>;
 
  using BlockSparseMatrix = BlockSparseMatrix<double>;
 
  using PreconditionSSOR = PreconditionSSOR<SparseMatrix>;
} // namespace LinearAlgebraDealII
 
 
#ifdef DEAL_II_WITH_PETSC
 
namespace LinearAlgebraPETSc
{
  using SolverCG = PETScWrappers::SolverCG;
 
  using SolverGMRES = PETScWrappers::SolverGMRES;
 
  namespace MPI
  {
    using Vector = PETScWrappers::MPI::Vector;
 
    using BlockVector = PETScWrappers::MPI::BlockVector;
 
    using SparseMatrix = PETScWrappers::MPI::SparseMatrix;
 
    using BlockSparseMatrix = PETScWrappers::MPI::BlockSparseMatrix;
 
    using BlockCompressedSparsityPattern = ::BlockDynamicSparsityPattern;
 
    using PreconditionAMG = PETScWrappers::PreconditionBoomerAMG;
 
    using PreconditionIC = PETScWrappers::PreconditionICC;
 
    using PreconditionILU = PETScWrappers::PreconditionILU;
 
    using PreconditionJacobi = PETScWrappers::PreconditionJacobi;
 
    using PreconditionSSOR = PETScWrappers::PreconditionSSOR;
 
    using PreconditionASM = PETScWrappers::PreconditionASM;
 
  } // namespace MPI
 
} // namespace LinearAlgebraPETSc
#endif // DEAL_II_WITH_PETSC
 
#ifdef DEAL_II_WITH_TRILINOS
 
namespace LinearAlgebraTrilinos
{
#  ifdef DEAL_II_TRILINOS_WITH_EPETRA
  using SolverCG = TrilinosWrappers::SolverCG;
 
  using SolverGMRES = TrilinosWrappers::SolverGMRES;
#  else
#    ifdef DEAL_II_TRILINOS_WITH_TPETRA
  using SolverCG = ::SolverCG<TrilinosWrappers::MPI::Vector>;
 
  using SolverGMRES = ::SolverGMRES<TrilinosWrappers::MPI::Vector>;
#    endif
#  endif
 
  namespace MPI
  {
    using Vector = TrilinosWrappers::MPI::Vector;
 
    using BlockVector = TrilinosWrappers::MPI::BlockVector;
 
    using SparseMatrix = TrilinosWrappers::SparseMatrix;
 
    using BlockSparseMatrix = TrilinosWrappers::BlockSparseMatrix;
 
#  ifdef DEAL_II_TRILINOS_WITH_EPETRA
    using BlockCompressedSparsityPattern =
      TrilinosWrappers::BlockSparsityPattern;
 
    using PreconditionAMG = TrilinosWrappers::PreconditionAMG;
 
    using PreconditionIC = TrilinosWrappers::PreconditionIC;
#  else
#    ifdef DEAL_II_TRILINOS_WITH_TPETRA_MUELU
    using PreconditionAMG =
      LinearAlgebra::TpetraWrappers::PreconditionAMGMueLu<double>;
#    endif
#  endif
 
    using PreconditionILU = TrilinosWrappers::PreconditionILU;
 
    using PreconditionJacobi = TrilinosWrappers::PreconditionJacobi;
 
    using PreconditionSSOR = TrilinosWrappers::PreconditionSSOR;
 
 
  } // namespace MPI
 
} // namespace LinearAlgebraTrilinos
 
 
#endif // DEAL_II_WITH_TRILINOS
 
DEAL_II_NAMESPACE_CLOSE
 
#endif