25 qm[0] -= dx[0]; qm[1] -= dx[1]; qm[2] -= dx[2]; qm[3] -= dx[3];
26 qm[4] -= dx[4]; qm[5] -= dx[5]; qm[6] -= dx[6]; qm[7] -= dx[7];
27 }
else if (*error > 0) {
28 qm[0] += dx[0]; qm[1] += dx[1]; qm[2] += dx[2]; qm[3] += dx[3];
29 qm[4] += dx[4]; qm[5] += dx[5]; qm[6] += dx[6]; qm[7] += dx[7];
32 round += dl[0] * qm[0] + dl[1] * qm[1] + dl[2] * qm[2] + dl[3] * qm[3] +
33 dl[4] * qm[4] + dl[5] * qm[5] + dl[6] * qm[6] + dl[7] * qm[7];
35 dx[0] = dx[1]; dx[1] = dx[2]; dx[2] = dx[3]; dx[3] = dx[4];
36 dl[0] = dl[1]; dl[1] = dl[2]; dl[2] = dl[3]; dl[3] = dl[4];
38 dx[4] = ((dl[4] >> 30) | 1);
39 dx[5] = ((dl[5] >> 30) | 2) & ~1;
40 dx[6] = ((dl[6] >> 30) | 2) & ~1;
41 dx[7] = ((dl[7] >> 30) | 4) & ~3;
44 *in += (round >>
shift);
46 dl[4] = -dl[5]; dl[5] = -dl[6];
47 dl[6] = *in - dl[7]; dl[7] = *
in;
48 dl[5] += dl[6]; dl[4] += dl[5];