I am running with release mode at the moment (but this should not affect the computation, am I right?). I passed the "update_inverse_jacobians" flag to FEValues. I have only one cell for simplicity and it is a QGauss Q4 element with linear shape functions, hence with const QGauss<dim> quadrature_formula (2).
Kind regards, S. A. Mohseni -- 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]. For more options, visit https://groups.google.com/d/optout.
