#include <stdio.h>
#include "misc.h"
#include "fft.h"
#include "rnd.h"
#include <math.h>
#define P (1<<17)
c d[P];
int main(){
R const tp = 2.*3.1415926535897932, r2 = 1./sqrt(2);
{int j=P; while(j--) d[j] = cmp(r2*rnd(), r2*rnd()) + 0.02*cos(60436*tp*j/P);}
   fft(d, P);
   {int j=P; while(j--) {R m = m2(d[j]);
      if (m > 1130000) {printf("j=%7d m2=%e ", j, m); pc(d[j]);}}}
return 0;}