Pd-cvs August 2004

pd-cvs@lists.iem.at

10 participants
171 discussions

externals/sc4pd config-pd-linux.txt,1.5,1.6 make-files.txt,1.18,1.19
by Tim Blechmann 07 Aug '04

07 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv24095 Modified Files: config-pd-linux.txt make-files.txt Log Message: sc4pd_dsp class and other changes Index: config-pd-linux.txt =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/config-pd-linux.txt,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** config-pd-linux.txt 3 Aug 2004 17:51:50 -0000 1.5 --- config-pd-linux.txt 7 Aug 2004 14:01:30 -0000 1.6 *************** *** 2,6 **** # your c++ compiler (if not g++) ! # CXX=icc --- 2,6 ---- # your c++ compiler (if not g++) ! CXX=icc Index: make-files.txt =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/make-files.txt,v retrieving revision 1.18 retrieving revision 1.19 diff -C2 -d -r1.18 -r1.19 *** make-files.txt 6 Aug 2004 08:34:41 -0000 1.18 --- make-files.txt 7 Aug 2004 14:01:30 -0000 1.19 *************** *** 13,15 **** Lag2.cpp Lag3.cpp LinExp.cpp DelayUnit.cpp DelayN.cpp DelayL.cpp \ DelayC.cpp CombN.cpp CombL.cpp CombC.cpp AllpassN.cpp AllpassL.cpp \ ! AllpassC.cpp PitchShift.cpp --- 13,15 ---- Lag2.cpp Lag3.cpp LinExp.cpp DelayUnit.cpp DelayN.cpp DelayL.cpp \ DelayC.cpp CombN.cpp CombL.cpp CombC.cpp AllpassN.cpp AllpassL.cpp \ ! AllpassC.cpp PitchShift.cpp Resonz.cpp OnePole.cpp

1 0

externals/sc4pd make-files.txt,1.17,1.18
by Tim Blechmann 06 Aug '04

06 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv16321 Modified Files: make-files.txt Log Message: checkin of PitchShift~ Index: make-files.txt =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/make-files.txt,v retrieving revision 1.17 retrieving revision 1.18 diff -C2 -d -r1.17 -r1.18 *** make-files.txt 4 Aug 2004 18:23:14 -0000 1.17 --- make-files.txt 6 Aug 2004 08:34:41 -0000 1.18 *************** *** 13,15 **** Lag2.cpp Lag3.cpp LinExp.cpp DelayUnit.cpp DelayN.cpp DelayL.cpp \ DelayC.cpp CombN.cpp CombL.cpp CombC.cpp AllpassN.cpp AllpassL.cpp \ ! AllpassC.cpp --- 13,15 ---- Lag2.cpp Lag3.cpp LinExp.cpp DelayUnit.cpp DelayN.cpp DelayL.cpp \ DelayC.cpp CombN.cpp CombL.cpp CombC.cpp AllpassN.cpp AllpassL.cpp \ ! AllpassC.cpp PitchShift.cpp

1 0

externals/sc4pd/source PitchShift.cpp,NONE,1.1 main.cpp,1.19,1.20
by Tim Blechmann 06 Aug '04

06 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/source In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv16321/source Modified Files: main.cpp Added Files: PitchShift.cpp Log Message: checkin of PitchShift~ Index: main.cpp =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/source/main.cpp,v retrieving revision 1.19 retrieving revision 1.20 diff -C2 -d -r1.19 -r1.20 *** main.cpp 4 Aug 2004 18:23:15 -0000 1.19 --- main.cpp 6 Aug 2004 08:34:42 -0000 1.20 *************** *** 63,67 **** " DelayL~, DelayC~, CombN~, CombL~, CombC~, AllpassN~, " "AllpassL~,\n" ! " AllpassC~" "\n" ); --- 63,67 ---- " DelayL~, DelayC~, CombN~, CombL~, CombC~, AllpassN~, " "AllpassL~,\n" ! " AllpassC~, PitchShift~" "\n" ); --- NEW FILE: PitchShift.cpp --- /* sc4pd PitchShift~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: */ #include "sc4pd.hpp" /* ------------------------ PitchShift~ -----------------------------*/ class PitchShift_ar :public flext_dsp { FLEXT_HEADER(PitchShift_ar,flext_dsp); public: PitchShift_ar(int argc,t_atom * argv); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out); void m_set_windowsize(float f) { m_windowsize = f; } void m_set_pitchratio (float f) { m_pitchratio = f; } void m_set_pitchdispersion (float f) { m_pitchdispersion = f; } void m_set_timedispersion (float f) { m_timedispersion = f; } private: float m_windowsize, m_pitchratio,m_pitchdispersion,m_timedispersion; RGen rgen; float *m_dlybuf; float m_dsamp1, m_dsamp1_slope, m_ramp1, m_ramp1_slope; float m_dsamp2, m_dsamp2_slope, m_ramp2, m_ramp2_slope; float m_dsamp3, m_dsamp3_slope, m_ramp3, m_ramp3_slope; float m_dsamp4, m_dsamp4_slope, m_ramp4, m_ramp4_slope; float m_fdelaylen, m_slope; long m_iwrphase, m_idelaylen, m_mask; long m_counter, m_stage, m_numoutput, m_framesize; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_set_windowsize); FLEXT_CALLBACK_F(m_set_pitchratio); FLEXT_CALLBACK_F(m_set_pitchdispersion); FLEXT_CALLBACK_F(m_set_timedispersion); }; FLEXT_LIB_DSP_V("PitchShift~",PitchShift_ar); PitchShift_ar::PitchShift_ar(int argc,t_atom * argv) { FLEXT_ADDMETHOD_(0,"windowSize",m_set_windowsize); FLEXT_ADDMETHOD_(0,"pitchRatio",m_set_pitchratio); FLEXT_ADDMETHOD_(0,"pitchDispersion",m_set_pitchdispersion); FLEXT_ADDMETHOD_(0,"timeDispersion",m_set_timedispersion); AtomList Args(argc,argv); if (Args.Count() != 4) { post("4 arguments needed"); return; } m_windowsize = sc_getfloatarg(Args,0); m_pitchratio = sc_getfloatarg(Args,1); m_pitchdispersion = sc_getfloatarg(Args,2); m_timedispersion = sc_getfloatarg(Args,3); rgen.init(timeseed()); AddOutSignal(); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); } void PitchShift_ar::m_dsp(int n, t_sample *const *in, t_sample *const *out) { /* initialization from PitchShift_Ctor(PitchShift *unit) */ long delaybufsize; float /* *out, *in,*/ *dlybuf; float winsize, pchratio; float fdelaylen, slope; long framesize, last; //out = ZOUT(0); //in = ZIN(0); pchratio = m_pitchratio; winsize = m_windowsize; delaybufsize = (long)ceil(winsize * SAMPLERATE * 3.f + 3.f); fdelaylen = delaybufsize - 3; delaybufsize = delaybufsize + BUFLENGTH; delaybufsize = NEXTPOWEROFTWO(delaybufsize); // round up to next power of two dlybuf = new float[delaybufsize]; //(float*)RTAlloc(unit->mWorld, delaybufsize * sizeof(float)); *dlybuf = 0; m_dlybuf = dlybuf; m_idelaylen = delaybufsize; m_fdelaylen = fdelaylen; m_iwrphase = 0; m_numoutput = 0; m_mask = last = (delaybufsize - 1); m_framesize = framesize = ((long)(winsize * SAMPLERATE) + 2) & ~3; m_slope = slope = 2.f / framesize; m_stage = 3; m_counter = framesize >> 2; m_ramp1 = 0.5; m_ramp2 = 1.0; m_ramp3 = 0.5; m_ramp4 = 0.0; m_ramp1_slope = -slope; m_ramp2_slope = -slope; m_ramp3_slope = slope; m_ramp4_slope = slope; dlybuf[last ] = 0.f; // put a few zeroes where we start the read heads dlybuf[last-1] = 0.f; dlybuf[last-2] = 0.f; m_numoutput = 0; // start all read heads 2 samples behind the write head m_dsamp1 = m_dsamp2 = m_dsamp3 = m_dsamp4 = 2.f; // pch ratio is initially zero for the read heads m_dsamp1_slope = m_dsamp2_slope = m_dsamp3_slope = m_dsamp4_slope = 1.f; } void PitchShift_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { float *nout, *nin, *dlybuf; float disppchratio, pchratio, pchratio1, value; float dsamp1, dsamp1_slope, ramp1, ramp1_slope; float dsamp2, dsamp2_slope, ramp2, ramp2_slope; float dsamp3, dsamp3_slope, ramp3, ramp3_slope; float dsamp4, dsamp4_slope, ramp4, ramp4_slope; float fdelaylen, d1, d2, frac, slope, samp_slope, startpos, winsize, pchdisp, timedisp; long remain, nsmps, idelaylen, irdphase, irdphaseb, iwrphase; long mask, idsamp; long counter, stage, framesize, numoutput; RGET; nout = *out; nin = *in; pchratio = m_pitchratio; winsize = m_windowsize; pchdisp = m_pitchdispersion; timedisp = m_timedispersion; timedisp = sc_clip(timedisp, 0.f, winsize) * SAMPLERATE; dlybuf = m_dlybuf; fdelaylen = m_fdelaylen; idelaylen = m_idelaylen; iwrphase = m_iwrphase; numoutput = m_numoutput; counter = m_counter; stage = m_stage; mask = m_mask; framesize = m_framesize; dsamp1 = m_dsamp1; dsamp2 = m_dsamp2; dsamp3 = m_dsamp3; dsamp4 = m_dsamp4; dsamp1_slope = m_dsamp1_slope; dsamp2_slope = m_dsamp2_slope; dsamp3_slope = m_dsamp3_slope; dsamp4_slope = m_dsamp4_slope; ramp1 = m_ramp1; ramp2 = m_ramp2; ramp3 = m_ramp3; ramp4 = m_ramp4; ramp1_slope = m_ramp1_slope; ramp2_slope = m_ramp2_slope; ramp3_slope = m_ramp3_slope; ramp4_slope = m_ramp4_slope; slope = m_slope; remain = n; while (remain) { if (counter <= 0) { counter = framesize >> 2; m_stage = stage = (stage + 1) & 3; disppchratio = pchratio; if (pchdisp != 0.f) { disppchratio += (pchdisp * frand2(s1,s2,s3)); } disppchratio = sc_clip(disppchratio, 0.f, 4.f); pchratio1 = disppchratio - 1.f; samp_slope = -pchratio1; startpos = pchratio1 < 0.f ? 2.f : framesize * pchratio1 + 2.f; startpos += (timedisp * frand(s1,s2,s3)); switch(stage) { case 0 : m_dsamp1_slope = dsamp1_slope = samp_slope; dsamp1 = startpos; ramp1 = 0.0; m_ramp1_slope = ramp1_slope = slope; m_ramp3_slope = ramp3_slope = -slope; break; case 1 : m_dsamp2_slope = dsamp2_slope = samp_slope; dsamp2 = startpos; ramp2 = 0.0; m_ramp2_slope = ramp2_slope = slope; m_ramp4_slope = ramp4_slope = -slope; break; case 2 : m_dsamp3_slope = dsamp3_slope = samp_slope; dsamp3 = startpos; ramp3 = 0.0; m_ramp3_slope = ramp3_slope = slope; m_ramp1_slope = ramp1_slope = -slope; break; case 3 : m_dsamp4_slope = dsamp4_slope = samp_slope; dsamp4 = startpos; ramp4 = 0.0; m_ramp2_slope = ramp2_slope = -slope; m_ramp4_slope = ramp4_slope = slope; break; } } nsmps = sc_min(remain, counter); remain -= nsmps; counter -= nsmps; while (nsmps--) { numoutput++; iwrphase = (iwrphase + 1) & mask; dsamp1 += dsamp1_slope; idsamp = (long)dsamp1; frac = dsamp1 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; if (numoutput < idelaylen) { if (irdphase > iwrphase) { value = 0.f; } else if (irdphaseb > iwrphase) { d1 = dlybuf[irdphase]; value = (d1 - frac * d1) * ramp1; } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value = (d1 + frac * (d2 - d1)) * ramp1; } } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value = (d1 + frac * (d2 - d1)) * ramp1; } ramp1 += ramp1_slope; dsamp2 += dsamp2_slope; idsamp = (long)dsamp2; frac = dsamp2 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; if (numoutput < idelaylen) { if (irdphase > iwrphase) { //value += 0.f; } else if (irdphaseb > iwrphase) { d1 = dlybuf[irdphase]; value += (d1 - frac * d1) * ramp2; } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp2; } } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp2; } ramp2 += ramp2_slope; dsamp3 += dsamp3_slope; idsamp = (long)dsamp3; frac = dsamp3 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; if (numoutput < idelaylen) { if (irdphase > iwrphase) { //value += 0.f; } else if (irdphaseb > iwrphase) { d1 = dlybuf[irdphase]; value += (d1 - frac * d1) * ramp3; } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp3; } } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp3; } ramp3 += ramp3_slope; dsamp4 += dsamp4_slope; idsamp = (long)dsamp4; frac = dsamp4 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; if (numoutput < idelaylen) { if (irdphase > iwrphase) { //value += 0.f; } else if (irdphaseb > iwrphase) { d1 = dlybuf[irdphase]; value += (d1 - frac * d1) * ramp4; } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp4; } } else { d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp4; } ramp4 += ramp4_slope; dlybuf[iwrphase] = ZXP(nin); ZXP(nout) = value *= 0.5; } } m_counter = counter; m_stage = stage; m_mask = mask; m_dsamp1 = dsamp1; m_dsamp2 = dsamp2; m_dsamp3 = dsamp3; m_dsamp4 = dsamp4; m_ramp1 = ramp1; m_ramp2 = ramp2; m_ramp3 = ramp3; m_ramp4 = ramp4; m_numoutput = numoutput; m_iwrphase = iwrphase; if (numoutput >= idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); } RPUT; } void PitchShift_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { float *nout, *nin, *dlybuf; float disppchratio, pchratio, pchratio1, value; float dsamp1, dsamp1_slope, ramp1, ramp1_slope; float dsamp2, dsamp2_slope, ramp2, ramp2_slope; float dsamp3, dsamp3_slope, ramp3, ramp3_slope; float dsamp4, dsamp4_slope, ramp4, ramp4_slope; float fdelaylen, d1, d2, frac, slope, samp_slope, startpos, winsize, pchdisp, timedisp; long remain, nsmps, idelaylen, irdphase, irdphaseb, iwrphase, mask, idsamp; long counter, stage, framesize; RGET; nout = *out; nin = *in; pchratio = m_pitchratio; winsize = m_windowsize; pchdisp = m_pitchdispersion; timedisp = m_timedispersion; timedisp = sc_clip(timedisp, 0.f, winsize) * SAMPLERATE; dlybuf = m_dlybuf; fdelaylen = m_fdelaylen; idelaylen = m_idelaylen; iwrphase = m_iwrphase; counter = m_counter; stage = m_stage; mask = m_mask; framesize = m_framesize; dsamp1 = m_dsamp1; dsamp2 = m_dsamp2; dsamp3 = m_dsamp3; dsamp4 = m_dsamp4; dsamp1_slope = m_dsamp1_slope; dsamp2_slope = m_dsamp2_slope; dsamp3_slope = m_dsamp3_slope; dsamp4_slope = m_dsamp4_slope; ramp1 = m_ramp1; ramp2 = m_ramp2; ramp3 = m_ramp3; ramp4 = m_ramp4; ramp1_slope = m_ramp1_slope; ramp2_slope = m_ramp2_slope; ramp3_slope = m_ramp3_slope; ramp4_slope = m_ramp4_slope; slope = m_slope; remain = n; while (remain) { if (counter <= 0) { counter = framesize >> 2; m_stage = stage = (stage + 1) & 3; disppchratio = pchratio; if (pchdisp != 0.f) { disppchratio += (pchdisp * frand2(s1,s2,s3)); } disppchratio = sc_clip(disppchratio, 0.f, 4.f); pchratio1 = disppchratio - 1.f; samp_slope = -pchratio1; startpos = pchratio1 < 0.f ? 2.f : framesize * pchratio1 + 2.f; startpos += (timedisp * frand(s1,s2,s3)); switch(stage) { case 0 : m_dsamp1_slope = dsamp1_slope = samp_slope; dsamp1 = startpos; ramp1 = 0.0; m_ramp1_slope = ramp1_slope = slope; m_ramp3_slope = ramp3_slope = -slope; break; case 1 : m_dsamp2_slope = dsamp2_slope = samp_slope; dsamp2 = startpos; ramp2 = 0.0; m_ramp2_slope = ramp2_slope = slope; m_ramp4_slope = ramp4_slope = -slope; break; case 2 : m_dsamp3_slope = dsamp3_slope = samp_slope; dsamp3 = startpos; ramp3 = 0.0; m_ramp3_slope = ramp3_slope = slope; m_ramp1_slope = ramp1_slope = -slope; break; case 3 : m_dsamp4_slope = dsamp4_slope = samp_slope; dsamp4 = startpos; ramp4 = 0.0; m_ramp2_slope = ramp2_slope = -slope; m_ramp4_slope = ramp4_slope = slope; break; } } nsmps = sc_min(remain, counter); remain -= nsmps; counter -= nsmps; for (int i = 0; i!= nsmps;++i) { iwrphase = (iwrphase + 1) & mask; dsamp1 += dsamp1_slope; idsamp = (long)dsamp1; frac = dsamp1 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value = (d1 + frac * (d2 - d1)) * ramp1; ramp1 += ramp1_slope; dsamp2 += dsamp2_slope; idsamp = (long)dsamp2; frac = dsamp2 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp2; ramp2 += ramp2_slope; dsamp3 += dsamp3_slope; idsamp = (long)dsamp3; frac = dsamp3 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp3; ramp3 += ramp3_slope; dsamp4 += dsamp4_slope; idsamp = (long)dsamp4; frac = dsamp4 - idsamp; irdphase = (iwrphase - idsamp) & mask; irdphaseb = (irdphase - 1) & mask; d1 = dlybuf[irdphase]; d2 = dlybuf[irdphaseb]; value += (d1 + frac * (d2 - d1)) * ramp4; ramp4 += ramp4_slope; dlybuf[iwrphase] = ZXP(nin); ZXP(nout) = value *= 0.5; } } m_counter = counter; m_dsamp1 = dsamp1; m_dsamp2 = dsamp2; m_dsamp3 = dsamp3; m_dsamp4 = dsamp4; m_ramp1 = ramp1; m_ramp2 = ramp2; m_ramp3 = ramp3; m_ramp4 = ramp4; m_iwrphase = iwrphase; RPUT; } /* todo: does a control rate PitchShift make sense? */

1 0

externals/sc4pd/pd combc.pd,NONE,1.1
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/pd In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv17161/pd Added Files: combc.pd Log Message: --- NEW FILE: combc.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 67 187 dac~; #X msg 185 88 delaytime \$1; #X obj 111 190 print~; #X obj 128 170 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X floatatom 316 83 5 0 0 0 - - -; #X msg 316 107 decaytime \$1; #X obj 79 57 WhiteNoise~; #X obj 77 138 CombN~ 1 0.001 0.2; #X connect 0 0 2 0; #X connect 2 0 8 0; #X connect 4 0 3 0; #X connect 5 0 6 0; #X connect 6 0 8 0; #X connect 7 0 8 0; #X connect 8 0 1 0; #X connect 8 0 1 1; #X connect 8 0 3 0;

1 0

externals/sc4pd/source CombC.cpp,NONE,1.1
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/source In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv17161/source Added Files: CombC.cpp Log Message: --- NEW FILE: CombC.cpp --- /* sc4pd CombC~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: Butcher/Charles/Doerner: The Contest Of Pleasures */ #include "sc4pd.hpp" #include "DelayUnit.hpp" class CombC_ar : public FeedbackDelay_ar { FLEXT_HEADER(CombC_ar,FeedbackDelay_ar); CombC_ar (int argc, t_atom *argv); ~CombC_ar (); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out) { delay_changed = decay_changed = false; FeedbackDelay_Reset(); } void m_delay(float f) { m_delaytime=f; delay_changed = true; } void m_decay(float f) { m_decaytime=f; decay_changed = true; } private: bool delay_changed, decay_changed; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_delay); FLEXT_CALLBACK_F(m_decay); }; FLEXT_LIB_DSP_V("CombC~",CombC_ar); CombC_ar::CombC_ar (int argc, t_atom *argv) { FLEXT_ADDMETHOD_(0,"delaytime",m_delay); FLEXT_ADDMETHOD_(0,"decaytime",m_decay); //parse arguments AtomList Args(argc,argv); if (Args.Count() != 3) { post("3 arguments are needed"); return; } m_maxdelaytime = sc_getfloatarg(Args,0); m_delaytime = sc_getfloatarg(Args,1); m_decaytime = sc_getfloatarg(Args,2); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); AddOutSignal(); } CombC_ar::~CombC_ar () { DelayUnit_Dtor(); } void CombC_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; float d0, d1, d2, d3; if (delay_changed || decay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for (int i = 0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; if (irdphase0 < 0) { dlybuf[iwrphase & mask] = ZXP(nin); ZXP(nout) = 0.f; } else { if (irdphase1 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; } else if (irdphase2 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; } else if (irdphase3 < 0) { d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; } else { d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; d3 = dlybuf[irdphase3 & mask]; } float value = cubicinterp(frac, d0, d1, d2, d3); dlybuf[iwrphase & mask] = ZXP(nin) + feedbk * value; ZXP(nout) = value; } feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for (int i = 0; i!= n;++i) { long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; if (irdphase0 < 0) { dlybuf[iwrphase & mask] = ZXP(nin); ZXP(nout) = 0.f; } else { if (irdphase1 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; } else if (irdphase2 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; } else if (irdphase3 < 0) { d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; } else { d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; d3 = dlybuf[irdphase3 & mask]; } float value = cubicinterp(frac, d0, d1, d2, d3); dlybuf[iwrphase & mask] = ZXP(nin) + feedbk * value; ZXP(nout) = value; } iwrphase++; } } m_iwrphase = iwrphase; m_numoutput += n; if (m_numoutput >= m_idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_)); } } void CombC_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; float d0, d1, d2, d3; if(delay_changed || decay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for(int i=0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; float d0 = dlybuf[irdphase0 & mask]; float d1 = dlybuf[irdphase1 & mask]; float d2 = dlybuf[irdphase2 & mask]; float d3 = dlybuf[irdphase3 & mask]; float value = cubicinterp(frac, d0, d1, d2, d3); dlybuf[iwrphase & mask] = ZXP(nin) + feedbk * value; ZXP(nout) = value; feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for(int i=0; i!= n;++i) { long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; float d0 = dlybuf[irdphase0 & mask]; float d1 = dlybuf[irdphase1 & mask]; float d2 = dlybuf[irdphase2 & mask]; float d3 = dlybuf[irdphase3 & mask]; float value = cubicinterp(frac, d0, d1, d2, d3); dlybuf[iwrphase & mask] = ZXP(nin) + feedbk * value; ZXP(nout) = value; iwrphase++; } } m_iwrphase = iwrphase; } /* todo: CombC for control rate ? */

1 0

externals/sc4pd make-files.txt,1.16,1.17 makefile.pd-linux,1.2,1.3
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv11639 Modified Files: make-files.txt makefile.pd-linux Log Message: some allpass delay lines Index: make-files.txt =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/make-files.txt,v retrieving revision 1.16 retrieving revision 1.17 diff -C2 -d -r1.16 -r1.17 *** make-files.txt 4 Aug 2004 16:43:25 -0000 1.16 --- make-files.txt 4 Aug 2004 18:23:14 -0000 1.17 *************** *** 12,14 **** LFPulse.cpp Impulse.cpp Integrator.cpp Decay.cpp Decay2.cpp Lag.cpp \ Lag2.cpp Lag3.cpp LinExp.cpp DelayUnit.cpp DelayN.cpp DelayL.cpp \ ! DelayC.cpp CombN.cpp CombL.cpp --- 12,15 ---- LFPulse.cpp Impulse.cpp Integrator.cpp Decay.cpp Decay2.cpp Lag.cpp \ Lag2.cpp Lag3.cpp LinExp.cpp DelayUnit.cpp DelayN.cpp DelayL.cpp \ ! DelayC.cpp CombN.cpp CombL.cpp CombC.cpp AllpassN.cpp AllpassL.cpp \ ! AllpassC.cpp Index: makefile.pd-linux =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/makefile.pd-linux,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** makefile.pd-linux 16 Jul 2004 16:43:47 -0000 1.2 --- makefile.pd-linux 4 Aug 2004 18:23:14 -0000 1.3 *************** *** 46,50 **** $(patsubst %,$(SRCDIR)/%,$(SRCS)): $(patsubst %,$(SRCDIR)/%,$(HDRS)) $(MAKEFILE) $(CONFIG) ! touch $(patsubst %,$(SRCDIR)/%,$(SRCS)) $(TARGDIR): --- 46,50 ---- $(patsubst %,$(SRCDIR)/%,$(SRCS)): $(patsubst %,$(SRCDIR)/%,$(HDRS)) $(MAKEFILE) $(CONFIG) ! # touch $(patsubst %,$(SRCDIR)/%,$(SRCS)) $(TARGDIR):

1 0

externals/sc4pd/source AllpassC.cpp,NONE,1.1 AllpassL.cpp,NONE,1.1 AllpassN.cpp,NONE,1.1 main.cpp,1.18,1.19
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/source In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv11639/source Modified Files: main.cpp Added Files: AllpassC.cpp AllpassL.cpp AllpassN.cpp Log Message: some allpass delay lines Index: main.cpp =================================================================== RCS file: /cvsroot/pure-data/externals/sc4pd/source/main.cpp,v retrieving revision 1.18 retrieving revision 1.19 diff -C2 -d -r1.18 -r1.19 *** main.cpp 4 Aug 2004 16:43:26 -0000 1.18 --- main.cpp 4 Aug 2004 18:23:15 -0000 1.19 *************** *** 61,65 **** " Integrator(~), Decay~, Decay2~, Lag~, Lag2~, LinExp(~), " "DelayN~,\n" ! " DelayL~, DelayC~, CombN~, CombL~" ); --- 61,68 ---- " Integrator(~), Decay~, Decay2~, Lag~, Lag2~, LinExp(~), " "DelayN~,\n" ! " DelayL~, DelayC~, CombN~, CombL~, CombC~, AllpassN~, " ! "AllpassL~,\n" ! " AllpassC~" ! "\n" ); *************** *** 218,221 **** --- 221,232 ---- FLEXT_DSP_SETUP(CombL_ar); + + FLEXT_DSP_SETUP(CombC_ar); + + FLEXT_DSP_SETUP(AllpassN_ar); + + FLEXT_DSP_SETUP(AllpassL_ar); + + FLEXT_DSP_SETUP(AllpassC_ar); } --- NEW FILE: AllpassN.cpp --- /* sc4pd AllpassN~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: Efzeg: Boogie */ #include "sc4pd.hpp" #include "DelayUnit.hpp" class AllpassN_ar : public FeedbackDelay_ar { FLEXT_HEADER(AllpassN_ar,FeedbackDelay_ar); AllpassN_ar (int argc, t_atom *argv); ~AllpassN_ar (); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out) { delay_changed = decay_changed = false; FeedbackDelay_Reset(); } void m_delay(float f) { m_delaytime=f; delay_changed = true; } void m_decay(float f) { m_decaytime=f; decay_changed = true; } private: bool delay_changed, decay_changed; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_delay); FLEXT_CALLBACK_F(m_decay); }; FLEXT_LIB_DSP_V("AllpassN~",AllpassN_ar); AllpassN_ar::AllpassN_ar (int argc, t_atom *argv) { FLEXT_ADDMETHOD_(0,"delaytime",m_delay); FLEXT_ADDMETHOD_(0,"decaytime",m_decay); //parse arguments AtomList Args(argc,argv); if (Args.Count() != 3) { post("3 arguments are needed"); return; } m_maxdelaytime = sc_getfloatarg(Args,0); m_delaytime = sc_getfloatarg(Args,1); m_decaytime = sc_getfloatarg(Args,2); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); AddOutSignal(); } AllpassN_ar::~AllpassN_ar () { DelayUnit_Dtor(); } void AllpassN_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; if (delay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for (int i = 0; i!= n;++i) { dsamp += dsamp_slope; long irdphase = iwrphase - (long)dsamp; if (irdphase < 0) { float dwr = ZXP(nin); dlybuf[iwrphase & mask] = dwr; ZXP(nout) = -feedbk * dwr; } else { float value = dlybuf[irdphase & mask]; float dwr = feedbk * value + ZXP(nin); dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long irdphase = iwrphase - (long)dsamp; float* dlybuf1 = dlybuf - ZOFF; float* dlyN = dlybuf1 + m_idelaylen; if (decay_changed) { float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); long remain = n; while (remain) { float* dlywr = dlybuf1 + (iwrphase & mask); float* dlyrd = dlybuf1 + (irdphase & mask); long rdspace = dlyN - dlyrd; long wrspace = dlyN - dlywr; long nsmps = sc_min(rdspace, wrspace); nsmps = sc_min(remain, nsmps); remain -= nsmps; if (irdphase < 0) { dlyrd += nsmps; for (int i = 0; i!= nsmps;++i) { float dwr = ZXP(nin); ZXP(dlywr) = dwr; ZXP(nout) = -feedbk * dwr; feedbk += feedbk_slope; } } else { for (int i = 0; i!= nsmps;++i) { float x1 = ZXP(dlyrd); float dwr = x1 * feedbk + ZXP(nin); ZXP(dlywr) = dwr; ZXP(nout) = x1 - feedbk * dwr; feedbk += feedbk_slope; } } iwrphase += nsmps; irdphase += nsmps; } decay_changed=false; } else { long remain = n; while (remain) { float* dlyrd = dlybuf1 + (irdphase & mask); float* dlywr = dlybuf1 + (iwrphase & mask); long rdspace = dlyN - dlyrd; long wrspace = dlyN - dlywr; long nsmps = sc_min(rdspace, wrspace); nsmps = sc_min(remain, nsmps); remain -= nsmps; if (irdphase < 0) { feedbk = -feedbk; for (int i = 0; i!= nsmps;++i) { float dwr = ZXP(nin); ZXP(dlywr) = dwr; ZXP(nout) = feedbk * dwr; } feedbk = -feedbk; } else { for (int i = 0; i!= nsmps;++i) { float x1 = ZXP(dlyrd); float dwr = x1 * feedbk + ZXP(nin); ZXP(dlywr) = dwr; ZXP(nout) = x1 - feedbk * dwr; } } iwrphase += nsmps; irdphase += nsmps; } m_feedbk = feedbk; } } m_iwrphase = iwrphase; m_numoutput += n; if (m_numoutput >= m_idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_)); } } void AllpassN_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; if(delay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for(int i=0; i!= n;++i) { dsamp += dsamp_slope; ++iwrphase; long irdphase = iwrphase - (long)dsamp; float value = dlybuf[irdphase & mask]; float dwr = value * feedbk + ZXP(nin); dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; feedbk += feedbk_slope; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long irdphase = iwrphase - (long)dsamp; float* dlybuf1 = dlybuf - ZOFF; float* dlyrd = dlybuf1 + (irdphase & mask); float* dlywr = dlybuf1 + (iwrphase & mask); float* dlyN = dlybuf1 + m_idelaylen; if(decay_changed) { float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); long remain = n; while (remain) { long rdspace = dlyN - dlyrd; long wrspace = dlyN - dlywr; long nsmps = sc_min(rdspace, wrspace); nsmps = sc_min(remain, nsmps); remain -= nsmps; for(int i=0; i!= nsmps;++i) { float value = ZXP(dlyrd); float dwr = value * feedbk + ZXP(nin); ZXP(dlywr) = dwr; ZXP(nout) = value - feedbk * dwr; feedbk += feedbk_slope; } if (dlyrd == dlyN) dlyrd = dlybuf1; if (dlywr == dlyN) dlywr = dlybuf1; } m_feedbk = feedbk; decay_changed = false; } else { long remain = n; while (remain) { long rdspace = dlyN - dlyrd; long wrspace = dlyN - dlywr; long nsmps = sc_min(rdspace, wrspace); nsmps = sc_min(remain, nsmps); remain -= nsmps; for (int i = 0; i!= nsmps; ++i) { float value = ZXP(dlyrd); float dwr = value * feedbk + ZXP(nin); ZXP(dlywr) = dwr; ZXP(nout) = value - feedbk * dwr; } if (dlyrd == dlyN) dlyrd = dlybuf1; if (dlywr == dlyN) dlywr = dlybuf1; } } iwrphase += n; } m_iwrphase = iwrphase; } /* todo: AllpassN for control rate ? */ --- NEW FILE: AllpassL.cpp --- /* sc4pd AllpassL~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: Rashid Ali & Frank Lowe: Duo Exchange */ #include "sc4pd.hpp" #include "DelayUnit.hpp" class AllpassL_ar : public FeedbackDelay_ar { FLEXT_HEADER(AllpassL_ar,FeedbackDelay_ar); AllpassL_ar (int argc, t_atom *argv); ~AllpassL_ar (); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out) { delay_changed = decay_changed = false; FeedbackDelay_Reset(); } void m_delay(float f) { m_delaytime=f; delay_changed = true; } void m_decay(float f) { m_decaytime=f; decay_changed = true; } private: bool delay_changed, decay_changed; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_delay); FLEXT_CALLBACK_F(m_decay); }; FLEXT_LIB_DSP_V("AllpassL~",AllpassL_ar); AllpassL_ar::AllpassL_ar (int argc, t_atom *argv) { FLEXT_ADDMETHOD_(0,"delaytime",m_delay); FLEXT_ADDMETHOD_(0,"decaytime",m_decay); //parse arguments AtomList Args(argc,argv); if (Args.Count() != 3) { post("3 arguments are needed"); return; } m_maxdelaytime = sc_getfloatarg(Args,0); m_delaytime = sc_getfloatarg(Args,1); m_decaytime = sc_getfloatarg(Args,2); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); AddOutSignal(); } AllpassL_ar::~AllpassL_ar () { DelayUnit_Dtor(); } void AllpassL_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; if (delay_changed || decay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for (int i = 0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; float zin = ZXP(nin); if (irdphase < 0) { dlybuf[iwrphase & mask] = zin; ZXP(nout) = - feedbk * zin; //check: probably a bug = 0? } else if (irdphaseb < 0) { float d1 = dlybuf[irdphase & mask]; float value = d1 - frac * d1; float dwr = zin + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } else { float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; float value = lininterp(frac, d1, d2); float dwr = zin + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; float zin = ZXP(nin); for (int i = 0; i!= n;++i) { long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; float zin = ZXP(nin); if (irdphase < 0) { dlybuf[iwrphase & mask] = zin; ZXP(nout) = - feedbk * zin; } else if (irdphaseb < 0) { float d1 = dlybuf[irdphase & mask]; float value = d1 - frac * d1; float dwr = zin + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } else { float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; float value = lininterp(frac, d1, d2); float dwr = zin + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } iwrphase++; } } m_iwrphase = iwrphase; m_numoutput += n; if (m_numoutput >= m_idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_)); } } void AllpassL_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; if(delay_changed || decay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for(int i=0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; float value = lininterp(frac, d1, d2); float dwr = ZXP(nin) + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for(int i=0; i!= n;++i) { long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; float value = lininterp(frac, d1, d2); float dwr = ZXP(nin) + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; iwrphase++; } } m_iwrphase = iwrphase; } /* todo: AllpassL for control rate ? */ --- NEW FILE: AllpassC.cpp --- /* sc4pd AllpassC~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: Rashid Ali & Frank Lowe: Duo Exchange */ #include "sc4pd.hpp" #include "DelayUnit.hpp" class AllpassC_ar : public FeedbackDelay_ar { FLEXT_HEADER(AllpassC_ar,FeedbackDelay_ar); AllpassC_ar (int argc, t_atom *argv); ~AllpassC_ar (); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out) { delay_changed = decay_changed = false; FeedbackDelay_Reset(); } void m_delay(float f) { m_delaytime=f; delay_changed = true; } void m_decay(float f) { m_decaytime=f; decay_changed = true; } private: bool delay_changed, decay_changed; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_delay); FLEXT_CALLBACK_F(m_decay); }; FLEXT_LIB_DSP_V("AllpassC~",AllpassC_ar); AllpassC_ar::AllpassC_ar (int argc, t_atom *argv) { FLEXT_ADDMETHOD_(0,"delaytime",m_delay); FLEXT_ADDMETHOD_(0,"decaytime",m_decay); //parse arguments AtomList Args(argc,argv); if (Args.Count() != 3) { post("3 arguments are needed"); return; } m_maxdelaytime = sc_getfloatarg(Args,0); m_delaytime = sc_getfloatarg(Args,1); m_decaytime = sc_getfloatarg(Args,2); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); AddOutSignal(); } AllpassC_ar::~AllpassC_ar () { DelayUnit_Dtor(); } void AllpassC_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; float d0, d1, d2, d3; if (delay_changed || decay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for (int i = 0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; if (irdphase0 < 0) { dlybuf[iwrphase & mask] = ZXP(nin); ZXP(nout) = 0.f; } else { if (irdphase1 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; } else if (irdphase2 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; } else if (irdphase3 < 0) { d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; } else { d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; d3 = dlybuf[irdphase3 & mask]; } float value = cubicinterp(frac, d0, d1, d2, d3); float dwr = ZXP(nin) + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for (int i = 0; i!= n;++i) { long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; if (irdphase0 < 0) { dlybuf[iwrphase & mask] = ZXP(nin); ZXP(nout) = 0.f; } else { if (irdphase1 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; } else if (irdphase2 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; } else if (irdphase3 < 0) { d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; } else { d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; d3 = dlybuf[irdphase3 & mask]; } float value = cubicinterp(frac, d0, d1, d2, d3); float dwr = ZXP(nin) + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; } iwrphase++; } } m_iwrphase = iwrphase; m_numoutput += n; if (m_numoutput >= m_idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_)); } } void AllpassC_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; float feedbk = m_feedbk; long mask = m_mask; float d0, d1, d2, d3; if(delay_changed || decay_changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); float next_feedbk = CalcFeedback(m_delaytime, m_decaytime); float feedbk_slope = CALCSLOPE(next_feedbk, feedbk); for(int i=0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; float d0 = dlybuf[irdphase0 & mask]; float d1 = dlybuf[irdphase1 & mask]; float d2 = dlybuf[irdphase2 & mask]; float d3 = dlybuf[irdphase3 & mask]; float value = cubicinterp(frac, d0, d1, d2, d3); float dwr = ZXP(nin) + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; feedbk += feedbk_slope; iwrphase++; } m_feedbk = feedbk; m_dsamp = dsamp; delay_changed = decay_changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for(int i=0; i!= n;++i) { long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; float d0 = dlybuf[irdphase0 & mask]; float d1 = dlybuf[irdphase1 & mask]; float d2 = dlybuf[irdphase2 & mask]; float d3 = dlybuf[irdphase3 & mask]; float value = cubicinterp(frac, d0, d1, d2, d3); float dwr = ZXP(nin) + feedbk * value; dlybuf[iwrphase & mask] = dwr; ZXP(nout) = value - feedbk * dwr; iwrphase++; } } m_iwrphase = iwrphase; } /* todo: AllpassC for control rate ? */

1 0

externals/sc4pd/pd allpassc.pd,NONE,1.1 allpassl.pd,NONE,1.1 allpassn.pd,NONE,1.1
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/pd In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv11639/pd Added Files: allpassc.pd allpassl.pd allpassn.pd Log Message: some allpass delay lines --- NEW FILE: allpassn.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 67 187 dac~; #X msg 185 88 delaytime \$1; #X obj 111 190 print~; #X obj 128 170 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X floatatom 316 83 5 0 0 0 - - -; #X msg 316 107 decaytime \$1; #X obj 77 138 AllpassN~ 1 0.001 0.2; #X obj 79 57 Dust~ 4; #X connect 0 0 2 0; #X connect 2 0 7 0; #X connect 4 0 3 0; #X connect 5 0 6 0; #X connect 6 0 7 0; #X connect 7 0 1 0; #X connect 7 0 1 1; #X connect 7 0 3 0; #X connect 8 0 7 0; --- NEW FILE: allpassc.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 67 187 dac~; #X msg 185 88 delaytime \$1; #X obj 111 190 print~; #X obj 128 170 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X floatatom 316 83 5 0 0 0 - - -; #X msg 316 107 decaytime \$1; #X obj 77 138 AllpassC~ 1 0.001 0.2; #X obj 79 57 Dust~ 55; #X connect 0 0 2 0; #X connect 2 0 7 0; #X connect 4 0 3 0; #X connect 5 0 6 0; #X connect 6 0 7 0; #X connect 7 0 1 0; #X connect 7 0 1 1; #X connect 7 0 3 0; #X connect 8 0 7 0; --- NEW FILE: allpassl.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 67 187 dac~; #X msg 185 88 delaytime \$1; #X obj 111 190 print~; #X obj 128 170 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X floatatom 316 83 5 0 0 0 - - -; #X msg 316 107 decaytime \$1; #X obj 79 57 Dust~ 4; #X obj 77 138 AllpassL~ 1 0.001 0.2; #X connect 0 0 2 0; #X connect 2 0 8 0; #X connect 4 0 3 0; #X connect 5 0 6 0; #X connect 6 0 8 0; #X connect 7 0 8 0; #X connect 8 0 1 0; #X connect 8 0 1 1; #X connect 8 0 3 0;

1 0

externals/sc4pd/source DelayC.cpp,NONE,1.1 DelayL.cpp,NONE,1.1
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/source In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv26400/source Added Files: DelayC.cpp DelayL.cpp Log Message: --- NEW FILE: DelayL.cpp --- /* sc4pd DelayL~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: Tom Cora: Halleluja, Anyway */ #include "sc4pd.hpp" #include "DelayUnit.hpp" class DelayL_ar : private DelayUnit_ar { FLEXT_HEADER(DelayL_ar,DelayUnit_ar); DelayL_ar (int argc, t_atom *argv); ~DelayL_ar (); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out) { changed = false; DelayUnit_Reset(); } void m_set(float f) { m_delaytime=f; changed = true; } private: bool changed; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_set); }; FLEXT_LIB_DSP_V("DelayL~",DelayL_ar); DelayL_ar::DelayL_ar (int argc, t_atom *argv) { FLEXT_ADDMETHOD_(0,"delaytime",m_set); //parse arguments AtomList Args(argc,argv); if (Args.Count() != 2) { post("2 arguments are needed"); return; } m_delaytime = sc_getfloatarg(Args,0); m_maxdelaytime = sc_getfloatarg(Args,1); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); AddOutSignal(); } DelayL_ar::~DelayL_ar () { DelayUnit_Dtor(); } void DelayL_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; long mask = m_mask; if (changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); for (int i = 0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; dlybuf[iwrphase & mask] = ZXP(nin); if (irdphase < 0) { ZXP(nout) = 0.f; } else if (irdphaseb < 0) { float d1 = dlybuf[irdphase & mask]; ZXP(nout) = d1 - frac * d1; } else { float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; ZXP(nout) = lininterp(frac, d1, d2); } iwrphase++; } m_dsamp = dsamp; changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for (int i = 0; i!= n;++i) { long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; dlybuf[iwrphase & mask] = ZXP(nin); if (irdphase < 0) { ZXP(nout) = 0.f; } else if (irdphaseb < 0) { float d1 = dlybuf[irdphase & mask]; ZXP(nout) = d1 - frac * d1; } else { float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; ZXP(nout) = lininterp(frac, d1, d2); } iwrphase++; } } m_iwrphase = iwrphase; m_numoutput += n; if (m_numoutput >= m_idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_)); } } void DelayL_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; long mask = m_mask; if (changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); for (int i = 0; i!= n;++i) { dlybuf[iwrphase & mask] = ZXP(nin); dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; ZXP(nout) = lininterp(frac, d1, d2); iwrphase++; } m_dsamp = dsamp; changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for (int i = 0; i!= n;++i) { dlybuf[iwrphase & mask] = ZXP(nin); long irdphase = iwrphase - idsamp; long irdphaseb = irdphase - 1; float d1 = dlybuf[irdphase & mask]; float d2 = dlybuf[irdphaseb & mask]; ZXP(nout) = lininterp(frac, d1, d2); iwrphase++; } } m_iwrphase = iwrphase; } /* todo: DelayL for control rate ? */ --- NEW FILE: DelayC.cpp --- /* sc4pd DelayC~ Copyright (c) 2004 Tim Blechmann. This code is derived from: SuperCollider real time audio synthesis system Copyright (c) 2002 James McCartney. All rights reserved. http://www.audiosynth.com This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. Based on: PureData by Miller Puckette and others. http://www.crca.ucsd.edu/~msp/software.html FLEXT by Thomas Grill http://www.parasitaere-kapazitaeten.net/ext SuperCollider by James McCartney http://www.audiosynth.com Coded while listening to: Tom Cora: Halleluja, Anyway */ #include "sc4pd.hpp" #include "DelayUnit.hpp" class DelayC_ar : private DelayUnit_ar { FLEXT_HEADER(DelayC_ar,DelayUnit_ar); DelayC_ar (int argc, t_atom *argv); ~DelayC_ar (); protected: virtual void m_signal(int n, t_sample *const *in, t_sample *const *out) { m_signal_fun(n,in,out); } virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out) { changed = false; DelayUnit_Reset(); } void m_set(float f) { m_delaytime=f; changed = true; } private: bool changed; DEFSIGCALL(m_signal_fun); DEFSIGFUN(m_signal_); DEFSIGFUN(m_signal_z); FLEXT_CALLBACK_F(m_set); }; FLEXT_LIB_DSP_V("DelayC~",DelayC_ar); DelayC_ar::DelayC_ar (int argc, t_atom *argv) { FLEXT_ADDMETHOD_(0,"delaytime",m_set); //parse arguments AtomList Args(argc,argv); if (Args.Count() != 2) { post("2 arguments are needed"); return; } m_delaytime = sc_getfloatarg(Args,0); m_maxdelaytime = sc_getfloatarg(Args,1); SETSIGFUN(m_signal_fun,SIGFUN(m_signal_z)); AddOutSignal(); } DelayC_ar::~DelayC_ar () { DelayUnit_Dtor(); } void DelayC_ar::m_signal_z(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; long mask = m_mask; float d0, d1, d2, d3; if (changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); for (int i = 0; i!= n;++i) { dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; dlybuf[iwrphase & mask] = ZXP(nin); if (irdphase0 < 0) { ZXP(nout) = 0.f; } else { if (irdphase1 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; } else if (irdphase2 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; } else if (irdphase3 < 0) { d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; } else { d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; d3 = dlybuf[irdphase3 & mask]; } ZXP(nout) = cubicinterp(frac, d0, d1, d2, d3); } iwrphase++; } m_dsamp = dsamp; changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for (int i = 0; i!= n;++i) { long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; dlybuf[iwrphase & mask] = ZXP(nin); if (irdphase0 < 0) { ZXP(nout) = 0.f; } else { if (irdphase1 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; } else if (irdphase2 < 0) { d1 = d2 = d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; } else if (irdphase3 < 0) { d3 = 0.f; d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; } else { d0 = dlybuf[irdphase0 & mask]; d1 = dlybuf[irdphase1 & mask]; d2 = dlybuf[irdphase2 & mask]; d3 = dlybuf[irdphase3 & mask]; } ZXP(nout) = cubicinterp(frac, d0, d1, d2, d3); } iwrphase++; } } m_iwrphase = iwrphase; m_numoutput += n; if (m_numoutput >= m_idelaylen) { SETSIGFUN(m_signal_fun,SIGFUN(m_signal_)); } } void DelayC_ar::m_signal_(int n, t_sample *const *in, t_sample *const *out) { t_sample *nin = *in; t_sample *nout = *out; float *dlybuf = m_dlybuf; long iwrphase = m_iwrphase; float dsamp = m_dsamp; long mask = m_mask; if (changed) { float next_dsamp = CalcDelay(m_delaytime); float dsamp_slope = CALCSLOPE(next_dsamp, dsamp); for (int i = 0; i!= n;++i) { dlybuf[iwrphase & mask] = ZXP(nin); dsamp += dsamp_slope; long idsamp = (long)dsamp; float frac = dsamp - idsamp; long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; float d0 = dlybuf[irdphase0 & mask]; float d1 = dlybuf[irdphase1 & mask]; float d2 = dlybuf[irdphase2 & mask]; float d3 = dlybuf[irdphase3 & mask]; ZXP(nout) = cubicinterp(frac, d0, d1, d2, d3); iwrphase++; } m_dsamp = dsamp; changed = false; } else { long idsamp = (long)dsamp; float frac = dsamp - idsamp; for (int i = 0; i!= n;++i) { dlybuf[iwrphase & mask] = ZXP(nin); long irdphase1 = iwrphase - idsamp; long irdphase2 = irdphase1 - 1; long irdphase3 = irdphase1 - 2; long irdphase0 = irdphase1 + 1; float d0 = dlybuf[irdphase0 & mask]; float d1 = dlybuf[irdphase1 & mask]; float d2 = dlybuf[irdphase2 & mask]; float d3 = dlybuf[irdphase3 & mask]; ZXP(nout) = cubicinterp(frac, d0, d1, d2, d3); iwrphase++; } } m_iwrphase = iwrphase; } /* todo: DelayC for control rate ? */

1 0

externals/sc4pd/pd combl.pd,NONE,1.1 delayc.pd,NONE,1.1 delayl.pd,NONE,1.1
by Tim Blechmann 04 Aug '04

04 Aug '04

Update of /cvsroot/pure-data/externals/sc4pd/pd In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv26400/pd Added Files: combl.pd delayc.pd delayl.pd Log Message: --- NEW FILE: combl.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 67 187 dac~; #X msg 185 88 delaytime \$1; #X obj 111 190 print~; #X obj 128 170 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X floatatom 316 83 5 0 0 0 - - -; #X msg 316 107 decaytime \$1; #X obj 79 57 WhiteNoise~; #X obj 77 138 CombL~ 1 0.001 0.2; #X connect 0 0 2 0; #X connect 2 0 8 0; #X connect 4 0 3 0; #X connect 5 0 6 0; #X connect 6 0 8 0; #X connect 7 0 8 0; #X connect 8 0 1 0; #X connect 8 0 1 1; #X connect 8 0 3 0; --- NEW FILE: delayl.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 90 233 dac~; #X msg 185 88 delaytime \$1; #X obj 79 57 Dust~ 4; #X obj 77 138 DelayL~ 1 3; #X connect 0 0 2 0; #X connect 2 0 4 0; #X connect 3 0 1 1; #X connect 3 0 4 0; #X connect 4 0 1 0; --- NEW FILE: delayc.pd --- #N canvas 0 0 450 300 10; #X floatatom 188 40 5 0 0 0 - - -; #X obj 67 187 dac~; #X msg 185 88 delaytime \$1; #X obj 77 138 CombN~ 0.01 0.003 3; #X obj 79 57 Dust~ 440; #X connect 0 0 2 0; #X connect 2 0 3 0; #X connect 3 0 1 0; #X connect 3 0 1 1; #X connect 4 0 3 0;

1 0

Jump to page:

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

2003

Pd-cvs August 2004