Dear all, sorry for the second message on the same day. I need interpolate the Jacobian matrix of at the center of an element. I wrote the following part but I am not sure whether it is right (because this is my first time to use FE_nothing)
...DoFHandler<dim>::active_cell_iterator &cell ... // update the Jacobian matrix at the center of the cell Point<dim> center_point; for (unsigned int i = 0; i < dim; ++i) center_point[i] = 0.5; Quadrature<dim> center_quadrature(center_point); FE_Nothing<dim> fe_nothing; MappingQ1<dim> mapping; FEValues<dim> temp_fe_values(mapping, fe_nothing, center_quadrature, update_jacobians); temp_fe_values.reinit(cell); DerivativeForm<1, dim, dim> center_jacobian = temp_fe_values.jacobian(0); double det_J0 = center_jacobian.determinant(); Best Regards, Leo -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion visit https://groups.google.com/d/msgid/dealii/ae83ae55-840a-42da-b291-da68aff04a84n%40googlegroups.com.
