126       tao->trust = tao->trust0;  in TaoBNKInitialize()
132       tao->trust = tao->trust0;  in TaoBNKInitialize()
148           PetscCall(VecAXPY(tao->solution, -tao->trust / bnk->gnorm, tao->gradient));  in TaoBNKInitialize()
162               sigma = -tao->trust / bnk->gnorm;  in TaoBNKInitialize()
179 …        prered = tao->trust * (bnk->gnorm - 0.5 * tao->trust * prered / (bnk->gnorm * bnk->gnorm));  in TaoBNKInitialize()
187 …tau_1   = bnk->theta_i * bnk->gnorm * tao->trust / (bnk->theta_i * bnk->gnorm * tao->trust + (1.0 …  in TaoBNKInitialize()
188 …tau_2   = bnk->theta_i * bnk->gnorm * tao->trust / (bnk->theta_i * bnk->gnorm * tao->trust - (1.0 …  in TaoBNKInitialize()
194               max_radius = PetscMax(max_radius, tao->trust);  in TaoBNKInitialize()
205               max_radius = PetscMax(max_radius, tao->trust);  in TaoBNKInitialize()
231           tao->trust = tau * 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()
265       tao->trust = PetscMin(tao->trust, bnk->max_radius);  in TaoBNKInitialize()
270       tao->trust = 0.0;  in TaoBNKInitialize()
444   PetscCall(KSPCGSetRadius(tao->ksp, tao->trust));  in TaoBNKComputeStep()
453     if (0.0 == tao->trust) {  in TaoBNKComputeStep()
456         tao->trust = bnk->dnorm;  in TaoBNKComputeStep()
459         tao->trust = PetscMax(tao->trust, bnk->min_radius);  in TaoBNKComputeStep()
460         tao->trust = PetscMin(tao->trust, bnk->max_radius);  in TaoBNKComputeStep()
464         tao->trust = tao->trust0;  in TaoBNKComputeStep()
467         tao->trust = PetscMax(tao->trust, bnk->min_radius);  in TaoBNKComputeStep()
468         tao->trust = PetscMin(tao->trust, bnk->max_radius);  in TaoBNKComputeStep()
470         PetscCall(KSPCGSetRadius(tao->ksp, tao->trust));  in TaoBNKComputeStep()
809         tao->trust = bnk->omega1 * PetscMin(bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
812         tao->trust = bnk->omega2 * PetscMin(bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
816           tao->trust = bnk->omega3 * PetscMin(bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
818           tao->trust = PetscMax(bnk->omega3 * bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
822         tao->trust = PetscMax(bnk->omega4 * bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
825         tao->trust = PetscMax(bnk->omega5 * bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
829       tao->trust = bnk->omega1 * PetscMin(bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
839         tao->trust = bnk->alpha1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
842           tao->trust = bnk->alpha1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
852             tao->trust = bnk->alpha1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
857             if (bnk->dnorm == tao->trust) {  in TaoBNKUpdateTrustRadius()
860                 tao->trust = bnk->alpha2 * tao->trust;  in TaoBNKUpdateTrustRadius()
863                 tao->trust = bnk->alpha3 * tao->trust;  in TaoBNKUpdateTrustRadius()
866                 tao->trust = bnk->alpha4 * tao->trust;  in TaoBNKUpdateTrustRadius()
869                 tao->trust = bnk->alpha5 * tao->trust;  in TaoBNKUpdateTrustRadius()
877       tao->trust = bnk->alpha1 * PetscMin(bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
887         tao->trust = bnk->gamma1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
890           tao->trust = bnk->gamma1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
908               tao->trust = PetscMax(tao->trust, bnk->gamma3 * bnk->dnorm);  in TaoBNKUpdateTrustRadius()
910               tao->trust = PetscMax(tao->trust, bnk->gamma4 * bnk->dnorm);  in TaoBNKUpdateTrustRadius()
912               tao->trust = PetscMax(tao->trust, tau_max * bnk->dnorm);  in TaoBNKUpdateTrustRadius()
918               tao->trust = bnk->gamma2 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
920               tao->trust = PetscMax(tao->trust, bnk->gamma3 * bnk->dnorm);  in TaoBNKUpdateTrustRadius()
922               tao->trust = tau_max * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
924               tao->trust = PetscMax(tao->trust, tau_max * bnk->dnorm);  in TaoBNKUpdateTrustRadius()
929               tao->trust = bnk->gamma2 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
931               tao->trust = bnk->gamma1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
933               tao->trust = bnk->gamma1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
935               tao->trust = tau_1 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
937               tao->trust = tau_2 * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
939               tao->trust = tau_max * PetscMin(tao->trust, bnk->dnorm);  in TaoBNKUpdateTrustRadius()
946       tao->trust = bnk->gamma1 * PetscMin(bnk->dnorm, tao->trust);  in TaoBNKUpdateTrustRadius()
951   tao->trust = PetscMin(tao->trust, bnk->max_radius);  in TaoBNKUpdateTrustRadius()
952   tao->trust = PetscMax(tao->trust, bnk->min_radius);  in TaoBNKUpdateTrustRadius()