MatRARt#

Creates the matrix product \(C = R * A * R^T\)

Synopsis#

#include "petscmat.h" 
PetscErrorCode MatRARt(Mat A, Mat R, MatReuse scall, PetscReal fill, Mat *C)

Neighbor-wise Collective

Input Parameters#

Output Parameter#

  • C - the product matrix

Notes#

C will be created and must be destroyed by the user with MatDestroy().

An alternative approach to this function is to use MatProductCreate() and set the desired options before the computation is done

This routine is currently only implemented for pairs of MATAIJ matrices and classes which inherit from MATAIJ. Due to PETSc sparse matrix block row distribution among processes, parallel MatRARt() is implemented via explicit transpose of R, which could be very expensive. We recommend using MatPtAP().

The deprecated PETSC_DEFAULT in fill also means use the current value

See Also#

Matrices, Mat, MatProductCreate(), MatMatMult(), MatPtAP()

Level#

intermediate

Location#

src/mat/interface/matrix.c


Index of all Mat routines
Table of Contents for all manual pages
Index of all manual pages