faust
faust copied to clipboard
Published
20 hours ago •
grame-cncm
grame-cncm
DSP on the Teensy with Faust not working when FaustZita object connected
Hello, we are a media art team working with Teensy and various sensors.
We have been using custom made headphones on which Teensy MicroMod (and various sensors) is embedded, so that output of Teensy is directly connected to the headphones' speakers. You may want to check out the artist website: http://www.byungjun.pe.kr/
For the next project we are trying to use Faust to create real-time audio signal processing program that can be uploaded to Teensy. Eventually we will try to get gyro and acceleration data from a sensor and use the data as serial inputs for the Teensy program to make interactive sound installations.
We read Faust homepage's introductory web document ( https://faustdoc.grame.fr/tutorials/teensy/ ) and the linked paper <Real Time Audio Digital Signal Processing with Faust and the Teensy> ( https://zenodo.org/record/3249282#.YdzaqtOA4-S ).
The <Band-Limited Sawtooth Oscillator on the Teensy> tutorial worked just fine with the source code provided ( https://faustdoc.grame.fr/tutorials/teensy/misc/teensy.zip ), and by following the provided procedure using the Faust online compilation service on the Faust Web IDE to generate .cpp and .h file from .dsp file and uploading from Teensyduino appropriate .ino code with the generated files in the same directory.
However, the <Audio Effect> tutorial does not work, either with the source code provided or following the provided procedure.
We have tried increased memory allocation size ( AudioMemory(100); ) and alternative connections ( connecting pink noise synth to faustZita as input rather than an AudioInputI2S object ), but including a FaustZita object always results in no sound at all whereas everything works fine if a FaustZita object is not included in the patchcord connection (declaring a FaustZita object is just fine as long it is not included in the patchcord connection).
For example, the following .ino code works, making pink noise only on the left channel.
#include <Audio.h>
#include "FaustZita.h"
FaustZita faustZita;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioSynthNoisePink pink;
AudioConnection patchCord0(pink, 0, out, 0);
void setup() {
AudioMemory(100);
audioShield.enable();
pink.amplitude(1);
audioShield.volume(0.1);
faustZita.setParamValue("level", 0);
}
void loop() {
}
But the following code does not work, making no sound.
#include <Audio.h>
#include "FaustZita.h"
FaustZita faustZita;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioSynthNoisePink pink;
AudioConnection patchCord0(pink, 0, faustZita, 0);
AudioConnection patchCord1(faustZita, 0, out, 0);
void setup() {
AudioMemory(100);
audioShield.enable();
pink.amplitude(1);
audioShield.volume(0.1);
faustZita.setParamValue("level", 0);
}
void loop() {
}
We are using Teensy 4.0 MicroMod and Teensyduino, and have been trying with several different computers running on Linux and Mac.
There is no problem with compiling and uploading to Teensy, but when we check the Serial Monitor from Teensyduino, the Teensy seems to be keep rebooting itself ( the top bar text turns to "[offline] (Teensy)" from "[/dev/cu.usbmodem84929401] (Teensy)" for approximately 1 second ), once every approximately 8 seconds. There is no sound at all, other than background electric noise.
The following are FaustZita.cpp and FaustZita.h, just in case.
// FaustZita.h
#ifndef faust_FaustZita_h_
#define faust_FaustZita_h_
#include <string>
#include "Arduino.h"
#include "AudioStream.h"
#include "Audio.h"
class dsp;
class MapUI;
class FaustZita : public AudioStream
{
public:
FaustZita();
~FaustZita();
virtual void update(void);
void setParamValue(const std::string& path, float value);
private:
float** fInChannel;
float** fOutChannel;
MapUI* UI;
dsp* fDSP;
};
#endif
// FaustZita.cpp
/* ------------------------------------------------------------
name: "FaustZita"
Code generated with Faust 2.15.5 (https://faust.grame.fr)
Compilation options: cpp, -scal -ftz 0
------------------------------------------------------------ */
#ifndef __mydsp_H__
#define __mydsp_H__
#include "FaustZita.h"
// IMPORTANT: in order for MapUI to work, the teensy linker must be g++
/************************************************************************
FAUST Architecture File
Copyright (C) 2003-2017 GRAME, Centre National de Creation Musicale
---------------------------------------------------------------------
This Architecture section 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 3 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, see <http://www.gnu.org/licenses/>.
EXCEPTION : As a special exception, you may create a larger work
that contains this FAUST architecture section and distribute
that work under terms of your choice, so long as this FAUST
architecture section is not modified.
************************************************************************/
#ifndef FAUST_MAPUI_H
#define FAUST_MAPUI_H
#include <vector>
#include <map>
#include <string>
/************************************************************************
FAUST Architecture File
Copyright (C) 2003-2017 GRAME, Centre National de Creation Musicale
---------------------------------------------------------------------
This Architecture section 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 3 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, see <http://www.gnu.org/licenses/>.
EXCEPTION : As a special exception, you may create a larger work
that contains this FAUST architecture section and distribute
that work under terms of your choice, so long as this FAUST
architecture section is not modified.
************************************************************************/
#ifndef __UI_H__
#define __UI_H__
#ifndef FAUSTFLOAT
#define FAUSTFLOAT float
#endif
/*******************************************************************************
* UI : Faust DSP User Interface
* User Interface as expected by the buildUserInterface() method of a DSP.
* This abstract class contains only the method that the Faust compiler can
* generate to describe a DSP user interface.
******************************************************************************/
struct Soundfile;
class UI
{
public:
UI() {}
virtual ~UI() {}
// -- widget's layouts
virtual void openTabBox(const char* label) = 0;
virtual void openHorizontalBox(const char* label) = 0;
virtual void openVerticalBox(const char* label) = 0;
virtual void closeBox() = 0;
// -- active widgets
virtual void addButton(const char* label, FAUSTFLOAT* zone) = 0;
virtual void addCheckButton(const char* label, FAUSTFLOAT* zone) = 0;
virtual void addVerticalSlider(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT init, FAUSTFLOAT min, FAUSTFLOAT max, FAUSTFLOAT step) = 0;
virtual void addHorizontalSlider(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT init, FAUSTFLOAT min, FAUSTFLOAT max, FAUSTFLOAT step) = 0;
virtual void addNumEntry(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT init, FAUSTFLOAT min, FAUSTFLOAT max, FAUSTFLOAT step) = 0;
// -- passive widgets
virtual void addHorizontalBargraph(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT min, FAUSTFLOAT max) = 0;
virtual void addVerticalBargraph(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT min, FAUSTFLOAT max) = 0;
// -- soundfiles
virtual void addSoundfile(const char* label, const char* filename, Soundfile** sf_zone) = 0;
// -- metadata declarations
virtual void declare(FAUSTFLOAT*, const char*, const char*) {}
};
#endif
/************************************************************************
FAUST Architecture File
Copyright (C) 2003-2017 GRAME, Centre National de Creation Musicale
---------------------------------------------------------------------
This Architecture section 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 3 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, see <http://www.gnu.org/licenses/>.
EXCEPTION : As a special exception, you may create a larger work
that contains this FAUST architecture section and distribute
that work under terms of your choice, so long as this FAUST
architecture section is not modified.
************************************************************************/
#ifndef FAUST_PATHBUILDER_H
#define FAUST_PATHBUILDER_H
#include <vector>
#include <string>
#include <algorithm>
/*******************************************************************************
* PathBuilder : Faust User Interface
* Helper class to build complete hierarchical path for UI items.
******************************************************************************/
class PathBuilder
{
protected:
std::vector<std::string> fControlsLevel;
public:
PathBuilder() {}
virtual ~PathBuilder() {}
std::string buildPath(const std::string& label)
{
std::string res = "/";
for (size_t i = 0; i < fControlsLevel.size(); i++) {
res += fControlsLevel[i];
res += "/";
}
res += label;
std::replace(res.begin(), res.end(), ' ', '_');
return res;
}
void pushLabel(const std::string& label) { fControlsLevel.push_back(label); }
void popLabel() { fControlsLevel.pop_back(); }
};
#endif // FAUST_PATHBUILDER_H
/*******************************************************************************
* MapUI : Faust User Interface
* This class creates a map of complete hierarchical path and zones for each UI items.
******************************************************************************/
class MapUI : public UI, public PathBuilder
{
protected:
// Complete path map
std::map<std::string, FAUSTFLOAT*> fPathZoneMap;
// Label zone map
std::map<std::string, FAUSTFLOAT*> fLabelZoneMap;
public:
MapUI() {};
virtual ~MapUI() {};
// -- widget's layouts
void openTabBox(const char* label)
{
fControlsLevel.push_back(label);
}
void openHorizontalBox(const char* label)
{
fControlsLevel.push_back(label);
}
void openVerticalBox(const char* label)
{
fControlsLevel.push_back(label);
}
void closeBox()
{
fControlsLevel.pop_back();
}
// -- active widgets
void addButton(const char* label, FAUSTFLOAT* zone)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
void addCheckButton(const char* label, FAUSTFLOAT* zone)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
void addVerticalSlider(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT init, FAUSTFLOAT fmin, FAUSTFLOAT fmax, FAUSTFLOAT step)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
void addHorizontalSlider(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT init, FAUSTFLOAT fmin, FAUSTFLOAT fmax, FAUSTFLOAT step)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
void addNumEntry(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT init, FAUSTFLOAT fmin, FAUSTFLOAT fmax, FAUSTFLOAT step)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
// -- passive widgets
void addHorizontalBargraph(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT fmin, FAUSTFLOAT fmax)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
void addVerticalBargraph(const char* label, FAUSTFLOAT* zone, FAUSTFLOAT fmin, FAUSTFLOAT fmax)
{
fPathZoneMap[buildPath(label)] = zone;
fLabelZoneMap[label] = zone;
}
// -- soundfiles
virtual void addSoundfile(const char* label, const char* filename, Soundfile** sf_zone) {}
// -- metadata declarations
void declare(FAUSTFLOAT* zone, const char* key, const char* val)
{}
// set/get
void setParamValue(const std::string& path, FAUSTFLOAT value)
{
if (fPathZoneMap.find(path) != fPathZoneMap.end()) {
*fPathZoneMap[path] = value;
} else if (fLabelZoneMap.find(path) != fLabelZoneMap.end()) {
*fLabelZoneMap[path] = value;
}
}
FAUSTFLOAT getParamValue(const std::string& path)
{
if (fPathZoneMap.find(path) != fPathZoneMap.end()) {
return *fPathZoneMap[path];
} else if (fLabelZoneMap.find(path) != fLabelZoneMap.end()) {
return *fLabelZoneMap[path];
} else {
return FAUSTFLOAT(0);
}
}
// map access
std::map<std::string, FAUSTFLOAT*>& getMap() { return fPathZoneMap; }
int getParamsCount() { return int(fPathZoneMap.size()); }
std::string getParamAddress(int index)
{
std::map<std::string, FAUSTFLOAT*>::iterator it = fPathZoneMap.begin();
while (index-- > 0 && it++ != fPathZoneMap.end()) {}
return (*it).first;
}
std::string getParamAddress(FAUSTFLOAT* zone)
{
std::map<std::string, FAUSTFLOAT*>::iterator it = fPathZoneMap.begin();
do {
if ((*it).second == zone) return (*it).first;
}
while (it++ != fPathZoneMap.end());
return "";
}
};
#endif // FAUST_MAPUI_H
/************************************************************************
FAUST Architecture File
Copyright (C) 2003-2017 GRAME, Centre National de Creation Musicale
---------------------------------------------------------------------
This Architecture section 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 3 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, see <http://www.gnu.org/licenses/>.
EXCEPTION : As a special exception, you may create a larger work
that contains this FAUST architecture section and distribute
that work under terms of your choice, so long as this FAUST
architecture section is not modified.
************************************************************************/
#ifndef __meta__
#define __meta__
struct Meta
{
virtual void declare(const char* key, const char* value) = 0;
virtual ~Meta() {};
};
#endif
/************************************************************************
FAUST Architecture File
Copyright (C) 2003-2017 GRAME, Centre National de Creation Musicale
---------------------------------------------------------------------
This Architecture section 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 3 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, see <http://www.gnu.org/licenses/>.
EXCEPTION : As a special exception, you may create a larger work
that contains this FAUST architecture section and distribute
that work under terms of your choice, so long as this FAUST
architecture section is not modified.
************************************************************************/
#ifndef __dsp__
#define __dsp__
#include <string>
#include <vector>
#ifndef FAUSTFLOAT
#define FAUSTFLOAT float
#endif
class UI;
struct Meta;
/**
* DSP memory manager.
*/
struct dsp_memory_manager {
virtual ~dsp_memory_manager() {}
virtual void* allocate(size_t size) = 0;
virtual void destroy(void* ptr) = 0;
};
/**
* Signal processor definition.
*/
class dsp {
public:
dsp() {}
virtual ~dsp() {}
/* Return instance number of audio inputs */
virtual int getNumInputs() = 0;
/* Return instance number of audio outputs */
virtual int getNumOutputs() = 0;
/**
* Trigger the ui_interface parameter with instance specific calls
* to 'addBtton', 'addVerticalSlider'... in order to build the UI.
*
* @param ui_interface - the user interface builder
*/
virtual void buildUserInterface(UI* ui_interface) = 0;
/* Returns the sample rate currently used by the instance */
virtual int getSampleRate() = 0;
/**
* Global init, calls the following methods:
* - static class 'classInit': static tables initialization
* - 'instanceInit': constants and instance state initialization
*
* @param samplingRate - the sampling rate in Hertz
*/
virtual void init(int samplingRate) = 0;
/**
* Init instance state
*
* @param samplingRate - the sampling rate in Hertz
*/
virtual void instanceInit(int samplingRate) = 0;
/**
* Init instance constant state
*
* @param samplingRate - the sampling rate in Hertz
*/
virtual void instanceConstants(int samplingRate) = 0;
/* Init default control parameters values */
virtual void instanceResetUserInterface() = 0;
/* Init instance state (delay lines...) */
virtual void instanceClear() = 0;
/**
* Return a clone of the instance.
*
* @return a copy of the instance on success, otherwise a null pointer.
*/
virtual dsp* clone() = 0;
/**
* Trigger the Meta* parameter with instance specific calls to 'declare' (key, value) metadata.
*
* @param m - the Meta* meta user
*/
virtual void metadata(Meta* m) = 0;
/**
* DSP instance computation, to be called with successive in/out audio buffers.
*
* @param count - the number of frames to compute
* @param inputs - the input audio buffers as an array of non-interleaved FAUSTFLOAT samples (eiher float, double or quad)
* @param outputs - the output audio buffers as an array of non-interleaved FAUSTFLOAT samples (eiher float, double or quad)
*
*/
virtual void compute(int count, FAUSTFLOAT** inputs, FAUSTFLOAT** outputs) = 0;
/**
* DSP instance computation: alternative method to be used by subclasses.
*
* @param date_usec - the timestamp in microsec given by audio driver.
* @param count - the number of frames to compute
* @param inputs - the input audio buffers as an array of non-interleaved FAUSTFLOAT samples (eiher float, double or quad)
* @param outputs - the output audio buffers as an array of non-interleaved FAUSTFLOAT samples (eiher float, double or quad)
*
*/
virtual void compute(double /*date_usec*/, int count, FAUSTFLOAT** inputs, FAUSTFLOAT** outputs) { compute(count, inputs, outputs); }
};
/**
* Generic DSP decorator.
*/
class decorator_dsp : public dsp {
protected:
dsp* fDSP;
public:
decorator_dsp(dsp* dsp = 0):fDSP(dsp) {}
virtual ~decorator_dsp() { delete fDSP; }
virtual int getNumInputs() { return fDSP->getNumInputs(); }
virtual int getNumOutputs() { return fDSP->getNumOutputs(); }
virtual void buildUserInterface(UI* ui_interface) { fDSP->buildUserInterface(ui_interface); }
virtual int getSampleRate() { return fDSP->getSampleRate(); }
virtual void init(int samplingRate) { fDSP->init(samplingRate); }
virtual void instanceInit(int samplingRate) { fDSP->instanceInit(samplingRate); }
virtual void instanceConstants(int samplingRate) { fDSP->instanceConstants(samplingRate); }
virtual void instanceResetUserInterface() { fDSP->instanceResetUserInterface(); }
virtual void instanceClear() { fDSP->instanceClear(); }
virtual decorator_dsp* clone() { return new decorator_dsp(fDSP->clone()); }
virtual void metadata(Meta* m) { fDSP->metadata(m); }
// Beware: subclasses usually have to overload the two 'compute' methods
virtual void compute(int count, FAUSTFLOAT** inputs, FAUSTFLOAT** outputs) { fDSP->compute(count, inputs, outputs); }
virtual void compute(double date_usec, int count, FAUSTFLOAT** inputs, FAUSTFLOAT** outputs) { fDSP->compute(date_usec, count, inputs, outputs); }
};
/**
* DSP factory class.
*/
class dsp_factory {
protected:
// So that to force sub-classes to use deleteDSPFactory(dsp_factory* factory);
virtual ~dsp_factory() {}
public:
virtual std::string getName() = 0;
virtual std::string getSHAKey() = 0;
virtual std::string getDSPCode() = 0;
virtual std::string getCompileOptions() = 0;
virtual std::vector<std::string> getLibraryList() = 0;
virtual std::vector<std::string> getIncludePathnames() = 0;
virtual dsp* createDSPInstance() = 0;
virtual void setMemoryManager(dsp_memory_manager* manager) = 0;
virtual dsp_memory_manager* getMemoryManager() = 0;
};
/**
* On Intel set FZ (Flush to Zero) and DAZ (Denormals Are Zero)
* flags to avoid costly denormals.
*/
#ifdef __SSE__
#include <xmmintrin.h>
#ifdef __SSE2__
#define AVOIDDENORMALS _mm_setcsr(_mm_getcsr() | 0x8040)
#else
#define AVOIDDENORMALS _mm_setcsr(_mm_getcsr() | 0x8000)
#endif
#else
#define AVOIDDENORMALS
#endif
#endif
#ifndef FAUSTFLOAT
#define FAUSTFLOAT float
#endif
#include <algorithm>
#include <cmath>
#include <math.h>
static float mydsp_faustpower2_f(float value) {
return (value * value);
}
#ifndef FAUSTCLASS
#define FAUSTCLASS mydsp
#endif
#ifdef __APPLE__
#define exp10f __exp10f
#define exp10 __exp10
#endif
class mydsp : public dsp {
private:
FAUSTFLOAT fVslider0;
float fRec0[2];
FAUSTFLOAT fVslider1;
float fRec1[2];
int IOTA;
float fVec0[4096];
int fSamplingFreq;
float fConst0;
float fConst1;
float fConst2;
float fConst3;
float fConst4;
float fConst5;
float fConst6;
float fConst7;
float fConst8;
float fConst9;
float fConst10;
float fConst11;
float fConst12;
float fConst13;
float fConst14;
float fConst15;
float fConst16;
float fConst17;
float fConst18;
float fConst19;
float fRec15[2];
float fRec14[2];
float fVec1[2048];
float fConst20;
int iConst21;
float fVec2[4096];
int iConst22;
float fVec3[2048];
int iConst23;
float fRec12[2];
float fConst24;
float fConst25;
float fConst26;
float fConst27;
float fConst28;
float fConst29;
float fConst30;
float fConst31;
float fConst32;
float fConst33;
float fConst34;
float fRec19[2];
float fRec18[2];
float fVec4[2048];
float fConst35;
int iConst36;
float fVec5[4096];
int iConst37;
float fRec16[2];
float fConst38;
float fConst39;
float fConst40;
float fConst41;
float fConst42;
float fConst43;
float fConst44;
float fConst45;
float fConst46;
float fConst47;
float fConst48;
float fRec23[2];
float fRec22[2];
float fVec6[2048];
float fConst49;
int iConst50;
float fVec7[4096];
int iConst51;
float fRec20[2];
float fConst52;
float fConst53;
float fConst54;
float fConst55;
float fConst56;
float fConst57;
float fConst58;
float fConst59;
float fConst60;
float fConst61;
float fConst62;
float fRec27[2];
float fRec26[2];
float fVec8[2048];
float fConst63;
int iConst64;
float fVec9[4096];
int iConst65;
float fRec24[2];
float fConst66;
float fConst67;
float fConst68;
float fConst69;
float fConst70;
float fConst71;
float fConst72;
float fConst73;
float fConst74;
float fConst75;
float fConst76;
float fRec31[2];
float fRec30[2];
float fVec10[2048];
float fConst77;
int iConst78;
float fVec11[2048];
int iConst79;
float fRec28[2];
float fConst80;
float fConst81;
float fConst82;
float fConst83;
float fConst84;
float fConst85;
float fConst86;
float fConst87;
float fConst88;
float fConst89;
float fConst90;
float fRec35[2];
float fRec34[2];
float fVec12[2048];
float fConst91;
int iConst92;
float fVec13[4096];
int iConst93;
float fRec32[2];
float fConst94;
float fConst95;
float fConst96;
float fConst97;
float fConst98;
float fConst99;
float fConst100;
float fConst101;
float fConst102;
float fConst103;
float fConst104;
float fRec39[2];
float fRec38[2];
float fVec14[2048];
float fConst105;
int iConst106;
float fVec15[4096];
int iConst107;
float fRec36[2];
float fConst108;
float fConst109;
float fConst110;
float fConst111;
float fConst112;
float fConst113;
float fConst114;
float fConst115;
float fConst116;
float fConst117;
float fConst118;
float fRec43[2];
float fRec42[2];
float fVec16[2048];
float fConst119;
int iConst120;
float fVec17[2048];
int iConst121;
float fRec40[2];
float fRec4[3];
float fRec5[3];
float fRec6[3];
float fRec7[3];
float fRec8[3];
float fRec9[3];
float fRec10[3];
float fRec11[3];
float fConst122;
float fRec3[3];
float fConst123;
float fRec2[3];
float fRec45[3];
float fRec44[3];
public:
void metadata(Meta* m) {
m->declare("basics.lib/name", "Faust Basic Element Library");
m->declare("basics.lib/version", "0.0");
m->declare("delays.lib/name", "Faust Delay Library");
m->declare("delays.lib/version", "0.1");
m->declare("filename", "FaustZita");
m->declare("filters.lib/name", "Faust Filters Library");
m->declare("filters.lib/version", "0.0");
m->declare("maths.lib/author", "GRAME");
m->declare("maths.lib/copyright", "GRAME");
m->declare("maths.lib/license", "LGPL with exception");
m->declare("maths.lib/name", "Faust Math Library");
m->declare("maths.lib/version", "2.1");
m->declare("name", "FaustZita");
m->declare("routes.lib/name", "Faust Signal Routing Library");
m->declare("routes.lib/version", "0.0");
m->declare("signals.lib/name", "Faust Signal Routing Library");
m->declare("signals.lib/version", "0.0");
}
virtual int getNumInputs() {
return 2;
}
virtual int getNumOutputs() {
return 2;
}
virtual int getInputRate(int channel) {
int rate;
switch (channel) {
case 0: {
rate = 1;
break;
}
case 1: {
rate = 1;
break;
}
default: {
rate = -1;
break;
}
}
return rate;
}
virtual int getOutputRate(int channel) {
int rate;
switch (channel) {
case 0: {
rate = 1;
break;
}
case 1: {
rate = 1;
break;
}
default: {
rate = -1;
break;
}
}
return rate;
}
static void classInit(int samplingFreq) {
}
virtual void instanceConstants(int samplingFreq) {
fSamplingFreq = samplingFreq;
fConst0 = std::min<float>(192000.0f, std::max<float>(1.0f, float(fSamplingFreq)));
fConst1 = (9424.77832f / fConst0);
fConst2 = ((1.0f - fConst1) / (fConst1 + 1.0f));
fConst3 = (1979.20337f / fConst0);
fConst4 = ((1.0f - fConst3) / (fConst3 + 1.0f));
fConst5 = std::cos((37699.1133f / fConst0));
fConst6 = std::floor(((0.219990999f * fConst0) + 0.5f));
fConst7 = ((0.0f - (6.90775537f * fConst6)) / fConst0);
fConst8 = std::exp((0.5f * fConst7));
fConst9 = mydsp_faustpower2_f(fConst8);
fConst10 = (1.0f - (fConst5 * fConst9));
fConst11 = (1.0f - fConst9);
fConst12 = (fConst10 / fConst11);
fConst13 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst10) / mydsp_faustpower2_f(fConst11)) + -1.0f)));
fConst14 = (fConst12 - fConst13);
fConst15 = ((fConst13 + (1.0f - fConst12)) * fConst8);
fConst16 = ((std::exp((0.333333343f * fConst7)) / fConst8) + -1.0f);
fConst17 = (1.0f / std::tan((628.318542f / fConst0)));
fConst18 = (1.0f / (fConst17 + 1.0f));
fConst19 = (1.0f - fConst17);
fConst20 = std::floor(((0.0191229992f * fConst0) + 0.5f));
iConst21 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst6 - fConst20))));
iConst22 = int(std::min<float>(2700.0f, std::max<float>(0.0f, (0.0599999987f * fConst0))));
iConst23 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst20 + -1.0f))));
fConst24 = std::floor(((0.256891012f * fConst0) + 0.5f));
fConst25 = ((0.0f - (6.90775537f * fConst24)) / fConst0);
fConst26 = std::exp((0.5f * fConst25));
fConst27 = mydsp_faustpower2_f(fConst26);
fConst28 = (1.0f - (fConst5 * fConst27));
fConst29 = (1.0f - fConst27);
fConst30 = (fConst28 / fConst29);
fConst31 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst28) / mydsp_faustpower2_f(fConst29)) + -1.0f)));
fConst32 = (fConst30 - fConst31);
fConst33 = ((fConst31 + (1.0f - fConst30)) * fConst26);
fConst34 = ((std::exp((0.333333343f * fConst25)) / fConst26) + -1.0f);
fConst35 = std::floor(((0.0273330007f * fConst0) + 0.5f));
iConst36 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst24 - fConst35))));
iConst37 = int(std::min<float>(2048.0f, std::max<float>(0.0f, (fConst35 + -1.0f))));
fConst38 = std::floor(((0.192303002f * fConst0) + 0.5f));
fConst39 = ((0.0f - (6.90775537f * fConst38)) / fConst0);
fConst40 = std::exp((0.5f * fConst39));
fConst41 = mydsp_faustpower2_f(fConst40);
fConst42 = (1.0f - (fConst5 * fConst41));
fConst43 = (1.0f - fConst41);
fConst44 = (fConst42 / fConst43);
fConst45 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst42) / mydsp_faustpower2_f(fConst43)) + -1.0f)));
fConst46 = (fConst44 - fConst45);
fConst47 = ((fConst45 + (1.0f - fConst44)) * fConst40);
fConst48 = ((std::exp((0.333333343f * fConst39)) / fConst40) + -1.0f);
fConst49 = std::floor(((0.0292910002f * fConst0) + 0.5f));
iConst50 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst38 - fConst49))));
iConst51 = int(std::min<float>(2048.0f, std::max<float>(0.0f, (fConst49 + -1.0f))));
fConst52 = std::floor(((0.210389003f * fConst0) + 0.5f));
fConst53 = ((0.0f - (6.90775537f * fConst52)) / fConst0);
fConst54 = std::exp((0.5f * fConst53));
fConst55 = mydsp_faustpower2_f(fConst54);
fConst56 = (1.0f - (fConst5 * fConst55));
fConst57 = (1.0f - fConst55);
fConst58 = (fConst56 / fConst57);
fConst59 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst56) / mydsp_faustpower2_f(fConst57)) + -1.0f)));
fConst60 = (fConst58 - fConst59);
fConst61 = ((fConst59 + (1.0f - fConst58)) * fConst54);
fConst62 = ((std::exp((0.333333343f * fConst53)) / fConst54) + -1.0f);
fConst63 = std::floor(((0.0244210009f * fConst0) + 0.5f));
iConst64 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst52 - fConst63))));
iConst65 = int(std::min<float>(2048.0f, std::max<float>(0.0f, (fConst63 + -1.0f))));
fConst66 = std::floor(((0.125f * fConst0) + 0.5f));
fConst67 = ((0.0f - (6.90775537f * fConst66)) / fConst0);
fConst68 = std::exp((0.5f * fConst67));
fConst69 = mydsp_faustpower2_f(fConst68);
fConst70 = (1.0f - (fConst5 * fConst69));
fConst71 = (1.0f - fConst69);
fConst72 = (fConst70 / fConst71);
fConst73 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst70) / mydsp_faustpower2_f(fConst71)) + -1.0f)));
fConst74 = (fConst72 - fConst73);
fConst75 = ((fConst73 + (1.0f - fConst72)) * fConst68);
fConst76 = ((std::exp((0.333333343f * fConst67)) / fConst68) + -1.0f);
fConst77 = std::floor(((0.0134579996f * fConst0) + 0.5f));
iConst78 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst66 - fConst77))));
iConst79 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst77 + -1.0f))));
fConst80 = std::floor(((0.127837002f * fConst0) + 0.5f));
fConst81 = ((0.0f - (6.90775537f * fConst80)) / fConst0);
fConst82 = std::exp((0.5f * fConst81));
fConst83 = mydsp_faustpower2_f(fConst82);
fConst84 = (1.0f - (fConst5 * fConst83));
fConst85 = (1.0f - fConst83);
fConst86 = (fConst84 / fConst85);
fConst87 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst84) / mydsp_faustpower2_f(fConst85)) + -1.0f)));
fConst88 = (fConst86 - fConst87);
fConst89 = ((fConst87 + (1.0f - fConst86)) * fConst82);
fConst90 = ((std::exp((0.333333343f * fConst81)) / fConst82) + -1.0f);
fConst91 = std::floor(((0.0316039994f * fConst0) + 0.5f));
iConst92 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst80 - fConst91))));
iConst93 = int(std::min<float>(2048.0f, std::max<float>(0.0f, (fConst91 + -1.0f))));
fConst94 = std::floor(((0.174713001f * fConst0) + 0.5f));
fConst95 = ((0.0f - (6.90775537f * fConst94)) / fConst0);
fConst96 = std::exp((0.5f * fConst95));
fConst97 = mydsp_faustpower2_f(fConst96);
fConst98 = (1.0f - (fConst5 * fConst97));
fConst99 = (1.0f - fConst97);
fConst100 = (fConst98 / fConst99);
fConst101 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst98) / mydsp_faustpower2_f(fConst99)) + -1.0f)));
fConst102 = (fConst100 - fConst101);
fConst103 = ((fConst101 + (1.0f - fConst100)) * fConst96);
fConst104 = ((std::exp((0.333333343f * fConst95)) / fConst96) + -1.0f);
fConst105 = std::floor(((0.0229039993f * fConst0) + 0.5f));
iConst106 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst94 - fConst105))));
iConst107 = int(std::min<float>(2048.0f, std::max<float>(0.0f, (fConst105 + -1.0f))));
fConst108 = std::floor(((0.153128996f * fConst0) + 0.5f));
fConst109 = ((0.0f - (6.90775537f * fConst108)) / fConst0);
fConst110 = std::exp((0.5f * fConst109));
fConst111 = mydsp_faustpower2_f(fConst110);
fConst112 = (1.0f - (fConst5 * fConst111));
fConst113 = (1.0f - fConst111);
fConst114 = (fConst112 / fConst113);
fConst115 = std::sqrt(std::max<float>(0.0f, ((mydsp_faustpower2_f(fConst112) / mydsp_faustpower2_f(fConst113)) + -1.0f)));
fConst116 = (fConst114 - fConst115);
fConst117 = ((fConst115 + (1.0f - fConst114)) * fConst110);
fConst118 = ((std::exp((0.333333343f * fConst109)) / fConst110) + -1.0f);
fConst119 = std::floor(((0.0203460008f * fConst0) + 0.5f));
iConst120 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst108 - fConst119))));
iConst121 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst119 + -1.0f))));
fConst122 = (0.0f - ((fConst4 + 1.0f) * std::cos(fConst3)));
fConst123 = (0.0f - ((fConst2 + 1.0f) * std::cos(fConst1)));
}
virtual void instanceResetUserInterface() {
fVslider0 = FAUSTFLOAT(-6.0f);
fVslider1 = FAUSTFLOAT(0.0f);
}
virtual void instanceClear() {
for (int l0 = 0; (l0 < 2); l0 = (l0 + 1)) {
fRec0[l0] = 0.0f;
}
for (int l1 = 0; (l1 < 2); l1 = (l1 + 1)) {
fRec1[l1] = 0.0f;
}
IOTA = 0;
for (int l2 = 0; (l2 < 4096); l2 = (l2 + 1)) {
fVec0[l2] = 0.0f;
}
for (int l3 = 0; (l3 < 2); l3 = (l3 + 1)) {
fRec15[l3] = 0.0f;
}
for (int l4 = 0; (l4 < 2); l4 = (l4 + 1)) {
fRec14[l4] = 0.0f;
}
for (int l5 = 0; (l5 < 2048); l5 = (l5 + 1)) {
fVec1[l5] = 0.0f;
}
for (int l6 = 0; (l6 < 4096); l6 = (l6 + 1)) {
fVec2[l6] = 0.0f;
}
for (int l7 = 0; (l7 < 2048); l7 = (l7 + 1)) {
fVec3[l7] = 0.0f;
}
for (int l8 = 0; (l8 < 2); l8 = (l8 + 1)) {
fRec12[l8] = 0.0f;
}
for (int l9 = 0; (l9 < 2); l9 = (l9 + 1)) {
fRec19[l9] = 0.0f;
}
for (int l10 = 0; (l10 < 2); l10 = (l10 + 1)) {
fRec18[l10] = 0.0f;
}
for (int l11 = 0; (l11 < 2048); l11 = (l11 + 1)) {
fVec4[l11] = 0.0f;
}
for (int l12 = 0; (l12 < 4096); l12 = (l12 + 1)) {
fVec5[l12] = 0.0f;
}
for (int l13 = 0; (l13 < 2); l13 = (l13 + 1)) {
fRec16[l13] = 0.0f;
}
for (int l14 = 0; (l14 < 2); l14 = (l14 + 1)) {
fRec23[l14] = 0.0f;
}
for (int l15 = 0; (l15 < 2); l15 = (l15 + 1)) {
fRec22[l15] = 0.0f;
}
for (int l16 = 0; (l16 < 2048); l16 = (l16 + 1)) {
fVec6[l16] = 0.0f;
}
for (int l17 = 0; (l17 < 4096); l17 = (l17 + 1)) {
fVec7[l17] = 0.0f;
}
for (int l18 = 0; (l18 < 2); l18 = (l18 + 1)) {
fRec20[l18] = 0.0f;
}
for (int l19 = 0; (l19 < 2); l19 = (l19 + 1)) {
fRec27[l19] = 0.0f;
}
for (int l20 = 0; (l20 < 2); l20 = (l20 + 1)) {
fRec26[l20] = 0.0f;
}
for (int l21 = 0; (l21 < 2048); l21 = (l21 + 1)) {
fVec8[l21] = 0.0f;
}
for (int l22 = 0; (l22 < 4096); l22 = (l22 + 1)) {
fVec9[l22] = 0.0f;
}
for (int l23 = 0; (l23 < 2); l23 = (l23 + 1)) {
fRec24[l23] = 0.0f;
}
for (int l24 = 0; (l24 < 2); l24 = (l24 + 1)) {
fRec31[l24] = 0.0f;
}
for (int l25 = 0; (l25 < 2); l25 = (l25 + 1)) {
fRec30[l25] = 0.0f;
}
for (int l26 = 0; (l26 < 2048); l26 = (l26 + 1)) {
fVec10[l26] = 0.0f;
}
for (int l27 = 0; (l27 < 2048); l27 = (l27 + 1)) {
fVec11[l27] = 0.0f;
}
for (int l28 = 0; (l28 < 2); l28 = (l28 + 1)) {
fRec28[l28] = 0.0f;
}
for (int l29 = 0; (l29 < 2); l29 = (l29 + 1)) {
fRec35[l29] = 0.0f;
}
for (int l30 = 0; (l30 < 2); l30 = (l30 + 1)) {
fRec34[l30] = 0.0f;
}
for (int l31 = 0; (l31 < 2048); l31 = (l31 + 1)) {
fVec12[l31] = 0.0f;
}
for (int l32 = 0; (l32 < 4096); l32 = (l32 + 1)) {
fVec13[l32] = 0.0f;
}
for (int l33 = 0; (l33 < 2); l33 = (l33 + 1)) {
fRec32[l33] = 0.0f;
}
for (int l34 = 0; (l34 < 2); l34 = (l34 + 1)) {
fRec39[l34] = 0.0f;
}
for (int l35 = 0; (l35 < 2); l35 = (l35 + 1)) {
fRec38[l35] = 0.0f;
}
for (int l36 = 0; (l36 < 2048); l36 = (l36 + 1)) {
fVec14[l36] = 0.0f;
}
for (int l37 = 0; (l37 < 4096); l37 = (l37 + 1)) {
fVec15[l37] = 0.0f;
}
for (int l38 = 0; (l38 < 2); l38 = (l38 + 1)) {
fRec36[l38] = 0.0f;
}
for (int l39 = 0; (l39 < 2); l39 = (l39 + 1)) {
fRec43[l39] = 0.0f;
}
for (int l40 = 0; (l40 < 2); l40 = (l40 + 1)) {
fRec42[l40] = 0.0f;
}
for (int l41 = 0; (l41 < 2048); l41 = (l41 + 1)) {
fVec16[l41] = 0.0f;
}
for (int l42 = 0; (l42 < 2048); l42 = (l42 + 1)) {
fVec17[l42] = 0.0f;
}
for (int l43 = 0; (l43 < 2); l43 = (l43 + 1)) {
fRec40[l43] = 0.0f;
}
for (int l44 = 0; (l44 < 3); l44 = (l44 + 1)) {
fRec4[l44] = 0.0f;
}
for (int l45 = 0; (l45 < 3); l45 = (l45 + 1)) {
fRec5[l45] = 0.0f;
}
for (int l46 = 0; (l46 < 3); l46 = (l46 + 1)) {
fRec6[l46] = 0.0f;
}
for (int l47 = 0; (l47 < 3); l47 = (l47 + 1)) {
fRec7[l47] = 0.0f;
}
for (int l48 = 0; (l48 < 3); l48 = (l48 + 1)) {
fRec8[l48] = 0.0f;
}
for (int l49 = 0; (l49 < 3); l49 = (l49 + 1)) {
fRec9[l49] = 0.0f;
}
for (int l50 = 0; (l50 < 3); l50 = (l50 + 1)) {
fRec10[l50] = 0.0f;
}
for (int l51 = 0; (l51 < 3); l51 = (l51 + 1)) {
fRec11[l51] = 0.0f;
}
for (int l52 = 0; (l52 < 3); l52 = (l52 + 1)) {
fRec3[l52] = 0.0f;
}
for (int l53 = 0; (l53 < 3); l53 = (l53 + 1)) {
fRec2[l53] = 0.0f;
}
for (int l54 = 0; (l54 < 3); l54 = (l54 + 1)) {
fRec45[l54] = 0.0f;
}
for (int l55 = 0; (l55 < 3); l55 = (l55 + 1)) {
fRec44[l55] = 0.0f;
}
}
virtual void init(int samplingFreq) {
classInit(samplingFreq);
instanceInit(samplingFreq);
}
virtual void instanceInit(int samplingFreq) {
instanceConstants(samplingFreq);
instanceResetUserInterface();
instanceClear();
}
virtual mydsp* clone() {
return new mydsp();
}
virtual int getSampleRate() {
return fSamplingFreq;
}
virtual void buildUserInterface(UI* ui_interface) {
ui_interface->openHorizontalBox("zita");
ui_interface->addVerticalSlider("dryWet", &fVslider1, 0.0f, -1.0f, 1.0f, 0.00999999978f);
ui_interface->addVerticalSlider("level", &fVslider0, -6.0f, -70.0f, 40.0f, 0.100000001f);
ui_interface->closeBox();
}
virtual void compute(int count, FAUSTFLOAT** inputs, FAUSTFLOAT** outputs) {
FAUSTFLOAT* input0 = inputs[0];
FAUSTFLOAT* input1 = inputs[1];
FAUSTFLOAT* output0 = outputs[0];
FAUSTFLOAT* output1 = outputs[1];
float fSlow0 = (0.00100000005f * std::pow(10.0f, (0.0500000007f * float(fVslider0))));
float fSlow1 = (0.00100000005f * float(fVslider1));
for (int i = 0; (i < count); i = (i + 1)) {
fRec0[0] = (fSlow0 + (0.999000013f * fRec0[1]));
fRec1[0] = (fSlow1 + (0.999000013f * fRec1[1]));
float fTemp0 = (fRec1[0] + 1.0f);
float fTemp1 = float(input0[i]);
fVec0[(IOTA & 4095)] = fTemp1;
float fTemp2 = (1.0f - (0.5f * fTemp0));
fRec15[0] = (0.0f - (fConst18 * ((fConst19 * fRec15[1]) - (fRec11[1] + fRec11[2]))));
fRec14[0] = ((fConst14 * fRec14[1]) + (fConst15 * (fRec11[1] + (fConst16 * fRec15[0]))));
fVec1[(IOTA & 2047)] = ((0.353553385f * fRec14[0]) + 9.99999968e-21f);
float fTemp3 = float(input1[i]);
fVec2[(IOTA & 4095)] = fTemp3;
float fTemp4 = (0.300000012f * fVec2[((IOTA - iConst22) & 4095)]);
float fTemp5 = (((0.600000024f * fRec12[1]) + fVec1[((IOTA - iConst21) & 2047)]) - fTemp4);
fVec3[(IOTA & 2047)] = fTemp5;
fRec12[0] = fVec3[((IOTA - iConst23) & 2047)];
float fRec13 = (0.0f - (0.600000024f * fTemp5));
fRec19[0] = (0.0f - (fConst18 * ((fConst19 * fRec19[1]) - (fRec7[1] + fRec7[2]))));
fRec18[0] = ((fConst32 * fRec18[1]) + (fConst33 * (fRec7[1] + (fConst34 * fRec19[0]))));
fVec4[(IOTA & 2047)] = ((0.353553385f * fRec18[0]) + 9.99999968e-21f);
float fTemp6 = (((0.600000024f * fRec16[1]) + fVec4[((IOTA - iConst36) & 2047)]) - fTemp4);
fVec5[(IOTA & 4095)] = fTemp6;
fRec16[0] = fVec5[((IOTA - iConst37) & 4095)];
float fRec17 = (0.0f - (0.600000024f * fTemp6));
fRec23[0] = (0.0f - (fConst18 * ((fConst19 * fRec23[1]) - (fRec9[1] + fRec9[2]))));
fRec22[0] = ((fConst46 * fRec22[1]) + (fConst47 * (fRec9[1] + (fConst48 * fRec23[0]))));
fVec6[(IOTA & 2047)] = ((0.353553385f * fRec22[0]) + 9.99999968e-21f);
float fTemp7 = (fVec6[((IOTA - iConst50) & 2047)] + (fTemp4 + (0.600000024f * fRec20[1])));
fVec7[(IOTA & 4095)] = fTemp7;
fRec20[0] = fVec7[((IOTA - iConst51) & 4095)];
float fRec21 = (0.0f - (0.600000024f * fTemp7));
fRec27[0] = (0.0f - (fConst18 * ((fConst19 * fRec27[1]) - (fRec5[1] + fRec5[2]))));
fRec26[0] = ((fConst60 * fRec26[1]) + (fConst61 * (fRec5[1] + (fConst62 * fRec27[0]))));
fVec8[(IOTA & 2047)] = ((0.353553385f * fRec26[0]) + 9.99999968e-21f);
float fTemp8 = (fTemp4 + ((0.600000024f * fRec24[1]) + fVec8[((IOTA - iConst64) & 2047)]));
fVec9[(IOTA & 4095)] = fTemp8;
fRec24[0] = fVec9[((IOTA - iConst65) & 4095)];
float fRec25 = (0.0f - (0.600000024f * fTemp8));
fRec31[0] = (0.0f - (fConst18 * ((fConst19 * fRec31[1]) - (fRec10[1] + fRec10[2]))));
fRec30[0] = ((fConst74 * fRec30[1]) + (fConst75 * (fRec10[1] + (fConst76 * fRec31[0]))));
fVec10[(IOTA & 2047)] = ((0.353553385f * fRec30[0]) + 9.99999968e-21f);
float fTemp9 = (0.300000012f * fVec0[((IOTA - iConst22) & 4095)]);
float fTemp10 = (fVec10[((IOTA - iConst78) & 2047)] - (fTemp9 + (0.600000024f * fRec28[1])));
fVec11[(IOTA & 2047)] = fTemp10;
fRec28[0] = fVec11[((IOTA - iConst79) & 2047)];
float fRec29 = (0.600000024f * fTemp10);
fRec35[0] = (0.0f - (fConst18 * ((fConst19 * fRec35[1]) - (fRec6[1] + fRec6[2]))));
fRec34[0] = ((fConst88 * fRec34[1]) + (fConst89 * (fRec6[1] + (fConst90 * fRec35[0]))));
fVec12[(IOTA & 2047)] = ((0.353553385f * fRec34[0]) + 9.99999968e-21f);
float fTemp11 = (fVec12[((IOTA - iConst92) & 2047)] - (fTemp9 + (0.600000024f * fRec32[1])));
fVec13[(IOTA & 4095)] = fTemp11;
fRec32[0] = fVec13[((IOTA - iConst93) & 4095)];
float fRec33 = (0.600000024f * fTemp11);
fRec39[0] = (0.0f - (fConst18 * ((fConst19 * fRec39[1]) - (fRec8[1] + fRec8[2]))));
fRec38[0] = ((fConst102 * fRec38[1]) + (fConst103 * (fRec8[1] + (fConst104 * fRec39[0]))));
fVec14[(IOTA & 2047)] = ((0.353553385f * fRec38[0]) + 9.99999968e-21f);
float fTemp12 = ((fTemp9 + fVec14[((IOTA - iConst106) & 2047)]) - (0.600000024f * fRec36[1]));
fVec15[(IOTA & 4095)] = fTemp12;
fRec36[0] = fVec15[((IOTA - iConst107) & 4095)];
float fRec37 = (0.600000024f * fTemp12);
fRec43[0] = (0.0f - (fConst18 * ((fConst19 * fRec43[1]) - (fRec4[1] + fRec4[2]))));
fRec42[0] = ((fConst116 * fRec42[1]) + (fConst117 * (fRec4[1] + (fConst118 * fRec43[0]))));
fVec16[(IOTA & 2047)] = ((0.353553385f * fRec42[0]) + 9.99999968e-21f);
float fTemp13 = ((fVec16[((IOTA - iConst120) & 2047)] + fTemp9) - (0.600000024f * fRec40[1]));
fVec17[(IOTA & 2047)] = fTemp13;
fRec40[0] = fVec17[((IOTA - iConst121) & 2047)];
float fRec41 = (0.600000024f * fTemp13);
float fTemp14 = (fRec41 + fRec37);
float fTemp15 = (fRec29 + (fRec33 + fTemp14));
fRec4[0] = (fRec12[1] + (fRec16[1] + (fRec20[1] + (fRec24[1] + (fRec28[1] + (fRec32[1] + (fRec36[1] + (fRec40[1] + (fRec13 + (fRec17 + (fRec21 + (fRec25 + fTemp15))))))))))));
fRec5[0] = ((fRec28[1] + (fRec32[1] + (fRec36[1] + (fRec40[1] + fTemp15)))) - (fRec12[1] + (fRec16[1] + (fRec20[1] + (fRec24[1] + (fRec13 + (fRec17 + (fRec25 + fRec21))))))));
float fTemp16 = (fRec33 + fRec29);
fRec6[0] = ((fRec20[1] + (fRec24[1] + (fRec36[1] + (fRec40[1] + (fRec21 + (fRec25 + fTemp14)))))) - (fRec12[1] + (fRec16[1] + (fRec28[1] + (fRec32[1] + (fRec13 + (fRec17 + fTemp16)))))));
fRec7[0] = ((fRec12[1] + (fRec16[1] + (fRec36[1] + (fRec40[1] + (fRec13 + (fRec17 + fTemp14)))))) - (fRec20[1] + (fRec24[1] + (fRec28[1] + (fRec32[1] + (fRec21 + (fRec25 + fTemp16)))))));
float fTemp17 = (fRec41 + fRec33);
float fTemp18 = (fRec37 + fRec29);
fRec8[0] = ((fRec16[1] + (fRec24[1] + (fRec32[1] + (fRec40[1] + (fRec17 + (fRec25 + fTemp17)))))) - (fRec12[1] + (fRec20[1] + (fRec28[1] + (fRec36[1] + (fRec13 + (fRec21 + fTemp18)))))));
fRec9[0] = ((fRec12[1] + (fRec20[1] + (fRec32[1] + (fRec40[1] + (fRec13 + (fRec21 + fTemp17)))))) - (fRec16[1] + (fRec24[1] + (fRec28[1] + (fRec36[1] + (fRec17 + (fRec25 + fTemp18)))))));
float fTemp19 = (fRec41 + fRec29);
float fTemp20 = (fRec37 + fRec33);
fRec10[0] = ((fRec12[1] + (fRec24[1] + (fRec28[1] + (fRec40[1] + (fRec13 + (fRec25 + fTemp19)))))) - (fRec16[1] + (fRec20[1] + (fRec32[1] + (fRec36[1] + (fRec17 + (fRec21 + fTemp20)))))));
fRec11[0] = ((fRec16[1] + (fRec20[1] + (fRec28[1] + (fRec40[1] + (fRec17 + (fRec21 + fTemp19)))))) - (fRec12[1] + (fRec24[1] + (fRec32[1] + (fRec36[1] + (fRec13 + (fRec25 + fTemp20)))))));
float fTemp21 = (0.370000005f * (fRec5[0] + fRec6[0]));
float fTemp22 = (fConst122 * fRec3[1]);
fRec3[0] = (fTemp21 - (fTemp22 + (fConst4 * fRec3[2])));
float fTemp23 = (fConst4 * fRec3[0]);
float fTemp24 = (0.5f * ((fTemp23 + (fRec3[2] + (fTemp21 + fTemp22))) + ((fTemp23 + (fTemp22 + fRec3[2])) - fTemp21)));
float fTemp25 = (fConst123 * fRec2[1]);
fRec2[0] = (fTemp24 - (fTemp25 + (fConst2 * fRec2[2])));
float fTemp26 = (fConst2 * fRec2[0]);
output0[i] = FAUSTFLOAT((0.5f * (fRec0[0] * ((fTemp0 * fTemp1) + (fTemp2 * ((fTemp26 + (fRec2[2] + (fTemp24 + fTemp25))) + ((fTemp26 + (fTemp25 + fRec2[2])) - fTemp24)))))));
float fTemp27 = (0.370000005f * (fRec5[0] - fRec6[0]));
float fTemp28 = (fConst122 * fRec45[1]);
fRec45[0] = (fTemp27 - (fTemp28 + (fConst4 * fRec45[2])));
float fTemp29 = (fConst4 * fRec45[0]);
float fTemp30 = (0.5f * ((fTemp29 + (fRec45[2] + (fTemp27 + fTemp28))) + ((fTemp29 + (fTemp28 + fRec45[2])) - fTemp27)));
float fTemp31 = (fConst123 * fRec44[1]);
fRec44[0] = (fTemp30 - (fTemp31 + (fConst2 * fRec44[2])));
float fTemp32 = (fConst2 * fRec44[0]);
output1[i] = FAUSTFLOAT((0.5f * (fRec0[0] * ((fTemp0 * fTemp3) + (fTemp2 * ((fTemp32 + (fRec44[2] + (fTemp30 + fTemp31))) + ((fTemp32 + (fTemp31 + fRec44[2])) - fTemp30)))))));
fRec0[1] = fRec0[0];
fRec1[1] = fRec1[0];
IOTA = (IOTA + 1);
fRec15[1] = fRec15[0];
fRec14[1] = fRec14[0];
fRec12[1] = fRec12[0];
fRec19[1] = fRec19[0];
fRec18[1] = fRec18[0];
fRec16[1] = fRec16[0];
fRec23[1] = fRec23[0];
fRec22[1] = fRec22[0];
fRec20[1] = fRec20[0];
fRec27[1] = fRec27[0];
fRec26[1] = fRec26[0];
fRec24[1] = fRec24[0];
fRec31[1] = fRec31[0];
fRec30[1] = fRec30[0];
fRec28[1] = fRec28[0];
fRec35[1] = fRec35[0];
fRec34[1] = fRec34[0];
fRec32[1] = fRec32[0];
fRec39[1] = fRec39[0];
fRec38[1] = fRec38[0];
fRec36[1] = fRec36[0];
fRec43[1] = fRec43[0];
fRec42[1] = fRec42[0];
fRec40[1] = fRec40[0];
fRec4[2] = fRec4[1];
fRec4[1] = fRec4[0];
fRec5[2] = fRec5[1];
fRec5[1] = fRec5[0];
fRec6[2] = fRec6[1];
fRec6[1] = fRec6[0];
fRec7[2] = fRec7[1];
fRec7[1] = fRec7[0];
fRec8[2] = fRec8[1];
fRec8[1] = fRec8[0];
fRec9[2] = fRec9[1];
fRec9[1] = fRec9[0];
fRec10[2] = fRec10[1];
fRec10[1] = fRec10[0];
fRec11[2] = fRec11[1];
fRec11[1] = fRec11[0];
fRec3[2] = fRec3[1];
fRec3[1] = fRec3[0];
fRec2[2] = fRec2[1];
fRec2[1] = fRec2[0];
fRec45[2] = fRec45[1];
fRec45[1] = fRec45[0];
fRec44[2] = fRec44[1];
fRec44[1] = fRec44[0];
}
}
};
#define MULT_16 2147483647
#define DIV_16 4.6566129e-10
FaustZita::FaustZita() : AudioStream((fDSP = new mydsp())->getNumInputs(), new audio_block_t*[fDSP->getNumInputs()])
{
UI = new MapUI();
fDSP->init(AUDIO_SAMPLE_RATE_EXACT);
fDSP->buildUserInterface(UI);
// allocating Faust inputs
if(fDSP->getNumInputs()>0){
fInChannel = new float*[fDSP->getNumInputs()];
for(int i=0; i<fDSP->getNumInputs(); i++){
fInChannel[i] = new float[AUDIO_BLOCK_SAMPLES];
}
}
// allocating Faust outputs
if(fDSP->getNumOutputs()>0){
fOutChannel = new float*[fDSP->getNumOutputs()];
for(int i=0; i<fDSP->getNumOutputs(); i++){
fOutChannel[i] = new float[AUDIO_BLOCK_SAMPLES];
}
}
}
void FaustZita::update(void)
{
audio_block_t *inBlock[fDSP->getNumInputs()], *outBlock[fDSP->getNumOutputs()];
int32_t val;
if(fDSP->getNumInputs()>0){
for(int channel=0; channel<fDSP->getNumInputs(); channel++){
inBlock[channel] = receiveReadOnly(channel);
for(int i=0; i < AUDIO_BLOCK_SAMPLES; i++) {
val = inBlock[channel]->data[i] << 16;
fInChannel[channel][i] = val*DIV_16;
}
release(inBlock[channel]);
}
}
fDSP->compute(AUDIO_BLOCK_SAMPLES,fInChannel,fOutChannel);
for(int channel=0; channel<fDSP->getNumOutputs(); channel++){
outBlock[channel] = allocate();
if(outBlock[channel]){
for(int i=0; i < AUDIO_BLOCK_SAMPLES; i++) {
val = fOutChannel[channel][i]*MULT_16;
outBlock[channel]->data[i] = val >> 16;
}
transmit(outBlock[channel],channel);
release(outBlock[channel]);
}
}
}
FaustZita::~FaustZita(){
delete fDSP;
delete UI;
for(int i=0; i<fDSP->getNumInputs(); i++) {
delete[] fInChannel[i];
}
delete [] fInChannel;
for (int i=0; i < fDSP->getNumOutputs(); i++) {
delete[] fOutChannel[i];
}
delete [] fOutChannel;
}
void FaustZita::setParamValue(const std::string& path, float value){
UI->setParamValue(path,value);
}
#endif
What do you think is the problem?
Best regards, Seoul Open Media
This could indeed be a memory size related problem:
- I see you're still using Faust 2.15.5, first you should use a more recent version like the latest published one
- secondly your could try to use the
-dltoption to change the way delays lines are generated. Try for instance-dlt 0. Read here for more info.
Hi @sletz, I am having trouble with this as well. I am using Teensy 4.1. I am working directly from the demo.
> cat FaustZita.dsp
import("stdfaust.lib");
process = dm.zita_light;
I generated the code with
faust2teensy -lib FaustZita.dsp -dlt 0
I compiled faust recently
> faust -v
FAUST Version 2.40.7
Embedded backends:
DSP to C
DSP to C++
DSP to CSharp
DSP to DLang
DSP to FIR
DSP to Interpreter
DSP to Java
DSP to Julia
DSP to LLVM IR
DSP to old C++
DSP to Rust
DSP to SOUL
DSP to WebAssembly (wast/wasm)
Build with LLVM version 10.0.0
Copyright (C) 2002-2022, GRAME - Centre National de Creation Musicale. All rights reserved.
My faustZita.ino is
#include <Audio.h>
#include "FaustZita.h"
FaustZita faustZita;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioConnection patchCord0(in,0,faustZita,0);
AudioConnection patchCord1(in,1,faustZita,1);
AudioConnection patchCord2(faustZita,0,out,0);
AudioConnection patchCord3(faustZita,1,out,1);
void setup() {
AudioMemory(6);
audioShield.enable();
faustZita.setParamValue("level",0.5);
faustZita.setParamValue("dryWet",0.5);
}
void loop() {
delay(10);
}
I can confirm that faust code from the other demo (the Band-Limited Sawtooth Oscillator) works great. I can also confirm that my audio input is working because this sketch works fine when bypassing faust.
#include <Audio.h>
#include "FaustZita.h"
FaustZita faustZita;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioConnection patchCord0(in,0,out,0);
AudioConnection patchCord1(in,1,out,1);
void setup() {
AudioMemory(6);
audioShield.enable();
audioShield.inputSelect(AUDIO_INPUT_LINEIN);
audioShield.volume(0.5);
// faustZita.setParamValue("level",0.5);
// faustZita.setParamValue("dryWet",0.5);
}
void loop() {
delay(10);
}
However, there might be a bug abound. Because if I uncomment those faustZita lines even though I'm not patching through faustZita, I will stop getting audio input directed to the output! E.g. the following code no longer works:
#include <Audio.h>
#include "FaustZita.h"
FaustZita faustZita;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioConnection patchCord0(in,0,out,0);
AudioConnection patchCord1(in,1,out,1);
void setup() {
AudioMemory(6);
audioShield.enable();
audioShield.inputSelect(AUDIO_INPUT_LINEIN);
audioShield.volume(0.5);
faustZita.setParamValue("level",0.5);
faustZita.setParamValue("dryWet",0.5);
}
void loop() {
delay(10);
}
There's no output in the serial command so I can't tell if its failing somehow. Just no audio and serial output whenever the setParamValue is enabled.
Actually I just saw issue https://github.com/grame-cncm/faust/issues/650 and I am noticing the same thing! When the setParamValue is defined there is a restart in the serial terminal like every 8-10 seconds. Could it be crashing like in that issue?
I would sincerely love some help on this and I promise not to ghost you like others in these issues...
Could it be simply a memory issue? Does the simpler freeverb works ? (see https://faustdoc.grame.fr/examples/reverb/#freeverb)
I'm not sure - I assumed the memory would be fine because it is listed in the tutorial. I tried using Freeverb and still got no input. Actually this one is weird too because I had to comment out the Freeverb to get my input patched to output.
This makes sound:
#include <Audio.h>
//#include "Freeverb.h"
//
//Freeverb fx;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioConnection patchCord0(in,0,out,0);
AudioConnection patchCord1(in,1,out,1);
void setup() {
AudioMemory(6);
audioShield.enable();
audioShield.inputSelect(AUDIO_INPUT_LINEIN);
audioShield.volume(0.5);
}
void loop() {
delay(100);
}
But this doesn't make sound even though the Freeverb is not hooked up:
#include <Audio.h>
#include "Freeverb.h"
Freeverb fx;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioConnection patchCord0(in,0,out,0);
AudioConnection patchCord1(in,1,out,1);
void setup() {
AudioMemory(6);
audioShield.enable();
audioShield.inputSelect(AUDIO_INPUT_LINEIN);
audioShield.volume(0.5);
}
void loop() {
delay(100);
}
Again in that second script te serial terminal is restarting constantly.
I also just tried using the source code published on the website (https://faustdoc.grame.fr/tutorials/teensy/misc/teensyEffect.zip) and it is not working either - same issue: no input sound coming out.
It does!
#include <Audio.h>
#include "Simple.h"
Simple fx;
AudioInputI2S in;
AudioOutputI2S out;
AudioControlSGTL5000 audioShield;
AudioConnection patchCord0(in,0,fx,0);
AudioConnection patchCord1(in,1,fx,1);
AudioConnection patchCord2(fx,0,out,0);
AudioConnection patchCord3(fx,1,out,1);
void setup() {
AudioMemory(20);
audioShield.enable();
audioShield.inputSelect(AUDIO_INPUT_LINEIN);
audioShield.volume(0.1);
}
void loop() {
delay(100);
}
works
So then the possible problem is the memory size of Freeverb or FaustZita object globally allocated, that is as a static variable. I guess you should look in this direction, possibly starting from simple FX then more complex ones.
I was able to get FaustZita working on Teensy by modifying the .cpp slightly. See MarkT's post here about null incoming blocks. https://forum.pjrc.com/threads/64311-Memory-issue-with-AudioLibrary-Teensyduino-since-1-52-1-54b?p=258857&viewfull=1#post258857
Can you prepare a PR to correct the architecture file https://github.com/grame-cncm/faust/blob/master-dev/architecture/teensy/teensy.cpp ? Or at least send me the modified file ? Thanks.
Correction: the pointer is already checked (https://github.com/grame-cncm/faust/blob/a7e111e0a77effb2ac7a778830cf0d365a6afafe/architecture/teensy/teensy.cpp#L173)...
So @MattKuebrich do you really have a working version with the current master-dev state?
I was using the files from the tutorial here. https://faustdoc.grame.fr/tutorials/teensy/misc/teensyEffect.zip
Here's that .cpp with the MarkT's fix added. FaustZita.zip
Hopefully that helps.
I can confirm that the tutorial code works with MarkT's fix, thanks @MattKuebrich.
But I can't transfer the fix from that code to any of my code generated with faust recently. It looks like the faust templates have changed a bit for teensy. From that fixed tutorial FaustZita.cpp the block around the snippet looks like
void FaustZita::update(void)
{
audio_block_t *inBlock[fDSP->getNumInputs()], *outBlock[fDSP->getNumOutputs()];
int32_t val;
if(fDSP->getNumInputs()>0){
for(int channel=0; channel<fDSP->getNumInputs(); channel++){
inBlock[channel] = receiveReadOnly(channel);
if (inBlock[channel] == NULL)
continue ;
for(int i=0; i < AUDIO_BLOCK_SAMPLES; i++) {
val = inBlock[channel]->data[i] << 16;
fInChannel[channel][i] = val*DIV_16;
}
release(inBlock[channel]);
}
}
But in the newly generated snippets it looks like:
void Volume::updateImp(void)
{
#if MIDICTRL
// Process the MIDI messages received by the Teensy
fMIDIHandler->processMidi();
// Synchronize all GUI controllers
GUI::updateAllGuis();
#endif
if (INPUTS > 0) {
audio_block_t* inBlock[INPUTS];
for (int channel = 0; channel < INPUTS; channel++) {
inBlock[channel] = receiveReadOnly(channel);
// THESE NEXT TWO LINES ADDED IN BY ME
if (inBlock[channel] == NULL)
continue ;
if (inBlock[channel]) {
for (int i = 0; i < AUDIO_BLOCK_SAMPLES; i++) {
int32_t val = inBlock[channel]->data[i] << 16;
fInChannel[channel][i] = val*DIV_16;
}
release(inBlock[channel]);
} else {
memset(fInChannel[channel], 0, AUDIO_BLOCK_SAMPLES * sizeof(float));
}
}
}
(generated for a simple gain control). The applied change doesn't do anything here, as @sletz says it looks like the pointer is already checked....
Do you know what version of faust was used to generate the tutorial? Maybe I can regress faust2teensy and try that.
In the generated file you get this kind of header
/* ------------------------------------------------------------
name: "FaustZita"
Code generated with Faust 2.15.5 (https://faust.grame.fr)
Compilation options: cpp, -scal -ftz 0
------------------------------------------------------------ */
But the annoying thing is that latest Faust + architecture files does not work anymore. My guess is that the generated DSP takes more memory. You could debug that using sizeof , printing sizeof(mydsp), and sizeof(MapUI) here https://github.com/grame-cncm/faust/blob/master-dev/architecture/teensy/teensy.cpp#L99
So comparing what the olde 2.15.5 generated version gives and what the current version gives.
Okay I got some progress. @sletz, you are of course right, it is a memory issue. I generated using 2.15.11 (could not find a tag for 2.15.5) and got a very similar file except all the numbers were higher still. This version also did not work. Even if I made the if (inBlock[channel] == NULL) fix. No audio input was processed.
But 2.15.11 wasn't exactly the same as 2.15.5, here are some example differences:
< = newly generated with 2.15.11
< iConst21 = int(std::min<float>(16384.0f, std::max<float>(0.0f, (fConst6 - fConst20))));
< iConst22 = int(std::min<float>(8192.0f, std::max<float>(0.0f, (0.0599999987f * fConst0))));
---
> working version from 2.15.5
> iConst21 = int(std::min<float>(1024.0f, std::max<float>(0.0f, (fConst6 - fConst20))));
> iConst22 = int(std::min<float>(2700.0f, std::max<float>(0.0f, (0.0599999987f * fConst0))));
So, it did work once I changed all the numbers to lower numbers, as in the working file (16384 -> 1024, 8192 -> 2700, etc. etc.).
I'm so very new to teensy and I'm a little naive at why Teensy4.1 needs this change. I did install the PSRAM so supposedly I should be able to allocate more memory for this? I guess I can just prefix all the buffers with EXTMEM? I'm not sure. But at least its working now. I'll go back to the latest version of Faust and try dropping the numbers down to see if this fix replicates.
I did not, I only did what was in that working example: -scal -ftz 0. I can try that when I get back to the device.
No luck on using the latest dev and reducing all the memory to low amounts (2048) (FAUST Version 2.40.9)
No luck using the -dlt 0 option on version 2.15.11.
But using 2.15.11, I did figure out that I can move some big arrays out of the dsp class into the global variables and set them to EXTMEM and its working (instead of reducing all the amounts of memory):
EXTMEM float fVec0[16384];
EXTMEM float fVec1[32768];
EXTMEM float fVec4[32768];
EXTMEM float fVec8[32768];
EXTMEM float fVec2[16384];
class mydsp : public dsp {
private:
...
No, its not working the same on 2.40.9. I tried exactly the same as above (fVec0,1,2,4, and 8 as EXTMEM) Also reducing the memory for those variables didn't make a difference either. I tried with and without -dlt 0 and neither worked with 2.40.9
It might take me a bit, but I could do a git bisection to find where it changes between 2.15.11 and 2.40.9. Unless there are some other things I could try first
Trying the technique showed here: https://github.com/modularev/faust/commit/f84a55c0393b3e3cfa87a73de3bc7c28e57225b2
Or https://join.slack.com/t/faustaudio/shared_invite/zt-a624szlz-fL4v2DTR~ZGlI7wARryT7g
Hi - Was there ever a resolution to this? I also can't get audio through the teensy using faust
