366 #define Mult_W_W picofftsg_mult_w_w  macro
1753             wd1r -= Mult_W_W(ss1, wk1i);  in cftb1st()
1754             wd1i += Mult_W_W(ss1, wk1r);  in cftb1st()
1755             wd3r -= Mult_W_W(ss3, wk3i);  in cftb1st()
1756             wd3i += Mult_W_W(ss3, wk3r);  in cftb1st()
1775             a[j2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftb1st()
1776             a[j2 + 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftb1st()
1779             a[j3] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i);  in cftb1st()
1780             a[j3 + 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r);  in cftb1st()
1795             a[j2 + 2] = Mult_W_W(wd1r, x0r) - Mult_W_W(wd1i, x0i);  in cftb1st()
1796             a[j2 + 3] = Mult_W_W(wd1r, x0i) + Mult_W_W(wd1i, x0r);  in cftb1st()
1799             a[j3 + 2] = Mult_W_W(wd3r, x0r) + Mult_W_W(wd3i, x0i);  in cftb1st()
1800             a[j3 + 3] = Mult_W_W(wd3r, x0i) - Mult_W_W(wd3i, x0r);  in cftb1st()
1819             a[j2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftb1st()
1820             a[j2 + 1] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftb1st()
1823             a[j3] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i);  in cftb1st()
1824             a[j3 + 1] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r);  in cftb1st()
1839             a[j2 - 2] = Mult_W_W(wd1i, x0r) - Mult_W_W(wd1r, x0i);  in cftb1st()
1840             a[j2 - 1] = Mult_W_W(wd1i, x0i) + Mult_W_W(wd1r, x0r);  in cftb1st()
1843             a[j3 - 2] = Mult_W_W(wd3i, x0r) + Mult_W_W(wd3r, x0i);  in cftb1st()
1844             a[j3 - 1] = Mult_W_W(wd3i, x0i) - Mult_W_W(wd3r, x0r);  in cftb1st()
1845             wk1r -= Mult_W_W(ss1, wd1i);  in cftb1st()
1846             wk1i += Mult_W_W(ss1, wd1r);  in cftb1st()
1847             wk3r -= Mult_W_W(ss3, wd3i);  in cftb1st()
1848             wk3i += Mult_W_W(ss3, wd3r);  in cftb1st()
1873     a[j2 - 2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftb1st()
1874     a[j2 - 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftb1st()
1877     a[j3 - 2] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i);  in cftb1st()
1878     a[j3 - 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r);  in cftb1st()
1913     a[j2 + 2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftb1st()
1914     a[j2 + 3] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftb1st()
1917     a[j3 + 2] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i);  in cftb1st()
1918     a[j3 + 3] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r);  in cftb1st()
2078             wd1r -= Mult_W_W(ss1, wk1i);  in cftmdl1()
2079             wd1i += Mult_W_W(ss1, wk1r);  in cftmdl1()
2080             wd3r -= Mult_W_W(ss3, wk3i);  in cftmdl1()
2081             wd3i += Mult_W_W(ss3, wk3r);  in cftmdl1()
2099             a[j2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftmdl1()
2100             a[j2 + 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftmdl1()
2103             a[j3] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i);  in cftmdl1()
2104             a[j3 + 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r);  in cftmdl1()
2119             a[j2 + 2] = Mult_W_W(wd1r, x0r) - Mult_W_W(wd1i, x0i);  in cftmdl1()
2120             a[j2 + 3] = Mult_W_W(wd1r, x0i) + Mult_W_W(wd1i, x0r);  in cftmdl1()
2123             a[j3 + 2] = Mult_W_W(wd3r, x0r) + Mult_W_W(wd3i, x0i);  in cftmdl1()
2124             a[j3 + 3] = Mult_W_W(wd3r, x0i) - Mult_W_W(wd3i, x0r);  in cftmdl1()
2143             a[j2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftmdl1()
2144             a[j2 + 1] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftmdl1()
2147             a[j3] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i);  in cftmdl1()
2148             a[j3 + 1] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r);  in cftmdl1()
2163             a[j2 - 2] = Mult_W_W(wd1i, x0r) - Mult_W_W(wd1r, x0i);  in cftmdl1()
2164             a[j2 - 1] = Mult_W_W(wd1i, x0i) + Mult_W_W(wd1r, x0r);  in cftmdl1()
2167             a[j3 - 2] = Mult_W_W(wd3i, x0r) + Mult_W_W(wd3r, x0i);  in cftmdl1()
2168             a[j3 - 1] = Mult_W_W(wd3i, x0i) - Mult_W_W(wd3r, x0r);  in cftmdl1()
2169             wk1r -= Mult_W_W(ss1, wd1i);  in cftmdl1()
2170             wk1i += Mult_W_W(ss1, wd1r);  in cftmdl1()
2171             wk3r -= Mult_W_W(ss3, wd3i);  in cftmdl1()
2172             wk3i += Mult_W_W(ss3, wd3r);  in cftmdl1()
2197     a[j2 - 2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftmdl1()
2198     a[j2 - 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftmdl1()
2201     a[j3 - 2] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i);  in cftmdl1()
2202     a[j3 - 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r);  in cftmdl1()
2237     a[j2 + 2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftmdl1()
2238     a[j2 + 3] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftmdl1()
2241     a[j3 + 2] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i);  in cftmdl1()
2242     a[j3 + 3] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r);  in cftmdl1()
2307             wl1r -= Mult_W_W(ss1, wk1i);  in cftmdl2()
2308             wl1i += Mult_W_W(ss1, wk1r);  in cftmdl2()
2309             wl3r -= Mult_W_W(ss3, wk3i);  in cftmdl2()
2310             wl3i += Mult_W_W(ss3, wk3r);  in cftmdl2()
2311             we1r -= Mult_W_W(ss1, wd1i);  in cftmdl2()
2312             we1i += Mult_W_W(ss1, wd1r);  in cftmdl2()
2313             we3r -= Mult_W_W(ss3, wd3i);  in cftmdl2()
2314             we3i += Mult_W_W(ss3, wd3r);  in cftmdl2()
2326             y0r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftmdl2()
2327             y0i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftmdl2()
2328             y2r = Mult_W_W(wd1r, x2r) - Mult_W_W(wd1i, x2i);  in cftmdl2()
2329             y2i = Mult_W_W(wd1r, x2i) + Mult_W_W(wd1i, x2r);  in cftmdl2()
2334             y0r = Mult_W_W(wk3r, x1r) + Mult_W_W(wk3i, x1i);  in cftmdl2()
2335             y0i = Mult_W_W(wk3r, x1i) - Mult_W_W(wk3i, x1r);  in cftmdl2()
2336             y2r = Mult_W_W(wd3r, x3r) + Mult_W_W(wd3i, x3i);  in cftmdl2()
2337             y2i = Mult_W_W(wd3r, x3i) - Mult_W_W(wd3i, x3r);  in cftmdl2()
2350             y0r = Mult_W_W(wl1r, x0r) - Mult_W_W(wl1i, x0i);  in cftmdl2()
2351             y0i = Mult_W_W(wl1r, x0i) + Mult_W_W(wl1i, x0r);  in cftmdl2()
2352             y2r = Mult_W_W(we1r, x2r) - Mult_W_W(we1i, x2i);  in cftmdl2()
2353             y2i = Mult_W_W(we1r, x2i) + Mult_W_W(we1i, x2r);  in cftmdl2()
2358             y0r = Mult_W_W(wl3r, x1r) + Mult_W_W(wl3i, x1i);  in cftmdl2()
2359             y0i = Mult_W_W(wl3r, x1i) - Mult_W_W(wl3i, x1r);  in cftmdl2()
2360             y2r = Mult_W_W(we3r, x3r) + Mult_W_W(we3i, x3i);  in cftmdl2()
2361             y2i = Mult_W_W(we3r, x3i) - Mult_W_W(we3i, x3r);  in cftmdl2()
2378             y0r = Mult_W_W(wd1i, x0r) - Mult_W_W(wd1r, x0i);  in cftmdl2()
2379             y0i = Mult_W_W(wd1i, x0i) + Mult_W_W(wd1r, x0r);  in cftmdl2()
2380             y2r = Mult_W_W(wk1i, x2r) - Mult_W_W(wk1r, x2i);  in cftmdl2()
2381             y2i = Mult_W_W(wk1i, x2i) + Mult_W_W(wk1r, x2r);  in cftmdl2()
2386             y0r = Mult_W_W(wd3i, x1r) + Mult_W_W(wd3r, x1i);  in cftmdl2()
2387             y0i = Mult_W_W(wd3i, x1i) - Mult_W_W(wd3r, x1r);  in cftmdl2()
2388             y2r = Mult_W_W(wk3i, x3r) + Mult_W_W(wk3r, x3i);  in cftmdl2()
2389             y2i = Mult_W_W(wk3i, x3i) - Mult_W_W(wk3r, x3r);  in cftmdl2()
2402             y0r = Mult_W_W(we1i, x0r) - Mult_W_W(we1r, x0i);  in cftmdl2()
2403             y0i = Mult_W_W(we1i, x0i) + Mult_W_W(we1r, x0r);  in cftmdl2()
2404             y2r = Mult_W_W(wl1i, x2r) - Mult_W_W(wl1r, x2i);  in cftmdl2()
2405             y2i = Mult_W_W(wl1i, x2i) + Mult_W_W(wl1r, x2r);  in cftmdl2()
2410             y0r = Mult_W_W(we3i, x1r) + Mult_W_W(we3r, x1i);  in cftmdl2()
2411             y0i = Mult_W_W(we3i, x1i) - Mult_W_W(we3r, x1r);  in cftmdl2()
2412             y2r = Mult_W_W(wl3i, x3r) + Mult_W_W(wl3r, x3i);  in cftmdl2()
2413             y2i = Mult_W_W(wl3i, x3i) - Mult_W_W(wl3r, x3r);  in cftmdl2()
2418             wk1r -= Mult_W_W(ss1, wl1i);  in cftmdl2()
2419             wk1i += Mult_W_W(ss1, wl1r);  in cftmdl2()
2420             wk3r -= Mult_W_W(ss3, wl3i);  in cftmdl2()
2421             wk3i += Mult_W_W(ss3, wl3r);  in cftmdl2()
2422             wd1r -= Mult_W_W(ss1, we1i);  in cftmdl2()
2423             wd1i += Mult_W_W(ss1, we1r);  in cftmdl2()
2424             wd3r -= Mult_W_W(ss3, we3i);  in cftmdl2()
2425             wd3i += Mult_W_W(ss3, we3r);  in cftmdl2()
2445     y0r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftmdl2()
2446     y0i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftmdl2()
2447     y2r = Mult_W_W(wd1r, x2r) - Mult_W_W(wd1i, x2i);  in cftmdl2()
2448     y2i = Mult_W_W(wd1r, x2i) + Mult_W_W(wd1i, x2r);  in cftmdl2()
2453     y0r = Mult_W_W(wk3r, x1r) + Mult_W_W(wk3i, x1i);  in cftmdl2()
2454     y0i = Mult_W_W(wk3r, x1i) - Mult_W_W(wk3i, x1r);  in cftmdl2()
2455     y2r = Mult_W_W(wd3r, x3r) + Mult_W_W(wd3i, x3i);  in cftmdl2()
2456     y2i = Mult_W_W(wd3r, x3i) - Mult_W_W(wd3i, x3r);  in cftmdl2()
2469     y0r = Mult_W_W(wl1r, x0r) - Mult_W_W(wl1i, x0i);  in cftmdl2()
2470     y0i = Mult_W_W(wl1r, x0i) + Mult_W_W(wl1i, x0r);  in cftmdl2()
2471     y2r = Mult_W_W(wl1i, x2r) - Mult_W_W(wl1r, x2i);  in cftmdl2()
2472     y2i = Mult_W_W(wl1i, x2i) + Mult_W_W(wl1r, x2r);  in cftmdl2()
2477     y0r = Mult_W_W(wl1i, x1r) - Mult_W_W(wl1r, x1i);  in cftmdl2()
2478     y0i = Mult_W_W(wl1i, x1i) + Mult_W_W(wl1r, x1r);  in cftmdl2()
2479     y2r = Mult_W_W(wl1r, x3r) - Mult_W_W(wl1i, x3i);  in cftmdl2()
2480     y2i = Mult_W_W(wl1r, x3i) + Mult_W_W(wl1i, x3r);  in cftmdl2()
2493     y0r = Mult_W_W(wd1i, x0r) - Mult_W_W(wd1r, x0i);  in cftmdl2()
2494     y0i = Mult_W_W(wd1i, x0i) + Mult_W_W(wd1r, x0r);  in cftmdl2()
2495     y2r = Mult_W_W(wk1i, x2r) - Mult_W_W(wk1r, x2i);  in cftmdl2()
2496     y2i = Mult_W_W(wk1i, x2i) + Mult_W_W(wk1r, x2r);  in cftmdl2()
2501     y0r = Mult_W_W(wd3i, x1r) + Mult_W_W(wd3r, x1i);  in cftmdl2()
2502     y0i = Mult_W_W(wd3i, x1i) - Mult_W_W(wd3r, x1r);  in cftmdl2()
2503     y2r = Mult_W_W(wk3i, x3r) + Mult_W_W(wk3r, x3i);  in cftmdl2()
2504     y2i = Mult_W_W(wk3i, x3i) - Mult_W_W(wk3r, x3r);  in cftmdl2()
2571     y9r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftf161()
2572     y9i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftf161()
2575     y13r = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftf161()
2576     y13i = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftf161()
2611     y11r = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftf161()
2612     y11i = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftf161()
2615     y15r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftf161()
2616     y15i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftf161()
2726     x1r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftf162()
2727     x1i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftf162()
2730     x2r = Mult_W_W(wk3i, x0r) - Mult_W_W(wk3r, x0i);  in cftf162()
2731     x2i = Mult_W_W(wk3i, x0i) + Mult_W_W(wk3r, x0r);  in cftf162()
2738     x1r = Mult_W_W(wk3r, x0r) - Mult_W_W(wk3i, x0i);  in cftf162()
2739     x1i = Mult_W_W(wk3r, x0i) + Mult_W_W(wk3i, x0r);  in cftf162()
2742     x2r = Mult_W_W(wk1r, x0r) + Mult_W_W(wk1i, x0i);  in cftf162()
2743     x2i = Mult_W_W(wk1r, x0i) - Mult_W_W(wk1i, x0r);  in cftf162()
2750     x1r = Mult_W_W(wk2r, x0r) - Mult_W_W(wk2i, x0i);  in cftf162()
2751     x1i = Mult_W_W(wk2r, x0i) + Mult_W_W(wk2i, x0r);  in cftf162()
2754     x2r = Mult_W_W(wk2i, x0r) - Mult_W_W(wk2r, x0i);  in cftf162()
2755     x2i = Mult_W_W(wk2i, x0i) + Mult_W_W(wk2r, x0r);  in cftf162()
2762     x1r = Mult_W_W(wk2i, x0r) - Mult_W_W(wk2r, x0i);  in cftf162()
2763     x1i = Mult_W_W(wk2i, x0i) + Mult_W_W(wk2r, x0r);  in cftf162()
2766     x2r = Mult_W_W(wk2r, x0r) - Mult_W_W(wk2i, x0i);  in cftf162()
2767     x2i = Mult_W_W(wk2r, x0i) + Mult_W_W(wk2i, x0r);  in cftf162()
2774     x1r = Mult_W_W(wk3r, x0r) - Mult_W_W(wk3i, x0i);  in cftf162()
2775     x1i = Mult_W_W(wk3r, x0i) + Mult_W_W(wk3i, x0r);  in cftf162()
2778     x2r = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftf162()
2779     x2i = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftf162()
2786     x1r = Mult_W_W(wk1i, x0r) + Mult_W_W(wk1r, x0i);  in cftf162()
2787     x1i = Mult_W_W(wk1i, x0i) - Mult_W_W(wk1r, x0r);  in cftf162()
2790     x2r = Mult_W_W(wk3i, x0r) - Mult_W_W(wk3r, x0i);  in cftf162()
2791     x2i = Mult_W_W(wk3i, x0i) + Mult_W_W(wk3r, x0r);  in cftf162()
2956     y4r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftf082()
2957     y4i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftf082()
2960     y5r = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftf082()
2961     y5i = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftf082()
2964     y6r = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i);  in cftf082()
2965     y6i = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r);  in cftf082()
2968     y7r = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i);  in cftf082()
2969     y7i = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r);  in cftf082()
3104             yr = Mult_W_W(wdr, xr) - Mult_W_W(wdi, xi);  in rftfsub()
3105             yi = Mult_W_W(wdr, xi) + Mult_W_W(wdi, xr);  in rftfsub()
3110             wkr += Mult_W_W(ss, wdi);  in rftfsub()
3114             yr = Mult_W_W(wkr, xr) - Mult_W_W(wki, xi);  in rftfsub()
3115             yi = Mult_W_W(wkr, xi) + Mult_W_W(wki, xr);  in rftfsub()
3120             wdr += Mult_W_W(ss, wki);  in rftfsub()
3130     yr = Mult_W_W(wdr, xr) - Mult_W_W(wdi, xi);  in rftfsub()
3131     yi = Mult_W_W(wdr, xi) + Mult_W_W(wdi, xr);  in rftfsub()
3164             yr = Mult_W_W(wdr, xr) + Mult_W_W(wdi, xi);  in rftbsub()
3165             yi = Mult_W_W(wdr, xi) - Mult_W_W(wdi, xr);  in rftbsub()
3170             wkr += Mult_W_W(ss, wdi);  in rftbsub()
3174             yr = Mult_W_W(wkr, xr) + Mult_W_W(wki, xi);  in rftbsub()
3175             yi = Mult_W_W(wkr, xi) - Mult_W_W(wki, xr);  in rftbsub()
3180             wdr += Mult_W_W(ss, wki);  in rftbsub()
3189     yr = Mult_W_W(wdr, xr) + Mult_W_W(wdi, xi);  in rftbsub()
3190     yi = Mult_W_W(wdr, xi) - Mult_W_W(wdi, xr);  in rftbsub()
3233             wkr -= Mult_W_W(ss, wdi);  in dctsub()
3234             wki += Mult_W_W(ss, wdr);  in dctsub()
3235             yr = Mult_W_W(wki, a[j]) - Mult_W_W(wkr, a[k]);  in dctsub()
3236             yi = Mult_W_W(wkr, a[j]) + Mult_W_W(wki, a[k]);  in dctsub()
3237             wdr -= Mult_W_W(ss, wki);  in dctsub()
3238             wdi += Mult_W_W(ss, wkr);  in dctsub()
3251     a[m] = Mult_W_W(WR5000, a[m]);  in dctsub()
3263         wdr = Mult_W_W(wki, WI2500);  in dctsub4()
3264         wdi = Mult_W_W(wki, WR2500);  in dctsub4()
3265         xr = Mult_W_W(wdi, a[1]) - Mult_W_W(wdr, a[3]);  in dctsub4()
3266         a[1] = Mult_W_W(wdr, a[1]) + Mult_W_W(wdi, a[3]);  in dctsub4()
3269     a[m] = Mult_W_W(wki, a[m]);  in dctsub4()