On Mon, Nov 02, 2009 at 04:31:10PM +0100, Jed Brown wrote: > Anders Logg wrote: > > > The default update procedure would be interpolation. > > Which interpolation? Perhaps the user needs it to be conservative or > preserve some other structure. What if they want to use different > schemes for different functions on that mesh?
The standard interpolation defined by the finite element space. So for example nodal evaluation for Lagrange elements. > You can't implement refinement in-place anyway. Providing an interface > that looks in-place just means that users have to rewrite larger pieces > of their code when they need more control. I would suggest that > refinement produces a new mesh and then provide methods to transfer > functions and solver objects to the new mesh. Perhaps, but it could require quite a bit of work (from the user) to keep track of which objects should be transferred: function spaces, functions, forms. -- Anders
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list DOLFIN-dev@fenics.org http://www.fenics.org/mailman/listinfo/dolfin-dev