/* * File: Srsly.cpp * * Version: 1.0 * * Created: 9/3/19 * * Copyright: Copyright © 2019 Airwindows, All Rights Reserved * * Disclaimer: IMPORTANT: This Apple software is supplied to you by Apple Computer, Inc. ("Apple") in * consideration of your agreement to the following terms, and your use, installation, modification * or redistribution of this Apple software constitutes acceptance of these terms. If you do * not agree with these terms, please do not use, install, modify or redistribute this Apple * software. * * In consideration of your agreement to abide by the following terms, and subject to these terms, * Apple grants you a personal, non-exclusive license, under Apple's copyrights in this * original Apple software (the "Apple Software"), to use, reproduce, modify and redistribute the * Apple Software, with or without modifications, in source and/or binary forms; provided that if you * redistribute the Apple Software in its entirety and without modifications, you must retain this * notice and the following text and disclaimers in all such redistributions of the Apple Software. * Neither the name, trademarks, service marks or logos of Apple Computer, Inc. may be used to * endorse or promote products derived from the Apple Software without specific prior written * permission from Apple. Except as expressly stated in this notice, no other rights or * licenses, express or implied, are granted by Apple herein, including but not limited to any * patent rights that may be infringed by your derivative works or by other works in which the * Apple Software may be incorporated. * * The Apple Software is provided by Apple on an "AS IS" basis. APPLE MAKES NO WARRANTIES, EXPRESS OR * IMPLIED, INCLUDING WITHOUT LIMITATION THE IMPLIED WARRANTIES OF NON-INFRINGEMENT, MERCHANTABILITY * AND FITNESS FOR A PARTICULAR PURPOSE, REGARDING THE APPLE SOFTWARE OR ITS USE AND OPERATION ALONE * OR IN COMBINATION WITH YOUR PRODUCTS. * * IN NO EVENT SHALL APPLE BE LIABLE FOR ANY SPECIAL, INDIRECT, INCIDENTAL OR CONSEQUENTIAL * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS * OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) ARISING IN ANY WAY OUT OF THE USE, * REPRODUCTION, MODIFICATION AND/OR DISTRIBUTION OF THE APPLE SOFTWARE, HOWEVER CAUSED AND WHETHER * UNDER THEORY OF CONTRACT, TORT (INCLUDING NEGLIGENCE), STRICT LIABILITY OR OTHERWISE, EVEN * IF APPLE HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. * */ /*============================================================================= Srsly.cpp =============================================================================*/ #include "Srsly.h" //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ COMPONENT_ENTRY(Srsly) //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::Srsly //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Srsly::Srsly(AudioUnit component) : AUEffectBase(component) { CreateElements(); Globals()->UseIndexedParameters(kNumberOfParameters); SetParameter(kParam_One, kDefaultValue_ParamOne ); SetParameter(kParam_Two, kDefaultValue_ParamTwo ); SetParameter(kParam_Three, kDefaultValue_ParamThree ); SetParameter(kParam_Four, kDefaultValue_ParamFour ); SetParameter(kParam_Five, kDefaultValue_ParamFive ); #if AU_DEBUG_DISPATCHER mDebugDispatcher = new AUDebugDispatcher (this); #endif } //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::GetParameterValueStrings //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ComponentResult Srsly::GetParameterValueStrings(AudioUnitScope inScope, AudioUnitParameterID inParameterID, CFArrayRef * outStrings) { return kAudioUnitErr_InvalidProperty; } //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::GetParameterInfo //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ComponentResult Srsly::GetParameterInfo(AudioUnitScope inScope, AudioUnitParameterID inParameterID, AudioUnitParameterInfo &outParameterInfo ) { ComponentResult result = noErr; outParameterInfo.flags = kAudioUnitParameterFlag_IsWritable | kAudioUnitParameterFlag_IsReadable; if (inScope == kAudioUnitScope_Global) { switch(inParameterID) { case kParam_One: AUBase::FillInParameterName (outParameterInfo, kParameterOneName, false); outParameterInfo.unit = kAudioUnitParameterUnit_Generic; outParameterInfo.minValue = 0.0; outParameterInfo.maxValue = 1.0; outParameterInfo.defaultValue = kDefaultValue_ParamOne; break; case kParam_Two: AUBase::FillInParameterName (outParameterInfo, kParameterTwoName, false); outParameterInfo.unit = kAudioUnitParameterUnit_Generic; outParameterInfo.minValue = 0.0; outParameterInfo.maxValue = 1.0; outParameterInfo.defaultValue = kDefaultValue_ParamTwo; break; case kParam_Three: AUBase::FillInParameterName (outParameterInfo, kParameterThreeName, false); outParameterInfo.unit = kAudioUnitParameterUnit_Generic; outParameterInfo.minValue = 0.0; outParameterInfo.maxValue = 1.0; outParameterInfo.defaultValue = kDefaultValue_ParamThree; break; case kParam_Four: AUBase::FillInParameterName (outParameterInfo, kParameterFourName, false); outParameterInfo.unit = kAudioUnitParameterUnit_Generic; outParameterInfo.minValue = 0.0; outParameterInfo.maxValue = 1.0; outParameterInfo.defaultValue = kDefaultValue_ParamFour; break; case kParam_Five: AUBase::FillInParameterName (outParameterInfo, kParameterFiveName, false); outParameterInfo.unit = kAudioUnitParameterUnit_Generic; outParameterInfo.minValue = 0.0; outParameterInfo.maxValue = 1.0; outParameterInfo.defaultValue = kDefaultValue_ParamFive; break; default: result = kAudioUnitErr_InvalidParameter; break; } } else { result = kAudioUnitErr_InvalidParameter; } return result; } //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::GetPropertyInfo //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ComponentResult Srsly::GetPropertyInfo (AudioUnitPropertyID inID, AudioUnitScope inScope, AudioUnitElement inElement, UInt32 & outDataSize, Boolean & outWritable) { return AUEffectBase::GetPropertyInfo (inID, inScope, inElement, outDataSize, outWritable); } //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // state that plugin supports only stereo-in/stereo-out processing UInt32 Srsly::SupportedNumChannels(const AUChannelInfo ** outInfo) { if (outInfo != NULL) { static AUChannelInfo info; info.inChannels = 2; info.outChannels = 2; *outInfo = &info; } return 1; } //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::GetProperty //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ComponentResult Srsly::GetProperty( AudioUnitPropertyID inID, AudioUnitScope inScope, AudioUnitElement inElement, void * outData ) { return AUEffectBase::GetProperty (inID, inScope, inElement, outData); } // Srsly::Initialize //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ComponentResult Srsly::Initialize() { ComponentResult result = AUEffectBase::Initialize(); if (result == noErr) Reset(kAudioUnitScope_Global, 0); return result; } #pragma mark ____SrslyEffectKernel //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::SrslyKernel::Reset() //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ComponentResult Srsly::Reset(AudioUnitScope inScope, AudioUnitElement inElement) { for (int x = 0; x < 11; x++) { biquadM2[x] = 0.0; biquadM7[x] = 0.0; biquadM10[x] = 0.0; biquadL3[x] = 0.0; biquadL7[x] = 0.0; biquadR3[x] = 0.0; biquadR7[x] = 0.0; biquadS3[x] = 0.0; biquadS5[x] = 0.0; } fpd = 17; return noErr; } //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ // Srsly::ProcessBufferLists //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ OSStatus Srsly::ProcessBufferLists(AudioUnitRenderActionFlags & ioActionFlags, const AudioBufferList & inBuffer, AudioBufferList & outBuffer, UInt32 inFramesToProcess) { Float32 * inputL = (Float32*)(inBuffer.mBuffers[0].mData); Float32 * inputR = (Float32*)(inBuffer.mBuffers[1].mData); Float32 * outputL = (Float32*)(outBuffer.mBuffers[0].mData); Float32 * outputR = (Float32*)(outBuffer.mBuffers[1].mData); UInt32 nSampleFrames = inFramesToProcess; Float64 sampleRate = GetSampleRate(); if (sampleRate < 22000) sampleRate = 22000; //keep biquads in range long double tempSample; biquadM2[0] = 2000 / sampleRate; //up biquadM7[0] = 7000 / sampleRate; //down biquadM10[0] = 10000 / sampleRate; //down biquadL3[0] = 3000 / sampleRate; //up biquadL7[0] = 7000 / sampleRate; //way up biquadR3[0] = 3000 / sampleRate; //up biquadR7[0] = 7000 / sampleRate; //way up biquadS3[0] = 3000 / sampleRate; //up biquadS5[0] = 5000 / sampleRate; //way down Float64 focusM = 15.0-(GetParameter( kParam_One )*10.0); Float64 focusS = 20.0-(GetParameter( kParam_Two )*15.0); Float64 Q = GetParameter( kParam_Four )+0.25; //add Q control: from half to double intensity biquadM2[1] = focusM*0.25*Q; //Q, mid 2K boost is much broader biquadM7[1] = focusM*Q; //Q biquadM10[1] = focusM*Q; //Q biquadS3[1] = focusM*Q; //Q biquadS5[1] = focusM*Q; //Q biquadL3[1] = focusS*Q; //Q biquadL7[1] = focusS*Q; //Q biquadR3[1] = focusS*Q; //Q biquadR7[1] = focusS*Q; //Q double K = tan(M_PI * biquadM2[0]); double norm = 1.0 / (1.0 + K / biquadM2[1] + K * K); biquadM2[2] = K / biquadM2[1] * norm; biquadM2[4] = -biquadM2[2]; biquadM2[5] = 2.0 * (K * K - 1.0) * norm; biquadM2[6] = (1.0 - K / biquadM2[1] + K * K) * norm; K = tan(M_PI * biquadM7[0]); norm = 1.0 / (1.0 + K / biquadM7[1] + K * K); biquadM7[2] = K / biquadM7[1] * norm; biquadM7[4] = -biquadM7[2]; biquadM7[5] = 2.0 * (K * K - 1.0) * norm; biquadM7[6] = (1.0 - K / biquadM7[1] + K * K) * norm; K = tan(M_PI * biquadM10[0]); norm = 1.0 / (1.0 + K / biquadM10[1] + K * K); biquadM10[2] = K / biquadM10[1] * norm; biquadM10[4] = -biquadM10[2]; biquadM10[5] = 2.0 * (K * K - 1.0) * norm; biquadM10[6] = (1.0 - K / biquadM10[1] + K * K) * norm; K = tan(M_PI * biquadL3[0]); norm = 1.0 / (1.0 + K / biquadL3[1] + K * K); biquadL3[2] = K / biquadL3[1] * norm; biquadL3[4] = -biquadL3[2]; biquadL3[5] = 2.0 * (K * K - 1.0) * norm; biquadL3[6] = (1.0 - K / biquadL3[1] + K * K) * norm; K = tan(M_PI * biquadL7[0]); norm = 1.0 / (1.0 + K / biquadL7[1] + K * K); biquadL7[2] = K / biquadL7[1] * norm; biquadL7[4] = -biquadL7[2]; biquadL7[5] = 2.0 * (K * K - 1.0) * norm; biquadL7[6] = (1.0 - K / biquadL7[1] + K * K) * norm; K = tan(M_PI * biquadR3[0]); norm = 1.0 / (1.0 + K / biquadR3[1] + K * K); biquadR3[2] = K / biquadR3[1] * norm; biquadR3[4] = -biquadR3[2]; biquadR3[5] = 2.0 * (K * K - 1.0) * norm; biquadR3[6] = (1.0 - K / biquadR3[1] + K * K) * norm; K = tan(M_PI * biquadR7[0]); norm = 1.0 / (1.0 + K / biquadR7[1] + K * K); biquadR7[2] = K / biquadR7[1] * norm; biquadR7[4] = -biquadR7[2]; biquadR7[5] = 2.0 * (K * K - 1.0) * norm; biquadR7[6] = (1.0 - K / biquadR7[1] + K * K) * norm; K = tan(M_PI * biquadS3[0]); norm = 1.0 / (1.0 + K / biquadS3[1] + K * K); biquadS3[2] = K / biquadS3[1] * norm; biquadS3[4] = -biquadS3[2]; biquadS3[5] = 2.0 * (K * K - 1.0) * norm; biquadS3[6] = (1.0 - K / biquadS3[1] + K * K) * norm; K = tan(M_PI * biquadS5[0]); norm = 1.0 / (1.0 + K / biquadS5[1] + K * K); biquadS5[2] = K / biquadS5[1] * norm; biquadS5[4] = -biquadS5[2]; biquadS5[5] = 2.0 * (K * K - 1.0) * norm; biquadS5[6] = (1.0 - K / biquadS5[1] + K * K) * norm; Float64 depthM = pow(GetParameter( kParam_One ),2)*2.0; //proportion to mix in the filtered stuff Float64 depthS = pow(GetParameter( kParam_Two ),2)*2.0; //proportion to mix in the filtered stuff Float64 level = GetParameter( kParam_Three ); //output pad Float64 wet = GetParameter( kParam_Five ); //dry/wet //biquad contains these values: //[0] is frequency: 0.000001 to 0.499999 is near-zero to near-Nyquist //[1] is resonance, 0.7071 is Butterworth. Also can't be zero //[2] is a0 but you need distinct ones for additional biquad instances so it's here //[3] is a1 but you need distinct ones for additional biquad instances so it's here //[4] is a2 but you need distinct ones for additional biquad instances so it's here //[5] is b1 but you need distinct ones for additional biquad instances so it's here //[6] is b2 but you need distinct ones for additional biquad instances so it's here //[7] is LEFT stored delayed sample (freq and res are stored so you can move them sample by sample) //[8] is LEFT stored delayed sample (you have to include the coefficient making code if you do that) //[9] is RIGHT stored delayed sample (freq and res are stored so you can move them sample by sample) //[10] is RIGHT stored delayed sample (you have to include the coefficient making code if you do that) while (nSampleFrames-- > 0) { long double inputSampleL = *inputL; long double inputSampleR = *inputR; if (fabs(inputSampleL)<1.18e-37) inputSampleL = fpd * 1.18e-37; if (fabs(inputSampleR)<1.18e-37) inputSampleR = fpd * 1.18e-37; long double drySampleL = inputSampleL; long double drySampleR = inputSampleR; inputSampleL = sin(inputSampleL); inputSampleR = sin(inputSampleR); //encode Console5: good cleanness long double mid = inputSampleL + inputSampleR; long double rawmid = mid * 0.5; //we'll use this to isolate L&R a little long double side = inputSampleL - inputSampleR; long double boostside = side * depthS; //assign mid and side.Between these sections, you can do mid/side processing tempSample = (mid * biquadM2[2]) + biquadM2[7]; biquadM2[7] = (-tempSample * biquadM2[5]) + biquadM2[8]; biquadM2[8] = (mid * biquadM2[4]) - (tempSample * biquadM2[6]); long double M2Sample = tempSample; //like mono AU, 7 and 8 store L channel tempSample = (mid * biquadM7[2]) + biquadM7[7]; biquadM7[7] = (-tempSample * biquadM7[5]) + biquadM7[8]; biquadM7[8] = (mid * biquadM7[4]) - (tempSample * biquadM7[6]); long double M7Sample = -tempSample*2.0; //like mono AU, 7 and 8 store L channel tempSample = (mid * biquadM10[2]) + biquadM10[7]; biquadM10[7] = (-tempSample * biquadM10[5]) + biquadM10[8]; biquadM10[8] = (mid * biquadM10[4]) - (tempSample * biquadM10[6]); long double M10Sample = -tempSample*2.0; //like mono AU, 7 and 8 store L channel //mid tempSample = (side * biquadS3[2]) + biquadS3[7]; biquadS3[7] = (-tempSample * biquadS3[5]) + biquadS3[8]; biquadS3[8] = (side * biquadS3[4]) - (tempSample * biquadS3[6]); long double S3Sample = tempSample*2.0; //like mono AU, 7 and 8 store L channel tempSample = (side * biquadS5[2]) + biquadS5[7]; biquadS5[7] = (-tempSample * biquadS5[5]) + biquadS5[8]; biquadS5[8] = (side * biquadS5[4]) - (tempSample * biquadS5[6]); long double S5Sample = -tempSample*5.0; //like mono AU, 7 and 8 store L channel mid = (M2Sample + M7Sample + M10Sample)*depthM; side = (S3Sample + S5Sample + boostside)*depthS; long double msOutSampleL = (mid+side)/2.0; long double msOutSampleR = (mid-side)/2.0; //unassign mid and side long double isoSampleL = inputSampleL-rawmid; long double isoSampleR = inputSampleR-rawmid; //trying to isolate L and R a little tempSample = (isoSampleL * biquadL3[2]) + biquadL3[7]; biquadL3[7] = (-tempSample * biquadL3[5]) + biquadL3[8]; biquadL3[8] = (isoSampleL * biquadL3[4]) - (tempSample * biquadL3[6]); long double L3Sample = tempSample; //like mono AU, 7 and 8 store L channel tempSample = (isoSampleR * biquadR3[2]) + biquadR3[9]; biquadR3[9] = (-tempSample * biquadR3[5]) + biquadR3[10]; biquadR3[10] = (isoSampleR * biquadR3[4]) - (tempSample * biquadR3[6]); long double R3Sample = tempSample; //note: 9 and 10 store the R channel tempSample = (isoSampleL * biquadL7[2]) + biquadL7[7]; biquadL7[7] = (-tempSample * biquadL7[5]) + biquadL7[8]; biquadL7[8] = (isoSampleL * biquadL7[4]) - (tempSample * biquadL7[6]); long double L7Sample = tempSample*3.0; //like mono AU, 7 and 8 store L channel tempSample = (isoSampleR * biquadR7[2]) + biquadR7[9]; biquadR7[9] = (-tempSample * biquadR7[5]) + biquadR7[10]; biquadR7[10] = (isoSampleR * biquadR7[4]) - (tempSample * biquadR7[6]); long double R7Sample = tempSample*3.0; //note: 9 and 10 store the R channel long double processingL = msOutSampleL + ((L3Sample + L7Sample)*depthS); long double processingR = msOutSampleR + ((R3Sample + R7Sample)*depthS); //done with making filters, now we apply them inputSampleL += processingL; inputSampleR += processingR; if (level < 1.0) { inputSampleL *= level; inputSampleR *= level; } if (inputSampleL > 1.0) inputSampleL = 1.0; if (inputSampleL < -1.0) inputSampleL = -1.0; if (inputSampleR > 1.0) inputSampleR = 1.0; if (inputSampleR < -1.0) inputSampleR = -1.0; //without this, you can get a NaN condition where it spits out DC offset at full blast! inputSampleL = asin(inputSampleL); inputSampleR = asin(inputSampleR); //amplitude aspect if (wet < 1.0) { inputSampleL = (inputSampleL * wet)+(drySampleL * (1.0-wet)); inputSampleR = (inputSampleR * wet)+(drySampleR * (1.0-wet)); } //begin 32 bit stereo floating point dither int expon; frexpf((float)inputSampleL, &expon); fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; inputSampleL += ((double(fpd)-uint32_t(0x7fffffff)) * 5.5e-36l * pow(2,expon+62)); frexpf((float)inputSampleR, &expon); fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; inputSampleR += ((double(fpd)-uint32_t(0x7fffffff)) * 5.5e-36l * pow(2,expon+62)); //end 32 bit stereo floating point dither *outputL = inputSampleL; *outputR = inputSampleR; //direct stereo out inputL += 1; inputR += 1; outputL += 1; outputR += 1; } return noErr; }