From: Roy S. <roy...@ic...> - 2012-08-17 16:27:09
|
On Fri, 17 Aug 2012, Subramanya Gautam Sadasiva wrote: > I am trying to write a small multiphysics code with Libmesh and part > of it needs a very dense mesh while other parts can do with a less > dense mesh. However, the same fields need to be accessible across > all the meshes. Is it possible to save both a refined and coarsened > mesh state at the same time and transfer results between them? You can use two separate Mesh objects, two EquationSystems, and a MeshFunction object to do this. The scaling won't be great, though - because of the point locator required by MeshFunction, your assembly would become O(N_E*log(N_E)) instead of just O(N_E). If you're running a non-adaptive time-dependent problem (or multiple forcing functions, or some other reason to do multiple assemblies on the same mesh) you could cache mesh-to-mesh mappings in a hash map to amortize the lookups. --- Roy |