from Brijesh Upadhaya <brijesh =dot= upadhaya =at= gmail =dot= com>

I don't know if this is necessary, but the scaling of a complex (sparse) matrix with a real vector, complex sparse matrix-vector multiplication, and unpacking(in place) a complex vector could be added. I had to somehow do the above-mentioned stuff while solving a time-harmonic FE problem.

There are alloc and calloc for the blocks and vectors. Is there any particular reason for not having anything to resize a block or a vector?

In gsl/spmatrix/getset_source.c, _set() replaces the matrix elements, could it be possible to assemble the matrix as

(BASE) ptr = (BASE) + x;

and copy the entire code to a new name.