impls/morethuente/morethuente.c

366   PetscReal         gamma1, p, q, r, s, sgnd, stpc, stpf, stpq, theta;
390     gamma1     = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s));
391     if (*stp < *stx) gamma1 = -gamma1;
393     p    = (gamma1 - *dx) + theta;
394     q    = ((gamma1 - *dx) + gamma1) + *dp;
413     gamma1     = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s));
414     if (*stp > *stx) gamma1 = -gamma1;
415     p    = (gamma1 - *dp) + theta;
416     q    = ((gamma1 - *dp) + gamma1) + *dx;
440     /* The case gamma1 = 0 only arises if the cubic does not tend
442     gamma1 = s * PetscSqrtScalar(PetscMax(0.0, PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s)));
443     if (*stp > *stx) gamma1 = -gamma1;
444     p = (gamma1 - *dp) + theta;
445     q = (gamma1 + (*dx - *dp)) + gamma1;
447     if (r < 0.0 && gamma1 != 0.0) stpc = *stp + r * (*stx - *stp);
471       gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dy / s) * (*dp / s));
472       if (*stp > *sty) gamma1 = -gamma1;
473       p    = (gamma1 - *dp) + theta;
474       q    = ((gamma1 - *dp) + gamma1) + *dy;