>> 2. What functionality would you like to see added to GSL? > Other interplants, particularly pchip. Interfaces to sparse direct linear > solvers like umfpack. >
The current git has routines for constructing sparse matrices (ie: gsl_spmatrix_{get,set}) for both triplet and compressed storage formats. This should make it fairly simple to pass the relevant arrays to other sparse solver libraries. I've personally interfaced gsl_spmatrix into the Lis, pastix, and mumps libraries (haven't tried umfpack). All thats needed is a simple wrapper routine to call the external library with gsl_spmatrix.