timowest / flauta Goto Github PK

hi,
I made a python script who compare the outputs. I don't include receptivity. Maybe the average is not a very good idea for obtain the error, any sugestion is welcome.

please run the script like this:
$ python comp.py

replace the constants with controllable parameters

make sure that resonator implementations are in sync

ideas

• notes: 70/75/80
• three pressure levels
  ?!?

test for Resonator.cpp

Add printing of parameters in faust tests

clean up resonator definition in faust

excitation, resonator and flauta

defined in jet (excitatiion.dsp)

Also fix k parameter in junction

Bottom up verification of the code

• Excitation
  • Blow
    • Envelope
    • Vibrato
    • Bernoulli
  • Jet
    • Receptivity
  • Sources
    • Jet drive
    • Turbulence
    • Vortex
• Resonator

The basic procedure for testing the code is drive the block in C++ and Faust with a common float typed WAV file and verify that the output WAV file matches.

We start with Bernoulli

Faust test is already available. Here is the contract

in : mouth pressure
out : uj, uj_steady, impulse

Developing a new C++ Testing program that allows us to Test the inputs and outputs of each Module (in C++).

 - Multiple options for input testing : Sinewave, Triangular wave, Impulse, Noise, Steps.
 - Multiple Inputs : 1, 2 or 3 inputs max.
 - Doesn't matter the length of the input. (using Vector from Stl)
 - All in form of a function in() and out() to avoid modifying the original C++ Code.
 - All the input examples are generated in MatLab and saved in CSV Files.

all these modifications where proposed so we can eliminate variables that can cause bad testing outputs.

Make sure that the filter declarations match those in the C++ code.

IIR filters are declared as tf2, tf3 and tf4 function invocations.

We need to verified which variables are inputs into the Receptivity Module, because i notice that the Faust version has 3 inputs and the C++ version has only 2 inputs (Acoustic Velocity and Hydrodynamic Feedback). the 3rd one is generated by 'Flow In' Module with this 2 inputs from Receptivity. (or at least i think it is, please someone take a look at this! ^^)

Add vibrato envelope to blow in excitation.dsp

The contract on the Faust side is

In: Vac, Pp
out : sources, impulse

it will be more like the C++ version and now the gate button doesn't works good

test version should contain constants instead of GUI declarations