dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Kent Andre wrote:
> Hi, 
> A.zero() is called both when reset_tensor is true and false
> in Assembler::init_global_tensor. 
>
> The code is as follows: 
>
> void Assembler::init_global_tensor(GenericTensor& A,
>                                    const Form& a,
>                                    UFC& ufc, bool reset_tensor)
> {
>   if (reset_tensor)
>   {
>     ... 
>   }
>   else 
>     A.zero();
> }
>
> Is this on purpose or is it a bug ? 

The code is correct. 'A.zero()' zeroes the entries but preserves the 
non-zero structure of the matrix.

Garth

> I thought 'A.zero()' should not happen when reset_tensor=False. 
>
> Kent
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev