Lines Matching refs:nmodelpoints

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