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.
139 Commits
1 Branch
3.6MB
Tree: 14effa86d8
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '14effa86d8'
${ noResults }
Repository/sources/plugins/airwindows-lv2/debian/patches/01_optimizations.patch

3201 lines
109KB
Raw Blame History 

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

  2. index 37cd72f..5f8e689 100644

  3. --- a/CMakeLists.txt

  4. +++ b/CMakeLists.txt

  5. @@ -95,6 +95,9 @@ add_plugin(UltrasonicLite)

  6.  add_plugin(UltrasonicMed)

  7.  add_plugin(Vibrato)

  8.  

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

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

  11. +

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

  13.  

  14.  if(APPLE OR UNIX)

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

  16. index e502a66..9f0c763 100644

  17. --- a/src/ADClip7/ADClip7.c

  18. +++ b/src/ADClip7/ADClip7.c

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

  20.  #include <math.h>

  21.  #include <stdint.h>

  22.  #include <stdlib.h>

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

  24.  

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

  26.  

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

  28.  	const float* enhance;

  29.  	const float* mode;

  30.  

  31. -	uint32_t fpdL;

  32. -	uint32_t fpdR;

  33.  	//default stuff

  34.  	double lastSampleL;

  35.  	double lastSampleR;

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

  37.  

  38.  static LV2_Handle instantiate(

  39.  	const LV2_Descriptor* descriptor,

  40. -	double rate,

  41. +	realdouble rate,

  42.  	const char* bundle_path,

  43.  	const LV2_Feature* const* features)

  44.  {

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

  46.  	adclip7->iirLowsAR = 0.0;

  47.  	adclip7->iirLowsBL = 0.0;

  48.  	adclip7->iirLowsBR = 0.0;

  49. -

  50. -	adclip7->fpdL = 1.0;

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

  52. -	adclip7->fpdR = 1.0;

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

  54.  }

  55.  

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

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

  58.  	while (sampleFrames-- > 0) {

  59.  		inputSampleL = *in1;

  60.  		inputSampleR = *in2;

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

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

  63. -

  64.  

  65.  		if (inputGain != 1.0) {

  66.  			inputSampleL *= inputGain;

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

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

  69.  		//final iron bar

  70.  

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

  72. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  84. -

  85.  		*out1 = (float) inputSampleL;

  86.  		*out2 = (float) inputSampleR;

  87.  

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

  89. index d6af6c5..f4366a0 100644

  90. --- a/src/Acceleration/Acceleration.c

  91. +++ b/src/Acceleration/Acceleration.c

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

  93.  #include <math.h>

  94.  #include <stdint.h>

  95.  #include <stdlib.h>

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

  97.  

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

  99.  

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

  101.  	const float* limit;

  102.  	const float* drywet;

  103.  

  104. -	uint32_t fpdL;

  105. -	uint32_t fpdR;

  106.  	//default stuff

  107.  	double ataLastOutL;

  108.  	double s1L;

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

  110.  

  111.  static LV2_Handle instantiate(

  112.  	const LV2_Descriptor* descriptor,

  113. -	double rate,

  114. +	realdouble rate,

  115.  	const char* bundle_path,

  116.  	const LV2_Feature* const* features)

  117.  {

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

  119.  	acceleration->m1R = 0.0;

  120.  	acceleration->m2R = 0.0;

  121.  	acceleration->desR = 0.0;

  122. -

  123. -	acceleration->fpdL = 1.0;

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

  125. -	acceleration->fpdR = 1.0;

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

  127.  }

  128.  

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

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

  131.  	while (sampleFrames-- > 0) {

  132.  		inputSampleL = *in1;

  133.  		inputSampleR = *in2;

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

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

  136.  		drySampleL = inputSampleL;

  137.  		drySampleR = inputSampleR;

  138.  

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

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

  141.  		}

  142.  

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

  144. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  156. -

  157.  		*out1 = (float) inputSampleL;

  158.  		*out2 = (float) inputSampleR;

  159.  

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

  161. index 3cdf047..88b5f03 100644

  162. --- a/src/Acceleration2/Acceleration2.c

  163. +++ b/src/Acceleration2/Acceleration2.c

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

  165.  #include <math.h>

  166.  #include <stdint.h>

  167.  #include <stdlib.h>

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

  169.  

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

  171.  

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

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

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

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

  176. -

  177. -	uint32_t fpdL;

  178. -	uint32_t fpdR;

  179.  } Acceleration2;

  180.  

  181.  static LV2_Handle instantiate(

  182.  	const LV2_Descriptor* descriptor,

  183. -	double rate,

  184. +	realdouble rate,

  185.  	const char* bundle_path,

  186.  	const LV2_Feature* const* features)

  187.  {

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

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

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

  191.  	}

  192. -

  193. -	acceleration2->fpdL = 1.0;

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

  195. -	acceleration2->fpdR = 1.0;

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

  197.  }

  198.  

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

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

  201.  	while (sampleFrames-- > 0) {

  202.  		double inputSampleL = *in1;

  203.  		double inputSampleR = *in2;

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

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

  206.  		double drySampleL = inputSampleL;

  207.  		double drySampleR = inputSampleR;

  208.  

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

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

  211.  		}

  212.  

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

  214. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  226. -

  227.  		*out1 = (float) inputSampleL;

  228.  		*out2 = (float) inputSampleR;

  229.  

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

  231. index 263b1ca..82d7f52 100644

  232. --- a/src/Baxandall/Baxandall.c

  233. +++ b/src/Baxandall/Baxandall.c

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

  235.  #include <stdbool.h>

  236.  #include <stdint.h>

  237.  #include <stdlib.h>

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

  239.  

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

  241.  

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

  243.  	const float* bass;

  244.  	const float* outputGain;

  245.  

  246. -	uint32_t fpdL;

  247. -	uint32_t fpdR;

  248.  	//default stuff

  249.  	double trebleAL[9];

  250.  	double trebleBL[9];

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

  252.  

  253.  static LV2_Handle instantiate(

  254.  	const LV2_Descriptor* descriptor,

  255. -	double rate,

  256. +	realdouble rate,

  257.  	const char* bundle_path,

  258.  	const LV2_Feature* const* features)

  259.  {

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

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

  262.  	}

  263.  	baxandall->flip = false;

  264. -	baxandall->fpdL = 1.0;

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

  266. -	baxandall->fpdR = 1.0;

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

  268.  }

  269.  

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

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

  272.  	while (sampleFrames-- > 0) {

  273.  		double inputSampleL = *in1;

  274.  		double inputSampleR = *in2;

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

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

  277.  

  278.  		if (output != 1.0) {

  279.  			inputSampleL *= output;

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

  281.  		inputSampleR = asin(inputSampleR);

  282.  		//amplitude aspect

  283.  

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

  285. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  297. -

  298.  		*out1 = (float) inputSampleL;

  299.  		*out2 = (float) inputSampleR;

  300.  

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

  302. index 4c5a5d1..1235099 100644

  303. --- a/src/BitShiftGain/BitShiftGain.c

  304. +++ b/src/BitShiftGain/BitShiftGain.c

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

  306.  

  307.  #include <stdint.h>

  308.  #include <stdlib.h>

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

  310.  

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

  312.  

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

  314.  

  315.  static LV2_Handle instantiate(

  316.  	const LV2_Descriptor* descriptor,

  317. -	double rate,

  318. +	realdouble rate,

  319.  	const char* bundle_path,

  320.  	const LV2_Feature* const* features)

  321.  {

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

  323. index 907c2d3..fe0d466 100644

  324. --- a/src/Capacitor/Capacitor.c

  325. +++ b/src/Capacitor/Capacitor.c

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

  327.  #include <math.h>

  328.  #include <stdint.h>

  329.  #include <stdlib.h>

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

  331.  

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

  333.  

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

  335.  	double lastLowpass;

  336.  	double lastHighpass;

  337.  	double lastWet;

  338. -

  339. -	uint32_t fpdL;

  340. -	uint32_t fpdR;

  341.  } Capacitor;

  342.  

  343.  static LV2_Handle instantiate(

  344.  	const LV2_Descriptor* descriptor,

  345. -	double rate,

  346. +	realdouble rate,

  347.  	const char* bundle_path,

  348.  	const LV2_Feature* const* features)

  349.  {

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

  351.  	capacitor->lastLowpass = 1000.0;

  352.  	capacitor->lastHighpass = 1000.0;

  353.  	capacitor->lastWet = 1000.0;

  354. -

  355. -	capacitor->fpdL = 1.0;

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

  357. -	capacitor->fpdR = 1.0;

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

  359.  }

  360.  

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

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

  363.  	while (sampleFrames-- > 0) {

  364.  		inputSampleL = *in1;

  365.  		inputSampleR = *in2;

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

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

  368.  		drySampleL = inputSampleL;

  369.  		drySampleR = inputSampleR;

  370.  

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

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

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

  374.  

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

  376. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  388. -

  389.  		*out1 = (float) inputSampleL;

  390.  		*out2 = (float) inputSampleR;

  391.  

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

  393. index 111ff8a..53050eb 100644

  394. --- a/src/Capacitor2/Capacitor2.c

  395. +++ b/src/Capacitor2/Capacitor2.c

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

  397.  #include <math.h>

  398.  #include <stdint.h>

  399.  #include <stdlib.h>

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

  401.  

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

  403.  

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

  405.  	double lastLowpass;

  406.  	double lastHighpass;

  407.  	double lastWet;

  408. -

  409. -	uint32_t fpdL;

  410. -	uint32_t fpdR;

  411.  } Capacitor2;

  412.  

  413.  static LV2_Handle instantiate(

  414.  	const LV2_Descriptor* descriptor,

  415. -	double rate,

  416. +	realdouble rate,

  417.  	const char* bundle_path,

  418.  	const LV2_Feature* const* features)

  419.  {

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

  421.  	capacitor2->lastLowpass = 1000.0;

  422.  	capacitor2->lastHighpass = 1000.0;

  423.  	capacitor2->lastWet = 1000.0;

  424. -

  425. -	capacitor2->fpdL = 1.0;

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

  427. -	capacitor2->fpdR = 1.0;

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

  429.  }

  430.  

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

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

  433.  	while (sampleFrames-- > 0) {

  434.  		double inputSampleL = *in1;

  435.  		double inputSampleR = *in2;

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

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

  438.  		double drySampleL = inputSampleL;

  439.  		double drySampleR = inputSampleR;

  440.  

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

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

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

  444.  

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

  446. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  458. -

  459.  		*out1 = (float) inputSampleL;

  460.  		*out2 = (float) inputSampleR;

  461.  

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

  463. index 5130655..80e75a7 100644

  464. --- a/src/Channel8/Channel8.c

  465. +++ b/src/Channel8/Channel8.c

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

  467.  #include <stdbool.h>

  468.  #include <stdint.h>

  469.  #include <stdlib.h>

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

  471.  

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

  473.  

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

  475.  	const float* drive;

  476.  	const float* outputGain;

  477.  

  478. -	uint32_t fpdL;

  479. -	uint32_t fpdR;

  480.  	//default stuff

  481.  	double iirSampleLA;

  482.  	double iirSampleRA;

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

  484.  

  485.  static LV2_Handle instantiate(

  486.  	const LV2_Descriptor* descriptor,

  487. -	double rate,

  488. +	realdouble rate,

  489.  	const char* bundle_path,

  490.  	const LV2_Feature* const* features)

  491.  {

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

  493.  {

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

  495.  

  496. -	channel8->fpdL = 1.0;

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

  498. -	channel8->fpdR = 1.0;

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

  500.  	channel8->iirSampleLA = 0.0;

  501.  	channel8->iirSampleRA = 0.0;

  502.  	channel8->iirSampleLB = 0.0;

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

  504.  	while (sampleFrames-- > 0) {

  505.  		double inputSampleL = *in1;

  506.  		double inputSampleR = *in2;

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

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

  509.  

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

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

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

  513.  			inputSampleR *= output;

  514.  		}

  515.  

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

  517. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  529. -

  530.  		*out1 = (float) inputSampleL;

  531.  		*out2 = (float) inputSampleR;

  532.  

  533. diff --git a/src/Channel9/Channel9.c b/src/Channel9/Channel9.c

  534. index 5e18790..da51ecd 100644

  535. --- a/src/Channel9/Channel9.c

  536. +++ b/src/Channel9/Channel9.c

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

  538.  #include <stdbool.h>

  539.  #include <stdint.h>

  540.  #include <stdlib.h>

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

  542.  

  543.  #define CHANNEL9_URI "https://hannesbraun.net/ns/lv2/airwindows/channel9"

  544.  

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

  546.  	const float* drive;

  547.  	const float* outputGain;

  548.  

  549. -	uint32_t fpdL;

  550. -	uint32_t fpdR;

  551.  	//default stuff

  552.  	double iirSampleLA;

  553.  	double iirSampleRA;

  554. @@ -48,7 +47,7 @@ typedef struct {

  555.  

  556.  static LV2_Handle instantiate(

  557.  	const LV2_Descriptor* descriptor,

  558. -	double rate,

  559. +	realdouble rate,

  560.  	const char* bundle_path,

  561.  	const LV2_Feature* const* features)

  562.  {

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

  564.  		channel9->biquadA[x] = 0.0;

  565.  		channel9->biquadB[x] = 0.0;

  566.  	}

  567. -	channel9->fpdL = 1.0;

  568. -	while (channel9->fpdL < 16386) channel9->fpdL = rand() * UINT32_MAX;

  569. -	channel9->fpdR = 1.0;

  570. -	while (channel9->fpdR < 16386) channel9->fpdR = rand() * UINT32_MAX;

  571.  	channel9->iirSampleLA = 0.0;

  572.  	channel9->iirSampleRA = 0.0;

  573.  	channel9->iirSampleLB = 0.0;

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

  575.  	while (sampleFrames-- > 0) {

  576.  		double inputSampleL = *in1;

  577.  		double inputSampleR = *in2;

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

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

  580.  		double tempSample;

  581.  

  582.  		if (channel9->biquadA[0] < 0.49999) {

  583. @@ -313,20 +306,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  584.  			channel9->biquadB[13] = inputSampleR; //DF1 right

  585.  		}

  586.  

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

  588. -		int expon;

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

  590. -		channel9->fpdL ^= channel9->fpdL << 13;

  591. -		channel9->fpdL ^= channel9->fpdL >> 17;

  592. -		channel9->fpdL ^= channel9->fpdL << 5;

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

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

  595. -		channel9->fpdR ^= channel9->fpdR << 13;

  596. -		channel9->fpdR ^= channel9->fpdR >> 17;

  597. -		channel9->fpdR ^= channel9->fpdR << 5;

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

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

  600. -

  601.  		*out1 = (float) inputSampleL;

  602.  		*out2 = (float) inputSampleR;

  603.  

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

  605. index 2d83605..94f0ef8 100644

  606. --- a/src/ClipOnly/ClipOnly.c

  607. +++ b/src/ClipOnly/ClipOnly.c

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

  609.  #include <stdbool.h>

  610.  #include <stdint.h>

  611.  #include <stdlib.h>

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

  613.  

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

  615.  

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

  617.  

  618.  static LV2_Handle instantiate(

  619.  	const LV2_Descriptor* descriptor,

  620. -	double rate,

  621. +	realdouble rate,

  622.  	const char* bundle_path,

  623.  	const LV2_Feature* const* features)

  624.  {

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

  626. index eec1237..7f28c8c 100644

  627. --- a/src/ClipOnly2/ClipOnly2.c

  628. +++ b/src/ClipOnly2/ClipOnly2.c

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

  630.  #include <stdbool.h>

  631.  #include <stdint.h>

  632.  #include <stdlib.h>

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

  634.  

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

  636.  

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

  638.  

  639.  static LV2_Handle instantiate(

  640.  	const LV2_Descriptor* descriptor,

  641. -	double rate,

  642. +	realdouble rate,

  643.  	const char* bundle_path,

  644.  	const LV2_Feature* const* features)

  645.  {

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

  647. index f6f006e..2cd7e8b 100644

  648. --- a/src/Console7Buss/Console7Buss.c

  649. +++ b/src/Console7Buss/Console7Buss.c

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

  651.  #include <math.h>

  652.  #include <stdint.h>

  653.  #include <stdlib.h>

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

  655.  

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

  657.  

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

  659.  

  660.  	double biquadA[15];

  661.  	double biquadB[15];

  662. -

  663. -	uint32_t fpdL;

  664. -	uint32_t fpdR;

  665.  } Console7Buss;

  666.  

  667.  static LV2_Handle instantiate(

  668.  	const LV2_Descriptor* descriptor,

  669. -	double rate,

  670. +	realdouble rate,

  671.  	const char* bundle_path,

  672.  	const LV2_Feature* const* features)

  673.  {

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

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

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

  677.  	}

  678. -	console7Buss->fpdL = 1.0;

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

  680. -	console7Buss->fpdR = 1.0;

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

  682.  }

  683.  

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

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

  686.  	while (sampleFrames-- > 0) {

  687.  		double inputSampleL = *in1;

  688.  		double inputSampleR = *in2;

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

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

  691.  

  692.  		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];

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

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

  695.  		}

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

  697.  

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

  699. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  711. -

  712.  		*out1 = (float) inputSampleL;

  713.  		*out2 = (float) inputSampleR;

  714.  

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

  716. index 2f5f8ee..e4fd69d 100644

  717. --- a/src/Console7Cascade/Console7Cascade.c

  718. +++ b/src/Console7Cascade/Console7Cascade.c

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

  720.  #include <math.h>

  721.  #include <stdint.h>

  722.  #include <stdlib.h>

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

  724.  

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

  726.  

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

  728.  	double biquadC[15];

  729.  	double biquadD[15];

  730.  	double biquadE[15];

  731. -

  732. -	uint32_t fpdL;

  733. -	uint32_t fpdR;

  734.  } Console7Cascade;

  735.  

  736.  static LV2_Handle instantiate(

  737.  	const LV2_Descriptor* descriptor,

  738. -	double rate,

  739. +	realdouble rate,

  740.  	const char* bundle_path,

  741.  	const LV2_Feature* const* features)

  742.  {

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

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

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

  746.  	}

  747. -	console7Cascade->fpdL = 1.0;

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

  749. -	console7Cascade->fpdR = 1.0;

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

  751.  }

  752.  

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

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

  755.  	while (sampleFrames-- > 0) {

  756.  		double inputSampleL = *in1;

  757.  		double inputSampleR = *in2;

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

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

  760.  

  761.  		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];

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

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

  764.  		}

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

  766.  

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

  768. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  780. -

  781.  		*out1 = (float) inputSampleL;

  782.  		*out2 = (float) inputSampleR;

  783.  

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

  785. index 9d531c4..e76936a 100644

  786. --- a/src/Console7Channel/Console7Channel.c

  787. +++ b/src/Console7Channel/Console7Channel.c

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

  789.  #include <math.h>

  790.  #include <stdint.h>

  791.  #include <stdlib.h>

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

  793.  

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

  795.  

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

  797.  	double chasespeed;

  798.  

  799.  	double biquadA[15];

  800. -

  801. -	uint32_t fpdL;

  802. -	uint32_t fpdR;

  803.  } Console7Channel;

  804.  

  805.  static LV2_Handle instantiate(

  806.  	const LV2_Descriptor* descriptor,

  807. -	double rate,

  808. +	realdouble rate,

  809.  	const char* bundle_path,

  810.  	const LV2_Feature* const* features)

  811.  {

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

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

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

  815.  	}

  816. -	console7Channel->fpdL = 1.0;

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

  818. -	console7Channel->fpdR = 1.0;

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

  820.  }

  821.  

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

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

  824.  	while (sampleFrames-- > 0) {

  825.  		double inputSampleL = *in1;

  826.  		double inputSampleR = *in2;

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

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

  829.  

  830.  		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];

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

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

  833.  		}

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

  835.  

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

  837. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  849. -

  850.  		*out1 = (float) inputSampleL;

  851.  		*out2 = (float) inputSampleR;

  852.  

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

  854. index 8f3049e..7d51bb7 100644

  855. --- a/src/Console7Crunch/Console7Crunch.c

  856. +++ b/src/Console7Crunch/Console7Crunch.c

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

  858.  #include <math.h>

  859.  #include <stdint.h>

  860.  #include <stdlib.h>

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

  862.  

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

  864.  

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

  866.  	double biquadC[15];

  867.  	double biquadD[15];

  868.  	double biquadE[15];

  869. -

  870. -	uint32_t fpdL;

  871. -	uint32_t fpdR;

  872.  } Console7Crunch;

  873.  

  874.  static LV2_Handle instantiate(

  875.  	const LV2_Descriptor* descriptor,

  876. -	double rate,

  877. +	realdouble rate,

  878.  	const char* bundle_path,

  879.  	const LV2_Feature* const* features)

  880.  {

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

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

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

  884.  	}

  885. -	console7Crunch->fpdL = 1.0;

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

  887. -	console7Crunch->fpdR = 1.0;

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

  889.  }

  890.  

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

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

  893.  	while (sampleFrames-- > 0) {

  894.  		double inputSampleL = *in1;

  895.  		double inputSampleR = *in2;

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

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

  898.  

  899.  		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];

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

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

  902.  		}

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

  904.  

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

  906. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  918. -

  919.  		*out1 = (float) inputSampleL;

  920.  		*out2 = (float) inputSampleR;

  921.  

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

  923. index 18f904e..8d17086 100644

  924. --- a/src/DCVoltage/DCVoltage.c

  925. +++ b/src/DCVoltage/DCVoltage.c

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

  927.  

  928.  #include <stdint.h>

  929.  #include <stdlib.h>

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

  931.  

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

  933.  

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

  935.  

  936.  static LV2_Handle instantiate(

  937.  	const LV2_Descriptor* descriptor,

  938. -	double rate,

  939. +	realdouble rate,

  940.  	const char* bundle_path,

  941.  	const LV2_Feature* const* features)

  942.  {

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

  944. index 54a84d6..85ec475 100644

  945. --- a/src/DeBess/DeBess.c

  946. +++ b/src/DeBess/DeBess.c

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

  948.  #include <stdbool.h>

  949.  #include <stdint.h>

  950.  #include <stdlib.h>

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

  952.  

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

  954.  

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

  956.  

  957.  

  958.  	bool flip;

  959. -	uint32_t fpdL;

  960. -	uint32_t fpdR;

  961.  } DeBess;

  962.  

  963.  static LV2_Handle instantiate(

  964.  	const LV2_Descriptor* descriptor,

  965. -	double rate,

  966. +	realdouble rate,

  967.  	const char* bundle_path,

  968.  	const LV2_Feature* const* features)

  969.  {

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

  971.  	debess->iirSampleBR = 0.0;

  972.  

  973.  	debess->flip = false;

  974. -	debess->fpdL = 1.0;

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

  976. -	debess->fpdR = 1.0;

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

  978.  }

  979.  

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

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

  982.  	while (sampleFrames-- > 0) {

  983.  		double inputSampleL = *in1;

  984.  		double inputSampleR = *in2;

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

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

  987.  

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

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

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

  991.  		}

  992.  		//sense monitoring

  993.  

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

  995. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1007. -

  1008.  		*out1 = (float) inputSampleL;

  1009.  		*out2 = (float) inputSampleR;

  1010.  

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

  1012. index ed028dd..cbbdf58 100644

  1013. --- a/src/Dyno/Dyno.c

  1014. +++ b/src/Dyno/Dyno.c

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

  1016.  #include <math.h>

  1017.  #include <stdint.h>

  1018.  #include <stdlib.h>

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

  1020.  

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

  1022.  

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

  1024.  	const float* input[2];

  1025.  	float* output[2];

  1026.  	const float* gain;

  1027. -

  1028. -	uint32_t fpdL;

  1029. -	uint32_t fpdR;

  1030.  } Dyno;

  1031.  

  1032.  static LV2_Handle instantiate(

  1033.  	const LV2_Descriptor* descriptor,

  1034. -	double rate,

  1035. +	realdouble rate,

  1036.  	const char* bundle_path,

  1037.  	const LV2_Feature* const* features)

  1038.  {

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

  1040.  	}

  1041.  }

  1042.  

  1043. -static void activate(LV2_Handle instance)

  1044. -{

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

  1046. -	dyno->fpdL = 1.0;

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

  1048. -	dyno->fpdR = 1.0;

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

  1050. -}

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

  1052.  

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

  1054.  {

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

  1056.  	while (sampleFrames-- > 0) {

  1057.  		double inputSampleL = *in1;

  1058.  		double inputSampleR = *in2;

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

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

  1061.  

  1062.  		if (gain != 1.0) {

  1063.  			inputSampleL *= gain;

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

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

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

  1067.  

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

  1069. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1081. -

  1082.  		*out1 = (float) inputSampleL;

  1083.  		*out2 = (float) inputSampleR;

  1084.  

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

  1086. index 9fb5f2d..3a2d982 100644

  1087. --- a/src/EdIsDim/EdIsDim.c

  1088. +++ b/src/EdIsDim/EdIsDim.c

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

  1090.  #include <math.h>

  1091.  #include <stdint.h>

  1092.  #include <stdlib.h>

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

  1094.  

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

  1096.  

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

  1098.  	const float* input[2];

  1099.  	float* output[2];

  1100.  	const float* edisdim;

  1101. -

  1102. -	uint32_t fpdL;

  1103. -	uint32_t fpdR;

  1104.  } EdIsDim;

  1105.  

  1106.  static LV2_Handle instantiate(

  1107.  	const LV2_Descriptor* descriptor,

  1108. -	double rate,

  1109. +	realdouble rate,

  1110.  	const char* bundle_path,

  1111.  	const LV2_Feature* const* features)

  1112.  {

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

  1114.  	}

  1115.  }

  1116.  

  1117. -static void activate(LV2_Handle instance)

  1118. -{

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

  1120. -	edisdim->fpdL = 1.0;

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

  1122. -	edisdim->fpdR = 1.0;

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

  1124. -}

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

  1126.  

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

  1128.  {

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

  1130.  	while (sampleFrames-- > 0) {

  1131.  		inputSampleL = *in1;

  1132.  		inputSampleR = *in2;

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

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

  1135.  

  1136.  		inputSampleL *= midgain;

  1137.  		inputSampleR *= sidegain;

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

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

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

  1141.  

  1142. -		//begin 32 bit floating point dither

  1143. -		int expon;

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

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

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

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

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

  1149. -		//end 32 bit floating point dither

  1150. -		//begin 32 bit floating point dither

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

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

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

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

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

  1156. -		//end 32 bit floating point dither

  1157. -

  1158.  		*out1 = (float) mid;

  1159.  		*out2 = (float) side;

  1160.  

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

  1162. index d22e1b9..58f8c57 100644

  1163. --- a/src/EveryTrim/EveryTrim.c

  1164. +++ b/src/EveryTrim/EveryTrim.c

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

  1166.  #include <math.h>

  1167.  #include <stdint.h>

  1168.  #include <stdlib.h>

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

  1170.  

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

  1172.  

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

  1174.  	const float* mid;

  1175.  	const float* side;

  1176.  	const float* master;

  1177. -

  1178. -	uint32_t fpdL;

  1179. -	uint32_t fpdR;

  1180.  } EveryTrim;

  1181.  

  1182.  static LV2_Handle instantiate(

  1183.  	const LV2_Descriptor* descriptor,

  1184. -	double rate,

  1185. +	realdouble rate,

  1186.  	const char* bundle_path,

  1187.  	const LV2_Feature* const* features)

  1188.  {

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

  1190.  	}

  1191.  }

  1192.  

  1193. -static void activate(LV2_Handle instance)

  1194. -{

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

  1196. -	everytrim->fpdL = 1.0;

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

  1198. -	everytrim->fpdR = 1.0;

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

  1200. -}

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

  1202.  

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

  1204.  {

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

  1206.  	while (sampleFrames-- > 0) {

  1207.  		inputSampleL = *in1;

  1208.  		inputSampleR = *in2;

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

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

  1211.  

  1212.  		mid = inputSampleL + inputSampleR;

  1213.  		side = inputSampleL - inputSampleR;

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

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

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

  1217.  

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

  1219. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1231. -

  1232.  		*out1 = (float) inputSampleL;

  1233.  		*out2 = (float) inputSampleR;

  1234.  

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

  1236. index 4629372..e38e881 100644

  1237. --- a/src/Galactic/Galactic.c

  1238. +++ b/src/Galactic/Galactic.c

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

  1240.  #include <math.h>

  1241.  #include <stdint.h>

  1242.  #include <stdlib.h>

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

  1244.  

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

  1246.  

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

  1248.  

  1249.  	double aML[3111];

  1250.  	double aMR[3111];

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

  1252. +	double vibML, vibMR, depthM;

  1253.  

  1254.  	double feedbackAL;

  1255.  	double feedbackBL;

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

  1257.  

  1258.  	double vibM;

  1259.  

  1260. -	uint32_t fpdL;

  1261. -	uint32_t fpdR;

  1262.  	//default stuff

  1263.  } Galactic;

  1264.  

  1265.  static LV2_Handle instantiate(

  1266.  	const LV2_Descriptor* descriptor,

  1267. -	double rate,

  1268. +	realdouble rate,

  1269.  	const char* bundle_path,

  1270.  	const LV2_Feature* const* features)

  1271.  {

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

  1273.  	galactic->cycle = 0;

  1274.  

  1275.  	galactic->vibM = 3.0;

  1276. -

  1277. -	galactic->oldfpd = 429496.7295;

  1278. -

  1279. -	galactic->fpdL = 1;

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

  1281. -	galactic->fpdR = 1;

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

  1283.  }

  1284.  

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

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

  1287.  	while (sampleFrames-- > 0) {

  1288.  		double inputSampleL = *in1;

  1289.  		double inputSampleR = *in2;

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

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

  1292.  		double drySampleL = inputSampleL;

  1293.  		double drySampleR = inputSampleR;

  1294.  

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

  1296. +		galactic->vibM += drift;

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

  1298.  			galactic->vibM = 0.0;

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

  1300.  		}

  1301.  

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

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

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

  1305.  		}

  1306.  

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

  1308. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1320. -

  1321.  		*out1 = (float) inputSampleL;

  1322.  		*out2 = (float) inputSampleR;

  1323.  

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

  1325. index 1f2d728..fd754f4 100644

  1326. --- a/src/HermeTrim/HermeTrim.c

  1327. +++ b/src/HermeTrim/HermeTrim.c

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

  1329.  #include <math.h>

  1330.  #include <stdint.h>

  1331.  #include <stdlib.h>

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

  1333.  

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

  1335.  

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

  1337.  	const float* mid;

  1338.  	const float* side;

  1339.  	const float* master;

  1340. -

  1341. -	uint32_t fpdL;

  1342. -	uint32_t fpdR;

  1343.  } HermeTrim;

  1344.  

  1345.  static LV2_Handle instantiate(

  1346.  	const LV2_Descriptor* descriptor,

  1347. -	double rate,

  1348. +	realdouble rate,

  1349.  	const char* bundle_path,

  1350.  	const LV2_Feature* const* features)

  1351.  {

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

  1353.  	}

  1354.  }

  1355.  

  1356. -static void activate(LV2_Handle instance)

  1357. -{

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

  1359. -	hermetrim->fpdL = 1.0;

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

  1361. -	hermetrim->fpdR = 1.0;

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

  1363. -}

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

  1365.  

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

  1367.  {

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

  1369.  	while (sampleFrames-- > 0) {

  1370.  		inputSampleL = *in1;

  1371.  		inputSampleR = *in2;

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

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

  1374.  

  1375.  		mid = inputSampleL + inputSampleR;

  1376.  		side = inputSampleL - inputSampleR;

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

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

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

  1380.  

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

  1382. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1394. -

  1395.  		*out1 = (float) inputSampleL;

  1396.  		*out2 = (float) inputSampleR;

  1397.  

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

  1399. index 68e8a7e..4f9daa7 100644

  1400. --- a/src/Hypersonic/Hypersonic.c

  1401. +++ b/src/Hypersonic/Hypersonic.c

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

  1403.  #include <math.h>

  1404.  #include <stdint.h>

  1405.  #include <stdlib.h>

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

  1407.  

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

  1409.  

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

  1411.  	double fixE[fix_total];

  1412.  	double fixF[fix_total];

  1413.  	double fixG[fix_total];

  1414. -

  1415. -	uint32_t fpdL;

  1416. -	uint32_t fpdR;

  1417.  	//default stuff

  1418.  } Hypersonic;

  1419.  

  1420.  static LV2_Handle instantiate(

  1421.  	const LV2_Descriptor* descriptor,

  1422. -	double rate,

  1423. +	realdouble rate,

  1424.  	const char* bundle_path,

  1425.  	const LV2_Feature* const* features)

  1426.  {

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

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

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

  1430.  	}

  1431. -	hypersonic->fpdL = 1.0;

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

  1433. -	hypersonic->fpdR = 1.0;

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

  1435.  }

  1436.  

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

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

  1439.  	while (sampleFrames-- > 0) {

  1440.  		double inputSampleL = *in1;

  1441.  		double inputSampleR = *in2;

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

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

  1444.  

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

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

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

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

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

  1450.  

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

  1452. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1464. -

  1465.  		*out1 = (float) inputSampleL;

  1466.  		*out2 = (float) inputSampleR;

  1467.  

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

  1469. index 82c2d71..7288316 100644

  1470. --- a/src/Infrasonic/Infrasonic.c

  1471. +++ b/src/Infrasonic/Infrasonic.c

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

  1473.  #include <math.h>

  1474.  #include <stdint.h>

  1475.  #include <stdlib.h>

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

  1477.  

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

  1479.  

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

  1481.  	double biquadC[15];

  1482.  	double biquadD[15];

  1483.  	double biquadE[15];

  1484. -	uint32_t fpdL;

  1485. -	uint32_t fpdR;

  1486.  } Infrasonic;

  1487.  

  1488.  static LV2_Handle instantiate(

  1489.  	const LV2_Descriptor* descriptor,

  1490. -	double rate,

  1491. +	realdouble rate,

  1492.  	const char* bundle_path,

  1493.  	const LV2_Feature* const* features)

  1494.  {

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

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

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

  1498.  	}

  1499. -	infrasonic->fpdL = 1.0;

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

  1501. -	infrasonic->fpdR = 1.0;

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

  1503.  }

  1504.  

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

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

  1507.  	while (sampleFrames-- > 0) {

  1508.  		double inputSampleL = *in1;

  1509.  		double inputSampleR = *in2;

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

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

  1512.  

  1513.  		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];

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

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

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

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

  1518.  

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

  1520. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1532. -

  1533.  		*out1 = (float) inputSampleL;

  1534.  		*out2 = (float) inputSampleR;

  1535.  

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

  1537. index 885ef83..2d4df77 100644

  1538. --- a/src/Interstage/Interstage.c

  1539. +++ b/src/Interstage/Interstage.c

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

  1541.  #include <stdbool.h>

  1542.  #include <stdint.h>

  1543.  #include <stdlib.h>

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

  1545.  

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

  1547.  

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

  1549.  	double iirSampleER;

  1550.  	double iirSampleFR;

  1551.  	double lastSampleR;

  1552. -	uint32_t fpdL;

  1553. -	uint32_t fpdR;

  1554.  	bool flip;

  1555.  } Interstage;

  1556.  

  1557.  static LV2_Handle instantiate(

  1558.  	const LV2_Descriptor* descriptor,

  1559. -	double rate,

  1560. +	realdouble rate,

  1561.  	const char* bundle_path,

  1562.  	const LV2_Feature* const* features)

  1563.  {

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

  1565.  	interstage->iirSampleER = 0.0;

  1566.  	interstage->iirSampleFR = 0.0;

  1567.  	interstage->lastSampleR = 0.0;

  1568. -	interstage->fpdL = 1.0;

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

  1570. -	interstage->fpdR = 1.0;

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

  1572.  	interstage->flip = true;

  1573.  }

  1574.  

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

  1576.  	while (sampleFrames-- > 0) {

  1577.  		double inputSampleL = *in1;

  1578.  		double inputSampleR = *in2;

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

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

  1581.  		double drySampleL = inputSampleL;

  1582.  		double drySampleR = inputSampleR;

  1583.  

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

  1585.  		interstage->lastSampleL = inputSampleL;

  1586.  		interstage->lastSampleR = inputSampleR;

  1587.  

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

  1589. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1601. -

  1602.  		*out1 = (float) inputSampleL;

  1603.  		*out2 = (float) inputSampleR;

  1604.  

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

  1606. index bd0b49d..8fedca5 100644

  1607. --- a/src/LeftoMono/LeftoMono.c

  1608. +++ b/src/LeftoMono/LeftoMono.c

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

  1610.  

  1611.  #include <stdint.h>

  1612.  #include <stdlib.h>

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

  1614.  

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

  1616.  

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

  1618.  

  1619.  static LV2_Handle instantiate(

  1620.  	const LV2_Descriptor* descriptor,

  1621. -	double rate,

  1622. +	realdouble rate,

  1623.  	const char* bundle_path,

  1624.  	const LV2_Feature* const* features)

  1625.  {

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

  1627. index 8d57a06..25f7b82 100644

  1628. --- a/src/MV/MV.c

  1629. +++ b/src/MV/MV.c

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

  1631.  #include <math.h>

  1632.  #include <stdint.h>

  1633.  #include <stdlib.h>

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

  1635.  

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

  1637.  

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

  1639.  	const float* gain;

  1640.  	const float* drywet;

  1641.  

  1642. -	uint32_t fpdL;

  1643. -	uint32_t fpdR;

  1644. -

  1645.  	double aAL[15150];

  1646.  	double aBL[14618];

  1647.  	double aCL[14358];

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

  1649.  

  1650.  static LV2_Handle instantiate(

  1651.  	const LV2_Descriptor* descriptor,

  1652. -	double rate,

  1653. +	realdouble rate,

  1654.  	const char* bundle_path,

  1655.  	const LV2_Feature* const* features)

  1656.  {

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

  1658.  

  1659.  	mv->feedbackL = 0.0;

  1660.  	mv->feedbackR = 0.0;

  1661. -

  1662. -	mv->fpdL = 1.0;

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

  1664. -	mv->fpdR = 1.0;

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

  1666.  }

  1667.  

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

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

  1670.  		double inputSampleL = *in1;

  1671.  		double inputSampleR = *in2;

  1672.  

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

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

  1675.  		double drySampleL = inputSampleL;

  1676.  		double drySampleR = inputSampleR;

  1677.  

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

  1679.  		}

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

  1681.  

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

  1683. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1695. -

  1696.  		*out1 = (float) inputSampleL;

  1697.  		*out2 = (float) inputSampleR;

  1698.  

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

  1700. index 5250dd8..933d3a7 100644

  1701. --- a/src/Mackity/Mackity.c

  1702. +++ b/src/Mackity/Mackity.c

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

  1704.  #include <math.h>

  1705.  #include <stdint.h>

  1706.  #include <stdlib.h>

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

  1708.  

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

  1710.  

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

  1712.  	double iirSampleBR;

  1713.  	double biquadA[15];

  1714.  	double biquadB[15];

  1715. -

  1716. -	uint32_t fpdL;

  1717. -	uint32_t fpdR;

  1718.  } Mackity;

  1719.  

  1720.  static LV2_Handle instantiate(

  1721.  	const LV2_Descriptor* descriptor,

  1722. -	double rate,

  1723. +	realdouble rate,

  1724.  	const char* bundle_path,

  1725.  	const LV2_Feature* const* features)

  1726.  {

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

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

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

  1730.  	}

  1731. -

  1732. -	mackity->fpdL = 1.0;

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

  1734. -	mackity->fpdR = 1.0;

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

  1736.  }

  1737.  

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

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

  1740.  	while (sampleFrames-- > 0) {

  1741.  		double inputSampleL = *in1;

  1742.  		double inputSampleR = *in2;

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

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

  1745.  

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

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

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

  1749.  			inputSampleR *= outPad;

  1750.  		}

  1751.  

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

  1753. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1765. -

  1766.  		*out1 = (float) inputSampleL;

  1767.  		*out2 = (float) inputSampleR;

  1768.  

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

  1770. index 4c27ce7..7ee0c90 100644

  1771. --- a/src/MidSide/MidSide.c

  1772. +++ b/src/MidSide/MidSide.c

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

  1774.  #include <math.h>

  1775.  #include <stdint.h>

  1776.  #include <stdlib.h>

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

  1778.  

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

  1780.  

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

  1782.  	const float* input[2];

  1783.  	float* output[2];

  1784.  	const float* midside;

  1785. -

  1786. -	uint32_t fpdL;

  1787. -	uint32_t fpdR;

  1788.  } MidSide;

  1789.  

  1790.  static LV2_Handle instantiate(

  1791.  	const LV2_Descriptor* descriptor,

  1792. -	double rate,

  1793. +	realdouble rate,

  1794.  	const char* bundle_path,

  1795.  	const LV2_Feature* const* features)

  1796.  {

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

  1798.  	}

  1799.  }

  1800.  

  1801. -static void activate(LV2_Handle instance)

  1802. -{

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

  1804. -	midside->fpdL = 1.0;

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

  1806. -	midside->fpdR = 1.0;

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

  1808. -}

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

  1810.  

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

  1812.  {

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

  1814.  	while (sampleFrames-- > 0) {

  1815.  		inputSampleL = *in1;

  1816.  		inputSampleR = *in2;

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

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

  1819.  

  1820.  		mid = inputSampleL + inputSampleR;

  1821.  		side = inputSampleL - inputSampleR;

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

  1823.  		mid *= midgain;

  1824.  		side *= sidegain;

  1825.  

  1826. -		//begin 32 bit floating point dither

  1827. -		int expon;

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

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

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

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

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

  1833. -		//end 32 bit floating point dither

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

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

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

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

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

  1839. -		//end 32 bit floating point dither

  1840. -

  1841.  		*out1 = (float) mid;

  1842.  		*out2 = (float) side;

  1843.  

  1844. diff --git a/src/MoNoam/MoNoam.c b/src/MoNoam/MoNoam.c

  1845. index 8c9a33b..b59a4bc 100644

  1846. --- a/src/MoNoam/MoNoam.c

  1847. +++ b/src/MoNoam/MoNoam.c

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

  1849.  #include <math.h>

  1850.  #include <stdint.h>

  1851.  #include <stdlib.h>

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

  1853.  

  1854.  #define MONOAM_URI "https://hannesbraun.net/ns/lv2/airwindows/monoam"

  1855.  

  1856. @@ -27,14 +28,11 @@ typedef struct {

  1857.  	const float* input[2];

  1858.  	float* output[2];

  1859.  	const float* monotor;

  1860. -

  1861. -	uint32_t fpdL;

  1862. -	uint32_t fpdR;

  1863.  } MoNoam;

  1864.  

  1865.  static LV2_Handle instantiate(

  1866.  	const LV2_Descriptor* descriptor,

  1867. -	double rate,

  1868. +	realdouble rate,

  1869.  	const char* bundle_path,

  1870.  	const LV2_Feature* const* features)

  1871.  {

  1872. @@ -65,14 +63,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  1873.  	}

  1874.  }

  1875.  

  1876. -static void activate(LV2_Handle instance)

  1877. -{

  1878. -	MoNoam* moNoam = (MoNoam*) instance;

  1879. -	moNoam->fpdL = 1.0;

  1880. -	while (moNoam->fpdL < 16386) moNoam->fpdL = rand() * UINT32_MAX;

  1881. -	moNoam->fpdR = 1.0;

  1882. -	while (moNoam->fpdR < 16386) moNoam->fpdR = rand() * UINT32_MAX;

  1883. -}

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

  1885.  

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

  1887.  {

  1888. @@ -93,8 +84,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1889.  	while (sampleFrames-- > 0) {

  1890.  		double inputSampleL = *in1;

  1891.  		double inputSampleR = *in2;

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

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

  1894.  

  1895.  		double mid;

  1896.  		mid = inputSampleL + inputSampleR;

  1897. @@ -117,20 +106,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  1898.  			inputSampleR = *in2;

  1899.  		}

  1900.  

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

  1902. -		int expon;

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

  1904. -		moNoam->fpdL ^= moNoam->fpdL << 13;

  1905. -		moNoam->fpdL ^= moNoam->fpdL >> 17;

  1906. -		moNoam->fpdL ^= moNoam->fpdL << 5;

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

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

  1909. -		moNoam->fpdR ^= moNoam->fpdR << 13;

  1910. -		moNoam->fpdR ^= moNoam->fpdR >> 17;

  1911. -		moNoam->fpdR ^= moNoam->fpdR << 5;

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

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

  1914. -

  1915.  		*out1 = (float) inputSampleL;

  1916.  		*out2 = (float) inputSampleR;

  1917.  

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

  1919. index ee39090..d8ef7cf 100644

  1920. --- a/src/Mojo/Mojo.c

  1921. +++ b/src/Mojo/Mojo.c

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

  1923.  #include <math.h>

  1924.  #include <stdint.h>

  1925.  #include <stdlib.h>

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

  1927.  

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

  1929.  

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

  1931.  	const float* input[2];

  1932.  	float* output[2];

  1933.  	const float* gain;

  1934. -

  1935. -	uint32_t fpdL;

  1936. -	uint32_t fpdR;

  1937.  } Mojo;

  1938.  

  1939.  static LV2_Handle instantiate(

  1940.  	const LV2_Descriptor* descriptor,

  1941. -	double rate,

  1942. +	realdouble rate,

  1943.  	const char* bundle_path,

  1944.  	const LV2_Feature* const* features)

  1945.  {

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

  1947.  	}

  1948.  }

  1949.  

  1950. -static void activate(LV2_Handle instance)

  1951. -{

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

  1953. -	mojo->fpdL = 1.0;

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

  1955. -	mojo->fpdR = 1.0;

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

  1957. -}

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

  1959.  

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

  1961.  {

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

  1963.  	while (sampleFrames-- > 0) {

  1964.  		double inputSampleL = *in1;

  1965.  		double inputSampleR = *in2;

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

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

  1968.  

  1969.  		if (gain != 1.0) {

  1970.  			inputSampleL *= gain;

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

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

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

  1974.  

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

  1976. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  1988. -

  1989.  		*out1 = (float) inputSampleL;

  1990.  		*out2 = (float) inputSampleR;

  1991.  

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

  1993. index 772fbe2..0a5a274 100644

  1994. --- a/src/Nikola/Nikola.c

  1995. +++ b/src/Nikola/Nikola.c

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

  1997.  #include <stdbool.h>

  1998.  #include <stdint.h>

  1999.  #include <stdlib.h>

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

  2001.  

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

  2003.  

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

  2005.  	int framenumberR;

  2006.  	bool wasNegativeR;

  2007.  	double outlevelR;

  2008. -

  2009. -	uint32_t fpdL;

  2010. -	uint32_t fpdR;

  2011.  } Nikola;

  2012.  

  2013.  static LV2_Handle instantiate(

  2014.  	const LV2_Descriptor* descriptor,

  2015. -	double rate,

  2016. +	realdouble rate,

  2017.  	const char* bundle_path,

  2018.  	const LV2_Feature* const* features)

  2019.  {

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

  2021.  	nikola->wasNegativeR = false;

  2022.  	nikola->outlevelR = 0.0;

  2023.  	nikola->framenumberR = 0;

  2024. -

  2025. -	nikola->fpdL = 1.0;

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

  2027. -	nikola->fpdR = 1.0;

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

  2029.  }

  2030.  

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

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

  2033.  	while (sampleFrames-- > 0) {

  2034.  		double inputSampleL = *in1;

  2035.  		double inputSampleR = *in2;

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

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

  2038.  		double drySampleL = inputSampleL;

  2039.  		double drySampleR = inputSampleR;

  2040.  

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

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

  2043.  		}

  2044.  

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

  2046. -		int expon;

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

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

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

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

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

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

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

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

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

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

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

  2058. -

  2059.  		*out1 = (float) inputSampleL;

  2060.  		*out2 = (float) inputSampleR;

  2061.  

  2062. diff --git a/src/PocketVerbs/PocketVerbs.c b/src/PocketVerbs/PocketVerbs.c

  2063. index 694321a..df4be1f 100644

  2064. --- a/src/PocketVerbs/PocketVerbs.c

  2065. +++ b/src/PocketVerbs/PocketVerbs.c

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

  2067.  #include <math.h>

  2068.  #include <stdint.h>

  2069.  #include <stdlib.h>

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

  2071.  

  2072.  #define POCKETVERBS_URI "https://hannesbraun.net/ns/lv2/airwindows/pocketverbs"

  2073.  

  2074. @@ -245,14 +246,11 @@ typedef struct {

  2075.  	int countdown;

  2076.  	double peakL;

  2077.  	double peakR;

  2078. -

  2079. -	uint32_t fpdL;

  2080. -	uint32_t fpdR;

  2081.  } PocketVerbs;

  2082.  

  2083.  static LV2_Handle instantiate(

  2084.  	const LV2_Descriptor* descriptor,

  2085. -	double rate,

  2086. +	realdouble rate,

  2087.  	const char* bundle_path,

  2088.  	const LV2_Feature* const* features)

  2089.  {

  2090. @@ -831,10 +829,6 @@ static void activate(LV2_Handle instance)

  2091.  	pocketverbs->countdown = -1;

  2092.  	pocketverbs->peakL = 1.0;

  2093.  	pocketverbs->peakR = 1.0;

  2094. -	pocketverbs->fpdL = 1.0;

  2095. -	while (pocketverbs->fpdL < 16386) pocketverbs->fpdL = rand() * UINT32_MAX;

  2096. -	pocketverbs->fpdR = 1.0;

  2097. -	while (pocketverbs->fpdR < 16386) pocketverbs->fpdR = rand() * UINT32_MAX;

  2098.  }

  2099.  

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

  2101. @@ -1778,9 +1772,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2102.  		double inputSampleL = *in1;

  2103.  		double inputSampleR = *in2;

  2104.  

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

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

  2107. -

  2108.  		double drySampleL = inputSampleL;

  2109.  		double drySampleR = inputSampleR;

  2110.  

  2111. @@ -13590,20 +13581,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2112.  		inputSampleR += drySampleR;

  2113.  		//here we combine the tanks with the dry signal

  2114.  

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

  2116. -		int expon;

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

  2118. -		pocketverbs->fpdL ^= pocketverbs->fpdL << 13;

  2119. -		pocketverbs->fpdL ^= pocketverbs->fpdL >> 17;

  2120. -		pocketverbs->fpdL ^= pocketverbs->fpdL << 5;

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

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

  2123. -		pocketverbs->fpdR ^= pocketverbs->fpdR << 13;

  2124. -		pocketverbs->fpdR ^= pocketverbs->fpdR >> 17;

  2125. -		pocketverbs->fpdR ^= pocketverbs->fpdR << 5;

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

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

  2128. -

  2129.  		*out1 = (float) inputSampleL;

  2130.  		*out2 = (float) inputSampleR;

  2131.  

  2132. diff --git a/src/Pressure5/Pressure5.c b/src/Pressure5/Pressure5.c

  2133. index 0ea5c79..e16c291 100644

  2134. --- a/src/Pressure5/Pressure5.c

  2135. +++ b/src/Pressure5/Pressure5.c

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

  2137.  #include <math.h>

  2138.  #include <stdint.h>

  2139.  #include <stdlib.h>

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

  2141.  

  2142.  #define PRESSURE5_URI "https://hannesbraun.net/ns/lv2/airwindows/pressure5"

  2143.  

  2144. @@ -70,14 +71,11 @@ typedef struct {

  2145.  	bool wasNegClipR; //Stereo ClipOnly2

  2146.  

  2147.  	double slewMax; //to adust mewiness

  2148. -

  2149. -	uint32_t fpdL;

  2150. -	uint32_t fpdR;

  2151.  } Pressure5;

  2152.  

  2153.  static LV2_Handle instantiate(

  2154.  	const LV2_Descriptor* descriptor,

  2155. -	double rate,

  2156. +	realdouble rate,

  2157.  	const char* bundle_path,

  2158.  	const LV2_Feature* const* features)

  2159.  {

  2160. @@ -149,11 +147,6 @@ static void activate(LV2_Handle instance)

  2161.  		pressure5->intermediateR[x] = 0.0;

  2162.  	}

  2163.  	pressure5->slewMax = 0.0;

  2164. -

  2165. -	pressure5->fpdL = 1.0;

  2166. -	while (pressure5->fpdL < 16386) pressure5->fpdL = rand() * UINT32_MAX;

  2167. -	pressure5->fpdR = 1.0;

  2168. -	while (pressure5->fpdR < 16386) pressure5->fpdR = rand() * UINT32_MAX;

  2169.  }

  2170.  

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

  2172. @@ -202,8 +195,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2173.  	while (sampleFrames-- > 0) {

  2174.  		double inputSampleL = *in1;

  2175.  		double inputSampleR = *in2;

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

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

  2178.  		double drySampleL = inputSampleL;

  2179.  		double drySampleR = inputSampleR;

  2180.  

  2181. @@ -373,20 +364,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2182.  		//end ClipOnly2 stereo as a little, compressed chunk that can be dropped into code

  2183.  		//final clip runs AFTER the Dry/Wet. It serves as a safety clip even if you're not full wet

  2184.  

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

  2186. -		int expon;

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

  2188. -		pressure5->fpdL ^= pressure5->fpdL << 13;

  2189. -		pressure5->fpdL ^= pressure5->fpdL >> 17;

  2190. -		pressure5->fpdL ^= pressure5->fpdL << 5;

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

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

  2193. -		pressure5->fpdR ^= pressure5->fpdR << 13;

  2194. -		pressure5->fpdR ^= pressure5->fpdR >> 17;

  2195. -		pressure5->fpdR ^= pressure5->fpdR << 5;

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

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

  2198. -

  2199.  		*out1 = (float) inputSampleL;

  2200.  		*out2 = (float) inputSampleR;

  2201.  

  2202. diff --git a/src/PurestConsole2Buss/PurestConsole2Buss.c b/src/PurestConsole2Buss/PurestConsole2Buss.c

  2203. index bd444ec..75c6163 100644

  2204. --- a/src/PurestConsole2Buss/PurestConsole2Buss.c

  2205. +++ b/src/PurestConsole2Buss/PurestConsole2Buss.c

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

  2207.  #include <math.h>

  2208.  #include <stdint.h>

  2209.  #include <stdlib.h>

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

  2211.  

  2212.  #define PURESTCONSOLE2BUSS_URI "https://hannesbraun.net/ns/lv2/airwindows/purestconsole2buss"

  2213.  

  2214. @@ -19,13 +20,11 @@ typedef struct {

  2215.  	float* output[2];

  2216.  

  2217.  	double biquadA[15];

  2218. -	uint32_t fpdL;

  2219. -	uint32_t fpdR;

  2220.  } PurestConsole2Buss;

  2221.  

  2222.  static LV2_Handle instantiate(

  2223.  	const LV2_Descriptor* descriptor,

  2224. -	double rate,

  2225. +	realdouble rate,

  2226.  	const char* bundle_path,

  2227.  	const LV2_Feature* const* features)

  2228.  {

  2229. @@ -60,10 +59,6 @@ static void activate(LV2_Handle instance)

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

  2231.  		purestConsole2Buss->biquadA[x] = 0.0;

  2232.  	}

  2233. -	purestConsole2Buss->fpdL = 1.0;

  2234. -	while (purestConsole2Buss->fpdL < 16386) purestConsole2Buss->fpdL = rand() * UINT32_MAX;

  2235. -	purestConsole2Buss->fpdR = 1.0;

  2236. -	while (purestConsole2Buss->fpdR < 16386) purestConsole2Buss->fpdR = rand() * UINT32_MAX;

  2237.  }

  2238.  

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

  2240. @@ -89,8 +84,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2241.  	while (sampleFrames-- > 0) {

  2242.  		double inputSampleL = *in1;

  2243.  		double inputSampleR = *in2;

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

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

  2246.  

  2247.  		if (purestConsole2Buss->biquadA[0] < 0.49999) {

  2248.  			double tempSample = purestConsole2Buss->biquadA[2] * inputSampleL + purestConsole2Buss->biquadA[3] * purestConsole2Buss->biquadA[7] + purestConsole2Buss->biquadA[4] * purestConsole2Buss->biquadA[8] - purestConsole2Buss->biquadA[5] * purestConsole2Buss->biquadA[9] - purestConsole2Buss->biquadA[6] * purestConsole2Buss->biquadA[10];

  2249. @@ -117,20 +110,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2250.  		inputSampleR = asin(inputSampleR);

  2251.  		//amplitude aspect

  2252.  

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

  2254. -		int expon;

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

  2256. -		purestConsole2Buss->fpdL ^= purestConsole2Buss->fpdL << 13;

  2257. -		purestConsole2Buss->fpdL ^= purestConsole2Buss->fpdL >> 17;

  2258. -		purestConsole2Buss->fpdL ^= purestConsole2Buss->fpdL << 5;

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

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

  2261. -		purestConsole2Buss->fpdR ^= purestConsole2Buss->fpdR << 13;

  2262. -		purestConsole2Buss->fpdR ^= purestConsole2Buss->fpdR >> 17;

  2263. -		purestConsole2Buss->fpdR ^= purestConsole2Buss->fpdR << 5;

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

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

  2266. -

  2267.  		*out1 = (float) inputSampleL;

  2268.  		*out2 = (float) inputSampleR;

  2269.  

  2270. diff --git a/src/PurestConsole2Channel/PurestConsole2Channel.c b/src/PurestConsole2Channel/PurestConsole2Channel.c

  2271. index 32c7d35..37ed37a 100644

  2272. --- a/src/PurestConsole2Channel/PurestConsole2Channel.c

  2273. +++ b/src/PurestConsole2Channel/PurestConsole2Channel.c

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

  2275.  #include <math.h>

  2276.  #include <stdint.h>

  2277.  #include <stdlib.h>

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

  2279.  

  2280.  #define PURESTCONSOLE2CHANNEL_URI "https://hannesbraun.net/ns/lv2/airwindows/purestconsole2channel"

  2281.  

  2282. @@ -19,13 +20,11 @@ typedef struct {

  2283.  	float* output[2];

  2284.  

  2285.  	double biquadA[15];

  2286. -	uint32_t fpdL;

  2287. -	uint32_t fpdR;

  2288.  } PurestConsole2Channel;

  2289.  

  2290.  static LV2_Handle instantiate(

  2291.  	const LV2_Descriptor* descriptor,

  2292. -	double rate,

  2293. +	realdouble rate,

  2294.  	const char* bundle_path,

  2295.  	const LV2_Feature* const* features)

  2296.  {

  2297. @@ -61,10 +60,6 @@ static void activate(LV2_Handle instance)

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

  2299.  		purestConsole2Channel->biquadA[x] = 0.0;

  2300.  	}

  2301. -	purestConsole2Channel->fpdL = 1.0;

  2302. -	while (purestConsole2Channel->fpdL < 16386) purestConsole2Channel->fpdL = rand() * UINT32_MAX;

  2303. -	purestConsole2Channel->fpdR = 1.0;

  2304. -	while (purestConsole2Channel->fpdR < 16386) purestConsole2Channel->fpdR = rand() * UINT32_MAX;

  2305.  }

  2306.  

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

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

  2309.  	while (sampleFrames-- > 0) {

  2310.  		double inputSampleL = *in1;

  2311.  		double inputSampleR = *in2;

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

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

  2314.  

  2315.  		if (purestConsole2Channel->biquadA[0] < 0.49999) {

  2316.  			double tempSample = purestConsole2Channel->biquadA[2] * inputSampleL + purestConsole2Channel->biquadA[3] * purestConsole2Channel->biquadA[7] + purestConsole2Channel->biquadA[4] * purestConsole2Channel->biquadA[8] - purestConsole2Channel->biquadA[5] * purestConsole2Channel->biquadA[9] - purestConsole2Channel->biquadA[6] * purestConsole2Channel->biquadA[10];

  2317. @@ -117,20 +110,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2318.  		inputSampleR = sin(inputSampleR);

  2319.  		//amplitude aspect

  2320.  

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

  2322. -		int expon;

  2323. -		frexpf((float)inputSampleL, &expon);

  2324. -		purestConsole2Channel->fpdL ^= purestConsole2Channel->fpdL << 13;

  2325. -		purestConsole2Channel->fpdL ^= purestConsole2Channel->fpdL >> 17;

  2326. -		purestConsole2Channel->fpdL ^= purestConsole2Channel->fpdL << 5;

  2327. -		inputSampleL += (((double)purestConsole2Channel->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2328. -		frexpf((float)inputSampleR, &expon);

  2329. -		purestConsole2Channel->fpdR ^= purestConsole2Channel->fpdR << 13;

  2330. -		purestConsole2Channel->fpdR ^= purestConsole2Channel->fpdR >> 17;

  2331. -		purestConsole2Channel->fpdR ^= purestConsole2Channel->fpdR << 5;

  2332. -		inputSampleR += (((double)purestConsole2Channel->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2333. -		//end 32 bit stereo floating point dither

  2334. -

  2335.  		*out1 = (float) inputSampleL;

  2336.  		*out2 = (float) inputSampleR;

  2337.  

  2338. diff --git a/src/PurestConsoleBuss/PurestConsoleBuss.c b/src/PurestConsoleBuss/PurestConsoleBuss.c

  2339. index 663ff27..a8cae6c 100644

  2340. --- a/src/PurestConsoleBuss/PurestConsoleBuss.c

  2341. +++ b/src/PurestConsoleBuss/PurestConsoleBuss.c

  2342. @@ -3,6 +3,7 @@

  2343.  #include <math.h>

  2344.  #include <stdint.h>

  2345.  #include <stdlib.h>

  2346. +#include "../optimizations.h"

  2347.  

  2348.  #define PURESTCONSOLEBUSS_URI "https://hannesbraun.net/ns/lv2/airwindows/purestconsolebuss"

  2349.  

  2350. @@ -16,14 +17,11 @@ typedef enum {

  2351.  typedef struct {

  2352.  	const float* input[2];

  2353.  	float* output[2];

  2354. -

  2355. -	uint32_t fpdL;

  2356. -	uint32_t fpdR;

  2357.  } PurestConsoleBuss;

  2358.  

  2359.  static LV2_Handle instantiate(

  2360.  	const LV2_Descriptor* descriptor,

  2361. -	double rate,

  2362. +	realdouble rate,

  2363.  	const char* bundle_path,

  2364.  	const LV2_Feature* const* features)

  2365.  {

  2366. @@ -51,14 +49,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  2367.  	}

  2368.  }

  2369.  

  2370. -static void activate(LV2_Handle instance)

  2371. -{

  2372. -	PurestConsoleBuss* purestConsoleBuss = (PurestConsoleBuss*) instance;

  2373. -	purestConsoleBuss->fpdL = 1.0;

  2374. -	while (purestConsoleBuss->fpdL < 16386) purestConsoleBuss->fpdL = rand() * UINT32_MAX;

  2375. -	purestConsoleBuss->fpdR = 1.0;

  2376. -	while (purestConsoleBuss->fpdR < 16386) purestConsoleBuss->fpdR = rand() * UINT32_MAX;

  2377. -}

  2378. +static void activate(LV2_Handle instance) {}

  2379.  

  2380.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2381.  {

  2382. @@ -75,8 +66,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2383.  	while (sampleFrames-- > 0) {

  2384.  		inputSampleL = *in1;

  2385.  		inputSampleR = *in2;

  2386. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = purestConsoleBuss->fpdL * 1.18e-17;

  2387. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = purestConsoleBuss->fpdR * 1.18e-17;

  2388.  

  2389.  		if (inputSampleL > 1.0) inputSampleL = 1.0;

  2390.  		if (inputSampleL < -1.0) inputSampleL = -1.0;

  2391. @@ -88,20 +77,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2392.  		inputSampleR = asin(inputSampleR);

  2393.  		//amplitude aspect

  2394.  

  2395. -		//begin 32 bit stereo floating point dither

  2396. -		int expon;

  2397. -		frexpf((float)inputSampleL, &expon);

  2398. -		purestConsoleBuss->fpdL ^= purestConsoleBuss->fpdL << 13;

  2399. -		purestConsoleBuss->fpdL ^= purestConsoleBuss->fpdL >> 17;

  2400. -		purestConsoleBuss->fpdL ^= purestConsoleBuss->fpdL << 5;

  2401. -		inputSampleL += (((double)purestConsoleBuss->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2402. -		frexpf((float)inputSampleR, &expon);

  2403. -		purestConsoleBuss->fpdR ^= purestConsoleBuss->fpdR << 13;

  2404. -		purestConsoleBuss->fpdR ^= purestConsoleBuss->fpdR >> 17;

  2405. -		purestConsoleBuss->fpdR ^= purestConsoleBuss->fpdR << 5;

  2406. -		inputSampleR += (((double)purestConsoleBuss->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2407. -		//end 32 bit stereo floating point dither

  2408. -

  2409.  		*out1 = (float) inputSampleL;

  2410.  		*out2 = (float) inputSampleR;

  2411.  

  2412. diff --git a/src/PurestConsoleChannel/PurestConsoleChannel.c b/src/PurestConsoleChannel/PurestConsoleChannel.c

  2413. index cee3dcb..29f7187 100644

  2414. --- a/src/PurestConsoleChannel/PurestConsoleChannel.c

  2415. +++ b/src/PurestConsoleChannel/PurestConsoleChannel.c

  2416. @@ -3,6 +3,7 @@

  2417.  #include <math.h>

  2418.  #include <stdint.h>

  2419.  #include <stdlib.h>

  2420. +#include "../optimizations.h"

  2421.  

  2422.  #define PURESTCONSOLECHANNEL_URI "https://hannesbraun.net/ns/lv2/airwindows/purestconsolechannel"

  2423.  

  2424. @@ -16,14 +17,11 @@ typedef enum {

  2425.  typedef struct {

  2426.  	const float* input[2];

  2427.  	float* output[2];

  2428. -

  2429. -	uint32_t fpdL;

  2430. -	uint32_t fpdR;

  2431.  } PurestConsoleChannel;

  2432.  

  2433.  static LV2_Handle instantiate(

  2434.  	const LV2_Descriptor* descriptor,

  2435. -	double rate,

  2436. +	realdouble rate,

  2437.  	const char* bundle_path,

  2438.  	const LV2_Feature* const* features)

  2439.  {

  2440. @@ -51,14 +49,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  2441.  	}

  2442.  }

  2443.  

  2444. -static void activate(LV2_Handle instance)

  2445. -{

  2446. -	PurestConsoleChannel* purestConsoleChannel = (PurestConsoleChannel*) instance;

  2447. -	purestConsoleChannel->fpdL = 1.0;

  2448. -	while (purestConsoleChannel->fpdL < 16386) purestConsoleChannel->fpdL = rand() * UINT32_MAX;

  2449. -	purestConsoleChannel->fpdR = 1.0;

  2450. -	while (purestConsoleChannel->fpdR < 16386) purestConsoleChannel->fpdR = rand() * UINT32_MAX;

  2451. -}

  2452. +static void activate(LV2_Handle instance) {}

  2453.  

  2454.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2455.  {

  2456. @@ -75,27 +66,11 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2457.  	while (sampleFrames-- > 0) {

  2458.  		inputSampleL = *in1;

  2459.  		inputSampleR = *in2;

  2460. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = purestConsoleChannel->fpdL * 1.18e-17;

  2461. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = purestConsoleChannel->fpdR * 1.18e-17;

  2462.  

  2463.  		inputSampleL = sin(inputSampleL);

  2464.  		inputSampleR = sin(inputSampleR);

  2465.  		//amplitude aspect

  2466.  

  2467. -		//begin 32 bit stereo floating point dither

  2468. -		int expon;

  2469. -		frexpf((float)inputSampleL, &expon);

  2470. -		purestConsoleChannel->fpdL ^= purestConsoleChannel->fpdL << 13;

  2471. -		purestConsoleChannel->fpdL ^= purestConsoleChannel->fpdL >> 17;

  2472. -		purestConsoleChannel->fpdL ^= purestConsoleChannel->fpdL << 5;

  2473. -		inputSampleL += (((double)purestConsoleChannel->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2474. -		frexpf((float)inputSampleR, &expon);

  2475. -		purestConsoleChannel->fpdR ^= purestConsoleChannel->fpdR << 13;

  2476. -		purestConsoleChannel->fpdR ^= purestConsoleChannel->fpdR >> 17;

  2477. -		purestConsoleChannel->fpdR ^= purestConsoleChannel->fpdR << 5;

  2478. -		inputSampleR += (((double)purestConsoleChannel->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2479. -		//end 32 bit stereo floating point dither

  2480. -

  2481.  		*out1 = (float) inputSampleL;

  2482.  		*out2 = (float) inputSampleR;

  2483.  

  2484. diff --git a/src/PurestGain/PurestGain.c b/src/PurestGain/PurestGain.c

  2485. index 97c9925..b06c780 100644

  2486. --- a/src/PurestGain/PurestGain.c

  2487. +++ b/src/PurestGain/PurestGain.c

  2488. @@ -3,6 +3,7 @@

  2489.  #include <math.h>

  2490.  #include <stdint.h>

  2491.  #include <stdlib.h>

  2492. +#include "../optimizations.h"

  2493.  

  2494.  #define PURESTGAIN_URI "https://hannesbraun.net/ns/lv2/airwindows/purestgain"

  2495.  

  2496. @@ -22,8 +23,6 @@ typedef struct {

  2497.  	const float* gain;

  2498.  	const float* slowfade;

  2499.  

  2500. -	uint32_t fpdL;

  2501. -	uint32_t fpdR;

  2502.  	//default stuff

  2503.  	double gainchase;

  2504.  	double settingchase;

  2505. @@ -33,7 +32,7 @@ typedef struct {

  2506.  

  2507.  static LV2_Handle instantiate(

  2508.  	const LV2_Descriptor* descriptor,

  2509. -	double rate,

  2510. +	realdouble rate,

  2511.  	const char* bundle_path,

  2512.  	const LV2_Feature* const* features)

  2513.  {

  2514. @@ -75,10 +74,6 @@ static void activate(LV2_Handle instance)

  2515.  	purestGain->settingchase = -90.0;

  2516.  	purestGain->gainBchase = -90.0;

  2517.  	purestGain->chasespeed = 350.0;

  2518. -	purestGain->fpdL = 1.0;

  2519. -	while (purestGain->fpdL < 16386) purestGain->fpdL = rand() * UINT32_MAX;

  2520. -	purestGain->fpdR = 1.0;

  2521. -	while (purestGain->fpdR < 16386) purestGain->fpdR = rand() * UINT32_MAX;

  2522.  }

  2523.  

  2524.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2525. @@ -143,8 +138,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2526.  

  2527.  		inputSampleL = *in1;

  2528.  		inputSampleR = *in2;

  2529. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = purestGain->fpdL * 1.18e-17;

  2530. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = purestGain->fpdR * 1.18e-17;

  2531.  

  2532.  		if (1.0 == outputgain) {

  2533.  			*out1 = *in1;

  2534. @@ -152,19 +145,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2535.  		} else {

  2536.  			inputSampleL *= outputgain;

  2537.  			inputSampleR *= outputgain;

  2538. -			//begin 32 bit stereo floating point dither

  2539. -			int expon;

  2540. -			frexpf((float)inputSampleL, &expon);

  2541. -			purestGain->fpdL ^= purestGain->fpdL << 13;

  2542. -			purestGain->fpdL ^= purestGain->fpdL >> 17;

  2543. -			purestGain->fpdL ^= purestGain->fpdL << 5;

  2544. -			inputSampleL += (((double)purestGain->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2545. -			frexpf((float)inputSampleR, &expon);

  2546. -			purestGain->fpdR ^= purestGain->fpdR << 13;

  2547. -			purestGain->fpdR ^= purestGain->fpdR >> 17;

  2548. -			purestGain->fpdR ^= purestGain->fpdR << 5;

  2549. -			inputSampleR += (((double)purestGain->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2550. -			//end 32 bit stereo floating point dither

  2551.  			*out1 = (float) inputSampleL;

  2552.  			*out2 = (float) inputSampleR;

  2553.  		}

  2554. diff --git a/src/RightoMono/RightoMono.c b/src/RightoMono/RightoMono.c

  2555. index 03acec9..e6eb262 100644

  2556. --- a/src/RightoMono/RightoMono.c

  2557. +++ b/src/RightoMono/RightoMono.c

  2558. @@ -2,6 +2,7 @@

  2559.  

  2560.  #include <stdint.h>

  2561.  #include <stdlib.h>

  2562. +#include "../optimizations.h"

  2563.  

  2564.  #define RIGHTOMONO_URI "https://hannesbraun.net/ns/lv2/airwindows/rightomono"

  2565.  

  2566. @@ -22,7 +23,7 @@ typedef struct {

  2567.  

  2568.  static LV2_Handle instantiate(

  2569.  	const LV2_Descriptor* descriptor,

  2570. -	double rate,

  2571. +	realdouble rate,

  2572.  	const char* bundle_path,

  2573.  	const LV2_Feature* const* features)

  2574.  {

  2575. diff --git a/src/Sidepass/Sidepass.c b/src/Sidepass/Sidepass.c

  2576. index 0e99ac2..b55f35e 100644

  2577. --- a/src/Sidepass/Sidepass.c

  2578. +++ b/src/Sidepass/Sidepass.c

  2579. @@ -4,6 +4,7 @@

  2580.  #include <stdbool.h>

  2581.  #include <stdint.h>

  2582.  #include <stdlib.h>

  2583. +#include "../optimizations.h"

  2584.  

  2585.  #define SIDEPASS_URI "https://hannesbraun.net/ns/lv2/airwindows/sidepass"

  2586.  

  2587. @@ -24,13 +25,11 @@ typedef struct {

  2588.  	double iirSampleA;

  2589.  	double iirSampleB;

  2590.  	bool flip;

  2591. -	uint32_t fpdL;

  2592. -	uint32_t fpdR;

  2593.  } Sidepass;

  2594.  

  2595.  static LV2_Handle instantiate(

  2596.  	const LV2_Descriptor* descriptor,

  2597. -	double rate,

  2598. +	realdouble rate,

  2599.  	const char* bundle_path,

  2600.  	const LV2_Feature* const* features)

  2601.  {

  2602. @@ -68,10 +67,6 @@ static void activate(LV2_Handle instance)

  2603.  	sidepass->iirSampleA = 0.0;

  2604.  	sidepass->iirSampleB = 0.0;

  2605.  	sidepass->flip = true;

  2606. -	sidepass->fpdL = 1.0;

  2607. -	while (sidepass->fpdL < 16386) sidepass->fpdL = rand() * UINT32_MAX;

  2608. -	sidepass->fpdR = 1.0;

  2609. -	while (sidepass->fpdR < 16386) sidepass->fpdR = rand() * UINT32_MAX;

  2610.  }

  2611.  

  2612.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2613. @@ -95,8 +90,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2614.  	while (sampleFrames-- > 0) {

  2615.  		inputSampleL = *in1;

  2616.  		inputSampleR = *in2;

  2617. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = sidepass->fpdL * 1.18e-17;

  2618. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = sidepass->fpdR * 1.18e-17;

  2619.  

  2620.  		mid = inputSampleL + inputSampleR;

  2621.  		side = inputSampleL - inputSampleR;

  2622. @@ -113,20 +106,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2623.  		inputSampleL = (mid + side) / 2.0;

  2624.  		inputSampleR = (mid - side) / 2.0;

  2625.  

  2626. -		//begin 32 bit stereo floating point dither

  2627. -		int expon;

  2628. -		frexpf((float)inputSampleL, &expon);

  2629. -		sidepass->fpdL ^= sidepass->fpdL << 13;

  2630. -		sidepass->fpdL ^= sidepass->fpdL >> 17;

  2631. -		sidepass->fpdL ^= sidepass->fpdL << 5;

  2632. -		inputSampleL += (((double)sidepass->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2633. -		frexpf((float)inputSampleR, &expon);

  2634. -		sidepass->fpdR ^= sidepass->fpdR << 13;

  2635. -		sidepass->fpdR ^= sidepass->fpdR >> 17;

  2636. -		sidepass->fpdR ^= sidepass->fpdR << 5;

  2637. -		inputSampleR += (((double)sidepass->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2638. -		//end 32 bit stereo floating point dither

  2639. -

  2640.  		*out1 = (float) inputSampleL;

  2641.  		*out2 = (float) inputSampleR;

  2642.  

  2643. diff --git a/src/Slew/Slew.c b/src/Slew/Slew.c

  2644. index e23ad33..97fe6fb 100644

  2645. --- a/src/Slew/Slew.c

  2646. +++ b/src/Slew/Slew.c

  2647. @@ -3,6 +3,7 @@

  2648.  #include <math.h>

  2649.  #include <stdint.h>

  2650.  #include <stdlib.h>

  2651. +#include "../optimizations.h"

  2652.  

  2653.  #define SLEW_URI "https://hannesbraun.net/ns/lv2/airwindows/slew"

  2654.  

  2655. @@ -22,13 +23,11 @@ typedef struct {

  2656.  

  2657.  	double lastSampleL;

  2658.  	double lastSampleR;

  2659. -	uint32_t fpdL;

  2660. -	uint32_t fpdR;

  2661.  } Slew;

  2662.  

  2663.  static LV2_Handle instantiate(

  2664.  	const LV2_Descriptor* descriptor,

  2665. -	double rate,

  2666. +	realdouble rate,

  2667.  	const char* bundle_path,

  2668.  	const LV2_Feature* const* features)

  2669.  {

  2670. @@ -91,8 +90,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2671.  	while (sampleFrames-- > 0) {

  2672.  		inputSampleL = *in1;

  2673.  		inputSampleR = *in2;

  2674. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = slew->fpdL * 1.18e-17;

  2675. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = slew->fpdR * 1.18e-17;

  2676.  

  2677.  		clamp = inputSampleL - slew->lastSampleL;

  2678.  		outputSampleL = inputSampleL;

  2679. diff --git a/src/SlewOnly/SlewOnly.c b/src/SlewOnly/SlewOnly.c

  2680. index d2de114..2148b40 100644

  2681. --- a/src/SlewOnly/SlewOnly.c

  2682. +++ b/src/SlewOnly/SlewOnly.c

  2683. @@ -3,6 +3,7 @@

  2684.  #include <math.h>

  2685.  #include <stdint.h>

  2686.  #include <stdlib.h>

  2687. +#include "../optimizations.h"

  2688.  

  2689.  #define SLEWONLY_URI "https://hannesbraun.net/ns/lv2/airwindows/slewonly"

  2690.  

  2691. @@ -19,13 +20,11 @@ typedef struct {

  2692.  

  2693.  	double lastSampleL;

  2694.  	double lastSampleR;

  2695. -	uint32_t fpdL;

  2696. -	uint32_t fpdR;

  2697.  } SlewOnly;

  2698.  

  2699.  static LV2_Handle instantiate(

  2700.  	const LV2_Descriptor* descriptor,

  2701. -	double rate,

  2702. +	realdouble rate,

  2703.  	const char* bundle_path,

  2704.  	const LV2_Feature* const* features)

  2705.  {

  2706. @@ -78,8 +77,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2707.  	while (sampleFrames-- > 0) {

  2708.  		inputSampleL = *in1;

  2709.  		inputSampleR = *in2;

  2710. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = slewonly->fpdL * 1.18e-17;

  2711. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = slewonly->fpdR * 1.18e-17;

  2712.  

  2713.  		outputSampleL = (inputSampleL - slewonly->lastSampleL) * trim;

  2714.  		outputSampleR = (inputSampleR - slewonly->lastSampleR) * trim;

  2715. diff --git a/src/Spiral/Spiral.c b/src/Spiral/Spiral.c

  2716. index 05a3361..619b60a 100644

  2717. --- a/src/Spiral/Spiral.c

  2718. +++ b/src/Spiral/Spiral.c

  2719. @@ -3,6 +3,7 @@

  2720.  #include <math.h>

  2721.  #include <stdint.h>

  2722.  #include <stdlib.h>

  2723. +#include "../optimizations.h"

  2724.  

  2725.  #define SPIRAL_URI "https://hannesbraun.net/ns/lv2/airwindows/spiral"

  2726.  

  2727. @@ -16,14 +17,11 @@ typedef enum {

  2728.  typedef struct {

  2729.  	const float* input[2];

  2730.  	float* output[2];

  2731. -

  2732. -	uint32_t fpdL;

  2733. -	uint32_t fpdR;

  2734.  } Spiral;

  2735.  

  2736.  static LV2_Handle instantiate(

  2737.  	const LV2_Descriptor* descriptor,

  2738. -	double rate,

  2739. +	realdouble rate,

  2740.  	const char* bundle_path,

  2741.  	const LV2_Feature* const* features)

  2742.  {

  2743. @@ -51,14 +49,7 @@ static void connect_port(LV2_Handle instance, uint32_t port, void* data)

  2744.  	}

  2745.  }

  2746.  

  2747. -static void activate(LV2_Handle instance)

  2748. -{

  2749. -	Spiral* spiral = (Spiral*) instance;

  2750. -	spiral->fpdL = 1.0;

  2751. -	while (spiral->fpdL < 16386) spiral->fpdL = rand() * UINT32_MAX;

  2752. -	spiral->fpdR = 1.0;

  2753. -	while (spiral->fpdR < 16386) spiral->fpdR = rand() * UINT32_MAX;

  2754. -}

  2755. +static void activate(LV2_Handle instance) {}

  2756.  

  2757.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2758.  {

  2759. @@ -73,28 +64,10 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2760.  		double inputSampleL = *in1;

  2761.  		double inputSampleR = *in2;

  2762.  

  2763. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = spiral->fpdL * 1.18e-17;

  2764. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = spiral->fpdR * 1.18e-17;

  2765. -

  2766. -

  2767.  		//clip to 1.2533141373155 to reach maximum output

  2768.  		inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((fabs(inputSampleL) == 0.0) ? 1 : fabs(inputSampleL));

  2769.  		inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((fabs(inputSampleR) == 0.0) ? 1 : fabs(inputSampleR));

  2770.  

  2771. -		//begin 32 bit stereo floating point dither

  2772. -		int expon;

  2773. -		frexpf((float)inputSampleL, &expon);

  2774. -		spiral->fpdL ^= spiral->fpdL << 13;

  2775. -		spiral->fpdL ^= spiral->fpdL >> 17;

  2776. -		spiral->fpdL ^= spiral->fpdL << 5;

  2777. -		inputSampleL += (((double)spiral->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2778. -		frexpf((float)inputSampleR, &expon);

  2779. -		spiral->fpdR ^= spiral->fpdR << 13;

  2780. -		spiral->fpdR ^= spiral->fpdR >> 17;

  2781. -		spiral->fpdR ^= spiral->fpdR << 5;

  2782. -		inputSampleR += (((double)spiral->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2783. -		//end 32 bit stereo floating point dither

  2784. -

  2785.  		*out1 = (float) inputSampleL;

  2786.  		*out2 = (float) inputSampleR;

  2787.  

  2788. diff --git a/src/StarChild/StarChild.c b/src/StarChild/StarChild.c

  2789. index c45f41d..e16f9ef 100644

  2790. --- a/src/StarChild/StarChild.c

  2791. +++ b/src/StarChild/StarChild.c

  2792. @@ -3,6 +3,7 @@

  2793.  #include <math.h>

  2794.  #include <stdint.h>

  2795.  #include <stdlib.h>

  2796. +#include "../optimizations.h"

  2797.  

  2798.  #define STARCHILD_URI "https://hannesbraun.net/ns/lv2/airwindows/starchild"

  2799.  

  2800. @@ -23,8 +24,6 @@ typedef struct {

  2801.  	const float* grain;

  2802.  	const float* drywet;

  2803.  

  2804. -	uint32_t fpdL;

  2805. -	uint32_t fpdR;

  2806.  	//default stuff

  2807.  

  2808.  	double d[45102];

  2809. @@ -47,7 +46,7 @@ typedef struct {

  2810.  

  2811.  static LV2_Handle instantiate(

  2812.  	const LV2_Descriptor* descriptor,

  2813. -	double rate,

  2814. +	realdouble rate,

  2815.  	const char* bundle_path,

  2816.  	const LV2_Feature* const* features)

  2817.  {

  2818. @@ -303,11 +302,6 @@ static void activate(LV2_Handle instance)

  2819.  	starchild->pitchCounter = 2;

  2820.  	starchild->increment = 1;

  2821.  	starchild->dutyCycle = 1;

  2822. -

  2823. -	starchild->fpdL = 1.0;

  2824. -	while (starchild->fpdL < 16386) starchild->fpdL = rand() * UINT32_MAX;

  2825. -	starchild->fpdR = 1.0;

  2826. -	while (starchild->fpdR < 16386) starchild->fpdR = rand() * UINT32_MAX;

  2827.  }

  2828.  

  2829.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2830. @@ -437,8 +431,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2831.  	while (sampleFrames-- > 0) {

  2832.  		inputSampleL = *in1;

  2833.  		inputSampleR = *in2;

  2834. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = starchild->fpdL * 1.18e-17;

  2835. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = starchild->fpdR * 1.18e-17;

  2836.  		drySampleL = inputSampleL;

  2837.  		drySampleR = inputSampleR;

  2838.  

  2839. @@ -1058,20 +1050,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2840.  		inputSampleR += drySampleR;

  2841.  		//here we combine the tanks with the dry signal

  2842.  

  2843. -		//begin 32 bit stereo floating point dither

  2844. -		int expon;

  2845. -		frexpf((float)inputSampleL, &expon);

  2846. -		starchild->fpdL ^= starchild->fpdL << 13;

  2847. -		starchild->fpdL ^= starchild->fpdL >> 17;

  2848. -		starchild->fpdL ^= starchild->fpdL << 5;

  2849. -		inputSampleL += (((double)starchild->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2850. -		frexpf((float)inputSampleR, &expon);

  2851. -		starchild->fpdR ^= starchild->fpdR << 13;

  2852. -		starchild->fpdR ^= starchild->fpdR >> 17;

  2853. -		starchild->fpdR ^= starchild->fpdR << 5;

  2854. -		inputSampleR += (((double)starchild->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2855. -		//end 32 bit stereo floating point dither

  2856. -

  2857.  		*out1 = (float) inputSampleL;

  2858.  		*out2 = (float) inputSampleR;

  2859.  

  2860. diff --git a/src/SubsOnly/SubsOnly.c b/src/SubsOnly/SubsOnly.c

  2861. index 886ae22..e649464 100644

  2862. --- a/src/SubsOnly/SubsOnly.c

  2863. +++ b/src/SubsOnly/SubsOnly.c

  2864. @@ -3,6 +3,7 @@

  2865.  #include <math.h>

  2866.  #include <stdint.h>

  2867.  #include <stdlib.h>

  2868. +#include "../optimizations.h"

  2869.  

  2870.  #define SUBSONLY_URI "https://hannesbraun.net/ns/lv2/airwindows/subsonly"

  2871.  

  2872. @@ -71,14 +72,11 @@ typedef struct {

  2873.  	double iirSampleXR;

  2874.  	double iirSampleYR;

  2875.  	double iirSampleZR;

  2876. -

  2877. -	uint32_t fpdL;

  2878. -	uint32_t fpdR;

  2879.  } SubsOnly;

  2880.  

  2881.  static LV2_Handle instantiate(

  2882.  	const LV2_Descriptor* descriptor,

  2883. -	double rate,

  2884. +	realdouble rate,

  2885.  	const char* bundle_path,

  2886.  	const LV2_Feature* const* features)

  2887.  {

  2888. @@ -188,8 +186,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2889.  	while (sampleFrames-- > 0) {

  2890.  		inputSampleL = *in1;

  2891.  		inputSampleR = *in2;

  2892. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = subsonly->fpdL * 1.18e-17;

  2893. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = subsonly->fpdR * 1.18e-17;

  2894.  

  2895.  		gain = gaintarget;

  2896.  

  2897. diff --git a/src/Ultrasonic/Ultrasonic.c b/src/Ultrasonic/Ultrasonic.c

  2898. index 57b2398..6aacb3b 100644

  2899. --- a/src/Ultrasonic/Ultrasonic.c

  2900. +++ b/src/Ultrasonic/Ultrasonic.c

  2901. @@ -3,6 +3,7 @@

  2902.  #include <math.h>

  2903.  #include <stdint.h>

  2904.  #include <stdlib.h>

  2905. +#include "../optimizations.h"

  2906.  

  2907.  #define ULTRASONIC_URI "https://hannesbraun.net/ns/lv2/airwindows/ultrasonic"

  2908.  

  2909. @@ -22,13 +23,11 @@ typedef struct {

  2910.  	double biquadC[15];

  2911.  	double biquadD[15];

  2912.  	double biquadE[15];

  2913. -	uint32_t fpdL;

  2914. -	uint32_t fpdR;

  2915.  } Ultrasonic;

  2916.  

  2917.  static LV2_Handle instantiate(

  2918.  	const LV2_Descriptor* descriptor,

  2919. -	double rate,

  2920. +	realdouble rate,

  2921.  	const char* bundle_path,

  2922.  	const LV2_Feature* const* features)

  2923.  {

  2924. @@ -66,10 +65,6 @@ static void activate(LV2_Handle instance)

  2925.  		ultrasonic->biquadD[x] = 0.0;

  2926.  		ultrasonic->biquadE[x] = 0.0;

  2927.  	}

  2928. -	ultrasonic->fpdL = 1.0;

  2929. -	while (ultrasonic->fpdL < 16386) ultrasonic->fpdL = rand() * UINT32_MAX;

  2930. -	ultrasonic->fpdR = 1.0;

  2931. -	while (ultrasonic->fpdR < 16386) ultrasonic->fpdR = rand() * UINT32_MAX;

  2932.  }

  2933.  

  2934.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  2935. @@ -136,8 +131,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2936.  	while (sampleFrames-- > 0) {

  2937.  		double inputSampleL = *in1;

  2938.  		double inputSampleR = *in2;

  2939. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = ultrasonic->fpdL * 1.18e-17;

  2940. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = ultrasonic->fpdR * 1.18e-17;

  2941.  

  2942.  		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];

  2943.  		ultrasonic->biquadA[8] = ultrasonic->biquadA[7];

  2944. @@ -209,20 +202,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  2945.  		ultrasonic->biquadE[14] = ultrasonic->biquadE[13];

  2946.  		ultrasonic->biquadE[13] = inputSampleR; // DF1 right

  2947.  

  2948. -		//begin 32 bit stereo floating point dither

  2949. -		int expon;

  2950. -		frexpf((float)inputSampleL, &expon);

  2951. -		ultrasonic->fpdL ^= ultrasonic->fpdL << 13;

  2952. -		ultrasonic->fpdL ^= ultrasonic->fpdL >> 17;

  2953. -		ultrasonic->fpdL ^= ultrasonic->fpdL << 5;

  2954. -		inputSampleL += (((double)ultrasonic->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2955. -		frexpf((float)inputSampleR, &expon);

  2956. -		ultrasonic->fpdR ^= ultrasonic->fpdR << 13;

  2957. -		ultrasonic->fpdR ^= ultrasonic->fpdR >> 17;

  2958. -		ultrasonic->fpdR ^= ultrasonic->fpdR << 5;

  2959. -		inputSampleR += (((double)ultrasonic->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  2960. -		//end 32 bit stereo floating point dither

  2961. -

  2962.  		*out1 = (float) inputSampleL;

  2963.  		*out2 = (float) inputSampleR;

  2964.  

  2965. diff --git a/src/UltrasonicLite/UltrasonicLite.c b/src/UltrasonicLite/UltrasonicLite.c

  2966. index f73d3f5..ef932f2 100644

  2967. --- a/src/UltrasonicLite/UltrasonicLite.c

  2968. +++ b/src/UltrasonicLite/UltrasonicLite.c

  2969. @@ -3,6 +3,7 @@

  2970.  #include <math.h>

  2971.  #include <stdint.h>

  2972.  #include <stdlib.h>

  2973. +#include "../optimizations.h"

  2974.  

  2975.  #define ULTRASONICLITE_URI "https://hannesbraun.net/ns/lv2/airwindows/ultrasoniclite"

  2976.  

  2977. @@ -18,13 +19,11 @@ typedef struct {

  2978.  	const float* input[2];

  2979.  	float* output[2];

  2980.  	double biquadA[15];

  2981. -	uint32_t fpdL;

  2982. -	uint32_t fpdR;

  2983.  } UltrasonicLite;

  2984.  

  2985.  static LV2_Handle instantiate(

  2986.  	const LV2_Descriptor* descriptor,

  2987. -	double rate,

  2988. +	realdouble rate,

  2989.  	const char* bundle_path,

  2990.  	const LV2_Feature* const* features)

  2991.  {

  2992. @@ -58,10 +57,6 @@ static void activate(LV2_Handle instance)

  2993.  	for (int x = 0; x < 15; x++) {

  2994.  		ultrasoniclite->biquadA[x] = 0.0;

  2995.  	}

  2996. -	ultrasoniclite->fpdL = 1.0;

  2997. -	while (ultrasoniclite->fpdL < 16386) ultrasoniclite->fpdL = rand() * UINT32_MAX;

  2998. -	ultrasoniclite->fpdR = 1.0;

  2999. -	while (ultrasoniclite->fpdR < 16386) ultrasoniclite->fpdR = rand() * UINT32_MAX;

  3000.  }

  3001.  

  3002.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  3003. @@ -90,8 +85,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  3004.  	while (sampleFrames-- > 0) {

  3005.  		double inputSampleL = *in1;

  3006.  		double inputSampleR = *in2;

  3007. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = ultrasoniclite->fpdL * 1.18e-17;

  3008. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = ultrasoniclite->fpdR * 1.18e-17;

  3009.  

  3010.  		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];

  3011.  		ultrasoniclite->biquadA[8] = ultrasoniclite->biquadA[7];

  3012. @@ -107,20 +100,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  3013.  		ultrasoniclite->biquadA[14] = ultrasoniclite->biquadA[13];

  3014.  		ultrasoniclite->biquadA[13] = inputSampleR; //DF1 right

  3015.  

  3016. -		//begin 32 bit stereo floating point dither

  3017. -		int expon;

  3018. -		frexpf((float)inputSampleL, &expon);

  3019. -		ultrasoniclite->fpdL ^= ultrasoniclite->fpdL << 13;

  3020. -		ultrasoniclite->fpdL ^= ultrasoniclite->fpdL >> 17;

  3021. -		ultrasoniclite->fpdL ^= ultrasoniclite->fpdL << 5;

  3022. -		inputSampleL += (((double)ultrasoniclite->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  3023. -		frexpf((float)inputSampleR, &expon);

  3024. -		ultrasoniclite->fpdR ^= ultrasoniclite->fpdR << 13;

  3025. -		ultrasoniclite->fpdR ^= ultrasoniclite->fpdR >> 17;

  3026. -		ultrasoniclite->fpdR ^= ultrasoniclite->fpdR << 5;

  3027. -		inputSampleR += (((double)ultrasoniclite->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  3028. -		//end 32 bit stereo floating point dither

  3029. -

  3030.  		*out1 = (float) inputSampleL;

  3031.  		*out2 = (float) inputSampleR;

  3032.  

  3033. diff --git a/src/UltrasonicMed/UltrasonicMed.c b/src/UltrasonicMed/UltrasonicMed.c

  3034. index 2231697..29d91be 100644

  3035. --- a/src/UltrasonicMed/UltrasonicMed.c

  3036. +++ b/src/UltrasonicMed/UltrasonicMed.c

  3037. @@ -3,6 +3,7 @@

  3038.  #include <math.h>

  3039.  #include <stdint.h>

  3040.  #include <stdlib.h>

  3041. +#include "../optimizations.h"

  3042.  

  3043.  #define ULTRASONICMED_URI "https://hannesbraun.net/ns/lv2/airwindows/ultrasonicmed"

  3044.  

  3045. @@ -19,13 +20,11 @@ typedef struct {

  3046.  	float* output[2];

  3047.  	double biquadA[15];

  3048.  	double biquadB[15];

  3049. -	uint32_t fpdL;

  3050. -	uint32_t fpdR;

  3051.  } UltrasonicMed;

  3052.  

  3053.  static LV2_Handle instantiate(

  3054.  	const LV2_Descriptor* descriptor,

  3055. -	double rate,

  3056. +	realdouble rate,

  3057.  	const char* bundle_path,

  3058.  	const LV2_Feature* const* features)

  3059.  {

  3060. @@ -61,10 +60,6 @@ static void activate(LV2_Handle instance)

  3061.  		ultrasonicMed->biquadA[x] = 0.0;

  3062.  		ultrasonicMed->biquadB[x] = 0.0;

  3063.  	}

  3064. -	ultrasonicMed->fpdL = 1.0;

  3065. -	while (ultrasonicMed->fpdL < 16386) ultrasonicMed->fpdL = rand() * UINT32_MAX;

  3066. -	ultrasonicMed->fpdR = 1.0;

  3067. -	while (ultrasonicMed->fpdR < 16386) ultrasonicMed->fpdR = rand() * UINT32_MAX;

  3068.  }

  3069.  

  3070.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  3071. @@ -105,8 +100,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  3072.  	while (sampleFrames-- > 0) {

  3073.  		double inputSampleL = *in1;

  3074.  		double inputSampleR = *in2;

  3075. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = ultrasonicMed->fpdL * 1.18e-17;

  3076. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = ultrasonicMed->fpdR * 1.18e-17;

  3077.  

  3078.  		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];

  3079.  		ultrasonicMed->biquadA[8] = ultrasonicMed->biquadA[7];

  3080. @@ -136,20 +129,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  3081.  		ultrasonicMed->biquadB[14] = ultrasonicMed->biquadB[13];

  3082.  		ultrasonicMed->biquadB[13] = inputSampleR; //DF1 right

  3083.  

  3084. -		//begin 32 bit stereo floating point dither

  3085. -		int expon;

  3086. -		frexpf((float)inputSampleL, &expon);

  3087. -		ultrasonicMed->fpdL ^= ultrasonicMed->fpdL << 13;

  3088. -		ultrasonicMed->fpdL ^= ultrasonicMed->fpdL >> 17;

  3089. -		ultrasonicMed->fpdL ^= ultrasonicMed->fpdL << 5;

  3090. -		inputSampleL += (((double)ultrasonicMed->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  3091. -		frexpf((float)inputSampleR, &expon);

  3092. -		ultrasonicMed->fpdR ^= ultrasonicMed->fpdR << 13;

  3093. -		ultrasonicMed->fpdR ^= ultrasonicMed->fpdR >> 17;

  3094. -		ultrasonicMed->fpdR ^= ultrasonicMed->fpdR << 5;

  3095. -		inputSampleR += (((double)ultrasonicMed->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  3096. -		//end 32 bit stereo floating point dither

  3097. -

  3098.  		*out1 = (float) inputSampleL;

  3099.  		*out2 = (float) inputSampleR;

  3100.  

  3101. diff --git a/src/Vibrato/Vibrato.c b/src/Vibrato/Vibrato.c

  3102. index 8f63069..551d6d1 100644

  3103. --- a/src/Vibrato/Vibrato.c

  3104. +++ b/src/Vibrato/Vibrato.c

  3105. @@ -4,6 +4,7 @@

  3106.  #include <stdbool.h>

  3107.  #include <stdint.h>

  3108.  #include <stdlib.h>

  3109. +#include "../optimizations.h"

  3110.  

  3111.  #define VIBRATO_URI "https://hannesbraun.net/ns/lv2/airwindows/vibrato"

  3112.  

  3113. @@ -44,13 +45,11 @@ typedef struct {

  3114.  	double airFactorR;

  3115.  

  3116.  	bool flip;

  3117. -	uint32_t fpdL;

  3118. -	uint32_t fpdR;

  3119.  } Vibrato;

  3120.  

  3121.  static LV2_Handle instantiate(

  3122.  	const LV2_Descriptor* descriptor,

  3123. -	double rate,

  3124. +	realdouble rate,

  3125.  	const char* bundle_path,

  3126.  	const LV2_Feature* const* features)

  3127.  {

  3128. @@ -115,11 +114,6 @@ static void activate(LV2_Handle instance)

  3129.  	vibrato->airFactorR = 0.0;

  3130.  

  3131.  	vibrato->flip = false;

  3132. -

  3133. -	vibrato->fpdL = 1.0;

  3134. -	while (vibrato->fpdL < 16386) vibrato->fpdL = rand() * UINT32_MAX;

  3135. -	vibrato->fpdR = 1.0;

  3136. -	while (vibrato->fpdR < 16386) vibrato->fpdR = rand() * UINT32_MAX;

  3137.  }

  3138.  

  3139.  static void run(LV2_Handle instance, uint32_t sampleFrames)

  3140. @@ -141,8 +135,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  3141.  	while (sampleFrames-- > 0) {

  3142.  		double inputSampleL = *in1;

  3143.  		double inputSampleR = *in2;

  3144. -		if (fabs(inputSampleL) < 1.18e-23) inputSampleL = vibrato->fpdL * 1.18e-17;

  3145. -		if (fabs(inputSampleR) < 1.18e-23) inputSampleR = vibrato->fpdR * 1.18e-17;

  3146.  		double drySampleL = inputSampleL;

  3147.  		double drySampleR = inputSampleR;

  3148.  

  3149. @@ -219,20 +211,6 @@ static void run(LV2_Handle instance, uint32_t sampleFrames)

  3150.  		}

  3151.  		//Inv/Dry/Wet control

  3152.  

  3153. -		//begin 32 bit stereo floating point dither

  3154. -		int expon;

  3155. -		frexpf((float)inputSampleL, &expon);

  3156. -		vibrato->fpdL ^= vibrato->fpdL << 13;

  3157. -		vibrato->fpdL ^= vibrato->fpdL >> 17;

  3158. -		vibrato->fpdL ^= vibrato->fpdL << 5;

  3159. -		inputSampleL += (((double)vibrato->fpdL - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  3160. -		frexpf((float)inputSampleR, &expon);

  3161. -		vibrato->fpdR ^= vibrato->fpdR << 13;

  3162. -		vibrato->fpdR ^= vibrato->fpdR >> 17;

  3163. -		vibrato->fpdR ^= vibrato->fpdR << 5;

  3164. -		inputSampleR += (((double)vibrato->fpdR - (uint32_t)0x7fffffff) * 5.5e-36l * pow(2, expon + 62));

  3165. -		//end 32 bit stereo floating point dither

  3166. -

  3167.  		*out1 = (float) inputSampleL;

  3168.  		*out2 = (float) inputSampleR;

  3169.  

  3170. diff --git a/src/optimizations.h b/src/optimizations.h

  3171. new file mode 100644

  3172. index 0000000..1f6516b

  3173. --- /dev/null

  3174. +++ b/src/optimizations.h

  3175. @@ -0,0 +1,25 @@

  3176. +#include <math.h>

  3177. +#include <stdint.h>

  3178. +#include <stdlib.h>

  3179. +

  3180. +// force doubles into floats

  3181. +typedef double realdouble;

  3182. +#define double float

  3183. +#define asin asinf

  3184. +#define atan2 atan2f

  3185. +#define cos cosf

  3186. +#define exp expf

  3187. +#define exp2 exp2f

  3188. +#define fabs fabsf

  3189. +#define log2 log2f

  3190. +#define log10 log10f

  3191. +#define pow powf

  3192. +#define sin sinf

  3193. +#define tan tanf

  3194. +

  3195. +// ensure denormal things is unused

  3196. +#define fpdL error-if-fpdL-is-used

  3197. +#define fpdR error-if-fpdR-is-used

  3198. +

  3199. +// ensure 32bit dithering is unused

  3200. +#define frexpf error-if-frexpf-is-used