p.s. the mapping is the same for both cases. On Monday, January 29, 2018 at 4:16:48 PM UTC+1, Denis Davydov wrote: > > Dear all, > > I am playing around with something in the following setup: > 1. I have a single p::d::Tria > 2. two different DoFHandler's using scalar FE_Q<dim> of different degrees > 3. The same quadrature rule (QGauss) used for both cases > 4. Two different MatrixFree objects > > I am not 100% sure, but it looks like the iteration over SIMD chunks of > cells via two different matrix-free > objects actually leads to the different sequence of cells. I see that when > I try to calculate some integral quantities > of my quadrature point data. Namely, i see the difference in the integral > > my_integral += quadrature_data(cell,q) * fe_eval.JxW(q); > > when I use two different MatrixFree objects to initialize FEEvaluation > (obviously with the same n_q_points_1d but different fe_degree). > The order of difference is far beyond any numerical roundoff error. > > > Is this something expected but not described, AFAICT, in documentation? > > > Regards, > Denis. >
-- 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.
