Hello, 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. Thank you, Alex -- 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/6ce3d2a3-5507-42e3-bd22-6be4484edd9do%40googlegroups.com.