46   am->dualres *= am->mu;  in TaoADMMToleranceUpdate()
101   am->muold = am->mu;  in AdaptiveADMMPenaltyUpdate()
103     am->mu = PetscSqrtReal(a_k * b_k);  in AdaptiveADMMPenaltyUpdate()
105     am->mu = a_k;  in AdaptiveADMMPenaltyUpdate()
107     am->mu = b_k;  in AdaptiveADMMPenaltyUpdate()
109   if (am->mu > am->muold) am->mu = am->muold;  in AdaptiveADMMPenaltyUpdate()
110   if (am->mu < am->mumin) am->mu = am->mumin;  in AdaptiveADMMPenaltyUpdate()
189   *f += temp + (am->mu / 2) * temp2;  in SubObjGradUpdate()
193   PetscCall(VecAXPY(g, am->mu, tempJR));  in SubObjGradUpdate()
217   *f += temp + (am->mu / 2) * temp2;  in RegObjGradUpdate()
221   PetscCall(VecAXPY(g, am->mu, tempJR));  in RegObjGradUpdate()
254     if (H && (am->muold != am->mu)) {  in ADMMInternalHessianUpdate()
256         PetscCall(MatAXPY(H, am->mu - am->muold, Constraint, DIFFERENT_NONZERO_PATTERN));  in ADMMInternalHessianUpdate()
258         PetscCall(MatShift(H, am->mu - am->muold));  in ADMMInternalHessianUpdate()
386         PetscCall(MatAXPY(am->subsolverX->hessian, am->mu, am->ATA, DIFFERENT_NONZERO_PATTERN));  in TaoSolve_ADMM()
388         PetscCall(MatShift(am->subsolverX->hessian, am->mu));  in TaoSolve_ADMM()
393         PetscCall(MatAXPY(am->subsolverZ->hessian, am->mu, am->BTB, DIFFERENT_NONZERO_PATTERN));  in TaoSolve_ADMM()
395         PetscCall(MatShift(am->subsolverZ->hessian, am->mu));  in TaoSolve_ADMM()
425       PetscCall(VecWAXPY(am->workJacobianRight, 1 / am->mu, am->y, am->Ax));  in TaoSolve_ADMM()
426 …PetscCall(TaoSoftThreshold(am->workJacobianRight, -am->lambda / am->mu, am->lambda / am->mu, am->s…  in TaoSolve_ADMM()
433     PetscCall(VecWAXPY(am->y, am->mu, am->residual, am->yold));  in TaoSolve_ADMM()
441       am->muold = am->mu;  in TaoSolve_ADMM()
449         PetscCall(VecWAXPY(am->yhatold, -am->mu, am->residual, am->yold));  in TaoSolve_ADMM()
452         am->muold = am->mu;  in TaoSolve_ADMM()
457         PetscCall(VecWAXPY(am->yhat, -am->mu, am->residual, am->yold));  in TaoSolve_ADMM()
464         am->muold = am->mu;  in TaoSolve_ADMM()
512 …tao_admm_spectral_penalty", "Constant for Augmented Lagrangian term.", "", am->mu, &am->mu, NULL));  in TaoSetFromOptions_ADMM()
717   am->mu              = 1.;  in TaoCreate_ADMM()
817 PetscErrorCode TaoADMMSetSpectralPenalty(Tao tao, PetscReal mu)  in TaoADMMSetSpectralPenalty()  argument
822   am->mu = mu;  in TaoADMMSetSpectralPenalty()
841 PetscErrorCode TaoADMMGetSpectralPenalty(Tao tao, PetscReal *mu)  in TaoADMMGetSpectralPenalty()  argument
847   PetscAssertPointer(mu, 2);  in TaoADMMGetSpectralPenalty()
848   *mu = am->mu;  in TaoADMMGetSpectralPenalty()
935 PetscErrorCode TaoADMMSetMinimumSpectralPenalty(Tao tao, PetscReal mu)  in TaoADMMSetMinimumSpectralPenalty()  argument
940   am->mumin = mu;  in TaoADMMSetMinimumSpectralPenalty()