264   PetscInt     i, j, k, num, np = mfqP->nmodelpoints;  in getquadpounders()
360   mfqP->nmodelpoints = mfqP->n + 1;  in morepoints()
361   while (mfqP->nmodelpoints < mfqP->npmax && point >= 0) {  in morepoints()
385     mfqP->M[(mfqP->n + 1) * mfqP->nmodelpoints] = 1.0;  in morepoints()
386 …for (j = 0; j < mfqP->n; j++) mfqP->M[j + 1 + ((mfqP->n + 1) * mfqP->nmodelpoints)] = (x[j] - mfqP…  in morepoints()
388 …&mfqP->M[1 + (mfqP->n + 1) * mfqP->nmodelpoints], mfqP->n, &mfqP->N[mfqP->n * (mfqP->n + 1) / 2 * …  in morepoints()
395     PetscCall(PetscBLASIntCast(mfqP->nmodelpoints + 1, &blasnp));  in morepoints()
416     PetscCall(PetscBLASIntCast(mfqP->nmodelpoints - mfqP->n, &blasint));  in morepoints()
424       mfqP->model_indices[mfqP->nmodelpoints] = point;  in morepoints()
428       mfqP->nmodelpoints++;  in morepoints()
429       PetscCall(PetscBLASIntCast(mfqP->nmodelpoints, &blasnp));  in morepoints()
437   PetscCall(PetscBLASIntCast(mfqP->nmodelpoints, &blasnp));  in morepoints()
451   if (mfqP->nmodelpoints == mfqP->n + 1) {  in morepoints()
479   mfqP->model_indices[mfqP->nmodelpoints] = mfqP->nHist;  in addpoint()
480   mfqP->nmodelpoints++;  in addpoint()
496   PetscCall(PetscBLASIntCast(mfqP->nmodelpoints, &blask));  in modelimprove()
508   for (i = mfqP->nmodelpoints; i < mfqP->n; i++) {  in modelimprove()
517     if (i == mfqP->nmodelpoints || mfqP->work[i] < minvalue) {  in modelimprove()
546       PetscCall(PetscBLASIntCast(PetscMax(mfqP->n - mfqP->nmodelpoints, 0), &blasj));  in affpoints()
549         PetscCall(PetscBLASIntCast(mfqP->nmodelpoints, &blask));  in affpoints()
553       PetscCallBLAS("BLASnrm2", proj = BLASnrm2_(&blasj, &mfqP->work2[mfqP->nmodelpoints], &ione));  in affpoints()
556         mfqP->model_indices[mfqP->nmodelpoints] = i;  in affpoints()
557         mfqP->nmodelpoints++;  in affpoints()
558 …BLAScopy_(&blasn, mfqP->work, &ione, &mfqP->Q_tmp[mfqP->npmax * (mfqP->nmodelpoints - 1)], &ione));  in affpoints()
559         PetscCall(PetscBLASIntCast(mfqP->npmax * (mfqP->nmodelpoints), &blask));  in affpoints()
561         PetscCall(PetscBLASIntCast(mfqP->nmodelpoints, &blask));  in affpoints()
567       if (mfqP->nmodelpoints == mfqP->n) break;  in affpoints()
754   mfqP->nmodelpoints = mfqP->n + 1;  in TaoSolve_POUNDERS()
811       mfqP->nmodelpoints = 0;  in TaoSolve_POUNDERS()
813       if (mfqP->nmodelpoints < mfqP->n) {  in TaoSolve_POUNDERS()
832     mfqP->nmodelpoints = 0;  in TaoSolve_POUNDERS()
835     if (mfqP->nmodelpoints == mfqP->n) {  in TaoSolve_POUNDERS()
841       if (mfqP->n > mfqP->nmodelpoints) {  in TaoSolve_POUNDERS()
843         PetscCall(modelimprove(tao, mfqP, mfqP->n - mfqP->nmodelpoints));  in TaoSolve_POUNDERS()
846     for (i = mfqP->nmodelpoints; i > 0; i--) mfqP->model_indices[i] = mfqP->model_indices[i - 1];  in TaoSolve_POUNDERS()
847     mfqP->nmodelpoints++;  in TaoSolve_POUNDERS()
850     for (i = 0; i < mfqP->nmodelpoints; i++) {  in TaoSolve_POUNDERS()
866 …o, "Get Quad, size: %" PetscInt_FMT ", points: %" PetscInt_FMT "\n", mfqP->n, mfqP->nmodelpoints));  in TaoSolve_POUNDERS()
896     if (mfqP->nmodelpoints == mfqP->last_nmodelpoints) {  in TaoSolve_POUNDERS()
901     for (i = 0; i < mfqP->nmodelpoints; i++) {  in TaoSolve_POUNDERS()
911     mfqP->last_nmodelpoints = mfqP->nmodelpoints;  in TaoSolve_POUNDERS()
1140 …PetscCall(PetscViewerASCIIPrintf(viewer, "model points: %" PetscInt_FMT "\n", mfqP->nmodelpoints));  in TaoView_POUNDERS()