From: Martin Kronbichler Date: Mon, 6 Aug 2018 15:52:12 +0000 (+0200) Subject: Add transform_from_q_points_to_basis in inverse mass matrix. X-Git-Tag: v9.1.0-rc1~841^2~3 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=aa3c5d75ad5c9f5844f3543fb42a1ecfd9d92f2f;p=dealii.git Add transform_from_q_points_to_basis in inverse mass matrix. --- diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h index 9f5d20f4a8..cd331f38fc 100644 --- a/include/deal.II/matrix_free/operators.h +++ b/include/deal.II/matrix_free/operators.h @@ -636,6 +636,44 @@ namespace MatrixFreeOperators const VectorizedArray * in_array, VectorizedArray * out_array) const; + /** + * This operation performs a projection from the data given in quadrature + * points to the actual basis underlying this object. This projection can + * also be interpreted as a change of the basis from the Lagrange + * interpolation polynomials in the quadrature points to the + * basis underlying the current `fe_eval` object. + * + * Calling this function on an array as + * @code + * inverse_mass.transform_from_q_points_to_basis(1, array, + * phi.begin_dof_values()); + * @endcode + * is equivalent to + * @code + * for (unsigned int q=0; q *in_array, + VectorizedArray *out_array) const; + /** * Fills the given array with the inverse of the JxW values, i.e., a mass * matrix with coefficient 1. Non-unit coefficients must be multiplied (in @@ -999,6 +1037,49 @@ namespace MatrixFreeOperators } } + + + template + inline void + CellwiseInverseMassMatrix:: + transform_from_q_points_to_basis(const unsigned int n_actual_components, + const VectorizedArray *in_array, + VectorizedArray *out_array) const + { + const unsigned int dofs_per_cell = + Utilities::fixed_int_power::value; + internal::EvaluatorTensorProduct> + evaluator(AlignedVector>(), + AlignedVector>(), + inverse_shape); + + for (unsigned int d = 0; d < n_actual_components; ++d) + { + const VectorizedArray *in = in_array + d * dofs_per_cell; + VectorizedArray * out = out_array + d * dofs_per_cell; + + if (dim == 3) + { + evaluator.template hessians<2, true, false>(in, out); + evaluator.template hessians<1, true, false>(out, out); + evaluator.template hessians<0, true, false>(out, out); + } + if (dim == 2) + { + evaluator.template hessians<1, true, false>(in, out); + evaluator.template hessians<0, true, false>(out, out); + } + if (dim == 1) + evaluator.template hessians<0, true, false>(in, out); + } + } + + + //----------------- Base operator ----------------------------- template Base::Base()