74     double *grad = QP->gradient; /* gradient at current x */  in QPBoundary()  local
94     double *D = QP->D; /* diagonal of quadratic */  in QPBoundary()  local
128             if (grad[k] > 0.0)  in QPBoundary()
181                 grad[Ei[p]] += delta_xk * ((Ex) ? Ex[p] : 1);  in QPBoundary()
183             grad[k] += delta_xk * D[k];  in QPBoundary()
220                 if (0.5 * D[k] + grad[k] >= 0) /* flip lowers cost */  in QPBoundary()
233                 if (grad[k] - 0.5 * D[k] <= 0) /* flip lowers cost */  in QPBoundary()
249                     grad[Ei[p]] += (Ex) ? Ex[p] : 1;  in QPBoundary()
251                 grad[k] += D[k];  in QPBoundary()
257                     grad[Ei[p]] -= (Ex) ? Ex[p] : 1;  in QPBoundary()
259                 grad[k] -= D[k];  in QPBoundary()
351                 double c1  = (grad[j] - .5 * D[j] * dxj) * dxj  in QPBoundary()
352                             + (grad[i] - .5 * D[i] * dxi) * dxi;  in QPBoundary()
367                 double c2 = (grad[j] - 0.5 * D[j] * dxj) * dxj  in QPBoundary()
368                             + (grad[i] - 0.5 * D[i] * dxi) * dxi;  in QPBoundary()
414                     grad[Ei[p]] -= ((Ex) ? Ex[p] : 1) * dxj;  in QPBoundary()
418                     grad[Ei[p]] -= ((Ex) ? Ex[p] : 1) * dxi;  in QPBoundary()
420                 grad[j] -= D[j] * dxj;  in QPBoundary()
421                 grad[i] -= D[i] * dxi;  in QPBoundary()
532         double dxj, dxi, s = grad[j] / aj - grad[i] / ai;  in QPBoundary()
579             grad[Ei[k]] -= ((Ex) ? Ex[k] : 1) * dxj;  in QPBoundary()
583             grad[Ei[k]] -= ((Ex) ? Ex[k] : 1) * dxi;  in QPBoundary()
585         grad[j] -= D[j] * dxj;  in QPBoundary()
586         grad[i] -= D[i] * dxi;  in QPBoundary()
671         double c1 = (grad[j] - 0.5 * D[j] * dxj) * dxj;  in QPBoundary()
672         double c2 = (grad[j] - 0.5 * D[j] * dxi) * dxi;  in QPBoundary()
725                 grad[Ei[p]] -= ((Ex) ? Ex[p] : 1) * dxj;  in QPBoundary()
727             grad[j] -= D[j] * dxj;  in QPBoundary()