2555   const PetscReal ds2 = 0.5 * dis;  in snapToCylinder()  local
2557   if ((PetscAbsScalar(u[0]) <= ds2) && (PetscAbsScalar(u[1]) <= ds2)) {  in snapToCylinder()
2569       dc = PetscAbsReal(ds2 / sinp);  in snapToCylinder()
2571       xc = ds2 * x / PetscAbsReal(y);  in snapToCylinder()
2572       yc = ds2 * PetscSignReal(y);  in snapToCylinder()
2574       dc = PetscAbsReal(ds2 / cosp);  in snapToCylinder()
2576       xc = ds2 * PetscSignReal(x);  in snapToCylinder()
2577       yc = ds2 * y / PetscAbsReal(x);  in snapToCylinder()
2805     const PetscReal ds2 = dis / 2.0;  in DMPlexCreateHexCylinderMesh_Internal()  local
2825       coords[0 * dim + 0]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2826       coords[0 * dim + 1]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2828       coords[1 * dim + 0]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2829       coords[1 * dim + 1]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2831       coords[2 * dim + 0]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2832       coords[2 * dim + 1]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2834       coords[3 * dim + 0]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2835       coords[3 * dim + 1]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2837       coords[4 * dim + 0]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2838       coords[4 * dim + 1]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2840       coords[5 * dim + 0]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2841       coords[5 * dim + 1]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2843       coords[6 * dim + 0]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2844       coords[6 * dim + 1]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2846       coords[7 * dim + 0]  = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2847       coords[7 * dim + 1]  = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2874         /* 15 31 19 */ coords[16 * dim + 0] = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2875         coords[16 * dim + 1]                = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2877         /* 16 32 22 */ coords[17 * dim + 0] = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2878         coords[17 * dim + 1]                = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2880         /* 17 33 21 */ coords[18 * dim + 0] = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2881         coords[18 * dim + 1]                = ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2883         /* 18 34 20 */ coords[19 * dim + 0] = -ds2;  in DMPlexCreateHexCylinderMesh_Internal()
2884         coords[19 * dim + 1]                = ds2;  in DMPlexCreateHexCylinderMesh_Internal()