#include <string.h>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <complex.h>
#include "h.h"
static int rm(int x, int X){return x+x<X?x:x-X;}
static R sq(R x){return x*x;}
static const R pi = 3.14159265358979324;
// The following is official C but missing from Apples complex.h
c CMPLX(R r, R i){c x; __real__ x = r; __imag__ x = i; return x;}
static int X, Y; static R rad[5];
extern c* a;


// No change
static void nop(){int k=X*Y; while(k--) a[k] = ~a[k];};


// Abrupt frequency cutoff (cold turkey)
static void ct(){R sx = rad[0]/X, sy = rad[0]/Y;
  {int k=X*Y; while(k--) a[k] = ~a[k];}
  {int x=X; while(x--){int y=Y; while(y--) {
  if((sq(rm(x, X)*sx) + sq(rm(y, Y)*sy)) > 1) a[y*X + x] = 0;}}}};


// Gaussian blur
static void Gauss(){R sx = rad[0]/X, sy = rad[0]/Y;
{int x=X; while(x--){int y=Y; while(y--) 
   a[y*X + x] = ~a[y*X + x] * exp(-(sq(rm(x, X)*sx) + sq(rm(y, Y)*sy)));}}}


// Convolve with uniform round disk (unfocus)
static void unf(){R s2 = rad[0]*rad[0];
  c * b = (c *)malloc(sizeof(c)*X*Y);
  {int x=X; while(x--){int y=Y; while(y--) {
  b[y*X + x] = ((sq(rm(x, X)) + sq(rm(y, Y))) < s2)?1:0;}}};
  ext2D(b);
  {int j=X*Y; while(j--) a[j] = ~a[j]*b[j];}
  free(b);}


// Auto correlation:
static void ac(){int k=X*Y; while(k--) a[k] *= ~a[k];};


// Vertical translation
static void up(){c ph = 1; R b = 2*pi/Y*rad[0];
  c f = CMPLX(cos(b), sin(b));
  int y; for(y=0; y<Y; ++y) {int Y = y*X;
    {int x; for(x = 0; x<X; ++x) a[Y+x] = ~a[Y+x]*ph;} ph*=f;}}


// Vertical and Horizontal translation, plus double exposure.
static void db(){
  R v = 2*pi/Y*rad[0]; c dv = CMPLX(cos(v), sin(v));
  R h = 2*pi/X*rad[1]; c dh = CMPLX(cos(h), sin(h));
  int y; c py = 1; for(y=0; y<Y; ++y) {int Y = y*X; c px = py;
    {int x; for(x = 0; x<X; ++x) {c*t=a+Y+x;
       *t = ~*t*((1-rad[2])*px + rad[2]); px*=dh;}} py*=dv;}}


void disp(int argcL, char * * argvL, int x, int y){
typedef void manip(void);
typedef struct{char * nm; manip *en;} cp;
X = x; Y = y;
cp d[] = {{"nop", nop}, {"ct", ct}, {"Gauss", Gauss}, {"up", up},
          {"db", db}, {"ac", ac}, {"unf", unf}};
{int k; for(k=3; k<argcL; ++k) rad[k-3] = atof(argvL[k]);}
{int j = sizeof(d)/sizeof(d[0]);
  while(j--) if(!strcmp(argvL[0], d[j].nm)) {d[j].en(); return;}
  exit(printf("No such sub-command as %s\n", argvL[3]));}}