Fix ooura

Author: kunitoki

modules/yup_dsp/fft/yup_OouraFFT8g.cpp

@@ -313,11 +313,13 @@ Appendix :
 namespace yup
 {
 
+void makewt(int nw, int *ip, float *w);
+void cftfsub(int n, float *a, int *ip, int nw, float *w);
+void cftbsub(int n, float *a, int *ip, int nw, float *w);
+
+
 void cdft(int n, int isgn, float *a, int *ip, float *w)
 {
-  void makewt(int nw, int *ip, float *w);
-  void cftfsub(int n, float *a, int *ip, int nw, float *w);
-  void cftbsub(int n, float *a, int *ip, int nw, float *w);
   int nw;
 
   nw = ip[0];
@@ -336,15 +338,15 @@ void cdft(int n, int isgn, float *a, int *ip, float *w)
   }
 }
 
+void makewt(int nw, int *ip, float *w);
+void makect(int nc, int *ip, float *c);
+void cftfsub(int n, float *a, int *ip, int nw, float *w);
+void cftbsub(int n, float *a, int *ip, int nw, float *w);
+void rftfsub(int n, float *a, int nc, float *c);
+void rftbsub(int n, float *a, int nc, float *c);
 
 void rdft(int n, int isgn, float *a, int *ip, float *w)
 {
-  void makewt(int nw, int *ip, float *w);
-  void makect(int nc, int *ip, float *c);
-  void cftfsub(int n, float *a, int *ip, int nw, float *w);
-  void cftbsub(int n, float *a, int *ip, int nw, float *w);
-  void rftfsub(int n, float *a, int nc, float *c);
-  void rftbsub(int n, float *a, int nc, float *c);
   int nw, nc;
   float xi;
 
@@ -391,16 +393,16 @@ void rdft(int n, int isgn, float *a, int *ip, float *w)
   }
 }
 
+void makewt(int nw, int *ip, float *w);
+void makect(int nc, int *ip, float *c);
+void cftfsub(int n, float *a, int *ip, int nw, float *w);
+void cftbsub(int n, float *a, int *ip, int nw, float *w);
+void rftfsub(int n, float *a, int nc, float *c);
+void rftbsub(int n, float *a, int nc, float *c);
+void dctsub(int n, float *a, int nc, float *c);
 
 void ddct(int n, int isgn, float *a, int *ip, float *w)
 {
-  void makewt(int nw, int *ip, float *w);
-  void makect(int nc, int *ip, float *c);
-  void cftfsub(int n, float *a, int *ip, int nw, float *w);
-  void cftbsub(int n, float *a, int *ip, int nw, float *w);
-  void rftfsub(int n, float *a, int nc, float *c);
-  void rftbsub(int n, float *a, int nc, float *c);
-  void dctsub(int n, float *a, int nc, float *c);
   int j, nw, nc;
   float xr;
 
@@ -459,16 +461,16 @@ void ddct(int n, int isgn, float *a, int *ip, float *w)
   }
 }
 
+void makewt(int nw, int *ip, float *w);
+void makect(int nc, int *ip, float *c);
+void cftfsub(int n, float *a, int *ip, int nw, float *w);
+void cftbsub(int n, float *a, int *ip, int nw, float *w);
+void rftfsub(int n, float *a, int nc, float *c);
+void rftbsub(int n, float *a, int nc, float *c);
+void dstsub(int n, float *a, int nc, float *c);
 
 void ddst(int n, int isgn, float *a, int *ip, float *w)
 {
-  void makewt(int nw, int *ip, float *w);
-  void makect(int nc, int *ip, float *c);
-  void cftfsub(int n, float *a, int *ip, int nw, float *w);
-  void cftbsub(int n, float *a, int *ip, int nw, float *w);
-  void rftfsub(int n, float *a, int nc, float *c);
-  void rftbsub(int n, float *a, int nc, float *c);
-  void dstsub(int n, float *a, int nc, float *c);
   int j, nw, nc;
   float xr;
 
@@ -527,14 +529,14 @@ void ddst(int n, int isgn, float *a, int *ip, float *w)
   }
 }
 
+void makewt(int nw, int *ip, float *w);
+void makect(int nc, int *ip, float *c);
+void cftfsub(int n, float *a, int *ip, int nw, float *w);
+void rftfsub(int n, float *a, int nc, float *c);
+void dctsub(int n, float *a, int nc, float *c);
 
 void dfct(int n, float *a, float *t, int *ip, float *w)
 {
-  void makewt(int nw, int *ip, float *w);
-  void makect(int nc, int *ip, float *c);
-  void cftfsub(int n, float *a, int *ip, int nw, float *w);
-  void rftfsub(int n, float *a, int nc, float *c);
-  void dctsub(int n, float *a, int nc, float *c);
   int j, k, l, m, mh, nw, nc;
   float xr, xi, yr, yi;
 
@@ -636,14 +638,14 @@ void dfct(int n, float *a, float *t, int *ip, float *w)
   }
 }
 
+void makewt(int nw, int *ip, float *w);
+void makect(int nc, int *ip, float *c);
+void cftfsub(int n, float *a, int *ip, int nw, float *w);
+void rftfsub(int n, float *a, int nc, float *c);
+void dstsub(int n, float *a, int nc, float *c);
 
 void dfst(int n, float *a, float *t, int *ip, float *w)
 {
-  void makewt(int nw, int *ip, float *w);
-  void makect(int nc, int *ip, float *c);
-  void cftfsub(int n, float *a, int *ip, int nw, float *w);
-  void rftfsub(int n, float *a, int nc, float *c);
-  void dstsub(int n, float *a, int nc, float *c);
   int j, k, l, m, mh, nw, nc;
   float xr, xi, yr, yi;
 
@@ -737,10 +739,10 @@ void dfst(int n, float *a, float *t, int *ip, float *w)
 
 /* -------- initializing routines -------- */
 
+void makeipt(int nw, int *ip);
 
 void makewt(int nw, int *ip, float *w)
 {
-  void makeipt(int nw, int *ip);
   int j, nwh, nw0, nw1;
   float delta, wn4r, wk1r, wk1i, wk3r, wk3i;
 
@@ -916,24 +918,23 @@ void makect(int nc, int *ip, float *c)
     }
 #endif /* USE_CDFT_WINTHREADS */
 
-
-void cftfsub(int n, float *a, int *ip, int nw, float *w)
-{
-  void bitrv2(int n, int *ip, float *a);
-  void bitrv216(float *a);
-  void bitrv208(float *a);
-  void cftf1st(int n, float *a, float *w);
-  void cftrec4(int n, float *a, int nw, float *w);
-  void cftleaf(int n, int isplt, float *a, int nw, float *w);
-  void cftfx41(int n, float *a, int nw, float *w);
-  void cftf161(float *a, float *w);
-  void cftf081(float *a, float *w);
-  void cftf040(float *a);
-  void cftx020(float *a);
+void bitrv2(int n, int *ip, float *a);
+void bitrv216(float *a);
+void bitrv208(float *a);
+void cftf1st(int n, float *a, float *w);
+void cftrec4(int n, float *a, int nw, float *w);
+void cftleaf(int n, int isplt, float *a, int nw, float *w);
+void cftfx41(int n, float *a, int nw, float *w);
+void cftf161(float *a, float *w);
+void cftf081(float *a, float *w);
+void cftf040(float *a);
+void cftx020(float *a);
 #ifdef USE_CDFT_THREADS
-  void cftrec4_th(int n, float *a, int nw, float *w);
+void cftrec4_th(int n, float *a, int nw, float *w);
 #endif /* USE_CDFT_THREADS */
 
+void cftfsub(int n, float *a, int *ip, int nw, float *w)
+{
   if (n > 8)
   {
     if (n > 32)
@@ -981,24 +982,23 @@ void cftfsub(int n, float *a, int *ip, int nw, float *w)
   }
 }
 
-
-void cftbsub(int n, float *a, int *ip, int nw, float *w)
-{
-  void bitrv2conj(int n, int *ip, float *a);
-  void bitrv216neg(float *a);
-  void bitrv208neg(float *a);
-  void cftb1st(int n, float *a, float *w);
-  void cftrec4(int n, float *a, int nw, float *w);
-  void cftleaf(int n, int isplt, float *a, int nw, float *w);
-  void cftfx41(int n, float *a, int nw, float *w);
-  void cftf161(float *a, float *w);
-  void cftf081(float *a, float *w);
-  void cftb040(float *a);
-  void cftx020(float *a);
+void bitrv2conj(int n, int *ip, float *a);
+void bitrv216neg(float *a);
+void bitrv208neg(float *a);
+void cftb1st(int n, float *a, float *w);
+void cftrec4(int n, float *a, int nw, float *w);
+void cftleaf(int n, int isplt, float *a, int nw, float *w);
+void cftfx41(int n, float *a, int nw, float *w);
+void cftf161(float *a, float *w);
+void cftf081(float *a, float *w);
+void cftb040(float *a);
+void cftx020(float *a);
 #ifdef USE_CDFT_THREADS
-  void cftrec4_th(int n, float *a, int nw, float *w);
+void cftrec4_th(int n, float *a, int nw, float *w);
 #endif /* USE_CDFT_THREADS */
 
+void cftbsub(int n, float *a, int *ip, int nw, float *w)
+{
   if (n > 8)
   {
     if (n > 32)
@@ -2368,11 +2368,11 @@ struct cdft_arg_st
 };
 typedef struct cdft_arg_st cdft_arg_t;
 
+void *cftrec1_th(void *p);
+void *cftrec2_th(void *p);
 
 void cftrec4_th(int n, float *a, int nw, float *w)
 {
-  void *cftrec1_th(void *p);
-  void *cftrec2_th(void *p);
   int i, idiv4, m, nthread;
   cdft_thread_t th[4];
   cdft_arg_t ag[4];
@@ -2408,12 +2408,12 @@ void cftrec4_th(int n, float *a, int nw, float *w)
   }
 }
 
+int cfttree(int n, int j, int k, float *a, int nw, float *w);
+void cftleaf(int n, int isplt, float *a, int nw, float *w);
+void cftmdl1(int n, float *a, float *w);
 
 void *cftrec1_th(void *p)
 {
-  int cfttree(int n, int j, int k, float *a, int nw, float *w);
-  void cftleaf(int n, int isplt, float *a, int nw, float *w);
-  void cftmdl1(int n, float *a, float *w);
   int isplt, j, k, m, n, n0, nw;
   float *a, *w;
 
@@ -2439,12 +2439,12 @@ void *cftrec1_th(void *p)
   return (void *)0;
 }
 
+int cfttree(int n, int j, int k, float *a, int nw, float *w);
+void cftleaf(int n, int isplt, float *a, int nw, float *w);
+void cftmdl2(int n, float *a, float *w);
 
 void *cftrec2_th(void *p)
 {
-  int cfttree(int n, int j, int k, float *a, int nw, float *w);
-  void cftleaf(int n, int isplt, float *a, int nw, float *w);
-  void cftmdl2(int n, float *a, float *w);
   int isplt, j, k, m, n, n0, nw;
   float *a, *w;
 
@@ -2473,12 +2473,12 @@ void *cftrec2_th(void *p)
 }
 #endif /* USE_CDFT_THREADS */
 
+int cfttree(int n, int j, int k, float *a, int nw, float *w);
+void cftleaf(int n, int isplt, float *a, int nw, float *w);
+void cftmdl1(int n, float *a, float *w);
 
 void cftrec4(int n, float *a, int nw, float *w)
 {
-  int cfttree(int n, int j, int k, float *a, int nw, float *w);
-  void cftleaf(int n, int isplt, float *a, int nw, float *w);
-  void cftmdl1(int n, float *a, float *w);
   int isplt, j, k, m;
 
   m = n;
@@ -2497,11 +2497,11 @@ void cftrec4(int n, float *a, int nw, float *w)
   }
 }
 
+void cftmdl1(int n, float *a, float *w);
+void cftmdl2(int n, float *a, float *w);
 
 int cfttree(int n, int j, int k, float *a, int nw, float *w)
 {
-  void cftmdl1(int n, float *a, float *w);
-  void cftmdl2(int n, float *a, float *w);
   int i, isplt, m;
 
   if ((k & 3) != 0)
@@ -2544,16 +2544,15 @@ int cfttree(int n, int j, int k, float *a, int nw, float *w)
   return isplt;
 }
 
+void cftmdl1(int n, float *a, float *w);
+void cftmdl2(int n, float *a, float *w);
+void cftf161(float *a, float *w);
+void cftf162(float *a, float *w);
+void cftf081(float *a, float *w);
+void cftf082(float *a, float *w);
 
 void cftleaf(int n, int isplt, float *a, int nw, float *w)
 {
-  void cftmdl1(int n, float *a, float *w);
-  void cftmdl2(int n, float *a, float *w);
-  void cftf161(float *a, float *w);
-  void cftf162(float *a, float *w);
-  void cftf081(float *a, float *w);
-  void cftf082(float *a, float *w);
-
   if (n == 512)
   {
     cftmdl1(128, a, &w[nw - 64]);
@@ -2864,14 +2863,13 @@ void cftmdl2(int n, float *a, float *w)
   a[j3 + 1] = y0i + y2i;
 }
 
+void cftf161(float *a, float *w);
+void cftf162(float *a, float *w);
+void cftf081(float *a, float *w);
+void cftf082(float *a, float *w);
 
 void cftfx41(int n, float *a, int nw, float *w)
 {
-  void cftf161(float *a, float *w);
-  void cftf162(float *a, float *w);
-  void cftf081(float *a, float *w);
-  void cftf082(float *a, float *w);
-
   if (n == 128)
   {
     cftf161(a, &w[nw - 8]);