comp_iwfac.cc

linear time－frequency toolbox

#include <octave/oct.h>

extern "C"
{
  void iwfac(double *gf, const int L, const int R,
	     const int a, const int M, double *g);

  int gcd(const int a, const int b,int *r, int *s);

}

DEFUN_DLD (comp_iwfac, args, ,
  "Computes inverse window factorization.\n\
  Usage: c=comp_iwfac(gf,L,a,M);\n\
  Yeah.")
{

  ComplexMatrix gf = args(0).complex_matrix_value();
  const int L = args(1).int_value();
  const int a = args(2).int_value();
  const int M = args(3).int_value();

  const int R=gf.rows()*gf.columns()/L;

  ComplexMatrix g(L,R);
  
  iwfac((double*)gf.data(),L,R,a,M,(double*)g.data());

  return octave_value (g);
}