7 #ifndef _amandus_matrix_darcy_h 8 #define _amandus_matrix_darcy_h 10 #include <amandus/integrator.h> 11 #include <deal.II/integrators/divergence.h> 12 #include <deal.II/integrators/l2.h> 13 #include <deal.II/meshworker/integration_info.h> 28 virtual void cell(dealii::MeshWorker::DoFInfo<dim>& dinfo,
29 dealii::MeshWorker::IntegrationInfo<dim>& info)
const;
30 virtual void boundary(dealii::MeshWorker::DoFInfo<dim>& dinfo,
31 dealii::MeshWorker::IntegrationInfo<dim>& info)
const;
32 virtual void face(dealii::MeshWorker::DoFInfo<dim>& dinfo1,
33 dealii::MeshWorker::DoFInfo<dim>& dinfo2,
34 dealii::MeshWorker::IntegrationInfo<dim>& info1,
35 dealii::MeshWorker::IntegrationInfo<dim>& info2)
const;
44 this->use_boundary =
false;
45 this->use_face =
false;
51 dealii::MeshWorker::IntegrationInfo<dim>& info)
const 53 AssertDimension(dinfo.n_matrices(), 4);
54 const unsigned int id = dinfo.cell->material_id();
58 L2::mass_matrix(dinfo.matrix(0,
false).matrix, info.fe_values(0), R);
60 dinfo.matrix(1,
false).matrix.copy_transposed(dinfo.matrix(2,
false).matrix);
66 typename dealii::MeshWorker::IntegrationInfo<dim>& )
const 73 dealii::MeshWorker::DoFInfo<dim>& ,
74 dealii::MeshWorker::IntegrationInfo<dim>& ,
75 dealii::MeshWorker::IntegrationInfo<dim>& )
const 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
virtual void face(dealii::MeshWorker::DoFInfo< dim > &dinfo1, dealii::MeshWorker::DoFInfo< dim > &dinfo2, dealii::MeshWorker::IntegrationInfo< dim > &info1, dealii::MeshWorker::IntegrationInfo< dim > &info2) const
Definition: darcy/matrix.h:72
std::vector< double > resistance
Definition: darcy/matrix.h:37
virtual void cell(dealii::MeshWorker::DoFInfo< dim > &dinfo, dealii::MeshWorker::IntegrationInfo< dim > &info) const
Definition: darcy/matrix.h:50
virtual void boundary(dealii::MeshWorker::DoFInfo< dim > &dinfo, dealii::MeshWorker::IntegrationInfo< dim > &info) const
Definition: darcy/matrix.h:65
Definition: estimator.h:33
Definition: darcy/matrix.h:24
Definition: integrator.h:29