116 lines
4.7 KiB
C++
116 lines
4.7 KiB
C++
#ifndef MAYER_H
|
|
#define MAYER_H
|
|
|
|
#define REAL float
|
|
|
|
//void mayer_realfft(int n, REAL* real);
|
|
//void mayer_realifft(int n, REAL* real);
|
|
class MayerFFT {
|
|
public:
|
|
void realfft(int n, REAL* real);
|
|
void realifft(int n, REAL* real);
|
|
|
|
MayerFFT() :halsec{
|
|
0,
|
|
0,
|
|
.54119610014619698439972320536638942006107206337801,
|
|
.50979557910415916894193980398784391368261849190893,
|
|
.50241928618815570551167011928012092247859337193963,
|
|
.50060299823519630134550410676638239611758632599591,
|
|
.50015063602065098821477101271097658495974913010340,
|
|
.50003765191554772296778139077905492847503165398345,
|
|
.50000941253588775676512870469186533538523133757983,
|
|
.50000235310628608051401267171204408939326297376426,
|
|
.50000058827484117879868526730916804925780637276181,
|
|
.50000014706860214875463798283871198206179118093251,
|
|
.50000003676714377807315864400643020315103490883972,
|
|
.50000000919178552207366560348853455333939112569380,
|
|
.50000000229794635411562887767906868558991922348920,
|
|
.50000000057448658687873302235147272458812263401372
|
|
},
|
|
costab{
|
|
.00000000000000000000000000000000000000000000000000,
|
|
.70710678118654752440084436210484903928483593768847,
|
|
.92387953251128675612818318939678828682241662586364,
|
|
.98078528040323044912618223613423903697393373089333,
|
|
.99518472667219688624483695310947992157547486872985,
|
|
.99879545620517239271477160475910069444320361470461,
|
|
.99969881869620422011576564966617219685006108125772,
|
|
.99992470183914454092164649119638322435060646880221,
|
|
.99998117528260114265699043772856771617391725094433,
|
|
.99999529380957617151158012570011989955298763362218,
|
|
.99999882345170190992902571017152601904826792288976,
|
|
.99999970586288221916022821773876567711626389934930,
|
|
.99999992646571785114473148070738785694820115568892,
|
|
.99999998161642929380834691540290971450507605124278,
|
|
.99999999540410731289097193313960614895889430318945,
|
|
.99999999885102682756267330779455410840053741619428
|
|
},
|
|
sintab{
|
|
1.0000000000000000000000000000000000000000000000000,
|
|
.70710678118654752440084436210484903928483593768846,
|
|
.38268343236508977172845998403039886676134456248561,
|
|
.19509032201612826784828486847702224092769161775195,
|
|
.09801714032956060199419556388864184586113667316749,
|
|
.04906767432741801425495497694268265831474536302574,
|
|
.02454122852291228803173452945928292506546611923944,
|
|
.01227153828571992607940826195100321214037231959176,
|
|
.00613588464915447535964023459037258091705788631738,
|
|
.00306795676296597627014536549091984251894461021344,
|
|
.00153398018628476561230369715026407907995486457522,
|
|
.00076699031874270452693856835794857664314091945205,
|
|
.00038349518757139558907246168118138126339502603495,
|
|
.00019174759731070330743990956198900093346887403385,
|
|
.00009587379909597734587051721097647635118706561284,
|
|
.00004793689960306688454900399049465887274686668768
|
|
},
|
|
coswrk{
|
|
.00000000000000000000000000000000000000000000000000,
|
|
.70710678118654752440084436210484903928483593768847,
|
|
.92387953251128675612818318939678828682241662586364,
|
|
.98078528040323044912618223613423903697393373089333,
|
|
.99518472667219688624483695310947992157547486872985,
|
|
.99879545620517239271477160475910069444320361470461,
|
|
.99969881869620422011576564966617219685006108125772,
|
|
.99992470183914454092164649119638322435060646880221,
|
|
.99998117528260114265699043772856771617391725094433,
|
|
.99999529380957617151158012570011989955298763362218,
|
|
.99999882345170190992902571017152601904826792288976,
|
|
.99999970586288221916022821773876567711626389934930,
|
|
.99999992646571785114473148070738785694820115568892,
|
|
.99999998161642929380834691540290971450507605124278,
|
|
.99999999540410731289097193313960614895889430318945,
|
|
.99999999885102682756267330779455410840053741619428
|
|
},
|
|
sinwrk{
|
|
1.0000000000000000000000000000000000000000000000000,
|
|
.70710678118654752440084436210484903928483593768846,
|
|
.38268343236508977172845998403039886676134456248561,
|
|
.19509032201612826784828486847702224092769161775195,
|
|
.09801714032956060199419556388864184586113667316749,
|
|
.04906767432741801425495497694268265831474536302574,
|
|
.02454122852291228803173452945928292506546611923944,
|
|
.01227153828571992607940826195100321214037231959176,
|
|
.00613588464915447535964023459037258091705788631738,
|
|
.00306795676296597627014536549091984251894461021344,
|
|
.00153398018628476561230369715026407907995486457522,
|
|
.00076699031874270452693856835794857664314091945205,
|
|
.00038349518757139558907246168118138126339502603495,
|
|
.00019174759731070330743990956198900093346887403385,
|
|
.00009587379909597734587051721097647635118706561284,
|
|
.00004793689960306688454900399049465887274686668768
|
|
}
|
|
|
|
{
|
|
}
|
|
private:
|
|
|
|
void fht(REAL* fz, int n);
|
|
REAL halsec[20];
|
|
REAL costab[20];
|
|
REAL sintab[20];
|
|
REAL coswrk[20];
|
|
REAL sinwrk[20];
|
|
};
|
|
|
|
#endif |