636   const PetscScalar *ptr_X, *ptr_Y;  in SolutionErrorNorms()  local
646   PetscCall(VecGetArrayRead(X, &ptr_X));  in SolutionErrorNorms()
649 …if (i < count_slow / 2 || i > count_slow / 2 + count_fast - 1) *nrm1 += hs * PetscAbs(ptr_X[i] - p…  in SolutionErrorNorms()
650     else *nrm1 += hf * PetscAbs(ptr_X[i] - ptr_Y[i]);  in SolutionErrorNorms()
652   PetscCall(VecRestoreArrayRead(X, &ptr_X));  in SolutionErrorNorms()
767     const PetscScalar *ptr_X, *ptr_X0;  in main()  local
777     PetscCall(DMDAVecGetArrayRead(da, X, (void *)&ptr_X));  in main()
782           mass_final   = mass_final + hs * ptr_X[i * dof + k];  in main()
787           mass_final   = mass_final + hf * ptr_X[i * dof + k];  in main()
792     PetscCall(DMDAVecRestoreArrayRead(da, X, (void *)&ptr_X));  in main()
815       PetscCall(VecGetArrayRead(X, &ptr_X));  in main()
818 …slow / 2 || i > count_slow / 2 + count_fast - 1) nrm1 = nrm1 + hs * PetscAbs(ptr_X[i] - ptr_XR[i]);  in main()
819         else nrm1 = nrm1 + hf * PetscAbs(ptr_X[i] - ptr_XR[i]);  in main()
821       PetscCall(VecRestoreArrayRead(X, &ptr_X));  in main()