For nonaffine geometries it would be much cleaner if the mesh coordinates were always represented by a Function.
In the initial iteration this can be made working by essentially ignoring the mesh.geometry() in the assembly, but to get that working it must be special cased in various places around dolfin. So it would be much nicer if MeshGeometry could simply be a Function. Or optionally have a Function, maybe that's easier to do. The problem is circular dependencies, which are in a way caused by the grouping of geometry and topology in the mesh. I don't really have a question at this point, just want to let everyone know that this is an issue that needs considering from various angles. Martin
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
