comp_idgt_fw.cc

来自「linear time－frequency toolbox」· CC 代码 · 共 32 行

#include <octave/oct.h>

extern "C"
{
  void idgt_fw(double *cin, double *gf, const int L, const int W,
	      const int R, const int a, const int M, double *f);

}

DEFUN_DLD (comp_idgt_fw, args, ,
  "This function calls the C-library\n\
  c=comp_idgt_fw(c,gf,L,a,M);\n\
  Yeah.\n")
{

  ComplexMatrix coef = args(0).complex_matrix_value();
  ComplexMatrix gf = args(1).complex_matrix_value();
  const int L = args(2).int_value();
  const int a = args(3).int_value();
  const int M = args(4).int_value();
  const int N = L/a;
  const int R = gf.rows()*gf.columns()/L;
  const int W = coef.rows()*coef.columns()/(M*N*R);

  ComplexMatrix f(L,W);  
  
  idgt_fw((double*)coef.data(),(double*)gf.data(),L,W,R,a,M,
	  (double*)f.data());  

  return octave_value (f);
}