17 #ifndef __deal2__integrators_l2_h
18 #define __deal2__integrators_l2_h
21 #include <deal.II/base/config.h>
23 #include <deal.II/base/quadrature.h>
24 #include <deal.II/lac/full_matrix.h>
25 #include <deal.II/fe/mapping.h>
26 #include <deal.II/fe/fe_values.h>
27 #include <deal.II/meshworker/dof_info.h>
60 const double factor = 1.)
63 const unsigned int n_components = fe.
get_fe().n_components();
67 const double dx = fe.
JxW(k) * factor;
69 for (
unsigned int i=0; i<n_dofs; ++i)
70 for (
unsigned int j=0; j<n_dofs; ++j)
71 for (
unsigned int d=0; d<n_components; ++d)
88 template <
int dim,
typename number>
92 const std::vector<double> &input,
93 const double factor = 1.)
101 for (
unsigned int i=0; i<n_dofs; ++i)
117 template <
int dim,
typename number>
121 const VectorSlice<
const std::vector<std::vector<double> > > &input,
122 const double factor = 1.)
125 const unsigned int fe_components = fe.
get_fe().n_components();
126 const unsigned int n_components = input.size();
132 for (
unsigned int i=0; i<n_dofs; ++i)
133 for (
unsigned int d=0; d<n_components; ++d)
161 const double factor1 = 1.,
162 const double factor2 = 1.)
166 const unsigned int n_components = fe1.
get_fe().n_components();
181 const double dx = fe1.
JxW(k);
183 for (
unsigned int i=0; i<n1_dofs; ++i)
184 for (
unsigned int j=0; j<n1_dofs; ++j)
185 for (
unsigned int d=0; d<n_components; ++d)
192 M11(i,j) += dx * u1*v1;
193 M12(i,j) += dx * u2*v1;
194 M21(i,j) += dx * u1*v2;
195 M22(i,j) += dx * u2*v2;
202 DEAL_II_NAMESPACE_CLOSE
#define AssertDimension(dim1, dim2)
const unsigned int dofs_per_cell
void jump_matrix(FullMatrix< double > &M11, FullMatrix< double > &M12, FullMatrix< double > &M21, FullMatrix< double > &M22, const FEValuesBase< dim > &fe1, const FEValuesBase< dim > &fe2, const double factor1=1., const double factor2=1.)
Library of integrals over cells and faces.
const double & shape_value(const unsigned int function_no, const unsigned int point_no) const
void mass_matrix(FullMatrix< double > &M, const FEValuesBase< dim > &fe, const double factor=1.)
double JxW(const unsigned int quadrature_point) const
#define Assert(cond, exc)
double shape_value_component(const unsigned int function_no, const unsigned int point_no, const unsigned int component) const
const FiniteElement< dim, spacedim > & get_fe() const
::ExceptionBase & ExcNotImplemented()
const unsigned int n_quadrature_points
void L2(Vector< number > &result, const FEValuesBase< dim > &fe, const std::vector< double > &input, const double factor=1.)