lasp_slm.h

Go to the documentation of this file.

#pragma once
#include "lasp_biquadbank.h"
#include "lasp_filter.h"
#include "lasp_thread.h"
#include <memory>
#include <optional>
 
class SLM {
  GlobalThreadPool _pool;
  std::unique_ptr<Filter> _pre_filter;
 
  std::vector<std::unique_ptr<Filter>> _bandpass;
  d _alpha = -1;
  vd _sp_storage;
 
  d Lrefsq;            
  us downsampling_fac; 
  us cur_offset = 0;   
public:
  vd Pm;
  vd Pmax; 
  vd Ppeak;
 
  us N = 0; 
 
  SLM(const d fs, const d Lref, const us downsampling_fac, const d tau,
      std::unique_ptr<Filter> pre_filter,
      std::vector<std::unique_ptr<Filter>> bandpass);
 
  static SLM fromBiquads(const d fs, const d Lref, const us downsampling_fac,
                         const d tau, const vd &pre_filter_coefs,
                         const dmat &bandpass_coefs);
  static SLM fromBiquads(const d fs, const d Lref, const us downsampling_fac,
                         const d tau, const dmat &bandpass_coefs);
 
  ~SLM();
 
  void reset();
 
  SLM(const SLM &o) = delete;
  SLM &operator=(const SLM &o) = delete;
  SLM(SLM &&o) = default;
 
  dmat run(const vd &input);
 
  vd Lpeak() const { return 10 * arma::log10(Ppeak / Lrefsq); };
  vd Leq() const { return 10 * arma::log10(Pm / Lrefsq); };
  vd Lmax() const { return 10 * arma::log10(Pmax / Lrefsq); };
 
  static us suggestedDownSamplingFac(const d fs,const d tw);
 
private:
  vd run_single(vd input, const us filter_no);
};