Sparsity pattern in BaSpaCho

[nikhil-vorticity]

❓ Questions and Help

Hello,

.Could you guide me to a reference about the sparsity pattern of the matrices and its factors?

[maurimo]

Hey @nikhil-vorticity
Sorry for the late reply!

Yes it would indeed be nice having such a function.
The best way to proceed is certainly from a duplicate of CoalescedBlockMatrixSkel::densify in baspacho/baspacho/CoalescedBlockMatrix.cpp. Basically we already have such a function (allowing to return the factor data either before/after factorization), but it only allows you to get the data as a filled dense matrix.
You can certainly modify it to add the written blocks to a sparse matrix, I don't have much time to do this myself but will be very happy to merge a diff adding this!

Best regards,
Maurizio