50   PetscReal          fmin, ftrial, prered, actred, kappa, sigma, beta;  in TaoSolve_NTR()  local
135           PetscCall(VecDot(tao->gradient, tao->stepdirection, &prered));  in TaoSolve_NTR()
137           prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm));  in TaoSolve_NTR()
139           if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) {  in TaoSolve_NTR()
142             kappa = actred / prered;  in TaoSolve_NTR()
145 … * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust + (1.0 - tr->theta_i) * prered - actred);  in TaoSolve_NTR()
146 … * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust - (1.0 + tr->theta_i) * prered + actred);  in TaoSolve_NTR()
289         PetscCall(KSPCGGetObjFcn(tao->ksp, &prered));  in TaoSolve_NTR()
290         if (prered >= 0.0) {  in TaoSolve_NTR()
306             prered = -prered;  in TaoSolve_NTR()
307 …          if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) {  in TaoSolve_NTR()
310               kappa = actred / prered;  in TaoSolve_NTR()
338         PetscCall(KSPCGGetObjFcn(tao->ksp, &prered));  in TaoSolve_NTR()
339         if (prered >= 0.0) {  in TaoSolve_NTR()
353             prered = -prered;  in TaoSolve_NTR()
354 …          if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) {  in TaoSolve_NTR()
357               kappa = actred / prered;  in TaoSolve_NTR()
360             tau_1   = tr->theta * beta / (tr->theta * beta - (1.0 - tr->theta) * prered + actred);  in TaoSolve_NTR()
361             tau_2   = tr->theta * beta / (tr->theta * beta + (1.0 + tr->theta) * prered - actred);  in TaoSolve_NTR()