comp_candual_fac.cc

来自「Matlab时频分析工具箱,希望能对大家有所帮助啊」· CC 代码 · 共 37 行

#include <octave/oct.h>
#include "config.h"
#include "../src/ltfat_blas.c"
#include "../src/ltfat_lapack.c"
#include "../src/candual_fac.c"
#include "../src/gcd.c"

DEFUN_DLD (comp_candual_fac, args, ,
  "This function calls the C-library\n\
  gdualf=comp_candual_fac(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;

  const int N = L/a;

  int h_a, h_m;

  const int c=gcd(a, M,&h_a, &h_m);
  const int p=a/c;
  const int q=M/c;
  const int d=N/q;

  ComplexMatrix gdualf(p*q*R,c*d);
  
  candual_fac((ltfat_complex*)gf.data(),L,R,a,M,
	 (ltfat_complex*)gdualf.data());  

  return octave_value (gdualf);
}