From 43f4fc1bbf342df74cefa161477afaf32b2cab75 Mon Sep 17 00:00:00 2001 From: kanschat Date: Thu, 26 Jun 2014 19:30:26 +0000 Subject: [PATCH] add grad-div residual git-svn-id: https://svn.dealii.org/trunk@33090 0785d39b-7218-0410-832d-ea1e28bc413d --- .../include/deal.II/integrators/divergence.h | 39 +++++++++++++++++++ 1 file changed, 39 insertions(+) diff --git a/deal.II/include/deal.II/integrators/divergence.h b/deal.II/include/deal.II/integrators/divergence.h index 79438a8fba..4f35e8f3e1 100644 --- a/deal.II/include/deal.II/integrators/divergence.h +++ b/deal.II/include/deal.II/integrators/divergence.h @@ -519,6 +519,45 @@ namespace LocalIntegrators } } + /** + * The weak form of the grad-div residual + * @f[ + * \int_Z \nabla\!\cdot\!u \nabla\!\cdot\!v \,dx + * @f] + * + * @author Guido Kanschat + * @date 2014 + */ + template + void grad_div_residual ( + Vector &result, + const FEValuesBase &fetest, + const VectorSlice > > > &input, + const double factor = 1.) + { + const unsigned int n_dofs = fetest.dofs_per_cell; + + AssertDimension(fetest.get_fe().n_components(), dim); + AssertVectorVectorDimension(input, dim, fetest.n_quadrature_points); + + for (unsigned int k=0; k