← Back to team overview

dolfin team mailing list archive

Re: [HG DOLFIN] Don't initialize vector in DiscreteFunction if already initialized.

 

On ma., 2008-09-15 at 11:19 +0100, Garth N. Wells wrote:
> 
> DOLFIN wrote:
> > One or more new changesets pushed to the primary dolfin repository.
> > A short summary of the last three changesets is included below.
> > 
> > changeset:   4751:df802d1b2b349e91a51b70f799ba6d0d3b825376
> > tag:         tip
> > user:        Anders Logg <logg@xxxxxxxxx>
> > date:        Sat Sep 13 08:55:34 2008 +0200
> > files:       dolfin/function/DiscreteFunction.cpp
> > description:
> > Don't initialize vector in DiscreteFunction if already initialized.
> >
> 
> Vector::init() only resizes the vector if necessary, but it also zeroes 
> the vector. Looks to me like it is desirable to zero the vector?
> 
> Garth
> 

The reason for this, I think, is initialization as such: 

solve(A, x, b, cg, amg)
u = Function(element, mesh, x)
plot(u)

Kent



Follow ups

References