step-6.h

Go to the documentation of this file.

 
#ifndef STEP_6_H
#define STEP_6_H
 
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/fe/fe_system.h>
#include <deal.II/grid/tria.h>
#include <deal.II/lac/affine_constraints.h>
#include <deal.II/lac/solver_control.h>
#include <deal.II/lac/sparse_matrix.h>
#include <deal.II/lac/vector.h>
#include <deal.II/base/point.h>
#include <deal.II/base/convergence_table.h>
#include <deal.II/base/parameter_handler.h>
 
using namespace dealii;
 
template <int dim>
class SaddlePointProblem
{
public:
    class Parameters
    {
    public:
        Parameters();
 
        static void declare_parameters(ParameterHandler &prm);
 
        void get_parameters(ParameterHandler &prm);
 
        unsigned int n_refinement_cycles;
 
        unsigned int fe_degree;
 
        double convergence_tolerance;
 
        std::string custom_rhs_0;
 
        std::string custom_rhs_1;
 
        std::string custom_rhs_gamma;
    };
 
    SaddlePointProblem(const Parameters &parameters);
 
    void run();
 
private:
    void setup_system();
 
    void assemble_system();
 
    void assemble_system_bak();
 
    void solve();
 
    void refine_grid();
 
    void output_results(const unsigned int cycle) const;
 
    void output_exact_results(const std::string &filename) const;
 
    double exact_solution(const Point<dim> &p) const;
 
    void compute_error(ConvergenceTable &mytable) const;
 
    void write_convergence_table(ConvergenceTable &mytable) const;
 
    const Parameters &parameters;
 
    Triangulation<dim> triangulation;
 
    FESystem<dim> fe;
 
    DoFHandler<dim> dof_handler;
 
    AffineConstraints<double> constraints;
 
    SparseMatrix<double> system_matrix;
 
    SparsityPattern sparsity_pattern;
 
    Vector<double> solution;
 
    Vector<double> system_rhs;
};
 
 
template <int dim>
double coefficient(const Point<dim> &p)
{
    if (p.square() < 1)
        return 1;
    else
        return 1;
}
 
 
#endif // STEP_6_H