This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: Please implement typical functionality for sparse matrices
- From: Patrick Alken <alken at colorado dot edu>
- To: David Cortes <david dot cortes dot rivera at gmail dot com>, "gsl-discuss at sourceware dot org" <gsl-discuss at sourceware dot org>
- Date: Thu, 12 Jul 2018 22:01:29 -0600
- Subject: Re: Please implement typical functionality for sparse matrices
- References: <00f77ff86fbe3066184b5cb827abd6ee5dbc637c.camel@gmail.com>
Moving this over to gsl-discuss since its not a bug.
I would like to expand support for sparse matrices also. Ideally I would
like to see GSL support the proposed Sparse BLAS standard
(https://math.nist.gov/spblas/) which would likely address the points
you raise below. I had some preliminary discussions a few months back
with the maintainer of librsb (a sparse matrix library which fully
implements the sparse BLAS standard, http://librsb.sourceforge.net/). In
principle we might be able to port a lot of that code since it is a LGPL
library. However the main issue is time.
I am currently the only active GSL developer and have quite limited time
these days. I would like to see more people work on these things and
contribute to GSL, but I see few people volunteer. If you can write good
code and know about sparse matrices, please consider working on these
issues yourself and contributing a patch. Otherwise it will likely take
me a long time to get to this :). But this is a relatively high priority
item for me since I work with sparse matrices frequently in my own work.
Thanks,
Patrick
On 07/12/2018 02:52 PM, David Cortes wrote:
> The function for matrix multiplication for sparse matrices
> (gsl_spblas_dgemm) currently only works when both inputs are in column-
> compressed format. It would be straightforward to expand it to cover
> also row-compressed matrices (without transposing them).
>
> As well, the spmatrix module supports convertion of
> coordinates/triplets formats to row- and column-compressed, but it
> doesn't support convertion from a compressed format to
> coordinates/triplets.
>
> This makes it also problematic to perform basic algebra on matrices
> which might later get multiplied.
>
>