#include "petscmat.h" PetscErrorCode MatCreateKAIJ(Mat A,PetscInt p,PetscInt q,const PetscScalar S[],const PetscScalar T[],Mat *kaij)[I \otimes S + A \otimes T]
where S is a dense (p \times q) matrix T is a dense (p \times q) matrix A is an AIJ (n \times n) matrix I is the identity matrix The resulting matrix is (np \times nq)
S and T are always stored independently on all processes as PetscScalar arrays in column-major format.
Collective
| A | - the AIJ matrix | |
| p | - number of rows in S and T | |
| q | - number of columns in S and T | |
| S | - the S matrix (can be PETSC_NULL), stored as a PetscScalar array (column-major) | |
| T | - the T matrix (can be PETSC_NULL), stored as a PetscScalar array (column-major) |
| kaij | - the new KAIJ matrix |