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