void fillrhs_aff(void) {
int i;
double *r1, *r2, *r3, *r4;
r1 = work.rhs;
r2 = work.rhs + 114;
r3 = work.rhs + 198;
r4 = work.rhs + 282;
/* r1 = -A^Ty - G^Tz - Px - q. */
multbymAT(r1, work.y);
multbymGT(work.buffer, work.z);
for (i = 0; i < 114; i++)
r1[i] += work.buffer[i];
multbyP(work.buffer, work.x);
for (i = 0; i < 114; i++)
r1[i] -= work.buffer[i] + work.q[i];
/* r2 = -z. */
for (i = 0; i < 84; i++)
r2[i] = -work.z[i];
/* r3 = -Gx - s + h. */
multbymG(r3, work.x);
for (i = 0; i < 84; i++)
r3[i] += -work.s[i] + work.h[i];
/* r4 = -Ax + b. */
multbymA(r4, work.x);
for (i = 0; i < 12; i++)
r4[i] += work.b[i];
}
double eval_objv(void) {
int i;
double objv;
/* Borrow space in work.rhs. */
multbyP(work.rhs, work.x);
objv = 0;
for (i = 0; i < 308; i++)
objv += work.x[i]*work.rhs[i];
objv *= 0.5;
for (i = 0; i < 308; i++)
objv += work.q[i]*work.x[i];
objv += 0;
return objv;
}
double eval_objv(void) {
int i;
double objv;
/* Borrow space in work.rhs. */
multbyP(work.rhs, work.x);
objv = 0;
for (i = 0; i < 114; i++)
objv += work.x[i]*work.rhs[i];
objv *= 0.5;
for (i = 0; i < 114; i++)
objv += work.q[i]*work.x[i];
objv += work.quad_919544893440[0];
return objv;
}
double eval_objv(void) {
int i;
double objv;
/* Borrow space in work.rhs. */
multbyP(work.rhs, work.x);
objv = 0;
for (i = 0; i < 9; i++)
objv += work.x[i]*work.rhs[i];
objv *= 0.5;
for (i = 0; i < 9; i++)
objv += work.q[i]*work.x[i];
objv += work.quad_153257975808[0];
return objv;
}
