petsc-3.7.3 2016-08-01

MatCreateComposite

Creates a matrix as the sum of zero or more matrices

#include "petscmat.h" 
PetscErrorCode  MatCreateComposite(MPI_Comm comm,PetscInt nmat,const Mat *mats,Mat *mat)

Collective on MPI_Comm

A -the matrix

Alternative construction

      MatCreate(comm,&mat);

      MatSetSizes(mat,m,n,M,N);

      MatSetType(mat,MATCOMPOSITE);

      MatCompositeAddMat(mat,mats[0]);

      ....

      MatCompositeAddMat(mat,mats[nmat-1]);

      MatAssemblyBegin(mat,MAT_FINAL_ASSEMBLY);

      MatAssemblyEnd(mat,MAT_FINAL_ASSEMBLY);

For the multiplicative form the product is mat[nmat-1]*mat[nmat-2]*....*mat[0]