exact-solution.cc

Go to the documentation of this file.

 
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/base/function.h>
 
#include <deal.II/lac/block_vector.h>
#include <deal.II/lac/full_matrix.h>
#include <deal.II/lac/block_sparse_matrix.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/precondition.h>
 
#include <deal.II/lac/linear_operator.h>
#include <deal.II/lac/packaged_operation.h>
 
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/dofs/dof_renumbering.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/fe/fe_dgq.h>
#include <deal.II/fe/fe_system.h>
#include <deal.II/fe/fe_values.h>
#include <deal.II/numerics/vector_tools.h>
#include <deal.II/numerics/data_out.h>
 
#include <fstream>
#include <iostream>
 
#include <deal.II/fe/fe_raviart_thomas.h>
 
#include <deal.II/base/tensor_function.h>
 
using namespace dealii;
 
 
template <int dim, int n_equations>
class ExactSolution : public Function<dim>
{
public:
    ExactSolution()
      : Function<dim>(n_equations)
    {}
 
    virtual void vector_value(const Point<dim> &p,
                              Vector<double>   &values) const override;
};
 
 
// Out-of-line implementation of ExactSolution::vector_value() (see the
// class declaration for documentation). Fills @p values with
// (sin(pi x) * sin(pi y), 0).
template <int dim, int n_equations>
void ExactSolution<dim, n_equations>::vector_value(const Point<dim> &p,
                                                   Vector<double>   &values) const
{
    //AssertDimension(values.size(), 2);
 
    double pi = numbers::PI;
    values(0) = std::sin(pi * p[0]) * std::sin(pi * p[1]);
    values(1) = 0;
    //values(2) = 0;
}