9 #ifndef __laplace_noforce_h 10 #define __laplace_noforce_h 12 #include <amandus/integrator.h> 13 #include <deal.II/base/polynomial.h> 14 #include <deal.II/base/tensor.h> 15 #include <deal.II/integrators/divergence.h> 16 #include <deal.II/integrators/l2.h> 17 #include <deal.II/integrators/laplace.h> 36 virtual void cell(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const;
37 virtual void boundary(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const;
38 virtual void face(DoFInfo<dim>& dinfo1, DoFInfo<dim>& dinfo2, IntegrationInfo<dim>& info1,
39 IntegrationInfo<dim>& info2)
const;
53 virtual void cell(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const;
54 virtual void boundary(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const;
55 virtual void face(DoFInfo<dim>& dinfo1, DoFInfo<dim>& dinfo2, IntegrationInfo<dim>& info1,
56 IntegrationInfo<dim>& info2)
const;
64 this->use_boundary =
false;
65 this->use_face =
false;
83 IntegrationInfo<dim>&)
const 92 this->input_vector_names.push_back(
"Newton iterate");
99 Assert(info.values.size() >= 1, ExcDimensionMismatch(info.values.size(), 1));
100 Assert(info.gradients.size() >= 1, ExcDimensionMismatch(info.values.size(), 1));
102 std::vector<std::vector<double>> rhs(1,
103 std::vector<double>(info.fe_values(0).n_quadrature_points));
111 std::vector<double> null(info.fe_values(0).n_quadrature_points, 0.);
113 const unsigned int deg = info.fe_values(0).get_fe().tensor_degree();
114 Laplace::nitsche_residual(dinfo.vector(0).block(0),
117 info.gradients[0][0],
119 Laplace::compute_penalty(dinfo, dinfo, deg, deg));
125 IntegrationInfo<dim>& info2)
const 127 const unsigned int deg = info1.fe_values(0).get_fe().tensor_degree();
128 Laplace::ip_residual(dinfo1.vector(0).block(0),
129 dinfo2.vector(0).block(0),
133 info1.gradients[0][0],
135 info2.gradients[0][0],
136 Laplace::compute_penalty(dinfo1, dinfo2, deg, deg));
Definition: laplace/noforce.h:48
Definition: laplace/eigen.h:19
Definition: laplace/noforce.h:31
void cell_residual(dealii::Vector< number > &result, 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: elasticity/integrators.h:23
Definition: integrator.h:29