home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
The C Users' Group Library 1994 August
/
wc-cdrom-cusersgrouplibrary-1994-08.iso
/
vol_200
/
245_01
/
lca31.c
< prev
next >
Wrap
Text File
|
1987-10-28
|
16KB
|
576 lines
/* (3,1) Linear Cellular Automaton */
/* Reference: */
/* */
/* Kenneth E. Perry */
/* Abstract Mathematical Art */
/* BYTE */
/* December, 1986 */
/* pages 181-192 */
/* Copyright (C) 1987 */
/* Harold V. McIntosh */
/* Gerardo Cisneros S. */
/* G. Cisneros, 4.3.87 */
/* 10 April 1987 - modified for (4,2) [HVM] */
/* 26 April 1987 - Multiple menus [HVM] */
/* 28 April 1987 - back modified to (4,1) [HVM] */
/* 28 April 1987 - version for XVI Feria de Puebla [HVM] */
/* 14 May 1987 - modified for (3,1) and general rule [HVM] */
# include <bdos.h>
# define COLGRAF 4 /* graph resolution */
# define T80X25 3 /* text resolution */
# define WHCYMAG 1 /* color quad for normal screen */
# define AL 320 /* array length (screen width) */
# define TS 7 /* distinct sums w/totalistic rule */
# define DS 27 /* number of distinct sums */
# define KK 3 /* number of states per cell */
# define NX 21 /* number of sample rules */
char xrule[NX][KK][KK][KK];
char ixrule[NX][DS]=
"001200111011011222201121100", /* interfaces of 2 vel */
"001200211221001010221220120",
"001210120212011110020110100", /* mottled background */
"002110012010121101201012211", /* (0*102211200)* still */
"002110012110021201001012221", /* symmetrized */
"002200011112011011012022200", /* threads on triangles */
"010002101012100210102011211", /* dendrites */
"010020202121101100222101120", /* gldrs, puffr trn */
"010022101200112021010000110", /* threads among trngs */
"010102020102020202020202020",
"010102022102022221022221210", /* open triangles */
"011200122100121212211112212",
"012100200110011022102202022", /* 2 comes and goes */
"012112212012012112012212210", /* 12 irresistible gldr */
"020101011101011112011112120", /* diagonal black gaps */
"020122100021210222210022111", /* black triangles */
"020201010201010102010102020", /* class iv, sparse */
"021111222210112122112122220", /* 102 ia a barrier */
"021211110211110102110102020", /* blue bgrnd cl iv */
"210021220202112222100220222", /* 0122 barrier */
"221212021202011210012120111" /* 2 glider on 1 bkgrnd */
;
char xx[4], rule[DS+1], ascrule[KK][KK][KK];
int binrule[KK][KK][KK], arule[DS], arr1[AL], arr2[AL];
char trule[TS]="0000000";
main()
{
int i, j, i0, i1, i2, jj, n;
int more = 'r';
char a, b, c;
for (i=0; i<NX; i++) { /* copy to 3-index array */
i0=0; i1=0; i2=0;
for (j=0; j<DS; j++) {
xrule[i][i0][i1][i2]=ixrule[i][j];
i2++;
if (i2>2) {i2=0; i1++;};
if (i1>2) {i1=0; i0++;};
if (i0>2) {i2=0; };
};};
videopalette(WHCYMAG); /* white/cyan/magenta */
tuto();
while (!kbdst()) jj=rand(); /* wait for keypress */
c=kbdin(); /* ignore it */
jj=rand()%NX;
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
ascrule[i0][i1][i2] = xrule[jj][i0][i1][i2]; /* random sample rule */
};};};
rule[DS]=0;
ranlin(); /* random initial array */
videomode(T80X25);
videoscroll(3,0,5,71,0,3); /* menu on blue background */
videoscroll(19,0,24,71,0,3);
xmenu(jj+1); /* show initial rule */
while (more!='n') { /* execute multiple runs */
rmenu();
lmenu();
while (0<1) { /* set up one run */
c=kbdin();
if (c=='g') break; /* go draw graph */
if (c=='q') more='n'; /* quit for good */
if (more=='n') break;
switch (c) {
case '@': /* numbered tot rule */
nutoto(numin(0));
totoasc();
rmenu();
videocursor(0,4,0);
break;
case '$': /* dozen totalistics */
j=numin(0);
for (i=0; i<12; i++) {
nutoto(j+i);
totoasc();
ranlin();
evolve(rule);
};
videomode(T80X25);
rmenu();
lmenu();
break;
case 't': /* totalistic rule */
xblnk();
tmenu();
edtrule();
totoasc();
for (i0=0; i0<3; i0++) {
for (i1=0; i1<3; i1++) {
for (i2=0; i2<3; i2++) {
ascrule[i0][i1][i2]=trule[i0+i1+i2];
};};};
videocursor(0,4,0);
rmenu();
xmenu(totonu(0));
break;
case 'r': /* edit rule */
xblnk();
edrule();
videocursor(0,4,0);
rmenu();
break;
case 'l': /* edit cell string */
xblnk();
edline();
videocursor(0,3,0);
lmenu();
break;
case '#': /* read stored rule */
xmenu(NX);
n=(i=lim(1,numin(0),NX))-1;
xmenu(i);
xtoasc(n);
rmenu();
break;
case 'D': /* run through samples */
for (i=0; i<NX; i++) {
xmenu(i);
xtoasc(i);
ranlin();
evolve(rule);
};
videomode(T80X25);
rmenu();
break;
case 'u': /* sparse init arry */
xblnk();
for (i=0; i<AL; i++) arr1[i]=0;
arr1[AL/4]=1;
arr1[AL/2]=2;
arr1[(3*AL)/4]=1;
arr1[(3*AL)/4+2]=2;
lmenu();
break;
case 'x': /* random rule */
xblnk();
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
if ((KK*(KK*i0+i1)+i2)%4 == 0) i=rand();
ascrule[i0][i1][i2]='0'+i%3;
i/=3;
};};};
rmenu();
break;
case 'y': /* random line */
xblnk();
ranlin();
lmenu();
break;
case 'Y': /* symmetrize rule */
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
ascrule[i2][i1][i0]=ascrule[i0][i1][i2];
};};};
break;
case 'B': /* begin barrier */
a=kbdin();
b=kbdin();
ascrule[0][a-'0'][b-'0']=a;
ascrule[1][a-'0'][b-'0']=a;
ascrule[2][a-'0'][b-'0']=a;
rmenu();
break;
case 'E': /* end barrier */
a=kbdin();
b=kbdin();
ascrule[a-'0'][b-'0'][0]=b;
ascrule[a-'0'][b-'0'][1]=b;
ascrule[a-'0'][b-'0'][2]=b;
rmenu();
break;
case 'L': /* left glider link */
a=kbdin();
b=kbdin();
c=kbdin();
ascrule[a-'0'][b-'0'][c-'0']=c;
rmenu();
break;
case 'R': /* left glider link */
a=kbdin();
b=kbdin();
c=kbdin();
ascrule[a-'0'][b-'0'][c-'0']=a;
rmenu();
break;
case 'S': /* still life link */
a=kbdin();
b=kbdin();
c=kbdin();
ascrule[a-'0'][b-'0'][c-'0']=b;
rmenu();
break;
case '=':
for (i=1; i<8; i++) {
for (j=0; j<40; j++) arr1[40*i+j]=arr1[j];};
lmenu();
break;
case '~':
for (i=1; i<16; i++) {
for (j=0; j<20; j++) arr1[20*i+j]=arr1[j];};
lmenu();
break;
default: break;
};
};
if (more=='n') break;
do {
evolve(rule);
videocursor(0,0,0);
scrstr("More?");
videocursor(0,0,34);
scrstr("y/n/cr");
more=kbdin();
} while (more=='\015');
videomode(T80X25); /* reset the screen */
if (more=='n') break;
};
videomode(T80X25);}
/* edit the rule */
edrule() {
char c;
int i, i0, i1, i2;
i=6; i0=0; i1=0, i2=0;
while (0<1) {
videocursor(0,3,i);
c = kbdin();
if (c == '\015') break; /* carriage return exits */
switch (c) {
case '0': case '1': case '2': /* state */
ascrule[i0][i1][i2] = c;
i2++;
if (i2>2) {i2=0; i1++;};
if (i1>2) {i1=0; i0++;};
if (i0>2) {i2=0; };
videocattr(0,c,3,1);
if (i<6+DS) i++;
break;
case ' ': case '\315': /* space = advance */
i2++;
if (i2>2) {i2=0; i1++;};
if (i1>2) {i1=0; i0++;};
if (i0>2) {i2=0; };
if (i<6+DS) i++;
break;
case '\010': case '\313': /* backspace */
if (i2!=0) i2--; else {i2=2;
if (i1!=0) i1--; else {i1=2;
if (i0!=0) i0--; else {i0=2;
};};};
if (i>6) i--;
break;
default: break;
};
};
}
/* edit totalistic rule */
edtrule() {char c; int i;
i=0;
while (i<TS) {
c=trule[i];
videocursor(0,3,56+i);
videocattr(0,c,3,1);
c = kbdin();
if (c == '\015') break;
switch (c) {
case '0': case '1': case '2': /* state */
trule[i]=c;
videocattr(0,c,5,1);
i++;
break;
case ' ': case '\315': /* space = advance */
i++;
break;
case '\010': case '\313': /* backspace */
if (i!=0) i--;
break;
default: break;
};
};
}
/* edit the line */
edline() {
char c;
int i, j, k, ii, jj;
videocursor(0,19,0);
scrstr("insert states using 0, 1, 2");
videocursor(0,20,0);
scrstr("move cursor with n(north), s(south), e(east), w(west), or");
videocursor(0,21,0);
scrstr("with keyboard arrows. Space, backspace move right and left.");
videocursor(0,22,0);
scrstr("( seeks left margin, < absolutely, { up one line, [ down one line");
videocursor(0,23,0);
scrstr(") seeks right margin, > absolutely, } up one line, ] down one line");
videocursor(0,24,0);
scrstr("carriage return exits");
jj=4;
ii=1;
while (0<1) {
ii=lim(1,ii,40);
jj=lim(1,jj,8);
j=jj-1;
i=ii-1;
videocursor(0,j+9,i);
c=kbdin();
if (c == '\015') {videoscroll(19,0,24,70,0,3); break;};
switch (c) {
case '0': case '1': case '2': /* enter state */
arr1[40*j+i]=c-'0'; ii++; break;
case 's': case '\012': case '\320': jj++; break; /* down - next line */
case 'n': case '\013': case '\310': jj--; break; /* up - last line */
case 'e': case '\014': case '\315': case ' ': ii++; break; /* space = advance */
case 'w': case '\010': case '\313': ii--; break; /* backspace */
case '<': ii=1; jj=1; break; /* absolute left */
case '{': ii=1; jj--; break; /* left one row up */
case '(': ii=1; break; /* left this row */
case '[': ii=1; jj++; break; /* left next row */
case '>': ii=40; jj=40; break; /* absolute right */
case '}': ii=40; jj--; break; /* right one row up */
case ')': ii=40; break; /* right this row */
case ']': ii=40; jj++; break; /* right next row */
default: break;
};
videocursor(0,j+9,0);
for (k=0; k<40; k++) videoputc('0'+arr1[40*j+k],1);
};
}
/* display a screen of evolution */
evolve(rule) char *rule; {
int i0, i1, i2, i, j;
videomode(COLGRAF); /* erase the screen */
videocursor(0,0,0); /* top text line */
scrstr("Rule: ");
scrrul();
for (i=0; i<DS; i++) {arule[i] = rule[i]-'0';};
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
binrule[i0][i1][i2]=ascrule[i0][i1][i2]-'0';
};};};
for (j=8; j<200; j++) videodot(j,AL-1,2);
for (j=8; j<200; j++) { /* evolve for 192 generations */
arr2[0] = binrule[arr1[AL-1]][arr1[0]][arr1[1]];
for (i=1; i<AL-1; i++) {
arr2[i] = binrule[arr1[i-1]][arr1[i]][arr1[i+1]];
};
arr2[AL-1] = binrule[arr1[AL-2]][arr1[AL-1]][arr1[0]];
for (i=0; i<AL; i++) {videodot(j,i,arr1[i]); arr1[i] = arr2[i];};
if (kbdst()) {kbdin(); break;};
}
}
/* generate a random line of cells in arr1 */
ranlin() {int i, c;
for (i=0; i<AL; i++) {
if (i%8 == 0) c=rand();
arr1[i]=c%3; c/=3;};
}
/* tutorial and Help screen */
tuto() {
videomode(T80X25);
videocursor(0,2,0);
scrstr("<Copyright (C) 1987 - H.V. McIntosh, G. Cisneros S.>");
videocursor(0,4,0);
scrstr(" *** LIFE in One Dimension ***");
videocursor(0,6,0);
scrstr("Three States - Black(0), Cyan(1), Magenta(2).");
videocursor(0,8,0);
scrstr("First neighbors - one on each side, three altogether.");
videocursor(0,10,0);
scrstr("Complete transition rule - random, edited, or stored.");
videocursor(0,12,0);
scrstr("Initial Cellular Array - random, edited, or patterned.");
videocursor(0,14,0);
scrstr("Some rules are fragile and require several tries before");
videocursor(0,15,0);
scrstr("manifesting an interesting evolutionary pattern.");
videocursor(0,17,0);
scrstr("Use any key to terminate a display in progress.");
videocursor(0,21,0);
scrstr("now, press any key to continue ...");
}
/* rule menu */
rmenu() {
videocursor(0,0,0);
scrstr(" 0........1........2........");
videocursor(0,1,0);
scrstr(" 0..1..2..0..1..2..0..1..2..");
videocursor(0,2,0);
scrstr(" 012012012012012012012012012");
videocursor(0,3,0);
scrstr("Rule: ");
scrrul();
if (istot()==1) {scrstr(" "); printf("%5d",totonu(0));};
videocursor(0,5,0);
scrstr(" r(rule), l(line), #nn(stored rule), g(graph), q(quit),");
videocursor(0,6,0);
scrstr(" x(random rule), y(random line), u(unit line),");
videocursor(0,7,0);
scrstr(" @nn(tot/rule), $nn(dzn tot/rules), t(ed tot/rule).");
videocursor(0,5,0);
}
/* totalistic rule menu*/
tmenu() {
videocursor(0,2,50);
scrstr(" 0..1..2");
videocursor(0,3,50);
scrstr("rule: ");
tscrrul();
videocursor(0,3,56);
}
/* line menu */
lmenu() {int i, j;
for (j=0; 40*j<AL; j++) {
videocursor(0,9+j,0);
for (i=0; i<40; i++) videoputc('0'+arr1[40*j+i],1);
}
videocursor(0,5,0); }
/* display rule number */
xmenu(n) int n;
{int i, nn;
nn=sprintf(xx,"%3d",n);
videocursor(0,1,40);
for (i=0; i<nn; i++) videoputc(xx[i],1);
videocursor(0,1,40); }
/* clear screen space for rule number */
xblnk() {
videocursor(0,1,40);
scrstr(" ");
videocursor(0,5,0); }
/* copy saved rule #n into active rule */
xtoasc(n) int n; { int i0, i1, i2;
xmenu(n);
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
ascrule[i0][i1][i2] = xrule[n][i0][i1][i2]; /* random sample rule */
};};}; }
/* change totalistic rule to general rule */
totoasc() {
int i0, i1, i2;
for (i0=0; i0<3; i0++) {
for (i1=0; i1<3; i1++) {
for (i2=0; i2<3; i2++) {
ascrule[i0][i1][i2]=trule[i0+i1+i2];
};};};
}
/* change decimal totalistic rule to sum values */
nutoto(x) int x; {int i;
for (i=0; i<TS; i++) {trule[i]=x%3+'0'; x/=3;};}
/* change sum values to decimal totalistic rule */
int totonu(i) int i; {int r;
if (i<TS) r=(trule[i]-'0')+3*totonu(i+1); else r=0;
return r; }
/* test whether a rule is totalistic */
int istot() {int i0, i1, i2, l;
l=1;
trule[0]=ascrule[0][0][0];
trule[1]=ascrule[1][0][0];
trule[2]=ascrule[2][0][0];
trule[3]=ascrule[2][1][0];
trule[4]=ascrule[2][2][0];
trule[5]=ascrule[2][2][1];
trule[6]=ascrule[2][2][2];
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
if (ascrule[i0][i1][i2]!=trule[i0+i1+i2]) l=0;
};};};
return l; }
/* limit j to interval (i,k) */
int lim(i,j,k) int i, j, k;
{if (i>=j) return i; if (k<=j) return k; return j;}
/* display the rule number on the screen */
scrrul() {int i0, i1, i2;
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
videoputc(ascrule[i0][i1][i2],1);
};};}; }
/* display totalistic rule number by sum */
tscrrul() {int i;
for (i=0; i<TS; i++) {videoputc(trule[i]); }; }
/* write a string in graphics mode */
scrstr(s) char *s;
{for (; *s != '\0'; s++) videoputc(*s,1); }
/* keyboard status */
int kbdst() {return(bdos(11) & 0xFF);}
/* direct keyboard input, no echo */
kbdin() {int c;
if ((c = bdos(7) & 0xFF) == '\0') c = (bdos(7) & 0xFF) | 0x80;
return(c);
}
/* read number from keyboard */
int numin(n) int n; {char c;
c=kbdin();
if (c>='0'&&c<='9') return(numin(10*n+(c-'0'))); else return(n);
}
/* end */