Re: Sparse submatrix view

[Patrick Alken]

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