Actual source code: ex133.c


  2: static char help[] = "Test saving SeqSBAIJ matrix that is missing diagonal entries.";

  4: #include <petscmat.h>

  6: int main(int argc, char **args)
  7: {
  8:   Mat         A;
  9:   PetscInt    bs = 3, m = 4, i, j, val = 10, row[2], col[3], rstart;
 10:   PetscMPIInt size;
 11:   PetscScalar x[6][9];

 14:   PetscInitialize(&argc, &args, (char *)0, help);
 15:   MPI_Comm_size(PETSC_COMM_WORLD, &size);
 17:   MatCreateSeqSBAIJ(PETSC_COMM_SELF, bs, m * bs, m * bs, 1, NULL, &A);
 18:   MatSetOption(A, MAT_IGNORE_LOWER_TRIANGULAR, PETSC_TRUE);
 19:   rstart = 0;

 21:   row[0] = rstart + 0;
 22:   row[1] = rstart + 2;
 23:   col[0] = rstart + 0;
 24:   col[1] = rstart + 1;
 25:   col[2] = rstart + 3;
 26:   for (i = 0; i < 6; i++) {
 27:     for (j = 0; j < 9; j++) x[i][j] = (PetscScalar)val++;
 28:   }
 29:   MatSetValuesBlocked(A, 2, row, 3, col, &x[0][0], INSERT_VALUES);
 30:   MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY);
 31:   MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY);
 32:   MatView(A, PETSC_VIEWER_BINARY_WORLD);

 34:   MatDestroy(&A);
 35:   PetscFinalize();
 36:   return 0;
 37: }