dolfin team mailing list archive

[Jed Brown <jed@xxxxxxxx>]

On Tue 2008-08-19 11:59, Anders Logg wrote:
> On Thu, Aug 14, 2008 at 10:10:03PM +0000, Jed Brown wrote:
> > One way to implement this is to allocate a vector for Dirichlet values,
> > a vector for Homogeneous values, and a Combined vector.  The Homogeneous
> > vector is the only one that is externally visible.
> 
> Isn't this problematic? I want the entire vector visible externally
> (and not the homogeneous part). It would make it difficult to plot
> solutions, saving to file etc.
> 
> Maybe the Function class could handle the wrapping but it would involve a
> complication.

Right, by `externally visible' I mean to the solution process, that is
time-stepping, nonlinear solver, linear solvers, preconditioners.  The
vector you are concerned about is the post-processed state which you can
get with zero communication.  It is inherently tied to the mesh and
anything you do with it likely needs to know mesh connectivity.  I don't
think it is advantageous to lump this in with the global state vector.

Jed

Attachment: pgp1lpFweYJEQ.pgp
Description: PGP signature