public inbox for gsl-discuss@sourceware.org
 help / color / mirror / Atom feed
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.
>
>

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