3 static PetscErrorCode TaoSetUp_SSILS(Tao tao)  in TaoSetUp_SSILS()  argument
5   TAO_SSLS *ssls = (TAO_SSLS *)tao->data;  in TaoSetUp_SSILS()
8   PetscCall(VecDuplicate(tao->solution, &tao->gradient));  in TaoSetUp_SSILS()
9   PetscCall(VecDuplicate(tao->solution, &tao->stepdirection));  in TaoSetUp_SSILS()
10   PetscCall(VecDuplicate(tao->solution, &ssls->ff));  in TaoSetUp_SSILS()
11   PetscCall(VecDuplicate(tao->solution, &ssls->dpsi));  in TaoSetUp_SSILS()
12   PetscCall(VecDuplicate(tao->solution, &ssls->da));  in TaoSetUp_SSILS()
13   PetscCall(VecDuplicate(tao->solution, &ssls->db));  in TaoSetUp_SSILS()
14   PetscCall(VecDuplicate(tao->solution, &ssls->t1));  in TaoSetUp_SSILS()
15   PetscCall(VecDuplicate(tao->solution, &ssls->t2));  in TaoSetUp_SSILS()
19 static PetscErrorCode TaoDestroy_SSILS(Tao tao)  in TaoDestroy_SSILS()  argument
21   TAO_SSLS *ssls = (TAO_SSLS *)tao->data;  in TaoDestroy_SSILS()
30   PetscCall(KSPDestroy(&tao->ksp));  in TaoDestroy_SSILS()
31   PetscCall(PetscFree(tao->data));  in TaoDestroy_SSILS()
35 static PetscErrorCode TaoSolve_SSILS(Tao tao)  in TaoSolve_SSILS()  argument
37   TAO_SSLS                    *ssls = (TAO_SSLS *)tao->data;  in TaoSolve_SSILS()
48   PetscCall(TaoComputeVariableBounds(tao));  in TaoSolve_SSILS()
49   PetscCall(VecMedian(tao->XL, tao->solution, tao->XU, tao->solution));  in TaoSolve_SSILS()
50 …etscCall(TaoLineSearchSetObjectiveAndGradientRoutine(tao->linesearch, Tao_SSLS_FunctionGradient, t…  in TaoSolve_SSILS()
51   PetscCall(TaoLineSearchSetObjectiveRoutine(tao->linesearch, Tao_SSLS_Function, tao));  in TaoSolve_SSILS()
55 …PetscCall(TaoLineSearchComputeObjectiveAndGradient(tao->linesearch, tao->solution, &psi, ssls->dps…  in TaoSolve_SSILS()
58   tao->reason = TAO_CONTINUE_ITERATING;  in TaoSolve_SSILS()
60 …PetscCall(PetscInfo(tao, "iter: %" PetscInt_FMT ", merit: %g, ndpsi: %g\n", tao->niter, (double)ss…  in TaoSolve_SSILS()
62     PetscCall(TaoLogConvergenceHistory(tao, ssls->merit, ndpsi, 0.0, tao->ksp_its));  in TaoSolve_SSILS()
63     PetscCall(TaoMonitor(tao, tao->niter, ssls->merit, ndpsi, 0.0, t));  in TaoSolve_SSILS()
64     PetscUseTypeMethod(tao, convergencetest, tao->cnvP);  in TaoSolve_SSILS()
65     if (tao->reason != TAO_CONTINUE_ITERATING) break;  in TaoSolve_SSILS()
68     PetscTryTypeMethod(tao, update, tao->niter, tao->user_update);  in TaoSolve_SSILS()
69     tao->niter++;  in TaoSolve_SSILS()
73     PetscCall(KSPSetOperators(tao->ksp, tao->jacobian, tao->jacobian_pre));  in TaoSolve_SSILS()
74     PetscCall(KSPSolve(tao->ksp, ssls->ff, tao->stepdirection));  in TaoSolve_SSILS()
75     PetscCall(KSPGetIterationNumber(tao->ksp, &tao->ksp_its));  in TaoSolve_SSILS()
76     tao->ksp_tot_its += tao->ksp_its;  in TaoSolve_SSILS()
77     PetscCall(VecNorm(tao->stepdirection, NORM_2, &normd));  in TaoSolve_SSILS()
78     PetscCall(VecDot(tao->stepdirection, ssls->dpsi, &innerd));  in TaoSolve_SSILS()
82       PetscCall(PetscInfo(tao, "newton direction not descent\n"));  in TaoSolve_SSILS()
83       PetscCall(VecCopy(ssls->dpsi, tao->stepdirection));  in TaoSolve_SSILS()
84       PetscCall(VecDot(tao->stepdirection, ssls->dpsi, &innerd));  in TaoSolve_SSILS()
87     PetscCall(VecScale(tao->stepdirection, -1.0));  in TaoSolve_SSILS()
90     PetscCall(TaoLineSearchSetInitialStepLength(tao->linesearch, 1.0));  in TaoSolve_SSILS()
91 …PetscCall(TaoLineSearchApply(tao->linesearch, tao->solution, &psi, ssls->dpsi, tao->stepdirection,…  in TaoSolve_SSILS()
107 PETSC_EXTERN PetscErrorCode TaoCreate_SSILS(Tao tao)  in TaoCreate_SSILS()  argument
114   tao->data                = (void *)ssls;  in TaoCreate_SSILS()
115   tao->ops->solve          = TaoSolve_SSILS;  in TaoCreate_SSILS()
116   tao->ops->setup          = TaoSetUp_SSILS;  in TaoCreate_SSILS()
117   tao->ops->view           = TaoView_SSLS;  in TaoCreate_SSILS()
118   tao->ops->setfromoptions = TaoSetFromOptions_SSLS;  in TaoCreate_SSILS()
119   tao->ops->destroy        = TaoDestroy_SSILS;  in TaoCreate_SSILS()
124   PetscCall(TaoLineSearchCreate(((PetscObject)tao)->comm, &tao->linesearch));  in TaoCreate_SSILS()
125   PetscCall(PetscObjectIncrementTabLevel((PetscObject)tao->linesearch, (PetscObject)tao, 1));  in TaoCreate_SSILS()
126   PetscCall(TaoLineSearchSetType(tao->linesearch, armijo_type));  in TaoCreate_SSILS()
127   PetscCall(TaoLineSearchSetOptionsPrefix(tao->linesearch, tao->hdr.prefix));  in TaoCreate_SSILS()
128   PetscCall(TaoLineSearchSetFromOptions(tao->linesearch));  in TaoCreate_SSILS()
130   PetscCall(KSPCreate(((PetscObject)tao)->comm, &tao->ksp));  in TaoCreate_SSILS()
131   PetscCall(PetscObjectIncrementTabLevel((PetscObject)tao->ksp, (PetscObject)tao, 1));  in TaoCreate_SSILS()
132   PetscCall(KSPSetOptionsPrefix(tao->ksp, tao->hdr.prefix));  in TaoCreate_SSILS()
135   PetscCall(TaoParametersInitialize(tao));  in TaoCreate_SSILS()
136   PetscObjectParameterSetDefault(tao, max_it, 2000);  in TaoCreate_SSILS()
137   PetscObjectParameterSetDefault(tao, max_funcs, 4000);  in TaoCreate_SSILS()
138   PetscObjectParameterSetDefault(tao, gttol, 0);  in TaoCreate_SSILS()
139   PetscObjectParameterSetDefault(tao, grtol, 0);  in TaoCreate_SSILS()
140   PetscObjectParameterSetDefault(tao, gatol, PetscDefined(USE_REAL_SINGLE) ? 1.0e-6 : 1.0e-16);  in TaoCreate_SSILS()
141   PetscObjectParameterSetDefault(tao, fmin, PetscDefined(USE_REAL_SINGLE) ? 1.0e-4 : 1.0e-8);  in TaoCreate_SSILS()