impls/pounders/pounders.c

3 static PetscErrorCode pounders_h(Tao subtao, Vec v, Mat H, Mat Hpre, PetscCtx ctx)
9 static PetscErrorCode pounders_fg(Tao subtao, Vec x, PetscReal *f, Vec g, PetscCtx ctx)
88     PetscCall(TaoResetStatistics(mfqP->subtao));
89     /* PetscCall(TaoSetTolerances(mfqP->subtao,*gnorm,*gnorm,PETSC_CURRENT)); */
121     PetscCall(TaoSolve(mfqP->subtao));
122     PetscCall(TaoGetSolutionStatus(mfqP->subtao, NULL, qmin, NULL, NULL, NULL, NULL));
1024     PetscCall(TaoCreate(PETSC_COMM_SELF, &mfqP->subtao));
1025     PetscCall(PetscObjectIncrementTabLevel((PetscObject)mfqP->subtao, (PetscObject)tao, 1));
1026     PetscCall(TaoSetType(mfqP->subtao, TAOBNTR));
1027     PetscCall(TaoSetOptionsPrefix(mfqP->subtao, "pounders_subsolver_"));
1028     PetscCall(TaoSetSolution(mfqP->subtao, mfqP->subx));
1029     PetscCall(TaoSetObjectiveAndGradient(mfqP->subtao, NULL, pounders_fg, (void *)mfqP));
1030     PetscCall(TaoSetMaximumIterations(mfqP->subtao, mfqP->gqt_maxits));
1031     PetscCall(TaoSetFromOptions(mfqP->subtao));
1032     PetscCall(TaoGetKSP(mfqP->subtao, &ksp));
1037     PetscCall(TaoSetVariableBounds(mfqP->subtao, mfqP->subxl, mfqP->subxu));
1039     PetscCall(TaoSetHessian(mfqP->subtao, mfqP->subH, mfqP->subH, pounders_h, (void *)mfqP));
1051     PetscCall(TaoDestroy(&mfqP->subtao));
1144       PetscCall(TaoView(mfqP->subtao, viewer));