From: Patrick Alken <alken@colorado.edu>
To: David Cortes <david.cortes.rivera@gmail.com>,
"gsl-discuss@sourceware.org" <gsl-discuss@sourceware.org>
Subject: Re: Please implement typical functionality for sparse matrices
Date: Fri, 13 Jul 2018 04:01:00 -0000 [thread overview]
Message-ID: <d8987205-8894-c2cd-165f-474caee3937e@colorado.edu> (raw)
In-Reply-To: <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.
>
>
parent reply other threads:[~2018-07-13 4:01 UTC|newest]
Thread overview: expand[flat|nested] mbox.gz Atom feed
[parent not found: <00f77ff86fbe3066184b5cb827abd6ee5dbc637c.camel@gmail.com>]
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=d8987205-8894-c2cd-165f-474caee3937e@colorado.edu \
--to=alken@colorado.edu \
--cc=david.cortes.rivera@gmail.com \
--cc=gsl-discuss@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
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).