KXStudio
/
Repository
mirror of https://github.com/KXStudio/Repository
1
0
Fork 0
Code Releases 2 Activity
External, Non-PPA KXStudio Repository
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.
101 Commits
1 Branch
2.9MB
Tree: 73434d89bc
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '73434d89bc'
${ noResults }
Repository/sources/plugins/airwindows-lv2/debian/patches/01_optimizations.patch

2925 lines
98KB
Raw Blame History 

  1. diff --git a/CMakeLists.txt b/CMakeLists.txt

  2. index 9636c5e..6696c38 100644

  3. --- a/CMakeLists.txt

  4. +++ b/CMakeLists.txt

  5. @@ -39,52 +39,64 @@ function(add_plugin PLUGIN_NAME)

  6.    configure_file(src/${PLUGIN_NAME}/${PLUGIN_NAME}.ttl ${CMAKE_LIBRARY_OUTPUT_DIRECTORY}/${PLUGIN_NAME}.ttl)

  7.  endfunction()

  8. 

  9. -add_plugin(Acceleration)

  10. +# Yay Distortions

  11. +add_plugin(Dyno)

  12. +add_plugin(Mojo)

  13. +add_plugin(Spiral)

  14. +

  15. +# Yay Reverbs

  16. +add_plugin(Galactic)

  17. +add_plugin(PocketVerbs)

  18. +add_plugin(StarChild)

  19. +

  20. +# Yay Misc

  21. +add_plugin(DeBess)

  22. +add_plugin(Baxandall)

  23. +add_plugin(Channel8)

  24. +add_plugin(EveryTrim)

  25. +add_plugin(Mackity)

  26. +

  27. +# Maybe

  28.  add_plugin(Acceleration2)

  29.  add_plugin(ADClip7)

  30. -add_plugin(Baxandall)

  31. -add_plugin(BitShiftGain)

  32.  add_plugin(Capacitor)

  33.  add_plugin(Capacitor2)

  34. -add_plugin(Channel8)

  35.  add_plugin(ClipOnly)

  36.  add_plugin(ClipOnly2)

  37.  add_plugin(Console7Buss)

  38.  add_plugin(Console7Cascade)

  39.  add_plugin(Console7Channel)

  40.  add_plugin(Console7Crunch)

  41. -add_plugin(DCVoltage)

  42. -add_plugin(DeBess)

  43. -add_plugin(Dyno)

  44.  add_plugin(EdIsDim)

  45. -add_plugin(EveryTrim)

  46. -add_plugin(Galactic)

  47. +add_plugin(MidSide)

  48. +add_plugin(MV)

  49. +add_plugin(Nikola)

  50. +add_plugin(Pressure5)

  51. +add_plugin(Sidepass)

  52. +

  53. +# Nay

  54. +add_plugin(Acceleration)

  55. +add_plugin(BitShiftGain)

  56. +add_plugin(DCVoltage)

  57.  add_plugin(HermeTrim)

  58.  add_plugin(Hypersonic)

  59.  add_plugin(Infrasonic)

  60.  add_plugin(Interstage)

  61.  add_plugin(LeftoMono)

  62. -add_plugin(Mackity)

  63. -add_plugin(MidSide)

  64. -add_plugin(Mojo)

  65. -add_plugin(MV)

  66. -add_plugin(Nikola)

  67. -add_plugin(PocketVerbs)

  68. -add_plugin(Pressure5)

  69.  add_plugin(PurestConsoleBuss)

  70.  add_plugin(PurestConsoleChannel)

  71.  add_plugin(PurestGain)

  72.  add_plugin(RightoMono)

  73. -add_plugin(Sidepass)

  74.  add_plugin(Slew)

  75.  add_plugin(SlewOnly)

  76. -add_plugin(Spiral)

  77. -add_plugin(StarChild)

  78.  add_plugin(SubsOnly)

  79.  add_plugin(Ultrasonic)

  80.  add_plugin(UltrasonicLite)

  81.  add_plugin(UltrasonicMed)

  82. 

  83. +set (CMAKE_C_FLAGS "${CMAKE_C_FLAGS} -fsingle-precision-constant -std=gnu99 -Wno-implicit-fallthrough -Werror=incompatible-pointer-types")

  84. +set (CMAKE_SHARED_LINKER_FLAGS "${CMAKE_SHARED_LINKER_FLAGS} -fsingle-precision-constant")

  85. +

  86.  configure_file(src/manifest.ttl.in ${CMAKE_LIBRARY_OUTPUT_DIRECTORY}/manifest.ttl)

  87. 

  88.  if(APPLE OR UNIX)

  89. diff --git a/src/ADClip7/ADClip7.c b/src/ADClip7/ADClip7.c

  90. index e502a66..9f0c763 100644

  91. --- a/src/ADClip7/ADClip7.c

  92. +++ b/src/ADClip7/ADClip7.c

  93. @@ -3,6 +3,7 @@

  94.  #include <math.h>

  95.  #include <stdint.h>

  96.  #include <stdlib.h>

  97. +#include "../optimizations.h"

  98. 

  99.  #define ADCLIP7_URI "https://hannesbraun.net/ns/lv2/airwindows/adclip7"

  100. 

  101. @@ -26,8 +27,6 @@ typedef struct {

  102.  	const float* enhance;

  103.  	const float* mode;

  104. 

  105. -	uint32_t fpdL;

  106. -	uint32_t fpdR;

  107.  	//default stuff

  108.  	double lastSampleL;

  109.  	double lastSampleR;

  110. @@ -46,7 +45,7 @@ typedef struct {

  111. 

  112.  static LV2_Handle instantiate(

  113.  	const LV2_Descriptor* descriptor,

  114. -	double rate,

  115. +	realdouble rate,

  116.  	const char* bundle_path,

  117.  	const LV2_Feature* const* features)

  118.  {

  119. @@ -106,11 +105,6 @@ static void activate(LV2_Handle instance)

  120.  	adclip7->iirLowsAR = 0.0;

  121.  	adclip7->iirLowsBL = 0.0;

  122.  	adclip7->iirLowsBR = 0.0;

  123. -

  124. -	adclip7->fpdL = 1.0;

  125. -	while (adclip7->fpdL < 16386) adclip7->fpdL = rand() * UINT32_MAX;

  126. -	adclip7->fpdR = 1.0;

  127. -	while (adclip7->fpdR < 16386) adclip7->fpdR = rand() * UINT32_MAX;

  128.  }

  129. 

  130.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  131. @@ -169,9 +163,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  132.  	while (sampleFrames-- > 0) {

  133.  		inputSampleL = *in1;

  134.  		inputSampleR = *in2;

  135. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = adclip7->fpdL * 1.18e-17;

  136. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = adclip7->fpdR * 1.18e-17;

  137. -

  138. 

  139.  		if (inputGain != 1.0) {

  140.  			inputSampleL *= inputGain;

  141. @@ -517,20 +508,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  142.  		if (inputSampleR < -adclip7->refclipR) inputSampleR = -adclip7->refclipR;

  143.  		//final iron bar

  144. 

  145. -		//begin 32 bit stereo floating point dither

  146. -		int expon;

  147. -		frexpf((float) inputSampleL, &expon);

  148. -		adclip7->fpdL ^= adclip7->fpdL << 13;

  149. -		adclip7->fpdL ^= adclip7->fpdL >> 17;

  150. -		adclip7->fpdL ^= adclip7->fpdL << 5;

  151. -		inputSampleL += (((double) adclip7->fpdL - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  152. -		frexpf((float) inputSampleR, &expon);

  153. -		adclip7->fpdR ^= adclip7->fpdR << 13;

  154. -		adclip7->fpdR ^= adclip7->fpdR >> 17;

  155. -		adclip7->fpdR ^= adclip7->fpdR << 5;

  156. -		inputSampleR += (((double) adclip7->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  157. -		//end 32 bit stereo floating point dither

  158. -

  159.  		*out1 = (float) inputSampleL;

  160.  		*out2 = (float) inputSampleR;

  161. 

  162. diff --git a/src/Acceleration/Acceleration.c b/src/Acceleration/Acceleration.c

  163. index d6af6c5..f4366a0 100644

  164. --- a/src/Acceleration/Acceleration.c

  165. +++ b/src/Acceleration/Acceleration.c

  166. @@ -3,6 +3,7 @@

  167.  #include <math.h>

  168.  #include <stdint.h>

  169.  #include <stdlib.h>

  170. +#include "../optimizations.h"

  171. 

  172.  #define ACCELERATION_URI "https://hannesbraun.net/ns/lv2/airwindows/acceleration"

  173. 

  174. @@ -22,8 +23,6 @@ typedef struct {

  175.  	const float* limit;

  176.  	const float* drywet;

  177. 

  178. -	uint32_t fpdL;

  179. -	uint32_t fpdR;

  180.  	//default stuff

  181.  	double ataLastOutL;

  182.  	double s1L;

  183. @@ -50,7 +49,7 @@ typedef struct {

  184. 

  185.  static LV2_Handle instantiate(

  186.  	const LV2_Descriptor* descriptor,

  187. -	double rate,

  188. +	realdouble rate,

  189.  	const char* bundle_path,

  190.  	const LV2_Feature* const* features)

  191.  {

  192. @@ -108,11 +107,6 @@ static void activate(LV2_Handle instance)

  193.  	acceleration->m1R = 0.0;

  194.  	acceleration->m2R = 0.0;

  195.  	acceleration->desR = 0.0;

  196. -

  197. -	acceleration->fpdL = 1.0;

  198. -	while (acceleration->fpdL < 16386) acceleration->fpdL = rand() * UINT32_MAX;

  199. -	acceleration->fpdR = 1.0;

  200. -	while (acceleration->fpdR < 16386) acceleration->fpdR = rand() * UINT32_MAX;

  201.  }

  202. 

  203.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  204. @@ -145,8 +139,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  205.  	while (sampleFrames-- > 0) {

  206.  		inputSampleL = *in1;

  207.  		inputSampleR = *in2;

  208. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = acceleration->fpdL * 1.18e-17;

  209. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = acceleration->fpdR * 1.18e-17;

  210.  		drySampleL = inputSampleL;

  211.  		drySampleR = inputSampleR;

  212. 

  213. @@ -205,20 +197,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  214.  			inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0 - wet));

  215.  		}

  216. 

  217. -		//begin 32 bit stereo floating point dither

  218. -		int expon;

  219. -		frexpf((float) inputSampleL, &expon);

  220. -		acceleration->fpdL ^= acceleration->fpdL << 13;

  221. -		acceleration->fpdL ^= acceleration->fpdL >> 17;

  222. -		acceleration->fpdL ^= acceleration->fpdL << 5;

  223. -		inputSampleL += (((double) acceleration->fpdL - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  224. -		frexpf((float) inputSampleR, &expon);

  225. -		acceleration->fpdR ^= acceleration->fpdR << 13;

  226. -		acceleration->fpdR ^= acceleration->fpdR >> 17;

  227. -		acceleration->fpdR ^= acceleration->fpdR << 5;

  228. -		inputSampleR += (((double) acceleration->fpdR - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  229. -		//end 32 bit stereo floating point dither

  230. -

  231.  		*out1 = (float) inputSampleL;

  232.  		*out2 = (float) inputSampleR;

  233. 

  234. diff --git a/src/Acceleration2/Acceleration2.c b/src/Acceleration2/Acceleration2.c

  235. index 3cdf047..88b5f03 100644

  236. --- a/src/Acceleration2/Acceleration2.c

  237. +++ b/src/Acceleration2/Acceleration2.c

  238. @@ -3,6 +3,7 @@

  239.  #include <math.h>

  240.  #include <stdint.h>

  241.  #include <stdlib.h>

  242. +#include "../optimizations.h"

  243. 

  244.  #define ACCELERATION2_URI "https://hannesbraun.net/ns/lv2/airwindows/acceleration2"

  245. 

  246. @@ -36,14 +37,11 @@ typedef struct {

  247.  	double biquadB[11]; //note that this stereo form doesn't require L and R forms!

  248.  	//This is because so much of it is coefficients etc. that are the same on both channels.

  249.  	//So the stored samples are in 7-8 and 9-10, and freq/res/coefficients serve both.

  250. -

  251. -	uint32_t fpdL;

  252. -	uint32_t fpdR;

  253.  } Acceleration2;

  254. 

  255.  static LV2_Handle instantiate(

  256.  	const LV2_Descriptor* descriptor,

  257. -	double rate,

  258. +	realdouble rate,

  259.  	const char* bundle_path,

  260.  	const LV2_Feature* const* features)

  261.  {

  262. @@ -95,11 +93,6 @@ static void activate(LV2_Handle instance)

  263.  		acceleration2->biquadA[x] = 0.0;

  264.  		acceleration2->biquadB[x] = 0.0;

  265.  	}

  266. -

  267. -	acceleration2->fpdL = 1.0;

  268. -	while (acceleration2->fpdL < 16386) acceleration2->fpdL = rand() * UINT32_MAX;

  269. -	acceleration2->fpdR = 1.0;

  270. -	while (acceleration2->fpdR < 16386) acceleration2->fpdR = rand() * UINT32_MAX;

  271.  }

  272. 

  273.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  274. @@ -145,8 +138,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  275.  	while (sampleFrames-- > 0) {

  276.  		double inputSampleL = *in1;

  277.  		double inputSampleR = *in2;

  278. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = acceleration2->fpdL * 1.18e-17;

  279. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = acceleration2->fpdR * 1.18e-17;

  280.  		double drySampleL = inputSampleL;

  281.  		double drySampleR = inputSampleR;

  282. 

  283. @@ -194,20 +185,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  284.  			inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0 - wet));

  285.  		}

  286. 

  287. -		//begin 32 bit stereo floating point dither

  288. -		int expon;

  289. -		frexpf((float)inputSampleL, &expon);

  290. -		acceleration2->fpdL ^= acceleration2->fpdL << 13;

  291. -		acceleration2->fpdL ^= acceleration2->fpdL >> 17;

  292. -		acceleration2->fpdL ^= acceleration2->fpdL << 5;

  293. -		inputSampleL += (((double)acceleration2->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  294. -		frexpf((float)inputSampleR, &expon);

  295. -		acceleration2->fpdR ^= acceleration2->fpdR << 13;

  296. -		acceleration2->fpdR ^= acceleration2->fpdR >> 17;

  297. -		acceleration2->fpdR ^= acceleration2->fpdR << 5;

  298. -		inputSampleR += (((double)acceleration2->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  299. -		//end 32 bit stereo floating point dither

  300. -

  301.  		*out1 = (float) inputSampleL;

  302.  		*out2 = (float) inputSampleR;

  303. 

  304. diff --git a/src/Baxandall/Baxandall.c b/src/Baxandall/Baxandall.c

  305. index 263b1ca..82d7f52 100644

  306. --- a/src/Baxandall/Baxandall.c

  307. +++ b/src/Baxandall/Baxandall.c

  308. @@ -4,6 +4,7 @@

  309.  #include <stdbool.h>

  310.  #include <stdint.h>

  311.  #include <stdlib.h>

  312. +#include "../optimizations.h"

  313. 

  314.  #define BAXANDALL_URI "https://hannesbraun.net/ns/lv2/airwindows/baxandall"

  315. 

  316. @@ -25,8 +26,6 @@ typedef struct {

  317.  	const float* bass;

  318.  	const float* outputGain;

  319. 

  320. -	uint32_t fpdL;

  321. -	uint32_t fpdR;

  322.  	//default stuff

  323.  	double trebleAL[9];

  324.  	double trebleBL[9];

  325. @@ -42,7 +41,7 @@ typedef struct {

  326. 

  327.  static LV2_Handle instantiate(

  328.  	const LV2_Descriptor* descriptor,

  329. -	double rate,

  330. +	realdouble rate,

  331.  	const char* bundle_path,

  332.  	const LV2_Feature* const* features)

  333.  {

  334. @@ -94,10 +93,6 @@ static void activate(LV2_Handle instance)

  335.  		baxandall->bassBR[x] = 0.0;

  336.  	}

  337.  	baxandall->flip = false;

  338. -	baxandall->fpdL = 1.0;

  339. -	while (baxandall->fpdL < 16386) baxandall->fpdL = rand() * UINT32_MAX;

  340. -	baxandall->fpdR = 1.0;

  341. -	while (baxandall->fpdR < 16386) baxandall->fpdR = rand() * UINT32_MAX;

  342.  }

  343. 

  344.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  345. @@ -183,8 +178,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  346.  	while (sampleFrames-- > 0) {

  347.  		double inputSampleL = *in1;

  348.  		double inputSampleR = *in2;

  349. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = baxandall->fpdL * 1.18e-17;

  350. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = baxandall->fpdR * 1.18e-17;

  351. 

  352.  		if (output != 1.0) {

  353.  			inputSampleL *= output;

  354. @@ -258,20 +251,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  355.  		inputSampleR = asin(inputSampleR);

  356.  		//amplitude aspect

  357. 

  358. -		//begin 32 bit stereo floating point dither

  359. -		int expon;

  360. -		frexpf((float) inputSampleL, &expon);

  361. -		baxandall->fpdL ^= baxandall->fpdL << 13;

  362. -		baxandall->fpdL ^= baxandall->fpdL >> 17;

  363. -		baxandall->fpdL ^= baxandall->fpdL << 5;

  364. -		inputSampleL += (((double) baxandall->fpdL - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  365. -		frexpf((float) inputSampleR, &expon);

  366. -		baxandall->fpdR ^= baxandall->fpdR << 13;

  367. -		baxandall->fpdR ^= baxandall->fpdR >> 17;

  368. -		baxandall->fpdR ^= baxandall->fpdR << 5;

  369. -		inputSampleR += (((double) baxandall->fpdR - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  370. -		//end 32 bit stereo floating point dither

  371. -

  372.  		*out1 = (float) inputSampleL;

  373.  		*out2 = (float) inputSampleR;

  374. 

  375. diff --git a/src/BitShiftGain/BitShiftGain.c b/src/BitShiftGain/BitShiftGain.c

  376. index 4c5a5d1..1235099 100644

  377. --- a/src/BitShiftGain/BitShiftGain.c

  378. +++ b/src/BitShiftGain/BitShiftGain.c

  379. @@ -2,6 +2,7 @@

  380. 

  381.  #include <stdint.h>

  382.  #include <stdlib.h>

  383. +#include "../optimizations.h"

  384. 

  385.  #define BITSHIFTGAIN_URI "https://hannesbraun.net/ns/lv2/airwindows/bitshiftgain"

  386. 

  387. @@ -21,7 +22,7 @@ typedef struct {

  388. 

  389.  static LV2_Handle instantiate(

  390.  	const LV2_Descriptor* descriptor,

  391. -	double rate,

  392. +	realdouble rate,

  393.  	const char* bundle_path,

  394.  	const LV2_Feature* const* features)

  395.  {

  396. diff --git a/src/Capacitor/Capacitor.c b/src/Capacitor/Capacitor.c

  397. index 907c2d3..fe0d466 100644

  398. --- a/src/Capacitor/Capacitor.c

  399. +++ b/src/Capacitor/Capacitor.c

  400. @@ -3,6 +3,7 @@

  401.  #include <math.h>

  402.  #include <stdint.h>

  403.  #include <stdlib.h>

  404. +#include "../optimizations.h"

  405. 

  406.  #define CAPACITOR_URI "https://hannesbraun.net/ns/lv2/airwindows/capacitor"

  407. 

  408. @@ -62,14 +63,11 @@ typedef struct {

  409.  	double lastLowpass;

  410.  	double lastHighpass;

  411.  	double lastWet;

  412. -

  413. -	uint32_t fpdL;

  414. -	uint32_t fpdR;

  415.  } Capacitor;

  416. 

  417.  static LV2_Handle instantiate(

  418.  	const LV2_Descriptor* descriptor,

  419. -	double rate,

  420. +	realdouble rate,

  421.  	const char* bundle_path,

  422.  	const LV2_Feature* const* features)

  423.  {

  424. @@ -145,11 +143,6 @@ static void activate(LV2_Handle instance)

  425.  	capacitor->lastLowpass = 1000.0;

  426.  	capacitor->lastHighpass = 1000.0;

  427.  	capacitor->lastWet = 1000.0;

  428. -

  429. -	capacitor->fpdL = 1.0;

  430. -	while (capacitor->fpdL < 16386) capacitor->fpdL = rand() * UINT32_MAX;

  431. -	capacitor->fpdR = 1.0;

  432. -	while (capacitor->fpdR < 16386) capacitor->fpdR = rand() * UINT32_MAX;

  433.  }

  434. 

  435.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  436. @@ -186,8 +179,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  437.  	while (sampleFrames-- > 0) {

  438.  		inputSampleL = *in1;

  439.  		inputSampleR = *in2;

  440. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = capacitor->fpdL * 1.18e-17;

  441. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = capacitor->fpdR * 1.18e-17;

  442.  		drySampleL = inputSampleL;

  443.  		drySampleR = inputSampleR;

  444. 

  445. @@ -365,20 +356,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  446.  		inputSampleL = (drySampleL * dry) + (inputSampleL * capacitor->wet);

  447.  		inputSampleR = (drySampleR * dry) + (inputSampleR * capacitor->wet);

  448. 

  449. -		//begin 32 bit stereo floating point dither

  450. -		int expon;

  451. -		frexpf((float)inputSampleL, &expon);

  452. -		capacitor->fpdL ^= capacitor->fpdL << 13;

  453. -		capacitor->fpdL ^= capacitor->fpdL >> 17;

  454. -		capacitor->fpdL ^= capacitor->fpdL << 5;

  455. -		inputSampleL += (((double)capacitor->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  456. -		frexpf((float)inputSampleR, &expon);

  457. -		capacitor->fpdR ^= capacitor->fpdR << 13;

  458. -		capacitor->fpdR ^= capacitor->fpdR >> 17;

  459. -		capacitor->fpdR ^= capacitor->fpdR << 5;

  460. -		inputSampleR += (((double)capacitor->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  461. -		//end 32 bit stereo floating point dither

  462. -

  463.  		*out1 = (float) inputSampleL;

  464.  		*out2 = (float) inputSampleR;

  465. 

  466. diff --git a/src/Capacitor2/Capacitor2.c b/src/Capacitor2/Capacitor2.c

  467. index 111ff8a..53050eb 100644

  468. --- a/src/Capacitor2/Capacitor2.c

  469. +++ b/src/Capacitor2/Capacitor2.c

  470. @@ -3,6 +3,7 @@

  471.  #include <math.h>

  472.  #include <stdint.h>

  473.  #include <stdlib.h>

  474. +#include "../optimizations.h"

  475. 

  476.  #define CAPACITOR2_URI "https://hannesbraun.net/ns/lv2/airwindows/capacitor2"

  477. 

  478. @@ -65,14 +66,11 @@ typedef struct {

  479.  	double lastLowpass;

  480.  	double lastHighpass;

  481.  	double lastWet;

  482. -

  483. -	uint32_t fpdL;

  484. -	uint32_t fpdR;

  485.  } Capacitor2;

  486. 

  487.  static LV2_Handle instantiate(

  488.  	const LV2_Descriptor* descriptor,

  489. -	double rate,

  490. +	realdouble rate,

  491.  	const char* bundle_path,

  492.  	const LV2_Feature* const* features)

  493.  {

  494. @@ -152,11 +150,6 @@ static void activate(LV2_Handle instance)

  495.  	capacitor2->lastLowpass = 1000.0;

  496.  	capacitor2->lastHighpass = 1000.0;

  497.  	capacitor2->lastWet = 1000.0;

  498. -

  499. -	capacitor2->fpdL = 1.0;

  500. -	while (capacitor2->fpdL < 16386) capacitor2->fpdL = rand() * UINT32_MAX;

  501. -	capacitor2->fpdR = 1.0;

  502. -	while (capacitor2->fpdR < 16386) capacitor2->fpdR = rand() * UINT32_MAX;

  503.  }

  504. 

  505.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  506. @@ -185,8 +178,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  507.  	while (sampleFrames-- > 0) {

  508.  		double inputSampleL = *in1;

  509.  		double inputSampleR = *in2;

  510. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = capacitor2->fpdL * 1.18e-17;

  511. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = capacitor2->fpdR * 1.18e-17;

  512.  		double drySampleL = inputSampleL;

  513.  		double drySampleR = inputSampleR;

  514. 

  515. @@ -377,20 +368,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  516.  		inputSampleL = (drySampleL * (1.0 - capacitor2->wet)) + (inputSampleL * nonLinTrim * capacitor2->wet);

  517.  		inputSampleR = (drySampleR * (1.0 - capacitor2->wet)) + (inputSampleR * nonLinTrim * capacitor2->wet);

  518. 

  519. -		//begin 32 bit stereo floating point dither

  520. -		int expon;

  521. -		frexpf((float)inputSampleL, &expon);

  522. -		capacitor2->fpdL ^= capacitor2->fpdL << 13;

  523. -		capacitor2->fpdL ^= capacitor2->fpdL >> 17;

  524. -		capacitor2->fpdL ^= capacitor2->fpdL << 5;

  525. -		inputSampleL += (((double)capacitor2->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  526. -		frexpf((float)inputSampleR, &expon);

  527. -		capacitor2->fpdR ^= capacitor2->fpdR << 13;

  528. -		capacitor2->fpdR ^= capacitor2->fpdR >> 17;

  529. -		capacitor2->fpdR ^= capacitor2->fpdR << 5;

  530. -		inputSampleR += (((double)capacitor2->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  531. -		//end 32 bit stereo floating point dither

  532. -

  533.  		*out1 = (float) inputSampleL;

  534.  		*out2 = (float) inputSampleR;

  535. 

  536. diff --git a/src/Channel8/Channel8.c b/src/Channel8/Channel8.c

  537. index 5130655..80e75a7 100644

  538. --- a/src/Channel8/Channel8.c

  539. +++ b/src/Channel8/Channel8.c

  540. @@ -4,6 +4,7 @@

  541.  #include <stdbool.h>

  542.  #include <stdint.h>

  543.  #include <stdlib.h>

  544. +#include "../optimizations.h"

  545. 

  546.  #define CHANNEL8_URI "https://hannesbraun.net/ns/lv2/airwindows/channel8"

  547. 

  548. @@ -25,8 +26,6 @@ typedef struct {

  549.  	const float* drive;

  550.  	const float* outputGain;

  551. 

  552. -	uint32_t fpdL;

  553. -	uint32_t fpdR;

  554.  	//default stuff

  555.  	double iirSampleLA;

  556.  	double iirSampleRA;

  557. @@ -45,7 +44,7 @@ typedef struct {

  558. 

  559.  static LV2_Handle instantiate(

  560.  	const LV2_Descriptor* descriptor,

  561. -	double rate,

  562. +	realdouble rate,

  563.  	const char* bundle_path,

  564.  	const LV2_Feature* const* features)

  565.  {

  566. @@ -87,10 +86,6 @@ static void activate(LV2_Handle instance)

  567.  {

  568.  	Channel8* channel8 = (Channel8*) instance;

  569. 

  570. -	channel8->fpdL = 1.0;

  571. -	while (channel8->fpdL < 16386) channel8->fpdL = rand() * UINT32_MAX;

  572. -	channel8->fpdR = 1.0;

  573. -	while (channel8->fpdR < 16386) channel8->fpdR = rand() * UINT32_MAX;

  574.  	channel8->iirSampleLA = 0.0;

  575.  	channel8->iirSampleRA = 0.0;

  576.  	channel8->iirSampleLB = 0.0;

  577. @@ -150,8 +145,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  578.  	while (sampleFrames-- > 0) {

  579.  		double inputSampleL = *in1;

  580.  		double inputSampleR = *in2;

  581. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = channel8->fpdL * 1.18e-17;

  582. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = channel8->fpdR * 1.18e-17;

  583. 

  584.  		double dielectricScaleL = fabs(2.0 - ((inputSampleL + nonLin) / nonLin));

  585.  		double dielectricScaleR = fabs(2.0 - ((inputSampleR + nonLin) / nonLin));

  586. @@ -236,20 +229,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  587.  			inputSampleR *= output;

  588.  		}

  589. 

  590. -		//begin 32 bit stereo floating point dither

  591. -		int expon;

  592. -		frexpf((float)inputSampleL, &expon);

  593. -		channel8->fpdL ^= channel8->fpdL << 13;

  594. -		channel8->fpdL ^= channel8->fpdL >> 17;

  595. -		channel8->fpdL ^= channel8->fpdL << 5;

  596. -		inputSampleL += (((double)channel8->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  597. -		frexpf((float)inputSampleR, &expon);

  598. -		channel8->fpdR ^= channel8->fpdR << 13;

  599. -		channel8->fpdR ^= channel8->fpdR >> 17;

  600. -		channel8->fpdR ^= channel8->fpdR << 5;

  601. -		inputSampleR += (((double)channel8->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  602. -		//end 32 bit stereo floating point dither

  603. -

  604.  		*out1 = (float) inputSampleL;

  605.  		*out2 = (float) inputSampleR;

  606. 

  607. diff --git a/src/ClipOnly/ClipOnly.c b/src/ClipOnly/ClipOnly.c

  608. index 2d83605..94f0ef8 100644

  609. --- a/src/ClipOnly/ClipOnly.c

  610. +++ b/src/ClipOnly/ClipOnly.c

  611. @@ -3,6 +3,7 @@

  612.  #include <stdbool.h>

  613.  #include <stdint.h>

  614.  #include <stdlib.h>

  615. +#include "../optimizations.h"

  616. 

  617.  #define CLIPONLY_URI "https://hannesbraun.net/ns/lv2/airwindows/cliponly"

  618. 

  619. @@ -28,7 +29,7 @@ typedef struct {

  620. 

  621.  static LV2_Handle instantiate(

  622.  	const LV2_Descriptor* descriptor,

  623. -	double rate,

  624. +	realdouble rate,

  625.  	const char* bundle_path,

  626.  	const LV2_Feature* const* features)

  627.  {

  628. diff --git a/src/ClipOnly2/ClipOnly2.c b/src/ClipOnly2/ClipOnly2.c

  629. index eec1237..7f28c8c 100644

  630. --- a/src/ClipOnly2/ClipOnly2.c

  631. +++ b/src/ClipOnly2/ClipOnly2.c

  632. @@ -4,6 +4,7 @@

  633.  #include <stdbool.h>

  634.  #include <stdint.h>

  635.  #include <stdlib.h>

  636. +#include "../optimizations.h"

  637. 

  638.  #define CLIPONLY2_URI "https://hannesbraun.net/ns/lv2/airwindows/cliponly2"

  639. 

  640. @@ -31,7 +32,7 @@ typedef struct {

  641. 

  642.  static LV2_Handle instantiate(

  643.  	const LV2_Descriptor* descriptor,

  644. -	double rate,

  645. +	realdouble rate,

  646.  	const char* bundle_path,

  647.  	const LV2_Feature* const* features)

  648.  {

  649. diff --git a/src/Console7Buss/Console7Buss.c b/src/Console7Buss/Console7Buss.c

  650. index f6f006e..2cd7e8b 100644

  651. --- a/src/Console7Buss/Console7Buss.c

  652. +++ b/src/Console7Buss/Console7Buss.c

  653. @@ -3,6 +3,7 @@

  654.  #include <math.h>

  655.  #include <stdint.h>

  656.  #include <stdlib.h>

  657. +#include "../optimizations.h"

  658. 

  659.  #define CONSOLE7BUSS_URI "https://hannesbraun.net/ns/lv2/airwindows/console7buss"

  660. 

  661. @@ -25,14 +26,11 @@ typedef struct {

  662. 

  663.  	double biquadA[15];

  664.  	double biquadB[15];

  665. -

  666. -	uint32_t fpdL;

  667. -	uint32_t fpdR;

  668.  } Console7Buss;

  669. 

  670.  static LV2_Handle instantiate(

  671.  	const LV2_Descriptor* descriptor,

  672. -	double rate,

  673. +	realdouble rate,

  674.  	const char* bundle_path,

  675.  	const LV2_Feature* const* features)

  676.  {

  677. @@ -74,10 +72,6 @@ static void activate(LV2_Handle instance)

  678.  		console7Buss->biquadA[x] = 0.0;

  679.  		console7Buss->biquadB[x] = 0.0;

  680.  	}

  681. -	console7Buss->fpdL = 1.0;

  682. -	while (console7Buss->fpdL < 16386) console7Buss->fpdL = rand() * UINT32_MAX;

  683. -	console7Buss->fpdR = 1.0;

  684. -	while (console7Buss->fpdR < 16386) console7Buss->fpdR = rand() * UINT32_MAX;

  685.  }

  686. 

  687.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  688. @@ -118,8 +112,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  689.  	while (sampleFrames-- > 0) {

  690.  		double inputSampleL = *in1;

  691.  		double inputSampleR = *in2;

  692. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = console7Buss->fpdL * 1.18e-17;

  693. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = console7Buss->fpdR * 1.18e-17;

  694. 

  695.  		double outSampleL = console7Buss->biquadA[2] * inputSampleL + console7Buss->biquadA[3] * console7Buss->biquadA[7] + console7Buss->biquadA[4] * console7Buss->biquadA[8] - console7Buss->biquadA[5] * console7Buss->biquadA[9] - console7Buss->biquadA[6] * console7Buss->biquadA[10];

  696.  		console7Buss->biquadA[8] = console7Buss->biquadA[7];

  697. @@ -176,20 +168,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  698.  		}

  699.  		//we re-amplify after the distortion relative to how much we cut back previously.

  700. 

  701. -		//begin 32 bit stereo floating point dither

  702. -		int expon;

  703. -		frexpf((float)inputSampleL, &expon);

  704. -		console7Buss->fpdL ^= console7Buss->fpdL << 13;

  705. -		console7Buss->fpdL ^= console7Buss->fpdL >> 17;

  706. -		console7Buss->fpdL ^= console7Buss->fpdL << 5;

  707. -		inputSampleL += (((double)console7Buss->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  708. -		frexpf((float)inputSampleR, &expon);

  709. -		console7Buss->fpdR ^= console7Buss->fpdR << 13;

  710. -		console7Buss->fpdR ^= console7Buss->fpdR >> 17;

  711. -		console7Buss->fpdR ^= console7Buss->fpdR << 5;

  712. -		inputSampleR += (((double)console7Buss->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  713. -		//end 32 bit stereo floating point dither

  714. -

  715.  		*out1 = (float) inputSampleL;

  716.  		*out2 = (float) inputSampleR;

  717. 

  718. diff --git a/src/Console7Cascade/Console7Cascade.c b/src/Console7Cascade/Console7Cascade.c

  719. index 2f5f8ee..e4fd69d 100644

  720. --- a/src/Console7Cascade/Console7Cascade.c

  721. +++ b/src/Console7Cascade/Console7Cascade.c

  722. @@ -3,6 +3,7 @@

  723.  #include <math.h>

  724.  #include <stdint.h>

  725.  #include <stdlib.h>

  726. +#include "../optimizations.h"

  727. 

  728.  #define CONSOLE7CASCADE_URI "https://hannesbraun.net/ns/lv2/airwindows/console7cascade"

  729. 

  730. @@ -28,14 +29,11 @@ typedef struct {

  731.  	double biquadC[15];

  732.  	double biquadD[15];

  733.  	double biquadE[15];

  734. -

  735. -	uint32_t fpdL;

  736. -	uint32_t fpdR;

  737.  } Console7Cascade;

  738. 

  739.  static LV2_Handle instantiate(

  740.  	const LV2_Descriptor* descriptor,

  741. -	double rate,

  742. +	realdouble rate,

  743.  	const char* bundle_path,

  744.  	const LV2_Feature* const* features)

  745.  {

  746. @@ -80,10 +78,6 @@ static void activate(LV2_Handle instance)

  747.  		console7Cascade->biquadD[x] = 0.0;

  748.  		console7Cascade->biquadE[x] = 0.0;

  749.  	}

  750. -	console7Cascade->fpdL = 1.0;

  751. -	while (console7Cascade->fpdL < 16386) console7Cascade->fpdL = rand() * UINT32_MAX;

  752. -	console7Cascade->fpdR = 1.0;

  753. -	while (console7Cascade->fpdR < 16386) console7Cascade->fpdR = rand() * UINT32_MAX;

  754.  }

  755. 

  756.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  757. @@ -150,8 +144,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  758.  	while (sampleFrames-- > 0) {

  759.  		double inputSampleL = *in1;

  760.  		double inputSampleR = *in2;

  761. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = console7Cascade->fpdL * 1.18e-17;

  762. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = console7Cascade->fpdR * 1.18e-17;

  763. 

  764.  		double outSampleL = console7Cascade->biquadA[2] * inputSampleL + console7Cascade->biquadA[3] * console7Cascade->biquadA[7] + console7Cascade->biquadA[4] * console7Cascade->biquadA[8] - console7Cascade->biquadA[5] * console7Cascade->biquadA[9] - console7Cascade->biquadA[6] * console7Cascade->biquadA[10];

  765.  		console7Cascade->biquadA[8] = console7Cascade->biquadA[7];

  766. @@ -292,20 +284,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  767.  		}

  768.  		//we re-amplify after the distortion relative to how much we cut back previously.

  769. 

  770. -		//begin 32 bit stereo floating point dither

  771. -		int expon;

  772. -		frexpf((float)inputSampleL, &expon);

  773. -		console7Cascade->fpdL ^= console7Cascade->fpdL << 13;

  774. -		console7Cascade->fpdL ^= console7Cascade->fpdL >> 17;

  775. -		console7Cascade->fpdL ^= console7Cascade->fpdL << 5;

  776. -		inputSampleL += (((double)console7Cascade->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  777. -		frexpf((float)inputSampleR, &expon);

  778. -		console7Cascade->fpdR ^= console7Cascade->fpdR << 13;

  779. -		console7Cascade->fpdR ^= console7Cascade->fpdR >> 17;

  780. -		console7Cascade->fpdR ^= console7Cascade->fpdR << 5;

  781. -		inputSampleR += (((double)console7Cascade->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  782. -		//end 32 bit stereo floating point dither

  783. -

  784.  		*out1 = (float) inputSampleL;

  785.  		*out2 = (float) inputSampleR;

  786. 

  787. diff --git a/src/Console7Channel/Console7Channel.c b/src/Console7Channel/Console7Channel.c

  788. index 9d531c4..e76936a 100644

  789. --- a/src/Console7Channel/Console7Channel.c

  790. +++ b/src/Console7Channel/Console7Channel.c

  791. @@ -3,6 +3,7 @@

  792.  #include <math.h>

  793.  #include <stdint.h>

  794.  #include <stdlib.h>

  795. +#include "../optimizations.h"

  796. 

  797.  #define CONSOLE7CHANNEL_URI "https://hannesbraun.net/ns/lv2/airwindows/console7channel"

  798. 

  799. @@ -24,14 +25,11 @@ typedef struct {

  800.  	double chasespeed;

  801. 

  802.  	double biquadA[15];

  803. -

  804. -	uint32_t fpdL;

  805. -	uint32_t fpdR;

  806.  } Console7Channel;

  807. 

  808.  static LV2_Handle instantiate(

  809.  	const LV2_Descriptor* descriptor,

  810. -	double rate,

  811. +	realdouble rate,

  812.  	const char* bundle_path,

  813.  	const LV2_Feature* const* features)

  814.  {

  815. @@ -71,10 +69,6 @@ static void activate(LV2_Handle instance)

  816.  	for (int x = 0; x < 15; x++) {

  817.  		console7Channel->biquadA[x] = 0.0;

  818.  	}

  819. -	console7Channel->fpdL = 1.0;

  820. -	while (console7Channel->fpdL < 16386) console7Channel->fpdL = rand() * UINT32_MAX;

  821. -	console7Channel->fpdR = 1.0;

  822. -	while (console7Channel->fpdR < 16386) console7Channel->fpdR = rand() * UINT32_MAX;

  823.  }

  824. 

  825.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  826. @@ -109,8 +103,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  827.  	while (sampleFrames-- > 0) {

  828.  		double inputSampleL = *in1;

  829.  		double inputSampleR = *in2;

  830. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = console7Channel->fpdL * 1.18e-17;

  831. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = console7Channel->fpdR * 1.18e-17;

  832. 

  833.  		double outSampleL = console7Channel->biquadA[2] * inputSampleL + console7Channel->biquadA[3] * console7Channel->biquadA[7] + console7Channel->biquadA[4] * console7Channel->biquadA[8] - console7Channel->biquadA[5] * console7Channel->biquadA[9] - console7Channel->biquadA[6] * console7Channel->biquadA[10];

  834.  		console7Channel->biquadA[8] = console7Channel->biquadA[7];

  835. @@ -155,20 +147,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  836.  		}

  837.  		//we re-amplify after the distortion relative to how much we cut back previously.

  838. 

  839. -		//begin 32 bit stereo floating point dither

  840. -		int expon;

  841. -		frexpf((float)inputSampleL, &expon);

  842. -		console7Channel->fpdL ^= console7Channel->fpdL << 13;

  843. -		console7Channel->fpdL ^= console7Channel->fpdL >> 17;

  844. -		console7Channel->fpdL ^= console7Channel->fpdL << 5;

  845. -		inputSampleL += (((double)console7Channel->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  846. -		frexpf((float)inputSampleR, &expon);

  847. -		console7Channel->fpdR ^= console7Channel->fpdR << 13;

  848. -		console7Channel->fpdR ^= console7Channel->fpdR >> 17;

  849. -		console7Channel->fpdR ^= console7Channel->fpdR << 5;

  850. -		inputSampleR += (((double)console7Channel->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  851. -		//end 32 bit stereo floating point dither

  852. -

  853.  		*out1 = (float) inputSampleL;

  854.  		*out2 = (float) inputSampleR;

  855. 

  856. diff --git a/src/Console7Crunch/Console7Crunch.c b/src/Console7Crunch/Console7Crunch.c

  857. index 8f3049e..7d51bb7 100644

  858. --- a/src/Console7Crunch/Console7Crunch.c

  859. +++ b/src/Console7Crunch/Console7Crunch.c

  860. @@ -3,6 +3,7 @@

  861.  #include <math.h>

  862.  #include <stdint.h>

  863.  #include <stdlib.h>

  864. +#include "../optimizations.h"

  865. 

  866.  #define CONSOLE7CRUNCH_URI "https://hannesbraun.net/ns/lv2/airwindows/console7crunch"

  867. 

  868. @@ -28,14 +29,11 @@ typedef struct {

  869.  	double biquadC[15];

  870.  	double biquadD[15];

  871.  	double biquadE[15];

  872. -

  873. -	uint32_t fpdL;

  874. -	uint32_t fpdR;

  875.  } Console7Crunch;

  876. 

  877.  static LV2_Handle instantiate(

  878.  	const LV2_Descriptor* descriptor,

  879. -	double rate,

  880. +	realdouble rate,

  881.  	const char* bundle_path,

  882.  	const LV2_Feature* const* features)

  883.  {

  884. @@ -80,10 +78,6 @@ static void activate(LV2_Handle instance)

  885.  		console7Crunch->biquadD[x] = 0.0;

  886.  		console7Crunch->biquadE[x] = 0.0;

  887.  	}

  888. -	console7Crunch->fpdL = 1.0;

  889. -	while (console7Crunch->fpdL < 16386) console7Crunch->fpdL = rand() * UINT32_MAX;

  890. -	console7Crunch->fpdR = 1.0;

  891. -	while (console7Crunch->fpdR < 16386) console7Crunch->fpdR = rand() * UINT32_MAX;

  892.  }

  893. 

  894.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  895. @@ -150,8 +144,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  896.  	while (sampleFrames-- > 0) {

  897.  		double inputSampleL = *in1;

  898.  		double inputSampleR = *in2;

  899. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = console7Crunch->fpdL * 1.18e-17;

  900. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = console7Crunch->fpdR * 1.18e-17;

  901. 

  902.  		double outSampleL = console7Crunch->biquadA[2] * inputSampleL + console7Crunch->biquadA[3] * console7Crunch->biquadA[7] + console7Crunch->biquadA[4] * console7Crunch->biquadA[8] - console7Crunch->biquadA[5] * console7Crunch->biquadA[9] - console7Crunch->biquadA[6] * console7Crunch->biquadA[10];

  903.  		console7Crunch->biquadA[8] = console7Crunch->biquadA[7];

  904. @@ -292,20 +284,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  905.  		}

  906.  		//we re-amplify after the distortion relative to how much we cut back previously.

  907. 

  908. -		//begin 32 bit stereo floating point dither

  909. -		int expon;

  910. -		frexpf((float)inputSampleL, &expon);

  911. -		console7Crunch->fpdL ^= console7Crunch->fpdL << 13;

  912. -		console7Crunch->fpdL ^= console7Crunch->fpdL >> 17;

  913. -		console7Crunch->fpdL ^= console7Crunch->fpdL << 5;

  914. -		inputSampleL += (((double)console7Crunch->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  915. -		frexpf((float)inputSampleR, &expon);

  916. -		console7Crunch->fpdR ^= console7Crunch->fpdR << 13;

  917. -		console7Crunch->fpdR ^= console7Crunch->fpdR >> 17;

  918. -		console7Crunch->fpdR ^= console7Crunch->fpdR << 5;

  919. -		inputSampleR += (((double)console7Crunch->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  920. -		//end 32 bit stereo floating point dither

  921. -

  922.  		*out1 = (float) inputSampleL;

  923.  		*out2 = (float) inputSampleR;

  924. 

  925. diff --git a/src/DCVoltage/DCVoltage.c b/src/DCVoltage/DCVoltage.c

  926. index 18f904e..8d17086 100644

  927. --- a/src/DCVoltage/DCVoltage.c

  928. +++ b/src/DCVoltage/DCVoltage.c

  929. @@ -2,6 +2,7 @@

  930. 

  931.  #include <stdint.h>

  932.  #include <stdlib.h>

  933. +#include "../optimizations.h"

  934. 

  935.  #define DCVOLTAGE_URI "https://hannesbraun.net/ns/lv2/airwindows/dcvoltage"

  936. 

  937. @@ -21,7 +22,7 @@ typedef struct {

  938. 

  939.  static LV2_Handle instantiate(

  940.  	const LV2_Descriptor* descriptor,

  941. -	double rate,

  942. +	realdouble rate,

  943.  	const char* bundle_path,

  944.  	const LV2_Feature* const* features)

  945.  {

  946. diff --git a/src/DeBess/DeBess.c b/src/DeBess/DeBess.c

  947. index 54a84d6..85ec475 100644

  948. --- a/src/DeBess/DeBess.c

  949. +++ b/src/DeBess/DeBess.c

  950. @@ -4,6 +4,7 @@

  951.  #include <stdbool.h>

  952.  #include <stdint.h>

  953.  #include <stdlib.h>

  954. +#include "../optimizations.h"

  955. 

  956.  #define DEBESS_URI "https://hannesbraun.net/ns/lv2/airwindows/debess"

  957. 

  958. @@ -43,13 +44,11 @@ typedef struct {

  959. 

  960. 

  961.  	bool flip;

  962. -	uint32_t fpdL;

  963. -	uint32_t fpdR;

  964.  } DeBess;

  965. 

  966.  static LV2_Handle instantiate(

  967.  	const LV2_Descriptor* descriptor,

  968. -	double rate,

  969. +	realdouble rate,

  970.  	const char* bundle_path,

  971.  	const LV2_Feature* const* features)

  972.  {

  973. @@ -115,10 +114,6 @@ static void activate(LV2_Handle instance)

  974.  	debess->iirSampleBR = 0.0;

  975. 

  976.  	debess->flip = false;

  977. -	debess->fpdL = 1.0;

  978. -	while (debess->fpdL < 16386) debess->fpdL = rand() * UINT32_MAX;

  979. -	debess->fpdR = 1.0;

  980. -	while (debess->fpdR < 16386) debess->fpdR = rand() * UINT32_MAX;

  981.  }

  982. 

  983.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  984. @@ -145,8 +140,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  985.  	while (sampleFrames-- > 0) {

  986.  		double inputSampleL = *in1;

  987.  		double inputSampleR = *in2;

  988. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = debess->fpdL * 1.18e-17;

  989. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = debess->fpdR * 1.18e-17;

  990. 

  991.  		debess->sL[0] = inputSampleL; //set up so both [0] and [1] will be input sample

  992.  		debess->sR[0] = inputSampleR; //set up so both [0] and [1] will be input sample

  993. @@ -208,20 +201,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  994.  		}

  995.  		//sense monitoring

  996. 

  997. -		//begin 32 bit stereo floating point dither

  998. -		int expon;

  999. -		frexpf((float)inputSampleL, &expon);

  1000. -		debess->fpdL ^= debess->fpdL << 13;

  1001. -		debess->fpdL ^= debess->fpdL >> 17;

  1002. -		debess->fpdL ^= debess->fpdL << 5;

  1003. -		inputSampleL += (((double)debess->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1004. -		frexpf((float)inputSampleR, &expon);

  1005. -		debess->fpdR ^= debess->fpdR << 13;

  1006. -		debess->fpdR ^= debess->fpdR >> 17;

  1007. -		debess->fpdR ^= debess->fpdR << 5;

  1008. -		inputSampleR += (((double)debess->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1009. -		//end 32 bit stereo floating point dither

  1010. -

  1011.  		*out1 = (float) inputSampleL;

  1012.  		*out2 = (float) inputSampleR;

  1013. 

  1014. diff --git a/src/Dyno/Dyno.c b/src/Dyno/Dyno.c

  1015. index ed028dd..cbbdf58 100644

  1016. --- a/src/Dyno/Dyno.c

  1017. +++ b/src/Dyno/Dyno.c

  1018. @@ -3,6 +3,7 @@

  1019.  #include <math.h>

  1020.  #include <stdint.h>

  1021.  #include <stdlib.h>

  1022. +#include "../optimizations.h"

  1023. 

  1024.  #define DYNO_URI "https://hannesbraun.net/ns/lv2/airwindows/dyno"

  1025. 

  1026. @@ -18,14 +19,11 @@ typedef struct {

  1027.  	const float* input[2];

  1028.  	float* output[2];

  1029.  	const float* gain;

  1030. -

  1031. -	uint32_t fpdL;

  1032. -	uint32_t fpdR;

  1033.  } Dyno;

  1034. 

  1035.  static LV2_Handle instantiate(

  1036.  	const LV2_Descriptor* descriptor,

  1037. -	double rate,

  1038. +	realdouble rate,

  1039.  	const char* bundle_path,

  1040.  	const LV2_Feature* const* features)

  1041.  {

  1042. @@ -56,14 +54,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1043.  	}

  1044.  }

  1045. 

  1046. -static void activate(LV2_Handle instance)

  1047. -{

  1048. -	Dyno* dyno = (Dyno*) instance;

  1049. -	dyno->fpdL = 1.0;

  1050. -	while (dyno->fpdL < 16386) dyno->fpdL = rand() * UINT32_MAX;

  1051. -	dyno->fpdR = 1.0;

  1052. -	while (dyno->fpdR < 16386) dyno->fpdR = rand() * UINT32_MAX;

  1053. -}

  1054. +static void activate(LV2_Handle instance) {}

  1055. 

  1056.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1057.  {

  1058. @@ -79,8 +70,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1059.  	while (sampleFrames-- > 0) {

  1060.  		double inputSampleL = *in1;

  1061.  		double inputSampleR = *in2;

  1062. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = dyno_instance->fpdL * 1.18e-17;

  1063. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = dyno_instance->fpdR * 1.18e-17;

  1064. 

  1065.  		if (gain != 1.0) {

  1066.  			inputSampleL *= gain;

  1067. @@ -94,20 +83,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1068.  		if (dyno > 0.0) inputSampleR = (sin(inputSampleR * dyno) / dyno) * 1.1654321;

  1069.  		//dyno is the one that tries to raise peak energy

  1070. 

  1071. -		//begin 32 bit stereo floating point dither

  1072. -		int expon;

  1073. -		frexpf((float)inputSampleL, &expon);

  1074. -		dyno_instance->fpdL ^= dyno_instance->fpdL << 13;

  1075. -		dyno_instance->fpdL ^= dyno_instance->fpdL >> 17;

  1076. -		dyno_instance->fpdL ^= dyno_instance->fpdL << 5;

  1077. -		inputSampleL += (((double)dyno_instance->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1078. -		frexpf((float)inputSampleR, &expon);

  1079. -		dyno_instance->fpdR ^= dyno_instance->fpdR << 13;

  1080. -		dyno_instance->fpdR ^= dyno_instance->fpdR >> 17;

  1081. -		dyno_instance->fpdR ^= dyno_instance->fpdR << 5;

  1082. -		inputSampleR += (((double)dyno_instance->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1083. -		//end 32 bit stereo floating point dither

  1084. -

  1085.  		*out1 = (float) inputSampleL;

  1086.  		*out2 = (float) inputSampleR;

  1087. 

  1088. diff --git a/src/EdIsDim/EdIsDim.c b/src/EdIsDim/EdIsDim.c

  1089. index 9fb5f2d..3a2d982 100644

  1090. --- a/src/EdIsDim/EdIsDim.c

  1091. +++ b/src/EdIsDim/EdIsDim.c

  1092. @@ -3,6 +3,7 @@

  1093.  #include <math.h>

  1094.  #include <stdint.h>

  1095.  #include <stdlib.h>

  1096. +#include "../optimizations.h"

  1097. 

  1098.  #define EDISDIM_URI "https://hannesbraun.net/ns/lv2/airwindows/edisdim"

  1099. 

  1100. @@ -18,14 +19,11 @@ typedef struct {

  1101.  	const float* input[2];

  1102.  	float* output[2];

  1103.  	const float* edisdim;

  1104. -

  1105. -	uint32_t fpdL;

  1106. -	uint32_t fpdR;

  1107.  } EdIsDim;

  1108. 

  1109.  static LV2_Handle instantiate(

  1110.  	const LV2_Descriptor* descriptor,

  1111. -	double rate,

  1112. +	realdouble rate,

  1113.  	const char* bundle_path,

  1114.  	const LV2_Feature* const* features)

  1115.  {

  1116. @@ -56,14 +54,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1117.  	}

  1118.  }

  1119. 

  1120. -static void activate(LV2_Handle instance)

  1121. -{

  1122. -	EdIsDim* edisdim = (EdIsDim*) instance;

  1123. -	edisdim->fpdL = 1.0;

  1124. -	while (edisdim->fpdL < 16386) edisdim->fpdL = rand() * UINT32_MAX;

  1125. -	edisdim->fpdR = 1.0;

  1126. -	while (edisdim->fpdR < 16386) edisdim->fpdR = rand() * UINT32_MAX;

  1127. -}

  1128. +static void activate(LV2_Handle instance) {}

  1129. 

  1130.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1131.  {

  1132. @@ -85,8 +76,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1133.  	while (sampleFrames-- > 0) {

  1134.  		inputSampleL = *in1;

  1135.  		inputSampleR = *in2;

  1136. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = edisdim->fpdL * 1.18e-17;

  1137. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = edisdim->fpdR * 1.18e-17;

  1138. 

  1139.  		inputSampleL *= midgain;

  1140.  		inputSampleR *= sidegain;

  1141. @@ -94,22 +83,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1142.  		mid = (inputSampleL + inputSampleR) / 2.0;

  1143.  		side = (inputSampleL - inputSampleR) / 2.0;

  1144. 

  1145. -		//begin 32 bit floating point dither

  1146. -		int expon;

  1147. -		frexpf((float)mid, &expon);

  1148. -		edisdim->fpdL ^= edisdim->fpdL << 13;

  1149. -		edisdim->fpdL ^= edisdim->fpdL >> 17;

  1150. -		edisdim->fpdL ^= edisdim->fpdL << 5;

  1151. -		mid += (((double)edisdim->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1152. -		//end 32 bit floating point dither

  1153. -		//begin 32 bit floating point dither

  1154. -		frexpf((float)side, &expon);

  1155. -		edisdim->fpdR ^= edisdim->fpdR << 13;

  1156. -		edisdim->fpdR ^= edisdim->fpdR >> 17;

  1157. -		edisdim->fpdR ^= edisdim->fpdR << 5;

  1158. -		side += (((double)edisdim->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1159. -		//end 32 bit floating point dither

  1160. -

  1161.  		*out1 = (float) mid;

  1162.  		*out2 = (float) side;

  1163. 

  1164. diff --git a/src/EveryTrim/EveryTrim.c b/src/EveryTrim/EveryTrim.c

  1165. index d22e1b9..58f8c57 100644

  1166. --- a/src/EveryTrim/EveryTrim.c

  1167. +++ b/src/EveryTrim/EveryTrim.c

  1168. @@ -3,6 +3,7 @@

  1169.  #include <math.h>

  1170.  #include <stdint.h>

  1171.  #include <stdlib.h>

  1172. +#include "../optimizations.h"

  1173. 

  1174.  #define EVERYTRIM_URI "https://hannesbraun.net/ns/lv2/airwindows/everytrim"

  1175. 

  1176. @@ -26,14 +27,11 @@ typedef struct {

  1177.  	const float* mid;

  1178.  	const float* side;

  1179.  	const float* master;

  1180. -

  1181. -	uint32_t fpdL;

  1182. -	uint32_t fpdR;

  1183.  } EveryTrim;

  1184. 

  1185.  static LV2_Handle instantiate(

  1186.  	const LV2_Descriptor* descriptor,

  1187. -	double rate,

  1188. +	realdouble rate,

  1189.  	const char* bundle_path,

  1190.  	const LV2_Feature* const* features)

  1191.  {

  1192. @@ -76,14 +74,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1193.  	}

  1194.  }

  1195. 

  1196. -static void activate(LV2_Handle instance)

  1197. -{

  1198. -	EveryTrim* everytrim = (EveryTrim*) instance;

  1199. -	everytrim->fpdL = 1.0;

  1200. -	while (everytrim->fpdL < 16386) everytrim->fpdL = rand() * UINT32_MAX;

  1201. -	everytrim->fpdR = 1.0;

  1202. -	while (everytrim->fpdR < 16386) everytrim->fpdR = rand() * UINT32_MAX;

  1203. -}

  1204. +static void activate(LV2_Handle instance) {}

  1205. 

  1206.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1207.  {

  1208. @@ -111,8 +102,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1209.  	while (sampleFrames-- > 0) {

  1210.  		inputSampleL = *in1;

  1211.  		inputSampleR = *in2;

  1212. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = everytrim->fpdL * 1.18e-17;

  1213. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = everytrim->fpdR * 1.18e-17;

  1214. 

  1215.  		mid = inputSampleL + inputSampleR;

  1216.  		side = inputSampleL - inputSampleR;

  1217. @@ -122,20 +111,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1218.  		inputSampleR = (mid - side) * rightgain;

  1219.  		//contains mastergain and the gain trim fixing the mid/side

  1220. 

  1221. -		//begin 32 bit stereo floating point dither

  1222. -		int expon;

  1223. -		frexpf((float)inputSampleL, &expon);

  1224. -		everytrim->fpdL ^= everytrim->fpdL << 13;

  1225. -		everytrim->fpdL ^= everytrim->fpdL >> 17;

  1226. -		everytrim->fpdL ^= everytrim->fpdL << 5;

  1227. -		inputSampleL += (((double)everytrim->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1228. -		frexpf((float)inputSampleR, &expon);

  1229. -		everytrim->fpdR ^= everytrim->fpdR << 13;

  1230. -		everytrim->fpdR ^= everytrim->fpdR >> 17;

  1231. -		everytrim->fpdR ^= everytrim->fpdR << 5;

  1232. -		inputSampleR += (((double)everytrim->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1233. -		//end 32 bit stereo floating point dither

  1234. -

  1235.  		*out1 = (float) inputSampleL;

  1236.  		*out2 = (float) inputSampleR;

  1237. 

  1238. diff --git a/src/Galactic/Galactic.c b/src/Galactic/Galactic.c

  1239. index 4629372..e38e881 100644

  1240. --- a/src/Galactic/Galactic.c

  1241. +++ b/src/Galactic/Galactic.c

  1242. @@ -3,6 +3,7 @@

  1243.  #include <math.h>

  1244.  #include <stdint.h>

  1245.  #include <stdlib.h>

  1246. +#include "../optimizations.h"

  1247. 

  1248.  #define GALACTIC_URI "https://hannesbraun.net/ns/lv2/airwindows/galactic"

  1249. 

  1250. @@ -48,7 +49,7 @@ typedef struct {

  1251. 

  1252.  	double aML[3111];

  1253.  	double aMR[3111];

  1254. -	double vibML, vibMR, depthM, oldfpd;

  1255. +	double vibML, vibMR, depthM;

  1256. 

  1257.  	double feedbackAL;

  1258.  	double feedbackBL;

  1259. @@ -101,14 +102,12 @@ typedef struct {

  1260. 

  1261.  	double vibM;

  1262. 

  1263. -	uint32_t fpdL;

  1264. -	uint32_t fpdR;

  1265.  	//default stuff

  1266.  } Galactic;

  1267. 

  1268.  static LV2_Handle instantiate(

  1269.  	const LV2_Descriptor* descriptor,

  1270. -	double rate,

  1271. +	realdouble rate,

  1272.  	const char* bundle_path,

  1273.  	const LV2_Feature* const* features)

  1274.  {

  1275. @@ -253,13 +252,6 @@ static void activate(LV2_Handle instance)

  1276.  	galactic->cycle = 0;

  1277. 

  1278.  	galactic->vibM = 3.0;

  1279. -

  1280. -	galactic->oldfpd = 429496.7295;

  1281. -

  1282. -	galactic->fpdL = 1;

  1283. -	while (galactic->fpdL < 16386) galactic->fpdL = rand() * UINT32_MAX;

  1284. -	galactic->fpdR = 1;

  1285. -	while (galactic->fpdR < 16386) galactic->fpdR = rand() * UINT32_MAX;

  1286.  }

  1287. 

  1288.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1289. @@ -308,15 +300,12 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1290.  	while (sampleFrames-- > 0) {

  1291.  		double inputSampleL = *in1;

  1292.  		double inputSampleR = *in2;

  1293. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = galactic->fpdL * 1.18e-17;

  1294. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = galactic->fpdR * 1.18e-17;

  1295.  		double drySampleL = inputSampleL;

  1296.  		double drySampleR = inputSampleR;

  1297. 

  1298. -		galactic->vibM += (galactic->oldfpd * drift);

  1299. +		galactic->vibM += drift;

  1300.  		if (galactic->vibM > (3.141592653589793238 * 2.0)) {

  1301.  			galactic->vibM = 0.0;

  1302. -			galactic->oldfpd = 0.4294967295 + (galactic->fpdL * 0.0000000000618);

  1303.  		}

  1304. 

  1305.  		galactic->aML[galactic->countM] = inputSampleL * attenuate;

  1306. @@ -496,20 +485,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1307.  			inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0 - wet));

  1308.  		}

  1309. 

  1310. -		//begin 32 bit stereo floating point dither

  1311. -		int expon;

  1312. -		frexpf((float)inputSampleL, &expon);

  1313. -		galactic->fpdL ^= galactic->fpdL << 13;

  1314. -		galactic->fpdL ^= galactic->fpdL >> 17;

  1315. -		galactic->fpdL ^= galactic->fpdL << 5;

  1316. -		inputSampleL += ((((double) galactic->fpdL) - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1317. -		frexpf((float)inputSampleR, &expon);

  1318. -		galactic->fpdR ^= galactic->fpdR << 13;

  1319. -		galactic->fpdR ^= galactic->fpdR >> 17;

  1320. -		galactic->fpdR ^= galactic->fpdR << 5;

  1321. -		inputSampleR += ((((double) galactic->fpdL) - (uint32_t) 0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1322. -		//end 32 bit stereo floating point dither

  1323. -

  1324.  		*out1 = (float) inputSampleL;

  1325.  		*out2 = (float) inputSampleR;

  1326. 

  1327. diff --git a/src/HermeTrim/HermeTrim.c b/src/HermeTrim/HermeTrim.c

  1328. index 1f2d728..fd754f4 100644

  1329. --- a/src/HermeTrim/HermeTrim.c

  1330. +++ b/src/HermeTrim/HermeTrim.c

  1331. @@ -3,6 +3,7 @@

  1332.  #include <math.h>

  1333.  #include <stdint.h>

  1334.  #include <stdlib.h>

  1335. +#include "../optimizations.h"

  1336. 

  1337.  #define HERMETRIM_URI "https://hannesbraun.net/ns/lv2/airwindows/hermetrim"

  1338. 

  1339. @@ -26,14 +27,11 @@ typedef struct {

  1340.  	const float* mid;

  1341.  	const float* side;

  1342.  	const float* master;

  1343. -

  1344. -	uint32_t fpdL;

  1345. -	uint32_t fpdR;

  1346.  } HermeTrim;

  1347. 

  1348.  static LV2_Handle instantiate(

  1349.  	const LV2_Descriptor* descriptor,

  1350. -	double rate,

  1351. +	realdouble rate,

  1352.  	const char* bundle_path,

  1353.  	const LV2_Feature* const* features)

  1354.  {

  1355. @@ -76,14 +74,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1356.  	}

  1357.  }

  1358. 

  1359. -static void activate(LV2_Handle instance)

  1360. -{

  1361. -	HermeTrim* hermetrim = (HermeTrim*) instance;

  1362. -	hermetrim->fpdL = 1.0;

  1363. -	while (hermetrim->fpdL < 16386) hermetrim->fpdL = rand() * UINT32_MAX;

  1364. -	hermetrim->fpdR = 1.0;

  1365. -	while (hermetrim->fpdR < 16386) hermetrim->fpdR = rand() * UINT32_MAX;

  1366. -}

  1367. +static void activate(LV2_Handle instance) {}

  1368. 

  1369.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1370.  {

  1371. @@ -111,8 +102,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1372.  	while (sampleFrames-- > 0) {

  1373.  		inputSampleL = *in1;

  1374.  		inputSampleR = *in2;

  1375. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = hermetrim->fpdL * 1.18e-17;

  1376. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = hermetrim->fpdR * 1.18e-17;

  1377. 

  1378.  		mid = inputSampleL + inputSampleR;

  1379.  		side = inputSampleL - inputSampleR;

  1380. @@ -122,20 +111,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1381.  		inputSampleR = (mid - side) * rightgain;

  1382.  		//contains mastergain and the gain trim fixing the mid/side

  1383. 

  1384. -		//begin 32 bit stereo floating point dither

  1385. -		int expon;

  1386. -		frexpf((float)inputSampleL, &expon);

  1387. -		hermetrim->fpdL ^= hermetrim->fpdL << 13;

  1388. -		hermetrim->fpdL ^= hermetrim->fpdL >> 17;

  1389. -		hermetrim->fpdL ^= hermetrim->fpdL << 5;

  1390. -		inputSampleL += (((double)hermetrim->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1391. -		frexpf((float)inputSampleR, &expon);

  1392. -		hermetrim->fpdR ^= hermetrim->fpdR << 13;

  1393. -		hermetrim->fpdR ^= hermetrim->fpdR >> 17;

  1394. -		hermetrim->fpdR ^= hermetrim->fpdR << 5;

  1395. -		inputSampleR += (((double)hermetrim->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1396. -		//end 32 bit stereo floating point dither

  1397. -

  1398.  		*out1 = (float) inputSampleL;

  1399.  		*out2 = (float) inputSampleR;

  1400. 

  1401. diff --git a/src/Hypersonic/Hypersonic.c b/src/Hypersonic/Hypersonic.c

  1402. index 68e8a7e..4f9daa7 100644

  1403. --- a/src/Hypersonic/Hypersonic.c

  1404. +++ b/src/Hypersonic/Hypersonic.c

  1405. @@ -3,6 +3,7 @@

  1406.  #include <math.h>

  1407.  #include <stdint.h>

  1408.  #include <stdlib.h>

  1409. +#include "../optimizations.h"

  1410. 

  1411.  #define HYPERSONIC_URI "https://hannesbraun.net/ns/lv2/airwindows/hypersonic"

  1412. 

  1413. @@ -40,15 +41,12 @@ typedef struct {

  1414.  	double fixE[fix_total];

  1415.  	double fixF[fix_total];

  1416.  	double fixG[fix_total];

  1417. -

  1418. -	uint32_t fpdL;

  1419. -	uint32_t fpdR;

  1420.  	//default stuff

  1421.  } Hypersonic;

  1422. 

  1423.  static LV2_Handle instantiate(

  1424.  	const LV2_Descriptor* descriptor,

  1425. -	double rate,

  1426. +	realdouble rate,

  1427.  	const char* bundle_path,

  1428.  	const LV2_Feature* const* features)

  1429.  {

  1430. @@ -89,10 +87,6 @@ static void activate(LV2_Handle instance)

  1431.  		hypersonic->fixF[x] = 0.0;

  1432.  		hypersonic->fixG[x] = 0.0;

  1433.  	}

  1434. -	hypersonic->fpdL = 1.0;

  1435. -	while (hypersonic->fpdL < 16386) hypersonic->fpdL = rand() * UINT32_MAX;

  1436. -	hypersonic->fpdR = 1.0;

  1437. -	while (hypersonic->fpdR < 16386) hypersonic->fpdR = rand() * UINT32_MAX;

  1438.  }

  1439. 

  1440.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1441. @@ -182,8 +176,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1442.  	while (sampleFrames-- > 0) {

  1443.  		double inputSampleL = *in1;

  1444.  		double inputSampleR = *in2;

  1445. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = hypersonic->fpdL * 1.18e-17;

  1446. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = hypersonic->fpdR * 1.18e-17;

  1447. 

  1448.  		double outSample = (inputSampleL * hypersonic->fixA[fix_a0]) + hypersonic->fixA[fix_sL1];

  1449.  		hypersonic->fixA[fix_sL1] = (inputSampleL * hypersonic->fixA[fix_a1]) - (outSample * hypersonic->fixA[fix_b1]) + hypersonic->fixA[fix_sL2];

  1450. @@ -248,20 +240,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1451.  		hypersonic->fixG[fix_sR2] = (inputSampleR * hypersonic->fixG[fix_a2]) - (outSample * hypersonic->fixG[fix_b2]);

  1452.  		inputSampleR = outSample; //fixed biquad filtering ultrasonics

  1453. 

  1454. -		//begin 32 bit stereo floating point dither

  1455. -		int expon;

  1456. -		frexpf((float)inputSampleL, &expon);

  1457. -		hypersonic->fpdL ^= hypersonic->fpdL << 13;

  1458. -		hypersonic->fpdL ^= hypersonic->fpdL >> 17;

  1459. -		hypersonic->fpdL ^= hypersonic->fpdL << 5;

  1460. -		inputSampleL += (((double)hypersonic->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1461. -		frexpf((float)inputSampleR, &expon);

  1462. -		hypersonic->fpdR ^= hypersonic->fpdR << 13;

  1463. -		hypersonic->fpdR ^= hypersonic->fpdR >> 17;

  1464. -		hypersonic->fpdR ^= hypersonic->fpdR << 5;

  1465. -		inputSampleR += (((double)hypersonic->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1466. -		//end 32 bit stereo floating point dither

  1467. -

  1468.  		*out1 = (float) inputSampleL;

  1469.  		*out2 = (float) inputSampleR;

  1470. 

  1471. diff --git a/src/Infrasonic/Infrasonic.c b/src/Infrasonic/Infrasonic.c

  1472. index 82c2d71..7288316 100644

  1473. --- a/src/Infrasonic/Infrasonic.c

  1474. +++ b/src/Infrasonic/Infrasonic.c

  1475. @@ -3,6 +3,7 @@

  1476.  #include <math.h>

  1477.  #include <stdint.h>

  1478.  #include <stdlib.h>

  1479. +#include "../optimizations.h"

  1480. 

  1481.  #define INFRASONIC_URI "https://hannesbraun.net/ns/lv2/airwindows/infrasonic"

  1482. 

  1483. @@ -22,13 +23,11 @@ typedef struct {

  1484.  	double biquadC[15];

  1485.  	double biquadD[15];

  1486.  	double biquadE[15];

  1487. -	uint32_t fpdL;

  1488. -	uint32_t fpdR;

  1489.  } Infrasonic;

  1490. 

  1491.  static LV2_Handle instantiate(

  1492.  	const LV2_Descriptor* descriptor,

  1493. -	double rate,

  1494. +	realdouble rate,

  1495.  	const char* bundle_path,

  1496.  	const LV2_Feature* const* features)

  1497.  {

  1498. @@ -66,10 +65,6 @@ static void activate(LV2_Handle instance)

  1499.  		infrasonic->biquadD[x] = 0.0;

  1500.  		infrasonic->biquadE[x] = 0.0;

  1501.  	}

  1502. -	infrasonic->fpdL = 1.0;

  1503. -	while (infrasonic->fpdL < 16386) infrasonic->fpdL = rand() * UINT32_MAX;

  1504. -	infrasonic->fpdR = 1.0;

  1505. -	while (infrasonic->fpdR < 16386) infrasonic->fpdR = rand() * UINT32_MAX;

  1506.  }

  1507. 

  1508.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1509. @@ -136,8 +131,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1510.  	while (sampleFrames-- > 0) {

  1511.  		double inputSampleL = *in1;

  1512.  		double inputSampleR = *in2;

  1513. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = infrasonic->fpdL * 1.18e-17;

  1514. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = infrasonic->fpdR * 1.18e-17;

  1515. 

  1516.  		double outSampleL = infrasonic->biquadA[2] * inputSampleL + infrasonic->biquadA[3] * infrasonic->biquadA[7] + infrasonic->biquadA[4] * infrasonic->biquadA[8] - infrasonic->biquadA[5] * infrasonic->biquadA[9] - infrasonic->biquadA[6] * infrasonic->biquadA[10];

  1517.  		infrasonic->biquadA[8] = infrasonic->biquadA[7];

  1518. @@ -209,20 +202,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1519.  		infrasonic->biquadE[14] = infrasonic->biquadE[13];

  1520.  		infrasonic->biquadE[13] = inputSampleR; //DF1 right

  1521. 

  1522. -		//begin 32 bit stereo floating point dither

  1523. -		int expon;

  1524. -		frexpf((float)inputSampleL, &expon);

  1525. -		infrasonic->fpdL ^= infrasonic->fpdL << 13;

  1526. -		infrasonic->fpdL ^= infrasonic->fpdL >> 17;

  1527. -		infrasonic->fpdL ^= infrasonic->fpdL << 5;

  1528. -		inputSampleL += (((double)infrasonic->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1529. -		frexpf((float)inputSampleR, &expon);

  1530. -		infrasonic->fpdR ^= infrasonic->fpdR << 13;

  1531. -		infrasonic->fpdR ^= infrasonic->fpdR >> 17;

  1532. -		infrasonic->fpdR ^= infrasonic->fpdR << 5;

  1533. -		inputSampleR += (((double)infrasonic->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1534. -		//end 32 bit stereo floating point dither

  1535. -

  1536.  		*out1 = (float) inputSampleL;

  1537.  		*out2 = (float) inputSampleR;

  1538. 

  1539. diff --git a/src/Interstage/Interstage.c b/src/Interstage/Interstage.c

  1540. index 885ef83..2d4df77 100644

  1541. --- a/src/Interstage/Interstage.c

  1542. +++ b/src/Interstage/Interstage.c

  1543. @@ -4,6 +4,7 @@

  1544.  #include <stdbool.h>

  1545.  #include <stdint.h>

  1546.  #include <stdlib.h>

  1547. +#include "../optimizations.h"

  1548. 

  1549.  #define INTERSTAGE_URI "https://hannesbraun.net/ns/lv2/airwindows/interstage"

  1550. 

  1551. @@ -33,14 +34,12 @@ typedef struct {

  1552.  	double iirSampleER;

  1553.  	double iirSampleFR;

  1554.  	double lastSampleR;

  1555. -	uint32_t fpdL;

  1556. -	uint32_t fpdR;

  1557.  	bool flip;

  1558.  } Interstage;

  1559. 

  1560.  static LV2_Handle instantiate(

  1561.  	const LV2_Descriptor* descriptor,

  1562. -	double rate,

  1563. +	realdouble rate,

  1564.  	const char* bundle_path,

  1565.  	const LV2_Feature* const* features)

  1566.  {

  1567. @@ -86,10 +85,6 @@ static void activate(LV2_Handle instance)

  1568.  	interstage->iirSampleER = 0.0;

  1569.  	interstage->iirSampleFR = 0.0;

  1570.  	interstage->lastSampleR = 0.0;

  1571. -	interstage->fpdL = 1.0;

  1572. -	while (interstage->fpdL < 16386) interstage->fpdL = rand() * UINT32_MAX;

  1573. -	interstage->fpdR = 1.0;

  1574. -	while (interstage->fpdR < 16386) interstage->fpdR = rand() * UINT32_MAX;

  1575.  	interstage->flip = true;

  1576.  }

  1577. 

  1578. @@ -113,8 +108,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1579.  	while (sampleFrames-- > 0) {

  1580.  		double inputSampleL = *in1;

  1581.  		double inputSampleR = *in2;

  1582. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = interstage->fpdL * 1.18e-17;

  1583. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = interstage->fpdR * 1.18e-17;

  1584.  		double drySampleL = inputSampleL;

  1585.  		double drySampleR = inputSampleR;

  1586. 

  1587. @@ -174,20 +167,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1588.  		interstage->lastSampleL = inputSampleL;

  1589.  		interstage->lastSampleR = inputSampleR;

  1590. 

  1591. -		//begin 32 bit stereo floating point dither

  1592. -		int expon;

  1593. -		frexpf((float)inputSampleL, &expon);

  1594. -		interstage->fpdL ^= interstage->fpdL << 13;

  1595. -		interstage->fpdL ^= interstage->fpdL >> 17;

  1596. -		interstage->fpdL ^= interstage->fpdL << 5;

  1597. -		inputSampleL += (((double)interstage->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1598. -		frexpf((float)inputSampleR, &expon);

  1599. -		interstage->fpdR ^= interstage->fpdR << 13;

  1600. -		interstage->fpdR ^= interstage->fpdR >> 17;

  1601. -		interstage->fpdR ^= interstage->fpdR << 5;

  1602. -		inputSampleR += (((double)interstage->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1603. -		//end 32 bit stereo floating point dither

  1604. -

  1605.  		*out1 = (float) inputSampleL;

  1606.  		*out2 = (float) inputSampleR;

  1607. 

  1608. diff --git a/src/LeftoMono/LeftoMono.c b/src/LeftoMono/LeftoMono.c

  1609. index bd0b49d..8fedca5 100644

  1610. --- a/src/LeftoMono/LeftoMono.c

  1611. +++ b/src/LeftoMono/LeftoMono.c

  1612. @@ -2,6 +2,7 @@

  1613. 

  1614.  #include <stdint.h>

  1615.  #include <stdlib.h>

  1616. +#include "../optimizations.h"

  1617. 

  1618.  #define LEFTOMONO_URI "https://hannesbraun.net/ns/lv2/airwindows/leftomono"

  1619. 

  1620. @@ -22,7 +23,7 @@ typedef struct {

  1621. 

  1622.  static LV2_Handle instantiate(

  1623.  	const LV2_Descriptor* descriptor,

  1624. -	double rate,

  1625. +	realdouble rate,

  1626.  	const char* bundle_path,

  1627.  	const LV2_Feature* const* features)

  1628.  {

  1629. diff --git a/src/MV/MV.c b/src/MV/MV.c

  1630. index 8d57a06..25f7b82 100644

  1631. --- a/src/MV/MV.c

  1632. +++ b/src/MV/MV.c

  1633. @@ -3,6 +3,7 @@

  1634.  #include <math.h>

  1635.  #include <stdint.h>

  1636.  #include <stdlib.h>

  1637. +#include "../optimizations.h"

  1638. 

  1639.  #define MV_URI "https://hannesbraun.net/ns/lv2/airwindows/mv"

  1640. 

  1641. @@ -27,9 +28,6 @@ typedef struct {

  1642.  	const float* gain;

  1643.  	const float* drywet;

  1644. 

  1645. -	uint32_t fpdL;

  1646. -	uint32_t fpdR;

  1647. -

  1648.  	double aAL[15150];

  1649.  	double aBL[14618];

  1650.  	double aCL[14358];

  1651. @@ -173,7 +171,7 @@ typedef struct {

  1652. 

  1653.  static LV2_Handle instantiate(

  1654.  	const LV2_Descriptor* descriptor,

  1655. -	double rate,

  1656. +	realdouble rate,

  1657.  	const char* bundle_path,

  1658.  	const LV2_Feature* const* features)

  1659.  {

  1660. @@ -433,11 +431,6 @@ static void activate(LV2_Handle instance)

  1661. 

  1662.  	mv->feedbackL = 0.0;

  1663.  	mv->feedbackR = 0.0;

  1664. -

  1665. -	mv->fpdL = 1.0;

  1666. -	while (mv->fpdL < 16386) mv->fpdL = rand() * UINT32_MAX;

  1667. -	mv->fpdR = 1.0;

  1668. -	while (mv->fpdR < 16386) mv->fpdR = rand() * UINT32_MAX;

  1669.  }

  1670. 

  1671.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1672. @@ -468,8 +461,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1673.  		double inputSampleL = *in1;

  1674.  		double inputSampleR = *in2;

  1675. 

  1676. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = mv->fpdL * 1.18e-17;

  1677. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = mv->fpdR * 1.18e-17;

  1678.  		double drySampleL = inputSampleL;

  1679.  		double drySampleR = inputSampleR;

  1680. 

  1681. @@ -1314,20 +1305,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1682.  		}

  1683.  		//Dry/Wet control, defaults to the last slider

  1684. 

  1685. -		//begin 32 bit stereo floating point dither

  1686. -		int expon;

  1687. -		frexpf((float)inputSampleL, &expon);

  1688. -		mv->fpdL ^= mv->fpdL << 13;

  1689. -		mv->fpdL ^= mv->fpdL >> 17;

  1690. -		mv->fpdL ^= mv->fpdL << 5;

  1691. -		inputSampleL += (((double)mv->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1692. -		frexpf((float)inputSampleR, &expon);

  1693. -		mv->fpdR ^= mv->fpdR << 13;

  1694. -		mv->fpdR ^= mv->fpdR >> 17;

  1695. -		mv->fpdR ^= mv->fpdR << 5;

  1696. -		inputSampleR += (((double)mv->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1697. -		//end 32 bit stereo floating point dither

  1698. -

  1699.  		*out1 = (float) inputSampleL;

  1700.  		*out2 = (float) inputSampleR;

  1701. 

  1702. diff --git a/src/Mackity/Mackity.c b/src/Mackity/Mackity.c

  1703. index 5250dd8..933d3a7 100644

  1704. --- a/src/Mackity/Mackity.c

  1705. +++ b/src/Mackity/Mackity.c

  1706. @@ -3,6 +3,7 @@

  1707.  #include <math.h>

  1708.  #include <stdint.h>

  1709.  #include <stdlib.h>

  1710. +#include "../optimizations.h"

  1711. 

  1712.  #define MACKITY_URI "https://hannesbraun.net/ns/lv2/airwindows/mackity"

  1713. 

  1714. @@ -28,14 +29,11 @@ typedef struct {

  1715.  	double iirSampleBR;

  1716.  	double biquadA[15];

  1717.  	double biquadB[15];

  1718. -

  1719. -	uint32_t fpdL;

  1720. -	uint32_t fpdR;

  1721.  } Mackity;

  1722. 

  1723.  static LV2_Handle instantiate(

  1724.  	const LV2_Descriptor* descriptor,

  1725. -	double rate,

  1726. +	realdouble rate,

  1727.  	const char* bundle_path,

  1728.  	const LV2_Feature* const* features)

  1729.  {

  1730. @@ -82,11 +80,6 @@ static void activate(LV2_Handle instance)

  1731.  		mackity->biquadA[x] = 0.0;

  1732.  		mackity->biquadB[x] = 0.0;

  1733.  	}

  1734. -

  1735. -	mackity->fpdL = 1.0;

  1736. -	while (mackity->fpdL < 16386) mackity->fpdL = rand() * UINT32_MAX;

  1737. -	mackity->fpdR = 1.0;

  1738. -	while (mackity->fpdR < 16386) mackity->fpdR = rand() * UINT32_MAX;

  1739.  }

  1740. 

  1741.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1742. @@ -132,8 +125,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1743.  	while (sampleFrames-- > 0) {

  1744.  		double inputSampleL = *in1;

  1745.  		double inputSampleR = *in2;

  1746. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = mackity->fpdL * 1.18e-17;

  1747. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = mackity->fpdR * 1.18e-17;

  1748. 

  1749.  		if (fabs(mackity->iirSampleAL) < 1.18e-37) mackity->iirSampleAL = 0.0;

  1750.  		mackity->iirSampleAL = (mackity->iirSampleAL * (1.0 - iirAmountA)) + (inputSampleL * iirAmountA);

  1751. @@ -194,20 +185,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1752.  			inputSampleR *= outPad;

  1753.  		}

  1754. 

  1755. -		//begin 32 bit stereo floating point dither

  1756. -		int expon;

  1757. -		frexpf((float)inputSampleL, &expon);

  1758. -		mackity->fpdL ^= mackity->fpdL << 13;

  1759. -		mackity->fpdL ^= mackity->fpdL >> 17;

  1760. -		mackity->fpdL ^= mackity->fpdL << 5;

  1761. -		inputSampleL += (((double)mackity->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1762. -		frexpf((float)inputSampleR, &expon);

  1763. -		mackity->fpdR ^= mackity->fpdR << 13;

  1764. -		mackity->fpdR ^= mackity->fpdR >> 17;

  1765. -		mackity->fpdR ^= mackity->fpdR << 5;

  1766. -		inputSampleR += (((double)mackity->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1767. -		//end 32 bit stereo floating point dither

  1768. -

  1769.  		*out1 = (float) inputSampleL;

  1770.  		*out2 = (float) inputSampleR;

  1771. 

  1772. diff --git a/src/MidSide/MidSide.c b/src/MidSide/MidSide.c

  1773. index 4c27ce7..7ee0c90 100644

  1774. --- a/src/MidSide/MidSide.c

  1775. +++ b/src/MidSide/MidSide.c

  1776. @@ -3,6 +3,7 @@

  1777.  #include <math.h>

  1778.  #include <stdint.h>

  1779.  #include <stdlib.h>

  1780. +#include "../optimizations.h"

  1781. 

  1782.  #define MIDSIDE_URI "https://hannesbraun.net/ns/lv2/airwindows/midside"

  1783. 

  1784. @@ -18,14 +19,11 @@ typedef struct {

  1785.  	const float* input[2];

  1786.  	float* output[2];

  1787.  	const float* midside;

  1788. -

  1789. -	uint32_t fpdL;

  1790. -	uint32_t fpdR;

  1791.  } MidSide;

  1792. 

  1793.  static LV2_Handle instantiate(

  1794.  	const LV2_Descriptor* descriptor,

  1795. -	double rate,

  1796. +	realdouble rate,

  1797.  	const char* bundle_path,

  1798.  	const LV2_Feature* const* features)

  1799.  {

  1800. @@ -56,14 +54,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1801.  	}

  1802.  }

  1803. 

  1804. -static void activate(LV2_Handle instance)

  1805. -{

  1806. -	MidSide* midside = (MidSide*) instance;

  1807. -	midside->fpdL = 1.0;

  1808. -	while (midside->fpdL < 16386) midside->fpdL = rand() * UINT32_MAX;

  1809. -	midside->fpdR = 1.0;

  1810. -	while (midside->fpdR < 16386) midside->fpdR = rand() * UINT32_MAX;

  1811. -}

  1812. +static void activate(LV2_Handle instance) {}

  1813. 

  1814.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1815.  {

  1816. @@ -86,8 +77,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1817.  	while (sampleFrames-- > 0) {

  1818.  		inputSampleL = *in1;

  1819.  		inputSampleR = *in2;

  1820. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = midside->fpdL * 1.18e-17;

  1821. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = midside->fpdR * 1.18e-17;

  1822. 

  1823.  		mid = inputSampleL + inputSampleR;

  1824.  		side = inputSampleL - inputSampleR;

  1825. @@ -95,21 +84,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1826.  		mid *= midgain;

  1827.  		side *= sidegain;

  1828. 

  1829. -		//begin 32 bit floating point dither

  1830. -		int expon;

  1831. -		frexpf((float)mid, &expon);

  1832. -		midside->fpdL ^= midside->fpdL << 13;

  1833. -		midside->fpdL ^= midside->fpdL >> 17;

  1834. -		midside->fpdL ^= midside->fpdL << 5;

  1835. -		mid += (((double)midside->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1836. -		//end 32 bit floating point dither

  1837. -		frexpf((float)side, &expon);

  1838. -		midside->fpdR ^= midside->fpdR << 13;

  1839. -		midside->fpdR ^= midside->fpdR >> 17;

  1840. -		midside->fpdR ^= midside->fpdR << 5;

  1841. -		side += (((double)midside->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1842. -		//end 32 bit floating point dither

  1843. -

  1844.  		*out1 = (float) mid;

  1845.  		*out2 = (float) side;

  1846. 

  1847. diff --git a/src/Mojo/Mojo.c b/src/Mojo/Mojo.c

  1848. index ee39090..d8ef7cf 100644

  1849. --- a/src/Mojo/Mojo.c

  1850. +++ b/src/Mojo/Mojo.c

  1851. @@ -3,6 +3,7 @@

  1852.  #include <math.h>

  1853.  #include <stdint.h>

  1854.  #include <stdlib.h>

  1855. +#include "../optimizations.h"

  1856. 

  1857.  #define MOJO_URI "https://hannesbraun.net/ns/lv2/airwindows/mojo"

  1858. 

  1859. @@ -18,14 +19,11 @@ typedef struct {

  1860.  	const float* input[2];

  1861.  	float* output[2];

  1862.  	const float* gain;

  1863. -

  1864. -	uint32_t fpdL;

  1865. -	uint32_t fpdR;

  1866.  } Mojo;

  1867. 

  1868.  static LV2_Handle instantiate(

  1869.  	const LV2_Descriptor* descriptor,

  1870. -	double rate,

  1871. +	realdouble rate,

  1872.  	const char* bundle_path,

  1873.  	const LV2_Feature* const* features)

  1874.  {

  1875. @@ -56,14 +54,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1876.  	}

  1877.  }

  1878. 

  1879. -static void activate(LV2_Handle instance)

  1880. -{

  1881. -	Mojo* mojo = (Mojo*) instance;

  1882. -	mojo->fpdL = 1.0;

  1883. -	while (mojo->fpdL < 16386) mojo->fpdL = rand() * UINT32_MAX;

  1884. -	mojo->fpdR = 1.0;

  1885. -	while (mojo->fpdR < 16386) mojo->fpdR = rand() * UINT32_MAX;

  1886. -}

  1887. +static void activate(LV2_Handle instance) {}

  1888. 

  1889.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1890.  {

  1891. @@ -79,8 +70,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1892.  	while (sampleFrames-- > 0) {

  1893.  		double inputSampleL = *in1;

  1894.  		double inputSampleR = *in2;

  1895. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = mojoInstance->fpdL * 1.18e-17;

  1896. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = mojoInstance->fpdR * 1.18e-17;

  1897. 

  1898.  		if (gain != 1.0) {

  1899.  			inputSampleL *= gain;

  1900. @@ -94,20 +83,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1901.  		if (mojo > 0.0) inputSampleR = (sin(inputSampleR * mojo * M_PI * 0.5) / mojo) * 0.987654321;

  1902.  		//mojo is the one that flattens WAAAAY out very softly before wavefolding

  1903. 

  1904. -		//begin 32 bit stereo floating point dither

  1905. -		int expon;

  1906. -		frexpf((float)inputSampleL, &expon);

  1907. -		mojoInstance->fpdL ^= mojoInstance->fpdL << 13;

  1908. -		mojoInstance->fpdL ^= mojoInstance->fpdL >> 17;

  1909. -		mojoInstance->fpdL ^= mojoInstance->fpdL << 5;

  1910. -		inputSampleL += (((double)mojoInstance->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1911. -		frexpf((float)inputSampleR, &expon);

  1912. -		mojoInstance->fpdR ^= mojoInstance->fpdR << 13;

  1913. -		mojoInstance->fpdR ^= mojoInstance->fpdR >> 17;

  1914. -		mojoInstance->fpdR ^= mojoInstance->fpdR << 5;

  1915. -		inputSampleR += (((double)mojoInstance->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1916. -		//end 32 bit stereo floating point dither

  1917. -

  1918.  		*out1 = (float) inputSampleL;

  1919.  		*out2 = (float) inputSampleR;

  1920. 

  1921. diff --git a/src/Nikola/Nikola.c b/src/Nikola/Nikola.c

  1922. index 772fbe2..0a5a274 100644

  1923. --- a/src/Nikola/Nikola.c

  1924. +++ b/src/Nikola/Nikola.c

  1925. @@ -4,6 +4,7 @@

  1926.  #include <stdbool.h>

  1927.  #include <stdint.h>

  1928.  #include <stdlib.h>

  1929. +#include "../optimizations.h"

  1930. 

  1931.  #define NIKOLA_URI "https://hannesbraun.net/ns/lv2/airwindows/nikola"

  1932. 

  1933. @@ -29,14 +30,11 @@ typedef struct {

  1934.  	int framenumberR;

  1935.  	bool wasNegativeR;

  1936.  	double outlevelR;

  1937. -

  1938. -	uint32_t fpdL;

  1939. -	uint32_t fpdR;

  1940.  } Nikola;

  1941. 

  1942.  static LV2_Handle instantiate(

  1943.  	const LV2_Descriptor* descriptor,

  1944. -	double rate,

  1945. +	realdouble rate,

  1946.  	const char* bundle_path,

  1947.  	const LV2_Feature* const* features)

  1948.  {

  1949. @@ -80,11 +78,6 @@ static void activate(LV2_Handle instance)

  1950.  	nikola->wasNegativeR = false;

  1951.  	nikola->outlevelR = 0.0;

  1952.  	nikola->framenumberR = 0;

  1953. -

  1954. -	nikola->fpdL = 1.0;

  1955. -	while (nikola->fpdL < 16386) nikola->fpdL = rand() * UINT32_MAX;

  1956. -	nikola->fpdR = 1.0;

  1957. -	while (nikola->fpdR < 16386) nikola->fpdR = rand() * UINT32_MAX;

  1958.  }

  1959. 

  1960.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  1961. @@ -103,8 +96,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1962.  	while (sampleFrames-- > 0) {

  1963.  		double inputSampleL = *in1;

  1964.  		double inputSampleR = *in2;

  1965. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = nikola->fpdL * 1.18e-17;

  1966. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = nikola->fpdR * 1.18e-17;

  1967.  		double drySampleL = inputSampleL;

  1968.  		double drySampleR = inputSampleR;

  1969. 

  1970. @@ -2965,20 +2956,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1971.  			inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0 - wet));

  1972.  		}

  1973. 

  1974. -		//begin 32 bit stereo floating point dither

  1975. -		int expon;

  1976. -		frexpf((float)inputSampleL, &expon);

  1977. -		nikola->fpdL ^= nikola->fpdL << 13;

  1978. -		nikola->fpdL ^= nikola->fpdL >> 17;

  1979. -		nikola->fpdL ^= nikola->fpdL << 5;

  1980. -		inputSampleL += (((double)nikola->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1981. -		frexpf((float)inputSampleR, &expon);

  1982. -		nikola->fpdR ^= nikola->fpdR << 13;

  1983. -		nikola->fpdR ^= nikola->fpdR >> 17;

  1984. -		nikola->fpdR ^= nikola->fpdR << 5;

  1985. -		inputSampleR += (((double)nikola->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  1986. -		//end 32 bit stereo floating point dither

  1987. -

  1988.  		*out1 = (float) inputSampleL;

  1989.  		*out2 = (float) inputSampleR;

  1990. 

  1991. diff --git a/src/PocketVerbs/PocketVerbs.c b/src/PocketVerbs/PocketVerbs.c

  1992. index 694321a..df4be1f 100644

  1993. --- a/src/PocketVerbs/PocketVerbs.c

  1994. +++ b/src/PocketVerbs/PocketVerbs.c

  1995. @@ -3,6 +3,7 @@

  1996.  #include <math.h>

  1997.  #include <stdint.h>

  1998.  #include <stdlib.h>

  1999. +#include "../optimizations.h"

  2000. 

  2001.  #define POCKETVERBS_URI "https://hannesbraun.net/ns/lv2/airwindows/pocketverbs"

  2002. 

  2003. @@ -245,14 +246,11 @@ typedef struct {

  2004.  	int countdown;

  2005.  	double peakL;

  2006.  	double peakR;

  2007. -

  2008. -	uint32_t fpdL;

  2009. -	uint32_t fpdR;

  2010.  } PocketVerbs;

  2011. 

  2012.  static LV2_Handle instantiate(

  2013.  	const LV2_Descriptor* descriptor,

  2014. -	double rate,

  2015. +	realdouble rate,

  2016.  	const char* bundle_path,

  2017.  	const LV2_Feature* const* features)

  2018.  {

  2019. @@ -831,10 +829,6 @@ static void activate(LV2_Handle instance)

  2020.  	pocketverbs->countdown = -1;

  2021.  	pocketverbs->peakL = 1.0;

  2022.  	pocketverbs->peakR = 1.0;

  2023. -	pocketverbs->fpdL = 1.0;

  2024. -	while (pocketverbs->fpdL < 16386) pocketverbs->fpdL = rand() * UINT32_MAX;

  2025. -	pocketverbs->fpdR = 1.0;

  2026. -	while (pocketverbs->fpdR < 16386) pocketverbs->fpdR = rand() * UINT32_MAX;

  2027.  }

  2028. 

  2029.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2030. @@ -1778,9 +1772,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2031.  		double inputSampleL = *in1;

  2032.  		double inputSampleR = *in2;

  2033. 

  2034. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = pocketverbs->fpdL * 1.18e-17;

  2035. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = pocketverbs->fpdR * 1.18e-17;

  2036. -

  2037.  		double drySampleL = inputSampleL;

  2038.  		double drySampleR = inputSampleR;

  2039. 

  2040. @@ -13590,20 +13581,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2041.  		inputSampleR += drySampleR;

  2042.  		//here we combine the tanks with the dry signal

  2043. 

  2044. -		//begin 32 bit stereo floating point dither

  2045. -		int expon;

  2046. -		frexpf((float)inputSampleL, &expon);

  2047. -		pocketverbs->fpdL ^= pocketverbs->fpdL << 13;

  2048. -		pocketverbs->fpdL ^= pocketverbs->fpdL >> 17;

  2049. -		pocketverbs->fpdL ^= pocketverbs->fpdL << 5;

  2050. -		inputSampleL += (((double)pocketverbs->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2051. -		frexpf((float)inputSampleR, &expon);

  2052. -		pocketverbs->fpdR ^= pocketverbs->fpdR << 13;

  2053. -		pocketverbs->fpdR ^= pocketverbs->fpdR >> 17;

  2054. -		pocketverbs->fpdR ^= pocketverbs->fpdR << 5;

  2055. -		inputSampleR += (((double)pocketverbs->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2056. -		//end 32 bit stereo floating point dither

  2057. -

  2058.  		*out1 = (float) inputSampleL;

  2059.  		*out2 = (float) inputSampleR;

  2060. 

  2061. diff --git a/src/Pressure5/Pressure5.c b/src/Pressure5/Pressure5.c

  2062. index 0ea5c79..e16c291 100644

  2063. --- a/src/Pressure5/Pressure5.c

  2064. +++ b/src/Pressure5/Pressure5.c

  2065. @@ -4,6 +4,7 @@

  2066.  #include <math.h>

  2067.  #include <stdint.h>

  2068.  #include <stdlib.h>

  2069. +#include "../optimizations.h"

  2070. 

  2071.  #define PRESSURE5_URI "https://hannesbraun.net/ns/lv2/airwindows/pressure5"

  2072. 

  2073. @@ -70,14 +71,11 @@ typedef struct {

  2074.  	bool wasNegClipR; //Stereo ClipOnly2

  2075. 

  2076.  	double slewMax; //to adust mewiness

  2077. -

  2078. -	uint32_t fpdL;

  2079. -	uint32_t fpdR;

  2080.  } Pressure5;

  2081. 

  2082.  static LV2_Handle instantiate(

  2083.  	const LV2_Descriptor* descriptor,

  2084. -	double rate,

  2085. +	realdouble rate,

  2086.  	const char* bundle_path,

  2087.  	const LV2_Feature* const* features)

  2088.  {

  2089. @@ -149,11 +147,6 @@ static void activate(LV2_Handle instance)

  2090.  		pressure5->intermediateR[x] = 0.0;

  2091.  	}

  2092.  	pressure5->slewMax = 0.0;

  2093. -

  2094. -	pressure5->fpdL = 1.0;

  2095. -	while (pressure5->fpdL < 16386) pressure5->fpdL = rand() * UINT32_MAX;

  2096. -	pressure5->fpdR = 1.0;

  2097. -	while (pressure5->fpdR < 16386) pressure5->fpdR = rand() * UINT32_MAX;

  2098.  }

  2099. 

  2100.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2101. @@ -202,8 +195,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2102.  	while (sampleFrames-- > 0) {

  2103.  		double inputSampleL = *in1;

  2104.  		double inputSampleR = *in2;

  2105. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = pressure5->fpdL * 1.18e-17;

  2106. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = pressure5->fpdR * 1.18e-17;

  2107.  		double drySampleL = inputSampleL;

  2108.  		double drySampleR = inputSampleR;

  2109. 

  2110. @@ -373,20 +364,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2111.  		//end ClipOnly2 stereo as a little, compressed chunk that can be dropped into code

  2112.  		//final clip runs AFTER the Dry/Wet. It serves as a safety clip even if you're not full wet

  2113. 

  2114. -		//begin 32 bit stereo floating point dither

  2115. -		int expon;

  2116. -		frexpf((float)inputSampleL, &expon);

  2117. -		pressure5->fpdL ^= pressure5->fpdL << 13;

  2118. -		pressure5->fpdL ^= pressure5->fpdL >> 17;

  2119. -		pressure5->fpdL ^= pressure5->fpdL << 5;

  2120. -		inputSampleL += (((double)pressure5->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2121. -		frexpf((float)inputSampleR, &expon);

  2122. -		pressure5->fpdR ^= pressure5->fpdR << 13;

  2123. -		pressure5->fpdR ^= pressure5->fpdR >> 17;

  2124. -		pressure5->fpdR ^= pressure5->fpdR << 5;

  2125. -		inputSampleR += (((double)pressure5->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2126. -		//end 32 bit stereo floating point dither

  2127. -

  2128.  		*out1 = (float) inputSampleL;

  2129.  		*out2 = (float) inputSampleR;

  2130. 

  2131. diff --git a/src/PurestConsoleBuss/PurestConsoleBuss.c b/src/PurestConsoleBuss/PurestConsoleBuss.c

  2132. index 663ff27..a8cae6c 100644

  2133. --- a/src/PurestConsoleBuss/PurestConsoleBuss.c

  2134. +++ b/src/PurestConsoleBuss/PurestConsoleBuss.c

  2135. @@ -3,6 +3,7 @@

  2136.  #include <math.h>

  2137.  #include <stdint.h>

  2138.  #include <stdlib.h>

  2139. +#include "../optimizations.h"

  2140. 

  2141.  #define PURESTCONSOLEBUSS_URI "https://hannesbraun.net/ns/lv2/airwindows/purestconsolebuss"

  2142. 

  2143. @@ -16,14 +17,11 @@ typedef enum {

  2144.  typedef struct {

  2145.  	const float* input[2];

  2146.  	float* output[2];

  2147. -

  2148. -	uint32_t fpdL;

  2149. -	uint32_t fpdR;

  2150.  } PurestConsoleBuss;

  2151. 

  2152.  static LV2_Handle instantiate(

  2153.  	const LV2_Descriptor* descriptor,

  2154. -	double rate,

  2155. +	realdouble rate,

  2156.  	const char* bundle_path,

  2157.  	const LV2_Feature* const* features)

  2158.  {

  2159. @@ -51,14 +49,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  2160.  	}

  2161.  }

  2162. 

  2163. -static void activate(LV2_Handle instance)

  2164. -{

  2165. -	PurestConsoleBuss* purestConsoleBuss = (PurestConsoleBuss*) instance;

  2166. -	purestConsoleBuss->fpdL = 1.0;

  2167. -	while (purestConsoleBuss->fpdL < 16386) purestConsoleBuss->fpdL = rand() * UINT32_MAX;

  2168. -	purestConsoleBuss->fpdR = 1.0;

  2169. -	while (purestConsoleBuss->fpdR < 16386) purestConsoleBuss->fpdR = rand() * UINT32_MAX;

  2170. -}

  2171. +static void activate(LV2_Handle instance) {}

  2172. 

  2173.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2174.  {

  2175. @@ -75,8 +66,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2176.  	while (sampleFrames-- > 0) {

  2177.  		inputSampleL = *in1;

  2178.  		inputSampleR = *in2;

  2179. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = purestConsoleBuss->fpdL * 1.18e-17;

  2180. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = purestConsoleBuss->fpdR * 1.18e-17;

  2181. 

  2182.  		if (inputSampleL > 1.0) inputSampleL = 1.0;

  2183.  		if (inputSampleL < -1.0) inputSampleL = -1.0;

  2184. @@ -88,20 +77,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2185.  		inputSampleR = asin(inputSampleR);

  2186.  		//amplitude aspect

  2187. 

  2188. -		//begin 32 bit stereo floating point dither

  2189. -		int expon;

  2190. -		frexpf((float)inputSampleL, &expon);

  2191. -		purestConsoleBuss->fpdL ^= purestConsoleBuss->fpdL << 13;

  2192. -		purestConsoleBuss->fpdL ^= purestConsoleBuss->fpdL >> 17;

  2193. -		purestConsoleBuss->fpdL ^= purestConsoleBuss->fpdL << 5;

  2194. -		inputSampleL += (((double)purestConsoleBuss->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2195. -		frexpf((float)inputSampleR, &expon);

  2196. -		purestConsoleBuss->fpdR ^= purestConsoleBuss->fpdR << 13;

  2197. -		purestConsoleBuss->fpdR ^= purestConsoleBuss->fpdR >> 17;

  2198. -		purestConsoleBuss->fpdR ^= purestConsoleBuss->fpdR << 5;

  2199. -		inputSampleR += (((double)purestConsoleBuss->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2200. -		//end 32 bit stereo floating point dither

  2201. -

  2202.  		*out1 = (float) inputSampleL;

  2203.  		*out2 = (float) inputSampleR;

  2204. 

  2205. diff --git a/src/PurestConsoleChannel/PurestConsoleChannel.c b/src/PurestConsoleChannel/PurestConsoleChannel.c

  2206. index cee3dcb..29f7187 100644

  2207. --- a/src/PurestConsoleChannel/PurestConsoleChannel.c

  2208. +++ b/src/PurestConsoleChannel/PurestConsoleChannel.c

  2209. @@ -3,6 +3,7 @@

  2210.  #include <math.h>

  2211.  #include <stdint.h>

  2212.  #include <stdlib.h>

  2213. +#include "../optimizations.h"

  2214. 

  2215.  #define PURESTCONSOLECHANNEL_URI "https://hannesbraun.net/ns/lv2/airwindows/purestconsolechannel"

  2216. 

  2217. @@ -16,14 +17,11 @@ typedef enum {

  2218.  typedef struct {

  2219.  	const float* input[2];

  2220.  	float* output[2];

  2221. -

  2222. -	uint32_t fpdL;

  2223. -	uint32_t fpdR;

  2224.  } PurestConsoleChannel;

  2225. 

  2226.  static LV2_Handle instantiate(

  2227.  	const LV2_Descriptor* descriptor,

  2228. -	double rate,

  2229. +	realdouble rate,

  2230.  	const char* bundle_path,

  2231.  	const LV2_Feature* const* features)

  2232.  {

  2233. @@ -51,14 +49,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  2234.  	}

  2235.  }

  2236. 

  2237. -static void activate(LV2_Handle instance)

  2238. -{

  2239. -	PurestConsoleChannel* purestConsoleChannel = (PurestConsoleChannel*) instance;

  2240. -	purestConsoleChannel->fpdL = 1.0;

  2241. -	while (purestConsoleChannel->fpdL < 16386) purestConsoleChannel->fpdL = rand() * UINT32_MAX;

  2242. -	purestConsoleChannel->fpdR = 1.0;

  2243. -	while (purestConsoleChannel->fpdR < 16386) purestConsoleChannel->fpdR = rand() * UINT32_MAX;

  2244. -}

  2245. +static void activate(LV2_Handle instance) {}

  2246. 

  2247.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2248.  {

  2249. @@ -75,27 +66,11 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2250.  	while (sampleFrames-- > 0) {

  2251.  		inputSampleL = *in1;

  2252.  		inputSampleR = *in2;

  2253. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = purestConsoleChannel->fpdL * 1.18e-17;

  2254. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = purestConsoleChannel->fpdR * 1.18e-17;

  2255. 

  2256.  		inputSampleL = sin(inputSampleL);

  2257.  		inputSampleR = sin(inputSampleR);

  2258.  		//amplitude aspect

  2259. 

  2260. -		//begin 32 bit stereo floating point dither

  2261. -		int expon;

  2262. -		frexpf((float)inputSampleL, &expon);

  2263. -		purestConsoleChannel->fpdL ^= purestConsoleChannel->fpdL << 13;

  2264. -		purestConsoleChannel->fpdL ^= purestConsoleChannel->fpdL >> 17;

  2265. -		purestConsoleChannel->fpdL ^= purestConsoleChannel->fpdL << 5;

  2266. -		inputSampleL += (((double)purestConsoleChannel->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2267. -		frexpf((float)inputSampleR, &expon);

  2268. -		purestConsoleChannel->fpdR ^= purestConsoleChannel->fpdR << 13;

  2269. -		purestConsoleChannel->fpdR ^= purestConsoleChannel->fpdR >> 17;

  2270. -		purestConsoleChannel->fpdR ^= purestConsoleChannel->fpdR << 5;

  2271. -		inputSampleR += (((double)purestConsoleChannel->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2272. -		//end 32 bit stereo floating point dither

  2273. -

  2274.  		*out1 = (float) inputSampleL;

  2275.  		*out2 = (float) inputSampleR;

  2276. 

  2277. diff --git a/src/PurestGain/PurestGain.c b/src/PurestGain/PurestGain.c

  2278. index 97c9925..b06c780 100644

  2279. --- a/src/PurestGain/PurestGain.c

  2280. +++ b/src/PurestGain/PurestGain.c

  2281. @@ -3,6 +3,7 @@

  2282.  #include <math.h>

  2283.  #include <stdint.h>

  2284.  #include <stdlib.h>

  2285. +#include "../optimizations.h"

  2286. 

  2287.  #define PURESTGAIN_URI "https://hannesbraun.net/ns/lv2/airwindows/purestgain"

  2288. 

  2289. @@ -22,8 +23,6 @@ typedef struct {

  2290.  	const float* gain;

  2291.  	const float* slowfade;

  2292. 

  2293. -	uint32_t fpdL;

  2294. -	uint32_t fpdR;

  2295.  	//default stuff

  2296.  	double gainchase;

  2297.  	double settingchase;

  2298. @@ -33,7 +32,7 @@ typedef struct {

  2299. 

  2300.  static LV2_Handle instantiate(

  2301.  	const LV2_Descriptor* descriptor,

  2302. -	double rate,

  2303. +	realdouble rate,

  2304.  	const char* bundle_path,

  2305.  	const LV2_Feature* const* features)

  2306.  {

  2307. @@ -75,10 +74,6 @@ static void activate(LV2_Handle instance)

  2308.  	purestGain->settingchase = -90.0;

  2309.  	purestGain->gainBchase = -90.0;

  2310.  	purestGain->chasespeed = 350.0;

  2311. -	purestGain->fpdL = 1.0;

  2312. -	while (purestGain->fpdL < 16386) purestGain->fpdL = rand() * UINT32_MAX;

  2313. -	purestGain->fpdR = 1.0;

  2314. -	while (purestGain->fpdR < 16386) purestGain->fpdR = rand() * UINT32_MAX;

  2315.  }

  2316. 

  2317.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2318. @@ -143,8 +138,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2319. 

  2320.  		inputSampleL = *in1;

  2321.  		inputSampleR = *in2;

  2322. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = purestGain->fpdL * 1.18e-17;

  2323. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = purestGain->fpdR * 1.18e-17;

  2324. 

  2325.  		if (1.0 == outputgain) {

  2326.  			*out1 = *in1;

  2327. @@ -152,19 +145,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2328.  		} else {

  2329.  			inputSampleL *= outputgain;

  2330.  			inputSampleR *= outputgain;

  2331. -			//begin 32 bit stereo floating point dither

  2332. -			int expon;

  2333. -			frexpf((float)inputSampleL, &expon);

  2334. -			purestGain->fpdL ^= purestGain->fpdL << 13;

  2335. -			purestGain->fpdL ^= purestGain->fpdL >> 17;

  2336. -			purestGain->fpdL ^= purestGain->fpdL << 5;

  2337. -			inputSampleL += (((double)purestGain->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2338. -			frexpf((float)inputSampleR, &expon);

  2339. -			purestGain->fpdR ^= purestGain->fpdR << 13;

  2340. -			purestGain->fpdR ^= purestGain->fpdR >> 17;

  2341. -			purestGain->fpdR ^= purestGain->fpdR << 5;

  2342. -			inputSampleR += (((double)purestGain->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2343. -			//end 32 bit stereo floating point dither

  2344.  			*out1 = (float) inputSampleL;

  2345.  			*out2 = (float) inputSampleR;

  2346.  		}

  2347. diff --git a/src/RightoMono/RightoMono.c b/src/RightoMono/RightoMono.c

  2348. index 03acec9..e6eb262 100644

  2349. --- a/src/RightoMono/RightoMono.c

  2350. +++ b/src/RightoMono/RightoMono.c

  2351. @@ -2,6 +2,7 @@

  2352. 

  2353.  #include <stdint.h>

  2354.  #include <stdlib.h>

  2355. +#include "../optimizations.h"

  2356. 

  2357.  #define RIGHTOMONO_URI "https://hannesbraun.net/ns/lv2/airwindows/rightomono"

  2358. 

  2359. @@ -22,7 +23,7 @@ typedef struct {

  2360. 

  2361.  static LV2_Handle instantiate(

  2362.  	const LV2_Descriptor* descriptor,

  2363. -	double rate,

  2364. +	realdouble rate,

  2365.  	const char* bundle_path,

  2366.  	const LV2_Feature* const* features)

  2367.  {

  2368. diff --git a/src/Sidepass/Sidepass.c b/src/Sidepass/Sidepass.c

  2369. index 0e99ac2..b55f35e 100644

  2370. --- a/src/Sidepass/Sidepass.c

  2371. +++ b/src/Sidepass/Sidepass.c

  2372. @@ -4,6 +4,7 @@

  2373.  #include <stdbool.h>

  2374.  #include <stdint.h>

  2375.  #include <stdlib.h>

  2376. +#include "../optimizations.h"

  2377. 

  2378.  #define SIDEPASS_URI "https://hannesbraun.net/ns/lv2/airwindows/sidepass"

  2379. 

  2380. @@ -24,13 +25,11 @@ typedef struct {

  2381.  	double iirSampleA;

  2382.  	double iirSampleB;

  2383.  	bool flip;

  2384. -	uint32_t fpdL;

  2385. -	uint32_t fpdR;

  2386.  } Sidepass;

  2387. 

  2388.  static LV2_Handle instantiate(

  2389.  	const LV2_Descriptor* descriptor,

  2390. -	double rate,

  2391. +	realdouble rate,

  2392.  	const char* bundle_path,

  2393.  	const LV2_Feature* const* features)

  2394.  {

  2395. @@ -68,10 +67,6 @@ static void activate(LV2_Handle instance)

  2396.  	sidepass->iirSampleA = 0.0;

  2397.  	sidepass->iirSampleB = 0.0;

  2398.  	sidepass->flip = true;

  2399. -	sidepass->fpdL = 1.0;

  2400. -	while (sidepass->fpdL < 16386) sidepass->fpdL = rand() * UINT32_MAX;

  2401. -	sidepass->fpdR = 1.0;

  2402. -	while (sidepass->fpdR < 16386) sidepass->fpdR = rand() * UINT32_MAX;

  2403.  }

  2404. 

  2405.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2406. @@ -95,8 +90,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2407.  	while (sampleFrames-- > 0) {

  2408.  		inputSampleL = *in1;

  2409.  		inputSampleR = *in2;

  2410. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = sidepass->fpdL * 1.18e-17;

  2411. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = sidepass->fpdR * 1.18e-17;

  2412. 

  2413.  		mid = inputSampleL + inputSampleR;

  2414.  		side = inputSampleL - inputSampleR;

  2415. @@ -113,20 +106,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2416.  		inputSampleL = (mid + side) / 2.0;

  2417.  		inputSampleR = (mid - side) / 2.0;

  2418. 

  2419. -		//begin 32 bit stereo floating point dither

  2420. -		int expon;

  2421. -		frexpf((float)inputSampleL, &expon);

  2422. -		sidepass->fpdL ^= sidepass->fpdL << 13;

  2423. -		sidepass->fpdL ^= sidepass->fpdL >> 17;

  2424. -		sidepass->fpdL ^= sidepass->fpdL << 5;

  2425. -		inputSampleL += (((double)sidepass->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2426. -		frexpf((float)inputSampleR, &expon);

  2427. -		sidepass->fpdR ^= sidepass->fpdR << 13;

  2428. -		sidepass->fpdR ^= sidepass->fpdR >> 17;

  2429. -		sidepass->fpdR ^= sidepass->fpdR << 5;

  2430. -		inputSampleR += (((double)sidepass->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2431. -		//end 32 bit stereo floating point dither

  2432. -

  2433.  		*out1 = (float) inputSampleL;

  2434.  		*out2 = (float) inputSampleR;

  2435. 

  2436. diff --git a/src/Slew/Slew.c b/src/Slew/Slew.c

  2437. index e23ad33..97fe6fb 100644

  2438. --- a/src/Slew/Slew.c

  2439. +++ b/src/Slew/Slew.c

  2440. @@ -3,6 +3,7 @@

  2441.  #include <math.h>

  2442.  #include <stdint.h>

  2443.  #include <stdlib.h>

  2444. +#include "../optimizations.h"

  2445. 

  2446.  #define SLEW_URI "https://hannesbraun.net/ns/lv2/airwindows/slew"

  2447. 

  2448. @@ -22,13 +23,11 @@ typedef struct {

  2449. 

  2450.  	double lastSampleL;

  2451.  	double lastSampleR;

  2452. -	uint32_t fpdL;

  2453. -	uint32_t fpdR;

  2454.  } Slew;

  2455. 

  2456.  static LV2_Handle instantiate(

  2457.  	const LV2_Descriptor* descriptor,

  2458. -	double rate,

  2459. +	realdouble rate,

  2460.  	const char* bundle_path,

  2461.  	const LV2_Feature* const* features)

  2462.  {

  2463. @@ -91,8 +90,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2464.  	while (sampleFrames-- > 0) {

  2465.  		inputSampleL = *in1;

  2466.  		inputSampleR = *in2;

  2467. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = slew->fpdL * 1.18e-17;

  2468. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = slew->fpdR * 1.18e-17;

  2469. 

  2470.  		clamp = inputSampleL - slew->lastSampleL;

  2471.  		outputSampleL = inputSampleL;

  2472. diff --git a/src/SlewOnly/SlewOnly.c b/src/SlewOnly/SlewOnly.c

  2473. index d2de114..2148b40 100644

  2474. --- a/src/SlewOnly/SlewOnly.c

  2475. +++ b/src/SlewOnly/SlewOnly.c

  2476. @@ -3,6 +3,7 @@

  2477.  #include <math.h>

  2478.  #include <stdint.h>

  2479.  #include <stdlib.h>

  2480. +#include "../optimizations.h"

  2481. 

  2482.  #define SLEWONLY_URI "https://hannesbraun.net/ns/lv2/airwindows/slewonly"

  2483. 

  2484. @@ -19,13 +20,11 @@ typedef struct {

  2485. 

  2486.  	double lastSampleL;

  2487.  	double lastSampleR;

  2488. -	uint32_t fpdL;

  2489. -	uint32_t fpdR;

  2490.  } SlewOnly;

  2491. 

  2492.  static LV2_Handle instantiate(

  2493.  	const LV2_Descriptor* descriptor,

  2494. -	double rate,

  2495. +	realdouble rate,

  2496.  	const char* bundle_path,

  2497.  	const LV2_Feature* const* features)

  2498.  {

  2499. @@ -78,8 +77,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2500.  	while (sampleFrames-- > 0) {

  2501.  		inputSampleL = *in1;

  2502.  		inputSampleR = *in2;

  2503. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = slewonly->fpdL * 1.18e-17;

  2504. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = slewonly->fpdR * 1.18e-17;

  2505. 

  2506.  		outputSampleL = (inputSampleL - slewonly->lastSampleL) * trim;

  2507.  		outputSampleR = (inputSampleR - slewonly->lastSampleR) * trim;

  2508. diff --git a/src/Spiral/Spiral.c b/src/Spiral/Spiral.c

  2509. index 05a3361..619b60a 100644

  2510. --- a/src/Spiral/Spiral.c

  2511. +++ b/src/Spiral/Spiral.c

  2512. @@ -3,6 +3,7 @@

  2513.  #include <math.h>

  2514.  #include <stdint.h>

  2515.  #include <stdlib.h>

  2516. +#include "../optimizations.h"

  2517. 

  2518.  #define SPIRAL_URI "https://hannesbraun.net/ns/lv2/airwindows/spiral"

  2519. 

  2520. @@ -16,14 +17,11 @@ typedef enum {

  2521.  typedef struct {

  2522.  	const float* input[2];

  2523.  	float* output[2];

  2524. -

  2525. -	uint32_t fpdL;

  2526. -	uint32_t fpdR;

  2527.  } Spiral;

  2528. 

  2529.  static LV2_Handle instantiate(

  2530.  	const LV2_Descriptor* descriptor,

  2531. -	double rate,

  2532. +	realdouble rate,

  2533.  	const char* bundle_path,

  2534.  	const LV2_Feature* const* features)

  2535.  {

  2536. @@ -51,14 +49,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  2537.  	}

  2538.  }

  2539. 

  2540. -static void activate(LV2_Handle instance)

  2541. -{

  2542. -	Spiral* spiral = (Spiral*) instance;

  2543. -	spiral->fpdL = 1.0;

  2544. -	while (spiral->fpdL < 16386) spiral->fpdL = rand() * UINT32_MAX;

  2545. -	spiral->fpdR = 1.0;

  2546. -	while (spiral->fpdR < 16386) spiral->fpdR = rand() * UINT32_MAX;

  2547. -}

  2548. +static void activate(LV2_Handle instance) {}

  2549. 

  2550.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2551.  {

  2552. @@ -73,28 +64,10 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2553.  		double inputSampleL = *in1;

  2554.  		double inputSampleR = *in2;

  2555. 

  2556. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = spiral->fpdL * 1.18e-17;

  2557. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = spiral->fpdR * 1.18e-17;

  2558. -

  2559. -

  2560.  		//clip to 1.2533141373155 to reach maximum output

  2561.  		inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((fabs(inputSampleL) == 0.0) ? 1 : fabs(inputSampleL));

  2562.  		inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((fabs(inputSampleR) == 0.0) ? 1 : fabs(inputSampleR));

  2563. 

  2564. -		//begin 32 bit stereo floating point dither

  2565. -		int expon;

  2566. -		frexpf((float)inputSampleL, &expon);

  2567. -		spiral->fpdL ^= spiral->fpdL << 13;

  2568. -		spiral->fpdL ^= spiral->fpdL >> 17;

  2569. -		spiral->fpdL ^= spiral->fpdL << 5;

  2570. -		inputSampleL += (((double)spiral->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2571. -		frexpf((float)inputSampleR, &expon);

  2572. -		spiral->fpdR ^= spiral->fpdR << 13;

  2573. -		spiral->fpdR ^= spiral->fpdR >> 17;

  2574. -		spiral->fpdR ^= spiral->fpdR << 5;

  2575. -		inputSampleR += (((double)spiral->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2576. -		//end 32 bit stereo floating point dither

  2577. -

  2578.  		*out1 = (float) inputSampleL;

  2579.  		*out2 = (float) inputSampleR;

  2580. 

  2581. diff --git a/src/StarChild/StarChild.c b/src/StarChild/StarChild.c

  2582. index c45f41d..e16f9ef 100644

  2583. --- a/src/StarChild/StarChild.c

  2584. +++ b/src/StarChild/StarChild.c

  2585. @@ -3,6 +3,7 @@

  2586.  #include <math.h>

  2587.  #include <stdint.h>

  2588.  #include <stdlib.h>

  2589. +#include "../optimizations.h"

  2590. 

  2591.  #define STARCHILD_URI "https://hannesbraun.net/ns/lv2/airwindows/starchild"

  2592. 

  2593. @@ -23,8 +24,6 @@ typedef struct {

  2594.  	const float* grain;

  2595.  	const float* drywet;

  2596. 

  2597. -	uint32_t fpdL;

  2598. -	uint32_t fpdR;

  2599.  	//default stuff

  2600. 

  2601.  	double d[45102];

  2602. @@ -47,7 +46,7 @@ typedef struct {

  2603. 

  2604.  static LV2_Handle instantiate(

  2605.  	const LV2_Descriptor* descriptor,

  2606. -	double rate,

  2607. +	realdouble rate,

  2608.  	const char* bundle_path,

  2609.  	const LV2_Feature* const* features)

  2610.  {

  2611. @@ -303,11 +302,6 @@ static void activate(LV2_Handle instance)

  2612.  	starchild->pitchCounter = 2;

  2613.  	starchild->increment = 1;

  2614.  	starchild->dutyCycle = 1;

  2615. -

  2616. -	starchild->fpdL = 1.0;

  2617. -	while (starchild->fpdL < 16386) starchild->fpdL = rand() * UINT32_MAX;

  2618. -	starchild->fpdR = 1.0;

  2619. -	while (starchild->fpdR < 16386) starchild->fpdR = rand() * UINT32_MAX;

  2620.  }

  2621. 

  2622.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2623. @@ -437,8 +431,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2624.  	while (sampleFrames-- > 0) {

  2625.  		inputSampleL = *in1;

  2626.  		inputSampleR = *in2;

  2627. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = starchild->fpdL * 1.18e-17;

  2628. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = starchild->fpdR * 1.18e-17;

  2629.  		drySampleL = inputSampleL;

  2630.  		drySampleR = inputSampleR;

  2631. 

  2632. @@ -1058,20 +1050,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2633.  		inputSampleR += drySampleR;

  2634.  		//here we combine the tanks with the dry signal

  2635. 

  2636. -		//begin 32 bit stereo floating point dither

  2637. -		int expon;

  2638. -		frexpf((float)inputSampleL, &expon);

  2639. -		starchild->fpdL ^= starchild->fpdL << 13;

  2640. -		starchild->fpdL ^= starchild->fpdL >> 17;

  2641. -		starchild->fpdL ^= starchild->fpdL << 5;

  2642. -		inputSampleL += (((double)starchild->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2643. -		frexpf((float)inputSampleR, &expon);

  2644. -		starchild->fpdR ^= starchild->fpdR << 13;

  2645. -		starchild->fpdR ^= starchild->fpdR >> 17;

  2646. -		starchild->fpdR ^= starchild->fpdR << 5;

  2647. -		inputSampleR += (((double)starchild->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2648. -		//end 32 bit stereo floating point dither

  2649. -

  2650.  		*out1 = (float) inputSampleL;

  2651.  		*out2 = (float) inputSampleR;

  2652. 

  2653. diff --git a/src/SubsOnly/SubsOnly.c b/src/SubsOnly/SubsOnly.c

  2654. index 886ae22..e649464 100644

  2655. --- a/src/SubsOnly/SubsOnly.c

  2656. +++ b/src/SubsOnly/SubsOnly.c

  2657. @@ -3,6 +3,7 @@

  2658.  #include <math.h>

  2659.  #include <stdint.h>

  2660.  #include <stdlib.h>

  2661. +#include "../optimizations.h"

  2662. 

  2663.  #define SUBSONLY_URI "https://hannesbraun.net/ns/lv2/airwindows/subsonly"

  2664. 

  2665. @@ -71,14 +72,11 @@ typedef struct {

  2666.  	double iirSampleXR;

  2667.  	double iirSampleYR;

  2668.  	double iirSampleZR;

  2669. -

  2670. -	uint32_t fpdL;

  2671. -	uint32_t fpdR;

  2672.  } SubsOnly;

  2673. 

  2674.  static LV2_Handle instantiate(

  2675.  	const LV2_Descriptor* descriptor,

  2676. -	double rate,

  2677. +	realdouble rate,

  2678.  	const char* bundle_path,

  2679.  	const LV2_Feature* const* features)

  2680.  {

  2681. @@ -188,8 +186,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2682.  	while (sampleFrames-- > 0) {

  2683.  		inputSampleL = *in1;

  2684.  		inputSampleR = *in2;

  2685. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = subsonly->fpdL * 1.18e-17;

  2686. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = subsonly->fpdR * 1.18e-17;

  2687. 

  2688.  		gain = gaintarget;

  2689. 

  2690. diff --git a/src/Ultrasonic/Ultrasonic.c b/src/Ultrasonic/Ultrasonic.c

  2691. index 57b2398..6aacb3b 100644

  2692. --- a/src/Ultrasonic/Ultrasonic.c

  2693. +++ b/src/Ultrasonic/Ultrasonic.c

  2694. @@ -3,6 +3,7 @@

  2695.  #include <math.h>

  2696.  #include <stdint.h>

  2697.  #include <stdlib.h>

  2698. +#include "../optimizations.h"

  2699. 

  2700.  #define ULTRASONIC_URI "https://hannesbraun.net/ns/lv2/airwindows/ultrasonic"

  2701. 

  2702. @@ -22,13 +23,11 @@ typedef struct {

  2703.  	double biquadC[15];

  2704.  	double biquadD[15];

  2705.  	double biquadE[15];

  2706. -	uint32_t fpdL;

  2707. -	uint32_t fpdR;

  2708.  } Ultrasonic;

  2709. 

  2710.  static LV2_Handle instantiate(

  2711.  	const LV2_Descriptor* descriptor,

  2712. -	double rate,

  2713. +	realdouble rate,

  2714.  	const char* bundle_path,

  2715.  	const LV2_Feature* const* features)

  2716.  {

  2717. @@ -66,10 +65,6 @@ static void activate(LV2_Handle instance)

  2718.  		ultrasonic->biquadD[x] = 0.0;

  2719.  		ultrasonic->biquadE[x] = 0.0;

  2720.  	}

  2721. -	ultrasonic->fpdL = 1.0;

  2722. -	while (ultrasonic->fpdL < 16386) ultrasonic->fpdL = rand() * UINT32_MAX;

  2723. -	ultrasonic->fpdR = 1.0;

  2724. -	while (ultrasonic->fpdR < 16386) ultrasonic->fpdR = rand() * UINT32_MAX;

  2725.  }

  2726. 

  2727.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2728. @@ -136,8 +131,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2729.  	while (sampleFrames-- > 0) {

  2730.  		double inputSampleL = *in1;

  2731.  		double inputSampleR = *in2;

  2732. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = ultrasonic->fpdL * 1.18e-17;

  2733. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = ultrasonic->fpdR * 1.18e-17;

  2734. 

  2735.  		double outSampleL = ultrasonic->biquadA[2] * inputSampleL + ultrasonic->biquadA[3] * ultrasonic->biquadA[7] + ultrasonic->biquadA[4] * ultrasonic->biquadA[8] - ultrasonic->biquadA[5] * ultrasonic->biquadA[9] - ultrasonic->biquadA[6] * ultrasonic->biquadA[10];

  2736.  		ultrasonic->biquadA[8] = ultrasonic->biquadA[7];

  2737. @@ -209,20 +202,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2738.  		ultrasonic->biquadE[14] = ultrasonic->biquadE[13];

  2739.  		ultrasonic->biquadE[13] = inputSampleR; // DF1 right

  2740. 

  2741. -		//begin 32 bit stereo floating point dither

  2742. -		int expon;

  2743. -		frexpf((float)inputSampleL, &expon);

  2744. -		ultrasonic->fpdL ^= ultrasonic->fpdL << 13;

  2745. -		ultrasonic->fpdL ^= ultrasonic->fpdL >> 17;

  2746. -		ultrasonic->fpdL ^= ultrasonic->fpdL << 5;

  2747. -		inputSampleL += (((double)ultrasonic->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2748. -		frexpf((float)inputSampleR, &expon);

  2749. -		ultrasonic->fpdR ^= ultrasonic->fpdR << 13;

  2750. -		ultrasonic->fpdR ^= ultrasonic->fpdR >> 17;

  2751. -		ultrasonic->fpdR ^= ultrasonic->fpdR << 5;

  2752. -		inputSampleR += (((double)ultrasonic->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2753. -		//end 32 bit stereo floating point dither

  2754. -

  2755.  		*out1 = (float) inputSampleL;

  2756.  		*out2 = (float) inputSampleR;

  2757. 

  2758. diff --git a/src/UltrasonicLite/UltrasonicLite.c b/src/UltrasonicLite/UltrasonicLite.c

  2759. index f73d3f5..ef932f2 100644

  2760. --- a/src/UltrasonicLite/UltrasonicLite.c

  2761. +++ b/src/UltrasonicLite/UltrasonicLite.c

  2762. @@ -3,6 +3,7 @@

  2763.  #include <math.h>

  2764.  #include <stdint.h>

  2765.  #include <stdlib.h>

  2766. +#include "../optimizations.h"

  2767. 

  2768.  #define ULTRASONICLITE_URI "https://hannesbraun.net/ns/lv2/airwindows/ultrasoniclite"

  2769. 

  2770. @@ -18,13 +19,11 @@ typedef struct {

  2771.  	const float* input[2];

  2772.  	float* output[2];

  2773.  	double biquadA[15];

  2774. -	uint32_t fpdL;

  2775. -	uint32_t fpdR;

  2776.  } UltrasonicLite;

  2777. 

  2778.  static LV2_Handle instantiate(

  2779.  	const LV2_Descriptor* descriptor,

  2780. -	double rate,

  2781. +	realdouble rate,

  2782.  	const char* bundle_path,

  2783.  	const LV2_Feature* const* features)

  2784.  {

  2785. @@ -58,10 +57,6 @@ static void activate(LV2_Handle instance)

  2786.  	for (int x = 0; x < 15; x++) {

  2787.  		ultrasoniclite->biquadA[x] = 0.0;

  2788.  	}

  2789. -	ultrasoniclite->fpdL = 1.0;

  2790. -	while (ultrasoniclite->fpdL < 16386) ultrasoniclite->fpdL = rand() * UINT32_MAX;

  2791. -	ultrasoniclite->fpdR = 1.0;

  2792. -	while (ultrasoniclite->fpdR < 16386) ultrasoniclite->fpdR = rand() * UINT32_MAX;

  2793.  }

  2794. 

  2795.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2796. @@ -90,8 +85,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2797.  	while (sampleFrames-- > 0) {

  2798.  		double inputSampleL = *in1;

  2799.  		double inputSampleR = *in2;

  2800. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = ultrasoniclite->fpdL * 1.18e-17;

  2801. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = ultrasoniclite->fpdR * 1.18e-17;

  2802. 

  2803.  		double outSampleL = ultrasoniclite->biquadA[2] * inputSampleL + ultrasoniclite->biquadA[3] * ultrasoniclite->biquadA[7] + ultrasoniclite->biquadA[4] * ultrasoniclite->biquadA[8] - ultrasoniclite->biquadA[5] * ultrasoniclite->biquadA[9] - ultrasoniclite->biquadA[6] * ultrasoniclite->biquadA[10];

  2804.  		ultrasoniclite->biquadA[8] = ultrasoniclite->biquadA[7];

  2805. @@ -107,20 +100,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2806.  		ultrasoniclite->biquadA[14] = ultrasoniclite->biquadA[13];

  2807.  		ultrasoniclite->biquadA[13] = inputSampleR; //DF1 right

  2808. 

  2809. -		//begin 32 bit stereo floating point dither

  2810. -		int expon;

  2811. -		frexpf((float)inputSampleL, &expon);

  2812. -		ultrasoniclite->fpdL ^= ultrasoniclite->fpdL << 13;

  2813. -		ultrasoniclite->fpdL ^= ultrasoniclite->fpdL >> 17;

  2814. -		ultrasoniclite->fpdL ^= ultrasoniclite->fpdL << 5;

  2815. -		inputSampleL += (((double)ultrasoniclite->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2816. -		frexpf((float)inputSampleR, &expon);

  2817. -		ultrasoniclite->fpdR ^= ultrasoniclite->fpdR << 13;

  2818. -		ultrasoniclite->fpdR ^= ultrasoniclite->fpdR >> 17;

  2819. -		ultrasoniclite->fpdR ^= ultrasoniclite->fpdR << 5;

  2820. -		inputSampleR += (((double)ultrasoniclite->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2821. -		//end 32 bit stereo floating point dither

  2822. -

  2823.  		*out1 = (float) inputSampleL;

  2824.  		*out2 = (float) inputSampleR;

  2825. 

  2826. diff --git a/src/UltrasonicMed/UltrasonicMed.c b/src/UltrasonicMed/UltrasonicMed.c

  2827. index 2231697..29d91be 100644

  2828. --- a/src/UltrasonicMed/UltrasonicMed.c

  2829. +++ b/src/UltrasonicMed/UltrasonicMed.c

  2830. @@ -3,6 +3,7 @@

  2831.  #include <math.h>

  2832.  #include <stdint.h>

  2833.  #include <stdlib.h>

  2834. +#include "../optimizations.h"

  2835. 

  2836.  #define ULTRASONICMED_URI "https://hannesbraun.net/ns/lv2/airwindows/ultrasonicmed"

  2837. 

  2838. @@ -19,13 +20,11 @@ typedef struct {

  2839.  	float* output[2];

  2840.  	double biquadA[15];

  2841.  	double biquadB[15];

  2842. -	uint32_t fpdL;

  2843. -	uint32_t fpdR;

  2844.  } UltrasonicMed;

  2845. 

  2846.  static LV2_Handle instantiate(

  2847.  	const LV2_Descriptor* descriptor,

  2848. -	double rate,

  2849. +	realdouble rate,

  2850.  	const char* bundle_path,

  2851.  	const LV2_Feature* const* features)

  2852.  {

  2853. @@ -61,10 +60,6 @@ static void activate(LV2_Handle instance)

  2854.  		ultrasonicMed->biquadA[x] = 0.0;

  2855.  		ultrasonicMed->biquadB[x] = 0.0;

  2856.  	}

  2857. -	ultrasonicMed->fpdL = 1.0;

  2858. -	while (ultrasonicMed->fpdL < 16386) ultrasonicMed->fpdL = rand() * UINT32_MAX;

  2859. -	ultrasonicMed->fpdR = 1.0;

  2860. -	while (ultrasonicMed->fpdR < 16386) ultrasonicMed->fpdR = rand() * UINT32_MAX;

  2861.  }

  2862. 

  2863.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2864. @@ -105,8 +100,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2865.  	while (sampleFrames-- > 0) {

  2866.  		double inputSampleL = *in1;

  2867.  		double inputSampleR = *in2;

  2868. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = ultrasonicMed->fpdL * 1.18e-17;

  2869. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = ultrasonicMed->fpdR * 1.18e-17;

  2870. 

  2871.  		double outSampleL = ultrasonicMed->biquadA[2] * inputSampleL + ultrasonicMed->biquadA[3] * ultrasonicMed->biquadA[7] + ultrasonicMed->biquadA[4] * ultrasonicMed->biquadA[8] - ultrasonicMed->biquadA[5] * ultrasonicMed->biquadA[9] - ultrasonicMed->biquadA[6] * ultrasonicMed->biquadA[10];

  2872.  		ultrasonicMed->biquadA[8] = ultrasonicMed->biquadA[7];

  2873. @@ -136,20 +129,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2874.  		ultrasonicMed->biquadB[14] = ultrasonicMed->biquadB[13];

  2875.  		ultrasonicMed->biquadB[13] = inputSampleR; //DF1 right

  2876. 

  2877. -		//begin 32 bit stereo floating point dither

  2878. -		int expon;

  2879. -		frexpf((float)inputSampleL, &expon);

  2880. -		ultrasonicMed->fpdL ^= ultrasonicMed->fpdL << 13;

  2881. -		ultrasonicMed->fpdL ^= ultrasonicMed->fpdL >> 17;

  2882. -		ultrasonicMed->fpdL ^= ultrasonicMed->fpdL << 5;

  2883. -		inputSampleL += (((double)ultrasonicMed->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2884. -		frexpf((float)inputSampleR, &expon);

  2885. -		ultrasonicMed->fpdR ^= ultrasonicMed->fpdR << 13;

  2886. -		ultrasonicMed->fpdR ^= ultrasonicMed->fpdR >> 17;

  2887. -		ultrasonicMed->fpdR ^= ultrasonicMed->fpdR << 5;

  2888. -		inputSampleR += (((double)ultrasonicMed->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2889. -		//end 32 bit stereo floating point dither

  2890. -

  2891.  		*out1 = (float) inputSampleL;

  2892.  		*out2 = (float) inputSampleR;

  2893. 

  2894. diff --git a/src/optimizations.h b/src/optimizations.h

  2895. new file mode 100644

  2896. index 0000000..8670aea

  2897. --- /dev/null

  2898. +++ b/src/optimizations.h

  2899. @@ -0,0 +1,25 @@

  2900. +#include <math.h>

  2901. +#include <stdint.h>

  2902. +#include <stdlib.h>

  2903. +

  2904. +// force doubles into floats

  2905. +typedef double realdouble;

  2906. +#define double float

  2907. +#define asin asinf

  2908. +#define atan2 atan2f

  2909. +#define cos cosf

  2910. +#define exp expf

  2911. +#define exp2 exp2f

  2912. +#define fabs fabsf

  2913. +#define log2 log2f

  2914. +#define log10 log10f

  2915. +#define pow powf

  2916. +#define sin sinf

  2917. +#define tan tanf

  2918. +

  2919. +// ensure denormal things is unused

  2920. +#define fpdL error-if-fpdL-is-used

  2921. +#define fpdR error-if-fpdR-is-used

  2922. +

  2923. +// ensure 32bit dithering is unused

  2924. +#define frexpf error-if-frexpf-is-used