1178   Mat       Jac = NULL;  /* Jacobian matrix                        */  in SolveODE()  local
1212     PetscCall(MatCreate(PETSC_COMM_WORLD, &Jac));  in SolveODE()
1213     PetscCall(MatSetSizes(Jac, PETSC_DECIDE, PETSC_DECIDE, N, N));  in SolveODE()
1214     PetscCall(MatSetFromOptions(Jac));  in SolveODE()
1215     PetscCall(MatSetUp(Jac));  in SolveODE()
1216     PetscCall(TSSetRHSJacobian(ts, Jac, Jac, RHSJacobian, &ptype[0]));  in SolveODE()
1221     PetscCall(MatCreate(PETSC_COMM_WORLD, &Jac));  in SolveODE()
1222     PetscCall(MatSetSizes(Jac, PETSC_DECIDE, PETSC_DECIDE, N, N));  in SolveODE()
1223     PetscCall(MatSetFromOptions(Jac));  in SolveODE()
1224     PetscCall(MatSetUp(Jac));  in SolveODE()
1225     PetscCall(TSSetIJacobian(ts, Jac, Jac, IJacobian, &ptype[0]));  in SolveODE()
1251   PetscCall(MatDestroy(&Jac));  in SolveODE()