From 656b72dd1a85069ad57351147553e07226a990b6 Mon Sep 17 00:00:00 2001 From: Tim Blechmann Date: Wed, 4 Aug 2004 18:23:15 +0000 Subject: some allpass delay lines svn path=/trunk/externals/tb/; revision=1909 --- sc4pd/make-files.txt | 3 +- sc4pd/makefile.pd-linux | 2 +- sc4pd/pd/allpassc.pd | 20 +++ sc4pd/pd/allpassl.pd | 20 +++ sc4pd/pd/allpassn.pd | 20 +++ sc4pd/source/AllpassC.cpp | 315 +++++++++++++++++++++++++++++++++++++++++ sc4pd/source/AllpassL.cpp | 276 ++++++++++++++++++++++++++++++++++++ sc4pd/source/AllpassN.cpp | 348 ++++++++++++++++++++++++++++++++++++++++++++++ sc4pd/source/main.cpp | 13 +- 9 files changed, 1014 insertions(+), 3 deletions(-) create mode 100644 sc4pd/pd/allpassc.pd create mode 100644 sc4pd/pd/allpassl.pd create mode 100644 sc4pd/pd/allpassn.pd create mode 100644 sc4pd/source/AllpassC.cpp create mode 100644 sc4pd/source/AllpassL.cpp create mode 100644 sc4pd/source/AllpassN.cpp diff --git a/sc4pd/make-files.txt b/sc4pd/make-files.txt index 99913e7..3f9038d 100755 --- a/sc4pd/make-files.txt +++ b/sc4pd/make-files.txt @@ -11,4 +11,5 @@ SRCS= \ difsqr.cpp sumsqr.cpp sqrdif.cpp sqrsum.cpp absdif.cpp LFSaw.cpp \ 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 + DelayC.cpp CombN.cpp CombL.cpp CombC.cpp AllpassN.cpp AllpassL.cpp \ + AllpassC.cpp diff --git a/sc4pd/makefile.pd-linux b/sc4pd/makefile.pd-linux index fc55ec6..99d62ef 100644 --- a/sc4pd/makefile.pd-linux +++ b/sc4pd/makefile.pd-linux @@ -45,7 +45,7 @@ TARGET=$(TARGDIR)/$(NAME).pd_linux all: $(TARGDIR) $(TARGET) $(patsubst %,$(SRCDIR)/%,$(SRCS)): $(patsubst %,$(SRCDIR)/%,$(HDRS)) $(MAKEFILE) $(CONFIG) - touch $(patsubst %,$(SRCDIR)/%,$(SRCS)) +# touch $(patsubst %,$(SRCDIR)/%,$(SRCS)) $(TARGDIR): -mkdir $(TARGDIR) diff --git a/sc4pd/pd/allpassc.pd b/sc4pd/pd/allpassc.pd new file mode 100644 index 0000000..297cb91 --- /dev/null +++ b/sc4pd/pd/allpassc.pd @@ -0,0 +1,20 @@ +#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; diff --git a/sc4pd/pd/allpassl.pd b/sc4pd/pd/allpassl.pd new file mode 100644 index 0000000..e2f3251 --- /dev/null +++ b/sc4pd/pd/allpassl.pd @@ -0,0 +1,20 @@ +#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; diff --git a/sc4pd/pd/allpassn.pd b/sc4pd/pd/allpassn.pd new file mode 100644 index 0000000..c0406b9 --- /dev/null +++ b/sc4pd/pd/allpassn.pd @@ -0,0 +1,20 @@ +#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; diff --git a/sc4pd/source/AllpassC.cpp b/sc4pd/source/AllpassC.cpp new file mode 100644 index 0000000..72c05f5 --- /dev/null +++ b/sc4pd/source/AllpassC.cpp @@ -0,0 +1,315 @@ +/* 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 ? */ diff --git a/sc4pd/source/AllpassL.cpp b/sc4pd/source/AllpassL.cpp new file mode 100644 index 0000000..2353df3 --- /dev/null +++ b/sc4pd/source/AllpassL.cpp @@ -0,0 +1,276 @@ +/* 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 ? */ diff --git a/sc4pd/source/AllpassN.cpp b/sc4pd/source/AllpassN.cpp new file mode 100644 index 0000000..a03130b --- /dev/null +++ b/sc4pd/source/AllpassN.cpp @@ -0,0 +1,348 @@ +/* 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 ? */ diff --git a/sc4pd/source/main.cpp b/sc4pd/source/main.cpp index 79e67a9..77a620c 100644 --- a/sc4pd/source/main.cpp +++ b/sc4pd/source/main.cpp @@ -60,7 +60,10 @@ void sc4pd_library_setup() " sqrsum(~), absdif(~), LFSaw(~), LFPulse(~), Impulse(~),\n" " Integrator(~), Decay~, Decay2~, Lag~, Lag2~, LinExp(~), " "DelayN~,\n" - " DelayL~, DelayC~, CombN~, CombL~" + " DelayL~, DelayC~, CombN~, CombL~, CombC~, AllpassN~, " + "AllpassL~,\n" + " AllpassC~" + "\n" ); //initialize objects @@ -217,6 +220,14 @@ void sc4pd_library_setup() FLEXT_DSP_SETUP(CombN_ar); 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); } FLEXT_LIB_SETUP(sc4pd,sc4pd_library_setup); -- cgit v1.2.1