Hello,
For dense matrices, submatrix views work because all the matrix
elements are stored in a nice consecutive, linear manner. This is not
the case for sparse matrices. GSL additionally supports three different
storage formats for sparse matrices (COO, CSR, CSC). So it is not
trivial to implement submatrix views for sparse matrices.
You could simply construct separate sparse matrices for each subsystem,
and then add them all together in the final step.
Patrick
On 12/9/22 01:00, Brijesh Upadhaya wrote:
Hi,
I am trying to solve a coupled electromagnetic field problem. I need to
combine several subsystems into one large system. Is there any idea to get
the submatrix view of the gsl_spmatrix?
It seems to me that the submatrix view is not yet implemented in the
development version.
BR,
Brijesh
