M_rdft
M_rdft(nn, 1, a);
M_rdft(nn, 1, b);
M_rdft(nn, -1, b);
extern void M_rdft(int, int, double *);