| /petsc/src/sys/utils/ |
| H A D | sorti.c | 68 #define TwoWayPartition1(X, pivot, t1, lo, hi, l, r) \ argument 71 r = hi; \ 101 #define TwoWayPartitionReverse1(X, pivot, t1, lo, hi, l, r) \ argument 104 r = hi; \ 118 #define TwoWayPartition2(X, Y, pivot, t1, t2, lo, hi, l, r) \ argument 121 r = hi; \ 135 #define TwoWayPartition3(X, Y, Z, pivot, t1, t2, t3, lo, hi, l, r) \ argument 138 r = hi; \ 155 PetscCount i, j, p, l, r, hi = n - 1; \ 167 p = MEDIAN(X, hi); \ [all …]
|
| H A D | sortd.c | 208 PetscInt lo = 0, hi; in PetscFindReal() local 217 PetscCall(PetscIntCast(n, &hi)); in PetscFindReal() 218 while (hi - lo > 1) { in PetscFindReal() 219 PetscInt mid = lo + (hi - lo) / 2; in PetscFindReal() 220 …id] <= t[hi - 1], PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Input array was not sorted: (%g, %g, %g)"… in PetscFindReal() 221 if (key < t[mid]) hi = mid; in PetscFindReal()
|
| H A D | sortso.c | 825 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 D | pmap.c | 619 PetscMPIInt lo = 0, hi, t; in PetscLayoutFindOwner() local 625 hi = map->size; in PetscLayoutFindOwner() 626 while (hi - lo > 1) { in PetscLayoutFindOwner() 627 t = lo + (hi - lo) / 2; in PetscLayoutFindOwner() 628 if (idx < map->range[t]) hi = t; in PetscLayoutFindOwner() 654 PetscMPIInt lo = 0, hi, t; in PetscLayoutFindOwnerIndex() local 659 hi = map->size; in PetscLayoutFindOwnerIndex() 660 while (hi - lo > 1) { in PetscLayoutFindOwnerIndex() 661 t = lo + (hi - lo) / 2; in PetscLayoutFindOwnerIndex() 662 if (idx < map->range[t]) hi = t; in PetscLayoutFindOwnerIndex()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | seqhashmat.h | 4 PetscHashIter hi; in MatCopyHashToXAIJ_Seq_Hash() local 39 PetscHashIterBegin(a->ht, hi); in MatCopyHashToXAIJ_Seq_Hash() 40 while (!PetscHashIterAtEnd(a->ht, hi)) { in MatCopyHashToXAIJ_Seq_Hash() 41 PetscHashIterGetKey(a->ht, hi, key); in MatCopyHashToXAIJ_Seq_Hash() 42 PetscHashIterGetVal(a->ht, hi, value); in MatCopyHashToXAIJ_Seq_Hash() 45 PetscHashIterNext(a->ht, hi); in MatCopyHashToXAIJ_Seq_Hash()
|
| /petsc/src/sys/classes/draw/impls/x/ |
| H A D | ximage.c | 47 int lo, mid, hi = PETSC_DRAW_MAXCOLOR; in PetscDrawXiPixelToColor() local 52 while (hi - lo > 1) { in PetscDrawXiPixelToColor() 53 mid = lo + (hi - lo) / 2; in PetscDrawXiPixelToColor() 54 if (pix < cmap[arg[mid]]) hi = mid; in PetscDrawXiPixelToColor()
|
| /petsc/src/dm/impls/swarm/tutorials/ |
| H A D | ex1.c | 17 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 … " lo = %g hi = %g n = %" PetscInt_FMT " h = %g\n", (double)lo[i], (double)hi[i], faces[i], (doubl… in main()
|
| /petsc/src/dm/impls/swarm/tutorials/output/ |
| H A D | ex1_bjacobi.out | 12 lo = -2. hi = 2. n = 4 h = 1. 13 lo = 0. hi = 2. n = 2 h = 1.
|
| H A D | ex1_0.out | 12 lo = -2. hi = 2. n = 4 h = 1. 13 lo = 0. hi = 2. n = 2 h = 1.
|
| /petsc/src/ksp/pc/impls/patch/ |
| H A D | pcpatch.c | 740 PetscHashIter hi; in PCPatchCompleteCellPatch() local 755 PetscHashIterBegin(ht, hi); in PCPatchCompleteCellPatch() 756 while (!PetscHashIterAtEnd(ht, hi)) { in PCPatchCompleteCellPatch() 757 PetscHashIterGetKey(ht, hi, point); in PCPatchCompleteCellPatch() 758 PetscHashIterNext(ht, hi); in PCPatchCompleteCellPatch() 840 PetscHashIter hi; in PCPatchGetPointDofs() local 868 PetscHashIterBegin(pts, hi); in PCPatchGetPointDofs() 869 while (!PetscHashIterAtEnd(pts, hi)) { in PCPatchGetPointDofs() 870 PetscHashIterGetKey(pts, hi, p); in PCPatchGetPointDofs() 871 PetscHashIterNext(pts, hi); in PCPatchGetPointDofs() [all …]
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/tutorials/ |
| H A D | ex2.c | 16 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/ts/tests/ |
| H A D | ex30.c | 556 … lo[3] = {-ctx->radius[grid], -ctx->radius[grid], -ctx->radius[grid]}, hi[3] = {ctx->radius[gr… in go() local 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 …, hi = %14.7e; hp = %14.7e, %14.7e; kT_m = %g; \n", (double)lo[1], (double)hi[1], (double)hp[0], (… in go() 588 …e && PetscSqrtReal(PetscSqr(xx_t[grid][tid][pp]) + PetscSqr(yy_t[grid][tid][pp])) > 0.92 * hi[0]) { in go() 612 … while (zz_t[grid][tid][pp] >= hi[2] || zz_t[grid][tid][pp] <= lo[2]) zz_t[grid][tid][pp] *= .9; in go() 614 … while (xx_t[grid][tid][pp] >= hi[0] || xx_t[grid][tid][pp] <= lo[0]) xx_t[grid][tid][pp] *= .9; in go() 615 … while (yy_t[grid][tid][pp] >= hi[1] || yy_t[grid][tid][pp] <= lo[1]) yy_t[grid][tid][pp] *= .9; in go() 618 …al(PetscSqr(xx_t[grid][tid][pp]) + PetscSqr(yy_t[grid][tid][pp])) > 0.92 * hi[0]) { // safety fact… in go() [all …]
|
| /petsc/src/mat/impls/preallocator/ |
| H A D | matpreallocator.c | 153 PetscHashIter hi; in MatPreallocatorPreallocate_Preallocator() local 171 PetscHashIterBegin(p->ht, hi); in MatPreallocatorPreallocate_Preallocator() 172 while (!PetscHashIterAtEnd(p->ht, hi)) { in MatPreallocatorPreallocate_Preallocator() 173 PetscHashIterGetKey(p->ht, hi, key); in MatPreallocatorPreallocate_Preallocator() 177 PetscHashIterNext(p->ht, hi); in MatPreallocatorPreallocate_Preallocator()
|
| /petsc/src/tao/pde_constrained/tutorials/ |
| H A D | hyperbolic.c | 100 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() 153 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi2 - lo2, hi + lo2, 1, &user.d_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() [all …]
|
| H A D | parabolic.c | 102 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() 158 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi2 - lo2, hi + lo2, 1, &user.d_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() [all …]
|
| H A D | elliptic.c | 548 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() 604 PetscCall(ISCreateStride(PETSC_COMM_SELF, hi2 - lo2, hi + lo2, 1, &user->d_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() [all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex7.c | 223 …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 …%g hi = %g n = %" PetscInt_FMT " h = %g hp = %g\n", (double)lo[i], (double)hi[i], faces[i], (doubl… in go()
|
| /petsc/src/mat/impls/adj/mpi/ |
| H A D | mpiadj.c | 509 PetscHashIter hi; in MatAssemblyEnd_MPIAdj() local 536 PetscHashIterBegin(ht, hi); in MatAssemblyEnd_MPIAdj() 537 for (; !PetscHashIterAtEnd(ht, hi);) { in MatAssemblyEnd_MPIAdj() 538 PetscHashIterGetKey(ht, hi, key); in MatAssemblyEnd_MPIAdj() 540 PetscHashIterNext(ht, hi); in MatAssemblyEnd_MPIAdj() 546 PetscHashIterBegin(ht, hi); in MatAssemblyEnd_MPIAdj() 547 for (; !PetscHashIterAtEnd(ht, hi);) { in MatAssemblyEnd_MPIAdj() 548 PetscHashIterGetKey(ht, hi, key); in MatAssemblyEnd_MPIAdj() 550 PetscHashIterNext(ht, hi); in MatAssemblyEnd_MPIAdj()
|
| /petsc/src/mat/impls/composite/ |
| H A D | mcomposite.c | 135 PetscInt i, j, k, n, nuniq, lo, hi, mid, *gindices, *buf, *tmp, tot; in MatMult_Composite() local 211 hi = nuniq; in MatMult_Composite() 212 while (hi - lo > 1) { in MatMult_Composite() 213 mid = lo + (hi - lo) / 2; in MatMult_Composite() 214 if (garray[i] < gindices[mid]) hi = mid; in MatMult_Composite()
|
| /petsc/src/tao/pde_constrained/impls/lcl/ |
| H A D | lcl.c | 96 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/composite/ |
| H A D | pack.c | 913 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() 923 if (suboff[t] > subi) hi = t; in DMCompositeGetISLocalToGlobalMappings()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexsfc.c | 194 PetscInt lo = 0, hi = n; in ZCodeFind() local 197 while (hi - lo > 1) { in ZCodeFind() 198 PetscInt mid = lo + (hi - lo) / 2; in ZCodeFind() 199 if (key < X[mid]) hi = mid; in ZCodeFind()
|
| /petsc/src/dm/impls/forest/p4est/ |
| H A D | pforest.h | 658 PetscInt lo, hi, guess; in DMPforestComputeOverlappingRanks() local 661 hi = size; in DMPforestComputeOverlappingRanks() 672 hi = guess; in DMPforestComputeOverlappingRanks() 676 guess = lo + (hi - lo) / 2; in DMPforestComputeOverlappingRanks() 680 hi = size; in DMPforestComputeOverlappingRanks() 690 hi = guess; 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() 3234 hi = guess; in DMPforestLabelsInitialize() [all …]
|
| /petsc/src/mat/impls/hypre/ |
| H A D | mhypre.c | 975 HYPRE_Int *hi, *hj; local 977 hi = hypre_CSRMatrixI(hdiag); 979 PetscCall(PetscFree2(hi, hj)); 981 hi = hypre_CSRMatrixI(hoffd); 983 PetscCall(PetscFree2(hi, hj));
|
| /petsc/src/dm/impls/plex/cgns/ |
| H A D | plexcgns2.c | 921 PetscMPIInt lo = 0, hi, t; in PetscLayoutFindOwnerIndex_Internal() local 931 hi = map->size; in PetscLayoutFindOwnerIndex_Internal() 932 while (hi - lo > 1) { in PetscLayoutFindOwnerIndex_Internal() 933 t = lo + (hi - lo) / 2; in PetscLayoutFindOwnerIndex_Internal() 934 if (idx < map->range[t]) hi = t; in PetscLayoutFindOwnerIndex_Internal()
|