From: kanschat Date: Thu, 26 Sep 2013 08:19:33 +0000 (+0000) Subject: weak residuals X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=10455003df248cce46de25b56f892a3a9fcf8133;p=dealii-svn.git weak residuals git-svn-id: https://svn.dealii.org/trunk@30941 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/include/deal.II/integrators/divergence.h b/deal.II/include/deal.II/integrators/divergence.h index b78c4777f0..ef99627c83 100644 --- a/deal.II/include/deal.II/integrators/divergence.h +++ b/deal.II/include/deal.II/integrators/divergence.h @@ -147,19 +147,58 @@ namespace LocalIntegrators result(i) += dx * input[d][k][d] * fetest.shape_value(i,k); } } + + /** - * Cell matrix for gradient. The derivative is on the trial - * function. + * The residual of the divergence operator in weak form. * \f[ - * \int_Z \nabla u \cdot \mathbf v\,dx + * - \int_Z \nabla v \cdot \mathbf u \,dx * \f] - * This is the strong gradient and the trial space - * should be at least in H1. The test functions can - * be discontinuous. + * This is the weak divergence operator and the test + * space should be at least H1. The trial functions + * may be discontinuous. + * + * @todo Verify: The function cell_matrix() is the Frechet derivative of this function with respect + * to the test functions. * * @author Guido Kanschat - * @date 2011 + * @date 2013 */ + template + void cell_residual( + Vector &result, + const FEValuesBase &fetest, + const VectorSlice > > &input, + const double factor = 1.) + { + AssertDimension(fetest.get_fe().n_components(), 1); + AssertVectorVectorDimension(input, dim, fetest.n_quadrature_points); + const unsigned int t_dofs = fetest.dofs_per_cell; + Assert (result.size() == t_dofs, ExcDimensionMismatch(result.size(), t_dofs)); + + for (unsigned int k=0; kH1. The test functions can be discontinuous. + * + * @author Guido Kanschat + * @date 2011 + */ template void gradient_matrix( FullMatrix &M, @@ -228,6 +267,42 @@ namespace LocalIntegrators } } + /** + * The residual of the gradient operator in weak form. + * \f[ + * -\int_Z \nabla\cdot \mathbf v u \,dx + * \f] + * This is the weak gradient operator and the test + * space should be at least Hdiv. The trial functions + * may be discontinuous. + * + * @todo Verify: The function gradient_matrix() is the Frechet derivative of this function with respect to the test functions. + * + * @author Guido Kanschat + * @date 2013 + */ + template + void gradient_residual( + Vector &result, + const FEValuesBase &fetest, + const std::vector &input, + const double factor = 1.) + { + AssertDimension(fetest.get_fe().n_components(), dim); + AssertDimension(input.size(), fetest.n_quadrature_points); + const unsigned int t_dofs = fetest.dofs_per_cell; + Assert (result.size() == t_dofs, ExcDimensionMismatch(result.size(), t_dofs)); + + for (unsigned int k=0; k