lac.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 by Luca Heltai
//
// This file is part of the FSI-suite platform, based on the deal.II library.
//
// The FSI-suite platform 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 3.0 of the License,
// or (at your option) any later version. The full text of the license can be
// found in the file LICENSE at the top level of the FSI-suite platform
// distribution.
//
// ---------------------------------------------------------------------
 
#ifndef fsi_lac_h
#define fsi_lac_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/sparse_direct.h>
#include <deal.II/lac/sparse_matrix.h>
#include <deal.II/lac/sparsity_pattern.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_block_vector.h>
#  include <deal.II/lac/petsc_solver.h>
#  include <deal.II/lac/petsc_sparse_matrix.h>
#  include <deal.II/lac/petsc_vector.h>
 
#endif
 
#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_solver.h>
#  include <deal.II/lac/trilinos_sparse_matrix.h>
#  include <deal.II/lac/trilinos_vector.h>
#endif
 
#include "parsed_lac/amg.h"
#include "parsed_lac/amg_petsc.h"
#include "parsed_lac/ilu.h"
 
namespace LAC
{
  struct LAdealii
  {
    using Vector               = dealii::Vector<double>;
    using BlockVector          = dealii::BlockVector<double>;
    using SparseMatrix         = dealii::SparseMatrix<double>;
    using BlockSparseMatrix    = dealii::BlockSparseMatrix<double>;
    using SparsityPattern      = dealii::SparsityPattern;
    using BlockSparsityPattern = dealii::BlockSparsityPattern;
#ifdef DEAL_II_WITH_TRILINOS
    using AMG = ParsedLAC::AMGPreconditioner;
#else
    // Just use UMFPACK if Trilinos is not available
    using AMG = dealii::SparseDirectUMFPACK;
#endif
    using DirectSolver = dealii::SparseDirectUMFPACK;
  };
 
 
#ifdef DEAL_II_WITH_PETSC
 
  struct LAPETSc
  {
    using Vector               = dealii::PETScWrappers::MPI::Vector;
    using SparseMatrix         = dealii::PETScWrappers::MPI::SparseMatrix;
    using SparsityPattern      = dealii::SparsityPattern;
    using BlockVector          = dealii::PETScWrappers::MPI::BlockVector;
    using BlockSparseMatrix    = dealii::PETScWrappers::MPI::BlockSparseMatrix;
    using BlockSparsityPattern = dealii::BlockSparsityPattern;
    using AMG                  = ParsedLAC::PETScAMGPreconditioner;
    using DirectSolver         = dealii::PETScWrappers::SparseDirectMUMPS;
  };
 
#endif // DEAL_II_WITH_PETSC
 
#ifdef DEAL_II_WITH_TRILINOS
  struct LATrilinos
  {
    using Vector               = dealii::TrilinosWrappers::MPI::Vector;
    using SparseMatrix         = dealii::TrilinosWrappers::SparseMatrix;
    using SparsityPattern      = dealii::TrilinosWrappers::SparsityPattern;
    using BlockVector          = dealii::TrilinosWrappers::MPI::BlockVector;
    using BlockSparseMatrix    = dealii::TrilinosWrappers::BlockSparseMatrix;
    using BlockSparsityPattern = dealii::TrilinosWrappers::BlockSparsityPattern;
    using AMG                  = ParsedLAC::AMGPreconditioner;
    using DirectSolver         = dealii::TrilinosWrappers::SolverDirect;
  };
 
#endif // DEAL_II_WITH_TRILINOS
} // namespace LAC
#endif