← Back to team overview

fenics team mailing list archive

sparse matrix/parallel

 

I have a sparse matrix that I'm constructing and beforehand I'm setting the sparsity pattern and then I will populate the matrix. Obviously dolfin is not going to handle everything so I ask this:

Will the "right" way to do this, in parallel, be to separately use the underlying PETSc routines for setting the values using MatGetOwnershipRange and MatSetValues?

Also, after the appropriate entries in sparsity_pattern are set, should I then call sparsity_pattern.apply() ?


-gideon


Follow ups