| /petsc/src/snes/interface/noise/ |
| H A D | snesdnest.c | 126 d__3 = PetscMax(est1, est2); in SNESNoise_dnest_() 129 d__1 = PetscMax(d__3, est3) - est4; in SNESNoise_dnest_() 131 err2 = PetscMax(d__1, d__2); in SNESNoise_dnest_() 151 f_max = PetscMax(d__1, d__2); in SNESNoise_dnest_() 168 scale = PetscMax(d__2, d__3); in SNESNoise_dnest_() 193 if (PetscMin(d__1, d__2) < 0. && PetscMax(d__3, d__4) > 0.) dsgn[j - 1] = TRUE_; in SNESNoise_dnest_() 215 d__1 = PetscMax(eps[4], eps[5]); in SNESNoise_dnest_() 216 emax = PetscMax(d__1, eps[6]); in SNESNoise_dnest_() 235 d__1 = PetscMax(eps[3], eps[4]); in SNESNoise_dnest_() 236 emax = PetscMax(d__1, eps[5]); in SNESNoise_dnest_()
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | gqt.c | 248 anorm = PetscMax(anorm, wa2[j]); in gqt() 258 pars = PetscMax(pars, -wa1[j]); in gqt() 259 parl = PetscMax(parl, wa1[j] + wa2[j]); in gqt() 260 paru = PetscMax(paru, -wa1[j] + wa2[j]); in gqt() 262 parl = PetscMax(bnorm / delta - parl, pars); in gqt() 263 parl = PetscMax(0.0, parl); in gqt() 264 paru = PetscMax(0.0, bnorm / delta + paru); in gqt() 269 par = PetscMax(par, parl); in gqt() 273 paru = PetscMax(paru, (1.0 + rtol) * parl); in gqt() 280 if (par <= pars && paru > 0) par = PetscMax(p001, PetscSqrtScalar(parl / paru)) * paru; in gqt() [all …]
|
| /petsc/src/snes/tutorials/network/power/ |
| H A D | power.h | 21 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 45 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 67 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 98 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 128 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 141 } PFDATA PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar)));
|
| /petsc/src/mat/utils/ |
| H A D | bandwidth.c | 36 lbw[0] = PetscMax(lbw[0], r - cols[0]); in MatComputeBandwidth() 37 lbw[1] = PetscMax(lbw[1], cols[ncols - 1] - r); in MatComputeBandwidth() 42 *bw = 2 * PetscMax(gbw[0], gbw[1]) + 1; in MatComputeBandwidth()
|
| /petsc/src/tao/unconstrained/impls/nls/ |
| H A D | nls.c | 91 tao->trust = PetscMax(tao->trust, nlsP->min_radius); in TaoSolve_NLS() 168 tau_max = PetscMax(tau_1, tau_2); in TaoSolve_NLS() 172 max_radius = PetscMax(max_radius, tao->trust); in TaoSolve_NLS() 187 max_radius = PetscMax(max_radius, tao->trust); in TaoSolve_NLS() 231 tao->trust = PetscMax(tao->trust, max_radius); in TaoSolve_NLS() 234 tao->trust = PetscMax(tao->trust, nlsP->min_radius); in TaoSolve_NLS() 290 tao->trust = PetscMax(tao->trust, nlsP->min_radius); in TaoSolve_NLS() 298 tao->trust = PetscMax(tao->trust, nlsP->min_radius); in TaoSolve_NLS() 350 pert = PetscMin(nlsP->imax, PetscMax(nlsP->imin, nlsP->imfac * gnorm)); in TaoSolve_NLS() 353 pert = PetscMax(pert, -e_min); in TaoSolve_NLS() [all …]
|
| /petsc/src/dm/impls/patch/ |
| H A D | patch.c | 63 blower.i = PetscMax(lower.i - halo, 0); in DMPatchZoom() 65 blower.j = PetscMax(lower.j - halo, 0); in DMPatchZoom() 67 blower.k = PetscMax(lower.k - halo, 0); in DMPatchZoom() 85 sxr = PetscMax(sxb, lower.i - blower.i); in DMPatchZoom() 86 syr = PetscMax(syb, lower.j - blower.j); in DMPatchZoom() 87 szr = PetscMax(szb, lower.k - blower.k); in DMPatchZoom() 91 …PetscCall(PetscMalloc2(dof * rM * rN * PetscMax(rP, 1), &localPoints, dof * rM * rN * PetscMax(rP,… in DMPatchZoom() 208 M = PetscMax(M, 1); in DMPatchSolve() 209 l = PetscMax(l, 1); in DMPatchSolve() 210 N = PetscMax(N, 1); in DMPatchSolve() [all …]
|
| /petsc/src/tao/linesearch/impls/morethuente/ |
| H A D | morethuente.c | 117 ls->stepmax = PetscMax(stx, sty); in TaoLineSearchApply_MT() 124 ls->step = PetscMax(ls->step, ls->stepmin); in TaoLineSearchApply_MT() 372 …PetscCheck(!mtP->bracket || (*stp > PetscMin(*stx, *sty) && *stp < PetscMax(*stx, *sty)), PETSC_CO… in Tao_mcstep() 388 s = PetscMax(PetscAbsReal(theta), PetscAbsReal(*dx)); in Tao_mcstep() 389 s = PetscMax(s, PetscAbsReal(*dp)); in Tao_mcstep() 411 s = PetscMax(PetscAbsReal(theta), PetscAbsReal(*dx)); in Tao_mcstep() 412 s = PetscMax(s, PetscAbsReal(*dp)); in Tao_mcstep() 437 s = PetscMax(PetscAbsReal(theta), PetscAbsReal(*dx)); in Tao_mcstep() 438 s = PetscMax(s, PetscAbsReal(*dp)); in Tao_mcstep() 442 …gamma1 = s * PetscSqrtScalar(PetscMax(0.0, PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s))… in Tao_mcstep() [all …]
|
| /petsc/include/petsc/private/ |
| H A D | dmnetworkimpl.h | 34 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 39 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 44 } DMNetworkComponent PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 130 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar)));
|
| /petsc/src/tao/unconstrained/impls/ntr/ |
| H A D | ntr.c | 72 tao->trust = PetscMax(tao->trust, tr->min_radius); in TaoSolve_NTR() 148 tau_max = PetscMax(tau_1, tau_2); in TaoSolve_NTR() 152 max_radius = PetscMax(max_radius, tao->trust); in TaoSolve_NTR() 163 max_radius = PetscMax(max_radius, tao->trust); in TaoSolve_NTR() 207 tao->trust = PetscMax(tao->trust, max_radius); in TaoSolve_NTR() 210 tao->trust = PetscMax(tao->trust, tr->min_radius); in TaoSolve_NTR() 257 tao->trust = PetscMax(tao->trust, tr->min_radius); in TaoSolve_NTR() 265 tao->trust = PetscMax(tao->trust, tr->min_radius); in TaoSolve_NTR() 327 tao->trust = PetscMax(tr->alpha4 * norm_d, tao->trust); in TaoSolve_NTR() 330 tao->trust = PetscMax(tr->alpha5 * norm_d, tao->trust); in TaoSolve_NTR() [all …]
|
| /petsc/src/tao/unconstrained/impls/cg/ |
| H A D | taocg.c | 42 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG() 46 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG() 69 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG() 73 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG() 96 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG() 100 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG() 164 beta = PetscMax((gnorm2 - ginner) / gnorm2_old, 0.0); in TaoSolve_CG() 190 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG()
|
| /petsc/src/tao/util/ |
| H A D | tao_util.c | 301 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDFischer() 315 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDFischer() 332 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDFischer() 347 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDFischer() 422 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDSFischer() 430 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDSFischer() 442 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDSFischer() 450 ai = PetscMax(PETSC_MACHINE_EPSILON, ai); in MatDSFischer() 475 return PetscMax(0, PetscRealPart(in) - ub) - PetscMax(0, -PetscRealPart(in) - PetscAbsReal(lb)); in ST_InternalPN() 480 …return PetscMax(0, PetscRealPart(in) + PetscAbsReal(ub)) - PetscMax(0, -PetscRealPart(in) - PetscA… in ST_InternalNN() [all …]
|
| /petsc/src/snes/tutorials/network/water/ |
| H A D | water.h | 53 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 92 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 106 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar)));
|
| /petsc/src/tao/unconstrained/impls/ntl/ |
| H A D | ntl.c | 66 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL() 142 tau_max = PetscMax(tau_1, tau_2); in TaoSolve_NTL() 146 max_radius = PetscMax(max_radius, tao->trust); in TaoSolve_NTL() 161 max_radius = PetscMax(max_radius, tao->trust); in TaoSolve_NTL() 205 tao->trust = PetscMax(tao->trust, max_radius); in TaoSolve_NTL() 208 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL() 256 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL() 264 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL() 333 tao->trust = PetscMax(tl->alpha4 * norm_d, tao->trust); in TaoSolve_NTL() 336 tao->trust = PetscMax(tl->alpha5 * norm_d, tao->trust); in TaoSolve_NTL() [all …]
|
| /petsc/src/mat/graphops/color/impls/minpack/ |
| H A D | dsm.c | 105 *mingrp = PetscMax(i__2, i__3); in MINPACKdsm() 118 *mingrp = PetscMax(*mingrp, maxclq); in MINPACKdsm() 129 *mingrp = PetscMax(*mingrp, maxclq); in MINPACKdsm()
|
| /petsc/src/dm/impls/da/ |
| H A D | fdda.c | 850 pstart = (bx == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -i)); in DMCreateMatrix_DA_2d_MPISELL() 856 lstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_2d_MPISELL() 888 pstart = (bx == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -i)); in DMCreateMatrix_DA_2d_MPISELL() 894 lstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_2d_MPISELL() 952 istart = (bx == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -i)); in DMCreateMatrix_DA_3d_MPISELL() 955 jstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_3d_MPISELL() 958 kstart = (bz == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -k)); in DMCreateMatrix_DA_3d_MPISELL() 994 istart = (bx == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -i)); in DMCreateMatrix_DA_3d_MPISELL() 997 jstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_3d_MPISELL() 1000 kstart = (bz == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -k)); in DMCreateMatrix_DA_3d_MPISELL() [all …]
|
| H A D | gr2.c | 331 …nt64)vec_size, PetscMin((PetscInt64)max_chunk_size, PetscMax((PetscInt64)avg_local_vec_size, Petsc… in VecGetHDF5ChunkSize() 333 …t)PetscMax(1, chunkDims[0]) * PetscMax(1, chunkDims[1]) * PetscMax(1, chunkDims[2]) * PetscMax(1, … in VecGetHDF5ChunkSize() 372 …t)PetscMax(1, chunkDims[0]) * PetscMax(1, chunkDims[1]) * PetscMax(1, chunkDims[2]) * PetscMax(1, … in VecGetHDF5ChunkSize() 403 …t)PetscMax(1, chunkDims[0]) * PetscMax(1, chunkDims[1]) * PetscMax(1, chunkDims[2]) * PetscMax(1, … in VecGetHDF5ChunkSize()
|
| H A D | grvtk.c | 99 maxnnodes = PetscMax(maxnnodes, nnodes); in DMDAVTKWriteAll_VTS() 116 maxbs = PetscMax(maxbs, bs); in DMDAVTKWriteAll_VTS() 148 …PetscCall(PetscMalloc2(maxnnodes * PetscMax(3, maxbs), &array, maxnnodes * PetscMax(3, maxbs), &ar… in DMDAVTKWriteAll_VTS() 313 maxnnodes = PetscMax(maxnnodes, nnodes); in DMDAVTKWriteAll_VTR() 333 maxbs = PetscMax(maxbs, bs); in DMDAVTKWriteAll_VTR() 364 …PetscCall(PetscMalloc2(PetscMax(maxnnodes * maxbs, info.xm + info.ym + info.zm), &array, PetscMax(… in DMDAVTKWriteAll_VTR()
|
| /petsc/src/ts/tutorials/network/ |
| H A D | wash.h | 43 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar))); 69 } PETSC_ATTRIBUTEALIGNED(PetscMax(sizeof(double), sizeof(PetscScalar)));
|
| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 17 oldest = PetscMax(0, ldb->k - ldb->sigma_hist); in SymBroydenRescaleUpdateScalar() 80 stDs = PetscMax(stDs, ldb->tol); in DiagonalUpdate() 133 oldest = PetscMax(0, ldb->k - ldb->sigma_hist); in SymBroydenRescaleUpdateDiagonal() 179 PetscInt start = PetscMax(oldest, lmvm->k - ldb->sigma_hist); in SymBroydenRescaleUpdateDiagonal() 257 ys_sum = PetscMax(ldb->tol, ys_sum); in SymBroydenRescaleUpdateDiagonal() 272 PetscReal denom = PetscMax(-b - sqrt_disc, ldb->tol); in SymBroydenRescaleUpdateDiagonal() 276 PetscReal denom = PetscMax(2 * a, ldb->tol); in SymBroydenRescaleUpdateDiagonal() 282 PetscReal denom = PetscMax(2 * a, ldb->tol); in SymBroydenRescaleUpdateDiagonal() 286 PetscReal denom = PetscMax(-b + sqrt_disc, ldb->tol); in SymBroydenRescaleUpdateDiagonal() 325 PetscInt new_oldest = PetscMax(0, next - ldb->sigma_hist); in SymBroydenRescaleUpdate() [all …]
|
| /petsc/src/ksp/ksp/impls/cg/pipelcg/ |
| H A D | pipelcg.c | 4 #define offset(j) PetscMax(((j) - (2 * l)), 0) 51 PetscCall(VecDuplicateVecs(plcg->p, PetscMax(3, l + 1), &plcg->Z)); in KSPSetUp_PIPELCG() 68 PetscCall(VecDestroyVecs(PetscMax(3, l + 1), &plcg->Z)); in KSPReset_PIPELCG() 168 temp = Z[PetscMax(l, 2)]; in KSPSolve_InnerLoop_PIPELCG() 169 for (i = PetscMax(l, 2); i > 0; --i) Z[i] = Z[i - 1]; in KSPSolve_InnerLoop_PIPELCG() 184 for (j = 0; j <= PetscMax(l, 2); ++j) PetscCall(VecScale(Z[j], 1.0 / beta)); in KSPSolve_InnerLoop_PIPELCG() 192 …for (j = PetscMax(0, it - 3 * l + 1); j <= it - 2 * l; j++) G(j, it - l + 1) = G(it - 2 * l + 1, j… in KSPSolve_InnerLoop_PIPELCG() 198 for (j = PetscMax(it - 3 * l + 1, 0); j <= it - l + 1; ++j) G(j, it - l + 1) *= invbeta2; in KSPSolve_InnerLoop_PIPELCG() 201 for (j = PetscMax(it - 2 * l + 2, 0); j <= it - l; ++j) { in KSPSolve_InnerLoop_PIPELCG() 203 …for (k = PetscMax(it - 3 * l + 1, 0); k <= j - 1; ++k) sum_dummy = sum_dummy + G(k, j) * G(k, it -… in KSPSolve_InnerLoop_PIPELCG() [all …]
|
| /petsc/src/tao/bound/impls/bnk/ |
| H A D | bnk.c | 190 tau_max = PetscMax(tau_1, tau_2); in TaoBNKInitialize() 194 max_radius = PetscMax(max_radius, tao->trust); in TaoBNKInitialize() 205 max_radius = PetscMax(max_radius, tao->trust); in TaoBNKInitialize() 261 tao->trust = PetscMax(tao->trust, max_radius); in TaoBNKInitialize() 264 tao->trust = PetscMax(tao->trust, bnk->min_radius); in TaoBNKInitialize() 459 tao->trust = PetscMax(tao->trust, bnk->min_radius); in TaoBNKComputeStep() 467 tao->trust = PetscMax(tao->trust, bnk->min_radius); in TaoBNKComputeStep() 574 bnk->pert = PetscMin(bnk->imax, PetscMax(bnk->imin, bnk->imfac * bnk->gnorm)); in TaoBNKSafeguardStep() 578 bnk->pert = PetscMax(bnk->pert, -e_min); in TaoBNKSafeguardStep() 582 bnk->pert = PetscMin(bnk->pmax, PetscMax(bnk->pgfac * bnk->pert, bnk->pmgfac * bnk->gnorm)); in TaoBNKSafeguardStep() [all …]
|
| /petsc/src/vec/vec/utils/tagger/impls/ |
| H A D | relative.c | 31 bxs[j].max = PetscMax(bxs[j].max, vArray[k]); in VecTaggerComputeBoxes_Relative() 34 …bxs[j].max = PetscCMPLX(PetscMax(PetscRealPart(bxs[j].max), PetscRealPart(vArray[k])), PetscMax(Pe… in VecTaggerComputeBoxes_Relative()
|
| H A D | cdf.c | 30 minCDF = PetscMax(0., bxs->min); in VecTaggerComputeBox_CDF_SortedArray() 35 boxes->max = cArray[PetscMax(minInd, maxInd - 1)]; in VecTaggerComputeBox_CDF_SortedArray() 137 B[i].max = PetscMax(A[i].max, B[i].max); in VecTaggerCDFStatsReduce() 196 stats[2].max = PetscMax(stats[2].max, val); in VecTaggerComputeBox_CDF_SortedArray_Iterative() 218 newBounds[i][j][0] = PetscMax(bounds[i][0], bounds[i][1]); in VecTaggerComputeBox_CDF_SortedArray_Iterative() 233 stats[i][section].max = PetscMax(stats[i][section].max, val); in VecTaggerComputeBox_CDF_SortedArray_Iterative() 238 newBounds[i][section][1] = PetscMax(newBounds[i][section][0], thisInd + 1); in VecTaggerComputeBox_CDF_SortedArray_Iterative() 262 maxDiff = PetscMax(maxDiff, diff); in VecTaggerComputeBox_CDF_SortedArray_Iterative()
|
| /petsc/src/sys/classes/draw/interface/ |
| H A D | dtext.c | 138 mlen = PetscMax(mlen, len); in PetscDrawStringBoxed() 152 draw->boundbox_xr = PetscMax(draw->boundbox_xr, right); in PetscDrawStringBoxed() 154 draw->boundbox_yr = PetscMax(draw->boundbox_yr, top); in PetscDrawStringBoxed()
|
| /petsc/src/mat/impls/aij/seq/lusol/ |
| H A D | lusol.c | 236 nnz = PetscMax(lusol->nnz, (int)(lusol->elbowroom * nz)); in MatLUFactorNumeric_LUSOL() 237 nnz = PetscMax(nnz, 5 * n); in MatLUFactorNumeric_LUSOL() 245 nnz = PetscMax(nnz, (int)(lusol->luroom * (lusol->luparm[22] + lusol->luparm[23]))); in MatLUFactorNumeric_LUSOL() 343 lusol->elbowroom = PetscMax(lusol->elbowroom, info->fill); in MatLUFactorSymbolic_LUSOL() 344 nnz = PetscMax((int)(lusol->elbowroom * nz), 5 * n); in MatLUFactorSymbolic_LUSOL()
|