9 #ifndef __darcy_noforce_h 10 #define __darcy_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 DarcyNoForceResidual();
37 virtual void cell(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const;
38 virtual void boundary(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const;
39 virtual void face(DoFInfo<dim>& dinfo1, DoFInfo<dim>& dinfo2, IntegrationInfo<dim>& info1,
40 IntegrationInfo<dim>& info2)
const;
46 inline DarcyNoForceResidual<dim>::DarcyNoForceResidual()
48 this->input_vector_names.push_back(
"Newton iterate");
53 DarcyNoForceResidual<dim>::cell(DoFInfo<dim>& dinfo, IntegrationInfo<dim>& info)
const 55 Assert(info.values.size() >= 1, ExcDimensionMismatch(info.values.size(), 1));
56 Assert(info.gradients.size() >= 1, ExcDimensionMismatch(info.values.size(), 1));
58 L2::L2(dinfo.vector(0).block(0), info.fe_values(0), make_slice(info.values[0], 0, dim));
59 Divergence::gradient_residual(
60 dinfo.vector(0).block(0), info.fe_values(0), info.values[0][dim], -1.);
62 dinfo.vector(0).block(1), info.fe_values(1), make_slice(info.gradients[0], 0, dim), 1.);
67 DarcyNoForceResidual<dim>::boundary(DoFInfo<dim>&, IntegrationInfo<dim>&)
const 73 DarcyNoForceResidual<dim>::face(DoFInfo<dim>&, DoFInfo<dim>&, IntegrationInfo<dim>&,
74 IntegrationInfo<dim>&)
const Definition: darcy/noforce.h:33
Definition: estimator.h:33
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