298   PetscReal   lambda;  in TaoSolve_BMRM()  local
313   lambda = bmrm->lambda;  in TaoSolve_BMRM()
372         df.Q[i][tao->niter] = df.Q[tao->niter][i] = reg / lambda;  in TaoSolve_BMRM()
387         PetscCall(VecAXPY(bmrm->local_w, -df.x[i] / lambda, pgrad->V));  in TaoSolve_BMRM()
392       reg = 0.5 * lambda * reg * reg;  in TaoSolve_BMRM()
460 …PetscCall(PetscOptionsReal("-tao_bmrm_lambda", "regulariser weight", "", 100, &bmrm->lambda, NULL)…  in TaoSetFromOptions_BMRM()
499   bmrm->lambda = 1.0;  in TaoCreate_BMRM()
671 static PetscReal phi(PetscReal *x, PetscInt n, PetscReal lambda, PetscReal *a, PetscReal b, PetscRe…  in phi()  argument
677     x[i] = -c[i] + lambda * a[i];  in phi()
695   PetscReal lambda, lambdal, lambdau, dlambda, lambda_new;  in project()  local
701   lambda    = 0;  in project()
715   r = phi(x, n, lambda, a, b, c, l, u);  in project()
726     lambdal = lambda;  in project()
728     lambda  = lambda + dlambda;  in project()
729     r       = phi(x, n, lambda, a, b, c, l, u);  in project()
731       lambdal = lambda;  in project()
735       lambda  = lambda + dlambda;  in project()
737       r       = phi(x, n, lambda, a, b, c, l, u);  in project()
739     lambdau = lambda;  in project()
742     lambdau = lambda;  in project()
744     lambda  = lambda - dlambda;  in project()
745     r       = phi(x, n, lambda, a, b, c, l, u);  in project()
747       lambdau = lambda;  in project()
751       lambda  = lambda - dlambda;  in project()
753       r       = phi(x, n, lambda, a, b, c, l, u);  in project()
755     lambdal = lambda;  in project()
766   lambda  = lambdau - dlambda;  in project()
767   r       = phi(x, n, lambda, a, b, c, l, u);  in project()
769 …while (PetscAbsReal(r) > TOL_R && dlambda > TOL_LAM * (1.0 + PetscAbsReal(lambda)) && innerIter < …  in project()
773         lambdau = lambda;  in project()
777         lambda  = lambdau - dlambda;  in project()
781         dlambda    = (lambdau - lambda) / s;  in project()
782         lambda_new = 0.75 * lambdal + 0.25 * lambda;  in project()
783         if (lambda_new < (lambda - dlambda)) lambda_new = lambda - dlambda;  in project()
784         lambdau = lambda;  in project()
786         lambda  = lambda_new;  in project()
787         s       = (lambdau - lambdal) / (lambdau - lambda);  in project()
791         lambdal = lambda;  in project()
795         lambda  = lambdau - dlambda;  in project()
799         dlambda    = (lambda - lambdal) / s;  in project()
800         lambda_new = 0.75 * lambdau + 0.25 * lambda;  in project()
801         if (lambda_new > (lambda + dlambda)) lambda_new = lambda + dlambda;  in project()
802         lambdal = lambda;  in project()
804         lambda  = lambda_new;  in project()
805         s       = (lambdau - lambdal) / (lambdau - lambda);  in project()
808     r = phi(x, n, lambda, a, b, c, l, u);  in project()
811   *lam_ext = lambda;  in project()