impls/ntr/ntr.c

51   PetscReal          tau, tau_1, tau_2, tau_max, tau_min, max_radius;
145           tau_1   = tr->theta_i * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust + (1.0 - tr->theta_i) * prered - actred);
147           tau_min = PetscMin(tau_1, tau_2);
148           tau_max = PetscMax(tau_1, tau_2);
180             } else if ((tau_1 >= tr->gamma1_i) && (tau_1 < 1.0) && ((tau_2 < tr->gamma1_i) || (tau_2 >= 1.0))) {
181               tau = tau_1;
182             } else if ((tau_2 >= tr->gamma1_i) && (tau_2 < 1.0) && ((tau_1 < tr->gamma1_i) || (tau_2 >= 1.0))) {
360             tau_1   = tr->theta * beta / (tr->theta * beta - (1.0 - tr->theta) * prered + actred);
362             tau_min = PetscMin(tau_1, tau_2);
363             tau_max = PetscMax(tau_1, tau_2);
396               } else if ((tau_1 >= tr->gamma1) && (tau_1 < 1.0) && ((tau_2 < tr->gamma1) || (tau_2 >= 1.0))) {
397                 tao->trust = tau_1 * PetscMin(tao->trust, norm_d);
398               } else if ((tau_2 >= tr->gamma1) && (tau_2 < 1.0) && ((tau_1 < tr->gamma1) || (tau_2 >= 1.0))) {