197   Mat          Ap; /* dfdp */  in main()  local
238   PetscCall(MatCreate(PETSC_COMM_WORLD, &Ap));  in main()
239   PetscCall(MatSetSizes(Ap, n, 1, PETSC_DETERMINE, PETSC_DETERMINE));  in main()
240   PetscCall(MatSetType(Ap, MATDENSE));  in main()
241   PetscCall(MatSetFromOptions(Ap));  in main()
242   PetscCall(MatSetUp(Ap));  in main()
243   PetscCall(MatZeroEntries(Ap)); /* initialize to zeros */  in main()
257   PetscCall(TSSetRHSJacobianP(ts, Ap, RHSJacobianP, &app));  in main()
302   PetscCall(MatCreateVecs(Ap, &mu[0], NULL));  in main()
303   PetscCall(MatCreateVecs(Ap, &mu[1], NULL));  in main()
330   PetscCall(MatDestroy(&Ap));  in main()