I am wondering if anybody has also found that the inverse of the Jacobian from FE Values, with MappingQGeneric does not satisfy the Geometric Conservation Law (GCL), in the sense of:

Kopriva, David A. "Metric identities and the discontinuous spectral element method on curvilinear meshes." /Journal of Scientific Computing/ 26.3 (2006): 301.

on curvilinear elements/manifolds in 3D.
That is:
\frac{\partial }{\partial \hat{x}_1} *det(J)* \frac{\partial \hat{x}_1 }{\partial x_1} + \frac{\partial }{\partial \hat{x}_2} *det(J)* \frac{\partial \hat{x}_2}{\partial x} + \frac{\partial }{\partial \hat{x}_3} * det(J)*\frac{\partial \hat{x}_3 }{\partial x_1} != 0 (GCL says it should =0, similarly for x_2 and x_3)

If so or if not, also, has anybody found a remedy to have the inverse of the Jacobian from FE Values with MappingQGeneric to satisfy the GCL.

I'm not sure any of us have ever thought about it. (I haven't -- but I really shouldn't speak for anyone else.) Can you explain what this equality represents? Why should it hold?

I'm also unsure whether we've ever checked whether it holds (exactly or approximately). Can you create a small test program that illustrates the behavior you are seeing?


Wolfgang Bangerth          email:       

The deal.II project is located at
For mailing list/forum options, see
--- 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 view this discussion on the web visit

Reply via email to