← Back to team overview

dolfin team mailing list archive

Re: [HG DOLFIN] Make rank() and dim() virtual (overloadable by user).

 

2007/9/17, Garth N. Wells <gnw20@xxxxxxxxx>:
> Alessio Quaglino wrote:
> >> Alessio Quaglino wrote:
> >>> With "the mesh remains constant" is meant that the vertices cannot move
> >>> or
> >>> that there should not be any connectivity change?
> >>>
> >> No connectivity changes.
> >
> > I see that in
> >
> > void Assembler::assemble(GenericTensor& A, const Form& form,
> >                          bool reset_tensor)
> > {
> >   assemble(A, form.form(), form.coefficients(), 0, 0, 0);
> > }
> >
> > reset_tensor is not used. Is this wanted? The initGlobalTensor takes 90%
> > of my runtime (I think it wasn't like this before the last release).
> >
>
>
> That's exactly why it's there. If the matrix has already been
> initialised, then there is not need to re-initialise. This will
> typically be useful for nonlinear problems.

But reset_tensor isn't used here, it's an obvious bug.

-- 
Martin


Follow ups

References