11   subroutine mymatmult(A, x, y, ierr)  argument
14     PetscErrorCode ierr
19   subroutine mymatmultadd(A, x, y, z, ierr)  argument
22     PetscErrorCode ierr
27   subroutine mymatmulttranspose(A, x, y, ierr)  argument
30     PetscErrorCode ierr
35   subroutine mymatmulthermitiantranspose(A, x, y, ierr)  argument
38     PetscErrorCode ierr
43   subroutine mymatmulttransposeadd(A, x, y, z, ierr)  argument
46     PetscErrorCode ierr
51   subroutine mymatmulthermitiantransposeadd(A, x, y, z, ierr)  argument
54     PetscErrorCode ierr
59   subroutine mymattranspose(A, reuse, B, ierr)  argument
62     PetscErrorCode ierr
67     PetscCallA(MatCreateShell(PETSC_COMM_SELF, i12, i12, i12, i12, i0, B, ierr))
68     PetscCallA(MatAssemblyBegin(B, MAT_FINAL_ASSEMBLY, ierr))
69     PetscCallA(MatAssemblyEnd(B, MAT_FINAL_ASSEMBLY, ierr))
74   subroutine mymatgetdiagonal(A, x, ierr)  argument
77     PetscErrorCode ierr
82   subroutine mymatdiagonalscale(A, x, y, ierr)  argument
85     PetscErrorCode ierr
90   subroutine mymatzeroentries(A, ierr)  argument
92     PetscErrorCode ierr
97   subroutine mymataxpy(A, alpha, B, str, ierr)  argument
101     PetscErrorCode ierr
106   subroutine mymatshift(A, alpha, ierr)  argument
109     PetscErrorCode ierr
114   subroutine mymatdiagonalset(A, x, ins, ierr)  argument
118     PetscErrorCode ierr
123   subroutine mymatdestroy(A, ierr)  argument
125     PetscErrorCode ierr
130   subroutine mymatview(A, viewer, ierr)  argument
133     PetscErrorCode ierr
138   subroutine mymatgetvecs(A, x, y, ierr)  argument
141     PetscErrorCode ierr
158   PetscErrorCode ierr
161   PetscCallA(PetscInitialize(ierr))
166   PetscCallA(VecCreateSeq(PETSC_COMM_SELF, i12, x, ierr))
167   PetscCallA(VecCreateSeq(PETSC_COMM_SELF, i12, y, ierr))
168   PetscCallA(VecCreateSeq(PETSC_COMM_SELF, i12, z, ierr))
169   PetscCallA(MatCreateShell(PETSC_COMM_SELF, i12, i12, i12, i12, i0, m, ierr))
170   PetscCallA(MatShellSetManageScalingShifts(m, ierr))
171   PetscCallA(MatAssemblyBegin(m, MAT_FINAL_ASSEMBLY, ierr))
172   PetscCallA(MatAssemblyEnd(m, MAT_FINAL_ASSEMBLY, ierr))
175   PetscCallA(MatShellSetOperation(m, op, mymatmult, ierr))
177   PetscCallA(MatShellSetOperation(m, op, mymatmultadd, ierr))
179   PetscCallA(MatShellSetOperation(m, op, mymatmulttranspose, ierr))
181   PetscCallA(MatShellSetOperation(m, op, mymatmulthermitiantranspose, ierr))
183   PetscCallA(MatShellSetOperation(m, op, mymatmulttransposeadd, ierr))
185   PetscCallA(MatShellSetOperation(m, op, mymatmulthermitiantransposeadd, ierr))
187   PetscCallA(MatShellSetOperation(m, op, mymattranspose, ierr))
189   PetscCallA(MatShellSetOperation(m, op, mymatgetdiagonal, ierr))
191   PetscCallA(MatShellSetOperation(m, op, mymatdiagonalscale, ierr))
193   PetscCallA(MatShellSetOperation(m, op, mymatzeroentries, ierr))
195   PetscCallA(MatShellSetOperation(m, op, mymataxpy, ierr))
197   PetscCallA(MatShellSetOperation(m, op, mymatshift, ierr))
199   PetscCallA(MatShellSetOperation(m, op, mymatdiagonalset, ierr))
201   PetscCallA(MatShellSetOperation(m, op, mymatdestroy, ierr))
203   PetscCallA(MatShellSetOperation(m, op, mymatview, ierr))
205   PetscCallA(MatShellSetOperation(m, op, mymatgetvecs, ierr))
207   PetscCallA(MatMult(m, x, y, ierr))
208   PetscCallA(MatMultAdd(m, x, y, z, ierr))
209   PetscCallA(MatMultTranspose(m, x, y, ierr))
210   PetscCallA(MatMultHermitianTranspose(m, x, y, ierr))
211   PetscCallA(MatMultTransposeAdd(m, x, y, z, ierr))
212   PetscCallA(MatMultHermitianTransposeAdd(m, x, y, z, ierr))
213   PetscCallA(MatTranspose(m, MAT_INITIAL_MATRIX, mt, ierr))
214   PetscCallA(MatGetDiagonal(m, x, ierr))
215   PetscCallA(MatDiagonalScale(m, x, y, ierr))
216   PetscCallA(MatZeroEntries(m, ierr))
218   PetscCallA(MatAXPY(m, a, mt, SAME_NONZERO_PATTERN, ierr))
219   PetscCallA(MatShift(m, a, ierr))
220   PetscCallA(MatDiagonalSet(m, x, INSERT_VALUES, ierr))
221   PetscCallA(MatView(m, viewer, ierr))
222   PetscCallA(MatCreateVecs(m, x, y, ierr))
223   PetscCallA(MatDestroy(m, ierr))
224   PetscCallA(MatDestroy(mt, ierr))
225   PetscCallA(VecDestroy(x, ierr))
226   PetscCallA(VecDestroy(y, ierr))
227   PetscCallA(VecDestroy(z, ierr))
229   PetscCallA(PetscFinalize(ierr))