impls/pounders/pounders.c

374       PetscCall(VecAXPY(mfqP->workxvec, -1.0, mfqP->Xhist[mfqP->minindex]));
468   PetscCall(VecAYPX(mfqP->Xhist[mfqP->nHist], mfqP->delta, mfqP->Xhist[mfqP->minindex]));
488   PetscInt     i, j, minindex = 0;
518       minindex = i;
523   if (!addallpoints) PetscCall(addpoint(tao, mfqP, minindex));
644   mfqP->minindex = 0;
659       mfqP->minindex = i;
663   PetscCall(VecCopy(mfqP->Xhist[mfqP->minindex], tao->solution));
664   PetscCall(VecCopy(mfqP->Fhist[mfqP->minindex], tao->ls_res));
671   /* Disp[i] = Xi-xmin, i=1,..,mfqP->minindex-1,mfqP->minindex+1,..,n */
672   /* Fdiff[i] = (Fi-Fmin)', i=1,..,mfqP->minindex-1,mfqP->minindex+1,..,n */
678     PetscCall(VecGetArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
680     PetscCall(VecRestoreArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
681     PetscCall(VecGetArrayRead(mfqP->Fhist[mfqP->minindex], &fmin));
683       if (i == mfqP->minindex) continue;
696     PetscCall(VecRestoreArrayRead(mfqP->Fhist[mfqP->minindex], &fmin));
699     PetscCall(VecScatterBegin(mfqP->scatterx, mfqP->Xhist[mfqP->minindex], mfqP->localxmin, INSERT_VALUES, SCATTER_FORWARD));
700     PetscCall(VecScatterEnd(mfqP->scatterx, mfqP->Xhist[mfqP->minindex], mfqP->localxmin, INSERT_VALUES, SCATTER_FORWARD));
706     PetscCall(VecScatterBegin(mfqP->scatterf, mfqP->Fhist[mfqP->minindex], mfqP->localfmin, INSERT_VALUES, SCATTER_FORWARD));
707     PetscCall(VecScatterEnd(mfqP->scatterf, mfqP->Fhist[mfqP->minindex], mfqP->localfmin, INSERT_VALUES, SCATTER_FORWARD));
710       if (i == mfqP->minindex) continue;
746   PetscCall(VecCopy(mfqP->Xhist[mfqP->minindex], tao->solution));
775     rho = (mfqP->Fres[mfqP->minindex] - mfqP->Fres[mfqP->nHist]) / mdec;
799       mfqP->minindex = mfqP->nHist - 1;
800       minnorm        = mfqP->Fres[mfqP->minindex];
801       PetscCall(VecCopy(mfqP->Fhist[mfqP->minindex], tao->ls_res));
803       PetscCall(VecGetArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
805       PetscCall(VecRestoreArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
848     mfqP->model_indices[0] = mfqP->minindex;
868     PetscCall(VecGetArrayRead(mfqP->Fhist[mfqP->minindex], &fmin));
885     PetscCall(VecCopy(mfqP->Xhist[mfqP->minindex], tao->solution));