File tree 4 files changed +1134
-1134
lines changed
4 files changed +1134
-1134
lines changed Original file line number Diff line number Diff line change @@ -23,40 +23,40 @@ minfft_mkaux_dft_1d (int N) {
23
23
if (N <=0 || N & (N - 1 ))
24
24
// error if N is negative or not a power of two
25
25
return NULL ;
26
- a = malloc (sizeof (minfft_aux ));
26
+ a = malloc (sizeof (minfft_aux ));
27
27
if (a == NULL )
28
28
goto err ;
29
- a -> N = N ;
29
+ a -> N = N ;
30
30
if (N >=16 ) {
31
- a -> t = malloc (N * sizeof (minfft_cmpl ));
31
+ a -> t = malloc (N * sizeof (minfft_cmpl ));
32
32
if (a -> t == NULL )
33
33
goto err ;
34
- a -> e = malloc (N * sizeof (minfft_cmpl ));
34
+ a -> e = malloc (N * sizeof (minfft_cmpl ));
35
35
if (a -> e == NULL )
36
36
goto err ;
37
- e = a -> e ;
37
+ e = a -> e ;
38
38
while (N >=16 ) {
39
39
for (n = 0 ; n < N /4 ; n += 2 ) {
40
- e1 = exp (-2 * pi * I * n /N );
41
- e1n = exp (-2 * pi * I * (n + 1 )/N );
42
- e3 = exp (-2 * pi * I * 3 * n /N );
43
- e3n = exp (-2 * pi * I * 3 * (n + 1 )/N );
44
- * e ++ = creal (e1 );
45
- * e ++ = creal (e1n );
46
- * e ++ = creal (e3 );
47
- * e ++ = creal (e3n );
48
- * e ++ = cimag (e1 );
49
- * e ++ = cimag (e1n );
50
- * e ++ = cimag (e3 );
51
- * e ++ = cimag (e3n );
40
+ e1 = exp (-2 * pi * I * n /N );
41
+ e1n = exp (-2 * pi * I * (n + 1 )/N );
42
+ e3 = exp (-2 * pi * I * 3 * n /N );
43
+ e3n = exp (-2 * pi * I * 3 * (n + 1 )/N );
44
+ * e ++ = creal (e1 );
45
+ * e ++ = creal (e1n );
46
+ * e ++ = creal (e3 );
47
+ * e ++ = creal (e3n );
48
+ * e ++ = cimag (e1 );
49
+ * e ++ = cimag (e1n );
50
+ * e ++ = cimag (e3 );
51
+ * e ++ = cimag (e3n );
52
52
}
53
53
N /= 2 ;
54
54
}
55
55
} else {
56
- a -> t = NULL ;
57
- a -> e = NULL ;
56
+ a -> t = NULL ;
57
+ a -> e = NULL ;
58
58
}
59
- a -> sub1 = a -> sub2 = NULL ;
59
+ a -> sub1 = a -> sub2 = NULL ;
60
60
return a ;
61
61
err : // memory allocation error
62
62
minfft_free_aux (a );
Original file line number Diff line number Diff line change @@ -22,32 +22,32 @@ minfft_mkaux_dft_1d (int N) {
22
22
if (N <=0 || N & (N - 1 ))
23
23
// error if N is negative or not a power of two
24
24
return NULL ;
25
- a = malloc (sizeof (minfft_aux ));
25
+ a = malloc (sizeof (minfft_aux ));
26
26
if (a == NULL )
27
27
goto err ;
28
- a -> N = N ;
28
+ a -> N = N ;
29
29
if (N >=16 ) {
30
- a -> t = malloc (N * sizeof (minfft_cmpl ));
30
+ a -> t = malloc (N * sizeof (minfft_cmpl ));
31
31
if (a -> t == NULL )
32
32
goto err ;
33
- a -> e = malloc (N * sizeof (minfft_cmpl ));
33
+ a -> e = malloc (N * sizeof (minfft_cmpl ));
34
34
if (a -> e == NULL )
35
35
goto err ;
36
- e = a -> e ;
36
+ e = a -> e ;
37
37
while (N >=16 ) {
38
38
for (n = 0 ; n < N /4 ; n += 2 ) {
39
- * e ++ = exp (-2 * pi * I * n /N );
40
- * e ++ = exp (-2 * pi * I * (n + 1 )/N );
41
- * e ++ = exp (-2 * pi * I * 3 * n /N );
42
- * e ++ = exp (-2 * pi * I * 3 * (n + 1 )/N );
39
+ * e ++ = exp (-2 * pi * I * n /N );
40
+ * e ++ = exp (-2 * pi * I * (n + 1 )/N );
41
+ * e ++ = exp (-2 * pi * I * 3 * n /N );
42
+ * e ++ = exp (-2 * pi * I * 3 * (n + 1 )/N );
43
43
}
44
44
N /= 2 ;
45
45
}
46
46
} else {
47
- a -> t = NULL ;
48
- a -> e = NULL ;
47
+ a -> t = NULL ;
48
+ a -> e = NULL ;
49
49
}
50
- a -> sub1 = a -> sub2 = NULL ;
50
+ a -> sub1 = a -> sub2 = NULL ;
51
51
return a ;
52
52
err : // memory allocation error
53
53
minfft_free_aux (a );
You can’t perform that action at this time.
0 commit comments