28   Vec           gauss, input;  in main()  local
47   PetscCall(DMCreateGlobalVector(da2D, &gauss));  in main()
48   PetscCall(PetscObjectSetName((PetscObject)gauss, "pressure"));  in main()
51   PetscCall(VecSet(gauss, value));  in main()
57   PetscCall(DMDAVecGetArray(da2D, gauss, &gauss_ptr));  in main()
61   PetscCall(DMDAVecRestoreArray(da2D, gauss, &gauss_ptr));  in main()
68   PetscCall(VecView(gauss, H5viewer));  in main()
73   PetscCall(VecDuplicate(gauss, &input));  in main()
74   PetscCall(PetscObjectGetName((PetscObject)gauss, &vecname));  in main()
83   PetscCall(VecAXPY(input, -1.0, gauss));  in main()
88   PetscCall(VecDestroy(&gauss));  in main()