?? hb_25.c
字號:
E T2W, T4d, T33, T4c, T2S, T32; T2S = T2Q - T2R; T2W = T2S - T2V; T4d = T2S + T2V; T32 = T30 - T31; T33 = T2Z + T32; T4c = T32 - T2Z; T34 = FNMS(KP998026728, T33, KP062790519 * T2W); T4p = FNMS(KP637423989, T4c, KP770513242 * T4d); T3l = FMA(KP062790519, T33, KP998026728 * T2W); T4e = FMA(KP770513242, T4c, KP637423989 * T4d); } } { E TM, TQ, T1U, T1L, T1N, T1Z, T1t, T1V, T1y, T1Y; { E TO, TP, T1B, T1M; TO = KP559016994 * (Ts - TL); TM = Ts + TL; TP = FNMS(KP250000000, TM, T9); TQ = TO + TP; T1U = TP - TO; T1B = KP559016994 * (T1z - T1A); T1L = T1z + T1A; T1M = FNMS(KP250000000, T1L, T1K); T1N = T1B + T1M; T1Z = T1M - T1B; } { E T19, T1s, T1w, T1x; T19 = TZ - T18; T1s = T1i - T1r; T1t = FMA(KP951056516, T19, KP587785252 * T1s); T1V = FNMS(KP951056516, T1s, KP587785252 * T19); T1w = Ti - Tr; T1x = TB - TK; T1y = FMA(KP951056516, T1w, KP587785252 * T1x); T1Y = FNMS(KP951056516, T1x, KP587785252 * T1w); } cr[0] = T9 + TM; ci[0] = T1K + T1L; { E T1u, T1O, TN, T1v; T1u = TQ - T1t; T1O = T1y + T1N; TN = W[8]; T1v = W[9]; cr[WS(rs, 5)] = FNMS(T1v, T1O, TN * T1u); ci[WS(rs, 5)] = FMA(T1v, T1u, TN * T1O); } { E T22, T24, T21, T23; T22 = T1U + T1V; T24 = T1Z - T1Y; T21 = W[28]; T23 = W[29]; cr[WS(rs, 15)] = FNMS(T23, T24, T21 * T22); ci[WS(rs, 15)] = FMA(T23, T22, T21 * T24); } { E T1W, T20, T1T, T1X; T1W = T1U - T1V; T20 = T1Y + T1Z; T1T = W[18]; T1X = W[19]; cr[WS(rs, 10)] = FNMS(T1X, T20, T1T * T1W); ci[WS(rs, 10)] = FMA(T1X, T1W, T1T * T20); } { E T1Q, T1S, T1P, T1R; T1Q = TQ + T1t; T1S = T1N - T1y; T1P = W[38]; T1R = W[39]; cr[WS(rs, 20)] = FNMS(T1R, T1S, T1P * T1Q); ci[WS(rs, 20)] = FMA(T1R, T1Q, T1P * T1S); } } { E T6H, T71, T6M, T74, T6i, T6x, T6y, T6z, T6Q, T6R, T6P, T6S; { E T6D, T6G, T6K, T6L; T6D = T6B + T6C; T6G = T6E - T6F; T6H = FMA(KP951056516, T6D, KP587785252 * T6G); T71 = FNMS(KP951056516, T6G, KP587785252 * T6D); T6K = T6l - T6o; T6L = T6v - T6s; T6M = FMA(KP951056516, T6K, KP587785252 * T6L); T74 = FNMS(KP951056516, T6L, KP587785252 * T6K); } { E T6p, T6w, T6N, T6O; T6i = T5v + T5w; T6p = T6l + T6o; T6w = T6s + T6v; T6x = T6p - T6w; T6y = FNMS(KP250000000, T6x, T6i); T6z = KP559016994 * (T6p + T6w); T6Q = T5R - T5Q; T6N = T6B - T6C; T6O = T6E + T6F; T6R = T6N + T6O; T6P = KP559016994 * (T6N - T6O); T6S = FNMS(KP250000000, T6R, T6Q); } { E T7c, T7e, T7b, T7d; T7c = T6i + T6x; T7e = T6Q + T6R; T7b = W[6]; T7d = W[7]; cr[WS(rs, 4)] = FNMS(T7d, T7e, T7b * T7c); ci[WS(rs, 4)] = FMA(T7d, T7c, T7b * T7e); } { E T72, T78, T76, T7a, T70, T75; T70 = T6y - T6z; T72 = T70 - T71; T78 = T70 + T71; T75 = T6S - T6P; T76 = T74 + T75; T7a = T75 - T74; { E T6Z, T73, T77, T79; T6Z = W[26]; T73 = W[27]; cr[WS(rs, 14)] = FNMS(T73, T76, T6Z * T72); ci[WS(rs, 14)] = FMA(T73, T72, T6Z * T76); T77 = W[36]; T79 = W[37]; cr[WS(rs, 19)] = FNMS(T79, T7a, T77 * T78); ci[WS(rs, 19)] = FMA(T79, T78, T77 * T7a); } } { E T6I, T6W, T6U, T6Y, T6A, T6T; T6A = T6y + T6z; T6I = T6A - T6H; T6W = T6A + T6H; T6T = T6P + T6S; T6U = T6M + T6T; T6Y = T6T - T6M; { E T6h, T6J, T6V, T6X; T6h = W[16]; T6J = W[17]; cr[WS(rs, 9)] = FNMS(T6J, T6U, T6h * T6I); ci[WS(rs, 9)] = FMA(T6J, T6I, T6h * T6U); T6V = W[46]; T6X = W[47]; cr[WS(rs, 24)] = FNMS(T6X, T6Y, T6V * T6W); ci[WS(rs, 24)] = FMA(T6X, T6W, T6V * T6Y); } } } { E T3n, T3N, T3s, T3Q, T3d, T3e, T36, T3f, T3C, T3D, T3v, T3E; { E T3j, T3m, T3q, T3r; T3j = T3h - T3i; T3m = T3k - T3l; T3n = FMA(KP951056516, T3j, KP587785252 * T3m); T3N = FNMS(KP951056516, T3m, KP587785252 * T3j); T3q = T2k + T2z; T3r = T2P - T34; T3s = FMA(KP951056516, T3q, KP587785252 * T3r); T3Q = FNMS(KP951056516, T3r, KP587785252 * T3q); } { E T2A, T35, T3t, T3u; T3d = T39 - T3c; T2A = T2k - T2z; T35 = T2P + T34; T3e = T2A + T35; T36 = KP559016994 * (T2A - T35); T3f = FNMS(KP250000000, T3e, T3d); T3C = T3y + T3B; T3t = T3h + T3i; T3u = T3k + T3l; T3D = T3t + T3u; T3v = KP559016994 * (T3t - T3u); T3E = FNMS(KP250000000, T3D, T3C); } { E T3Y, T40, T3X, T3Z; T3Y = T3d + T3e; T40 = T3C + T3D; T3X = W[2]; T3Z = W[3]; cr[WS(rs, 2)] = FNMS(T3Z, T40, T3X * T3Y); ci[WS(rs, 2)] = FMA(T3Z, T3Y, T3X * T40); } { E T3O, T3U, T3S, T3W, T3M, T3R; T3M = T3f - T36; T3O = T3M - T3N; T3U = T3M + T3N; T3R = T3E - T3v; T3S = T3Q + T3R; T3W = T3R - T3Q; { E T3L, T3P, T3T, T3V; T3L = W[22]; T3P = W[23]; cr[WS(rs, 12)] = FNMS(T3P, T3S, T3L * T3O); ci[WS(rs, 12)] = FMA(T3P, T3O, T3L * T3S); T3T = W[32]; T3V = W[33]; cr[WS(rs, 17)] = FNMS(T3V, T3W, T3T * T3U); ci[WS(rs, 17)] = FMA(T3V, T3U, T3T * T3W); } } { E T3o, T3I, T3G, T3K, T3g, T3F; T3g = T36 + T3f; T3o = T3g - T3n; T3I = T3g + T3n; T3F = T3v + T3E; T3G = T3s + T3F; T3K = T3F - T3s; { E T25, T3p, T3H, T3J; T25 = W[12]; T3p = W[13]; cr[WS(rs, 7)] = FNMS(T3p, T3G, T25 * T3o); ci[WS(rs, 7)] = FMA(T3p, T3o, T25 * T3G); T3H = W[42]; T3J = W[43]; cr[WS(rs, 22)] = FNMS(T3J, T3K, T3H * T3I); ci[WS(rs, 22)] = FMA(T3J, T3I, T3H * T3K); } } } { E T4r, T4L, T4w, T4O, T4h, T4i, T4g, T4j, T4A, T4B, T4z, T4C; { E T4n, T4q, T4u, T4v; T4n = T4l - T4m; T4q = T4o - T4p; T4r = FMA(KP951056516, T4n, KP587785252 * T4q); T4L = FNMS(KP951056516, T4q, KP587785252 * T4n); T4u = T44 + T47; T4v = T4b + T4e; T4w = FMA(KP951056516, T4u, KP587785252 * T4v); T4O = FNMS(KP951056516, T4v, KP587785252 * T4u); } { E T48, T4f, T4x, T4y; T4h = T39 + T3c; T48 = T44 - T47; T4f = T4b - T4e; T4i = T48 + T4f; T4g = KP559016994 * (T48 - T4f); T4j = FNMS(KP250000000, T4i, T4h); T4A = T3B - T3y; T4x = T4l + T4m; T4y = T4o + T4p; T4B = T4x + T4y; T4z = KP559016994 * (T4x - T4y); T4C = FNMS(KP250000000, T4B, T4A); } { E T4W, T4Y, T4V, T4X; T4W = T4h + T4i; T4Y = T4A + T4B; T4V = W[4]; T4X = W[5]; cr[WS(rs, 3)] = FNMS(T4X, T4Y, T4V * T4W); ci[WS(rs, 3)] = FMA(T4X, T4W, T4V * T4Y); } { E T4M, T4S, T4Q, T4U, T4K, T4P; T4K = T4j - T4g; T4M = T4K - T4L; T4S = T4K + T4L; T4P = T4C - T4z; T4Q = T4O + T4P; T4U = T4P - T4O; { E T4J, T4N, T4R, T4T; T4J = W[24]; T4N = W[25]; cr[WS(rs, 13)] = FNMS(T4N, T4Q, T4J * T4M); ci[WS(rs, 13)] = FMA(T4N, T4M, T4J * T4Q); T4R = W[34]; T4T = W[35]; cr[WS(rs, 18)] = FNMS(T4T, T4U, T4R * T4S); ci[WS(rs, 18)] = FMA(T4T, T4S, T4R * T4U); } } { E T4s, T4G, T4E, T4I, T4k, T4D; T4k = T4g + T4j; T4s = T4k - T4r; T4G = T4k + T4r; T4D = T4z + T4C; T4E = T4w + T4D; T4I = T4D - T4w; { E T41, T4t, T4F, T4H; T41 = W[14]; T4t = W[15]; cr[WS(rs, 8)] = FNMS(T4t, T4E, T41 * T4s); ci[WS(rs, 8)] = FMA(T4t, T4s, T41 * T4E); T4F = W[44]; T4H = W[45]; cr[WS(rs, 23)] = FNMS(T4H, T4I, T4F * T4G); ci[WS(rs, 23)] = FMA(T4H, T4G, T4F * T4I); } } } { E T5H, T63, T5M, T66, T5x, T5y, T5u, T5z, T5S, T5T, T5P, T5U; { E T5D, T5G, T5K, T5L; T5D = T5B - T5C; T5G = T5E - T5F; T5H = FMA(KP951056516, T5D, KP587785252 * T5G); T63 = FNMS(KP951056516, T5G, KP587785252 * T5D); T5K = T56 - T5d; T5L = T5l - T5s; T5M = FMA(KP951056516, T5K, KP587785252 * T5L); T66 = FNMS(KP951056516, T5L, KP587785252 * T5K); } { E T5e, T5t, T5N, T5O; T5x = T5v - T5w; T5e = T56 + T5d; T5t = T5l + T5s; T5y = T5e + T5t; T5u = KP559016994 * (T5e - T5t); T5z = FNMS(KP250000000, T5y, T5x); T5S = T5Q + T5R; T5N = T5B + T5C; T5O = T5E + T5F; T5T = T5N + T5O; T5P = KP559016994 * (T5N - T5O); T5U = FNMS(KP250000000, T5T, T5S); } { E T6e, T6g, T6d, T6f; T6e = T5x + T5y; T6g = T5S + T5T; T6d = W[0]; T6f = W[1]; cr[WS(rs, 1)] = FNMS(T6f, T6g, T6d * T6e); ci[WS(rs, 1)] = FMA(T6f, T6e, T6d * T6g); } { E T64, T6a, T68, T6c, T62, T67; T62 = T5z - T5u; T64 = T62 - T63; T6a = T62 + T63; T67 = T5U - T5P; T68 = T66 + T67; T6c = T67 - T66; { E T61, T65, T69, T6b; T61 = W[20]; T65 = W[21]; cr[WS(rs, 11)] = FNMS(T65, T68, T61 * T64); ci[WS(rs, 11)] = FMA(T65, T64, T61 * T68); T69 = W[30]; T6b = W[31]; cr[WS(rs, 16)] = FNMS(T6b, T6c, T69 * T6a); ci[WS(rs, 16)] = FMA(T6b, T6a, T69 * T6c); } } { E T5I, T5Y, T5W, T60, T5A, T5V; T5A = T5u + T5z; T5I = T5A - T5H; T5Y = T5A + T5H; T5V = T5P + T5U; T5W = T5M + T5V; T60 = T5V - T5M; { E T4Z, T5J, T5X, T5Z; T4Z = W[10]; T5J = W[11]; cr[WS(rs, 6)] = FNMS(T5J, T5W, T4Z * T5I); ci[WS(rs, 6)] = FMA(T5J, T5I, T4Z * T5W); T5X = W[40]; T5Z = W[41]; cr[WS(rs, 21)] = FNMS(T5Z, T60, T5X * T5Y); ci[WS(rs, 21)] = FMA(T5Z, T5Y, T5X * T60); } } } }}static const tw_instr twinstr[] = { {TW_FULL, 1, 25}, {TW_NEXT, 1, 0}};static const hc2hc_desc desc = { 25, "hb_25", twinstr, &GENUS, {260, 140, 140, 0} };void X(codelet_hb_25) (planner *p) { X(khc2hc_register) (p, hb_25, &desc);}#endif /* HAVE_FMA */
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -