DISTRHO
/
plugin-examples
mirror of git://github.com/DISTRHO/plugin-examples
1
0
Fork 0
Code Releases 0 Activity
DPF Plugin examples
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
19 Commits
1 Branch
132KB
Tree: 5794133242
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '5794133242'
${ noResults }
plugin-examples/plugins/openav-roomy/dsp/dsp_reverb.hxx

512 lines
21KB
Raw Blame History 

  1. /*

  2.  * Author: Harry van Haaren 2013

  3.  *         harryhaaren@gmail.com

  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,

  18.  * MA 02110-1301, USA.

  19.  */

  20. 

  21. /*

  22.   The DSP code in this file was generated using FAUST.

  23.   The FAUST code used is under the following licenses:

  24. 

  25.   effect.lib/name       Faust Audio Effect Library

  26.   effect.lib/author     Julius O. Smith (jos at ccrma.stanford.edu)

  27.   effect.lib/copyright  Julius O. Smith III

  28.   effect.lib/version    1.33

  29.   effect.lib/license    STK-4.3

  30.   effect.lib/reference  https://ccrma.stanford.edu/realsimple/faust_strings/

  31. 

  32.   filter.lib/author     Julius O. Smith (jos at ccrma.stanford.edu)

  33.   filter.lib/name       Faust Filter Library

  34.   filter.lib/copyright  Julius O. Smith III

  35.   filter.lib/version    1.29

  36.   filter.lib/license    STK-4.3

  37.   filter.lib/reference  https://ccrma.stanford.edu/~jos/filters/

  38. 

  39.   music.lib/name        Music Library

  40.   music.lib/author      GRAME

  41.   music.lib/copyright   GRAME

  42.   music.lib/version     1.0

  43.   music.lib/license     LGPL

  44. 

  45.   math.lib/name         Math Library

  46.   math.lib/author       GRAME

  47.   math.lib/copyright    GRAME

  48.   math.lib/version      1.0

  49.   math.lib/license      LGPL

  50. */

  51. 

  52. 

  53. #ifndef OPENAV_DSP_REVERB_H

  54. #define OPENAV_DSP_REVERB_H

  55. 

  56. #include <math.h>

  57. #include <cmath>

  58. #include <cstdlib>

  59. 

  60. template <class T> const T& max (const T& a, const T& b) {

  61.   return (a<b)?b:a;

  62. }

  63. 

  64. template <class T> const T& min (const T& a, const T& b) {

  65.   return !(b<a)?a:b;

  66. }

  67. 

  68. template <int N> inline float faustpower(float x) 		{ return powf(x,N); }

  69. template <int N> inline double faustpower(double x) 	{ return pow(x,N); }

  70. template <int N> inline int faustpower(int x) 			{ return faustpower<N/2>(x) * faustpower<N-N/2>(x); }

  71. template <> 	 inline int faustpower<0>(int x) 		{ return 1; }

  72. template <> 	 inline int faustpower<1>(int x) 		{ return x; }

  73. 

  74. class Reverb // : Effect

  75. {

  76.   public:

  77.     Reverb(int sr)

  78.     {

  79.       init( sr );

  80.     }

  81. 

  82.     int getNumInputs() { return 2; }

  83.     int getNumOutputs(){ return 2; }

  84. 

  85.     void damping(float d)

  86.     {

  87.       if( d > 1.0 ) d = 1.0f;

  88.       if( d < 0.0 ) d = 0.0f;

  89. 

  90.       fslider1 = (1-d) * 18500 + 1500.f;

  91.     }

  92. 

  93.     void rt60(float rt)

  94.     {

  95.       if( rt > 1.0 ) rt = 1.0f;

  96.       if( rt < 0.0 ) rt = 0.0f;

  97. 

  98.       fslider0 = 1 + rt * 5;

  99.     }

  100. 

  101.     void dryWet(float dw)

  102.     {

  103.       if( dw > 1.0 ) dw = 1.0f;

  104.       if( dw < 0.0 ) dw = 0.0f;

  105.       _dryWet = dw;

  106.     }

  107. 

  108.     void process (int count, const float** input, float** output)

  109.     {

  110.       float 	fSlow0 = fslider0;

  111.       float 	fSlow1 = expf((fConst2 / fSlow0));

  112.       float 	fSlow2 = faustpower<2>(fSlow1);

  113.       float 	fSlow3 = (1.0f - fSlow2);

  114.       float 	fSlow4 = cosf((fConst3 * fslider1));

  115.       float 	fSlow5 = (1.0f - (fSlow4 * fSlow2));

  116.       float 	fSlow6 = sqrtf(max(0.f, ((faustpower<2>(fSlow5) / faustpower<2>(fSlow3)) - 1.0f)));

  117.       float 	fSlow7 = (fSlow5 / fSlow3);

  118.       float 	fSlow8 = (fSlow7 - fSlow6);

  119.       float 	fSlow9 = (fSlow0 - 0.5f);

  120.       float 	fSlow10 = ((expf((fConst2 / fSlow9)) / fSlow1) - 1);

  121.       float 	fSlow11 = (fSlow1 * ((1.0f + fSlow6) - fSlow7));

  122.       float 	fSlow12 = expf((fConst12 / fSlow0));

  123.       float 	fSlow13 = faustpower<2>(fSlow12);

  124.       float 	fSlow14 = (1.0f - fSlow13);

  125.       float 	fSlow15 = (1.0f - (fSlow4 * fSlow13));

  126.       float 	fSlow16 = sqrtf(max(0.f, ((faustpower<2>(fSlow15) / faustpower<2>(fSlow14)) - 1.0f)));

  127.       float 	fSlow17 = (fSlow15 / fSlow14);

  128.       float 	fSlow18 = (fSlow17 - fSlow16);

  129.       float 	fSlow19 = ((expf((fConst12 / fSlow9)) / fSlow12) - 1);

  130.       float 	fSlow20 = (fSlow12 * ((1.0f + fSlow16) - fSlow17));

  131.       float 	fSlow21 = expf((fConst17 / fSlow0));

  132.       float 	fSlow22 = faustpower<2>(fSlow21);

  133.       float 	fSlow23 = (1.0f - fSlow22);

  134.       float 	fSlow24 = (1.0f - (fSlow4 * fSlow22));

  135.       float 	fSlow25 = sqrtf(max(0.f, ((faustpower<2>(fSlow24) / faustpower<2>(fSlow23)) - 1.0f)));

  136.       float 	fSlow26 = (fSlow24 / fSlow23);

  137.       float 	fSlow27 = (fSlow26 - fSlow25);

  138.       float 	fSlow28 = ((expf((fConst17 / fSlow9)) / fSlow21) - 1);

  139.       float 	fSlow29 = (fSlow21 * ((1.0f + fSlow25) - fSlow26));

  140.       float 	fSlow30 = expf((fConst22 / fSlow0));

  141.       float 	fSlow31 = faustpower<2>(fSlow30);

  142.       float 	fSlow32 = (1.0f - fSlow31);

  143.       float 	fSlow33 = (1.0f - (fSlow4 * fSlow31));

  144.       float 	fSlow34 = sqrtf(max(0.f, ((faustpower<2>(fSlow33) / faustpower<2>(fSlow32)) - 1.0f)));

  145.       float 	fSlow35 = (fSlow33 / fSlow32);

  146.       float 	fSlow36 = (fSlow35 - fSlow34);

  147.       float 	fSlow37 = ((expf((fConst22 / fSlow9)) / fSlow30) - 1);

  148.       float 	fSlow38 = (fSlow30 * ((1.0f + fSlow34) - fSlow35));

  149.       float 	fSlow39 = expf((fConst27 / fSlow0));

  150.       float 	fSlow40 = faustpower<2>(fSlow39);

  151.       float 	fSlow41 = (1.0f - fSlow40);

  152.       float 	fSlow42 = (1.0f - (fSlow4 * fSlow40));

  153.       float 	fSlow43 = sqrtf(max(0.f, ((faustpower<2>(fSlow42) / faustpower<2>(fSlow41)) - 1.0f)));

  154.       float 	fSlow44 = (fSlow42 / fSlow41);

  155.       float 	fSlow45 = (fSlow44 - fSlow43);

  156.       float 	fSlow46 = ((expf((fConst27 / fSlow9)) / fSlow39) - 1);

  157.       float 	fSlow47 = (fSlow39 * ((1.0f + fSlow43) - fSlow44));

  158.       float 	fSlow48 = expf((fConst32 / fSlow0));

  159.       float 	fSlow49 = faustpower<2>(fSlow48);

  160.       float 	fSlow50 = (1.0f - fSlow49);

  161.       float 	fSlow51 = (1.0f - (fSlow4 * fSlow49));

  162.       float 	fSlow52 = sqrtf(max(0.f, ((faustpower<2>(fSlow51) / faustpower<2>(fSlow50)) - 1.0f)));

  163.       float 	fSlow53 = (fSlow51 / fSlow50);

  164.       float 	fSlow54 = (fSlow53 - fSlow52);

  165.       float 	fSlow55 = ((expf((fConst32 / fSlow9)) / fSlow48) - 1);

  166.       float 	fSlow56 = (fSlow48 * ((1.0f + fSlow52) - fSlow53));

  167.       float 	fSlow57 = expf((fConst37 / fSlow0));

  168.       float 	fSlow58 = faustpower<2>(fSlow57);

  169.       float 	fSlow59 = (1.0f - fSlow58);

  170.       float 	fSlow60 = (1.0f - (fSlow4 * fSlow58));

  171.       float 	fSlow61 = sqrtf(max(0.f, ((faustpower<2>(fSlow60) / faustpower<2>(fSlow59)) - 1.0f)));

  172.       float 	fSlow62 = (fSlow60 / fSlow59);

  173.       float 	fSlow63 = (fSlow62 - fSlow61);

  174.       float 	fSlow64 = ((expf((fConst37 / fSlow9)) / fSlow57) - 1);

  175.       float 	fSlow65 = (fSlow57 * ((1.0f + fSlow61) - fSlow62));

  176.       float 	fSlow66 = expf((fConst42 / fSlow0));

  177.       float 	fSlow67 = faustpower<2>(fSlow66);

  178.       float 	fSlow68 = (1.0f - fSlow67);

  179.       float 	fSlow69 = (1.0f - (fSlow67 * fSlow4));

  180.       float 	fSlow70 = sqrtf(max(0.f, ((faustpower<2>(fSlow69) / faustpower<2>(fSlow68)) - 1.0f)));

  181.       float 	fSlow71 = (fSlow69 / fSlow68);

  182.       float 	fSlow72 = (fSlow71 - fSlow70);

  183.       float 	fSlow73 = ((expf((fConst42 / fSlow9)) / fSlow66) - 1);

  184.       float 	fSlow74 = (fSlow66 * ((1.0f + fSlow70) - fSlow71));

  185.       const float* input0 = input[0];

  186.       const float* input1 = input[1];

  187.       float* output0 = output[0];

  188.       float* output1 = output[1];

  189.       for (int i=0; i<count; i++) {

  190.         fRec11[0] = ((fConst7 * (fRec4[1] + fRec4[2])) + (fConst6 * fRec11[1]));

  191.         fRec10[0] = ((fSlow11 * (fRec4[1] + (fSlow10 * fRec11[0]))) + (fSlow8 * fRec10[1]));

  192.         fVec0[IOTA&8191] = (1e-20f + (0.35355339059327373f * fRec10[0]));

  193.         float fTemp0 = (0.3f * (float)input0[i]);

  194.         float fTemp1 = ((fTemp0 + fVec0[(IOTA-iConst9)&8191]) - (0.6f * fRec8[1]));

  195.         fVec1[IOTA&2047] = fTemp1;

  196.         fRec8[0] = fVec1[(IOTA-iConst10)&2047];

  197.         float 	fRec9 = (0.6f * fVec1[IOTA&2047]);

  198.         fRec15[0] = ((fConst7 * (fRec0[1] + fRec0[2])) + (fConst6 * fRec15[1]));

  199.         fRec14[0] = ((fSlow20 * (fRec0[1] + (fSlow19 * fRec15[0]))) + (fSlow18 * fRec14[1]));

  200.         fVec2[IOTA&8191] = (1e-20f + (0.35355339059327373f * fRec14[0]));

  201.         float fTemp2 = ((fTemp0 + fVec2[(IOTA-iConst14)&8191]) - (0.6f * fRec12[1]));

  202.         fVec3[IOTA&1023] = fTemp2;

  203.         fRec12[0] = fVec3[(IOTA-iConst15)&1023];

  204.         float 	fRec13 = (0.6f * fVec3[IOTA&1023]);

  205.         float fTemp3 = (fRec13 + fRec9);

  206.         fRec19[0] = ((fConst7 * (fRec2[1] + fRec2[2])) + (fConst6 * fRec19[1]));

  207.         fRec18[0] = ((fSlow29 * (fRec2[1] + (fSlow28 * fRec19[0]))) + (fSlow27 * fRec18[1]));

  208.         fVec4[IOTA&8191] = (1e-20f + (0.35355339059327373f * fRec18[0]));

  209.         float fTemp4 = (fVec4[(IOTA-iConst19)&8191] - (fTemp0 + (0.6f * fRec16[1])));

  210.         fVec5[IOTA&2047] = fTemp4;

  211.         fRec16[0] = fVec5[(IOTA-iConst20)&2047];

  212.         float 	fRec17 = (0.6f * fVec5[IOTA&2047]);

  213.         fRec23[0] = ((fConst7 * (fRec6[1] + fRec6[2])) + (fConst6 * fRec23[1]));

  214.         fRec22[0] = ((fSlow38 * (fRec6[1] + (fSlow37 * fRec23[0]))) + (fSlow36 * fRec22[1]));

  215.         fVec6[IOTA&8191] = (1e-20f + (0.35355339059327373f * fRec22[0]));

  216.         float fTemp5 = (fVec6[(IOTA-iConst24)&8191] - (fTemp0 + (0.6f * fRec20[1])));

  217.         fVec7[IOTA&1023] = fTemp5;

  218.         fRec20[0] = fVec7[(IOTA-iConst25)&1023];

  219.         float 	fRec21 = (0.6f * fVec7[IOTA&1023]);

  220.         float fTemp6 = (fRec21 + (fRec17 + fTemp3));

  221.         fRec27[0] = ((fConst7 * (fRec1[1] + fRec1[2])) + (fConst6 * fRec27[1]));

  222.         fRec26[0] = ((fSlow47 * (fRec1[1] + (fSlow46 * fRec27[0]))) + (fSlow45 * fRec26[1]));

  223.         fVec8[IOTA&16383] = (1e-20f + (0.35355339059327373f * fRec26[0]));

  224.         float fTemp7 = (0.3f * (float)input1[i]);

  225.         float fTemp8 = ((fTemp7 + fVec8[(IOTA-iConst29)&16383]) + (0.6f * fRec24[1]));

  226.         fVec9[IOTA&2047] = fTemp8;

  227.         fRec24[0] = fVec9[(IOTA-iConst30)&2047];

  228.         float 	fRec25 = (0 - (0.6f * fVec9[IOTA&2047]));

  229.         fRec31[0] = ((fConst7 * (fRec5[1] + fRec5[2])) + (fConst6 * fRec31[1]));

  230.         fRec30[0] = ((fSlow56 * (fRec5[1] + (fSlow55 * fRec31[0]))) + (fSlow54 * fRec30[1]));

  231.         fVec10[IOTA&8191] = (1e-20f + (0.35355339059327373f * fRec30[0]));

  232.         float fTemp9 = ((fTemp7 + fVec10[(IOTA-iConst34)&8191]) + (0.6f * fRec28[1]));

  233.         fVec11[IOTA&2047] = fTemp9;

  234.         fRec28[0] = fVec11[(IOTA-iConst35)&2047];

  235.         float 	fRec29 = (0 - (0.6f * fVec11[IOTA&2047]));

  236.         fRec35[0] = ((fConst7 * (fRec3[1] + fRec3[2])) + (fConst6 * fRec35[1]));

  237.         fRec34[0] = ((fSlow65 * (fRec3[1] + (fSlow64 * fRec35[0]))) + (fSlow63 * fRec34[1]));

  238.         fVec12[IOTA&16383] = (1e-20f + (0.35355339059327373f * fRec34[0]));

  239.         float fTemp10 = ((fVec12[(IOTA-iConst39)&16383] + (0.6f * fRec32[1])) - fTemp7);

  240.         fVec13[IOTA&2047] = fTemp10;

  241.         fRec32[0] = fVec13[(IOTA-iConst40)&2047];

  242.         float 	fRec33 = (0 - (0.6f * fVec13[IOTA&2047]));

  243.         fRec39[0] = ((fConst7 * (fRec7[1] + fRec7[2])) + (fConst6 * fRec39[1]));

  244.         fRec38[0] = ((fSlow74 * (fRec7[1] + (fSlow73 * fRec39[0]))) + (fSlow72 * fRec38[1]));

  245.         fVec14[IOTA&16383] = (1e-20f + (0.35355339059327373f * fRec38[0]));

  246.         float fTemp11 = ((fVec14[(IOTA-iConst44)&16383] + (0.6f * fRec36[1])) - fTemp7);

  247.         fVec15[IOTA&1023] = fTemp11;

  248.         fRec36[0] = fVec15[(IOTA-iConst45)&1023];

  249.         float 	fRec37 = (0 - (0.6f * fVec15[IOTA&1023]));

  250.         fRec0[0] = (fRec12[1] + (fRec8[1] + (fRec16[1] + (fRec20[1] + (fRec24[1] + (fRec28[1] + (fRec32[1] + (fRec36[1] + (fRec37 + (fRec33 + (fRec29 + (fRec25 + fTemp6))))))))))));

  251.         fRec1[0] = (0 - ((fRec24[1] + (fRec28[1] + (fRec32[1] + (fRec36[1] + (fRec37 + (fRec33 + (fRec25 + fRec29))))))) - (fRec12[1] + (fRec8[1] + (fRec16[1] + (fRec20[1] + fTemp6))))));

  252.         float fTemp12 = (fRec17 + fRec21);

  253.         fRec2[0] = (0 - ((fRec16[1] + (fRec20[1] + (fRec32[1] + (fRec36[1] + (fRec37 + (fRec33 + fTemp12)))))) - (fRec12[1] + (fRec8[1] + (fRec24[1] + (fRec28[1] + (fRec29 + (fRec25 + fTemp3))))))));

  254.         fRec3[0] = (0 - ((fRec16[1] + (fRec20[1] + (fRec24[1] + (fRec28[1] + (fRec29 + (fRec25 + fTemp12)))))) - (fRec12[1] + (fRec8[1] + (fRec32[1] + (fRec36[1] + (fRec37 + (fRec33 + fTemp3))))))));

  255.         float fTemp13 = (fRec13 + fRec17);

  256.         float fTemp14 = (fRec9 + fRec21);

  257.         fRec4[0] = (0 - ((fRec8[1] + (fRec20[1] + (fRec28[1] + (fRec36[1] + (fRec37 + (fRec29 + fTemp14)))))) - (fRec12[1] + (fRec16[1] + (fRec24[1] + (fRec32[1] + (fRec33 + (fRec25 + fTemp13))))))));

  258.         fRec5[0] = (0 - ((fRec8[1] + (fRec20[1] + (fRec24[1] + (fRec32[1] + (fRec33 + (fRec25 + fTemp14)))))) - (fRec12[1] + (fRec16[1] + (fRec28[1] + (fRec36[1] + (fRec37 + (fRec29 + fTemp13))))))));

  259.         float fTemp15 = (fRec13 + fRec21);

  260.         float fTemp16 = (fRec9 + fRec17);

  261.         fRec6[0] = (0 - ((fRec8[1] + (fRec16[1] + (fRec28[1] + (fRec32[1] + (fRec33 + (fRec29 + fTemp16)))))) - (fRec12[1] + (fRec20[1] + (fRec24[1] + (fRec36[1] + (fRec37 + (fRec25 + fTemp15))))))));

  262.         fRec7[0] = (0 - ((fRec8[1] + (fRec16[1] + (fRec24[1] + (fRec36[1] + (fRec37 + (fRec25 + fTemp16)))))) - (fRec12[1] + (fRec20[1] + (fRec28[1] + (fRec32[1] + (fRec33 + (fRec29 + fTemp15))))))));

  263. 

  264. 

  265.         float reverb0 = (float)(0.37f * (fRec1[0] + fRec2[0]));

  266.         float reverb1 = (float)(0.37f * (fRec1[0] - fRec2[0]));

  267. 

  268.         output0[i] = (input0[i] * (1-_dryWet)) + (reverb0 * _dryWet );

  269.         output1[i] = (input1[i] * (1-_dryWet)) + (reverb1 * _dryWet );

  270. 

  271.         // post processing

  272.         fRec7[2] = fRec7[1]; fRec7[1] = fRec7[0];

  273.         fRec6[2] = fRec6[1]; fRec6[1] = fRec6[0];

  274.         fRec5[2] = fRec5[1]; fRec5[1] = fRec5[0];

  275.         fRec4[2] = fRec4[1]; fRec4[1] = fRec4[0];

  276.         fRec3[2] = fRec3[1]; fRec3[1] = fRec3[0];

  277.         fRec2[2] = fRec2[1]; fRec2[1] = fRec2[0];

  278.         fRec1[2] = fRec1[1]; fRec1[1] = fRec1[0];

  279.         fRec0[2] = fRec0[1]; fRec0[1] = fRec0[0];

  280.         fRec36[1] = fRec36[0];

  281.         fRec38[1] = fRec38[0];

  282.         fRec39[1] = fRec39[0];

  283.         fRec32[1] = fRec32[0];

  284.         fRec34[1] = fRec34[0];

  285.         fRec35[1] = fRec35[0];

  286.         fRec28[1] = fRec28[0];

  287.         fRec30[1] = fRec30[0];

  288.         fRec31[1] = fRec31[0];

  289.         fRec24[1] = fRec24[0];

  290.         fRec26[1] = fRec26[0];

  291.         fRec27[1] = fRec27[0];

  292.         fRec20[1] = fRec20[0];

  293.         fRec22[1] = fRec22[0];

  294.         fRec23[1] = fRec23[0];

  295.         fRec16[1] = fRec16[0];

  296.         fRec18[1] = fRec18[0];

  297.         fRec19[1] = fRec19[0];

  298.         fRec12[1] = fRec12[0];

  299.         fRec14[1] = fRec14[0];

  300.         fRec15[1] = fRec15[0];

  301.         fRec8[1] = fRec8[0];

  302.         IOTA = IOTA+1;

  303.         fRec10[1] = fRec10[0];

  304.         fRec11[1] = fRec11[0];

  305.       }

  306.     }

  307. 

  308.   private:

  309.     float _dryWet;

  310.     float fslider0;

  311.     int iConst0;

  312.     float fConst1;

  313.     float fConst2;

  314.     float fslider1;

  315.     float fConst3;

  316.     float fConst4;

  317.     float fConst5;

  318.     float fConst6;

  319.     float fConst7;

  320.     float fRec11[2];

  321.     float fRec10[2];

  322.     int   IOTA;

  323.     float fVec0[8192];

  324.     float fConst8;

  325.     int iConst9;

  326.     float fVec1[2048];

  327.     int iConst10;

  328.     float fRec8[2];

  329.     float fConst11;

  330.     float fConst12;

  331.     float fRec15[2];

  332.     float fRec14[2];

  333.     float fVec2[8192];

  334.     float fConst13;

  335.     int iConst14;

  336.     float fVec3[1024];

  337.     int iConst15;

  338.     float fRec12[2];

  339.     float fConst16;

  340.     float fConst17;

  341.     float fRec19[2];

  342.     float fRec18[2];

  343.     float fVec4[8192];

  344.     float fConst18;

  345.     int iConst19;

  346.     float fVec5[2048];

  347.     int iConst20;

  348.     float fRec16[2];

  349.     float fConst21;

  350.     float fConst22;

  351.     float fRec23[2];

  352.     float fRec22[2];

  353.     float fVec6[8192];

  354.     float fConst23;

  355.     int iConst24;

  356.     float fVec7[1024];

  357.     int iConst25;

  358.     float fRec20[2];

  359.     float fConst26;

  360.     float fConst27;

  361.     float fRec27[2];

  362.     float fRec26[2];

  363.     float fVec8[16384];

  364.     float fConst28;

  365.     int iConst29;

  366.     float fVec9[2048];

  367.     int iConst30;

  368.     float fRec24[2];

  369.     float fConst31;

  370.     float fConst32;

  371.     float fRec31[2];

  372.     float fRec30[2];

  373.     float fVec10[8192];

  374.     float fConst33;

  375.     int iConst34;

  376.     float fVec11[2048];

  377.     int iConst35;

  378.     float fRec28[2];

  379.     float fConst36;

  380.     float fConst37;

  381.     float fRec35[2];

  382.     float fRec34[2];

  383.     float fVec12[16384];

  384.     float fConst38;

  385.     int iConst39;

  386.     float fVec13[2048];

  387.     int iConst40;

  388.     float fRec32[2];

  389.     float fConst41;

  390.     float fConst42;

  391.     float fRec39[2];

  392.     float fRec38[2];

  393.     float fVec14[16384];

  394.     float fConst43;

  395.     int iConst44;

  396.     float fVec15[1024];

  397.     int iConst45;

  398.     float fRec36[2];

  399.     float fRec0[3];

  400.     float fRec1[3];

  401.     float fRec2[3];

  402.     float fRec3[3];

  403.     float fRec4[3];

  404.     float fRec5[3];

  405.     float fRec6[3];

  406.     float fRec7[3];

  407. 

  408.     /// Long nasty function setting initial values

  409.     void init(int samplingFreq)

  410.     {

  411.       fslider0 = 3.0f;

  412.       iConst0 = min(192000, max(1, samplingFreq));

  413.       fConst1 = floorf((0.5f + (0.174713f * iConst0)));

  414.       fConst2 = ((0 - (6.907755278982138f * fConst1)) / iConst0);

  415.       fslider1 = 6e+03f;

  416.       fConst3 = (6.283185307179586f / float(iConst0));

  417.       fConst4 = (1.0f / tanf((1256.6370614359173f / iConst0)));

  418.       fConst5 = (1 + fConst4);

  419.       fConst6 = (0 - ((1 - fConst4) / fConst5));

  420.       fConst7 = (1.0f / fConst5);

  421.       for (int i=0; i<2; i++) fRec11[i] = 0;

  422.       for (int i=0; i<2; i++) fRec10[i] = 0;

  423.       IOTA = 0;

  424.       for (int i=0; i<8192; i++) fVec0[i] = 0;

  425.       fConst8 = floorf((0.5f + (0.022904f * iConst0)));

  426.       iConst9 = int((int((fConst1 - fConst8)) & 8191));

  427.       for (int i=0; i<2048; i++) fVec1[i] = 0;

  428.       iConst10 = int((int((fConst8 - 1)) & 2047));

  429.       for (int i=0; i<2; i++) fRec8[i] = 0;

  430.       fConst11 = floorf((0.5f + (0.153129f * iConst0)));

  431.       fConst12 = ((0 - (6.907755278982138f * fConst11)) / iConst0);

  432.       for (int i=0; i<2; i++) fRec15[i] = 0;

  433.       for (int i=0; i<2; i++) fRec14[i] = 0;

  434.       for (int i=0; i<8192; i++) fVec2[i] = 0;

  435.       fConst13 = floorf((0.5f + (0.020346f * iConst0)));

  436.       iConst14 = int((int((fConst11 - fConst13)) & 8191));

  437.       for (int i=0; i<1024; i++) fVec3[i] = 0;

  438.       iConst15 = int((int((fConst13 - 1)) & 1023));

  439.       for (int i=0; i<2; i++) fRec12[i] = 0;

  440.       fConst16 = floorf((0.5f + (0.127837f * iConst0)));

  441.       fConst17 = ((0 - (6.907755278982138f * fConst16)) / iConst0);

  442.       for (int i=0; i<2; i++) fRec19[i] = 0;

  443.       for (int i=0; i<2; i++) fRec18[i] = 0;

  444.       for (int i=0; i<8192; i++) fVec4[i] = 0;

  445.       fConst18 = floorf((0.5f + (0.031604f * iConst0)));

  446.       iConst19 = int((int((fConst16 - fConst18)) & 8191));

  447.       for (int i=0; i<2048; i++) fVec5[i] = 0;

  448.       iConst20 = int((int((fConst18 - 1)) & 2047));

  449.       for (int i=0; i<2; i++) fRec16[i] = 0;

  450.       fConst21 = floorf((0.5f + (0.125f * iConst0)));

  451.       fConst22 = ((0 - (6.907755278982138f * fConst21)) / iConst0);

  452.       for (int i=0; i<2; i++) fRec23[i] = 0;

  453.       for (int i=0; i<2; i++) fRec22[i] = 0;

  454.       for (int i=0; i<8192; i++) fVec6[i] = 0;

  455.       fConst23 = floorf((0.5f + (0.013458f * iConst0)));

  456.       iConst24 = int((int((fConst21 - fConst23)) & 8191));

  457.       for (int i=0; i<1024; i++) fVec7[i] = 0;

  458.       iConst25 = int((int((fConst23 - 1)) & 1023));

  459.       for (int i=0; i<2; i++) fRec20[i] = 0;

  460.       fConst26 = floorf((0.5f + (0.210389f * iConst0)));

  461.       fConst27 = ((0 - (6.907755278982138f * fConst26)) / iConst0);

  462.       for (int i=0; i<2; i++) fRec27[i] = 0;

  463.       for (int i=0; i<2; i++) fRec26[i] = 0;

  464.       for (int i=0; i<16384; i++) fVec8[i] = 0;

  465.       fConst28 = floorf((0.5f + (0.024421f * iConst0)));

  466.       iConst29 = int((int((fConst26 - fConst28)) & 16383));

  467.       for (int i=0; i<2048; i++) fVec9[i] = 0;

  468.       iConst30 = int((int((fConst28 - 1)) & 2047));

  469.       for (int i=0; i<2; i++) fRec24[i] = 0;

  470.       fConst31 = floorf((0.5f + (0.192303f * iConst0)));

  471.       fConst32 = ((0 - (6.907755278982138f * fConst31)) / iConst0);

  472.       for (int i=0; i<2; i++) fRec31[i] = 0;

  473.       for (int i=0; i<2; i++) fRec30[i] = 0;

  474.       for (int i=0; i<8192; i++) fVec10[i] = 0;

  475.       fConst33 = floorf((0.5f + (0.029291f * iConst0)));

  476.       iConst34 = int((int((fConst31 - fConst33)) & 8191));

  477.       for (int i=0; i<2048; i++) fVec11[i] = 0;

  478.       iConst35 = int((int((fConst33 - 1)) & 2047));

  479.       for (int i=0; i<2; i++) fRec28[i] = 0;

  480.       fConst36 = floorf((0.5f + (0.256891f * iConst0)));

  481.       fConst37 = ((0 - (6.907755278982138f * fConst36)) / iConst0);

  482.       for (int i=0; i<2; i++) fRec35[i] = 0;

  483.       for (int i=0; i<2; i++) fRec34[i] = 0;

  484.       for (int i=0; i<16384; i++) fVec12[i] = 0;

  485.       fConst38 = floorf((0.5f + (0.027333f * iConst0)));

  486.       iConst39 = int((int((fConst36 - fConst38)) & 16383));

  487.       for (int i=0; i<2048; i++) fVec13[i] = 0;

  488.       iConst40 = int((int((fConst38 - 1)) & 2047));

  489.       for (int i=0; i<2; i++) fRec32[i] = 0;

  490.       fConst41 = floorf((0.5f + (0.219991f * iConst0)));

  491.       fConst42 = ((0 - (6.907755278982138f * fConst41)) / iConst0);

  492.       for (int i=0; i<2; i++) fRec39[i] = 0;

  493.       for (int i=0; i<2; i++) fRec38[i] = 0;

  494.       for (int i=0; i<16384; i++) fVec14[i] = 0;

  495.       fConst43 = floorf((0.5f + (0.019123f * iConst0)));

  496.       iConst44 = int((int((fConst41 - fConst43)) & 16383));

  497.       for (int i=0; i<1024; i++) fVec15[i] = 0;

  498.       iConst45 = int((int((fConst43 - 1)) & 1023));

  499.       for (int i=0; i<2; i++) fRec36[i] = 0;

  500.       for (int i=0; i<3; i++) fRec0[i] = 0;

  501.       for (int i=0; i<3; i++) fRec1[i] = 0;

  502.       for (int i=0; i<3; i++) fRec2[i] = 0;

  503.       for (int i=0; i<3; i++) fRec3[i] = 0;

  504.       for (int i=0; i<3; i++) fRec4[i] = 0;

  505.       for (int i=0; i<3; i++) fRec5[i] = 0;

  506.       for (int i=0; i<3; i++) fRec6[i] = 0;

  507.       for (int i=0; i<3; i++) fRec7[i] = 0;

  508.     }

  509. };

  510. 

  511. #endif // OPENAV_DSP_REVERB_H