Floating Point Dither For All

1 files changed, 10 insertions, 29 deletions

diff --git a/plugins/MacAU/Wider/Wider.cpp b/plugins/MacAU/Wider/Wider.cpp
index 050f43a..3063f97 100755
--- a/plugins/MacAU/Wider/Wider.cpp
+++ b/plugins/MacAU/Wider/Wider.cpp
@@ -189,11 +189,8 @@ ComponentResult		Wider::Reset(AudioUnitScope inScope, AudioUnitElement inElement
 {
 	for(int fcount = 0; fcount < 4098; fcount++) {p[fcount] = 0.0;}
 	count = 0;
-	flip = false;
-	fpNShapeAL = 0.0;
-	fpNShapeBL = 0.0;
-	fpNShapeAR = 0.0;
-	fpNShapeBR = 0.0;
+	fpNShapeL = 0.0;
+	fpNShapeR = 0.0;
 	return noErr;
 }
 
@@ -215,9 +212,6 @@ OSStatus		Wider::ProcessBufferLists(AudioUnitRenderActionFlags & ioActionFlags,
 	Float64 overallscale = 1.0;
 	overallscale /= 44100.0;
 	overallscale *= GetSampleRate();
-	Float32 fpTemp;
-	Float64 fpOld = 0.618033988749894848204586; //golden ratio!
-	Float64 fpNew = 1.0 - fpOld;
 	long double inputSampleL;
 	long double inputSampleR;
 	Float64 drySampleL;
@@ -338,27 +332,14 @@ OSStatus		Wider::ProcessBufferLists(AudioUnitRenderActionFlags & ioActionFlags,
 		inputSampleL = (drySampleL * dry) + ((mid+side) * wet);
 		inputSampleR = (drySampleR * dry) + ((mid-side) * wet);
 		
-		//noise shaping to 32-bit floating point
-		if (flip) {
-			fpTemp = inputSampleL;
-			fpNShapeAL = (fpNShapeAL*fpOld)+((inputSampleL-fpTemp)*fpNew);
-			inputSampleL += fpNShapeAL;
-			
-			fpTemp = inputSampleR;
-			fpNShapeAR = (fpNShapeAR*fpOld)+((inputSampleR-fpTemp)*fpNew);
-			inputSampleR += fpNShapeAR;
-		}
-		else {
-			fpTemp = inputSampleL;
-			fpNShapeBL = (fpNShapeBL*fpOld)+((inputSampleL-fpTemp)*fpNew);
-			inputSampleL += fpNShapeBL;
-			
-			fpTemp = inputSampleR;
-			fpNShapeBR = (fpNShapeBR*fpOld)+((inputSampleR-fpTemp)*fpNew);
-			inputSampleR += fpNShapeBR;
-		}
-		flip = !flip;
-		//end noise shaping on 32 bit output
+		//stereo 32 bit dither, made small and tidy.
+		int expon; frexpf((Float32)inputSampleL, &expon);
+		long double dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62);
+		inputSampleL += (dither-fpNShapeL); fpNShapeL = dither;
+		frexpf((Float32)inputSampleR, &expon);
+		dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62);
+		inputSampleR += (dither-fpNShapeR); fpNShapeR = dither;
+		//end 32 bit dither
 
 		*outputL = inputSampleL;
 		*outputR = inputSampleR;