impls/pounders/pounders.c

351     PetscCall(VecGetArrayRead(mfqP->Xhist[mfqP->model_indices[i]], &x));
354     PetscCall(VecRestoreArrayRead(mfqP->Xhist[mfqP->model_indices[i]], &x));
373       PetscCall(VecCopy(mfqP->Xhist[point], mfqP->workxvec));
374       PetscCall(VecAXPY(mfqP->workxvec, -1.0, mfqP->Xhist[mfqP->minindex]));
384     PetscCall(VecGetArrayRead(mfqP->Xhist[point], &x));
387     PetscCall(VecRestoreArrayRead(mfqP->Xhist[point], &x));
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]));
476   PetscCall(pounders_feval(tao, mfqP->Xhist[mfqP->nHist], mfqP->Fhist[mfqP->nHist], &mfqP->Fres[mfqP->nHist]));
540     PetscCall(VecGetArrayRead(mfqP->Xhist[i], &x));
542     PetscCall(VecRestoreArrayRead(mfqP->Xhist[i], &x));
612     PetscCall(VecCopy(tao->solution, mfqP->Xhist[0]));
613     PetscCall(VecAXPY(mfqP->Xhist[0], -1.0, tao->XU));
614     PetscCall(VecMax(mfqP->Xhist[0], NULL, &val));
618     PetscCall(VecCopy(tao->XL, mfqP->Xhist[0]));
619     PetscCall(VecAXPY(mfqP->Xhist[0], -1.0, tao->solution));
620     PetscCall(VecMax(mfqP->Xhist[0], NULL, &val));
625     PetscCall(VecSet(mfqP->Xhist[0], mfqP->delta));
626     PetscCall(VecAXPY(mfqP->Xhist[0], 1.0, tao->solution));
627     PetscCall(VecAXPY(mfqP->Xhist[0], -1.0, tao->XU));
628     PetscCall(VecMax(mfqP->Xhist[0], NULL, &val));
637   PetscCall(VecCopy(tao->solution, mfqP->Xhist[0]));
643   PetscCall(pounders_feval(tao, mfqP->Xhist[0], mfqP->Fhist[0], &mfqP->Fres[0]));
647   PetscCall(VecGetOwnershipRange(mfqP->Xhist[0], &low, &high));
649     PetscCall(VecCopy(mfqP->Xhist[0], mfqP->Xhist[i]));
652       PetscCall(VecGetArray(mfqP->Xhist[i], &x));
654       PetscCall(VecRestoreArray(mfqP->Xhist[i], &x));
657     PetscCall(pounders_feval(tao, mfqP->Xhist[i], mfqP->Fhist[i], &mfqP->Fres[i]));
663   PetscCall(VecCopy(mfqP->Xhist[mfqP->minindex], tao->solution));
678     PetscCall(VecGetArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
680     PetscCall(VecRestoreArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
685       PetscCall(VecGetArray(mfqP->Xhist[i], &x));
687       PetscCall(VecRestoreArray(mfqP->Xhist[i], &x));
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));
712       PetscCall(VecScatterBegin(mfqP->scatterx, mfqP->Xhist[ii], mfqP->localx, INSERT_VALUES, SCATTER_FORWARD));
713       PetscCall(VecScatterEnd(mfqP->scatterx, mfqP->Xhist[ii], mfqP->localx, INSERT_VALUES, SCATTER_FORWARD));
746   PetscCall(VecCopy(mfqP->Xhist[mfqP->minindex], tao->solution));
767     PetscCall(VecDuplicate(tao->solution, &mfqP->Xhist[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]));
803       PetscCall(VecGetArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
805       PetscCall(VecRestoreArrayRead(mfqP->Xhist[mfqP->minindex], &xmint));
851       PetscCall(VecGetArray(mfqP->Xhist[mfqP->model_indices[i]], &x));
853       PetscCall(VecRestoreArray(mfqP->Xhist[mfqP->model_indices[i]], &x));
885     PetscCall(VecCopy(mfqP->Xhist[mfqP->minindex], tao->solution));
936   PetscCall(PetscMalloc1(tao->max_funcs + 100, &mfqP->Xhist));
939     PetscCall(VecDuplicate(tao->solution, &mfqP->Xhist[i]));
1098     PetscCall(VecDestroy(&mfqP->Xhist[i]));
1103   PetscCall(PetscFree(mfqP->Xhist));