dolfin team mailing list archive

["Alessio Quaglino" <alessio@xxxxxx>]

When I assemble a matrix then my program crashes if I multiply it with a
vector. This *seems* not to happen if the matrix mantains the same
structure after the assembly (if for example it was diagonal and remains
diagonal), but only when its structure is modified by the assembler. Just
doing:

A.init(N,M);
assembler.assemble(A, a, true);
A.mult(x,y);

I get the crash. While if I do:

A.init(N,M);
B.init(N,M)
assembler.assemble(A, a, true);
B = A;
B.mult(x,y);

everything is fine.

Alessio