Lines Matching refs:da

33 PetscErrorCode DMDAGetLogicalCoordinate(DM da, PetscScalar x, PetscScalar y, PetscScalar z, PetscIn…  in DMDAGetLogicalCoordinate()  argument
43 …PetscCheck(da->dim != 1, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Cannot get point from 1… in DMDAGetLogicalCoordinate()
44 …PetscCheck(da->dim != 3, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Cannot get point from 3… in DMDAGetLogicalCoordinate()
49 PetscCall(DMGetCoordinateDM(da, &dacoors)); in DMDAGetLogicalCoordinate()
51 PetscCall(DMGetCoordinates(da, &coors)); in DMDAGetLogicalCoordinate()
63 PetscCallMPI(MPIU_Allreduce(&D, &Dv, 1, MPIU_REAL, MPIU_MIN, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
71 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)da), &rank)); in DMDAGetLogicalCoordinate()
74 PetscCallMPI(MPIU_Allreduce(&rank, &root, 1, MPI_INT, MPI_SUM, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
76 PetscCallMPI(MPI_Bcast(X, 1, MPIU_SCALAR, root, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
77 PetscCallMPI(MPI_Bcast(Y, 1, MPIU_SCALAR, root, PetscObjectComm((PetscObject)da))); in DMDAGetLogicalCoordinate()
103 PetscErrorCode DMDAGetRay(DM da, DMDirection dir, PetscInt gp, Vec *newvec, VecScatter *scatter) in DMDAGetRay() argument
106 DM_DA *dd = (DM_DA *)da->data; in DMDAGetRay()
113 …PetscCheck(da->dim != 3, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Cannot get slice from 3… in DMDAGetRay()
114 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)da), &rank)); in DMDAGetRay()
115 PetscCall(DMDAGetAO(da, &ao)); in DMDAGetRay()
117 if (da->dim == 1) { in DMDAGetRay()
129 …PetscCheck(dir != DM_Y, PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "Cannot get Y slice from … in DMDAGetRay()
130 SETERRQ(PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Unknown DMDirection"); in DMDAGetRay()
164 PetscCall(DMGetGlobalVector(da, &vec)); in DMDAGetRay()
166 PetscCall(DMRestoreGlobalVector(da, &vec)); in DMDAGetRay()
199 PetscErrorCode DMDAGetProcessorSubset(DM da, DMDirection dir, PetscInt gp, MPI_Comm *comm) in DMDAGetProcessorSubset() argument
205 DM_DA *dd = (DM_DA *)da->data; in DMDAGetProcessorSubset()
208 PetscValidHeaderSpecificType(da, DM_CLASSID, 1, DMDA); in DMDAGetProcessorSubset()
210 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAGetProcessorSubset()
211 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)da), &size)); in DMDAGetProcessorSubset()
213 …PetscCheck(da->dim >= 3, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "DM_Z invalid… in DMDAGetProcessorSubset()
217 …PetscCheck(da->dim != 1, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "DM_Y invalid… in DMDAGetProcessorSubset()
223 } else SETERRQ(PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Invalid direction"); in DMDAGetProcessorSubset()
226 …cCallMPI(MPI_Allgather(&flag, 1, MPIU_INT, owners, 1, MPIU_INT, PetscObjectComm((PetscObject)da))); in DMDAGetProcessorSubset()
228 …PetscCall(PetscInfo(da, "DMDAGetProcessorSubset: dim=%" PetscInt_FMT ", direction=%d, procs: ", da in DMDAGetProcessorSubset()
233 PetscCall(PetscInfo(da, "%d", i)); /* checkbadSource \n */ in DMDAGetProcessorSubset()
236 PetscCall(PetscInfo(da, "\n")); in DMDAGetProcessorSubset()
237 PetscCallMPI(MPI_Comm_group(PetscObjectComm((PetscObject)da), &group)); in DMDAGetProcessorSubset()
239 PetscCallMPI(MPI_Comm_create(PetscObjectComm((PetscObject)da), subgroup, comm)); in DMDAGetProcessorSubset()
269 PetscErrorCode DMDAGetProcessorSubsets(DM da, DMDirection dir, MPI_Comm *subcomm) in DMDAGetProcessorSubsets() argument
279 PetscValidHeaderSpecificType(da, DM_CLASSID, 1, DMDA); in DMDAGetProcessorSubsets()
280 PetscCall(PetscObjectGetComm((PetscObject)da, &comm)); in DMDAGetProcessorSubsets()
281 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAGetProcessorSubsets()
284 PetscCheck(da->dim >= 3, comm, PETSC_ERR_ARG_OUTOFRANGE, "DM_Z invalid for DMDA dim < 3"); in DMDAGetProcessorSubsets()
287 PetscCheck(da->dim != 1, comm, PETSC_ERR_ARG_OUTOFRANGE, "DM_Y invalid for DMDA dim = 1"); in DMDAGetProcessorSubsets()
295 …PetscCall(PetscInfo(da, "DMDAGetProcessorSubset: dim=%" PetscInt_FMT ", direction=%d, procs: ", da in DMDAGetProcessorSubsets()
299 PetscCall(PetscInfo(da, "%d", p)); /* checkbadSource \n */ in DMDAGetProcessorSubsets()
302 PetscCall(PetscInfo(da, "\n")); in DMDAGetProcessorSubsets()