comp_dgt_fw.cc

linear time－frequency toolbox

#include <octave/oct.h>

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

}

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

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

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

  const int N = L/a;

  ComplexMatrix cout(M,N*W*R);  
  
  dgt_fw((double*)f.data(),(double*)gf.data(),L,W,R,a,M,
	 (double*)cout.data());  

  return octave_value (cout);
}