| | | | | | sampler.cpp // // SYNTOPIA. See http://Syntopia.sourceforge.net for details and documentation. // // Author of this file: Mikael Hvidtfeldt Christensen (mikaelc@users.sourceforge.net) // #include "Sampler.h" #include "math.h" #include <iostream> #include "SynthVoice.h" // FMOSC ---------------------------------------------------------------------------------------------------------- namespace SynthCore { void Sampler::setMidiNote(int note) { sample = mySampleMap->getSamples(note); } Module * Sampler::newInstance(SynthVoice * S1) { return new Sampler(S1); } Module * Sampler::sharedClone() { Sampler * S = new Sampler(hostVoice); S->ModID = ModID; S->firstVoice = false; S->freq = freq; S->interpolationType = interpolationType; if (S->mySampleMap != 0) delete S->mySampleMap; S->mySampleMap = mySampleMap; return S; } Sampler::Sampler(SynthVoice * S1) { freqs = new FreqTable(); hostVoice = S1; Output1 = new Output(this); freq = 440.0; pos = 0; interpolationType = spline5Point; factoryWaveform = sawtooth; mySampleMap = new SampleMap(); setSawtooth(1024); // Init with triangle waveform. sample = mySampleMap->getSamples(40); MAX = MIN = 0.0f; // For debugging } Sampler::~Sampler() { delete Output1; delete mySampleMap; delete (freqs); } std::string Sampler::getXML(int indent ) { std::string S = Utils::space(indent) + "<Module ID=\""+Utils::intToString((int)this)+"\" >"+ Utils::newline(); S += Utils::space(indent+2) + "<TypeID>" + typeid(*this).name() + "</TypeID>" + Utils::newline(); S += saveParameter("interpolationType",interpolationType, indent+2); S += saveParameter("factoryWaveform",factoryWaveform, indent+2); S += Output1->getXML("Output1", indent+2); S += saveID(indent+2); S += Utils::space(indent) + "</Module>" + Utils::newline(); return S; } void Sampler::loadXML(XMLNode * n, bool firstPass) { cout << "Parsing Sampler XML" << endl; string TypeID = n->getTagText("TypeID") ; if (typeid(*this).name() != TypeID) { throw parseError("Trying to put "+TypeID+" into "+typeid(*this).name()); } if (firstPass) { loadParameter("interpolationType",interpolationType, n); loadParameter("factoryWaveform",factoryWaveform, n); Output1->loadXML("Output1",n); } }; void Sampler::setTriangle(int length) { factoryWaveform = triangle; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::triangle(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setSawtooth(int length) { factoryWaveform = sawtooth; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::sawtooth(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setSquare(int length) { factoryWaveform = square; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::square(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setSine(int length) { factoryWaveform = sine; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::sine(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setWave4(int length) { factoryWaveform = wave4; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::wave4(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setWave1(int length) { factoryWaveform = wave1; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::wave1(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setWave2(int length) { factoryWaveform = wave2; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::wave2(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::setWave3(int length) { factoryWaveform = wave4; mySampleMap->clean(); for (int n=21; n<109; n++) { // find no harmonics; int maxPartials = freqs->getMaxPartialsFromMidi(n); mySampleMap->map[n] = Fourier::wave3(length,maxPartials,6); } sample = mySampleMap->getSamples(40); } void Sampler::resetPhase() { pos = 0.0f; } void Sampler::newMidiEvent() { freq = hostVoice->getMidiFreq(); setMidiNote(hostVoice->getMidiNote()); } void Sampler::update() { // Add phase to phase position: pos += freq / (48000.0f/1024.0f); // Check wrap-around. if (pos > 1024.0f) { pos -= 1024.0f; } //fractpart = (pos - (float) n); float half = 0.5; float fractpart ; float posi = pos; int n; // fast float->int conversion _asm{ fld posi fsub half fistp n ; fist / fistp does NOT round value stored on stack, must reload it. fild n fsubr posi fstp fractpart } switch (interpolationType) { case cubicHermite: // now do the interpolation Output1->write( ((( ( ( ( ( 3.0f* ( sample[n+1] - sample[n+2])) - sample[n] ) + sample[n+3] ) * 0.5f * fractpart) + ((2.0f* sample[n+2]) + sample[n]) - (((5.0f* sample[n+1])+ sample[n+3])*0.5f) ) * fractpart) + ((sample[n+2]-sample[n])*0.5f))*fractpart + sample[n+1] ); break; case spline5Point: { float p0 = sample[n]; float p1 = sample[n+1]; float p2 = sample[n+2]; float p3 = sample[n+3]; float p4 = sample[n+4]; float p5 = sample[n+5]; Output1->write( p2 + 0.04166666666*fractpart*((p3-p1)*16.0+(p0-p4)*2.0 + fractpart *((p3+p1)*16.0-p0-p2*30.0- p4 + fractpart *(p3*66.0-p2*70.0-p4*33.0+p1*39.0+ p5*7.0- p0*9.0 + fractpart *( p2*126.0-p3*124.0+p4*61.0-p1*64.0- p5*12.0+p0*13.0 + fractpart *((p3-p2)*50.0+(p1-p4)*25.0+(p5-p0)*5.0))))) ); } break; case none: Output1->write( sample[n] ); break; case linear: Output1->write( fractpart*(sample[n+1]+sample[n])+sample[n] ); break; case pureSine: Output1->write( (float) cos( 2*(pos/1024.0f)*3.14159276) ); break; default: break; } if (Output1->read() < MIN) { MIN = Output1->read(); } if (Output1->read() > MAX) { MAX = Output1->read(); std::cout << "--------------------MAX FOUND: " << MAX << " MIN = " << MIN << " AT n = " << n << " | SAMPLE[n] = " << sample[n] << " pos = " << pos << " posi = " << posi << " fractpart = " << fractpart << std::endl; } } }; // end of namespace: SynthCore Syntopia Project. Visit the web page, or the SourceForge page. Docs made by Doxygen. Email: Mikael Christensen