440   PetscReal tmp = 1.0, ynorm, ynorm2 = 1.0, snorm = 1.0, dk_yk = 1.0, gd;  in TaoBNCGStepDirectionUpdate()  local
450     ynorm2 = ynorm * ynorm;  in TaoBNCGStepDirectionUpdate()
497   } else if (ynorm2 > PETSC_MACHINE_EPSILON) {  in TaoBNCGStepDirectionUpdate()
503           PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->sts, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
520         PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->sts, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
535       ynorm2 = ynorm * ynorm;  in TaoBNCGStepDirectionUpdate()
538 …PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, step * step * dnorm * dnorm, &tau_k, c…  in TaoBNCGStepDirectionUpdate()
555 …    PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, snorm * snorm, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
570       ynorm2 = ynorm * ynorm;  in TaoBNCGStepDirectionUpdate()
574 …    PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, snorm * snorm, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
594         PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->yts, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
624           PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
627           beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk));  in TaoBNCGStepDirectionUpdate()
634           cg->yty = ynorm2;  in TaoBNCGStepDirectionUpdate()
668         PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
671 …beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk) + gd / (dnorm * dnorm)) - step * g…  in TaoBNCGStepDirectionUpdate()
677         cg->yty = ynorm2;  in TaoBNCGStepDirectionUpdate()
714           PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha));  in TaoBNCGStepDirectionUpdate()
715           beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk)) - step * gd / dk_yk;  in TaoBNCGStepDirectionUpdate()
730           cg->yty = ynorm2;  in TaoBNCGStepDirectionUpdate()
773       cg->yty = ynorm2;  in TaoBNCGStepDirectionUpdate()
802       cg->yty = ynorm2;  in TaoBNCGStepDirectionUpdate()
833       cg->yty = ynorm2;  in TaoBNCGStepDirectionUpdate()