Dear Luca I wanted to use MappingManifold for sphere, but it looks like this is not possible either with PolarManifold or SphericalManifold.
I am doing this for a shallow water model on sphere, where I want to do long time integration, and look at conservation of mass, energy, etc. Having an exact mapping is important for conservation. If I use *SphericalManifold + MappingQ*, I seem to have loss of mass conservation over long time integration. In this case, can you tell me about these questions ? (1) If phi_i are the basis functions, does grad(phi_i) lie exactly on the sphere ? (I guess not) (2) If (x_q,y_q,z_q) is some quadrature point, is sum_i grad(phi_i(x_q,y_q,z_q)) = 0 upto machine precision ? The last property is necessary for mass conservation in DG schemes for conservation laws. Best praveen -- 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.
