static void
increment_state (ranlxs_state_t * state)
{
int k, kmax, m;
double x, y1, y2, y3;
float *xflt = state->xflt;
double *xdbl = state->xdbl;
double *ydbl = state->ydbl;
double carry = state->carry;
unsigned int ir = state->ir;
unsigned int jr = state->jr;
for (k = 0; ir > 0; ++k)
{
y1 = xdbl[jr] - xdbl[ir];
y2 = y1 - carry;
if (y2 < 0)
{
carry = one_bit;
y2 += 1;
}
else
{
carry = 0;
}
xdbl[ir] = y2;
ir = next[ir];
jr = next[jr];
}
kmax = state->pr - 12;
for (; k <= kmax; k += 12)
{
y1 = xdbl[7] - xdbl[0];
y1 -= carry;
RANLUX_STEP (y2, y1, 8, 1, 0);
RANLUX_STEP (y3, y2, 9, 2, 1);
RANLUX_STEP (y1, y3, 10, 3, 2);
RANLUX_STEP (y2, y1, 11, 4, 3);
RANLUX_STEP (y3, y2, 0, 5, 4);
RANLUX_STEP (y1, y3, 1, 6, 5);
RANLUX_STEP (y2, y1, 2, 7, 6);
RANLUX_STEP (y3, y2, 3, 8, 7);
RANLUX_STEP (y1, y3, 4, 9, 8);
RANLUX_STEP (y2, y1, 5, 10, 9);
RANLUX_STEP (y3, y2, 6, 11, 10);
if (y3 < 0)
{
carry = one_bit;
y3 += 1;
}
else
{
carry = 0;
}
xdbl[11] = y3;
}
kmax = state->pr;
for (; k < kmax; ++k)
{
y1 = xdbl[jr] - xdbl[ir];
y2 = y1 - carry;
if (y2 < 0)
{
carry = one_bit;
y2 += 1;
}
else
{
carry = 0;
}
xdbl[ir] = y2;
ydbl[ir] = y2 + shift;
ir = next[ir];
jr = next[jr];
}
ydbl[ir] = xdbl[ir] + shift;
for (k = next[ir]; k > 0;)
{
ydbl[k] = xdbl[k] + shift;
k = next[k];
}
for (k = 0, m = 0; k < 12; ++k)
{
x = xdbl[k];
y2 = ydbl[k] - shift;
if (y2 > x)
y2 -= sone_bit;
y1 = (x - y2) * sbase;
xflt[m++] = (float) y1;
xflt[m++] = (float) y2;
}
state->ir = ir;
state->is = 2 * ir;
state->is_old = 2 * ir;
state->jr = jr;
state->carry = carry;
}
static inline void
increment_state (ranlxd_state_t * state)
{
int k, kmax;
double y1, y2, y3;
double *xdbl = state->xdbl;
double carry = state->carry;
unsigned int ir = state->ir;
unsigned int jr = state->jr;
for (k = 0; ir > 0; ++k)
{
y1 = xdbl[jr] - xdbl[ir];
y2 = y1 - carry;
if (y2 < 0)
{
carry = one_bit;
y2 += 1;
}
else
{
carry = 0;
}
xdbl[ir] = y2;
ir = next[ir];
jr = next[jr];
}
kmax = state->pr - 12;
for (; k <= kmax; k += 12)
{
y1 = xdbl[7] - xdbl[0];
y1 -= carry;
RANLUX_STEP (y2, y1, 8, 1, 0);
RANLUX_STEP (y3, y2, 9, 2, 1);
RANLUX_STEP (y1, y3, 10, 3, 2);
RANLUX_STEP (y2, y1, 11, 4, 3);
RANLUX_STEP (y3, y2, 0, 5, 4);
RANLUX_STEP (y1, y3, 1, 6, 5);
RANLUX_STEP (y2, y1, 2, 7, 6);
RANLUX_STEP (y3, y2, 3, 8, 7);
RANLUX_STEP (y1, y3, 4, 9, 8);
RANLUX_STEP (y2, y1, 5, 10, 9);
RANLUX_STEP (y3, y2, 6, 11, 10);
if (y3 < 0)
{
carry = one_bit;
y3 += 1;
}
else
{
carry = 0;
}
xdbl[11] = y3;
}
kmax = state->pr;
for (; k < kmax; ++k)
{
y1 = xdbl[jr] - xdbl[ir];
y2 = y1 - carry;
if (y2 < 0)
{
carry = one_bit;
y2 += 1;
}
else
{
carry = 0;
}
xdbl[ir] = y2;
ir = next[ir];
jr = next[jr];
}
state->ir = ir;
state->ir_old = ir;
state->jr = jr;
state->carry = carry;
}
