#include <stdio.h>
#define n1 4096
#define n2 4096
float a[n2*n1],b[n2*n1],c[n2*n1];
char head[72][80],head0[72][80];
FILE *fp;
main(int ac, char **av)
{
  int i,j,k;
  j=n1*n2;
  if(ac<4){ printf("\n\t *** 4k_float_fit div 4k_float *****\n");
            printf("\n\t div_fit f1 f2 f3   (f3=f1/f2)\n\n");
            exit(0); }
  fp=fopen(av[1],"rb"); if(fp==0){ printf("%s not found!\n",av[1]); exit(0);}
  fread(head,72,80,fp); k=fread(a,4,j,fp);
  if(k!=j){ printf("%s not a 4k*4k_float file!\n",av[1]); exit(0); } 
  swap4(a,4*k); fclose(fp);
  fp=fopen(av[2],"rb"); if(fp==0){ printf("%s not found!\n",av[2]); exit(0);}
  fread(head0,72,80,fp); k=fread(b,4,j,fp);
  if(k!=j){ printf("%s not a 4k*4k_float file!\n",av[2]); exit(0); } 
  swap4(b,4*k); fclose(fp);
  for(i=0;i<k;i++)c[i]=a[i]/b[i];
  fp=fopen(av[3],"wb"); fwrite(head,72,80,fp); swap4(c,4*k); 
  fwrite(c,4,j,fp); fclose(fp);
}