Adding sparse matrices

[rebpoli]

Hi

I am creating an additional matrix in a system, to avoid rebuilding the elastic part of the system every timestep.

I do that by:
my_mat = sys->add_matrix("elastic");
{ ... add to the matrix ... }
my_mat->close();

However, when I try to add to the global jacobian by:
jacobian->add( 1, *my_mat )

I get the following error:
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Object is in wrong state
[0]PETSC ERROR: Not for unassembled matrix (Y)

Shouldn't the "close()" call have prepared the matrix for addition?

[rebpoli]

The issue was closing the jacobian matrix. Solved.