On Fri, 22 Nov 2019, Renato Poli wrote: > I am alternatively considering a global L2 Glaerkin mapping. > > Something like - find f_i where: > f(x,y)=f_i \phi_i (new mesh) > g(x,y) (function evaluated in the old mesh) > \int f(x,y) \phi dA = \int g(x,y) \phi dA > > It make take more time,
It may take more CPU time but I bet it'll be *much* cheaper in programmer-hours. > but it looks safe ... doesn't it? That depends on your problem and your formulation, I think. Global projections can be a great way to maintain conservation, but they're prone to "ringing" artifacts, and either one of those factors can be critically important or nearly inconsequential depending on what you're solving and how. --- Roy _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
