568 static PetscErrorCode THIInitializePrm(THI thi, DM da2prm, Vec prm)  in THIInitializePrm()  argument
576   PetscCall(DMDAVecGetArray(da2prm, prm, &p));  in THIInitializePrm()
583   PetscCall(DMDAVecRestoreArray(da2prm, prm, &p));  in THIInitializePrm()
650 static PetscErrorCode THIDAGetPrm(DM da, PrmNode ***prm)  in THIDAGetPrm()  argument
660   PetscCall(DMDAVecGetArray(da2prm, X, prm));  in THIDAGetPrm()
664 static PetscErrorCode THIDARestorePrm(DM da, PrmNode ***prm)  in THIDARestorePrm()  argument
674   PetscCall(DMDAVecRestoreArray(da2prm, X, prm));  in THIDARestorePrm()
683   PrmNode **prm;  in THIInitial()  local
693   PetscCall(THIDAGetPrm(da, &prm));  in THIInitial()
699 …prm[i + 1][j].b + prm[i + 1][j].h - prm[i - 1][j].b - prm[i - 1][j].h) / (2 * hx), drivingy = thi-…  in THIInitial()
700         x[i][j][k].u = 0. * drivingx * prm[i][j].h * (PetscScalar)k / zm1;  in THIInitial()
701         x[i][j][k].v = 0. * drivingy * prm[i][j].h * (PetscScalar)k / zm1;  in THIInitial()
706   PetscCall(THIDARestorePrm(da, &prm));  in THIInitial()
756   PrmNode **prm;  in THIFunctionLocal()  local
772   PetscCall(THIDAGetPrm(info->da, &prm));  in THIFunctionLocal()
777       QuadExtract(prm, i, j, pn);  in THIFunctionLocal()
849   PetscCall(THIDARestorePrm(info->da, &prm));  in THIFunctionLocal()
969   PrmNode **prm;  in THIJacobianLocal_2D()  local
980   PetscCall(THIDAGetPrm(info->da, &prm));  in THIJacobianLocal_2D()
987       QuadExtract(prm, i, j, pn);  in THIJacobianLocal_2D()
1034   PetscCall(THIDARestorePrm(info->da, &prm));  in THIJacobianLocal_2D()
1047   PrmNode **prm;  local
1060   PetscCall(THIDAGetPrm(info->da, &prm));
1065       QuadExtract(prm, i, j, pn);
1224   PetscCall(THIDARestorePrm(info->da, &prm));