- for (i = 1; i < TABLE_4_3_SIZE; i++) {
- double f, fm;
- int e, m;
- double value = i / 4;
- if ((i & 3) == 0)
- pow43_val = value / IMDCT_SCALAR * cbrt(value);
- f = pow43_val * exp2_lut[i & 3];
- fm = frexp(f, &e);
- m = llrint(fm * (1LL << 31));
- e += FRAC_BITS - 31 + 5 - 100;
-
- /* normalized to FRAC_BITS */
- table_4_3_value[i] = m;
- table_4_3_exp[i] = -e;
- }