impls/pounders/pounders.c

359   point              = mfqP->nHist - 1;
464   PetscCall(VecDuplicate(mfqP->Xhist[0], &mfqP->Xhist[mfqP->nHist]));
465   PetscCall(VecSetValues(mfqP->Xhist[mfqP->nHist], mfqP->n, mfqP->indices, &mfqP->Q_tmp[index * mfqP->npmax], INSERT_VALUES));
466   PetscCall(VecAssemblyBegin(mfqP->Xhist[mfqP->nHist]));
467   PetscCall(VecAssemblyEnd(mfqP->Xhist[mfqP->nHist]));
468   PetscCall(VecAYPX(mfqP->Xhist[mfqP->nHist], mfqP->delta, mfqP->Xhist[mfqP->minindex]));
471   if (tao->XU && tao->XL) PetscCall(VecMedian(mfqP->Xhist[mfqP->nHist], tao->XL, tao->XU, mfqP->Xhist[mfqP->nHist]));
474   PetscCall(VecDuplicate(mfqP->Fhist[0], &mfqP->Fhist[mfqP->nHist]));
476   PetscCall(pounders_feval(tao, mfqP->Xhist[mfqP->nHist], mfqP->Fhist[mfqP->nHist], &mfqP->Fres[mfqP->nHist]));
479   mfqP->model_indices[mfqP->nmodelpoints] = mfqP->nHist;
481   mfqP->nHist++;
539   for (i = mfqP->nHist - 1; i >= 0; i--) {
753   mfqP->nHist        = mfqP->n + 1;
767     PetscCall(VecDuplicate(tao->solution, &mfqP->Xhist[mfqP->nHist]));
768     PetscCall(VecDuplicate(tao->ls_res, &mfqP->Fhist[mfqP->nHist]));
769     PetscCall(VecSetValues(mfqP->Xhist[mfqP->nHist], mfqP->n, mfqP->indices, mfqP->work, INSERT_VALUES));
770     PetscCall(VecAssemblyBegin(mfqP->Xhist[mfqP->nHist]));
771     PetscCall(VecAssemblyEnd(mfqP->Xhist[mfqP->nHist]));
773     PetscCall(pounders_feval(tao, mfqP->Xhist[mfqP->nHist], mfqP->Fhist[mfqP->nHist], &mfqP->Fres[mfqP->nHist]));
775     rho = (mfqP->Fres[mfqP->minindex] - mfqP->Fres[mfqP->nHist]) / mdec;
776     mfqP->nHist++;
799       mfqP->minindex = mfqP->nHist - 1;
944   mfqP->nHist = 0;
1097   for (i = 0; i < mfqP->nHist; i++) {