public inbox for gsl-discuss@sourceware.org
 help / color / mirror / Atom feed
* Re: Please implement typical functionality for sparse matrices
       [not found] <00f77ff86fbe3066184b5cb827abd6ee5dbc637c.camel@gmail.com>
@ 2018-07-13  4:01 ` Patrick Alken
  0 siblings, 0 replies; only message in thread
From: Patrick Alken @ 2018-07-13  4:01 UTC (permalink / raw)
  To: David Cortes, gsl-discuss

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.
>
>

^ permalink raw reply	[flat|nested] only message in thread

only message in thread, other threads:[~2018-07-13  4:01 UTC | newest]

Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
     [not found] <00f77ff86fbe3066184b5cb827abd6ee5dbc637c.camel@gmail.com>
2018-07-13  4:01 ` Please implement typical functionality for sparse matrices Patrick Alken

This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).