10 …l *dx, PetscReal *sty, PetscReal *fy, PetscReal *dy, PetscReal *stp, PetscReal *fp, PetscReal *dp);
363 …al *dx, PetscReal *sty, PetscReal *fy, PetscReal *dy, PetscReal *stp, PetscReal *fp, PetscReal *dp)  in Tao_mcstep()  argument
377   sgnd = *dp * (*dx / PetscAbsReal(*dx));  in Tao_mcstep()
387     theta      = 3 * (*fx - *fp) / (*stp - *stx) + *dx + *dp;  in Tao_mcstep()
389     s          = PetscMax(s, PetscAbsReal(*dp));  in Tao_mcstep()
390     gamma1     = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s));  in Tao_mcstep()
394     q    = ((gamma1 - *dx) + gamma1) + *dp;  in Tao_mcstep()
410     theta      = 3 * (*fx - *fp) / (*stp - *stx) + *dx + *dp;  in Tao_mcstep()
412     s          = PetscMax(s, PetscAbsReal(*dp));  in Tao_mcstep()
413     gamma1     = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s));  in Tao_mcstep()
415     p    = (gamma1 - *dp) + theta;  in Tao_mcstep()
416     q    = ((gamma1 - *dp) + gamma1) + *dx;  in Tao_mcstep()
419     stpq = *stp + (*dp / (*dp - *dx)) * (*stx - *stp);  in Tao_mcstep()
424   } else if (PetscAbsReal(*dp) < PetscAbsReal(*dx)) {  in Tao_mcstep()
436     theta      = 3 * (*fx - *fp) / (*stp - *stx) + *dx + *dp;  in Tao_mcstep()
438     s          = PetscMax(s, PetscAbsReal(*dp));  in Tao_mcstep()
442 …gamma1 = s * PetscSqrtScalar(PetscMax(0.0, PetscPowScalar(theta / s, 2.0) - (*dx / s) * (*dp / s))…  in Tao_mcstep()
444     p = (gamma1 - *dp) + theta;  in Tao_mcstep()
445     q = (gamma1 + (*dx - *dp)) + gamma1;  in Tao_mcstep()
450     stpq = *stp + (*dp / (*dp - *dx)) * (*stx - *stp);  in Tao_mcstep()
468       theta  = 3 * (*fp - *fy) / (*sty - *stp) + *dy + *dp;  in Tao_mcstep()
470       s      = PetscMax(s, PetscAbsReal(*dp));  in Tao_mcstep()
471       gamma1 = s * PetscSqrtScalar(PetscPowScalar(theta / s, 2.0) - (*dy / s) * (*dp / s));  in Tao_mcstep()
473       p    = (gamma1 - *dp) + theta;  in Tao_mcstep()
474       q    = ((gamma1 - *dp) + gamma1) + *dy;  in Tao_mcstep()
491     *dy  = *dp;  in Tao_mcstep()
500     *dx  = *dp;  in Tao_mcstep()