1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
#include "ltfat.h"
#include "ltfat_types.h"
LTFAT_EXTERN void
LTFAT_NAME(idgt_long)(const LTFAT_COMPLEX *cin, const LTFAT_COMPLEX *g,
const ltfatInt L, const ltfatInt W,
const ltfatInt a, const ltfatInt M,
const dgt_phasetype ptype,
LTFAT_COMPLEX *f)
{
LTFAT_COMPLEX *gf = ltfat_malloc(L * sizeof * gf);
LTFAT_NAME_COMPLEX(wfac)(g, L, 1, a, M, gf);
LTFAT_NAME(idgt_fac)(cin, (const LTFAT_COMPLEX*) gf, L, W, a, M, ptype, f);
ltfat_free(gf);
}
LTFAT_EXTERN void
LTFAT_NAME(idgtreal_long)(const LTFAT_COMPLEX *cin, const LTFAT_REAL *g,
const ltfatInt L, const ltfatInt W,
const ltfatInt a, const ltfatInt M,
const dgt_phasetype ptype, LTFAT_REAL *f)
{
// TO DO: Is it possible to use wfacreal? wfacreal_size(L,a,M)
LTFAT_COMPLEX *gf = ltfat_malloc(L * sizeof * gf);
LTFAT_NAME(wfac)(g, L, 1, a, M, gf);
LTFAT_NAME(idgtreal_fac)(cin, (const LTFAT_COMPLEX*) gf, L, W, a, M,
ptype, f);
ltfat_free(gf);
}
|
