/src/fftw3/rdft/scalar/r2cf/r2cfII_3.c
Line | Count | Source (jump to first uncovered line) |
1 | | /* |
2 | | * Copyright (c) 2003, 2007-14 Matteo Frigo |
3 | | * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology |
4 | | * |
5 | | * This program is free software; you can redistribute it and/or modify |
6 | | * it under the terms of the GNU General Public License as published by |
7 | | * the Free Software Foundation; either version 2 of the License, or |
8 | | * (at your option) any later version. |
9 | | * |
10 | | * This program is distributed in the hope that it will be useful, |
11 | | * but WITHOUT ANY WARRANTY; without even the implied warranty of |
12 | | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
13 | | * GNU General Public License for more details. |
14 | | * |
15 | | * You should have received a copy of the GNU General Public License |
16 | | * along with this program; if not, write to the Free Software |
17 | | * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
18 | | * |
19 | | */ |
20 | | |
21 | | /* This file was automatically generated --- DO NOT EDIT */ |
22 | | /* Generated on Sun Sep 8 06:41:29 UTC 2024 */ |
23 | | |
24 | | #include "rdft/codelet-rdft.h" |
25 | | |
26 | | #if defined(ARCH_PREFERS_FMA) || defined(ISA_EXTENSION_PREFERS_FMA) |
27 | | |
28 | | /* Generated by: ../../../genfft/gen_r2cf.native -fma -compact -variables 4 -pipeline-latency 4 -n 3 -name r2cfII_3 -dft-II -include rdft/scalar/r2cfII.h */ |
29 | | |
30 | | /* |
31 | | * This function contains 4 FP additions, 2 FP multiplications, |
32 | | * (or, 3 additions, 1 multiplications, 1 fused multiply/add), |
33 | | * 7 stack variables, 2 constants, and 6 memory accesses |
34 | | */ |
35 | | #include "rdft/scalar/r2cfII.h" |
36 | | |
37 | | static void r2cfII_3(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs) |
38 | | { |
39 | | DK(KP500000000, +0.500000000000000000000000000000000000000000000); |
40 | | DK(KP866025403, +0.866025403784438646763723170752936183471402627); |
41 | | { |
42 | | INT i; |
43 | | for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(12, rs), MAKE_VOLATILE_STRIDE(12, csr), MAKE_VOLATILE_STRIDE(12, csi)) { |
44 | | E T3, T1, T2, T4; |
45 | | T3 = R0[0]; |
46 | | T1 = R1[0]; |
47 | | T2 = R0[WS(rs, 1)]; |
48 | | T4 = T2 - T1; |
49 | | Ci[0] = -(KP866025403 * (T1 + T2)); |
50 | | Cr[0] = FNMS(KP500000000, T4, T3); |
51 | | Cr[WS(csr, 1)] = T3 + T4; |
52 | | } |
53 | | } |
54 | | } |
55 | | |
56 | | static const kr2c_desc desc = { 3, "r2cfII_3", { 3, 1, 1, 0 }, &GENUS }; |
57 | | |
58 | | void X(codelet_r2cfII_3) (planner *p) { X(kr2c_register) (p, r2cfII_3, &desc); |
59 | | } |
60 | | |
61 | | #else |
62 | | |
63 | | /* Generated by: ../../../genfft/gen_r2cf.native -compact -variables 4 -pipeline-latency 4 -n 3 -name r2cfII_3 -dft-II -include rdft/scalar/r2cfII.h */ |
64 | | |
65 | | /* |
66 | | * This function contains 4 FP additions, 2 FP multiplications, |
67 | | * (or, 3 additions, 1 multiplications, 1 fused multiply/add), |
68 | | * 7 stack variables, 2 constants, and 6 memory accesses |
69 | | */ |
70 | | #include "rdft/scalar/r2cfII.h" |
71 | | |
72 | | static void r2cfII_3(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs) |
73 | 0 | { |
74 | 0 | DK(KP500000000, +0.500000000000000000000000000000000000000000000); |
75 | 0 | DK(KP866025403, +0.866025403784438646763723170752936183471402627); |
76 | 0 | { |
77 | 0 | INT i; |
78 | 0 | for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(12, rs), MAKE_VOLATILE_STRIDE(12, csr), MAKE_VOLATILE_STRIDE(12, csi)) { |
79 | 0 | E T1, T2, T3, T4; |
80 | 0 | T1 = R0[0]; |
81 | 0 | T2 = R1[0]; |
82 | 0 | T3 = R0[WS(rs, 1)]; |
83 | 0 | T4 = T2 - T3; |
84 | 0 | Cr[WS(csr, 1)] = T1 - T4; |
85 | 0 | Ci[0] = -(KP866025403 * (T2 + T3)); |
86 | 0 | Cr[0] = FMA(KP500000000, T4, T1); |
87 | 0 | } |
88 | 0 | } |
89 | 0 | } |
90 | | |
91 | | static const kr2c_desc desc = { 3, "r2cfII_3", { 3, 1, 1, 0 }, &GENUS }; |
92 | | |
93 | 1 | void X(codelet_r2cfII_3) (planner *p) { X(kr2c_register) (p, r2cfII_3, &desc); |
94 | 1 | } |
95 | | |
96 | | #endif |