Home
last modified time | relevance | path

Searched refs:lo (Results 1 – 25 of 25) sorted by relevance

/petsc/src/sys/classes/draw/impls/x/
H A Dximage.c47 int lo, mid, hi = PETSC_DRAW_MAXCOLOR; in PetscDrawXiPixelToColor() local
49 for (lo = 0; lo < 8; lo++) in PetscDrawXiPixelToColor()
50 if (pix == cmap[lo]) return lo; in PetscDrawXiPixelToColor()
52 while (hi - lo > 1) { in PetscDrawXiPixelToColor()
53 mid = lo + (hi - lo) / 2; in PetscDrawXiPixelToColor()
55 else lo = mid; in PetscDrawXiPixelToColor()
57 return arg[lo]; in PetscDrawXiPixelToColor()
/petsc/src/sys/utils/
H A Dsorti.c68 #define TwoWayPartition1(X, pivot, t1, lo, hi, l, r) \ argument
70 l = lo; \
101 #define TwoWayPartitionReverse1(X, pivot, t1, lo, hi, l, r) \ argument
103 l = lo; \
118 #define TwoWayPartition2(X, Y, pivot, t1, t2, lo, hi, l, r) \ argument
120 l = lo; \
135 #define TwoWayPartition3(X, Y, Z, pivot, t1, t2, t3, lo, hi, l, r) \ argument
137 l = lo; \
537 PetscInt lo = 0, hi; in PetscFindInt() local
547 while (hi - lo > 1) { in PetscFindInt()
[all …]
H A Dsortd.c208 PetscInt lo = 0, hi; in PetscFindReal() local
218 while (hi - lo > 1) { in PetscFindReal()
219 PetscInt mid = lo + (hi - lo) / 2; in PetscFindReal()
220 …cAssert(t[lo] <= t[mid] && t[mid] <= t[hi - 1], PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Input array… in PetscFindReal()
222 else lo = mid; in PetscFindReal()
224 *loc = (PetscAbsReal(key - t[lo]) < eps) ? lo : -(lo + (key > t[lo]) + 1); in PetscFindReal()
H A Dsortso.c825 PetscInt lo = runstart, hi = ri; in PetscTimSortBuildRun_Private() local
826 for (; lo < hi; ++lo, --hi) COPYSWAPPY(arr + lo * size, arr + hi * size, tarr, size); in PetscTimSortBuildRun_Private()
867 PetscInt lo = runstart, hi = ri; in PetscTimSortBuildRunWithArray_Private() local
868 …for (; lo < hi; ++lo, --hi) COPYSWAPPY2(arr + lo * asize, arr + hi * asize, asize, barr + lo * bsi… in PetscTimSortBuildRunWithArray_Private()
/petsc/src/vec/is/utils/
H A Dpmap.c619 PetscMPIInt lo = 0, hi, t; in PetscLayoutFindOwner() local
626 while (hi - lo > 1) { in PetscLayoutFindOwner()
627 t = lo + (hi - lo) / 2; in PetscLayoutFindOwner()
629 else lo = t; in PetscLayoutFindOwner()
631 *owner = lo; in PetscLayoutFindOwner()
654 PetscMPIInt lo = 0, hi, t; in PetscLayoutFindOwnerIndex() local
660 while (hi - lo > 1) { in PetscLayoutFindOwnerIndex()
661 t = lo + (hi - lo) / 2; in PetscLayoutFindOwnerIndex()
663 else lo = t; in PetscLayoutFindOwnerIndex()
665 if (owner) *owner = lo; in PetscLayoutFindOwnerIndex()
[all …]
/petsc/src/vec/vec/tests/
H A Dex52.c17 PetscInt n, lo, N = 10; in main() local
36 PetscCall(VecGetOwnershipRange(x, &lo, NULL)); in main()
38 for (PetscInt i = 0; i < n; ++i) array[i] = (PetscScalar)(lo + i); in main()
/petsc/src/dm/impls/swarm/tutorials/
H A Dex1.c17 PetscReal lo[3], hi[3], h[3]; in main() local
34 PetscCall(DMGetBoundingBox(dm, lo, hi)); in main()
36 h[i] = (hi[i] - lo[i]) / faces[i]; in main()
37 …ETSC_COMM_SELF, " lo = %g hi = %g n = %" PetscInt_FMT " h = %g\n", (double)lo[i], (double)hi[i], f… in main()
/petsc/src/dm/impls/swarm/tutorials/output/
H A Dex1_bjacobi.out12 lo = -2. hi = 2. n = 4 h = 1.
13 lo = 0. hi = 2. n = 2 h = 1.
H A Dex1_0.out12 lo = -2. hi = 2. n = 4 h = 1.
13 lo = 0. hi = 2. n = 2 h = 1.
/petsc/src/tao/pde_constrained/tutorials/
H A Dhyperbolic.c100 PetscInt lo, hi, hi2, lo2, ksp_old; in main() local
146 PetscCall(VecGetOwnershipRange(user.y, &lo, &hi)); in main()
149 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo, 1, &is_allstate)); in main()
150 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo + lo2, 1, &user.s_is)); in main()
155 PetscCall(VecSetSizes(x, hi - lo + hi2 - lo2, user.n)); in main()
675 PetscInt n, i, j, linear_index, istart, iend, iblock, lo, hi; in HyperbolicInitialize() local
883 PetscCall(VecGetOwnershipRange(user->yi[i], &lo, &hi)); in HyperbolicInitialize()
884 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo, 1, &is_to_yi)); in HyperbolicInitialize()
885 … PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo + i * user->mx * user->mx, 1, &is_from_y)); in HyperbolicInitialize()
911 PetscCall(VecGetOwnershipRange(user->uxi[i], &lo, &hi)); in HyperbolicInitialize()
[all …]
H A Dparabolic.c102 PetscInt lo, hi, hi2, lo2, ksp_old; in main() local
151 PetscCall(VecGetOwnershipRange(user.y, &lo, &hi)); in main()
154 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo, 1, &is_allstate)); in main()
155 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo + lo2, 1, &user.s_is)); in main()
160 PetscCall(VecSetSizes(x, hi - lo + hi2 - lo2, user.n)); in main()
658 PetscInt m, n, i, j, k, linear_index, istart, iend, iblock, lo, hi, lo2, hi2; in ParabolicInitialize() local
975 PetscCall(VecGetOwnershipRange(user->yi[i], &lo, &hi)); in ParabolicInitialize()
976 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo, 1, &is_to_yi)); in ParabolicInitialize()
977 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo2 + istart, 1, &is_from_y)); in ParabolicInitialize()
979 istart = istart + hi - lo; in ParabolicInitialize()
[all …]
H A Delliptic.c548 PetscInt lo, hi, hi2, lo2, ysubnlocal, dsubnlocal; in EllipticInitialize() local
598 PetscCall(VecGetOwnershipRange(user->y, &lo, &hi)); in EllipticInitialize()
601 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo, 1, &is_allstate)); in EllipticInitialize()
602 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo + lo2, 1, &user->s_is)); in EllipticInitialize()
606 PetscCall(VecSetSizes(user->x, hi - lo + hi2 - lo2, user->n)); in EllipticInitialize()
625 PetscCall(VecGetOwnershipRange(user->suby, &lo, &hi)); in EllipticInitialize()
626 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo2 + istart, 1, &is_from_y)); in EllipticInitialize()
628 istart = istart + hi - lo; in EllipticInitialize()
648 PetscCall(VecGetOwnershipRange(user->subd, &lo, &hi)); in EllipticInitialize()
649 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi - lo, lo2 + istart, 1, &is_from_d)); in EllipticInitialize()
[all …]
/petsc/src/ts/tests/
H A Dex30.c556 …PetscReal lo[3] = {-ctx->radius[grid], -ctx->radius[grid], -ctx->radius[grid]}, hi[3] = {ctx… in go() local
563 if (dim == 2) lo[0] = 0; // Landau coordinate (r,z) in go()
570 hp[0] = (hi[0] - lo[0]) / Npi; in go()
571 hp[1] = (hi[1] - lo[1]) / Npj; in go()
572 hp[2] = (hi[2] - lo[2]) / Npk; in go()
574 …, " lo = %14.7e, hi = %14.7e; hp = %14.7e, %14.7e; kT_m = %g; \n", (double)lo[1], (double)hi[1], (… in go()
583 xx_t[grid][tid][pp] = lo[0] + hp[0] / 2.0 + pi * hp[0]; in go()
584 yy_t[grid][tid][pp] = lo[1] + hp[1] / 2.0 + pj * hp[1]; in go()
585 if (dim == 3) zz_t[grid][tid][pp] = lo[2] + hp[2] / 2.0 + pk * hp[2]; in go()
607 … if (dim == 2 && xx_t[grid][tid][pp] < lo[0]) xx_t[grid][tid][pp] = -xx_t[grid][tid][pp]; in go()
[all …]
/petsc/src/dm/dt/dualspace/impls/lagrange/tutorials/
H A Dex2.c16 PetscReal radius = 2, lo[] = {-radius, -radius, -radius}, hi[] = {radius, radius, radius}; in main() local
28 …PetscCall(DMPlexCreateBoxMesh(PETSC_COMM_SELF, dim, PETSC_FALSE, cells, lo, hi, periodicity, PETSC… in main()
/petsc/src/dm/impls/swarm/tests/
H A Dex7.c223 …PetscReal lo[3] = {-5, 0, -5}, hi[3] = {5, 5, 5}, h[3], hp[3], *xx_t[MAX_NUM_THRDS], *yy_t[MAX… in go() local
266 PetscCall(DMGetBoundingBox(dm_t[tid], lo, hi)); in go()
270 h[i] = (hi[i] - lo[i]) / faces[i]; in go()
271 hp[i] = (hi[i] - lo[i]) / Np[i]; in go()
272 vol *= (hi[i] - lo[i]); in go()
273 …t[tid], " lo = %g hi = %g n = %" PetscInt_FMT " h = %g hp = %g\n", (double)lo[i], (double)hi[i], f… in go()
288 xx_t[tid][pp] = lo[0] + hp[0] / 2. + pi * hp[0]; in go()
289 yy_t[tid][pp] = lo[1] + hp[1] / 2. + pj * hp[1]; in go()
/petsc/src/mat/impls/composite/
H A Dmcomposite.c135 PetscInt i, j, k, n, nuniq, lo, hi, mid, *gindices, *buf, *tmp, tot; in MatMult_Composite() local
209 lo = 0; in MatMult_Composite()
212 while (hi - lo > 1) { in MatMult_Composite()
213 mid = lo + (hi - lo) / 2; in MatMult_Composite()
215 else lo = mid; in MatMult_Composite()
217 shell->location[tot + i] = lo; /* gindices[lo] = garray[i] */ in MatMult_Composite()
218 lo++; /* Since garray[i+1] > garray[i], we can safely advance lo */ in MatMult_Composite()
/petsc/src/tao/pde_constrained/impls/lcl/
H A Dlcl.c96 PetscInt lo, hi, nlocalstate, nlocaldesign; in TaoSetup_LCL() local
155 PetscCall(VecGetOwnershipRange(lclP->U, &lo, &hi)); in TaoSetup_LCL()
156 PetscCall(ISCreateStride(((PetscObject)lclP->U)->comm, hi - lo, lo, 1, &is_state)); in TaoSetup_LCL()
157 PetscCall(VecGetOwnershipRange(lclP->V, &lo, &hi)); in TaoSetup_LCL()
164 PetscCall(ISCreateStride(((PetscObject)lclP->V)->comm, hi - lo, lo, 1, &is_design)); in TaoSetup_LCL()
/petsc/src/dm/impls/plex/
H A Dplexsfc.c194 PetscInt lo = 0, hi = n; in ZCodeFind() local
197 while (hi - lo > 1) { in ZCodeFind()
198 PetscInt mid = lo + (hi - lo) / 2; in ZCodeFind()
200 else lo = mid; in ZCodeFind()
202 return key == X[lo] ? lo : -(lo + (key > X[lo]) + 1); in ZCodeFind()
/petsc/src/dm/impls/composite/
H A Dpack.c913 PetscInt subi = indices[i], lo = 0, hi = size, t; in DMCompositeGetISLocalToGlobalMappings() local
921 while (hi - lo > 1) { in DMCompositeGetISLocalToGlobalMappings()
922 t = lo + (hi - lo) / 2; in DMCompositeGetISLocalToGlobalMappings()
924 else lo = t; in DMCompositeGetISLocalToGlobalMappings()
926 idx[i] = subi - suboff[lo] + next->grstarts[lo]; in DMCompositeGetISLocalToGlobalMappings()
/petsc/src/dm/impls/forest/p4est/
H A Dpforest.h658 PetscInt lo, hi, guess; in DMPforestComputeOverlappingRanks() local
660 lo = 0; in DMPforestComputeOverlappingRanks()
674 lo = guess + 1; in DMPforestComputeOverlappingRanks()
676 guess = lo + (hi - lo) / 2; in DMPforestComputeOverlappingRanks()
679 lo = 0; in DMPforestComputeOverlappingRanks()
692 lo = guess + 1; in DMPforestComputeOverlappingRanks()
694 guess = lo + (hi - lo) / 2; in DMPforestComputeOverlappingRanks()
3223 PetscInt lo = 0, hi = num_trees; in DMPforestLabelsInitialize() local
3231 …PetscCheck(guess >= lo && guess < num_trees && lo < hi, PETSC_COMM_SELF, PETSC_ERR_PLIB, "failed b… in DMPforestLabelsInitialize()
3240 lo = guess + 1; in DMPforestLabelsInitialize()
[all …]
/petsc/src/dm/impls/plex/cgns/
H A Dplexcgns2.c921 PetscMPIInt lo = 0, hi, t; in PetscLayoutFindOwnerIndex_Internal() local
932 while (hi - lo > 1) { in PetscLayoutFindOwnerIndex_Internal()
933 t = lo + (hi - lo) / 2; in PetscLayoutFindOwnerIndex_Internal()
935 else lo = t; in PetscLayoutFindOwnerIndex_Internal()
937 if (owner) *owner = lo; in PetscLayoutFindOwnerIndex_Internal()
938 if (lidx) *lidx = idx - map->range[lo]; in PetscLayoutFindOwnerIndex_Internal()
/petsc/lib/petsc/conf/
H A Drules15 @-${RM} ${CLEANFILES} ${TESTS} *.o *.lo *~ \
/petsc/src/ts/utils/dmplexlandau/
H A Dplexland.c597 …PetscReal lo[] = {-perp_radius, -par_radius, -par_radius}, hi[] = {perp_radius, par_radius, p… in LandauDMCreateVMeshes() local
599 if (dim == 2) lo[0] = 0; in LandauDMCreateVMeshes()
601 lo[1] = -perp_radius; in LandauDMCreateVMeshes()
604 …PetscCall(DMPlexCreateBoxMesh(comm_self, dim, PETSC_FALSE, ctx->cells0, lo, hi, periodicity, PETSC… in LandauDMCreateVMeshes()
/petsc/src/dm/impls/plex/hdf5/
H A Dplexhdf5.c2499 PetscInt lo, n, size; in DMPlexTopologyBuildFromLayers_Private() local
2501 PetscCall(PetscSectionGetChart(l->coneSizesSection, &lo, &n)); in DMPlexTopologyBuildFromLayers_Private()
2503 …PetscCheck(lo == 0, comm, PETSC_ERR_PLIB, "starting point should be 0 in coneSizesSection %" Petsc… in DMPlexTopologyBuildFromLayers_Private()
/petsc/src/ksp/pc/impls/bddc/
H A Dbddcprivate.c1962 PetscBool gl, lo; in PCBDDCComputeLocalTopologyInfo() local
2000 lo = (PetscBool)(l2l && corners); in PCBDDCComputeLocalTopologyInfo()
2001 … PetscCallMPI(MPIU_Allreduce(&lo, &gl, 1, MPI_C_BOOL, MPI_LAND, PetscObjectComm((PetscObject)pc))); in PCBDDCComputeLocalTopologyInfo()