The short answer is: assemble the sparse matrix in triplet format, then use 
`sparse` to convert to CCS format. Assembling in CCS format (what you are 
doing) is almost generically slow.

Christoph


Reply via email to