restore *.orig.c to each entry

Under each yyyy/dir there is a file of the form '*.orig.c'
containing the original source code, or as close as we can
get to the original source code.

2006/stewart/stewart.orig.c

@@ -0,0 +1,65 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#define P printf
+#define F(x,y) ((FILE*)Z(x,y))
+#define U(x,y) ((unsigned char*)Z(x,y))
+#define X(x,y,z) ((double*)x+y+7*z)
+#define Y(x,y) ((double*)x+y)
+#define Z(x,y) *(long*)Y(x,y)
+#define R(x,y) (x=y,y^=(y<<11),y^=(y>>15),y^=(y<<3),(x+y)&0x7fffffff)
+
+int main(int c,char**o)
+{
+   void* v=malloc(0170UL);
+   memset(v,0,0170UL);
+   if(c!=4)P("bad args\n"),exit(1);
+   Z(v,2)=atoi(o[1]);
+   Z(v,1)=atoi(o[2]);
+   Z(v,0)=Z(v,2)*Z(v,2)>>3;
+   Z(v,6)=(long)malloc(Z(v,0));
+   memset((void*)Z(v,6),0,Z(v,0));
+   Z(v,13)=(long)fopen(o[3],"r");
+   fscanf(F(v,13),"%ld",(long*)Y(v,3));
+   Z(v,8)=(long)malloc(050+Z(v,3)*070);
+   Z(Z(v,8),0)=Z(v,3);
+   fscanf(F(v,13),"%lf %lf %lf %lf",Y(Z(v,8),1),Y(Z(v,8),2),Y(Z(v,8),3),Y(Z(v,8),4));
+   for(Z(v,3)=0;Z(v,3)<Z(Z(v,8),0);
+       fscanf(F(v,13),"%lf %lf %lf %lf %lf %lf %lf",X(Z(v,8),5,Z(v,3)),
+              X(Z(v,8),6,Z(v,3)),X(Z(v,8),7,Z(v,3)),X(Z(v,8),8,Z(v,3)),
+              X(Z(v,8),9,Z(v,3)),X(Z(v,8),10,Z(v,3)),X(Z(v,8),11,Z(v,3))),
+           ++Z(v,3));
+   fclose(F(v,13));
+   Z(v,14)=8675309;
+   for(Z(v,3)=0;Z(v,3)<Z(v,1);++Z(v,3))
+   {
+      *Y(v,7)=(R(Z(v,4),Z(v,14))/(double)0x7fffffff),Z(v,5)=0;
+      for(Z(v,4)=0;Z(v,4)<Z(Z(v,8),0);++Z(v,4),++Z(v,5))
+         if(*Y(v,7)<*X(Z(v,8),11,Z(v,5)))
+            break;
+         else
+            *Y(v,7)-=*X(Z(v,8),11,Z(v,5));
+      *Y(v,11)=*X(Z(v,8),5,Z(v,5))**Y(v,9)+*X(Z(v,8),6,Z(v,5))**Y(v,10)+*X(Z(v,8),9,Z(v,5));
+      *Y(v,12)=*X(Z(v,8),7,Z(v,5))**Y(v,9)+*X(Z(v,8),8,Z(v,5))**Y(v,10)+*X(Z(v,8),10,Z(v,5));
+      *Y(v,9)=*Y(v,11);
+      *Y(v,10)=*Y(v,12);
+      if(Z(v,3)>10)
+         if((*Y(v,9)>=*Y(Z(v,8),1))&&(*Y(v,9)<=*Y(Z(v,8),3))&&(*Y(v,10)>=*Y(Z(v,8),2))&&(*Y(v,10)<=*Y(Z(v,8),4)))
+            *Y(v,11)=(*Y(v,9)+-*Y(Z(v,8),1))/(*Y(Z(v,8),3)-*Y(Z(v,8),1)),
+      *Y(v,12)=(*Y(v,10)+-*Y(Z(v,8),2))/(*Y(Z(v,8),4)-*Y(Z(v,8),2)),
+      *(U(v,6)+((((long)(Z(v,2)-*Y(v,12)*Z(v,2))*Z(v,2))+(long)(*Y(v,11)*Z(v,2)))>>3))|=
+      1<<((((long)(Z(v,2)-*Y(v,12)*Z(v,2))*Z(v,2))+(long)(*Y(v,11)*Z(v,2)))&7);
+   }
+   P("#define x_width %ld\n#define x_height %ld\nstatic char x_bits[] = {\n",Z(v,2),Z(v,2));
+   for(Z(v,3)=0;Z(v,3)<Z(v,0);++Z(v,3))
+   {
+      if(Z(v,3))
+         P(",");
+      P("0x%02x",*(U(v,6)+Z(v,3)));
+   }
+   P("};\n");
+   free(Y(Z(v,6),0));
+   free(Y(Z(v,8),0));
+   free(v);
+   return 0;
+}