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
> [email protected]
> http://www.fenics.org/mailman/listinfo/dolfin-dev
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
