41 template <int dim, int fe_degree, int n_q_points_1d, int n_components, typename Number>
46   operator()(Portable::FEEvaluation<dim, fe_degree, n_q_points_1d, n_components, Number> *fe_eval,  in operator()
55 template <int dim, int fe_degree, int n_q_points_1d, int n_components, typename Number>
60   operator()(Portable::FEEvaluation<dim, fe_degree, n_q_points_1d, n_components, Number> *fe_eval,  in operator()
69 template <int dim, int fe_degree, int n_q_points_1d, int n_components, typename Number>
74   operator()(const typename Portable::MatrixFree<dim, Number>::Data *data,  in operator()
78     Portable::FEEvaluation<dim, fe_degree, n_q_points_1d, n_components, Number> fe_eval(data);  in operator()
82       OperatorDealiiMassQuad<dim, fe_degree, n_q_points_1d, n_components, Number>());  in operator()
87   static const unsigned int n_local_dofs = Utilities::pow(fe_degree + 1, dim) * n_components;
88   static const unsigned int n_q_points   = Utilities::pow(n_q_points_1d, dim);
93 template <int dim, int fe_degree, int n_q_points_1d, int n_components, typename Number>
98   operator()(const typename Portable::MatrixFree<dim, Number>::Data *data,  in operator()
102     Portable::FEEvaluation<dim, fe_degree, n_q_points_1d, n_components, Number> fe_eval(data);  in operator()
106       OperatorDealiiLaplaceQuad<dim, fe_degree, n_q_points_1d, n_components, Number>());  in operator()
111   static const unsigned int n_local_dofs = Utilities::pow(fe_degree + 1, dim) * n_components;
112   static const unsigned int n_q_points   = Utilities::pow(n_q_points_1d, dim);
120 template <int dim, typename Number>
129   OperatorDealii(const Mapping<dim>              &mapping,  in OperatorDealii()
130                  const DoFHandler<dim>           &dof_handler,  in OperatorDealii()
132                  const Quadrature<dim>           &quadrature,  in OperatorDealii()
155     typename Portable::MatrixFree<dim, Number>::AdditionalData additional_data;  in reinit()
184     else if (n_components == dim && fe_degree == 1 && n_q_points_1d == 2)  in vmult()
185       this->vmult_internal<dim, 1, 2>(dst, src);  in vmult()
186     else if (n_components == dim && fe_degree == 2 && n_q_points_1d == 3)  in vmult()
187       this->vmult_internal<dim, 2, 3>(dst, src);  in vmult()
192     else if (n_components == dim && fe_degree == 1 && n_q_points_1d == 3)  in vmult()
193       this->vmult_internal<dim, 1, 3>(dst, src);  in vmult()
194     else if (n_components == dim && fe_degree == 2 && n_q_points_1d == 4)  in vmult()
195       this->vmult_internal<dim, 2, 4>(dst, src);  in vmult()
227     else if (n_components == dim && fe_degree == 1 && n_q_points_1d == 2)  in compute_inverse_diagonal()
228       this->compute_inverse_diagonal_internal<dim, 1, 2>(diagonal);  in compute_inverse_diagonal()
229     else if (n_components == dim && fe_degree == 2 && n_q_points_1d == 3)  in compute_inverse_diagonal()
230       this->compute_inverse_diagonal_internal<dim, 2, 3>(diagonal);  in compute_inverse_diagonal()
235     else if (n_components == dim && fe_degree == 1 && n_q_points_1d == 3)  in compute_inverse_diagonal()
236       this->compute_inverse_diagonal_internal<dim, 1, 3>(diagonal);  in compute_inverse_diagonal()
237     else if (n_components == dim && fe_degree == 2 && n_q_points_1d == 4)  in compute_inverse_diagonal()
238       this->compute_inverse_diagonal_internal<dim, 2, 4>(diagonal);  in compute_inverse_diagonal()
253         OperatorDealiiMassLocal<dim, fe_degree, n_q_points_1d, n_components, Number> mass_operator;  in vmult_internal()
258         OperatorDealiiLaplaceLocal<dim, fe_degree, n_q_points_1d, n_components, Number>  in vmult_internal()
273         OperatorDealiiMassQuad<dim, fe_degree, n_q_points_1d, n_components, Number> op_quad;  in compute_inverse_diagonal_internal()
275         MatrixFreeTools::compute_diagonal<dim, fe_degree, n_q_points_1d, n_components, Number>(  in compute_inverse_diagonal_internal()
280         OperatorDealiiLaplaceQuad<dim, fe_degree, n_q_points_1d, n_components, Number> op_quad;  in compute_inverse_diagonal_internal()
282         MatrixFreeTools::compute_diagonal<dim, fe_degree, n_q_points_1d, n_components, Number>(  in compute_inverse_diagonal_internal()
299   const Mapping<dim> &mapping;
304   const DoFHandler<dim> &dof_handler;
314   const Quadrature<dim> &quadrature;
324   Portable::MatrixFree<dim, Number> matrix_free;