@@ -7,13 +7,13 @@ public static class DoubleArrayFT
77{
88 private static ( double Sin , double Cos ) [ ] GetCoefficients ( int N , bool IsInverse = false )
99 {
10- var w = new ( double Sin , double Cos ) [ N ] ;
11- var w0 = IsInverse ? - Consts . pi2 / N : Consts . pi2 / N ;
10+ var w = new ( double Sin , double Cos ) [ N ] ;
11+ var w0 = IsInverse ? - Consts . pi2 / N : Consts . pi2 / N ;
1212 var arg = 0.0 ;
1313 for ( var i = 0 ; i < N ; i ++ )
1414 {
15- w [ i ] = ( Math . Sin ( arg ) , Math . Cos ( arg ) ) ;
16- arg += w0 ;
15+ w [ i ] = ( Math . Sin ( arg ) , Math . Cos ( arg ) ) ;
16+ arg += w0 ;
1717 }
1818 return w ;
1919 }
@@ -29,17 +29,17 @@ public static Spectrum GetFourierTransformation(this double[] s, bool IsInverse
2929
3030 return m =>
3131 {
32- var P = 0.0 ;
33- var Q = 0.0 ;
34- for ( var n = 0 ; n < N ; n ++ )
32+ var p = 0.0 ;
33+ var q = 0.0 ;
34+ for ( var n = 0 ; n < N ; n ++ )
3535 {
3636 var val = s [ n ] ;
3737 var ( sin , cos ) = w [ n * m % N ] ;
38- P += val * cos ;
39- Q += val * sin ;
38+ p += val * cos ;
39+ q += val * sin ;
4040 }
4141
42- return new ( P / N , Q / N ) ;
42+ return new ( p / N , q / N ) ;
4343 } ;
4444 }
4545
@@ -54,17 +54,17 @@ public static Spectrum GetFourierTransformation(this Complex[] s, bool IsInverse
5454
5555 return m =>
5656 {
57- var P = 0.0 ;
58- var Q = 0.0 ;
59- for ( var n = 0 ; n < N ; n ++ )
57+ var p = 0.0 ;
58+ var q = 0.0 ;
59+ for ( var n = 0 ; n < N ; n ++ )
6060 {
6161 var ( re , im ) = s [ n ] ;
6262 var ( sin , cos ) = w [ n * m % N ] ;
63- P += re * cos - im * sin ;
64- Q += im * cos + re * sin ;
63+ p += re * cos - im * sin ;
64+ q += im * cos + re * sin ;
6565 }
6666
67- return new ( P / N , Q / N ) ;
67+ return new ( p / N , q / N ) ;
6868 } ;
6969 }
7070}
0 commit comments