impls/pounders/pounders.c

71     PetscCall(VecSetValues(mfqP->subb, mfqP->n, mfqP->indices, mfqP->Gres, INSERT_VALUES));
141     PetscCall(gqt(mfqP->n, mfqP->Hres, mfqP->n, mfqP->Gres, 1.0, mfqP->gqt_rtol, atol, mfqP->gqt_maxits, gnorm, qmin, mfqP->Xsubproblem, &info, &its, mfqP->work, mfqP->work2, mfqP->work3));
158   for (i = 0; i < mfqP->n; i++) mfqP->Gres[i] = 0;
161   /* Compute Gres= sum_ij[wij * (cjgi + cigj)] */
167       PetscCallBLAS("BLASaxpy", BLASaxpy_(&blasn, &factor, &mfqP->Fdiff[blasn * i], &ione, mfqP->Gres, &ione));
183     /* General case: .5 * Gres= sum_ij[wij * (cjgi + cigj)] */
189       PetscCallBLAS("BLASaxpy", BLASaxpy_(&blasn, &factor, &mfqP->Fdiff[blasn * row], &ione, mfqP->Gres, &ione));
191       PetscCallBLAS("BLASaxpy", BLASaxpy_(&blasn, &factor, &mfqP->Fdiff[blasn * col], &ione, mfqP->Gres, &ione));
219     /* Default: Gres= sum_i[cigi] = G*c' */
221     PetscCallBLAS("BLASgemv", BLASgemv_("N", &blasn, &blasm, &one, mfqP->Fdiff, &blasn, mfqP->C, &ione, &zero, mfqP->Gres, &ione));
509     PetscCallBLAS("BLASdot", dp = BLASdot_(&blasn, &mfqP->Q_tmp[i * mfqP->npmax], &blas1, mfqP->Gres, &blas1));
513     /* mfqP->work[i] = Cres+Modeld(i,:)*(Gres+.5*Hres*Modeld(i,:)') */
514     for (j = 0; j < mfqP->n; j++) mfqP->work2[j] = mfqP->Gres[j];
741   PetscCall(VecSetValues(tao->gradient, mfqP->n, mfqP->indices, mfqP->Gres, INSERT_VALUES));
794       /* Cres += work*Gres + .5*work*Hres*work';
795        Gres += Hres*work'; */
798       for (i = 0; i < mfqP->n; i++) mfqP->Gres[i] += mfqP->work2[i];
886     PetscCall(VecSetValues(tao->gradient, mfqP->n, mfqP->indices, mfqP->Gres, INSERT_VALUES));
974   PetscCall(PetscMalloc1(mfqP->n, &mfqP->Gres));
1086   PetscCall(PetscFree(mfqP->Gres));