impls/pounders/pounders.c

333   /* Assumes mfqP->model_indices[0]  is minimum index
334    Finishes adding points to mfqP->model_indices (up to npmax)
351     PetscCall(VecGetArrayRead(mfqP->Xhist[mfqP->model_indices[i]], &x));
354     PetscCall(VecRestoreArrayRead(mfqP->Xhist[mfqP->model_indices[i]], &x));
365       if (point == mfqP->model_indices[j]) {
424       mfqP->model_indices[mfqP->nmodelpoints] = point;
479   mfqP->model_indices[mfqP->nmodelpoints] = mfqP->nHist;
556         mfqP->model_indices[mfqP->nmodelpoints] = i;
693       mfqP->model_indices[ii++] = i;
724       mfqP->model_indices[ii++] = i;
846     for (i = mfqP->nmodelpoints; i > 0; i--) mfqP->model_indices[i] = mfqP->model_indices[i - 1];
848     mfqP->model_indices[0] = mfqP->minindex;
851       PetscCall(VecGetArray(mfqP->Xhist[mfqP->model_indices[i]], &x));
853       PetscCall(VecRestoreArray(mfqP->Xhist[mfqP->model_indices[i]], &x));
854       PetscCall(VecGetArray(mfqP->Fhist[mfqP->model_indices[i]], &f));
862       PetscCall(VecRestoreArray(mfqP->Fhist[mfqP->model_indices[i]], &f));
903         if (mfqP->model_indices[i] == mfqP->last_model_indices[i]) {
909       mfqP->last_model_indices[i] = mfqP->model_indices[i];
977   PetscCall(PetscMalloc1(mfqP->npmax, &mfqP->model_indices));
1089   PetscCall(PetscFree(mfqP->model_indices));