generic_linear_algebra.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2008 - 2017 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------

#ifndef dealii_generic_linear_algebra_h
#define dealii_generic_linear_algebra_h

#include <deal.II/base/config.h>

#include <deal.II/lac/block_vector.h>
#include <deal.II/lac/precondition.h>
#include <deal.II/lac/sparse_matrix.h>
#include <deal.II/lac/vector.h>


DEAL_II_NAMESPACE_OPEN


namespace LinearAlgebraDealII
{
  using Vector      = Vector<double>;
  using BlockVector = BlockVector<double>;

  using SparseMatrix = SparseMatrix<double>;

  using PreconditionSSOR = PreconditionSSOR<SparseMatrix>;
} // namespace LinearAlgebraDealII


DEAL_II_NAMESPACE_CLOSE


#ifdef DEAL_II_WITH_PETSC

#  include <deal.II/lac/block_sparsity_pattern.h>
#  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>

DEAL_II_NAMESPACE_OPEN

namespace LinearAlgebraPETSc
{
  using namespace dealii;

  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;

  } // namespace MPI

} // namespace LinearAlgebraPETSc
DEAL_II_NAMESPACE_CLOSE


#endif // DEAL_II_WITH_PETSC

#ifdef DEAL_II_WITH_TRILINOS

#  include <deal.II/lac/block_sparsity_pattern.h>
#  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>

DEAL_II_NAMESPACE_OPEN

namespace LinearAlgebraTrilinos
{
  using namespace dealii;

  using SolverCG    = TrilinosWrappers::SolverCG;
  using SolverGMRES = TrilinosWrappers::SolverGMRES;

  namespace MPI
  {
    using Vector = TrilinosWrappers::MPI::Vector;

    using BlockVector = TrilinosWrappers::MPI::BlockVector;

    using SparseMatrix = TrilinosWrappers::SparseMatrix;

    using BlockSparseMatrix = TrilinosWrappers::BlockSparseMatrix;

    using BlockCompressedSparsityPattern =
      TrilinosWrappers::BlockSparsityPattern;

    using PreconditionAMG = TrilinosWrappers::PreconditionAMG;

    using PreconditionIC = TrilinosWrappers::PreconditionIC;

    using PreconditionILU = TrilinosWrappers::PreconditionILU;

    using PreconditionJacobi = TrilinosWrappers::PreconditionJacobi;

    using PreconditionSSOR = TrilinosWrappers::PreconditionSSOR;


  } // namespace MPI

} // namespace LinearAlgebraTrilinos

DEAL_II_NAMESPACE_CLOSE


#endif // DEAL_II_WITH_TRILINOS



#endif