Using which architecture?

Thanks.

Stéphane 

> Le 7 juil. 2021 à 23:34, Julius Smith <julius.sm...@gmail.com> a écrit :
> 
> Hmmm, this test points to a hbargraph gating behavior in its display only, 
> although I did not run it for a long time:
> 
> import("stdfaust.lib");
> Tg = 0.4;
> zi = an.ms_envelope_rect(Tg);
> //gain = hslider("Gain [unit:dB]",-10,-70,0,0.1) : ba.db2linear;
> gain = 0.01; // -40 dB - ok in Octave (I get -44.744 dB)
> sig = no.noise * gain;
> level = zi : ba.linear2db : *(0.5);
> //process = level, attach(sig, (sig : level : hbargraph("test",-70,0)));
> process = sig : level <: _, hbargraph("test",-70,0);
> 
> On Wed, Jul 7, 2021 at 12:50 PM Juan Carlos Blancas <lav...@gmail.com> wrote:
> Hi Klaus, 
> 
> Same here, it seems there is something with ms and rms_envelope.
> 
> Best,
> Juan Carlos
> 
> import("stdfaust.lib");
> 
> rms(n) = _ : square : mean(n) : sqrt;
> square(x) = x * x;
> mean(n) = float2fix : integrate(n) : fix2float : /(n); 
> integrate(n,x) = x - x@n : +~_;
> float2fix(x) = int(x*(1<<20));      
> fix2float(x) = float(x)/(1<<20);    
> 
> Tg = 0.4;
> zi  = an.ms_envelope_rect(Tg);
> ziR = an.rms_envelope_rect(Tg);
> 
> process = no.noise*1.737 * 
> ba.db2linear(hslider("[0]g[unit:dB]",-20,-95,-10,0.1)) <:
> attach(_, rms(ma.SR*Tg) : ba.linear2db : hbargraph("[1]rms",-95,0)),
> attach(_, sqrt(zi) : ba.linear2db : 
> hbargraph("[2]sqrt(ms_envelope_rect)",-95,0)),
> attach(_, ziR : ba.linear2db : hbargraph("[3]rms_envelope_rect",-95,0));
> 
> 
> > El 7 jul 2021, a las 9:59, Klaus Scheuermann <kla...@posteo.de> escribió:
> > 
> > Hi all,
> > I did some testing and
> > 
> > an.ms_envelope_rect()
> > 
> > seems to show some strange behaviour (at least to me). Here is a video
> > of the test:
> > https://cloud.4ohm.de/s/64caEPBqxXeRMt5
> > 
> > The audio is white noise and the testing code is:
> > 
> > import("stdfaust.lib");
> > Tg = 0.4;
> > zi = an.ms_envelope_rect(Tg);
> > process = _ : zi : ba.linear2db : hbargraph("test",-95,0);
> > 
> > Could you please verify?
> > 
> > Thanks, Klaus
> > 
> > 
> > 
> > On 05.07.21 20:16, Julius Smith wrote:
> >> Hmmm, '!' means "block the signal", but attach should save the bargraph
> >> from being optimized away as a result.  Maybe I misremembered the
> >> argument order to attach?  While it's very simple in concept, it can be
> >> confusing in practice.
> >> 
> >> I chose not to have a gate at all, but you can grab one from
> >> misceffects.lib if you like.  Low volume should not give -infinity,
> >> that's a bug, but zero should, and zero should become MIN as I mentioned
> >> so -infinity should never happen.
> >> 
> >> Cheers,
> >> Julius
> >> 
> >> 
> >> On Mon, Jul 5, 2021 at 10:39 AM Klaus Scheuermann <kla...@posteo.de
> >> <mailto:kla...@posteo.de>> wrote:
> >> 
> >>    Cheers Julius,
> >> 
> >> 
> >> 
> >>    At least I understood the 'attach' primitive now ;) Thanks.
> >> 
> >> 
> >> 
> >>    This does not show any meter here...
> >>    process(x,y) = x,y <: (_,_), attach(x, (Lk2 : vbargraph("LUFS",-90,0)))
> >>    : _,_,!;
> >> 
> >>    But this does for some reason (although the output is 3-channel then):
> >>    process(x,y) = x,y <: (_,_), attach(x, (Lk2 : vbargraph("LUFS",-90,0)))
> >>    : _,_,_;
> >> 
> >>    What does the '!' do?
> >> 
> >> 
> >> 
> >>    I still don't quite get the gating topic. In my understanding, the meter
> >>    should hold the current value if the input signal drops below a
> >>    threshold. In your version, the meter drops to -infinity when very low
> >>    volume content is played.
> >> 
> >>    Which part of your code does the gating?
> >> 
> >>    Many thanks,
> >>    Klaus
> >> 
> >> 
> >> 
> >>    On 05.07.21 18:06, Julius Smith wrote:
> >>> Hi Klaus,
> >>> 
> >>> Yes, I agree the filters are close enough.  I bet that the shelf is
> >>> exactly correct if we determined the exact transition frequency, and
> >>> that the Butterworth highpass is close enough to the
> >>    Bessel-or-whatever
> >>> that is inexplicably not specified as a filter type, leaving it
> >>> sample-rate dependent.  I would bet large odds that the differences
> >>> cannot be reliably detected in listening tests.
> >>> 
> >>> Yes, I just looked again, and there are "gating blocks" defined,
> >>    each Tg
> >>> = 0.4 sec long, so that only ungated blocks are averaged to form a
> >>> longer term level-estimate.  What I wrote gives a "sliding gating
> >>> block", which can be lowpass filtered further, and/or gated, etc.  
> >>> Instead of a gate, I would simply replace 0 by ma.EPSILON so that the
> >>> log always works (good for avoiding denormals as well).
> >>> 
> >>> I believe stereo is supposed to be handled like this:
> >>> 
> >>> Lk2 = _,0,_,0,0 : Lk5;
> >>> process(x,y) = Lk2(x,y);
> >>> 
> >>> or
> >>> 
> >>> Lk2 = Lk(0),Lk(2) :> 10 * log10 : -(0.691);
> >>> 
> >>> but since the center channel is processed identically to left
> >>    and right,
> >>> your solution also works.
> >>> 
> >>> Bypassing is normal Faust, e.g.,
> >>> 
> >>> process(x,y) = x,y <: (_,_), attach(x, (Lk2 :
> >>    vbargraph("LUFS",-90,0)))
> >>> : _,_,!;
> >>> 
> >>> Cheers,
> >>> Julius
> >>> 
> >>> 
> >>> On Mon, Jul 5, 2021 at 1:56 AM Klaus Scheuermann <kla...@posteo.de
> >>    <mailto:kla...@posteo.de>
> >>> <mailto:kla...@posteo.de <mailto:kla...@posteo.de>>> wrote:
> >>> 
> >>> 
> >>>      > I can never resist these things!   Faust makes it too
> >>    enjoyable :-)
> >>> 
> >>>      Glad you can't ;)
> >>> 
> >>>      I understood you approximate the filters with standard faust
> >>    filters.
> >>>      That is probably close enough for me :)
> >>> 
> >>>      I also get the part with the sliding window envelope. If I
> >>    wanted to
> >>>      make the meter follow slowlier, I would just widen the window
> >>    with Tg.
> >>> 
> >>>      The 'gating' part I don't understand for lack of mathematical
> >>    knowledge,
> >>>      but I suppose it is meant differently. When the input signal
> >>    falls below
> >>>      the gate threshold, the meter should stay at the current
> >>    value, not drop
> >>>      to -infinity, right? This is so 'silent' parts are not taken into
> >>>      account.
> >>> 
> >>>      If I wanted to make a stereo version it would be something like
> >>>      this, right?
> >>> 
> >>>      Lk2 = par(i,2, Lk(i)) :> 10 * log10 : -(0.691);
> >>>      process = _,_ : Lk2 : vbargraph("LUFS",-90,0);
> >>> 
> >>>      Probably very easy, but how do I attach this to a stereo
> >>    signal (passing
> >>>      through the stereo signal)?
> >>> 
> >>>      Thanks again!
> >>>      Klaus
> >>> 
> >>> 
> >>> 
> >>>      >
> >>>      > I made a pass, but there is a small scaling error.  I think
> >>    it can be
> >>>      > fixed by reducing boostFreqHz until the sine_test is nailed.
> >>>      > The highpass is close (and not a source of the scale error),
> >>    but I'm
> >>>      > using Butterworth instead of whatever they used.
> >>>      > I glossed over the discussion of "gating" in the spec, and
> >>    may have
> >>>      > missed something important there, but
> >>>      > I simply tried to make a sliding rectangular window, instead
> >>    of 75%
> >>>      > overlap, etc.
> >>>      >
> >>>      > If useful, let me know and I'll propose it for analyzers.lib!
> >>>      >
> >>>      > Cheers,
> >>>      > Julius
> >>>      >
> >>>      > import("stdfaust.lib");
> >>>      >
> >>>      > // Highpass:
> >>>      > // At 48 kHz, this is the right highpass filter (maybe a
> >>    Bessel or
> >>>      > Thiran filter?):
> >>>      > A48kHz = ( /* 1.0, */ -1.99004745483398, 0.99007225036621); 
> >>>      > B48kHz = (1.0, -2.0, 1.0); 
> >>>      > highpass48kHz = fi.iir(B48kHz,A48kHz);
> >>>      > highpass = fi.highpass(2, 40); // Butterworth highpass:
> >>    roll-off is a
> >>>      > little too sharp
> >>>      >
> >>>      > // High Shelf:
> >>>      > boostDB = 4;
> >>>      > boostFreqHz = 1430; // a little too high - they should give
> >>    us this!
> >>>      > highshelf = fi.high_shelf(boostDB, boostFreqHz); // Looks
> >>    very close,
> >>>      > but 1 kHz gain has to be nailed
> >>>      >
> >>>      > kfilter = highshelf : highpass;
> >>>      >
> >>>      > // Power sum:
> >>>      > Tg = 0.4; // spec calls for 75% overlap of successive
> >>    rectangular
> >>>      > windows - we're overlapping MUCH more (sliding window)
> >>>      > zi = an.ms_envelope_rect(Tg); // mean square: average power =
> >>>      energy/Tg
> >>>      > = integral of squared signal / Tg
> >>>      >
> >>>      > // Gain vector Gv = (GL,GR,GC,GLs,GRs):
> >>>      > N = 5;
> >>>      > Gv = (1, 1, 1, 1.41, 1.41); // left GL(-30deg), right GR
> >>    (30), center
> >>>      > GC(0), left surround GLs(-110), right surr. GRs(110)
> >>>      > G(i) = *(ba.take(i+1,Gv));
> >>>      > Lk(i) = kfilter : zi : G(i); // one channel, before summing
> >>    and before
> >>>      > taking dB and offsetting
> >>>      > LkDB(i) = Lk(i) : 10 * log10 : -(0.691); // Use this for a mono
> >>>      input signal
> >>>      >
> >>>      > // Five-channel surround input:
> >>>      > Lk5 = par(i,5,Lk(i)) :> 10 * log10 : -(0.691);
> >>>      >
> >>>      > // sine_test = os.oscrs(1000); // should give –3.01 LKFS, with
> >>>      > GL=GR=GC=1 (0dB) and GLs=GRs=1.41 (~1.5 dB)
> >>>      > sine_test = os.osc(1000);
> >>>      >            
> >>>      > process = sine_test : LkDB(0); // should read -3.01 LKFS -
> >>    high-shelf
> >>>      > gain at 1 kHz is critical
> >>>      > // process = 0,sine_test,0,0,0 : Lk5; // should read -3.01
> >>    LKFS for
> >>>      > left, center, and right
> >>>      > // Highpass test: process = 1-1' <: highpass, highpass48kHz;
> >>    // fft in
> >>>      > Octave
> >>>      > // High shelf test: process = 1-1' : highshelf; // fft in Octave
> >>>      >
> >>>      > On Sat, Jul 3, 2021 at 1:08 AM Klaus Scheuermann
> >>    <kla...@posteo.de <mailto:kla...@posteo.de>
> >>>      <mailto:kla...@posteo.de <mailto:kla...@posteo.de>>
> >>>      > <mailto:kla...@posteo.de <mailto:kla...@posteo.de>
> >>    <mailto:kla...@posteo.de <mailto:kla...@posteo.de>>>> wrote:
> >>>      >
> >>>      >     Hello everyone :)
> >>>      >
> >>>      >     Would someone be up for helping me implement an LUFS
> >>    loudness
> >>>      analyser
> >>>      >     in faust?
> >>>      >
> >>>      >     Or has someone done it already?
> >>>      >
> >>>      >     LUFS (aka LKFS) is becoming more and more the standard for
> >>>      loudness
> >>>      >     measurement in the audio industry. Youtube, Spotify and
> >>    broadcast
> >>>      >     stations use the concept to normalize loudness. A very
> >>>      positive side
> >>>      >     effect is, that loudness-wars are basically over.
> >>>      >
> >>>      >     I looked into it, but my programming skills clearly
> >>    don't match
> >>>      >     the level for implementing this.
> >>>      >
> >>>      >     Here is some resource about the topic:
> >>>      >
> >>>      >     https://en.wikipedia.org/wiki/LKFS
> >>    <https://en.wikipedia.org/wiki/LKFS>
> >>>      <https://en.wikipedia.org/wiki/LKFS
> >>    <https://en.wikipedia.org/wiki/LKFS>>
> >>>      <https://en.wikipedia.org/wiki/LKFS
> >>    <https://en.wikipedia.org/wiki/LKFS>
> >>>      <https://en.wikipedia.org/wiki/LKFS
> >>    <https://en.wikipedia.org/wiki/LKFS>>>
> >>>      >
> >>>      >     Specifications (in Annex 1):
> >>>      >   
> >>>    
> >>      
> >> https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf
> >>    
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf>
> >>>    
> >>     
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf
> >>    
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf>>
> >>>      >   
> >>>    
> >>      
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf
> >>    
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf>
> >>>    
> >>     
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf
> >>    
> >> <https://www.itu.int/dms_pubrec/itu-r/rec/bs/R-REC-BS.1770-3-201208-S!!PDF-E.pdf>>>
> >>>      >
> >>>      >     An implementation by 'klangfreund' in JUCE / C:
> >>>      >     https://github.com/klangfreund/LUFSMeter
> >>    <https://github.com/klangfreund/LUFSMeter>
> >>>      <https://github.com/klangfreund/LUFSMeter
> >>    <https://github.com/klangfreund/LUFSMeter>>
> >>>      >     <https://github.com/klangfreund/LUFSMeter
> >>    <https://github.com/klangfreund/LUFSMeter>
> >>>      <https://github.com/klangfreund/LUFSMeter
> >>    <https://github.com/klangfreund/LUFSMeter>>>
> >>>      >
> >>>      >     There is also a free LUFS Meter in JS / Reaper by
> >>    Geraint Luff.
> >>>      >     (The code can be seen in reaper, but I don't know if I
> >>    should
> >>>      paste it
> >>>      >     here.)
> >>>      >
> >>>      >     Please let me know if you are up for it!
> >>>      >
> >>>      >     Take care,
> >>>      >     Klaus
> >>>      >
> >>>      >
> >>>      >     _______________________________________________
> >>>      >     Faudiostream-users mailing list
> >>>      >     Faudiostream-users@lists.sourceforge.net
> >>    <mailto:Faudiostream-users@lists.sourceforge.net>
> >>>      <mailto:Faudiostream-users@lists.sourceforge.net
> >>    <mailto:Faudiostream-users@lists.sourceforge.net>>
> >>>      >     <mailto:Faudiostream-users@lists.sourceforge.net
> >>    <mailto:Faudiostream-users@lists.sourceforge.net>
> >>>      <mailto:Faudiostream-users@lists.sourceforge.net
> >>    <mailto:Faudiostream-users@lists.sourceforge.net>>>
> >>>      >   
> >>>    
> >>      https://lists.sourceforge.net/lists/listinfo/faudiostream-users
> >>    <https://lists.sourceforge.net/lists/listinfo/faudiostream-users>
> >>>    
> >>     <https://lists.sourceforge.net/lists/listinfo/faudiostream-users
> >>    <https://lists.sourceforge.net/lists/listinfo/faudiostream-users>>
> >>>      >   
> >>>    
> >>      <https://lists.sourceforge.net/lists/listinfo/faudiostream-users
> >>    <https://lists.sourceforge.net/lists/listinfo/faudiostream-users>
> >>>    
> >>     <https://lists.sourceforge.net/lists/listinfo/faudiostream-users
> >>    <https://lists.sourceforge.net/lists/listinfo/faudiostream-users>>>
> >>>      >
> >>>      >
> >>>      >
> >>>      > --
> >>>      > "Anybody who knows all about nothing knows everything" --
> >>    Leonard
> >>>      Susskind
> >>> 
> >>> 
> >>> 
> >>> --
> >>> "Anybody who knows all about nothing knows everything" -- Leonard
> >>    Susskind
> >> 
> >> 
> >> 
> >> -- 
> >> "Anybody who knows all about nothing knows everything" -- Leonard Susskind
> > 
> > 
> > _______________________________________________
> > Faudiostream-users mailing list
> > Faudiostream-users@lists.sourceforge.net
> > https://lists.sourceforge.net/lists/listinfo/faudiostream-users
> 
> 
> 
> _______________________________________________
> Faudiostream-users mailing list
> Faudiostream-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/faudiostream-users
> 
> 
> -- 
> "Anybody who knows all about nothing knows everything" -- Leonard Susskind
> _______________________________________________
> Faudiostream-users mailing list
> Faudiostream-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/faudiostream-users



_______________________________________________
Faudiostream-users mailing list
Faudiostream-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/faudiostream-users

Reply via email to