home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
The Fred Fish Collection 1.5
/
ffcollection-1-5-1992-11.iso
/
ff_disks
/
001-099
/
ff005.lzh
/
MandelbrotSrc
/
mand.c
< prev
next >
Wrap
C/C++ Source or Header
|
1986-01-11
|
22KB
|
723 lines
/*
MAND.C - Command parsing
Mandelbrot's Self-Squared Dragon Generator
For the Commodore Amiga
Version 2.01
Inspired by Scientific American, August/1985
Corrections and improvements suggested by
The Fractal Geometry of Nature
By Benoit Mandelbrot, W.H. Freeman and Company, 1983
(Used to be Z=Z^2+C, now is Z=Z^2-u, etc.)
Copyright (C) 1985, Robert S. French
Vastly Enhanced by =RJ Mical= 1985/86
Copyright (C) 1986, =Robert J. Mical=
Placed in the Public Domain
This program may be distributed free of charge as long as the above
notice is retained.
Please see the accompanying documentation for more information.
Send ANY problems or suggestions to the address in the <I>nformation
section. Thank you!
Programs should be compiled with:
1> lc -i:include/ -i:include/lattice/ mand.c
1> lc -i:include/ -i:include/lattice/ mand1.c
1> lc -i:include/ -i:include/lattice/ mand2.c
1> lc -i:include/ -i:include/lattice/ mand3.c
1> lc -i:include/ -i:include/lattice/ mand4.c
1> lc -i:include/ -i:include/lattice/ mand5.c
1> lc -i:include/ -i:include/lattice/ mand6.c
1> alink :lib/Lstartup.obj+mand.o+mand1.o+mand2.o+mand3.o+mand4.o+mand5.o+mand6.o to Mandelbrot lib :lib/lc.lib+:lib/amiga.lib
*/
/* ===========================================================================
===========================================================================
From: "french robert%d.mfenet"@LLL-MFE.ARPA
By the way, on a historical note, the "Mandelbrot Set" as related in
Scientific American is a little different from the u-Map implemented
here (as discussed in Mandelbrot's "The Fractal Geometry of Nature").
The main difference is the use of Z <- Z^2+C in Scientific American vs.
Z <- Z^2-u in the book and this program. The primary effect is the reversal
of the picture across the Y axis.
Enjoy!
Robert French
French#Robert%d@LLL-MFE.ARPA
========================================================================
======================================================================== */
#include "mand.h"
int MathBase,MathTransBase;
/*----------------------*/
/* Graphics definitions */
struct GfxBase *GfxBase;
struct IntuitionBase *IntuitionBase;
struct IconBase *IconBase;
SHORT Color0, Color1, Color2;
SHORT UserPalette[29];
/*----------------------------------*/
/* Miscellaneous Global Definitions */
union kludge {
float f;
int i;
} start_r,end_r,start_i,end_i;
int max_x,max_y,max_mem_y;
int max_count,color_inc,color_offset,color_set,color_mode,color_div;
int color_inset,func_num;
int v_starty,max_mem;
long v_offset;
UWORD *color_table,*v_mand_store;
int modified,want_read;
FILE *console,*v_fp = NULL,*redir_fp = NULL;
SHORT ZoomCenterX, ZoomCenterY, ZoomBoxSizeX, ZoomBoxSizeY;
SHORT ZoomBoxStartX, ZoomBoxStartY;
int cur_resource = NULL;
/*-------------*/
/* Here we go! */
main()
{
FILE *fopen();
char *fgets(),*stpblk();
float cnvf();
void anal_mand(),wait_close();
int temp, y, i;
char *cmd,inp_buf[80],secchar,*argpos;
union kludge scale;
FILE *fp;
max_mem_y = MAXMY;
color_table = (UWORD *)malloc(4096*sizeof(UWORD));
if (color_table == NULL)
abort("Can't allocate memory for color table");
cur_resource |= F_COLORTAB;
v_mand_store = (UWORD *)malloc(MAXX*max_mem_y*sizeof(UWORD));
if (v_mand_store == NULL)
abort("Can't allocate memory for set storage");
cur_resource |= F_SETSTORE;
MathBase = OpenLibrary("mathffp.library",0);
if (MathBase < 1)
abort("Can't open mathffp.library");
cur_resource |= F_MATH;
MathTransBase = OpenLibrary("mathtrans.library",0);
if (MathTransBase < 1)
abort("Can't open mathtrans.library");
cur_resource |= F_MATHTRANS;
GfxBase = (struct GfxBase *)OpenLibrary("graphics.library",0);
if (GfxBase == NULL)
abort("Can't open graphics.library");
cur_resource |= F_GRAPHICS;
IntuitionBase = (struct IntuitionBase *)OpenLibrary("intuition.library",0);
if (IntuitionBase == NULL)
abort("Can't open intuition.library");
cur_resource |= F_INTUITION;
IconBase = (struct IconBase *)OpenLibrary("icon.library",0);
if (IconBase == NULL)
abort("Can't open icon.library");
cur_resource |= F_ICON;
console = fopen("con:0/0/640/200/Mandelbrot Commands","w+");
if (console == NULL)
abort("Can't open console window");
cur_resource |= F_CONSOLE;
fprintf(console,
"Mandelbrot Self-Squared Dragon Generator - Version %s\n",VERSION);
fputs("Copyright (C) 1985, Robert S. French\n",console);
fputs(
"Vastly Enhanced (with Intuition!) by =RJ Mical= 1985/86\n",console);
fputs("Copyright (C) 1986, =Robert J. Mical=\n",console);
fputs("\n",console);
fputs(
"Placed in the public domain -- Read the source and learn!\n",console);
fputs("\n",console);
fputs("Type '?' or 'H' or <RETURN> for help\n\n",console);
SetPresets(0);
Color0 = 0x22F;
Color1 = 0xFFF;
Color2 = 0x000;
for (i = 0; i <= 28; i++)
UserPalette[i] = ((i & 0xF) << 8) | ((i & 0xF) << 4) | ((i & 0xF) << 0);
for (EVER) {
command:
fprintf(console,"Command: ");
rewind(console);
if (redir_fp)
if (Chk_Abort()) {
fputs("\nRedirected input aborted!\n",console);
fclose(redir_fp);
redir_fp = NULL;
}
else {
cmd = fgets(inp_buf,78,redir_fp);
if (cmd == NULL) {
fputs("End of file encountered.\n",console);
fclose(redir_fp);
redir_fp = NULL;
cmd = fgets(inp_buf,40,console);
}
else
fputs(inp_buf,console);
}
else
cmd = fgets(inp_buf,40,console);
*(cmd+strlen(cmd)-1) = '\0';
secchar = toupper(*(cmd+1));
argpos = stpblk(cmd+2);
rewind(console);
switch (toupper(*cmd)) {
case 'I':
sscanf(argpos,"%d",&temp);
Information(temp);
goto command;
case 'S':
if (secchar == 'R') {
sscanf(argpos,"%f",&start_r.f);
start_r.i = SPFieee(start_r.i);
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
if (secchar == 'I') {
sscanf(argpos,"%f",&start_i.i);
start_i.i = SPFieee(start_i.i);
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
if (secchar == 'H') {
fputs("Current settings:\n",console);
CurrentSettings();
goto command;
}
if (secchar == 'A') {
fp = fopen(argpos,"w");
if (fp == NULL) {
fprintf(console,"Cannot open file '%s'\n",argpos);
goto command;
}
putc((char) 1,fp);
fwrite(&start_r,sizeof(start_r),1,fp);
fwrite(&start_i,sizeof(start_i),1,fp);
fwrite(&end_i,sizeof(end_i),1,fp);
fwrite(&max_x,sizeof(max_x),1,fp);
fwrite(&max_y,sizeof(max_y),1,fp);
want_read = TRUE;
for (y=0;y<max_y;y++) {
v_pos_line(y);
if(fwrite(v_mand_store+(y-v_starty)*max_x,
max_x*sizeof(UWORD),1,fp) != 1) {
fputs("Error while writing to file\n",console);
break;
}
}
fclose(fp);
goto command;
}
ill_cmd();
goto command;
case 'E':
if (secchar == 'R') {
sscanf(argpos,"%f",&end_r.i);
end_r.i = SPFieee(end_r.i);
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
if (secchar == 'I') {
sscanf(argpos,"%f",&end_i.i);
end_i.i = SPFieee(end_i.i);
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
ill_cmd();
goto command;
case 'M':
if (secchar == 'X') {
sscanf(argpos,"%d",&temp);
if (temp < 5 ||
(temp > 320 && !(color_mode & 4)) ||
(temp > 640 && (color_mode & 4))) {
fputs("Illegal parameter!\n",console);
goto command;
}
max_x = temp;
max_mem = max_mem_y * MAXX;
max_mem /= max_x;
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
if (secchar == 'Y') {
sscanf(argpos,"%d",&temp);
if (temp < 5 ||
(temp > 200-STARTY && !(color_mode & 2)) ||
(temp > 400-STARTY && (color_mode & 2))) {
fputs("Illegal parameter!\n",console);
goto command;
}
max_y = temp;
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
if (secchar == 'C') {
sscanf(argpos,"%d",&temp);
if (temp * color_inc + color_offset > 4095) {
fputs("More than 4096 colors!\n",console);
goto command;
}
if (temp < 2) {
fputs("MaxCount must be greater than 1\n",console);
goto command;
}
max_count = temp;
goto command;
}
if (secchar == 'M') {
sscanf(argpos,"%d",&temp);
if (temp < 5) {
fputs("Number of lines must be >4!\n",console);
goto command;
}
free(v_mand_store);
v_mand_store = (UWORD *)malloc(MAXX*temp*sizeof(UWORD));
if (v_mand_store == NULL) {
fputs("Can't allocate that much memory for set storage",console);
v_mand_store = (UWORD *)malloc(MAXX*max_mem_y*sizeof(UWORD));
if (v_mand_store == NULL)
abort("Can't reallocate memory!!!");
goto command;
}
max_mem_y = temp;
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
max_mem = MAXX * max_mem_y;
max_mem /= max_x;
goto command;
}
ill_cmd();
goto command;
case 'L':
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
v_fp = fopen(stpblk(cmd+1),"r");
if (v_fp == NULL) {
fprintf(console,"Cannot open file '%s'\n",stpblk(cmd+1));
goto command;
}
if (getc(v_fp) != 1) {
fputs("File is not in proper format\n",console);
fclose(v_fp);
v_fp = NULL;
goto command;
}
fread(&start_r,sizeof(start_r),1,v_fp);
fread(&end_r,sizeof(end_r),1,v_fp);
fread(&start_i,sizeof(start_i),1,v_fp);
fread(&end_i,sizeof(end_i),1,v_fp);
fread(&max_x,sizeof(max_x),1,v_fp);
fread(&max_y,sizeof(max_y),1,v_fp);
v_offset = 25L;
modified = FALSE;
v_starty = max_mem+1;
want_read = TRUE;
v_pos_line(0);
goto command;
case 'X':
if (secchar == 'R') {
sscanf(argpos,"%f",&scale.f);
scale.i = SPFieee(scale.i);
start_r.i = SPAdd(start_r.i,scale.i);
end_r.i = SPAdd(end_r.i,scale.i);
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
if (secchar == 'I') {
sscanf(argpos,"%f",&scale.f);
scale.i = SPFieee(scale.i);
start_i.i = SPAdd(start_i.i,scale.i);
end_i.i = SPAdd(end_i.i,scale.i);
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
goto command;
}
ill_cmd();
goto command;
case 'Z':
if (secchar != 'R' && secchar != 'I' && secchar != 'B') {
ill_cmd();
goto command;
}
if (v_fp) {
fclose(v_fp);
v_fp = NULL;
}
if (secchar != 'I' ) {
/* scale along the real axis */
sscanf(argpos,"%f",&scale.f);
ZoomAlongDarling(SPFieee(scale.i), SPFlt(1));
}
if (secchar != 'R') {
/* scale along the complex axis */
sscanf(argpos,"%f",&scale.f);
ZoomAlongDarling(SPFlt(1), SPFieee(scale.i));
}
goto command;
case 'C':
if (secchar == 'I') {
sscanf(argpos,"%d",&temp);
if (max_count * temp + color_offset > 4095) {
fputs("More than 4096 colors!\n",console);
goto command;
}
if (temp < 1) {
fputs("Increment must be greater than 0!\n",console);
goto command;
}
color_inc = temp;
goto command;
}
if (secchar == 'O') {
sscanf(argpos,"%d",&temp);
if (max_count * color_inc + temp > 4095) {
fputs("More than 4096 colors!\n",console);
goto command;
}
if (temp < 0) {
fputs("Negative offset illegal!\n",console);
goto command;
}
color_offset = temp;
goto command;
}
if (secchar == 'S') {
sscanf(argpos,"%d",&temp);
if (temp < 0 || temp > 2) {
fputs("Illegal color set!\n",console);
goto command;
}
color_set = temp;
init_colors();
goto command;
}
if (secchar == 'M') {
sscanf(argpos,"%d",&temp);
if (temp < 0 || temp > 7 || ((temp&4) && !(temp&1))) {
fputs("Illegal graphics mode!\n",console);
goto command;
}
color_mode = temp;
if (!(color_mode & 0x2))
if (max_y > 200-STARTY)
max_y = 200-STARTY;
if (!(color_mode & 0x4))
if (max_x > 320)
max_x = 320;
goto command;
}
if (secchar == 'D') {
sscanf(argpos,"%d",&temp);
if (temp < 1) {
fputs("Divisor must be greater than 0!\n",console);
goto command;
}
color_div = temp;
goto command;
}
if (secchar == 'T') {
sscanf(argpos,"%d",&temp);
if (temp < 0 || temp > 4095) {
fputs("Color must be between 0 and 4095!\n",console);
goto command;
}
color_inset = temp;
goto command;
}
ill_cmd();
goto command;
case 'F':
sscanf(stpblk(cmd+1),"%d",&temp);
if (temp < 0 || temp > 1) {
fputs("Function number must be 0 or 1!\n",console);
goto command;
}
func_num = temp;
goto command;
case 'P':
sscanf(stpblk(cmd+1),"%d",&temp);
SetPresets(temp);
goto command;
case 'D':
if (v_fp == NULL) {
fputs("Must <G>enerate or <L>oad first!\n",console);
goto command;
}
if (disp_mand() == NULL)
wait_close();
goto command;
case 'G':
if (gen_mand() == NULL)
wait_close();
goto command;
case 'A':
if (v_fp == NULL) {
fputs("Must <G>enerate or <L>oad first!\n",console);
goto command;
}
if (!(color_mode & 1)) {
fputs("Cannot be in hold and modify!\n",console);
goto command;
}
anal_mand();
goto command;
case ';':
goto command; /* Lattice will complain about this line! */
case '<':
if (redir_fp) {
fclose(redir_fp);
redir_fp = NULL;
}
redir_fp = fopen(stpblk(cmd+1),"r");
if (redir_fp == NULL)
fprintf(console,"Can't open file %s!\n",stpblk(cmd+1));
goto command;
case 'Q':
abort("Bye!");
case '?':
case 'H':
default:
AvailableCommands();
goto command;
}
}
}
ill_cmd()
{
fputs("Unknown command!\n",console);
}
float cnvf(i)
int i;
{
union kludge n;
n.i = i;
n.i = SPTieee(n.i);
return (n.f);
}
/*-----------------------------------------*/
/* Non-intelligent virtual memory handling */
v_pos_line(l)
int l;
{
if (l < v_starty) {
if (modified)
v_flush();
v_starty = 0;
fseek(v_fp,v_offset,0);
if (want_read)
fread(v_mand_store,max_x*sizeof(UWORD),max_mem,v_fp);
}
if (l-v_starty > max_mem-1) {
if (modified)
v_flush();
v_starty += max_mem;
fseek(v_fp,(long)(v_starty*max_x*sizeof(UWORD)+v_offset),0);
if (want_read)
fread(v_mand_store,max_x*sizeof(UWORD),max_mem,v_fp);
}
}
v_flush()
{
fseek(v_fp,(long)(v_starty*max_x*sizeof(UWORD)+v_offset),0);
fwrite(v_mand_store,max_x*sizeof(UWORD),max_mem,v_fp);
modified = FALSE;
}
abort(s)
char *s;
{
if (cur_resource & F_MATHTRANS)
CloseLibrary(MathTransBase);
if (cur_resource & F_MATH)
CloseLibrary(MathBase);
if (cur_resource & F_CONSOLE)
fclose(console);
if (v_fp)
fclose(v_fp);
if (redir_fp)
fclose(redir_fp);
if (cur_resource & F_GRAPHICS)
CloseLibrary(GfxBase);
if (cur_resource & F_INTUITION)
CloseLibrary(IntuitionBase);
if (cur_resource & F_ICON)
CloseLibrary(IconBase);
if (cur_resource & F_SETSTORE)
free(v_mand_store);
if (cur_resource & F_COLORTAB)
free(color_table);
CloseDisplay();
puts(s);
exit(0);
}
SetPresets(preset)
{
/* these are some common defaults, preset here to avoid repetitive
* assignments in the case statements below. Feel free to override
* any of these defaults with your own in the case statements below
*/
max_x = 320;
max_y = 200;
max_count = 29;
color_inc = 1;
color_set = 1;
color_mode = 1;
color_div = 1;
func_num = 0;
switch (preset)
{
case 0:
fputs("Mandelbrot's Set\n", console);
start_r.f = -2.85;
end_r.f = 2.85;
start_i.f = -2.05;
end_i.f = 2.05;
max_count = 15;
color_offset = 91;
break;
case 1:
fputs("Unbounded Rhythms\n", console);
start_r.f = 1.703418;
end_r.f = 2.016930;
start_i.f = -0.169450;
end_i.f = 0.169450;
color_offset = 31;
break;
case 2:
fputs("RJ's Gangliactic\n", console);
start_r.f = 1.937821;
end_r.f = 1.945044;
start_i.f = -0.00259;
end_i.f = 0.00259;
color_offset = 31;
break;
case 3:
fputs("Mandelbrot Recursion\n", console);
start_r.f = -0.294473;
end_r.f = -0.288444;
start_i.f = 0.012677;
end_i.f = 0.014839;
color_offset = 26;
max_count = 90;
break;
case 4:
fputs("Crackle\n", console);
start_r.f = 0.225;
end_r.f = 0.275;
start_i.f = 0.8425;
end_i.f = 0.8525;
color_offset = 61;
max_count = 31;
break;
case 5:
fputs("Connections\n", console);
start_r.f = 0.115206;
end_r.f = 0.168190;
start_i.f = -1.036059;
end_i.f = -0.985386;
color_offset = 76;
color_inc = 2;
break;
default:
fputs("SORRY: there's not that many presets!\n",console);
break;
}
ZoomCenterX = max_x >> 1;
ZoomCenterY = max_y >> 1;
ZoomBoxSizeX = max_x;
ZoomBoxSizeY = max_y;
start_r.i = SPFieee(start_r.i);
end_r.i = SPFieee(end_r.i);
start_i.i = SPFieee(start_i.i);
end_i.i = SPFieee(end_i.i);
max_mem = max_mem_y * MAXX / max_x;
init_colors();
}