dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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

Attachment: signature.asc
Description: Digital signature