impls/pounders/gqt.c

248     anorm = PetscMax(anorm, wa2[j]);
258     pars = PetscMax(pars, -wa1[j]);
259     parl = PetscMax(parl, wa1[j] + wa2[j]);
260     paru = PetscMax(paru, -wa1[j] + wa2[j]);
262   parl = PetscMax(bnorm / delta - parl, pars);
263   parl = PetscMax(0.0, parl);
264   paru = PetscMax(0.0, bnorm / delta + paru);
269   par = PetscMax(par, parl);
273   paru = PetscMax(paru, (1.0 + rtol) * parl);
280     if (par <= pars && paru > 0) par = PetscMax(p001, PetscSqrtScalar(parl / paru)) * paru;
316       pars = PetscMax(pars, par - rznorm * rznorm);
355         parl = PetscMax(parl, par);
386       pars = PetscMax(pars, par + parc);
392       paru = PetscMax(paru, (1.0 + rtol) * pars);
396     parl = PetscMax(parl, pars);
425     par = PetscMax(parl, par + parc);