343   for (PetscInt i = 0; i < dm->periodic.num_affines; i++) {  in DMCoordAddPeriodicOffsets_Private()
344 …PetscCall(VecScatterBegin(dm->periodic.affine_to_local[i], dm->periodic.affine[i], l, ADD_VALUES, …  in DMCoordAddPeriodicOffsets_Private()
345 …PetscCall(VecScatterEnd(dm->periodic.affine_to_local[i], dm->periodic.affine[i], l, ADD_VALUES, SC…  in DMCoordAddPeriodicOffsets_Private()
561   if (plex->periodic.num_face_sfs > 0 && sfNatural_old) {  in DMPlexCorrectOrientationForIsoperiodic()
594   for (PetscInt f = 0; f < plex->periodic.num_face_sfs; f++) {  in DMPlexCorrectOrientationForIsoperiodic()
595     PetscSF face_sf                   = plex->periodic.face_sfs[f];  in DMPlexCorrectOrientationForIsoperiodic()
596     const PetscScalar (*transform)[4] = (const PetscScalar (*)[4])plex->periodic.transform[f];  in DMPlexCorrectOrientationForIsoperiodic()
1000 …periodic.composed_sf) PetscCall(DMPlexCreateIsoperiodicPointSF_Private(dm, plex->periodic.num_face…  in DMGetIsoperiodicPointSF_Plex()
1001   if (sf) *sf = plex->periodic.composed_sf;  in DMGetIsoperiodicPointSF_Plex()
1012   if (!plex->periodic.face_sfs) PetscFunctionReturn(PETSC_SUCCESS);  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1015   PetscCall(PetscMalloc1(plex->periodic.num_face_sfs, &new_face_sfs));  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1017   for (PetscInt f = 0; f < plex->periodic.num_face_sfs; f++) {  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1023 …PetscCall(PetscSFGetGraph(plex->periodic.face_sfs[f], &old_npoints, &old_nleaf, &old_local, &old_r…  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1048 …PetscCall(PetscSFBcastBegin(plex->periodic.face_sfs[f], MPIU_SF_NODE, rootdata, leafdata, MPI_REPL…  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1049 …PetscCall(PetscSFBcastEnd(plex->periodic.face_sfs[f], MPIU_SF_NODE, rootdata, leafdata, MPI_REPLAC…  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1077   PetscCall(DMPlexSetIsoperiodicFaceSF(dm, plex->periodic.num_face_sfs, new_face_sfs));  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1078 …cCall(DMPlexSetIsoperiodicFaceTransform(dm, plex->periodic.num_face_sfs, (PetscScalar *)plex->peri…  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1079 …for (PetscInt f = 0; f < plex->periodic.num_face_sfs; f++) PetscCall(PetscSFDestroy(&new_face_sfs[…  in DMPlexMigrateIsoperiodicFaceSF_Internal()
1092   if (!plex->periodic.face_sfs) PetscFunctionReturn(PETSC_SUCCESS);  in DMPeriodicCoordinateSetUp_Internal()
1093 …PetscCheck(plex->periodic.periodic_points, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONGSTATE, "Isoperiodic…  in DMPeriodicCoordinateSetUp_Internal()
1096   dm->periodic.num_affines = plex->periodic.num_face_sfs;  in DMPeriodicCoordinateSetUp_Internal()
1097   PetscCall(PetscFree2(dm->periodic.affine_to_local, dm->periodic.affine));  in DMPeriodicCoordinateSetUp_Internal()
1098 …PetscCall(PetscMalloc2(dm->periodic.num_affines, &dm->periodic.affine_to_local, dm->periodic.num_a…  in DMPeriodicCoordinateSetUp_Internal()
1100   for (PetscInt f = 0; f < plex->periodic.num_face_sfs; f++) {  in DMPeriodicCoordinateSetUp_Internal()
1103     IS              is = plex->periodic.periodic_points[f];  in DMPeriodicCoordinateSetUp_Internal()
1147       for (PetscInt j = 0; j < dim; j++) x[i * dim + j] = plex->periodic.transform[f][j][3];  in DMPeriodicCoordinateSetUp_Internal()
1151     dm->periodic.affine_to_local[f] = scatter;  in DMPeriodicCoordinateSetUp_Internal()
1152     dm->periodic.affine[f]          = P;  in DMPeriodicCoordinateSetUp_Internal()
1457 …if (num_face_sfs == plex->periodic.num_face_sfs && (num_face_sfs == 0 || face_sfs == plex->periodi…  in DMPlexSetIsoperiodicFaceSF()
1461   if (plex->periodic.num_face_sfs > 0) {  in DMPlexSetIsoperiodicFaceSF()
1462 …for (PetscInt i = 0; i < plex->periodic.num_face_sfs; i++) PetscCall(PetscSFDestroy(&plex->periodi…  in DMPlexSetIsoperiodicFaceSF()
1463     PetscCall(PetscFree(plex->periodic.face_sfs));  in DMPlexSetIsoperiodicFaceSF()
1466   plex->periodic.num_face_sfs = num_face_sfs;  in DMPlexSetIsoperiodicFaceSF()
1467   PetscCall(PetscCalloc1(num_face_sfs, &plex->periodic.face_sfs));  in DMPlexSetIsoperiodicFaceSF()
1468   for (PetscInt i = 0; i < num_face_sfs; i++) plex->periodic.face_sfs[i] = face_sfs[i];  in DMPlexSetIsoperiodicFaceSF()
1473     if (face_sfs) cdm->periodic.setup = DMPeriodicCoordinateSetUp_Internal;  in DMPlexSetIsoperiodicFaceSF()
1503     if (face_sfs) *face_sfs = plex->periodic.face_sfs;  in DMPlexGetIsoperiodicFaceSF()
1504     if (num_face_sfs) *num_face_sfs = plex->periodic.num_face_sfs;  in DMPlexGetIsoperiodicFaceSF()
1544 …periodic.num_face_sfs, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_OUTOFRANGE, "Number of tran…  in DMPlexSetIsoperiodicFaceTransform()
1546   PetscCall(PetscFree(plex->periodic.transform));  in DMPlexSetIsoperiodicFaceTransform()
1547   PetscCall(PetscMalloc1(n, &plex->periodic.transform));  in DMPlexSetIsoperiodicFaceTransform()
1552         plex->periodic.transform[i][j][k] = t[i * 16 + j * 4 + k];  in DMPlexSetIsoperiodicFaceTransform()