7 #ifndef __matrix_laplace_h 8 #define __matrix_laplace_h 10 #include <amandus/integrator.h> 11 #include <deal.II/integrators/divergence.h> 12 #include <deal.II/integrators/l2.h> 13 #include <deal.II/integrators/laplace.h> 14 #include <deal.II/meshworker/integration_info.h> 39 virtual void cell(MeshWorker::DoFInfo<dim>& dinfo, MeshWorker::IntegrationInfo<dim>& info)
const;
40 virtual void boundary(MeshWorker::DoFInfo<dim>& dinfo,
41 MeshWorker::IntegrationInfo<dim>& info)
const;
42 virtual void face(MeshWorker::DoFInfo<dim>& dinfo1, MeshWorker::DoFInfo<dim>& dinfo2,
43 MeshWorker::IntegrationInfo<dim>& info1,
44 MeshWorker::IntegrationInfo<dim>& info2)
const;
49 Matrix<dim>::cell(MeshWorker::DoFInfo<dim>& dinfo, MeshWorker::IntegrationInfo<dim>& info)
const 51 AssertDimension(dinfo.n_matrices(), 1);
58 typename MeshWorker::IntegrationInfo<dim>& info)
const 60 if (info.fe_values(0).get_fe().conforms(FiniteElementData<dim>::H1))
63 const unsigned int deg = info.fe_values(0).get_fe().tensor_degree();
64 Laplace::nitsche_matrix(dinfo.matrix(0,
false).matrix,
66 Laplace::compute_penalty(dinfo, dinfo, deg, deg));
72 MeshWorker::IntegrationInfo<dim>& info1,
73 MeshWorker::IntegrationInfo<dim>& info2)
const 75 if (info1.fe_values(0).get_fe().conforms(FiniteElementData<dim>::H1))
78 const unsigned int deg = info1.fe_values(0).get_fe().tensor_degree();
79 Laplace::ip_matrix(dinfo1.matrix(0,
false).matrix,
80 dinfo1.matrix(0,
true).matrix,
81 dinfo2.matrix(0,
true).matrix,
82 dinfo2.matrix(0,
false).matrix,
85 Laplace::compute_penalty(dinfo1, dinfo2, deg, deg));
Integrator for the matrix of the Laplace operator.
Definition: laplace/matrix.h:36
void cell_matrix(dealii::FullMatrix< double > &M, const dealii::FEValuesBase< dim > &fe, const dealii::VectorSlice< const std::vector< std::vector< dealii::Tensor< 1, dim >>>> &input, double lambda=0., double mu=1.)
Definition: matrix_integrators.h:23
Definition: laplace/eigen.h:19
Definition: integrator.h:29