dolfin team mailing list archive

[Andy R Terrel <question129481@xxxxxxxxxxxxxxxxxxxxx>]

New question #129481 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/129481

I need to write out some assembled matrices in Matrix Market format, hopefully from the Python interface.  Easily for small matrices I can do:

from dolfin import *
from scipy.io import mmio
from scipy.sparse import csr_matrix
...
       A = Matrix()
        b = Vector()

        assemble_system(a, L, bcs=bc0, A_tensor=A, b_tensor=b)

        mmio.mmwrite("A.mtx", csr_matrix(A.array()))
        mmio.mmwrite("b.mtx", csr_matrix(b.array()))

But once A gets large calling "A.array()" takes forever.  Is there a way in python that I can just query the sparsity structure of A?

Any other ideas on how to do this quickly?



-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.