[petsc-users] matrix-vector and matrix-matrix products in PETSc

[Zhang, Hong]

"Sparse Matrix-Matrix Products Executed Through Coloring" paper is a sequential algorithm for C=A*B^T.

We do not have paper discussing the parallel mat-mat product algorithms used in PETSc.
For source code, you may look at MatMatMultSymbolic_MPIAIJ_MPIAIJ() and MatMatMultNumeric_MPIAIJ_MPIAIJ()
in petsc/src/mat/impls/aij/mpi/mpimatmatmult.c.
Hong
________________________________
From: petsc-users <petsc-users-bounces at mcs.anl.gov> on behalf of Jose E. Roman <jroman at dsic.upv.es>
Sent: Thursday, November 11, 2021 7:33 AM
To: Matthew Knepley <knepley at gmail.com>
Cc: PETSc <petsc-users at mcs.anl.gov>; Majid Rasouli <maj.rasouli at gmail.com>
Subject: Re: [petsc-users] matrix-vector and matrix-matrix products in PETSc



> El 11 nov 2021, a las 13:41, Matthew Knepley <knepley at gmail.com> escribió:
>
> On Thu, Nov 11, 2021 at 3:00 AM Majid Rasouli <maj.rasouli at gmail.com> wrote:
> Dear all,
>
> I'm trying to find a brief explanation about the parallel (MPI) implementation of matrix-vector multiplication and matrix-matrix multiplication in PETSc. I have used them for some experiments for a research paper and I need to provide a brief explanation about the implementation of them. Is there a paper that I can read about that? If not, can anyone please help me find the source code that includes the implementation of those two functions?
>
> I think this describes our SpMV:
>
>   https://dl.acm.org/doi/abs/10.1145/3225058.3225100?casa_token=Pz_zMwXCoz8AAAAA%3AH8iNvXSY4opzDsJsEnoYVU0CccQjOK9pIH4Rs_Yb7a_ebchRpgAJ5GHoH2humWfkFYS54m2fO84s
>
> This paper describes the SpMM algorithms we use, and makes some notes on performance
>
>   https://epubs.siam.org/doi/abs/10.1137/15M1028807
>
> but I could not find a paper on our implementation.

Sparse Matrix-Matrix Products Executed Through Coloring
Michael McCourt, Barry Smith, and Hong Zhang
https://doi.org/10.1137/13093426X



>
>   Thanks,
>
>      Matt
>
> Thanks a lot,
> Majid.
>
>
> --
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211111/e12aff95/attachment.html>