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
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