From 966f2d253cd2ee6ce140ad68095a20a9d2b63052 Mon Sep 17 00:00:00 2001 From: Chris Johnson Date: Sun, 27 Jan 2019 21:13:54 -0500 Subject: Floating Point Dither For All --- .../MacVST/uLawDecode/source/uLawDecodeProc.cpp | 54 ++++++++-------------- .../christopherjohnson.pbxuser | 4 +- .../christopherjohnson.perspectivev3 | 17 ++++--- 3 files changed, 28 insertions(+), 47 deletions(-) (limited to 'plugins/MacVST/uLawDecode') diff --git a/plugins/MacVST/uLawDecode/source/uLawDecodeProc.cpp b/plugins/MacVST/uLawDecode/source/uLawDecodeProc.cpp index a78b3c2..1835623 100755 --- a/plugins/MacVST/uLawDecode/source/uLawDecodeProc.cpp +++ b/plugins/MacVST/uLawDecode/source/uLawDecodeProc.cpp @@ -88,18 +88,14 @@ void uLawDecode::processReplacing(float **inputs, float **outputs, VstInt32 samp inputSampleR = (inputSampleR * wet) + (drySampleR * dry); } - //noise shaping to 32-bit floating point - float fpTemp = inputSampleL; - fpNShapeL += (inputSampleL-fpTemp); - inputSampleL += fpNShapeL; - //if this confuses you look at the wordlength for fpTemp :) - fpTemp = inputSampleR; - fpNShapeR += (inputSampleR-fpTemp); - inputSampleR += fpNShapeR; - //for deeper space and warmth, we try a non-oscillating noise shaping - //that is kind of ruthless: it will forever retain the rounding errors - //except we'll dial it back a hair at the end of every buffer processed - //end noise shaping on 32 bit output + //stereo 32 bit dither, made small and tidy. + int expon; frexpf((float)inputSampleL, &expon); + long double dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62); + inputSampleL += (dither-fpNShapeL); fpNShapeL = dither; + frexpf((float)inputSampleR, &expon); + dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62); + inputSampleR += (dither-fpNShapeR); fpNShapeR = dither; + //end 32 bit dither *out1 = inputSampleL; *out2 = inputSampleR; @@ -109,12 +105,6 @@ void uLawDecode::processReplacing(float **inputs, float **outputs, VstInt32 samp *out1++; *out2++; } - fpNShapeL *= 0.999999; - fpNShapeR *= 0.999999; - //we will just delicately dial back the FP noise shaping, not even every sample - //this is a good place to put subtle 'no runaway' calculations, though bear in mind - //that it will be called more often when you use shorter sample buffers in the DAW. - //So, very low latency operation will call these calculations more often. } void uLawDecode::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) @@ -198,18 +188,16 @@ void uLawDecode::processDoubleReplacing(double **inputs, double **outputs, VstIn inputSampleR = (inputSampleR * wet) + (drySampleR * dry); } - //noise shaping to 64-bit floating point - double fpTemp = inputSampleL; - fpNShapeL += (inputSampleL-fpTemp); - inputSampleL += fpNShapeL; - //if this confuses you look at the wordlength for fpTemp :) - fpTemp = inputSampleR; - fpNShapeR += (inputSampleR-fpTemp); - inputSampleR += fpNShapeR; - //for deeper space and warmth, we try a non-oscillating noise shaping - //that is kind of ruthless: it will forever retain the rounding errors - //except we'll dial it back a hair at the end of every buffer processed - //end noise shaping on 64 bit output + //stereo 64 bit dither, made small and tidy. + int expon; frexp((double)inputSampleL, &expon); + long double dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62); + dither /= 536870912.0; //needs this to scale to 64 bit zone + inputSampleL += (dither-fpNShapeL); fpNShapeL = dither; + frexp((double)inputSampleR, &expon); + dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62); + dither /= 536870912.0; //needs this to scale to 64 bit zone + inputSampleR += (dither-fpNShapeR); fpNShapeR = dither; + //end 64 bit dither *out1 = inputSampleL; *out2 = inputSampleR; @@ -219,10 +207,4 @@ void uLawDecode::processDoubleReplacing(double **inputs, double **outputs, VstIn *out1++; *out2++; } - fpNShapeL *= 0.999999; - fpNShapeR *= 0.999999; - //we will just delicately dial back the FP noise shaping, not even every sample - //this is a good place to put subtle 'no runaway' calculations, though bear in mind - //that it will be called more often when you use shorter sample buffers in the DAW. - //So, very low latency operation will call these calculations more often. } diff --git a/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.pbxuser b/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.pbxuser index a221ade..78e474f 100755 --- a/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.pbxuser +++ b/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.pbxuser @@ -75,8 +75,8 @@ PBXFileDataSource_Warnings_ColumnID, ); }; - PBXPerProjectTemplateStateSaveDate = 561241927; - PBXWorkspaceStateSaveDate = 561241927; + PBXPerProjectTemplateStateSaveDate = 569773408; + PBXWorkspaceStateSaveDate = 569773408; }; sourceControlManager = 8B02375E1D42B1C400E1E8C8 /* Source Control */; userBuildSettings = { diff --git a/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.perspectivev3 b/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.perspectivev3 index 11e6474..45637a5 100755 --- a/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.perspectivev3 +++ b/plugins/MacVST/uLawDecode/uLawDecode.xcodeproj/christopherjohnson.perspectivev3 @@ -323,7 +323,7 @@ 185 RubberWindowFrame - 24 356 810 487 0 0 1440 878 + 426 311 810 487 0 0 1440 878 Module PBXSmartGroupTreeModule @@ -362,7 +362,7 @@ Frame {{0, 0}, {603, 0}} RubberWindowFrame - 24 356 810 487 0 0 1440 878 + 426 311 810 487 0 0 1440 878 Module PBXNavigatorGroup @@ -387,7 +387,7 @@ Frame {{10, 27}, {603, 414}} RubberWindowFrame - 24 356 810 487 0 0 1440 878 + 426 311 810 487 0 0 1440 878 Module XCDetailModule @@ -469,11 +469,11 @@ TableOfContents - 8B2721E12173DF5E00396442 + 8B79426F21F60D67006E9731 1CA23ED40692098700951B8B - 8B2721E22173DF5E00396442 + 8B79427021F60D67006E9731 8B0237581D42B1C400E1E8C8 - 8B2721E32173DF5E00396442 + 8B79427121F60D67006E9731 1CA23EDF0692099D00951B8B 1CA23EE00692099D00951B8B 1CA23EE10692099D00951B8B @@ -646,7 +646,7 @@ StatusbarIsVisible TimeStamp - 561241950.77135801 + 569773415.70611095 ToolbarConfigUserDefaultsMinorVersion 2 ToolbarDisplayMode @@ -663,11 +663,10 @@ 5 WindowOrderList - 8B2721E42173DF5E00396442 /Users/christopherjohnson/Desktop/MacVST/uLawDecode/uLawDecode.xcodeproj WindowString - 24 356 810 487 0 0 1440 878 + 426 311 810 487 0 0 1440 878 WindowToolsV3 -- cgit v1.2.3