some filters

svn path=/trunk/externals/tb/; revision=1921

6 files changed, 824 insertions, 0 deletions

diff --git a/sc4pd/source/BPZ2.cpp b/sc4pd/source/BPZ2.cpp
new file mode 100644
index 0000000..a18a9dc
--- /dev/null
+++ b/sc4pd/source/BPZ2.cpp
@@ -0,0 +1,140 @@
+/* sc4pd 
+   BPZ2~
+
+   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: William Parker: Compassion Seizes Bed-Stuy
+
+*/
+
+#include "sc4pd.hpp"
+
+/* ------------------------ BPZ2~ -------------------------------*/
+
+class BPZ2_ar:
+    public sc4pd_dsp
+{
+    FLEXT_HEADER(BPZ2_ar,sc4pd_dsp);
+    
+public:
+    BPZ2_ar(int argc, t_atom *argv);
+    
+protected:
+    virtual void m_signal(int n, t_sample *const *in, t_sample *const *out);
+    virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out)
+    {
+	mFilterLoops=sc_filterloops();
+	mFilterRemain=sc_filterremain();
+    }
+
+private:
+    float m_x1, m_x2;
+    int mFilterLoops, mFilterRemain;
+};
+
+FLEXT_LIB_DSP_V("BPZ2~",BPZ2_ar);
+
+BPZ2_ar::BPZ2_ar(int argc, t_atom *argv)
+{
+    AddOutSignal();
+
+    m_x1=0;
+    m_x2=0;
+}
+
+void BPZ2_ar::m_signal(int n, t_sample *const *in, 
+			  t_sample *const *out)
+{
+    t_sample *nin = *in;
+    t_sample *nout = *out;
+
+    float x0;
+    float x1 = m_x1;
+    float x2 = m_x2;
+
+    for (int i = 0; i!=mFilterLoops ;++i)
+    {
+	x0 = ZXP(nin);
+	ZXP(nout) = (x0 - x2) * 0.5f;
+	x2 = ZXP(nin);
+	ZXP(nout) = (x2 - x1) * 0.5f;
+	x1 = ZXP(nin);
+	ZXP(nout) = (x1 - x0) * 0.5f;
+    }
+    
+    for (int i = 0; i!= mFilterRemain;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = (x0 - x2) * 0.5f;
+	x2 = x1;
+	x1 = x0;
+    }
+    m_x1 = x1;
+    m_x2 = x2;
+}
+
+/* ------------------------ BPZ2 -------------------------------*/
+
+class BPZ2_kr:
+    public flext_base
+{
+    FLEXT_HEADER(BPZ2_kr,flext_base);
+    
+public:
+    BPZ2_kr(int argc, t_atom *argv);
+    
+protected:
+    void m_perform (float f);
+    
+private:
+    float m_x1;
+    float m_x2;
+    FLEXT_CALLBACK_F(m_perform);
+};
+
+FLEXT_LIB_V("BPZ2",BPZ2_kr);
+
+BPZ2_kr::BPZ2_kr(int argc, t_atom *argv)
+{
+    FLEXT_ADDMETHOD(0,m_perform);
+
+    AddOutFloat();
+
+    m_x2 = m_x1 = 0;
+}
+
+void BPZ2_kr::m_perform(float f)
+{
+    ToOutFloat(0,(f - m_x2) * 0.5f);
+    m_x2=m_x1;
+    m_x1=f;
+}
+
diff --git a/sc4pd/source/BRZ2.cpp b/sc4pd/source/BRZ2.cpp
new file mode 100644
index 0000000..4c14fa2
--- /dev/null
+++ b/sc4pd/source/BRZ2.cpp
@@ -0,0 +1,140 @@
+/* sc4pd 
+   BRZ2~
+
+   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: William Parker: Compassion Seizes Bed-Stuy
+
+*/
+
+#include "sc4pd.hpp"
+
+/* ------------------------ BRZ2~ -------------------------------*/
+
+class BRZ2_ar:
+    public sc4pd_dsp
+{
+    FLEXT_HEADER(BRZ2_ar,sc4pd_dsp);
+    
+public:
+    BRZ2_ar(int argc, t_atom *argv);
+    
+protected:
+    virtual void m_signal(int n, t_sample *const *in, t_sample *const *out);
+    virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out)
+    {
+	mFilterLoops=sc_filterloops();
+	mFilterRemain=sc_filterremain();
+    }
+
+private:
+    float m_x1, m_x2;
+    int mFilterLoops, mFilterRemain;
+};
+
+FLEXT_LIB_DSP_V("BRZ2~",BRZ2_ar);
+
+BRZ2_ar::BRZ2_ar(int argc, t_atom *argv)
+{
+    AddOutSignal();
+
+    m_x1=0;
+    m_x2=0;
+}
+
+void BRZ2_ar::m_signal(int n, t_sample *const *in, 
+			  t_sample *const *out)
+{
+    t_sample *nin = *in;
+    t_sample *nout = *out;
+
+    float x0;
+    float x1 = m_x1;
+    float x2 = m_x2;
+
+    for (int i = 0; i!=mFilterLoops ;++i)
+    {
+	x0 = ZXP(nin);
+	ZXP(nout) = (x0 + x2) * 0.5f;
+	x2 = ZXP(nin);
+	ZXP(nout) = (x2 + x1) * 0.5f;
+	x1 = ZXP(nin);
+	ZXP(nout) = (x1 + x0) * 0.5f;
+    }
+    
+    for (int i = 0; i!= mFilterRemain;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = (x0 + x2) * 0.5f;
+	x2 = x1;
+	x1 = x0;
+    }
+    m_x1 = x1;
+    m_x2 = x2;
+}
+
+/* ------------------------ BRZ2 -------------------------------*/
+
+class BRZ2_kr:
+    public flext_base
+{
+    FLEXT_HEADER(BRZ2_kr,flext_base);
+    
+public:
+    BRZ2_kr(int argc, t_atom *argv);
+    
+protected:
+    void m_perform (float f);
+    
+private:
+    float m_x1;
+    float m_x2;
+    FLEXT_CALLBACK_F(m_perform);
+};
+
+FLEXT_LIB_V("BRZ2",BRZ2_kr);
+
+BRZ2_kr::BRZ2_kr(int argc, t_atom *argv)
+{
+    FLEXT_ADDMETHOD(0,m_perform);
+
+    AddOutFloat();
+
+    m_x2 = m_x1 = 0;
+}
+
+void BRZ2_kr::m_perform(float f)
+{
+    ToOutFloat(0,(f + m_x2) * 0.5f);
+    m_x2=m_x1;
+    m_x1=f;
+}
+
diff --git a/sc4pd/source/HPZ1.cpp b/sc4pd/source/HPZ1.cpp
new file mode 100644
index 0000000..0ab5f91
--- /dev/null
+++ b/sc4pd/source/HPZ1.cpp
@@ -0,0 +1,132 @@
+/* sc4pd 
+   HPZ1~
+
+   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: William Parker: Compassion Seizes Bed-Stuy
+
+*/
+
+#include "sc4pd.hpp"
+
+/* ------------------------ HPZ1~ -------------------------------*/
+
+class HPZ1_ar:
+    public sc4pd_dsp
+{
+    FLEXT_HEADER(HPZ1_ar,sc4pd_dsp);
+    
+public:
+    HPZ1_ar(int argc, t_atom *argv);
+    
+protected:
+    virtual void m_signal(int n, t_sample *const *in, t_sample *const *out);
+
+private:
+    float m_x1;
+};
+
+FLEXT_LIB_DSP_V("HPZ1~",HPZ1_ar);
+
+HPZ1_ar::HPZ1_ar(int argc, t_atom *argv)
+{
+    AddOutSignal();
+
+    m_x1=0;
+}
+
+void HPZ1_ar::m_signal(int n, t_sample *const *in, 
+			  t_sample *const *out)
+{
+    t_sample *nin = *in;
+    t_sample *nout = *out;
+
+    float x0;
+    float x1 = m_x1;
+
+    int t = n >> 2;
+    for (int i = 0; i!= t;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x0 - x1);
+	x1 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x1 - x0);
+	x0 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x0 - x1);
+	x1 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x1 - x0);
+    }
+    
+    t = n & 3;
+    for (int i = 0; i!= t;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x0 - x1);
+	x1 = x0;
+    }
+    m_x1 = x1;
+}
+
+/* ------------------------ HPZ1 -------------------------------*/
+
+class HPZ1_kr:
+    public flext_base
+{
+    FLEXT_HEADER(HPZ1_kr,flext_base);
+    
+public:
+    HPZ1_kr(int argc, t_atom *argv);
+    
+protected:
+    void m_perform (float f);
+    
+private:
+    float m_x1;
+    FLEXT_CALLBACK_F(m_perform);
+};
+
+FLEXT_LIB_V("HPZ1",HPZ1_kr);
+
+HPZ1_kr::HPZ1_kr(int argc, t_atom *argv)
+{
+    FLEXT_ADDMETHOD(0,m_perform);
+
+    AddOutFloat();
+
+    m_x1=0;
+}
+
+void HPZ1_kr::m_perform(float f)
+{
+    ToOutFloat(0,0.5f * (f - m_x1));
+    m_x1=f;
+}
+
diff --git a/sc4pd/source/HPZ2.cpp b/sc4pd/source/HPZ2.cpp
new file mode 100644
index 0000000..c561455
--- /dev/null
+++ b/sc4pd/source/HPZ2.cpp
@@ -0,0 +1,140 @@
+/* sc4pd 
+   HPZ2~
+
+   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: William Parker: Compassion Seizes Bed-Stuy
+
+*/
+
+#include "sc4pd.hpp"
+
+/* ------------------------ HPZ2~ -------------------------------*/
+
+class HPZ2_ar:
+    public sc4pd_dsp
+{
+    FLEXT_HEADER(HPZ2_ar,sc4pd_dsp);
+    
+public:
+    HPZ2_ar(int argc, t_atom *argv);
+    
+protected:
+    virtual void m_signal(int n, t_sample *const *in, t_sample *const *out);
+    virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out)
+    {
+	mFilterLoops=sc_filterloops();
+	mFilterRemain=sc_filterremain();
+    }
+
+private:
+    float m_x1, m_x2;
+    int mFilterLoops, mFilterRemain;
+};
+
+FLEXT_LIB_DSP_V("HPZ2~",HPZ2_ar);
+
+HPZ2_ar::HPZ2_ar(int argc, t_atom *argv)
+{
+    AddOutSignal();
+
+    m_x1=0;
+    m_x2=0;
+}
+
+void HPZ2_ar::m_signal(int n, t_sample *const *in, 
+			  t_sample *const *out)
+{
+    t_sample *nin = *in;
+    t_sample *nout = *out;
+
+    float x0;
+    float x1 = m_x1;
+    float x2 = m_x2;
+
+    for (int i = 0; i!=mFilterLoops ;++i)
+    {
+	x0 = ZXP(nin);
+	ZXP(nout) = (x0 - 2.f * x1 + x2) * 0.25f;
+	x2 = ZXP(nin);
+	ZXP(nout) = (x2 - 2.f * x0 + x1) * 0.25f;
+	x1 = ZXP(nin);
+	ZXP(nout) = (x1 - 2.f * x2 + x0) * 0.25f;
+    }
+    
+    for (int i = 0; i!= mFilterRemain;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = (x0 - 2.f * x1 + x2) * 0.25f;
+	x2 = x1;
+	x1 = x0;
+    }
+    m_x1 = x1;
+    m_x2 = x2;
+}
+
+/* ------------------------ HPZ2 -------------------------------*/
+
+class HPZ2_kr:
+    public flext_base
+{
+    FLEXT_HEADER(HPZ2_kr,flext_base);
+    
+public:
+    HPZ2_kr(int argc, t_atom *argv);
+    
+protected:
+    void m_perform (float f);
+    
+private:
+    float m_x1;
+    float m_x2;
+    FLEXT_CALLBACK_F(m_perform);
+};
+
+FLEXT_LIB_V("HPZ2",HPZ2_kr);
+
+HPZ2_kr::HPZ2_kr(int argc, t_atom *argv)
+{
+    FLEXT_ADDMETHOD(0,m_perform);
+
+    AddOutFloat();
+
+    m_x2 = m_x1 = 0;
+}
+
+void HPZ2_kr::m_perform(float f)
+{
+    ToOutFloat(0,(f - 2.f * m_x1 + m_x2) * 0.25f);
+    m_x2=m_x1;
+    m_x1=f;
+}
+
diff --git a/sc4pd/source/LPZ1.cpp b/sc4pd/source/LPZ1.cpp
new file mode 100644
index 0000000..70e495b
--- /dev/null
+++ b/sc4pd/source/LPZ1.cpp
@@ -0,0 +1,132 @@
+/* sc4pd 
+   LPZ1~
+
+   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: William Parker: Compassion Seizes Bed-Stuy
+
+*/
+
+#include "sc4pd.hpp"
+
+/* ------------------------ LPZ1~ -------------------------------*/
+
+class LPZ1_ar:
+    public sc4pd_dsp
+{
+    FLEXT_HEADER(LPZ1_ar,sc4pd_dsp);
+    
+public:
+    LPZ1_ar(int argc, t_atom *argv);
+    
+protected:
+    virtual void m_signal(int n, t_sample *const *in, t_sample *const *out);
+
+private:
+    float m_x1;
+};
+
+FLEXT_LIB_DSP_V("LPZ1~",LPZ1_ar);
+
+LPZ1_ar::LPZ1_ar(int argc, t_atom *argv)
+{
+    AddOutSignal();
+
+    m_x1=0;
+}
+
+void LPZ1_ar::m_signal(int n, t_sample *const *in, 
+			  t_sample *const *out)
+{
+    t_sample *nin = *in;
+    t_sample *nout = *out;
+
+    float x0;
+    float x1 = m_x1;
+
+    int t = n >> 2;
+    for (int i = 0; i!= t;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x0 + x1);
+	x1 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x1 + x0);
+	x0 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x0 + x1);
+	x1 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x1 + x0);
+    }
+    
+    t = n & 3;
+    for (int i = 0; i!= t;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = 0.5f * (x0 + x1);
+	x1 = x0;
+    }
+    m_x1 = x1;
+}
+
+/* ------------------------ LPZ1 -------------------------------*/
+
+class LPZ1_kr:
+    public flext_base
+{
+    FLEXT_HEADER(LPZ1_kr,flext_base);
+    
+public:
+    LPZ1_kr(int argc, t_atom *argv);
+    
+protected:
+    void m_perform (float f);
+    
+private:
+    float m_x1;
+    FLEXT_CALLBACK_F(m_perform);
+};
+
+FLEXT_LIB_V("LPZ1",LPZ1_kr);
+
+LPZ1_kr::LPZ1_kr(int argc, t_atom *argv)
+{
+    FLEXT_ADDMETHOD(0,m_perform);
+
+    AddOutFloat();
+
+    m_x1=0;
+}
+
+void LPZ1_kr::m_perform(float f)
+{
+    ToOutFloat(0,0.5f * (f + m_x1));
+    m_x1=f;
+}
+
diff --git a/sc4pd/source/LPZ2.cpp b/sc4pd/source/LPZ2.cpp
new file mode 100644
index 0000000..4165690
--- /dev/null
+++ b/sc4pd/source/LPZ2.cpp
@@ -0,0 +1,140 @@
+/* sc4pd 
+   LPZ2~
+
+   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: William Parker: Compassion Seizes Bed-Stuy
+
+*/
+
+#include "sc4pd.hpp"
+
+/* ------------------------ LPZ2~ -------------------------------*/
+
+class LPZ2_ar:
+    public sc4pd_dsp
+{
+    FLEXT_HEADER(LPZ2_ar,sc4pd_dsp);
+    
+public:
+    LPZ2_ar(int argc, t_atom *argv);
+    
+protected:
+    virtual void m_signal(int n, t_sample *const *in, t_sample *const *out);
+    virtual void m_dsp(int n, t_sample *const *in, t_sample *const *out)
+    {
+	mFilterLoops=sc_filterloops();
+	mFilterRemain=sc_filterremain();
+    }
+
+private:
+    float m_x1, m_x2;
+    int mFilterLoops, mFilterRemain;
+};
+
+FLEXT_LIB_DSP_V("LPZ2~",LPZ2_ar);
+
+LPZ2_ar::LPZ2_ar(int argc, t_atom *argv)
+{
+    AddOutSignal();
+
+    m_x1=0;
+    m_x2=0;
+}
+
+void LPZ2_ar::m_signal(int n, t_sample *const *in, 
+			  t_sample *const *out)
+{
+    t_sample *nin = *in;
+    t_sample *nout = *out;
+
+    float x0;
+    float x1 = m_x1;
+    float x2 = m_x2;
+
+    for (int i = 0; i!=mFilterLoops ;++i)
+    {
+	x0 = ZXP(nin);
+	ZXP(nout) = (x0 + 2.f * x1 + x2) * 0.25f;
+	x2 = ZXP(nin);
+	ZXP(nout) = (x2 + 2.f * x0 + x1) * 0.25f;
+	x1 = ZXP(nin);
+	ZXP(nout) = (x1 + 2.f * x2 + x0) * 0.25f;
+    }
+    
+    for (int i = 0; i!= mFilterRemain;++i)
+    {
+	x0 = ZXP(nin); 
+	ZXP(nout) = (x0 + 2.f * x1 + x2) * 0.25f;
+	x2 = x1;
+	x1 = x0;
+    }
+    m_x1 = x1;
+    m_x2 = x2;
+}
+
+/* ------------------------ LPZ2 -------------------------------*/
+
+class LPZ2_kr:
+    public flext_base
+{
+    FLEXT_HEADER(LPZ2_kr,flext_base);
+    
+public:
+    LPZ2_kr(int argc, t_atom *argv);
+    
+protected:
+    void m_perform (float f);
+    
+private:
+    float m_x1;
+    float m_x2;
+    FLEXT_CALLBACK_F(m_perform);
+};
+
+FLEXT_LIB_V("LPZ2",LPZ2_kr);
+
+LPZ2_kr::LPZ2_kr(int argc, t_atom *argv)
+{
+    FLEXT_ADDMETHOD(0,m_perform);
+
+    AddOutFloat();
+
+    m_x2 = m_x1 = 0;
+}
+
+void LPZ2_kr::m_perform(float f)
+{
+    ToOutFloat(0,(f + 2.f * m_x1 + m_x2) * 0.25f);
+    m_x2=m_x1;
+    m_x1=f;
+}
+