Alexander,

I am wondering if anybody has also found that the inverse of the Jacobian fromFE Values, with MappingQGeneric does not satisfy the Geometric ConservationLaw (GCL), in the sense of:Kopriva, David A. "Metric identities and the discontinuous spectral elementmethod 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 theJacobian 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?`

Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/d8421b9f-0fa9-9769-b63c-281500879687%40colostate.edu.