poisson.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 pdes_mpi_poisson_h
#define pdes_mpi_poisson_h
 
#include "pdes/linear_problem.h"
 
namespace PDEs
{
  namespace MPI
  {
    using namespace dealii;
 
    template <int dim, int spacedim = dim>
    class Poisson : public LinearProblem<dim, spacedim, LAC::LATrilinos>
    {
    public:
      Poisson();
 
      virtual ~Poisson() = default;
 
      virtual void
      custom_estimator(dealii::Vector<float> &error_per_cell) const override;
 
      using ScratchData =
        typename LinearProblem<dim, spacedim, LAC::LATrilinos>::ScratchData;
 
      using CopyData =
        typename LinearProblem<dim, spacedim, LAC::LATrilinos>::CopyData;
 
      using VectorType =
        typename LinearProblem<dim, spacedim, LAC::LATrilinos>::VectorType;
 
    protected:
      virtual void
      assemble_system_one_cell(
        const typename DoFHandler<dim, spacedim>::active_cell_iterator &cell,
        ScratchData                                                    &scratch,
        CopyData &copy) override;
 
      virtual void
      solve() override;
 
      ParsedTools::Function<spacedim> coefficient;
    };
  } // namespace MPI
} // namespace PDEs
#endif