15#ifndef dealii_integrators_grad_div_h
16#define dealii_integrators_grad_div_h
63 const double dx = factor * fe.
JxW(k);
64 for (
unsigned int i = 0; i < n_dofs; ++i)
65 for (
unsigned int j = 0; j < n_dofs; ++j)
72 M(i, j) += dx * divu * divv;
83 template <
int dim,
typename number>
88 const double factor = 1.)
97 const double dx = factor * fetest.
JxW(k);
98 for (
unsigned int i = 0; i < n_dofs; ++i)
103 for (
unsigned int d = 0; d < dim; ++d)
104 du += input[d][k][d];
106 result(i) += dx * du * divv;
134 const double dx = factor * fe.
JxW(k);
136 for (
unsigned int i = 0; i < n_dofs; ++i)
137 for (
unsigned int j = 0; j < n_dofs; ++j)
143 double un = 0., vn = 0.;
144 for (
unsigned int d = 0; d < dim; ++d)
150 M(i, j) += dx * 2. * penalty * un * vn;
151 M(i, j) -= dx * (divu * vn + divv * un);
175 const ArrayView<
const std::vector<double>> &input,
177 const ArrayView<
const std::vector<double>> &data,
189 const double dx = factor * fe.
JxW(k);
194 for (
unsigned int d = 0; d < dim; ++d)
196 umgn += (input[d][k] - data[d][k]) * n[d];
197 divu += Dinput[d][k][d];
200 for (
unsigned int i = 0; i < n_dofs; ++i)
205 for (
unsigned int d = 0; d < dim; ++d)
209 dx * (2. * penalty * umgn * vn - divv * umgn - divu * vn);
229 double factor2 = -1.)
241 const double fi = factor1;
242 const double fe = (factor2 < 0) ? factor1 : factor2;
243 const double f = .5 * (fi + fe);
247 const double dx = fe1.
JxW(k);
249 for (
unsigned int i = 0; i < n_dofs; ++i)
250 for (
unsigned int j = 0; j < n_dofs; ++j)
265 for (
unsigned int d = 0; d < dim; ++d)
273 dx * (-.5 * fi * divvi * uni - .5 * fi * divui * vni +
274 f * penalty * uni * vni);
276 dx * (.5 * fi * divvi * une - .5 * fe * divue * vni -
277 f * penalty * vni * une);
279 dx * (-.5 * fe * divve * uni + .5 * fi * divui * vne -
280 f * penalty * uni * vne);
282 dx * (.5 * fe * divve * une + .5 * fe * divue * vne +
283 f * penalty * une * vne);
305 const ArrayView<
const std::vector<double>> &input1,
307 const ArrayView<
const std::vector<double>> &input2,
310 double int_factor = 1.,
311 double ext_factor = -1.)
321 const double fi = int_factor;
322 const double fe = (ext_factor < 0) ? int_factor : ext_factor;
323 const double penalty = .5 * pen * (fi + fe);
328 const double dx = fe1.
JxW(k);
334 for (
unsigned int d = 0; d < dim; ++d)
336 uni += input1[d][k] * n[d];
337 une += input2[d][k] * n[d];
338 divui += Dinput1[d][k][d];
339 divue += Dinput2[d][k][d];
342 for (
unsigned int i = 0; i < n1; ++i)
350 for (
unsigned int d = 0; d < dim; ++d)
356 result1(i) += dx * (-.5 * fi * divvi * uni -
357 .5 * fi * divui * vni + penalty * uni * vni);
358 result1(i) += dx * (.5 * fi * divvi * une -
359 .5 * fe * divue * vni - penalty * vni * une);
360 result2(i) += dx * (-.5 * fe * divve * uni +
361 .5 * fi * divui * vne - penalty * uni * vne);
362 result2(i) += dx * (.5 * fe * divve * une +
363 .5 * fe * divue * vne + penalty * une * vne);
const unsigned int dofs_per_cell
const Tensor< 1, spacedim > & normal_vector(const unsigned int q_point) const
double shape_value_component(const unsigned int i, const unsigned int q_point, const unsigned int component) const
const unsigned int n_quadrature_points
const FiniteElement< dim, spacedim > & get_fe() const
double JxW(const unsigned int q_point) const
unsigned int n_components() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define AssertVectorVectorDimension(VEC, DIM1, DIM2)
#define AssertDimension(dim1, dim2)
void nitsche_residual(Vector< double > &result, const FEValuesBase< dim > &fe, const ArrayView< const std::vector< double > > &input, const ArrayView< const std::vector< Tensor< 1, dim > > > &Dinput, const ArrayView< const std::vector< double > > &data, double penalty, double factor=1.)
void nitsche_matrix(FullMatrix< double > &M, const FEValuesBase< dim > &fe, double penalty, double factor=1.)
void cell_matrix(FullMatrix< double > &M, const FEValuesBase< dim > &fe, double factor=1.)
void ip_residual(Vector< double > &result1, Vector< double > &result2, const FEValuesBase< dim > &fe1, const FEValuesBase< dim > &fe2, const ArrayView< const std::vector< double > > &input1, const ArrayView< const std::vector< Tensor< 1, dim > > > &Dinput1, const ArrayView< const std::vector< double > > &input2, const ArrayView< const std::vector< Tensor< 1, dim > > > &Dinput2, double pen, double int_factor=1., double ext_factor=-1.)
void ip_matrix(FullMatrix< double > &M11, FullMatrix< double > &M12, FullMatrix< double > &M21, FullMatrix< double > &M22, const FEValuesBase< dim > &fe1, const FEValuesBase< dim > &fe2, double penalty, double factor1=1., double factor2=-1.)
void cell_residual(Vector< number > &result, const FEValuesBase< dim > &fetest, const ArrayView< const std::vector< Tensor< 1, dim > > > &input, const double factor=1.)
Library of integrals over cells and faces.