/* Resident odds and ends that don't fit anywhere else. */ #include #include #include #include #include #include #include #ifdef __TURBOC__ #include #endif #include "fractint.h" #include "fractype.h" #include "helpdefs.h" /* routines in this module */ extern void restore_active_ovly(void ); extern void findpath(char *filename,char *fullpathname); extern void notdiskmsg(void ); extern int cvtcentermag(double *Xctr,double *Yctr,double *Magnification); extern void updatesavename(char *name); extern int check_writefile(char *name,char *ext); extern int check_key(void ); extern int timer(int timertype,int(*subrtn)(),...); extern void showtrig(char *buf); extern int set_trig_array(int k,char *name); extern void set_trig_pointers(int which); extern int tab_display(void ); extern int endswithslash(char *fl); extern int ifsload(void); extern int find_file_item(char *filename,char *itemname,FILE **infile); extern int file_gets(char *buf,int maxlen,FILE *infile); extern void roundfloatd(double *); static void trigdetails(char *); int active_ovly = -1; long timer_start,timer_interval; /* timer(...) start & total */ extern char IFSFileName[80]; extern char IFSName[40]; extern float far *ifs_defn; extern int ifs_changed; extern int ifs_type; extern char temp[], temp1[256]; /* temporary strings */ extern int active_ovly; extern int xdots, ydots; extern int dotmode; extern int show_orbit; extern int timerflag; extern int debugflag; extern int maxit; extern int fractype; extern int got_status,curpass,totpasses,currow,curcol; extern int fileydots; extern int xxstart,xxstop,yystart,yystop; extern int display3d; extern char overwrite; /* call next when returning from resident routine and unsure whether caller is an overlay which has been displaced */ void restore_active_ovly() { switch (active_ovly) { case OVLY_MISCOVL: miscovl_overlay(); break; case OVLY_CMDFILES: cmdfiles_overlay(); break; case OVLY_HELP: help_overlay(); break; case OVLY_PROMPTS: prompts_overlay(); break; case OVLY_LOADFILE: loadfile_overlay(); break; case OVLY_ROTATE: rotate_overlay(); break; case OVLY_PRINTER: printer_overlay(); break; case OVLY_LINE3D: line3d_overlay(); break; case OVLY_ENCODER: encoder_overlay(); break; case OVLY_CALCFRAC: calcfrac_overlay(); break; case OVLY_INTRO: intro_overlay(); break; } } void findpath(char *filename, char *fullpathname) /* return full pathnames */ { if (filename[0] == '\\' || (filename[0] && filename[1] == ':')) { strcpy(fullpathname,filename); return; } fullpathname[0] = 0; /* indicate none found */ #ifdef __TURBOC__ /* look for the file */ strcpy(fullpathname,searchpath(filename)); #else _searchenv(filename,"PATH",fullpathname); #endif if (fullpathname[0] != 0) /* found it! */ if (strncmp(&fullpathname[2],"\\\\",2) == 0) /* stupid klooge! */ strcpy(&fullpathname[3],filename); } void notdiskmsg() { static char far sorrymsg[]={"\ I'm sorry, but because of its random-screen-access algorithms, this\n\ type cannot be created using a real-disk based 'video' mode."}; stopmsg(1,sorrymsg); } /* convert corners to center/mag */ int cvtcentermag(double *Xctr, double *Yctr, double *Magnification) { extern double xxmax,xxmin,yymax,yymin,xx3rd,yy3rd; double Width, Height, Radius, Ratio; Width = xxmax - xxmin; Height = yymax - yymin; Ratio = Height / Width; if(xx3rd != xxmin || yy3rd != yymin || Width < 0 || (Width > 1e-8 && (Ratio <= 0.74 || Ratio >= 0.76)) || Ratio < 0.66 || Ratio > 0.84) return(0); /* calculate center and magnification */ Radius = Height / 2.0; *Xctr = xxmin + (Width / 2.0); *Yctr = yymin + Radius; *Magnification = 1.0 / Radius; return(1); } void updatesavename(char *filename) /* go to the next file name */ { char *save, *hold; char name[80],suffix[80]; char *dotptr; strcpy(name,filename); suffix[0] = 0; if ((dotptr = strrchr(name,'.')) != NULL && dotptr > strrchr(name,'\\')) { strcpy(suffix,dotptr); *dotptr = 0; } hold = name + strlen(name) - 1; /* start at the end */ while(hold >= name && (*hold == ' ' || isdigit(*hold))) /* skip backwards */ hold--; hold++; /* recover first digit */ while (*hold == '0') /* skip leading zeros */ hold++; save = hold; while (*save) { /* check for all nines */ if (*save != '9') break; save++; } if (!*save) /* if the whole thing is nines then back */ save = hold - 1; /* up one place. Note that this will eat */ /* your last letter if you go to far. */ else save = hold; sprintf(save,"%d",atoi(hold)+1); /* increment the number */ strcpy(filename,name); strcat(filename,suffix); } int check_writefile(char *name,char *ext) { /* after v16 release, change encoder.c to also use this routine */ char openfile[80]; char opentype[20]; int i; nextname: strcpy(openfile,name); strcpy(opentype,ext); for (i = 0; i < strlen(openfile); i++) if (openfile[i] == '.') { strcpy(opentype,&openfile[i]); openfile[i] = 0; } strcat(openfile,opentype); if (access(openfile,0) != 0) /* file doesn't exist */ return 0; /* file already exists */ if (overwrite == 0) { updatesavename(name); goto nextname; } return 1; } int check_key() { int key; if((key = keypressed()) != 0) { if(key != 'o' && key != 'O') return(-1); getakey(); if (dotmode != 11) show_orbit = 1 - show_orbit; } return(0); } /* timer function: timer(0,(*fractal)()) fractal engine timer(1,NULL,int width) decoder timer(2) encoder */ int timer(int timertype,int(*subrtn)(),...) { va_list arg_marker; /* variable arg list */ char *timestring; time_t ltime; FILE *fp; int out; int i; int do_bench; va_start(arg_marker,subrtn); do_bench = timerflag; /* record time? */ if (timertype == 2) /* encoder, record time only if debug=200 */ do_bench = (debugflag == 200); if(do_bench) fp=fopen("bench","a"); timer_start = clock(); switch(timertype) { case 0: out = (*subrtn)(); break; case 1: i = va_arg(arg_marker,int); out = decoder(i); /* not indirect, safer with overlays */ break; case 2: out = encoder(); /* not indirect, safer with overlays */ break; } /* next assumes CLK_TCK is 10^n, n>=2 */ timer_interval = (clock() - timer_start) / (CLK_TCK/100); if(do_bench) { time(<ime); timestring = ctime(<ime); timestring[24] = 0; /*clobber newline in time string */ switch(timertype) { case 1: fprintf(fp,"decode "); break; case 2: fprintf(fp,"encode "); break; } fprintf(fp,"%s type=%s resolution = %dx%d maxiter=%d", timestring, curfractalspecific->name, xdots, ydots, maxit); fprintf(fp," time= %ld.%02ld secs\n",timer_interval/100,timer_interval%100); if(fp != NULL) fclose(fp); } return(out); } extern void lStkSin(void) ,dStkSin(void) ,mStkSin(void) ; extern void lStkCos(void) ,dStkCos(void) ,mStkCos(void) ; extern void lStkSinh(void),dStkSinh(void),mStkSinh(void); extern void lStkCosh(void),dStkCosh(void),mStkCosh(void); extern void lStkExp(void) ,dStkExp(void) ,mStkExp(void) ; extern void lStkLog(void) ,dStkLog(void) ,mStkLog(void) ; extern void lStkSqr(void) ,dStkSqr(void) ,mStkSqr(void) ; extern void lStkRecip(void) ,dStkRecip(void) ,mStkRecip(void) ; extern void StkIdent(void); extern void lStkTan(void) ,dStkTan(void) ,mStkTan(void) ; extern void lStkTanh(void),dStkTanh(void),mStkTanh(void); extern void lStkCoTan(void),dStkCoTan(void),mStkCoTan(void); extern void lStkCoTanh(void),dStkCoTanh(void),mStkCoTanh(void); extern void lStkCosXX(void) ,dStkCosXX(void) ,mStkCosXX(void) ; unsigned char trigndx[] = {SIN,SQR,SINH,COSH}; void (*ltrig0)() = lStkSin; void (*ltrig1)() = lStkSqr; void (*ltrig2)() = lStkSinh; void (*ltrig3)() = lStkCosh; void (*dtrig0)() = dStkSin; void (*dtrig1)() = dStkSqr; void (*dtrig2)() = dStkSinh; void (*dtrig3)() = dStkCosh; void (*mtrig0)() = mStkSin; void (*mtrig1)() = mStkSqr; void (*mtrig2)() = mStkSinh; void (*mtrig3)() = mStkCosh; struct trig_funct_lst trigfn[] = /* changing the order of these alters meaning of *.fra file */ /* maximum 6 characters in function names or recheck all related code */ { {"sin", lStkSin, dStkSin, mStkSin }, {"cosxx", lStkCosXX, dStkCosXX, mStkCosXX }, {"sinh", lStkSinh, dStkSinh, mStkSinh }, {"cosh", lStkCosh, dStkCosh, mStkCosh }, {"exp", lStkExp, dStkExp, mStkExp }, {"log", lStkLog, dStkLog, mStkLog }, {"sqr", lStkSqr, dStkSqr, mStkSqr }, {"recip", lStkRecip, dStkRecip, mStkRecip }, /* from recip on new in v16 */ {"ident", StkIdent, StkIdent, StkIdent }, {"cos", lStkCos, dStkCos, mStkCos }, {"tan", lStkTan, dStkTan, mStkTan }, {"tanh", lStkTanh, dStkTanh, mStkTanh }, {"cotan", lStkCoTan, dStkCoTan, mStkCoTan }, {"cotanh",lStkCoTanh,dStkCoTanh,mStkCoTanh}, }; int numtrigfn = sizeof(trigfn)/sizeof(struct trig_funct_lst); void showtrig(char *buf) /* return display form of active trig functions */ { char tmpbuf[30]; *buf = 0; /* null string if none */ trigdetails(tmpbuf); if (tmpbuf[0]) sprintf(buf," function=%s",tmpbuf); } static void trigdetails(char *buf) { extern char maxfn; int i, numfn; char tmpbuf[20]; numfn = (curfractalspecific->flags >> 6) & 7; if(curfractalspecific == &fractalspecific[FORMULA] || curfractalspecific == &fractalspecific[FFORMULA] ) numfn = maxfn; *buf = 0; /* null string if none */ if (numfn>0) { strcpy(buf,trigfn[trigndx[0]].name); i = 0; while(++i < numfn) { sprintf(tmpbuf,"/%s",trigfn[trigndx[i]].name); strcat(buf,tmpbuf); } } } /* set array of trig function indices according to "function=" command */ int set_trig_array(int k, char *name) { char trigname[10]; int i; char *slash; strncpy(trigname,name,6); trigname[6] = 0; /* safety first */ if ((slash = strchr(trigname,'/'))) *slash = 0; strlwr(trigname); for(i=0;i=2 */ calctime += (clock() - timer_start) / (CLK_TCK/100); stackscreen(); helptitle(); setattr(1,0,C_GENERAL_MED,24*80); /* init rest to background */ row = 2; putstring(row,2,C_GENERAL_MED,"Fractal type:"); if (display3d > 0) putstring(row,16,C_GENERAL_HI,"3D Transform"); else { putstring(row,16,C_GENERAL_HI, curfractalspecific->name[0] == '*' ? &curfractalspecific->name[1] : curfractalspecific->name); i = 0; if (fractype == FORMULA || fractype == FFORMULA) { putstring(row+1,16,C_GENERAL_HI,FormName); i = strlen(FormName)+1; } trigdetails(msg); putstring(row+1,16+i,C_GENERAL_HI,msg); if (fractype == LSYSTEM) putstring(row+1,16,C_GENERAL_HI,LName); if (fractype == IFS || fractype == IFS3D) putstring(row+1,16,C_GENERAL_HI,IFSName); } switch (calc_status) { case 0: msgptr = "Parms chgd since generated"; break; case 1: msgptr = "Still being generated"; break; case 2: msgptr = "Interrupted, resumable"; break; case 3: msgptr = "Interrupted, non-resumable"; break; case 4: msgptr = "Image completed"; break; default: msgptr = ""; } putstring(row,45,C_GENERAL_HI,msgptr); if (helpmode == HELPCYCLING) putstring(row+1,45,C_GENERAL_HI,"You are in color-cycling mode"); row += 2; i = j = 0; if (display3d > 0) { if (usr_floatflag) j = 1; } else if (floatflag) j = (usr_floatflag) ? 1 : 2; if (j) { putstring(row,45,C_GENERAL_HI,"Floating-point"); putstring(-1,-1,C_GENERAL_HI,(j == 1) ? " flag is activated" : " in use (required)"); i = 1; } if (calc_status == 1 || calc_status == 2) if (curfractalspecific->flags&INFCALC) { putstring(row,2,C_GENERAL_HI,"Note: this type runs forever."); i = 1; } row += i; if (calc_status == 1 || calc_status == 2) if (curfractalspecific->flags&NORESUME) putstring(row++,2,C_GENERAL_HI,"Note: can't resume this type after interrupts other than and "); ++row; if (got_status >= 0 && (calc_status == 1 || calc_status == 2)) { switch (got_status) { case 0: sprintf(msg,"%d Pass Mode",totpasses); putstring(row,2,C_GENERAL_HI,msg); break; case 1: putstring(row,2,C_GENERAL_HI,"Solid Guessing"); break; case 2: putstring(row,2,C_GENERAL_HI,"Boundary Tracing"); break; case 3: sprintf(msg,"Processing row %d (of %d) of input image",currow,fileydots); putstring(row,2,C_GENERAL_HI,msg); break; } ++row; if (got_status != 3) { sprintf(msg,"Working on block (y,x) [%d,%d]...[%d,%d], ", yystart,xxstart,yystop,xxstop); putstring(row,2,C_GENERAL_MED,msg); if (got_status == 2) { /* btm */ putstring(-1,-1,C_GENERAL_MED,"at "); sprintf(msg,"[%d,%d]",currow,curcol); putstring(-1,-1,C_GENERAL_HI,msg); } else { if (totpasses > 1) { putstring(-1,-1,C_GENERAL_MED,"pass "); sprintf(msg,"%d",curpass); putstring(-1,-1,C_GENERAL_HI,msg); putstring(-1,-1,C_GENERAL_MED," of "); sprintf(msg,"%d",totpasses); putstring(-1,-1,C_GENERAL_HI,msg); putstring(-1,-1,C_GENERAL_MED,", "); } putstring(-1,-1,C_GENERAL_MED,"at row "); sprintf(msg,"%d",currow); putstring(-1,-1,C_GENERAL_HI,msg); } ++row; } } putstring(row,2,C_GENERAL_MED,"Calculation time:"); sprintf(msg,"%3ld:%02ld:%02ld.%02ld", calctime/360000, (calctime%360000)/6000, (calctime%6000)/100, calctime%100); putstring(-1,-1,C_GENERAL_HI,msg); row += 2; if (videoentry.xdots) { sprintf(msg,"Video: %dx%dx%d %s %s", videoentry.xdots, videoentry.ydots, videoentry.colors, videoentry.name, videoentry.comment); putstring(row,2,C_GENERAL_MED,msg); } ++row; putstring(row,2,C_GENERAL_MED,"Corners: X Y"); putstring(++row,3,C_GENERAL_MED,"top-left"); sprintf(msg,"%20.16f %20.16f",xxmin,yymax); putstring(-1,17,C_GENERAL_HI,msg); putstring(++row,3,C_GENERAL_MED,"bottom-right"); sprintf(msg,"%20.16f %20.16f",xxmax,yymin); putstring(-1,17,C_GENERAL_HI,msg); adjust_corner(); /* make bottom left exact if very near exact */ if (cvtcentermag(&Xctr, &Yctr, &Magnification)) { putstring(row+=2,2,C_GENERAL_MED,"Center: "); sprintf(msg,"%20.16f %20.16f",Xctr,Yctr); putstring(-1,-1,C_GENERAL_HI,msg); putstring(-1,-1,C_GENERAL_MED," Mag: "); if (Magnification < 1e6) sprintf(msg,"%20.14f",Magnification); else if (Magnification < 1e12) sprintf(msg,"%20.8f",Magnification); else sprintf(msg,"%20.2f",Magnification); putstring(-1,-1,C_GENERAL_HI,msg); } else if (xxmin != xx3rd || yymin != yy3rd) { putstring(++row,3,C_GENERAL_MED,"bottom-left"); sprintf(msg,"%20.16f %20.16f",xx3rd,yy3rd); putstring(-1,17,C_GENERAL_HI,msg); } putstring(row+=2,2,C_GENERAL_MED,"Params,"); for (i = 0; i < 4; i++) { sprintf(msg,"%3d: ",i+1); putstring(-1,-1,C_GENERAL_MED,msg); sprintf(msg,"%12.9f",param[i]); putstring(-1,-1,C_GENERAL_HI,msg); } putstring(row+=2,2,C_GENERAL_MED,"Iteration maximum: "); sprintf(msg,"%d",maxit); putstring(-1,-1,C_GENERAL_HI,msg); putstring(-1,-1,C_GENERAL_MED," Effective bailout: "); sprintf(msg,"%f",rqlim); putstring(-1,-1,C_GENERAL_HI,msg); if (fractype == PLASMA) { putstring(++row,2,C_GENERAL_MED,"Current 'rseed': "); sprintf(msg,"%d",rseed); putstring(-1,-1,C_GENERAL_HI,msg); } if(invert) { extern double f_radius,f_xcenter,f_ycenter; putstring(++row,2,C_GENERAL_MED,"Inversion radius: "); sprintf(msg,"%12.9f",f_radius); putstring(-1,-1,C_GENERAL_HI,msg); putstring(-1,-1,C_GENERAL_MED," xcenter: "); sprintf(msg,"%12.9f",f_xcenter); putstring(-1,-1,C_GENERAL_HI,msg); putstring(-1,-1,C_GENERAL_MED," ycenter: "); sprintf(msg,"%12.9f",f_ycenter); putstring(-1,-1,C_GENERAL_HI,msg); } if ((row += 2) < 23) ++row; putstringcenter(row,0,80,C_GENERAL_LO,"...Press any key to continue..."); movecursor(25,80); getakeynohelp(); unstackscreen(); timer_start = clock(); /* tab display was "time out" */ return(0); } int endswithslash(char *fl) { int len; len = strlen(fl); if(len) if(fl[--len]=='\\') return(1); return(0); } /* --------------------------------------------------------------------- */ int ifsload() /* read in IFS parameters */ { int i; FILE *ifsfile; char buf[201]; char *bufptr; extern float dstack[]; /* shared temp */ int ret,rowsize; if (ifs_defn) { /* release prior parms */ farmemfree((char far *)ifs_defn); ifs_defn = NULL; } ifs_changed = ifs_type = 0; rowsize = IFSPARM; if (find_file_item(IFSFileName,IFSName,&ifsfile) < 0) return(-1); fgets(buf,200,ifsfile); strlwr(buf); bufptr = &buf[0]; while (*bufptr) { if (strncmp(bufptr,"(3d)",4) == 0) { ifs_type = 1; rowsize = IFS3DPARM; } ++bufptr; } for (i = 0; i < (NUMIFS+1)*IFS3DPARM; ++i) dstack[i] = 0.0; i = ret = 0; while (fscanf(ifsfile," %f ",&dstack[i])) { if (++i >= NUMIFS*rowsize) { static char far msg[]={"IFS definition has too many lines"}; stopmsg(0,msg); ret = -1; break; } } if ((i % rowsize) != 0 || getc(ifsfile) != '}') { static char far msg[]={"invalid IFS definition"}; stopmsg(0,msg); ret = -1; } if (i == 0 && ret == 0) { static char far msg[]={"Empty IFS definition"}; stopmsg(0,msg); ret = -1; } fclose(ifsfile); if (ret == 0) if ((ifs_defn = (float far *)farmemalloc( (long)((NUMIFS+1)*IFS3DPARM*sizeof(float)))) == NULL) { static char far msg[]={"Insufficient memory for IFS"}; stopmsg(0,msg); ret = -1; } else for (i = 0; i < (NUMIFS+1)*IFS3DPARM; ++i) ifs_defn[i] = dstack[i]; return(ret); } int find_file_item(char *filename,char *itemname,FILE **infile) { char tmpname[41]; long notepoint; char buf[201]; int c; if ((*infile = fopen(filename,"rt")) == NULL) { sprintf(buf,"Can't open %s",filename); stopmsg(0,buf); return(-1); } while (1) { while ((c = getc(*infile)) == ' ' || c == '\t' || c == '\n') { } if (c == EOF) break; if (c == ';') { while ((c = fgetc(*infile)) != '\n' && c != EOF) { } if (c == EOF) break; continue; } notepoint = ftell(*infile) - 1; ungetc(c,*infile); if (fscanf(*infile," %40[^ \n\t({]",tmpname) == EOF) break; while ((c = getc(*infile)) != EOF && c != '{' && c != '\n') { } if (c == EOF) break; if (c == '{') { if (stricmp(tmpname,itemname) == 0) { fseek(*infile,notepoint,SEEK_SET); return(0); } while ((c = getc(*infile)) != '}' && c != EOF) { } if (c == EOF) break; } } fclose(*infile); sprintf(buf,"'%s' definition not found",itemname); stopmsg(0,buf); return(-1); } int file_gets(char *buf,int maxlen,FILE *infile) { int len,c; /* similar to 'fgets', but file may be in either text or binary mode */ /* returns -1 at eof, length of string otherwise */ if (feof(infile)) return -1; len = 0; while (len < maxlen) { if ((c = getc(infile)) == EOF || c == '\x1a') { if (len) break; return -1; } if (c == '\n') break; /* linefeed is end of line */ if (c != '\r') buf[len++] = c; /* ignore c/r */ } buf[len] = 0; return len; } int matherr( struct exception *except ) { extern int debugflag; if(debugflag) { static int ct = 0; static FILE *fp=NULL; if(fp==NULL) fp = fopen("matherr","w"); if(ct++ < 100) { fprintf(fp,"err: %d\nname: %s\narg: %le\n", except->type, except->name, except->arg1); fflush(fp); } } if( except->type == DOMAIN ) { char buf[40]; sprintf(buf,"%le",except->arg1); /* This test may be unnecessary - from my experiments if the argument is too large or small the error is TLOSS not DOMAIN */ if(strstr(buf,"IN")||strstr(buf,"NAN")) /* trashed arg? */ /* "IND" with MSC, "INF" with BC++ */ { if( strcmp( except->name, "sin" ) == 0 ) { except->retval = 0.0; return(1); } else if( strcmp( except->name, "cos" ) == 0 ) { except->retval = 1.0; return(1); } } } return (0); } void roundfloatd(double *x) /* make double converted from float look ok */ { char buf[30]; sprintf(buf,"%-10.7g",*x); *x = atof(buf); }