Home » Archimedes archive » Archimedes World » AW-1996-04-Disc2.adf » !AcornAns_AcornAns » !Scape2/c/main
!Scape2/c/main
This website contains an archive of files for the Acorn Electron, BBC Micro, Acorn Archimedes, Commodore 16 and Commodore 64 computers, which Dominic Ford has rescued from his private collection of floppy disks and cassettes.
Some of these files were originally commercial releases in the 1980s and 1990s, but they are now widely available online. I assume that copyright over them is no longer being asserted. If you own the copyright and would like files to be removed, please contact me.
| Tape/disk: | Home » Archimedes archive » Archimedes World » AW-1996-04-Disc2.adf » !AcornAns_AcornAns |
| Filename: | !Scape2/c/main |
| Read OK: | ✔ |
| File size: | 65DD bytes |
| Load address: | 0000 |
| Exec address: | 0000 |
File contents
#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>
#include <string.h>
#include <time.h>
#include <math.h>
#include "wimpt.h"
#include "werr.h"
#include "bbc.h"
#include "os.h"
#include "flex.h"
#include "sprite.h"
#include "akbd.h"
#include "swis.h"
#define one (1<<16)
#define dimsn(f) ((int)(f)/5.0+2.0)
typedef enum {multi, nonmulti} monitor;
monitor monitortype;
int linesep;
typedef enum {
left=136,
right,
down,
up
} scrolldir;
typedef struct {
int re;
int im;
} complex;
int N;
int limit;
void *data = 0;
int tilt=one/3; /*in range 0 to one*/
int range;
int sealevel;
int powscale;
int linscale;
int orient, scrollx, scrolly;
int screen=0;
char *imname = "scape";
sprite_area *sa32=0;
int im32_yl;
int im32_height;
int dz;
int lx,ly,lz;
int vx,vy,vz;
int hx,hy,hz;
sprite_area *ditharea=0;
char lut[32768];
unsigned int pal256def[256];
unsigned int earth[256];
BOOL save32 = TRUE;
BOOL save8 = TRUE;
typedef struct {
signed char er;
signed char eg;
signed char eb;
signed char byte4; /*reminder that struct will be word sized*/
} compact_dither_error;
BOOL gen_256withnopalette(flex_ptr sa, char *name);
extern void g256wnp_internal(int x, int y, char *dp, char *sp, compact_dither_error *de, char *lut,
unsigned int *pal256def);
extern int rbbcinc(int r, int k);
extern int div_frac16(int number, int divisor);
extern int mul_frac16(int x, int a);
extern int mul_frac16c(int x, int a);
extern int sqrt_frac16(unsigned int x);
extern int gauss16(void);
extern void sgauss16(int seed);
extern int cos16(int a);
extern int sin16(int a);
extern int exp16(int a);
extern int ln16(int a);
extern int pow16(int a, int b);
/************************************** fourier transform functions ****************************************/
BOOL invfft_2d(complex *data, int limit)
{
int i, j;
int n = 1<<limit;
complex t1, *d1, *d1e, *d2;
int k, l;
int loop;
int rotangle, angle;
int cos, sin;
if (limit<3 || limit>32) return FALSE;
for (i=j=0; j<n; i=rbbcinc(i, limit), j++) if (i>j) {
for (d1=data+(i<<limit), d2=data+(j<<limit), d1e=d1+n; d1<d1e;) {
t1 = *d1;
*d1++ = *d2;
*d2++ = t1;
}
}
rotangle = -2*one;
for (k=loop=1; loop<=limit; loop++, k<<=1, rotangle/=2) {
for (angle=0, sin=0, cos=one, i=1; i<=k; i++, angle+=rotangle, sin=sin16(angle), cos=cos16(angle)) {
for (j=i-1; j<n; j+=k<<1) {
l = j+k;
for (d1=data+(l<<limit), d2=data+(j<<limit), d1e=d1+n; d1<d1e;) {
t1.re = (mul_frac16c(d1->re, cos)-mul_frac16c(d1->im, sin));
t1.im = (mul_frac16c(d1->re, sin)+mul_frac16c(d1->im, cos));
d1->re = d2->re-t1.re;
d1++->im = d2->im-t1.im;
d2->re += t1.re;
d2++->im += t1.im;
}
}
}
}
for (i=j=0; j<n; i=rbbcinc(i, limit), j++) if (i>j) {
for (d1=data+i, d2=data+j, d1e=d1+n*n; d1<d1e; d1+=n, d2+=n) {
t1 = *d1;
*d1 = *d2;
*d2 = t1;
}
}
rotangle = -2*one;
for (k=loop=1; loop<=limit; loop++, k<<=1, rotangle/=2) {
for (angle=0, sin=0, cos=one, i=1; i<=k; i++, angle+=rotangle, sin=sin16(angle), cos=cos16(angle)) {
for (j=i-1; j<n; j+=k<<1) {
l = j+k;
for (d1=data+l, d2=data+j, d1e=d1+n*n; d1<d1e; d1+=n, d2+=n) {
t1.re = (mul_frac16c(d1->re, cos)-mul_frac16c(d1->im, sin));
t1.im = (mul_frac16c(d1->re, sin)+mul_frac16c(d1->im, cos));
d1->re = d2->re-t1.re;
d1->im = d2->im-t1.im;
d2->re += t1.re;
d2->im += t1.im;
}
}
}
}
return TRUE;
}
/*********************************** data generation/manipulation functions *******************************/
/*The following assigns values to 2d data array with suitable properties.*/
/*It does however do so by scanning the array in a non-obvious order.
This has been chosen so that if n is increased by doubling, then those random values used for the
last grid will be now written to the low frequency quadrant of the new grid (if same seed used),
suitably re filtered. This should have effect of giving a degree of continuity to the change in n,
since former frequencies are carried over to the low frequencies of the higher res grid.
*/
BOOL assign_ft(complex *data, int limit, double dim, int seed)
{
int n=1<<limit;
int i, ic, j, i2, j2, k, l;
complex *d1, *d2, *d;
int phase, rad;
int im, re;
int mbpobt;
srand(-seed);
sgauss16(-seed);
mbpobt = (int)(-(4-dim)*one);
data->re = data->im = 0;
for (i=1; i<=n/2; i++) {
j2 = i2 = i*i;
k = n-i;
d1 = data+i*n;
d2 = data+k*n;
for (j=0; j<=i; j++) {
l = j==0 ? 0 : n-j;
phase = rand()>>13;
rad = mul_frac16(pow16(sqrt_frac16((i2+j*j)*(one/4))*2, mbpobt), gauss16());
re = mul_frac16c(rad, cos16(phase));
im = mul_frac16c(rad, sin16(phase));
d = d1+j;
d->re = re;
d->im = im;
d = d2+l;
d->re = re;
d->im = -im;
}
j = i;
l = n-j;
d1 = data+j;
d2 = data+l;
for (ic=i-1; ic>=0; ic--) {
k = ic==0 ? 0 : n-ic;
phase = rand()>>13;
rad = mul_frac16(pow16(sqrt_frac16((ic*ic+j2)*(one/4))*2, mbpobt), gauss16());
re = mul_frac16c(rad, cos16(phase));
im = mul_frac16c(rad, sin16(phase));
d = d1+ic*n;
d->re = re;
d->im = im;
d = d2+k*n;
d->re = re;
d->im = -im;
}
if (i<n/2) {
ic = n-i;
k = i;
d1 = data+ic*n;
d2 = data+k*n;
for (j=1; j<=i; j++) {
l = n-j;
phase = rand()>>13;
rad = mul_frac16(pow16(sqrt_frac16((i2+j*j)*(one/4))*2, mbpobt), gauss16());
re = mul_frac16c(rad, cos16(phase));
im = mul_frac16c(rad, sin16(phase));
d = d1+j;
d->re = re;
d->im = im;
d = d2+l;
d->re = re;
d->im = -im;
}
j=i;
l = n-j;
d1 = data+j;
d2 = data+l;
for (ic++; ic<n; ic++) {
k = n-ic;
phase = rand()>>13;
rad = mul_frac16(pow16(sqrt_frac16((k*k+j2)*(one/4))*2, mbpobt), gauss16());
re = mul_frac16c(rad, cos16(phase));
im = mul_frac16c(rad, sin16(phase));
d = d1+ic*n;
d->re = re;
d->im = im;
d = d2+k*n;
d->re = re;
d->im = -im;
}
}
}
(data+n*n/2)->im = (data+n/2)->im = (data+n*(n+1)/2)->im = 0;
return TRUE;
}
BOOL normalise_and_compact(void **data, int limit)
{
int n=1<<limit;
int min=1<<30;
int max=-min;
int scale, v;
complex *d1, *d1e;
int *d2;
for (d1=*data, d1e=d1+n*n; d1<d1e; d1++) {
if (d1->re<min) min=d1->re;
if (d1->re>max) max=d1->re;
}
scale = div_frac16(one, max-min);
for (d2=*data, d1=*data, d1e=d1+n*n; d1<d1e; d1++, d2++) {
v = mul_frac16(d1->re-min, scale);
if (v>=one) v=one-1;
if (v<0) v=0;
v = pow16(v, powscale);
if (v<sealevel) v=sealevel;
if (linscale!=one) v=mul_frac16c(v, linscale);
if (v>=one) v=one-1;
if (v<0) v=0;
*d2 = v;
}
/*NB given range on data, could pack each value into 2 chars, but is slower to read & write
half words than whole words, so just reduce from 2 ints down to 1 int.*/
return flex_extend((flex_ptr)data, n*n*sizeof(int));
}
BOOL clip_ft(complex *data, int limit, int rl)
{
int n=1<<limit;
int i, j;
int rh=n-rl;
for (i=0; i<n; i++)
for (j=0; j<n; j++, data++)
if (i>rl && i<rh || j>rl && j<rh) data->re=data->im=0;
return TRUE;
}
BOOL scroll(int *data, int limit, scrolldir key, int offset)
{
int n=1<<limit;
int i, j;
int *p, *p1, *p2;
int buf[512];
int si=offset<<(9-limit);
switch (key) {
case left:
scrollx += orient==0 ? -si : (orient==2 ? si : 0);
scrolly += orient==1 ? -si : (orient==3 ? si : 0);
break;
case right:
scrollx += orient==2 ? -si : (orient==0 ? si : 0);
scrolly += orient==3 ? -si : (orient==1 ? si : 0);
offset=n-offset;
key=left;
break;
case down:
scrollx += orient==3 ? -si : (orient==1 ? si : 0);
scrolly += orient==0 ? -si : (orient==2 ? si : 0);
break;
case up:
scrollx += orient==1 ? -si : (orient==3 ? si : 0);
scrolly += orient==2 ? -si : (orient==0 ? si : 0);
offset=n-offset;
key=down;
break;
}
scrollx &= 511;
scrolly &= 511;
offset &= n-1;
if (key==left) {
for (i=0, p=data; i<n; i++, p+=n) {
for (j=n-offset, p1=p+j, p2=buf; j<n; j++, p1++, p2++) *p2=*p1;
for (j=n-offset-1, p1=p+j, p2=p+j+offset; j>=0; j--, p1--, p2--) *p2=*p1;
for (j=0, p1=buf, p2=p; j<offset; j++, p1++, p2++) *p2=*p1;
}
}
else {
for (j=0, p=data; j<n; j++, p++) {
for (i=n-offset, p1=p+i*n, p2=buf; i<n; i++, p1+=n, p2++) *p2=*p1;
for (i=n-offset-1, p1=p+i*n, p2=p+(i+offset)*n; i>=0; i--, p1-=n, p2-=n) *p2=*p1;
for (i=0, p1=buf, p2=p; i<offset; i++, p1++, p2+=n) *p2=*p1;
}
}
return TRUE;
}
BOOL rotate(int *data, int limit, BOOL cw)
{
int n=1<<limit;
int i, j;
int k = n-1;
int l = n/2;
int *q1, *q2, *q3, *q4;
int t;
if (cw) {
orient++;
for (i=0, q1=data, q2=data+n*k, q3=data+(n+1)*k, q4=data+k; i<l;
i++, q1+=l, q2+=l*n+1, q3-=l, q4-=l*n+1)
for (j=0; j<l; j++, q1++, q2-=n, q3--, q4+=n) {
t=*q1;
*q1=*q4;
*q4=*q3;
*q3=*q2;
*q2=t;
}
}
else {
orient--;
for (i=0, q1=data, q2=data+n*k, q3=data+(n+1)*k, q4=data+k; i<l;
i++, q1+=l, q2+=l*n+1, q3-=l, q4-=l*n+1)
for (j=0; j<l; j++, q1++, q2-=n, q3--, q4+=n) {
t=*q1;
*q1=*q2;
*q2=*q3;
*q3=*q4;
*q4=t;
}
}
orient &= 3;
return TRUE;
}
/*********************************** colour rendering functions ********************************************/
BOOL create_image32(int x, int y)
{
int areasize, spritesize, imagesize;
sprite_header *sph;
if (sa32) flex_free((flex_ptr)&sa32);
imagesize = x*y*4;
spritesize = sizeof(sprite_header)+imagesize;
areasize = sizeof(sprite_area)+spritesize;
if (!flex_alloc((flex_ptr)&sa32, areasize)) return FALSE;
sa32->size = areasize;
sa32->number = 1;
sa32->sproff = 16;
sa32->freeoff = areasize;
sph = (sprite_header *)(sa32+1);
sph->next = spritesize;
memset(sph->name, 0, 12);
strcpy(sph->name, imname);
sph->width = x-1;
sph->height = y-1;
sph->lbit = 0;
sph->rbit = 31;
sph->image = sizeof(sprite_header);
sph->mask = sizeof(sprite_header);
sph->mode = 1 | (90<<1) | (90<<14) | (6<<27);
memset((char *)sph+sph->image, 0, imagesize);
return TRUE;
}
BOOL discard_image32(void)
{
if (sa32) flex_free((flex_ptr)&sa32);
return TRUE;
}
int calc_intens(int *data, int limit, int i, int j)
{
int n=1<<limit;
int nx,ny,nz,r,p;
int z;
int in=0;
z=data[i+j*n];
nx=z-data[((i+1)&(n-1))+j*n];
ny=data[i+((j-1)&(n-1))*n]-z;
r = sqrt_frac16(nx*nx+ny*ny+dz*dz)>>8;
nx = div_frac16(nx, r);
ny = div_frac16(ny, r);
nz = div_frac16(dz, r); /*surface normal*/
r = mul_frac16c(nx,vx)+mul_frac16c(ny,vy)+mul_frac16c(nz,vz);
if (r<0) return in;
in+=(int)(0.05*one);
r = mul_frac16c(nx,lx)+mul_frac16c(ny,ly)+mul_frac16c(nz,lz);
if (r<0) return in;
in += mul_frac16c(r, (int)(0.75*one));
r = mul_frac16c(nx,hx)+mul_frac16c(ny,hy)+mul_frac16c(nz,hz);
p = mul_frac16c(r,r);
p = mul_frac16c(p,p); /*^4, for ^3 use r,p instead of p,p*/
in += mul_frac16c(p, (int)(0.45*one));
return in;
}
BOOL render_image32(int *data, int limit, int yl, int xm, int ym)
{
int n=1<<limit;
int i,j;
int horizon;
int p,p0,p1;
int cos = cos16(tilt);
int sin = sin16(tilt);
int lastintens;
int lastr, lastg, lastb, newr, newg, newb;
BOOL lastvalid;
sprite_header *sph = (sprite_header *)((char *)sa32+sa32->sproff);
int *intens = (int *)((char *)sph+sph->image);
int in;
int deltar, deltag, deltab;
int r;
int incl, bear;
unsigned int colour;
int cr,cg,cb;
dz = 65536/(n-1);
vx = 0;
vy = -cos;
vz = sin;
incl = (int)(0.50*one);
bear = (int)(1.50*one);
r = cos16(incl);
lx = mul_frac16c(cos16(bear), r);
ly = mul_frac16c(sin16(bear),-r);
lz = sin16(incl);
hx = (lx+vx)/2;
hy = (ly+vy)/2;
hz = (lz+vz)/2;
r = sqrt_frac16(mul_frac16c(hx,hx)+mul_frac16c(hy,hy)+mul_frac16c(hz,hz));
hx = div_frac16(hx, r);
hy = div_frac16(hy, r);
hz = div_frac16(hz, r);
for (i=0; i<n; i++) {
p0 = ( (mul_frac16(data[i], cos)+mul_frac16(0<<(16-limit), sin))*3/256 - yl) >> (10-limit);
colour = earth[data[i]/256];
cr = colour & 255;
cg = (colour>>8) & 255;
cb = (colour>>16) & 255;
in = calc_intens(data, limit, i, 0);
cr = (cr*in)/256;
cg = (cg*in)/256;
cb = (cb*in)/256;
if (cr<0) cr=0;
if (cg<0) cg=0;
if (cb<0) cb=0;
if (cr>65535) cr=65535;
if (cg>65535) cg=65535;
if (cb>65535) cb=65535;
lastr = cr;
lastg = cg;
lastb = cb;
lastvalid=TRUE;
if (p0>=0 && p0<ym) intens[i+xm*(ym-1-p0)] = (cr>>8)+(cg&(255<<8))+((cb>>8)<<16);
horizon = p0;
for (j=1; j<n; j++) {
p1 = ( (mul_frac16(data[i+j*n], cos)+mul_frac16(j<<(16-limit), sin))*3/256 - yl) >> (10-limit);
if (p1>horizon) {
colour = earth[data[i+j*n]/256];
cr = colour & 255;
cg = (colour>>8) & 255;
cb = (colour>>16) & 255;
in = calc_intens(data, limit, i, j);
cr = (cr*in)/256;
cg = (cg*in)/256;
cb = (cb*in)/256;
if (cr<0) cr=0;
if (cg<0) cg=0;
if (cb<0) cb=0;
if (cr>65535) cr=65535;
if (cg>65535) cg=65535;
if (cb>65535) cb=65535;
newr = cr;
newg = cg;
newb = cb;
if (p1>=0 && p1<ym) intens[i+xm*(ym-1-p1)] = (cr>>8)+(cg&(255<<8))+((cb>>8)<<16);
p = p1-1;
if (p>horizon) {
if (!lastvalid) {
colour = earth[data[i+(j-1)*n]/256];
lastr = colour & 255;
lastg = (colour>>8) & 255;
lastb = (colour>>16) & 255;
lastintens = calc_intens(data, limit, i, j-1);
lastr = (lastr*lastintens)/256;
lastg = (lastg*lastintens)/256;
lastb = (lastb*lastintens)/256;
if (lastr<0) lastr=0;
if (lastg<0) lastg=0;
if (lastb<0) lastb=0;
if (lastr>65535) lastr=65535;
if (lastg>65535) lastg=65535;
if (lastb>65535) lastb=65535;
}
deltar = (lastr-newr)/(p0-p1);
deltag = (lastg-newg)/(p0-p1);
deltab = (lastb-newb)/(p0-p1);
while (p>horizon) {
cr-=deltar;
cg-=deltag;
cb-=deltab;
if (p>=0 && p<ym) intens[i+xm*(ym-1-p)] = (cr>>8)+(cg&(255<<8))+((cb>>8)<<16);
p-=1;
}
}
horizon = p1;
lastr=newr;
lastg=newg;
lastb=newb;
lastvalid = TRUE;
}
else {
lastvalid=FALSE;
}
p0=p1;
}
}
return TRUE;
}
BOOL dump_image32(int *data, int limit, int yl, int height, char *filename)
{
int n=1<<limit;
BOOL retcode=TRUE;
int x=n;
int y=height;
char buf[256];
if (!create_image32(x, y)) return FALSE;
if (!render_image32(data, limit, yl, x, y)) {
discard_image32();
return FALSE;
}
if (save32) {
if (sprite_area_save(sa32, filename)) retcode=FALSE;
}
if (save8) {
if (ditharea) flex_free((flex_ptr)&ditharea);
if (!flex_alloc((flex_ptr)&ditharea, sizeof(sprite_area))) {
werr(0, "ditharea init failed");
retcode=FALSE;
}
else {
sprite_area_initialise(ditharea, flex_size((flex_ptr)&ditharea));
if (gen_256withnopalette((flex_ptr)&sa32, imname)) {
strcat(strcpy(buf, filename), "d");
if (sprite_area_save(ditharea, buf)) retcode=FALSE;
}
else retcode = FALSE;
flex_free((flex_ptr)&ditharea);
}
}
discard_image32();
return retcode;
}
/******************************** dithering ********************************/
BOOL gen_256withnopalette(flex_ptr sa, char *name)
{
int newsize;
os_error *e;
compact_dither_error *de;
int x, y;
sprite_id si;
sprite_header *sh, *dh;
char *sp, *dp;
unsigned int modeword, type, xres, yres;
int sprsize, extra;
si.tag = sprite_id_name;
si.s.name = name;
if (wimpt_complain(sprite_select_rp((sprite_area *)*sa, &si, (sprite_ptr *)&sh))) return FALSE;
modeword = sh->mode;
type = modeword>>27;
xres = (modeword>> 1)&0x1fff;
yres = (modeword>>14)&0x1fff;
if (type!=6) {
werr(0, "Sprite is not new type 6 (32bpp) - I only understand these!");
return FALSE;
}
x = 1+sh->width;
y = sh->height + 1;
sprsize = sizeof(sprite_header)+y*((x+3)&~3);
extra = sprsize - (ditharea->size - ditharea->freeoff);
if (extra>0) {
newsize = ditharea->size + extra;
if (!flex_extend((flex_ptr)&ditharea, newsize)) {
werr(0, "No room for destination sprite");
return FALSE;
}
ditharea->size = newsize;
/*above flex will likely have invalidated sh, so recalc*/
if (wimpt_complain(sprite_select_rp((sprite_area *)*sa, &si, (sprite_ptr *)&sh))) return FALSE;
}
de = malloc(x*sizeof(compact_dither_error));
if (de==0) {
werr(0, "No room for dither error line buffer");
return FALSE;
}
e=sprite_create_rp(ditharea, name, sprite_nopalette, x, y, xres==2*yres ? 15 : 28, (sprite_ptr *)&dh);
if (e) {
free(de);
wimpt_complain(e);
return FALSE;
}
dp = (char *)dh + dh->image;
sp = (char *)sh + sh->image;
g256wnp_internal(x, y, dp, sp, de, lut, pal256def);
free(de);
return TRUE;
}
/****************************** initialisation ******************************/
/*--- Initialise the program, returning TRUE if it was all OK. ---*/
static BOOL initialise(int argc, char *argv[])
{
os_filestr f;
int i;
char pt[1536], *p;
if (argc>=2) {
i=atoi(argv[1]);
save32 = i&1 ? TRUE : FALSE;
save8 = i&2 ? TRUE : FALSE;
}
f.action=255;
f.name="<Scape2$Dir>.resources.lut";
f.loadaddr=(int)lut;
f.execaddr=0;
if (wimpt_complain(os_file(&f))) return FALSE;
if (f.start!=32768) { /*f.start is r4, here holding file length*/
werr(0, "Internal: resources.lut file corrupted - should be 32768 bytes long");
return FALSE;
}
f.action=255;
f.name="<Scape2$Dir>.resources.pal256def";
f.loadaddr=(int)pal256def;
f.execaddr=0;
if (wimpt_complain(os_file(&f))) return FALSE;
if (f.start!=1024) { /*f.start is r4, here holding file length*/
werr(0, "Internal: resources.pal256def file corrupted - should be 1024 bytes long");
return FALSE;
}
f.action=255;
f.name="<Scape2$Dir>.resources.earth";
f.loadaddr=(int)pt;
f.execaddr=0;
if (wimpt_complain(os_file(&f))) return FALSE;
if (f.start!=1536) { /*f.start is r4, here holding file length*/
werr(0, "Internal: resources.earth file corrupted - should be 1536 bytes long");
return FALSE;
}
for (i=0, p=pt; i<256; i++, p+=6) earth[i]=p[3]+(p[4]<<8)+(p[5]<<16);
/* All went ok */
return TRUE;
}
/***************************************** main functions **************************************************/
BOOL graph(int *data, int limit)
{
int n=1<<limit;
int i, j, k;
int pl, pr;
int vl, vr, v, v2;
int horizon[1026];
int horizon2[1026];
int cos = cos16(tilt);
int sin = sin16(tilt);
int image_yl=1024;
int image_yh=0;
for (i=0; i<1026; i++) {
horizon[i]=0;
horizon2[i]=1024;
}
im32_yl=1024;
bbc_origin(640,0);
for (i=0; i<n; i++) {
for (j=1, pl=-512, vl = data[0+n*i]; j<n; j++, pl=pr, vl=vr) {
for (pr=pl+(1024>>limit), vr=data[j+n*i], k=pl, v=vl; k<pr; k+=2, v+=(vr-vl)>>(9-limit)) {
v2 = (mul_frac16(v, cos)+mul_frac16(i<<(16-limit), sin))*3/256;
if (i==0) {
horizon2[k+512]=v2;
if (v2<im32_yl) im32_yl=v2;
}
if (horizon[k+512]<v2) {
bbc_gcol(0, 4+6*i/n);
bbc_plot(69, k, v2);
horizon[k+512]=v2;
}
else if (horizon2[k+512]>v2) {
bbc_gcol(0, 10+6*i/n);
bbc_plot(69, k, v2);
horizon2[k+512]=v2;
}
}
}
v2 = (mul_frac16(vr, cos)+mul_frac16(i<<(16-limit), sin))*3/256;
if (i==0) horizon2[k+512]=v2;
if (horizon[k+512]<v2) {
bbc_gcol(0, 4+6*i/n);
bbc_plot(69, k, v2);
horizon[k+512]=v2;
}
else if (horizon2[k+512]>v2) {
bbc_gcol(0, 10+6*i/n);
bbc_plot(69, k, v2);
horizon2[k+512]=v2;
}
}
for (i=0; i<1026; i++) {
if (horizon[i] >image_yh) image_yh=horizon[i];
if (horizon2[i]<image_yl) image_yl=horizon2[i];
}
bbc_move(-512, image_yl);
bbc_move(512, image_yh);
bbc_plot(189, -512, 512+(image_yl-image_yh)/2);
im32_yl = im32_yl-2;
im32_height = 3 + ( (image_yh-im32_yl-2) >> (10-limit) );
return TRUE;
}
monitor read_monitor_type(void)
{
int r2;
os_swi3r(6, 161,133,0, 0,0,&r2);
if ((r2/4 & 3) == 1) return multi;
else return nonmulti;
}
int main(int argc, char *argv[])
{
int seed, i;
char key;
double dim=2.2;
double d;
char buf[1024];
int torient, tscrollx, tscrolly;
monitortype = read_monitor_type();
flex_init();
if (!initialise(argc, argv)) return 0;
for (;;) {
bbc_mode(monitortype==multi ? 20 : 12);
printf("Please select which landscape to construct (any integer) ");
scanf("%i", &seed);
do {
printf("Please enter scape resolution (suggest 5...8, though 3...9 permissible) ");
scanf("%i", &limit);
} while (limit<3 || limit>9);
N = 1<<limit;
do {
printf("Please enter scape dimension (suggest 2.15, 2-3 permissible) ");
scanf("%Lf", &dim);
} while (dim<2 || dim>3);
do {
printf("Please enter frequency range (suggest 1, 0-1 permissible) ");
scanf("%Lf", &d);
} while (d<0 || d>1);
range = (int)(0.5+d*(double)(int)(N/2));
if (range<0) range=0;
if (range>N/2) range=N/2;
do {
printf("Please enter sea-level (suggest 0, 0-1 permissible) ");
scanf("%Lf", &d);
} while (d<0 || d>1);
sealevel = (int)(d*one);
do {
printf("Please enter power scaling law (suggest 1, or 0.3-3.0, though >= 0 permissible) ");
scanf("%Lf", &d);
} while (d<0);
powscale = (int)(d*one);
do {
printf("Please enter linear scaling factor (suggest 1, 0-1 permissible) ");
scanf("%Lf", &d);
} while (d<0 || d>1);
linscale = (int)(d*one);
do {
printf("Please enter orientation (suggest 0, 0...3 permissible) ");
scanf("%i", &torient);
} while (torient<0 || torient>3);
do {
printf("Please enter x scroll value (suggest 0, 0...511 permissible) ");
scanf("%i", &tscrollx);
} while (tscrollx<0 || tscrollx>=512);
do {
printf("Please enter y scroll value (suggest 0, 0...511 permissible) ");
scanf("%i", &tscrolly);
} while (tscrolly<0 || tscrolly>=512);
do {
printf("Please enter viewing angle (suggest 0.3, or 0.3-0.7, though 0-1 permissible) ");
scanf("%Lf", &d);
} while (d<0 || d>1);
tilt = (int)(d*one);
if (!flex_alloc((flex_ptr)&data, N*N*sizeof(complex))) werr(1, "Insufficient memory");
/*note, further flex calls may shift data, so be careful with ptrs into data*/
bbc_mode(monitortype==multi ? 20 : 12);
bbc_vduq(28,72,31,79,0);
bbc_palette(1, 16, 51,51,51);
bbc_palette(2, 16, 255,255,255);
bbc_palette(0, 16, 192,192,192);
for (i=0; i<6; i++) {
bbc_palette( 4+i, 16, 32+20*i,136+7*i,32+20*i);
bbc_palette(10+i, 16, 128+8*i,96+12*i,64+16*i);
}
linesep = monitortype==multi ? 20 : 40;
printf("\nPlease wait ...\n");
assign_ft(data, limit, dim, seed);
if (range<N/2) clip_ft(data, limit, range);
invfft_2d(data, limit);
if (!normalise_and_compact(&data, limit)) werr(1, "Can't compact data");
orient=scrollx=scrolly=0;
tscrollx>>=9-limit;
tscrolly>>=9-limit;
if (tscrollx>0) scroll(data, limit, right, tscrollx);
if (tscrollx<0) scroll(data, limit, left, -tscrollx);
if (tscrolly>0) scroll(data, limit, up, tscrolly);
if (tscrolly<0) scroll(data, limit, down, -tscrolly);
switch (torient) {
case 1:
rotate(data, limit, TRUE);
break;
case 2:
rotate(data, limit, TRUE);
rotate(data, limit, TRUE);
break;
case 3:
rotate(data, limit, FALSE);
break;
}
do {
bbc_clg();
graph(data, limit);
bbc_vduq(28,0,31,7,0);
printf(" scape\n%7i\n\n res\n%7i\n\n dim\n%7.2f\n\n range\n%7.2f\n\n sea\n%7.2f\n\n pscale\n%7.2f\n\n lscale\n%7.2f\n\n orient\n%7i\n\n x\n%7i\n\n y\n%7i\n\n tilt\n%7.2f", seed, limit, dim, (int)range/(double)(int)(N/2), (int)sealevel/65536.0, (int)powscale/65536.0, (int)linscale/65536.0, orient, scrollx, scrolly, (int)tilt/65536.0);
bbc_vduq(28,72,31,79,0);
bbc_cls();
printf("\n\n\n\n\n\n\n\nPress a key ...\n\n");
printf("[ '<'\n or\n '>'\n to\nrotate\n scape\n\ncursors\n (+\x8b)\n to\nscroll\n scape\n\n copy\n to\n save\n image\n\n\x8b copy\n to\n save\ncolour\n image\n ]\n");
do {
os_cli("fx 15 1");
bbc_cursor(0);
os_swi2(OS_Byte,4,1);
key = bbc_get();
os_swi2(OS_Byte,4,0);
if (key==135) {
screen+=1;
sprintf(buf, "<Scape2$Dir>.^.Image%02i", screen);
if (akbd_pollsh()) {
dump_image32(data, limit, im32_yl, im32_height, buf);
}
else {
wimpt_complain(sprite_screensave(buf, TRUE));
}
}
} while (key==135);
if (key==44 || key==46) rotate(data, limit, key==46);
if (key==136 || key==137 || key==138 || key==139) {
i=limit-2;
if (akbd_pollsh()) i-=2;
if (i<0) i=0;
scroll(data, limit, key, 1<<i);
}
} while (key==44 || key==46 || key==136 || key==137 || key==138 || key==139);
flex_free((flex_ptr)&data);
bbc_cls();
}
return 0;
}
00000000 23 69 6e 63 6c 75 64 65 20 3c 73 74 64 69 6f 2e |#include <stdio.|
00000010 68 3e 0a 23 69 6e 63 6c 75 64 65 20 3c 73 74 64 |h>.#include <std|
00000020 6c 69 62 2e 68 3e 0a 23 69 6e 63 6c 75 64 65 20 |lib.h>.#include |
00000030 3c 73 74 64 61 72 67 2e 68 3e 0a 23 69 6e 63 6c |<stdarg.h>.#incl|
00000040 75 64 65 20 3c 73 74 72 69 6e 67 2e 68 3e 0a 23 |ude <string.h>.#|
00000050 69 6e 63 6c 75 64 65 20 3c 74 69 6d 65 2e 68 3e |include <time.h>|
00000060 0a 23 69 6e 63 6c 75 64 65 20 3c 6d 61 74 68 2e |.#include <math.|
00000070 68 3e 0a 0a 23 69 6e 63 6c 75 64 65 20 22 77 69 |h>..#include "wi|
00000080 6d 70 74 2e 68 22 0a 23 69 6e 63 6c 75 64 65 20 |mpt.h".#include |
00000090 22 77 65 72 72 2e 68 22 0a 23 69 6e 63 6c 75 64 |"werr.h".#includ|
000000a0 65 20 22 62 62 63 2e 68 22 0a 23 69 6e 63 6c 75 |e "bbc.h".#inclu|
000000b0 64 65 20 22 6f 73 2e 68 22 0a 23 69 6e 63 6c 75 |de "os.h".#inclu|
000000c0 64 65 20 22 66 6c 65 78 2e 68 22 0a 23 69 6e 63 |de "flex.h".#inc|
000000d0 6c 75 64 65 20 22 73 70 72 69 74 65 2e 68 22 0a |lude "sprite.h".|
000000e0 23 69 6e 63 6c 75 64 65 20 22 61 6b 62 64 2e 68 |#include "akbd.h|
000000f0 22 0a 0a 23 69 6e 63 6c 75 64 65 20 22 73 77 69 |"..#include "swi|
00000100 73 2e 68 22 0a 23 64 65 66 69 6e 65 20 6f 6e 65 |s.h".#define one|
00000110 20 28 31 3c 3c 31 36 29 0a 0a 23 64 65 66 69 6e | (1<<16)..#defin|
00000120 65 20 64 69 6d 73 6e 28 66 29 20 28 28 69 6e 74 |e dimsn(f) ((int|
00000130 29 28 66 29 2f 35 2e 30 2b 32 2e 30 29 0a 0a 74 |)(f)/5.0+2.0)..t|
00000140 79 70 65 64 65 66 20 65 6e 75 6d 20 7b 6d 75 6c |ypedef enum {mul|
00000150 74 69 2c 20 6e 6f 6e 6d 75 6c 74 69 7d 20 6d 6f |ti, nonmulti} mo|
00000160 6e 69 74 6f 72 3b 0a 6d 6f 6e 69 74 6f 72 20 6d |nitor;.monitor m|
00000170 6f 6e 69 74 6f 72 74 79 70 65 3b 0a 69 6e 74 20 |onitortype;.int |
00000180 6c 69 6e 65 73 65 70 3b 0a 0a 74 79 70 65 64 65 |linesep;..typede|
00000190 66 20 65 6e 75 6d 20 7b 0a 20 20 6c 65 66 74 3d |f enum {. left=|
000001a0 31 33 36 2c 0a 20 20 72 69 67 68 74 2c 0a 20 20 |136,. right,. |
000001b0 64 6f 77 6e 2c 0a 20 20 75 70 0a 7d 20 73 63 72 |down,. up.} scr|
000001c0 6f 6c 6c 64 69 72 3b 0a 0a 74 79 70 65 64 65 66 |olldir;..typedef|
000001d0 20 73 74 72 75 63 74 20 7b 0a 20 20 69 6e 74 20 | struct {. int |
000001e0 72 65 3b 0a 20 20 69 6e 74 20 69 6d 3b 0a 7d 20 |re;. int im;.} |
000001f0 63 6f 6d 70 6c 65 78 3b 0a 0a 69 6e 74 20 4e 3b |complex;..int N;|
00000200 0a 69 6e 74 20 6c 69 6d 69 74 3b 0a 76 6f 69 64 |.int limit;.void|
00000210 20 2a 64 61 74 61 20 3d 20 30 3b 0a 69 6e 74 20 | *data = 0;.int |
00000220 74 69 6c 74 3d 6f 6e 65 2f 33 3b 20 2f 2a 69 6e |tilt=one/3; /*in|
00000230 20 72 61 6e 67 65 20 30 20 74 6f 20 6f 6e 65 2a | range 0 to one*|
00000240 2f 0a 69 6e 74 20 72 61 6e 67 65 3b 0a 69 6e 74 |/.int range;.int|
00000250 20 73 65 61 6c 65 76 65 6c 3b 0a 69 6e 74 20 70 | sealevel;.int p|
00000260 6f 77 73 63 61 6c 65 3b 0a 69 6e 74 20 6c 69 6e |owscale;.int lin|
00000270 73 63 61 6c 65 3b 0a 69 6e 74 20 6f 72 69 65 6e |scale;.int orien|
00000280 74 2c 20 73 63 72 6f 6c 6c 78 2c 20 73 63 72 6f |t, scrollx, scro|
00000290 6c 6c 79 3b 0a 0a 69 6e 74 20 73 63 72 65 65 6e |lly;..int screen|
000002a0 3d 30 3b 0a 0a 63 68 61 72 20 2a 69 6d 6e 61 6d |=0;..char *imnam|
000002b0 65 20 3d 20 22 73 63 61 70 65 22 3b 0a 73 70 72 |e = "scape";.spr|
000002c0 69 74 65 5f 61 72 65 61 20 2a 73 61 33 32 3d 30 |ite_area *sa32=0|
000002d0 3b 0a 69 6e 74 20 69 6d 33 32 5f 79 6c 3b 0a 69 |;.int im32_yl;.i|
000002e0 6e 74 20 69 6d 33 32 5f 68 65 69 67 68 74 3b 0a |nt im32_height;.|
000002f0 69 6e 74 20 64 7a 3b 0a 69 6e 74 20 6c 78 2c 6c |int dz;.int lx,l|
00000300 79 2c 6c 7a 3b 0a 69 6e 74 20 76 78 2c 76 79 2c |y,lz;.int vx,vy,|
00000310 76 7a 3b 0a 69 6e 74 20 68 78 2c 68 79 2c 68 7a |vz;.int hx,hy,hz|
00000320 3b 0a 0a 73 70 72 69 74 65 5f 61 72 65 61 20 2a |;..sprite_area *|
00000330 64 69 74 68 61 72 65 61 3d 30 3b 0a 63 68 61 72 |ditharea=0;.char|
00000340 20 6c 75 74 5b 33 32 37 36 38 5d 3b 0a 75 6e 73 | lut[32768];.uns|
00000350 69 67 6e 65 64 20 69 6e 74 20 70 61 6c 32 35 36 |igned int pal256|
00000360 64 65 66 5b 32 35 36 5d 3b 0a 75 6e 73 69 67 6e |def[256];.unsign|
00000370 65 64 20 69 6e 74 20 65 61 72 74 68 5b 32 35 36 |ed int earth[256|
00000380 5d 3b 0a 42 4f 4f 4c 20 73 61 76 65 33 32 20 3d |];.BOOL save32 =|
00000390 20 54 52 55 45 3b 0a 42 4f 4f 4c 20 73 61 76 65 | TRUE;.BOOL save|
000003a0 38 20 3d 20 54 52 55 45 3b 0a 74 79 70 65 64 65 |8 = TRUE;.typede|
000003b0 66 20 73 74 72 75 63 74 20 7b 0a 20 20 73 69 67 |f struct {. sig|
000003c0 6e 65 64 20 63 68 61 72 20 65 72 3b 0a 20 20 73 |ned char er;. s|
000003d0 69 67 6e 65 64 20 63 68 61 72 20 65 67 3b 0a 20 |igned char eg;. |
000003e0 20 73 69 67 6e 65 64 20 63 68 61 72 20 65 62 3b | signed char eb;|
000003f0 0a 20 20 73 69 67 6e 65 64 20 63 68 61 72 20 62 |. signed char b|
00000400 79 74 65 34 3b 20 20 2f 2a 72 65 6d 69 6e 64 65 |yte4; /*reminde|
00000410 72 20 74 68 61 74 20 73 74 72 75 63 74 20 77 69 |r that struct wi|
00000420 6c 6c 20 62 65 20 77 6f 72 64 20 73 69 7a 65 64 |ll be word sized|
00000430 2a 2f 0a 7d 20 63 6f 6d 70 61 63 74 5f 64 69 74 |*/.} compact_dit|
00000440 68 65 72 5f 65 72 72 6f 72 3b 0a 0a 42 4f 4f 4c |her_error;..BOOL|
00000450 20 67 65 6e 5f 32 35 36 77 69 74 68 6e 6f 70 61 | gen_256withnopa|
00000460 6c 65 74 74 65 28 66 6c 65 78 5f 70 74 72 20 73 |lette(flex_ptr s|
00000470 61 2c 20 63 68 61 72 20 2a 6e 61 6d 65 29 3b 0a |a, char *name);.|
00000480 65 78 74 65 72 6e 20 76 6f 69 64 20 67 32 35 36 |extern void g256|
00000490 77 6e 70 5f 69 6e 74 65 72 6e 61 6c 28 69 6e 74 |wnp_internal(int|
000004a0 20 78 2c 20 69 6e 74 20 79 2c 20 63 68 61 72 20 | x, int y, char |
000004b0 2a 64 70 2c 20 63 68 61 72 20 2a 73 70 2c 20 63 |*dp, char *sp, c|
000004c0 6f 6d 70 61 63 74 5f 64 69 74 68 65 72 5f 65 72 |ompact_dither_er|
000004d0 72 6f 72 20 2a 64 65 2c 20 63 68 61 72 20 2a 6c |ror *de, char *l|
000004e0 75 74 2c 0a 20 20 20 20 20 20 20 20 20 20 20 20 |ut,. |
000004f0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000500 20 75 6e 73 69 67 6e 65 64 20 69 6e 74 20 2a 70 | unsigned int *p|
00000510 61 6c 32 35 36 64 65 66 29 3b 0a 0a 65 78 74 65 |al256def);..exte|
00000520 72 6e 20 69 6e 74 20 72 62 62 63 69 6e 63 28 69 |rn int rbbcinc(i|
00000530 6e 74 20 72 2c 20 69 6e 74 20 6b 29 3b 0a 65 78 |nt r, int k);.ex|
00000540 74 65 72 6e 20 69 6e 74 20 64 69 76 5f 66 72 61 |tern int div_fra|
00000550 63 31 36 28 69 6e 74 20 6e 75 6d 62 65 72 2c 20 |c16(int number, |
00000560 69 6e 74 20 64 69 76 69 73 6f 72 29 3b 0a 65 78 |int divisor);.ex|
00000570 74 65 72 6e 20 69 6e 74 20 6d 75 6c 5f 66 72 61 |tern int mul_fra|
00000580 63 31 36 28 69 6e 74 20 78 2c 20 69 6e 74 20 61 |c16(int x, int a|
00000590 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 6d 75 |);.extern int mu|
000005a0 6c 5f 66 72 61 63 31 36 63 28 69 6e 74 20 78 2c |l_frac16c(int x,|
000005b0 20 69 6e 74 20 61 29 3b 0a 65 78 74 65 72 6e 20 | int a);.extern |
000005c0 69 6e 74 20 73 71 72 74 5f 66 72 61 63 31 36 28 |int sqrt_frac16(|
000005d0 75 6e 73 69 67 6e 65 64 20 69 6e 74 20 78 29 3b |unsigned int x);|
000005e0 0a 65 78 74 65 72 6e 20 69 6e 74 20 67 61 75 73 |.extern int gaus|
000005f0 73 31 36 28 76 6f 69 64 29 3b 0a 65 78 74 65 72 |s16(void);.exter|
00000600 6e 20 76 6f 69 64 20 73 67 61 75 73 73 31 36 28 |n void sgauss16(|
00000610 69 6e 74 20 73 65 65 64 29 3b 0a 65 78 74 65 72 |int seed);.exter|
00000620 6e 20 69 6e 74 20 63 6f 73 31 36 28 69 6e 74 20 |n int cos16(int |
00000630 61 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 73 |a);.extern int s|
00000640 69 6e 31 36 28 69 6e 74 20 61 29 3b 0a 65 78 74 |in16(int a);.ext|
00000650 65 72 6e 20 69 6e 74 20 65 78 70 31 36 28 69 6e |ern int exp16(in|
00000660 74 20 61 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 |t a);.extern int|
00000670 20 6c 6e 31 36 28 69 6e 74 20 61 29 3b 0a 65 78 | ln16(int a);.ex|
00000680 74 65 72 6e 20 69 6e 74 20 70 6f 77 31 36 28 69 |tern int pow16(i|
00000690 6e 74 20 61 2c 20 69 6e 74 20 62 29 3b 0a 0a 2f |nt a, int b);../|
000006a0 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
*
000006c0 2a 2a 2a 2a 2a 2a 20 66 6f 75 72 69 65 72 20 74 |****** fourier t|
000006d0 72 61 6e 73 66 6f 72 6d 20 66 75 6e 63 74 69 6f |ransform functio|
000006e0 6e 73 20 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |ns *************|
000006f0 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00000700 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2f 0a 0a 42 4f |***********/..BO|
00000710 4f 4c 20 69 6e 76 66 66 74 5f 32 64 28 63 6f 6d |OL invfft_2d(com|
00000720 70 6c 65 78 20 2a 64 61 74 61 2c 20 69 6e 74 20 |plex *data, int |
00000730 6c 69 6d 69 74 29 0a 7b 0a 20 20 69 6e 74 20 69 |limit).{. int i|
00000740 2c 20 6a 3b 0a 20 20 69 6e 74 20 6e 20 3d 20 31 |, j;. int n = 1|
00000750 3c 3c 6c 69 6d 69 74 3b 0a 20 20 63 6f 6d 70 6c |<<limit;. compl|
00000760 65 78 20 74 31 2c 20 2a 64 31 2c 20 2a 64 31 65 |ex t1, *d1, *d1e|
00000770 2c 20 2a 64 32 3b 0a 20 20 69 6e 74 20 6b 2c 20 |, *d2;. int k, |
00000780 6c 3b 0a 20 20 69 6e 74 20 6c 6f 6f 70 3b 0a 20 |l;. int loop;. |
00000790 20 69 6e 74 20 72 6f 74 61 6e 67 6c 65 2c 20 61 | int rotangle, a|
000007a0 6e 67 6c 65 3b 0a 20 20 69 6e 74 20 63 6f 73 2c |ngle;. int cos,|
000007b0 20 73 69 6e 3b 0a 20 20 69 66 20 28 6c 69 6d 69 | sin;. if (limi|
000007c0 74 3c 33 20 7c 7c 20 6c 69 6d 69 74 3e 33 32 29 |t<3 || limit>32)|
000007d0 20 72 65 74 75 72 6e 20 46 41 4c 53 45 3b 0a 0a | return FALSE;..|
000007e0 20 20 66 6f 72 20 28 69 3d 6a 3d 30 3b 20 6a 3c | for (i=j=0; j<|
000007f0 6e 3b 20 69 3d 72 62 62 63 69 6e 63 28 69 2c 20 |n; i=rbbcinc(i, |
00000800 6c 69 6d 69 74 29 2c 20 6a 2b 2b 29 20 69 66 20 |limit), j++) if |
00000810 28 69 3e 6a 29 20 7b 0a 20 20 20 20 66 6f 72 20 |(i>j) {. for |
00000820 28 64 31 3d 64 61 74 61 2b 28 69 3c 3c 6c 69 6d |(d1=data+(i<<lim|
00000830 69 74 29 2c 20 64 32 3d 64 61 74 61 2b 28 6a 3c |it), d2=data+(j<|
00000840 3c 6c 69 6d 69 74 29 2c 20 64 31 65 3d 64 31 2b |<limit), d1e=d1+|
00000850 6e 3b 20 64 31 3c 64 31 65 3b 29 20 7b 0a 20 20 |n; d1<d1e;) {. |
00000860 20 20 20 20 74 31 20 3d 20 2a 64 31 3b 0a 20 20 | t1 = *d1;. |
00000870 20 20 20 20 2a 64 31 2b 2b 20 3d 20 2a 64 32 3b | *d1++ = *d2;|
00000880 0a 20 20 20 20 20 20 2a 64 32 2b 2b 20 3d 20 74 |. *d2++ = t|
00000890 31 3b 0a 20 20 20 20 7d 0a 20 20 7d 0a 20 20 72 |1;. }. }. r|
000008a0 6f 74 61 6e 67 6c 65 20 3d 20 2d 32 2a 6f 6e 65 |otangle = -2*one|
000008b0 3b 0a 20 20 66 6f 72 20 28 6b 3d 6c 6f 6f 70 3d |;. for (k=loop=|
000008c0 31 3b 20 6c 6f 6f 70 3c 3d 6c 69 6d 69 74 3b 20 |1; loop<=limit; |
000008d0 6c 6f 6f 70 2b 2b 2c 20 6b 3c 3c 3d 31 2c 20 72 |loop++, k<<=1, r|
000008e0 6f 74 61 6e 67 6c 65 2f 3d 32 29 20 7b 0a 20 20 |otangle/=2) {. |
000008f0 20 20 66 6f 72 20 28 61 6e 67 6c 65 3d 30 2c 20 | for (angle=0, |
00000900 73 69 6e 3d 30 2c 20 63 6f 73 3d 6f 6e 65 2c 20 |sin=0, cos=one, |
00000910 69 3d 31 3b 20 69 3c 3d 6b 3b 20 69 2b 2b 2c 20 |i=1; i<=k; i++, |
00000920 61 6e 67 6c 65 2b 3d 72 6f 74 61 6e 67 6c 65 2c |angle+=rotangle,|
00000930 20 73 69 6e 3d 73 69 6e 31 36 28 61 6e 67 6c 65 | sin=sin16(angle|
00000940 29 2c 20 63 6f 73 3d 63 6f 73 31 36 28 61 6e 67 |), cos=cos16(ang|
00000950 6c 65 29 29 20 7b 0a 20 20 20 20 20 20 66 6f 72 |le)) {. for|
00000960 20 28 6a 3d 69 2d 31 3b 20 6a 3c 6e 3b 20 6a 2b | (j=i-1; j<n; j+|
00000970 3d 6b 3c 3c 31 29 20 7b 0a 20 20 20 20 20 20 20 |=k<<1) {. |
00000980 20 6c 20 3d 20 6a 2b 6b 3b 0a 20 20 20 20 20 20 | l = j+k;. |
00000990 20 20 66 6f 72 20 28 64 31 3d 64 61 74 61 2b 28 | for (d1=data+(|
000009a0 6c 3c 3c 6c 69 6d 69 74 29 2c 20 64 32 3d 64 61 |l<<limit), d2=da|
000009b0 74 61 2b 28 6a 3c 3c 6c 69 6d 69 74 29 2c 20 64 |ta+(j<<limit), d|
000009c0 31 65 3d 64 31 2b 6e 3b 20 64 31 3c 64 31 65 3b |1e=d1+n; d1<d1e;|
000009d0 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 74 31 |) {. t1|
000009e0 2e 72 65 20 3d 20 28 6d 75 6c 5f 66 72 61 63 31 |.re = (mul_frac1|
000009f0 36 63 28 64 31 2d 3e 72 65 2c 20 63 6f 73 29 2d |6c(d1->re, cos)-|
00000a00 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e |mul_frac16c(d1->|
00000a10 69 6d 2c 20 73 69 6e 29 29 3b 0a 20 20 20 20 20 |im, sin));. |
00000a20 20 20 20 20 20 74 31 2e 69 6d 20 3d 20 28 6d 75 | t1.im = (mu|
00000a30 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e 72 65 |l_frac16c(d1->re|
00000a40 2c 20 73 69 6e 29 2b 6d 75 6c 5f 66 72 61 63 31 |, sin)+mul_frac1|
00000a50 36 63 28 64 31 2d 3e 69 6d 2c 20 63 6f 73 29 29 |6c(d1->im, cos))|
00000a60 3b 0a 20 20 20 20 20 20 20 20 20 20 64 31 2d 3e |;. d1->|
00000a70 72 65 20 20 20 20 3d 20 64 32 2d 3e 72 65 2d 74 |re = d2->re-t|
00000a80 31 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 20 20 |1.re;. |
00000a90 64 31 2b 2b 2d 3e 69 6d 20 20 3d 20 64 32 2d 3e |d1++->im = d2->|
00000aa0 69 6d 2d 74 31 2e 69 6d 3b 0a 20 20 20 20 20 20 |im-t1.im;. |
00000ab0 20 20 20 20 64 32 2d 3e 72 65 20 20 20 2b 3d 20 | d2->re += |
00000ac0 74 31 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 20 |t1.re;. |
00000ad0 20 64 32 2b 2b 2d 3e 69 6d 20 2b 3d 20 74 31 2e | d2++->im += t1.|
00000ae0 69 6d 3b 0a 20 20 20 20 20 20 20 20 7d 0a 20 20 |im;. }. |
00000af0 20 20 20 20 7d 0a 20 20 20 20 7d 0a 20 20 7d 0a | }. }. }.|
00000b00 0a 20 20 66 6f 72 20 28 69 3d 6a 3d 30 3b 20 6a |. for (i=j=0; j|
00000b10 3c 6e 3b 20 69 3d 72 62 62 63 69 6e 63 28 69 2c |<n; i=rbbcinc(i,|
00000b20 20 6c 69 6d 69 74 29 2c 20 6a 2b 2b 29 20 69 66 | limit), j++) if|
00000b30 20 28 69 3e 6a 29 20 7b 0a 20 20 20 20 66 6f 72 | (i>j) {. for|
00000b40 20 28 64 31 3d 64 61 74 61 2b 69 2c 20 64 32 3d | (d1=data+i, d2=|
00000b50 64 61 74 61 2b 6a 2c 20 64 31 65 3d 64 31 2b 6e |data+j, d1e=d1+n|
00000b60 2a 6e 3b 20 64 31 3c 64 31 65 3b 20 64 31 2b 3d |*n; d1<d1e; d1+=|
00000b70 6e 2c 20 64 32 2b 3d 6e 29 20 7b 0a 20 20 20 20 |n, d2+=n) {. |
00000b80 20 20 74 31 20 3d 20 2a 64 31 3b 0a 20 20 20 20 | t1 = *d1;. |
00000b90 20 20 2a 64 31 20 3d 20 2a 64 32 3b 0a 20 20 20 | *d1 = *d2;. |
00000ba0 20 20 20 2a 64 32 20 3d 20 74 31 3b 0a 20 20 20 | *d2 = t1;. |
00000bb0 20 7d 0a 20 20 7d 0a 20 20 72 6f 74 61 6e 67 6c | }. }. rotangl|
00000bc0 65 20 3d 20 2d 32 2a 6f 6e 65 3b 0a 20 20 66 6f |e = -2*one;. fo|
00000bd0 72 20 28 6b 3d 6c 6f 6f 70 3d 31 3b 20 6c 6f 6f |r (k=loop=1; loo|
00000be0 70 3c 3d 6c 69 6d 69 74 3b 20 6c 6f 6f 70 2b 2b |p<=limit; loop++|
00000bf0 2c 20 6b 3c 3c 3d 31 2c 20 72 6f 74 61 6e 67 6c |, k<<=1, rotangl|
00000c00 65 2f 3d 32 29 20 7b 0a 20 20 20 20 66 6f 72 20 |e/=2) {. for |
00000c10 28 61 6e 67 6c 65 3d 30 2c 20 73 69 6e 3d 30 2c |(angle=0, sin=0,|
00000c20 20 63 6f 73 3d 6f 6e 65 2c 20 69 3d 31 3b 20 69 | cos=one, i=1; i|
00000c30 3c 3d 6b 3b 20 69 2b 2b 2c 20 61 6e 67 6c 65 2b |<=k; i++, angle+|
00000c40 3d 72 6f 74 61 6e 67 6c 65 2c 20 73 69 6e 3d 73 |=rotangle, sin=s|
00000c50 69 6e 31 36 28 61 6e 67 6c 65 29 2c 20 63 6f 73 |in16(angle), cos|
00000c60 3d 63 6f 73 31 36 28 61 6e 67 6c 65 29 29 20 7b |=cos16(angle)) {|
00000c70 0a 20 20 20 20 20 20 66 6f 72 20 28 6a 3d 69 2d |. for (j=i-|
00000c80 31 3b 20 6a 3c 6e 3b 20 6a 2b 3d 6b 3c 3c 31 29 |1; j<n; j+=k<<1)|
00000c90 20 7b 0a 20 20 20 20 20 20 20 20 6c 20 3d 20 6a | {. l = j|
00000ca0 2b 6b 3b 0a 20 20 20 20 20 20 20 20 66 6f 72 20 |+k;. for |
00000cb0 28 64 31 3d 64 61 74 61 2b 6c 2c 20 64 32 3d 64 |(d1=data+l, d2=d|
00000cc0 61 74 61 2b 6a 2c 20 64 31 65 3d 64 31 2b 6e 2a |ata+j, d1e=d1+n*|
00000cd0 6e 3b 20 64 31 3c 64 31 65 3b 20 64 31 2b 3d 6e |n; d1<d1e; d1+=n|
00000ce0 2c 20 64 32 2b 3d 6e 29 20 7b 0a 20 20 20 20 20 |, d2+=n) {. |
00000cf0 20 20 20 20 20 74 31 2e 72 65 20 3d 20 28 6d 75 | t1.re = (mu|
00000d00 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e 72 65 |l_frac16c(d1->re|
00000d10 2c 20 63 6f 73 29 2d 6d 75 6c 5f 66 72 61 63 31 |, cos)-mul_frac1|
00000d20 36 63 28 64 31 2d 3e 69 6d 2c 20 73 69 6e 29 29 |6c(d1->im, sin))|
00000d30 3b 0a 20 20 20 20 20 20 20 20 20 20 74 31 2e 69 |;. t1.i|
00000d40 6d 20 3d 20 28 6d 75 6c 5f 66 72 61 63 31 36 63 |m = (mul_frac16c|
00000d50 28 64 31 2d 3e 72 65 2c 20 73 69 6e 29 2b 6d 75 |(d1->re, sin)+mu|
00000d60 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e 69 6d |l_frac16c(d1->im|
00000d70 2c 20 63 6f 73 29 29 3b 0a 20 20 20 20 20 20 20 |, cos));. |
00000d80 20 20 20 64 31 2d 3e 72 65 20 20 3d 20 64 32 2d | d1->re = d2-|
00000d90 3e 72 65 2d 74 31 2e 72 65 3b 0a 20 20 20 20 20 |>re-t1.re;. |
00000da0 20 20 20 20 20 64 31 2d 3e 69 6d 20 20 3d 20 64 | d1->im = d|
00000db0 32 2d 3e 69 6d 2d 74 31 2e 69 6d 3b 0a 20 20 20 |2->im-t1.im;. |
00000dc0 20 20 20 20 20 20 20 64 32 2d 3e 72 65 20 2b 3d | d2->re +=|
00000dd0 20 74 31 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 | t1.re;. |
00000de0 20 20 64 32 2d 3e 69 6d 20 2b 3d 20 74 31 2e 69 | d2->im += t1.i|
00000df0 6d 3b 0a 20 20 20 20 20 20 20 20 7d 0a 20 20 20 |m;. }. |
00000e00 20 20 20 7d 0a 20 20 20 20 7d 0a 20 20 7d 0a 0a | }. }. }..|
00000e10 20 20 72 65 74 75 72 6e 20 54 52 55 45 3b 0a 7d | return TRUE;.}|
00000e20 0a 0a 2f 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |../*************|
00000e30 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00000e40 2a 2a 2a 2a 2a 2a 20 64 61 74 61 20 67 65 6e 65 |****** data gene|
00000e50 72 61 74 69 6f 6e 2f 6d 61 6e 69 70 75 6c 61 74 |ration/manipulat|
00000e60 69 6f 6e 20 66 75 6e 63 74 69 6f 6e 73 20 20 2a |ion functions *|
00000e70 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00000e80 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2f 0a |**************/.|
00000e90 0a 2f 2a 54 68 65 20 66 6f 6c 6c 6f 77 69 6e 67 |./*The following|
00000ea0 20 61 73 73 69 67 6e 73 20 76 61 6c 75 65 73 20 | assigns values |
00000eb0 74 6f 20 32 64 20 64 61 74 61 20 61 72 72 61 79 |to 2d data array|
00000ec0 20 77 69 74 68 20 73 75 69 74 61 62 6c 65 20 70 | with suitable p|
00000ed0 72 6f 70 65 72 74 69 65 73 2e 2a 2f 0a 2f 2a 49 |roperties.*/./*I|
00000ee0 74 20 64 6f 65 73 20 68 6f 77 65 76 65 72 20 64 |t does however d|
00000ef0 6f 20 73 6f 20 62 79 20 73 63 61 6e 6e 69 6e 67 |o so by scanning|
00000f00 20 74 68 65 20 61 72 72 61 79 20 69 6e 20 61 20 | the array in a |
00000f10 6e 6f 6e 2d 6f 62 76 69 6f 75 73 20 6f 72 64 65 |non-obvious orde|
00000f20 72 2e 0a 20 20 54 68 69 73 20 68 61 73 20 62 65 |r.. This has be|
00000f30 65 6e 20 63 68 6f 73 65 6e 20 73 6f 20 74 68 61 |en chosen so tha|
00000f40 74 20 69 66 20 6e 20 69 73 20 69 6e 63 72 65 61 |t if n is increa|
00000f50 73 65 64 20 62 79 20 64 6f 75 62 6c 69 6e 67 2c |sed by doubling,|
00000f60 20 74 68 65 6e 20 74 68 6f 73 65 20 72 61 6e 64 | then those rand|
00000f70 6f 6d 20 76 61 6c 75 65 73 20 75 73 65 64 20 66 |om values used f|
00000f80 6f 72 20 74 68 65 0a 20 20 6c 61 73 74 20 67 72 |or the. last gr|
00000f90 69 64 20 77 69 6c 6c 20 62 65 20 6e 6f 77 20 77 |id will be now w|
00000fa0 72 69 74 74 65 6e 20 74 6f 20 74 68 65 20 6c 6f |ritten to the lo|
00000fb0 77 20 66 72 65 71 75 65 6e 63 79 20 71 75 61 64 |w frequency quad|
00000fc0 72 61 6e 74 20 6f 66 20 74 68 65 20 6e 65 77 20 |rant of the new |
00000fd0 67 72 69 64 20 28 69 66 20 73 61 6d 65 20 73 65 |grid (if same se|
00000fe0 65 64 20 75 73 65 64 29 2c 0a 20 20 73 75 69 74 |ed used),. suit|
00000ff0 61 62 6c 79 20 72 65 20 66 69 6c 74 65 72 65 64 |ably re filtered|
00001000 2e 20 54 68 69 73 20 73 68 6f 75 6c 64 20 68 61 |. This should ha|
00001010 76 65 20 65 66 66 65 63 74 20 6f 66 20 67 69 76 |ve effect of giv|
00001020 69 6e 67 20 61 20 64 65 67 72 65 65 20 6f 66 20 |ing a degree of |
00001030 63 6f 6e 74 69 6e 75 69 74 79 20 74 6f 20 74 68 |continuity to th|
00001040 65 20 63 68 61 6e 67 65 20 69 6e 20 6e 2c 0a 20 |e change in n,. |
00001050 20 73 69 6e 63 65 20 66 6f 72 6d 65 72 20 66 72 | since former fr|
00001060 65 71 75 65 6e 63 69 65 73 20 61 72 65 20 63 61 |equencies are ca|
00001070 72 72 69 65 64 20 6f 76 65 72 20 74 6f 20 74 68 |rried over to th|
00001080 65 20 6c 6f 77 20 66 72 65 71 75 65 6e 63 69 65 |e low frequencie|
00001090 73 20 6f 66 20 74 68 65 20 68 69 67 68 65 72 20 |s of the higher |
000010a0 72 65 73 20 67 72 69 64 2e 0a 2a 2f 0a 42 4f 4f |res grid..*/.BOO|
000010b0 4c 20 61 73 73 69 67 6e 5f 66 74 28 63 6f 6d 70 |L assign_ft(comp|
000010c0 6c 65 78 20 2a 64 61 74 61 2c 20 69 6e 74 20 6c |lex *data, int l|
000010d0 69 6d 69 74 2c 20 64 6f 75 62 6c 65 20 64 69 6d |imit, double dim|
000010e0 2c 20 69 6e 74 20 73 65 65 64 29 0a 7b 0a 20 20 |, int seed).{. |
000010f0 69 6e 74 20 6e 3d 31 3c 3c 6c 69 6d 69 74 3b 0a |int n=1<<limit;.|
00001100 20 20 69 6e 74 20 69 2c 20 69 63 2c 20 6a 2c 20 | int i, ic, j, |
00001110 69 32 2c 20 6a 32 2c 20 6b 2c 20 6c 3b 0a 20 20 |i2, j2, k, l;. |
00001120 63 6f 6d 70 6c 65 78 20 2a 64 31 2c 20 2a 64 32 |complex *d1, *d2|
00001130 2c 20 2a 64 3b 0a 20 20 69 6e 74 20 70 68 61 73 |, *d;. int phas|
00001140 65 2c 20 72 61 64 3b 0a 20 20 69 6e 74 20 69 6d |e, rad;. int im|
00001150 2c 20 72 65 3b 0a 20 20 69 6e 74 20 6d 62 70 6f |, re;. int mbpo|
00001160 62 74 3b 0a 0a 20 20 73 72 61 6e 64 28 2d 73 65 |bt;.. srand(-se|
00001170 65 64 29 3b 0a 20 20 73 67 61 75 73 73 31 36 28 |ed);. sgauss16(|
00001180 2d 73 65 65 64 29 3b 0a 0a 20 20 6d 62 70 6f 62 |-seed);.. mbpob|
00001190 74 20 3d 20 28 69 6e 74 29 28 2d 28 34 2d 64 69 |t = (int)(-(4-di|
000011a0 6d 29 2a 6f 6e 65 29 3b 0a 20 20 64 61 74 61 2d |m)*one);. data-|
000011b0 3e 72 65 20 3d 20 64 61 74 61 2d 3e 69 6d 20 3d |>re = data->im =|
000011c0 20 30 3b 0a 20 20 66 6f 72 20 28 69 3d 31 3b 20 | 0;. for (i=1; |
000011d0 69 3c 3d 6e 2f 32 3b 20 69 2b 2b 29 20 7b 0a 20 |i<=n/2; i++) {. |
000011e0 20 20 20 6a 32 20 3d 20 69 32 20 3d 20 69 2a 69 | j2 = i2 = i*i|
000011f0 3b 0a 20 20 20 20 6b 20 20 3d 20 6e 2d 69 3b 0a |;. k = n-i;.|
00001200 20 20 20 20 64 31 20 3d 20 64 61 74 61 2b 69 2a | d1 = data+i*|
00001210 6e 3b 0a 20 20 20 20 64 32 20 3d 20 64 61 74 61 |n;. d2 = data|
00001220 2b 6b 2a 6e 3b 0a 20 20 20 20 66 6f 72 20 28 6a |+k*n;. for (j|
00001230 3d 30 3b 20 6a 3c 3d 69 3b 20 6a 2b 2b 29 20 7b |=0; j<=i; j++) {|
00001240 0a 20 20 20 20 20 20 6c 20 3d 20 6a 3d 3d 30 20 |. l = j==0 |
00001250 3f 20 30 20 3a 20 6e 2d 6a 3b 0a 20 20 20 20 20 |? 0 : n-j;. |
00001260 20 70 68 61 73 65 20 3d 20 72 61 6e 64 28 29 3e | phase = rand()>|
00001270 3e 31 33 3b 0a 20 20 20 20 20 20 72 61 64 20 3d |>13;. rad =|
00001280 20 6d 75 6c 5f 66 72 61 63 31 36 28 70 6f 77 31 | mul_frac16(pow1|
00001290 36 28 73 71 72 74 5f 66 72 61 63 31 36 28 28 69 |6(sqrt_frac16((i|
000012a0 32 2b 6a 2a 6a 29 2a 28 6f 6e 65 2f 34 29 29 2a |2+j*j)*(one/4))*|
000012b0 32 2c 20 6d 62 70 6f 62 74 29 2c 20 67 61 75 73 |2, mbpobt), gaus|
000012c0 73 31 36 28 29 29 3b 0a 20 20 20 20 20 20 72 65 |s16());. re|
000012d0 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 28 72 | = mul_frac16c(r|
000012e0 61 64 2c 20 63 6f 73 31 36 28 70 68 61 73 65 29 |ad, cos16(phase)|
000012f0 29 3b 0a 20 20 20 20 20 20 69 6d 20 3d 20 6d 75 |);. im = mu|
00001300 6c 5f 66 72 61 63 31 36 63 28 72 61 64 2c 20 73 |l_frac16c(rad, s|
00001310 69 6e 31 36 28 70 68 61 73 65 29 29 3b 0a 20 20 |in16(phase));. |
00001320 20 20 20 20 64 20 3d 20 64 31 2b 6a 3b 0a 20 20 | d = d1+j;. |
00001330 20 20 20 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a | d->re = re;.|
00001340 20 20 20 20 20 20 64 2d 3e 69 6d 20 3d 20 69 6d | d->im = im|
00001350 3b 0a 20 20 20 20 20 20 64 20 3d 20 64 32 2b 6c |;. d = d2+l|
00001360 3b 0a 20 20 20 20 20 20 64 2d 3e 72 65 20 3d 20 |;. d->re = |
00001370 72 65 3b 0a 20 20 20 20 20 20 64 2d 3e 69 6d 20 |re;. d->im |
00001380 3d 20 2d 69 6d 3b 0a 20 20 20 20 7d 0a 20 20 20 |= -im;. }. |
00001390 20 6a 20 3d 20 69 3b 0a 20 20 20 20 6c 20 3d 20 | j = i;. l = |
000013a0 6e 2d 6a 3b 0a 20 20 20 20 64 31 20 3d 20 64 61 |n-j;. d1 = da|
000013b0 74 61 2b 6a 3b 0a 20 20 20 20 64 32 20 3d 20 64 |ta+j;. d2 = d|
000013c0 61 74 61 2b 6c 3b 0a 20 20 20 20 66 6f 72 20 28 |ata+l;. for (|
000013d0 69 63 3d 69 2d 31 3b 20 69 63 3e 3d 30 3b 20 69 |ic=i-1; ic>=0; i|
000013e0 63 2d 2d 29 20 7b 0a 20 20 20 20 20 20 6b 20 3d |c--) {. k =|
000013f0 20 69 63 3d 3d 30 20 3f 20 30 20 3a 20 6e 2d 69 | ic==0 ? 0 : n-i|
00001400 63 3b 0a 20 20 20 20 20 20 70 68 61 73 65 20 3d |c;. phase =|
00001410 20 72 61 6e 64 28 29 3e 3e 31 33 3b 0a 20 20 20 | rand()>>13;. |
00001420 20 20 20 72 61 64 20 3d 20 6d 75 6c 5f 66 72 61 | rad = mul_fra|
00001430 63 31 36 28 70 6f 77 31 36 28 73 71 72 74 5f 66 |c16(pow16(sqrt_f|
00001440 72 61 63 31 36 28 28 69 63 2a 69 63 2b 6a 32 29 |rac16((ic*ic+j2)|
00001450 2a 28 6f 6e 65 2f 34 29 29 2a 32 2c 20 6d 62 70 |*(one/4))*2, mbp|
00001460 6f 62 74 29 2c 20 67 61 75 73 73 31 36 28 29 29 |obt), gauss16())|
00001470 3b 0a 20 20 20 20 20 20 72 65 20 3d 20 6d 75 6c |;. re = mul|
00001480 5f 66 72 61 63 31 36 63 28 72 61 64 2c 20 63 6f |_frac16c(rad, co|
00001490 73 31 36 28 70 68 61 73 65 29 29 3b 0a 20 20 20 |s16(phase));. |
000014a0 20 20 20 69 6d 20 3d 20 6d 75 6c 5f 66 72 61 63 | im = mul_frac|
000014b0 31 36 63 28 72 61 64 2c 20 73 69 6e 31 36 28 70 |16c(rad, sin16(p|
000014c0 68 61 73 65 29 29 3b 0a 20 20 20 20 20 20 64 20 |hase));. d |
000014d0 3d 20 64 31 2b 69 63 2a 6e 3b 0a 20 20 20 20 20 |= d1+ic*n;. |
000014e0 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 20 20 | d->re = re;. |
000014f0 20 20 20 64 2d 3e 69 6d 20 3d 20 69 6d 3b 0a 20 | d->im = im;. |
00001500 20 20 20 20 20 64 20 3d 20 64 32 2b 6b 2a 6e 3b | d = d2+k*n;|
00001510 0a 20 20 20 20 20 20 64 2d 3e 72 65 20 3d 20 72 |. d->re = r|
00001520 65 3b 0a 20 20 20 20 20 20 64 2d 3e 69 6d 20 3d |e;. d->im =|
00001530 20 2d 69 6d 3b 0a 20 20 20 20 7d 0a 20 20 20 20 | -im;. }. |
00001540 69 66 20 28 69 3c 6e 2f 32 29 20 7b 0a 20 20 20 |if (i<n/2) {. |
00001550 20 20 20 69 63 20 3d 20 6e 2d 69 3b 0a 20 20 20 | ic = n-i;. |
00001560 20 20 20 6b 20 3d 20 69 3b 0a 20 20 20 20 20 20 | k = i;. |
00001570 64 31 20 3d 20 64 61 74 61 2b 69 63 2a 6e 3b 0a |d1 = data+ic*n;.|
00001580 20 20 20 20 20 20 64 32 20 3d 20 64 61 74 61 2b | d2 = data+|
00001590 6b 2a 6e 3b 0a 20 20 20 20 20 20 66 6f 72 20 28 |k*n;. for (|
000015a0 6a 3d 31 3b 20 6a 3c 3d 69 3b 20 6a 2b 2b 29 20 |j=1; j<=i; j++) |
000015b0 7b 0a 20 20 20 20 20 20 20 20 6c 20 3d 20 6e 2d |{. l = n-|
000015c0 6a 3b 0a 20 20 20 20 20 20 20 20 70 68 61 73 65 |j;. phase|
000015d0 20 3d 20 72 61 6e 64 28 29 3e 3e 31 33 3b 0a 20 | = rand()>>13;. |
000015e0 20 20 20 20 20 20 20 72 61 64 20 3d 20 6d 75 6c | rad = mul|
000015f0 5f 66 72 61 63 31 36 28 70 6f 77 31 36 28 73 71 |_frac16(pow16(sq|
00001600 72 74 5f 66 72 61 63 31 36 28 28 69 32 2b 6a 2a |rt_frac16((i2+j*|
00001610 6a 29 2a 28 6f 6e 65 2f 34 29 29 2a 32 2c 20 6d |j)*(one/4))*2, m|
00001620 62 70 6f 62 74 29 2c 20 67 61 75 73 73 31 36 28 |bpobt), gauss16(|
00001630 29 29 3b 0a 20 20 20 20 20 20 20 20 72 65 20 3d |));. re =|
00001640 20 6d 75 6c 5f 66 72 61 63 31 36 63 28 72 61 64 | mul_frac16c(rad|
00001650 2c 20 63 6f 73 31 36 28 70 68 61 73 65 29 29 3b |, cos16(phase));|
00001660 0a 20 20 20 20 20 20 20 20 69 6d 20 3d 20 6d 75 |. im = mu|
00001670 6c 5f 66 72 61 63 31 36 63 28 72 61 64 2c 20 73 |l_frac16c(rad, s|
00001680 69 6e 31 36 28 70 68 61 73 65 29 29 3b 0a 20 20 |in16(phase));. |
00001690 20 20 20 20 20 20 64 20 3d 20 64 31 2b 6a 3b 0a | d = d1+j;.|
000016a0 20 20 20 20 20 20 20 20 64 2d 3e 72 65 20 3d 20 | d->re = |
000016b0 72 65 3b 0a 20 20 20 20 20 20 20 20 64 2d 3e 69 |re;. d->i|
000016c0 6d 20 3d 20 69 6d 3b 0a 20 20 20 20 20 20 20 20 |m = im;. |
000016d0 64 20 3d 20 64 32 2b 6c 3b 0a 20 20 20 20 20 20 |d = d2+l;. |
000016e0 20 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 20 | d->re = re;. |
000016f0 20 20 20 20 20 20 64 2d 3e 69 6d 20 3d 20 2d 69 | d->im = -i|
00001700 6d 3b 0a 20 20 20 20 20 20 7d 0a 20 20 20 20 20 |m;. }. |
00001710 20 6a 3d 69 3b 0a 20 20 20 20 20 20 6c 20 3d 20 | j=i;. l = |
00001720 6e 2d 6a 3b 0a 20 20 20 20 20 20 64 31 20 3d 20 |n-j;. d1 = |
00001730 64 61 74 61 2b 6a 3b 0a 20 20 20 20 20 20 64 32 |data+j;. d2|
00001740 20 3d 20 64 61 74 61 2b 6c 3b 0a 20 20 20 20 20 | = data+l;. |
00001750 20 66 6f 72 20 28 69 63 2b 2b 3b 20 69 63 3c 6e | for (ic++; ic<n|
00001760 3b 20 69 63 2b 2b 29 20 7b 0a 20 20 20 20 20 20 |; ic++) {. |
00001770 20 20 6b 20 3d 20 6e 2d 69 63 3b 0a 20 20 20 20 | k = n-ic;. |
00001780 20 20 20 20 70 68 61 73 65 20 3d 20 72 61 6e 64 | phase = rand|
00001790 28 29 3e 3e 31 33 3b 0a 20 20 20 20 20 20 20 20 |()>>13;. |
000017a0 72 61 64 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 |rad = mul_frac16|
000017b0 28 70 6f 77 31 36 28 73 71 72 74 5f 66 72 61 63 |(pow16(sqrt_frac|
000017c0 31 36 28 28 6b 2a 6b 2b 6a 32 29 2a 28 6f 6e 65 |16((k*k+j2)*(one|
000017d0 2f 34 29 29 2a 32 2c 20 6d 62 70 6f 62 74 29 2c |/4))*2, mbpobt),|
000017e0 20 67 61 75 73 73 31 36 28 29 29 3b 0a 20 20 20 | gauss16());. |
000017f0 20 20 20 20 20 72 65 20 3d 20 6d 75 6c 5f 66 72 | re = mul_fr|
00001800 61 63 31 36 63 28 72 61 64 2c 20 63 6f 73 31 36 |ac16c(rad, cos16|
00001810 28 70 68 61 73 65 29 29 3b 0a 20 20 20 20 20 20 |(phase));. |
00001820 20 20 69 6d 20 3d 20 6d 75 6c 5f 66 72 61 63 31 | im = mul_frac1|
00001830 36 63 28 72 61 64 2c 20 73 69 6e 31 36 28 70 68 |6c(rad, sin16(ph|
00001840 61 73 65 29 29 3b 0a 20 20 20 20 20 20 20 20 64 |ase));. d|
00001850 20 3d 20 64 31 2b 69 63 2a 6e 3b 0a 20 20 20 20 | = d1+ic*n;. |
00001860 20 20 20 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a | d->re = re;.|
00001870 20 20 20 20 20 20 20 20 64 2d 3e 69 6d 20 3d 20 | d->im = |
00001880 69 6d 3b 0a 20 20 20 20 20 20 20 20 64 20 3d 20 |im;. d = |
00001890 64 32 2b 6b 2a 6e 3b 0a 20 20 20 20 20 20 20 20 |d2+k*n;. |
000018a0 64 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 20 20 20 |d->re = re;. |
000018b0 20 20 20 20 64 2d 3e 69 6d 20 3d 20 2d 69 6d 3b | d->im = -im;|
000018c0 0a 20 20 20 20 20 20 7d 0a 20 20 20 20 7d 0a 20 |. }. }. |
000018d0 20 7d 0a 20 20 28 64 61 74 61 2b 6e 2a 6e 2f 32 | }. (data+n*n/2|
000018e0 29 2d 3e 69 6d 20 3d 20 28 64 61 74 61 2b 6e 2f |)->im = (data+n/|
000018f0 32 29 2d 3e 69 6d 20 3d 20 28 64 61 74 61 2b 6e |2)->im = (data+n|
00001900 2a 28 6e 2b 31 29 2f 32 29 2d 3e 69 6d 20 3d 20 |*(n+1)/2)->im = |
00001910 30 3b 0a 0a 20 20 72 65 74 75 72 6e 20 54 52 55 |0;.. return TRU|
00001920 45 3b 0a 7d 0a 0a 42 4f 4f 4c 20 6e 6f 72 6d 61 |E;.}..BOOL norma|
00001930 6c 69 73 65 5f 61 6e 64 5f 63 6f 6d 70 61 63 74 |lise_and_compact|
00001940 28 76 6f 69 64 20 2a 2a 64 61 74 61 2c 20 69 6e |(void **data, in|
00001950 74 20 6c 69 6d 69 74 29 0a 7b 0a 20 20 69 6e 74 |t limit).{. int|
00001960 20 6e 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 | n=1<<limit;. i|
00001970 6e 74 20 6d 69 6e 3d 31 3c 3c 33 30 3b 0a 20 20 |nt min=1<<30;. |
00001980 69 6e 74 20 6d 61 78 3d 2d 6d 69 6e 3b 0a 20 20 |int max=-min;. |
00001990 69 6e 74 20 73 63 61 6c 65 2c 20 76 3b 0a 20 20 |int scale, v;. |
000019a0 63 6f 6d 70 6c 65 78 20 2a 64 31 2c 20 2a 64 31 |complex *d1, *d1|
000019b0 65 3b 0a 20 20 69 6e 74 20 2a 64 32 3b 0a 0a 20 |e;. int *d2;.. |
000019c0 20 66 6f 72 20 28 64 31 3d 2a 64 61 74 61 2c 20 | for (d1=*data, |
000019d0 64 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 3c 64 |d1e=d1+n*n; d1<d|
000019e0 31 65 3b 20 64 31 2b 2b 29 20 7b 0a 20 20 20 20 |1e; d1++) {. |
000019f0 69 66 20 28 64 31 2d 3e 72 65 3c 6d 69 6e 29 20 |if (d1->re<min) |
00001a00 6d 69 6e 3d 64 31 2d 3e 72 65 3b 0a 20 20 20 20 |min=d1->re;. |
00001a10 69 66 20 28 64 31 2d 3e 72 65 3e 6d 61 78 29 20 |if (d1->re>max) |
00001a20 6d 61 78 3d 64 31 2d 3e 72 65 3b 0a 20 20 7d 0a |max=d1->re;. }.|
00001a30 20 20 73 63 61 6c 65 20 3d 20 64 69 76 5f 66 72 | scale = div_fr|
00001a40 61 63 31 36 28 6f 6e 65 2c 20 6d 61 78 2d 6d 69 |ac16(one, max-mi|
00001a50 6e 29 3b 0a 20 20 66 6f 72 20 28 64 32 3d 2a 64 |n);. for (d2=*d|
00001a60 61 74 61 2c 20 64 31 3d 2a 64 61 74 61 2c 20 64 |ata, d1=*data, d|
00001a70 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 3c 64 31 |1e=d1+n*n; d1<d1|
00001a80 65 3b 20 64 31 2b 2b 2c 20 64 32 2b 2b 29 20 7b |e; d1++, d2++) {|
00001a90 0a 20 20 20 20 76 20 3d 20 6d 75 6c 5f 66 72 61 |. v = mul_fra|
00001aa0 63 31 36 28 64 31 2d 3e 72 65 2d 6d 69 6e 2c 20 |c16(d1->re-min, |
00001ab0 73 63 61 6c 65 29 3b 0a 20 20 20 20 69 66 20 28 |scale);. if (|
00001ac0 76 3e 3d 6f 6e 65 29 20 76 3d 6f 6e 65 2d 31 3b |v>=one) v=one-1;|
00001ad0 0a 20 20 20 20 69 66 20 28 76 3c 30 29 20 20 20 |. if (v<0) |
00001ae0 20 76 3d 30 3b 0a 20 20 20 20 76 20 3d 20 70 6f | v=0;. v = po|
00001af0 77 31 36 28 76 2c 20 70 6f 77 73 63 61 6c 65 29 |w16(v, powscale)|
00001b00 3b 0a 20 20 20 20 69 66 20 28 76 3c 73 65 61 6c |;. if (v<seal|
00001b10 65 76 65 6c 29 20 76 3d 73 65 61 6c 65 76 65 6c |evel) v=sealevel|
00001b20 3b 0a 20 20 20 20 69 66 20 28 6c 69 6e 73 63 61 |;. if (linsca|
00001b30 6c 65 21 3d 6f 6e 65 29 20 76 3d 6d 75 6c 5f 66 |le!=one) v=mul_f|
00001b40 72 61 63 31 36 63 28 76 2c 20 6c 69 6e 73 63 61 |rac16c(v, linsca|
00001b50 6c 65 29 3b 0a 20 20 20 20 69 66 20 28 76 3e 3d |le);. if (v>=|
00001b60 6f 6e 65 29 20 76 3d 6f 6e 65 2d 31 3b 0a 20 20 |one) v=one-1;. |
00001b70 20 20 69 66 20 28 76 3c 30 29 20 20 20 20 76 3d | if (v<0) v=|
00001b80 30 3b 0a 20 20 20 20 2a 64 32 20 3d 20 76 3b 0a |0;. *d2 = v;.|
00001b90 20 20 7d 0a 20 20 2f 2a 4e 42 20 67 69 76 65 6e | }. /*NB given|
00001ba0 20 72 61 6e 67 65 20 6f 6e 20 64 61 74 61 2c 20 | range on data, |
00001bb0 63 6f 75 6c 64 20 70 61 63 6b 20 65 61 63 68 20 |could pack each |
00001bc0 76 61 6c 75 65 20 69 6e 74 6f 20 32 20 63 68 61 |value into 2 cha|
00001bd0 72 73 2c 20 62 75 74 20 69 73 20 73 6c 6f 77 65 |rs, but is slowe|
00001be0 72 20 74 6f 20 72 65 61 64 20 26 20 77 72 69 74 |r to read & writ|
00001bf0 65 0a 20 20 20 20 68 61 6c 66 20 77 6f 72 64 73 |e. half words|
00001c00 20 74 68 61 6e 20 77 68 6f 6c 65 20 77 6f 72 64 | than whole word|
00001c10 73 2c 20 73 6f 20 6a 75 73 74 20 72 65 64 75 63 |s, so just reduc|
00001c20 65 20 66 72 6f 6d 20 32 20 69 6e 74 73 20 64 6f |e from 2 ints do|
00001c30 77 6e 20 74 6f 20 31 20 69 6e 74 2e 2a 2f 0a 20 |wn to 1 int.*/. |
00001c40 20 72 65 74 75 72 6e 20 66 6c 65 78 5f 65 78 74 | return flex_ext|
00001c50 65 6e 64 28 28 66 6c 65 78 5f 70 74 72 29 64 61 |end((flex_ptr)da|
00001c60 74 61 2c 20 6e 2a 6e 2a 73 69 7a 65 6f 66 28 69 |ta, n*n*sizeof(i|
00001c70 6e 74 29 29 3b 0a 7d 0a 0a 42 4f 4f 4c 20 63 6c |nt));.}..BOOL cl|
00001c80 69 70 5f 66 74 28 63 6f 6d 70 6c 65 78 20 2a 64 |ip_ft(complex *d|
00001c90 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 2c 20 |ata, int limit, |
00001ca0 69 6e 74 20 72 6c 29 0a 7b 0a 20 20 69 6e 74 20 |int rl).{. int |
00001cb0 6e 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e |n=1<<limit;. in|
00001cc0 74 20 69 2c 20 6a 3b 0a 20 20 69 6e 74 20 72 68 |t i, j;. int rh|
00001cd0 3d 6e 2d 72 6c 3b 0a 0a 20 20 66 6f 72 20 28 69 |=n-rl;.. for (i|
00001ce0 3d 30 3b 20 69 3c 6e 3b 20 69 2b 2b 29 0a 20 20 |=0; i<n; i++). |
00001cf0 20 20 66 6f 72 20 28 6a 3d 30 3b 20 6a 3c 6e 3b | for (j=0; j<n;|
00001d00 20 6a 2b 2b 2c 20 64 61 74 61 2b 2b 29 0a 20 20 | j++, data++). |
00001d10 20 20 20 20 69 66 20 28 69 3e 72 6c 20 26 26 20 | if (i>rl && |
00001d20 69 3c 72 68 20 7c 7c 20 6a 3e 72 6c 20 26 26 20 |i<rh || j>rl && |
00001d30 6a 3c 72 68 29 20 64 61 74 61 2d 3e 72 65 3d 64 |j<rh) data->re=d|
00001d40 61 74 61 2d 3e 69 6d 3d 30 3b 0a 0a 20 20 72 65 |ata->im=0;.. re|
00001d50 74 75 72 6e 20 54 52 55 45 3b 0a 7d 0a 0a 42 4f |turn TRUE;.}..BO|
00001d60 4f 4c 20 73 63 72 6f 6c 6c 28 69 6e 74 20 2a 64 |OL scroll(int *d|
00001d70 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 2c 20 |ata, int limit, |
00001d80 73 63 72 6f 6c 6c 64 69 72 20 6b 65 79 2c 20 69 |scrolldir key, i|
00001d90 6e 74 20 6f 66 66 73 65 74 29 0a 7b 0a 20 20 69 |nt offset).{. i|
00001da0 6e 74 20 6e 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 |nt n=1<<limit;. |
00001db0 20 69 6e 74 20 69 2c 20 6a 3b 0a 20 20 69 6e 74 | int i, j;. int|
00001dc0 20 2a 70 2c 20 2a 70 31 2c 20 2a 70 32 3b 0a 20 | *p, *p1, *p2;. |
00001dd0 20 69 6e 74 20 62 75 66 5b 35 31 32 5d 3b 0a 20 | int buf[512];. |
00001de0 20 69 6e 74 20 73 69 3d 6f 66 66 73 65 74 3c 3c | int si=offset<<|
00001df0 28 39 2d 6c 69 6d 69 74 29 3b 0a 0a 20 20 73 77 |(9-limit);.. sw|
00001e00 69 74 63 68 20 28 6b 65 79 29 20 7b 0a 20 20 20 |itch (key) {. |
00001e10 20 63 61 73 65 20 6c 65 66 74 3a 0a 20 20 20 20 | case left:. |
00001e20 20 20 73 63 72 6f 6c 6c 78 20 2b 3d 20 6f 72 69 | scrollx += ori|
00001e30 65 6e 74 3d 3d 30 20 3f 20 2d 73 69 20 3a 20 28 |ent==0 ? -si : (|
00001e40 6f 72 69 65 6e 74 3d 3d 32 20 3f 20 73 69 20 3a |orient==2 ? si :|
00001e50 20 30 29 3b 0a 20 20 20 20 20 20 73 63 72 6f 6c | 0);. scrol|
00001e60 6c 79 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d 31 20 |ly += orient==1 |
00001e70 3f 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e 74 3d |? -si : (orient=|
00001e80 3d 33 20 3f 20 73 69 20 3a 20 30 29 3b 0a 20 20 |=3 ? si : 0);. |
00001e90 20 20 20 20 62 72 65 61 6b 3b 0a 20 20 20 20 63 | break;. c|
00001ea0 61 73 65 20 72 69 67 68 74 3a 0a 20 20 20 20 20 |ase right:. |
00001eb0 20 73 63 72 6f 6c 6c 78 20 2b 3d 20 6f 72 69 65 | scrollx += orie|
00001ec0 6e 74 3d 3d 32 20 3f 20 2d 73 69 20 3a 20 28 6f |nt==2 ? -si : (o|
00001ed0 72 69 65 6e 74 3d 3d 30 20 3f 20 73 69 20 3a 20 |rient==0 ? si : |
00001ee0 30 29 3b 0a 20 20 20 20 20 20 73 63 72 6f 6c 6c |0);. scroll|
00001ef0 79 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d 33 20 3f |y += orient==3 ?|
00001f00 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e 74 3d 3d | -si : (orient==|
00001f10 31 20 3f 20 73 69 20 3a 20 30 29 3b 0a 20 20 20 |1 ? si : 0);. |
00001f20 20 20 20 6f 66 66 73 65 74 3d 6e 2d 6f 66 66 73 | offset=n-offs|
00001f30 65 74 3b 0a 20 20 20 20 20 20 6b 65 79 3d 6c 65 |et;. key=le|
00001f40 66 74 3b 0a 20 20 20 20 20 20 62 72 65 61 6b 3b |ft;. break;|
00001f50 0a 20 20 20 20 63 61 73 65 20 64 6f 77 6e 3a 0a |. case down:.|
00001f60 20 20 20 20 20 20 73 63 72 6f 6c 6c 78 20 2b 3d | scrollx +=|
00001f70 20 6f 72 69 65 6e 74 3d 3d 33 20 3f 20 2d 73 69 | orient==3 ? -si|
00001f80 20 3a 20 28 6f 72 69 65 6e 74 3d 3d 31 20 3f 20 | : (orient==1 ? |
00001f90 73 69 20 3a 20 30 29 3b 0a 20 20 20 20 20 20 73 |si : 0);. s|
00001fa0 63 72 6f 6c 6c 79 20 2b 3d 20 6f 72 69 65 6e 74 |crolly += orient|
00001fb0 3d 3d 30 20 3f 20 2d 73 69 20 3a 20 28 6f 72 69 |==0 ? -si : (ori|
00001fc0 65 6e 74 3d 3d 32 20 3f 20 73 69 20 3a 20 30 29 |ent==2 ? si : 0)|
00001fd0 3b 0a 20 20 20 20 20 20 62 72 65 61 6b 3b 0a 20 |;. break;. |
00001fe0 20 20 20 63 61 73 65 20 75 70 3a 0a 20 20 20 20 | case up:. |
00001ff0 20 20 73 63 72 6f 6c 6c 78 20 2b 3d 20 6f 72 69 | scrollx += ori|
00002000 65 6e 74 3d 3d 31 20 3f 20 2d 73 69 20 3a 20 28 |ent==1 ? -si : (|
00002010 6f 72 69 65 6e 74 3d 3d 33 20 3f 20 73 69 20 3a |orient==3 ? si :|
00002020 20 30 29 3b 0a 20 20 20 20 20 20 73 63 72 6f 6c | 0);. scrol|
00002030 6c 79 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d 32 20 |ly += orient==2 |
00002040 3f 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e 74 3d |? -si : (orient=|
00002050 3d 30 20 3f 20 73 69 20 3a 20 30 29 3b 0a 20 20 |=0 ? si : 0);. |
00002060 20 20 20 20 6f 66 66 73 65 74 3d 6e 2d 6f 66 66 | offset=n-off|
00002070 73 65 74 3b 0a 20 20 20 20 20 20 6b 65 79 3d 64 |set;. key=d|
00002080 6f 77 6e 3b 0a 20 20 20 20 20 20 62 72 65 61 6b |own;. break|
00002090 3b 0a 20 20 7d 0a 20 20 73 63 72 6f 6c 6c 78 20 |;. }. scrollx |
000020a0 26 3d 20 35 31 31 3b 0a 20 20 73 63 72 6f 6c 6c |&= 511;. scroll|
000020b0 79 20 26 3d 20 35 31 31 3b 0a 20 20 6f 66 66 73 |y &= 511;. offs|
000020c0 65 74 20 20 26 3d 20 6e 2d 31 3b 0a 0a 20 20 69 |et &= n-1;.. i|
000020d0 66 20 28 6b 65 79 3d 3d 6c 65 66 74 29 20 7b 0a |f (key==left) {.|
000020e0 20 20 20 20 66 6f 72 20 28 69 3d 30 2c 20 70 3d | for (i=0, p=|
000020f0 64 61 74 61 3b 20 69 3c 6e 3b 20 69 2b 2b 2c 20 |data; i<n; i++, |
00002100 70 2b 3d 6e 29 20 7b 0a 20 20 20 20 20 20 66 6f |p+=n) {. fo|
00002110 72 20 28 6a 3d 6e 2d 6f 66 66 73 65 74 2c 20 70 |r (j=n-offset, p|
00002120 31 3d 70 2b 6a 2c 20 70 32 3d 62 75 66 3b 20 6a |1=p+j, p2=buf; j|
00002130 3c 6e 3b 20 6a 2b 2b 2c 20 70 31 2b 2b 2c 20 70 |<n; j++, p1++, p|
00002140 32 2b 2b 29 20 2a 70 32 3d 2a 70 31 3b 0a 20 20 |2++) *p2=*p1;. |
00002150 20 20 20 20 66 6f 72 20 28 6a 3d 6e 2d 6f 66 66 | for (j=n-off|
00002160 73 65 74 2d 31 2c 20 70 31 3d 70 2b 6a 2c 20 70 |set-1, p1=p+j, p|
00002170 32 3d 70 2b 6a 2b 6f 66 66 73 65 74 3b 20 6a 3e |2=p+j+offset; j>|
00002180 3d 30 3b 20 6a 2d 2d 2c 20 70 31 2d 2d 2c 20 70 |=0; j--, p1--, p|
00002190 32 2d 2d 29 20 2a 70 32 3d 2a 70 31 3b 0a 20 20 |2--) *p2=*p1;. |
000021a0 20 20 20 20 66 6f 72 20 28 6a 3d 30 2c 20 70 31 | for (j=0, p1|
000021b0 3d 62 75 66 2c 20 70 32 3d 70 3b 20 6a 3c 6f 66 |=buf, p2=p; j<of|
000021c0 66 73 65 74 3b 20 6a 2b 2b 2c 20 70 31 2b 2b 2c |fset; j++, p1++,|
000021d0 20 70 32 2b 2b 29 20 2a 70 32 3d 2a 70 31 3b 0a | p2++) *p2=*p1;.|
000021e0 20 20 20 20 7d 0a 20 20 7d 0a 20 20 65 6c 73 65 | }. }. else|
000021f0 20 7b 0a 20 20 20 20 66 6f 72 20 28 6a 3d 30 2c | {. for (j=0,|
00002200 20 70 3d 64 61 74 61 3b 20 6a 3c 6e 3b 20 6a 2b | p=data; j<n; j+|
00002210 2b 2c 20 70 2b 2b 29 20 7b 0a 20 20 20 20 20 20 |+, p++) {. |
00002220 66 6f 72 20 28 69 3d 6e 2d 6f 66 66 73 65 74 2c |for (i=n-offset,|
00002230 20 70 31 3d 70 2b 69 2a 6e 2c 20 70 32 3d 62 75 | p1=p+i*n, p2=bu|
00002240 66 3b 20 69 3c 6e 3b 20 69 2b 2b 2c 20 70 31 2b |f; i<n; i++, p1+|
00002250 3d 6e 2c 20 70 32 2b 2b 29 20 2a 70 32 3d 2a 70 |=n, p2++) *p2=*p|
00002260 31 3b 0a 20 20 20 20 20 20 66 6f 72 20 28 69 3d |1;. for (i=|
00002270 6e 2d 6f 66 66 73 65 74 2d 31 2c 20 70 31 3d 70 |n-offset-1, p1=p|
00002280 2b 69 2a 6e 2c 20 70 32 3d 70 2b 28 69 2b 6f 66 |+i*n, p2=p+(i+of|
00002290 66 73 65 74 29 2a 6e 3b 20 69 3e 3d 30 3b 20 69 |fset)*n; i>=0; i|
000022a0 2d 2d 2c 20 70 31 2d 3d 6e 2c 20 70 32 2d 3d 6e |--, p1-=n, p2-=n|
000022b0 29 20 2a 70 32 3d 2a 70 31 3b 0a 20 20 20 20 20 |) *p2=*p1;. |
000022c0 20 66 6f 72 20 28 69 3d 30 2c 20 70 31 3d 62 75 | for (i=0, p1=bu|
000022d0 66 2c 20 70 32 3d 70 3b 20 69 3c 6f 66 66 73 65 |f, p2=p; i<offse|
000022e0 74 3b 20 69 2b 2b 2c 20 70 31 2b 2b 2c 20 70 32 |t; i++, p1++, p2|
000022f0 2b 3d 6e 29 20 2a 70 32 3d 2a 70 31 3b 0a 20 20 |+=n) *p2=*p1;. |
00002300 20 20 7d 0a 20 20 7d 0a 0a 20 20 72 65 74 75 72 | }. }.. retur|
00002310 6e 20 54 52 55 45 3b 0a 7d 0a 0a 42 4f 4f 4c 20 |n TRUE;.}..BOOL |
00002320 72 6f 74 61 74 65 28 69 6e 74 20 2a 64 61 74 61 |rotate(int *data|
00002330 2c 20 69 6e 74 20 6c 69 6d 69 74 2c 20 42 4f 4f |, int limit, BOO|
00002340 4c 20 63 77 29 0a 7b 0a 20 20 69 6e 74 20 6e 3d |L cw).{. int n=|
00002350 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 20 |1<<limit;. int |
00002360 69 2c 20 6a 3b 0a 20 20 69 6e 74 20 6b 20 3d 20 |i, j;. int k = |
00002370 6e 2d 31 3b 0a 20 20 69 6e 74 20 6c 20 3d 20 6e |n-1;. int l = n|
00002380 2f 32 3b 0a 20 20 69 6e 74 20 2a 71 31 2c 20 2a |/2;. int *q1, *|
00002390 71 32 2c 20 2a 71 33 2c 20 2a 71 34 3b 0a 20 20 |q2, *q3, *q4;. |
000023a0 69 6e 74 20 74 3b 0a 0a 20 20 69 66 20 28 63 77 |int t;.. if (cw|
000023b0 29 20 7b 0a 20 20 20 20 6f 72 69 65 6e 74 2b 2b |) {. orient++|
000023c0 3b 0a 20 20 20 20 66 6f 72 20 28 69 3d 30 2c 20 |;. for (i=0, |
000023d0 71 31 3d 64 61 74 61 2c 20 71 32 3d 64 61 74 61 |q1=data, q2=data|
000023e0 2b 6e 2a 6b 2c 20 71 33 3d 64 61 74 61 2b 28 6e |+n*k, q3=data+(n|
000023f0 2b 31 29 2a 6b 2c 20 71 34 3d 64 61 74 61 2b 6b |+1)*k, q4=data+k|
00002400 3b 20 69 3c 6c 3b 0a 20 20 09 09 09 20 20 20 69 |; i<l;. ... i|
00002410 2b 2b 2c 20 71 31 2b 3d 6c 2c 20 71 32 2b 3d 6c |++, q1+=l, q2+=l|
00002420 2a 6e 2b 31 2c 20 71 33 2d 3d 6c 2c 20 71 34 2d |*n+1, q3-=l, q4-|
00002430 3d 6c 2a 6e 2b 31 29 0a 20 20 20 20 20 20 66 6f |=l*n+1). fo|
00002440 72 20 28 6a 3d 30 3b 20 6a 3c 6c 3b 20 6a 2b 2b |r (j=0; j<l; j++|
00002450 2c 20 71 31 2b 2b 2c 20 71 32 2d 3d 6e 2c 20 71 |, q1++, q2-=n, q|
00002460 33 2d 2d 2c 20 71 34 2b 3d 6e 29 20 7b 0a 20 20 |3--, q4+=n) {. |
00002470 20 20 20 20 20 20 74 3d 2a 71 31 3b 0a 20 20 20 | t=*q1;. |
00002480 20 20 20 20 20 2a 71 31 3d 2a 71 34 3b 0a 20 20 | *q1=*q4;. |
00002490 20 20 20 20 20 20 2a 71 34 3d 2a 71 33 3b 0a 20 | *q4=*q3;. |
000024a0 20 20 20 20 20 20 20 2a 71 33 3d 2a 71 32 3b 0a | *q3=*q2;.|
000024b0 20 20 20 20 20 20 20 20 2a 71 32 3d 74 3b 0a 20 | *q2=t;. |
000024c0 20 20 20 20 20 7d 0a 20 20 7d 0a 20 20 65 6c 73 | }. }. els|
000024d0 65 20 7b 0a 20 20 20 20 6f 72 69 65 6e 74 2d 2d |e {. orient--|
000024e0 3b 0a 20 20 20 20 66 6f 72 20 28 69 3d 30 2c 20 |;. for (i=0, |
000024f0 71 31 3d 64 61 74 61 2c 20 71 32 3d 64 61 74 61 |q1=data, q2=data|
00002500 2b 6e 2a 6b 2c 20 71 33 3d 64 61 74 61 2b 28 6e |+n*k, q3=data+(n|
00002510 2b 31 29 2a 6b 2c 20 71 34 3d 64 61 74 61 2b 6b |+1)*k, q4=data+k|
00002520 3b 20 69 3c 6c 3b 0a 20 20 09 09 09 20 20 20 69 |; i<l;. ... i|
00002530 2b 2b 2c 20 71 31 2b 3d 6c 2c 20 71 32 2b 3d 6c |++, q1+=l, q2+=l|
00002540 2a 6e 2b 31 2c 20 71 33 2d 3d 6c 2c 20 71 34 2d |*n+1, q3-=l, q4-|
00002550 3d 6c 2a 6e 2b 31 29 0a 20 20 20 20 20 20 66 6f |=l*n+1). fo|
00002560 72 20 28 6a 3d 30 3b 20 6a 3c 6c 3b 20 6a 2b 2b |r (j=0; j<l; j++|
00002570 2c 20 71 31 2b 2b 2c 20 71 32 2d 3d 6e 2c 20 71 |, q1++, q2-=n, q|
00002580 33 2d 2d 2c 20 71 34 2b 3d 6e 29 20 7b 0a 20 20 |3--, q4+=n) {. |
00002590 20 20 20 20 20 20 74 3d 2a 71 31 3b 0a 20 20 20 | t=*q1;. |
000025a0 20 20 20 20 20 2a 71 31 3d 2a 71 32 3b 0a 20 20 | *q1=*q2;. |
000025b0 20 20 20 20 20 20 2a 71 32 3d 2a 71 33 3b 0a 20 | *q2=*q3;. |
000025c0 20 20 20 20 20 20 20 2a 71 33 3d 2a 71 34 3b 0a | *q3=*q4;.|
000025d0 20 20 20 20 20 20 20 20 2a 71 34 3d 74 3b 0a 20 | *q4=t;. |
000025e0 20 20 20 20 20 7d 0a 20 20 7d 0a 20 20 6f 72 69 | }. }. ori|
000025f0 65 6e 74 20 26 3d 20 33 3b 0a 0a 20 20 72 65 74 |ent &= 3;.. ret|
00002600 75 72 6e 20 54 52 55 45 3b 0a 7d 0a 0a 2f 2a 2a |urn TRUE;.}../**|
00002610 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
*
00002630 2a 20 63 6f 6c 6f 75 72 20 72 65 6e 64 65 72 69 |* colour renderi|
00002640 6e 67 20 66 75 6e 63 74 69 6f 6e 73 20 2a 2a 2a |ng functions ***|
00002650 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
*
00002670 2a 2a 2a 2a 2a 2a 2a 2a 2a 2f 0a 0a 42 4f 4f 4c |*********/..BOOL|
00002680 20 63 72 65 61 74 65 5f 69 6d 61 67 65 33 32 28 | create_image32(|
00002690 69 6e 74 20 78 2c 20 69 6e 74 20 79 29 0a 7b 0a |int x, int y).{.|
000026a0 20 20 69 6e 74 20 61 72 65 61 73 69 7a 65 2c 20 | int areasize, |
000026b0 73 70 72 69 74 65 73 69 7a 65 2c 20 69 6d 61 67 |spritesize, imag|
000026c0 65 73 69 7a 65 3b 0a 20 20 73 70 72 69 74 65 5f |esize;. sprite_|
000026d0 68 65 61 64 65 72 20 2a 73 70 68 3b 0a 0a 20 20 |header *sph;.. |
000026e0 69 66 20 28 73 61 33 32 29 20 66 6c 65 78 5f 66 |if (sa32) flex_f|
000026f0 72 65 65 28 28 66 6c 65 78 5f 70 74 72 29 26 73 |ree((flex_ptr)&s|
00002700 61 33 32 29 3b 0a 20 20 69 6d 61 67 65 73 69 7a |a32);. imagesiz|
00002710 65 09 3d 20 78 2a 79 2a 34 3b 0a 20 20 73 70 72 |e.= x*y*4;. spr|
00002720 69 74 65 73 69 7a 65 09 3d 20 73 69 7a 65 6f 66 |itesize.= sizeof|
00002730 28 73 70 72 69 74 65 5f 68 65 61 64 65 72 29 2b |(sprite_header)+|
00002740 69 6d 61 67 65 73 69 7a 65 3b 0a 20 20 61 72 65 |imagesize;. are|
00002750 61 73 69 7a 65 09 3d 20 73 69 7a 65 6f 66 28 73 |asize.= sizeof(s|
00002760 70 72 69 74 65 5f 61 72 65 61 29 2b 73 70 72 69 |prite_area)+spri|
00002770 74 65 73 69 7a 65 3b 0a 20 20 69 66 20 28 21 66 |tesize;. if (!f|
00002780 6c 65 78 5f 61 6c 6c 6f 63 28 28 66 6c 65 78 5f |lex_alloc((flex_|
00002790 70 74 72 29 26 73 61 33 32 2c 20 61 72 65 61 73 |ptr)&sa32, areas|
000027a0 69 7a 65 29 29 20 72 65 74 75 72 6e 20 46 41 4c |ize)) return FAL|
000027b0 53 45 3b 0a 20 20 73 61 33 32 2d 3e 73 69 7a 65 |SE;. sa32->size|
000027c0 09 3d 20 61 72 65 61 73 69 7a 65 3b 0a 20 20 73 |.= areasize;. s|
000027d0 61 33 32 2d 3e 6e 75 6d 62 65 72 09 3d 20 31 3b |a32->number.= 1;|
000027e0 0a 20 20 73 61 33 32 2d 3e 73 70 72 6f 66 66 09 |. sa32->sproff.|
000027f0 3d 20 31 36 3b 0a 20 20 73 61 33 32 2d 3e 66 72 |= 16;. sa32->fr|
00002800 65 65 6f 66 66 09 3d 20 61 72 65 61 73 69 7a 65 |eeoff.= areasize|
00002810 3b 0a 20 20 73 70 68 09 09 3d 20 28 73 70 72 69 |;. sph..= (spri|
00002820 74 65 5f 68 65 61 64 65 72 20 2a 29 28 73 61 33 |te_header *)(sa3|
00002830 32 2b 31 29 3b 0a 20 20 73 70 68 2d 3e 6e 65 78 |2+1);. sph->nex|
00002840 74 09 3d 20 73 70 72 69 74 65 73 69 7a 65 3b 0a |t.= spritesize;.|
00002850 20 20 6d 65 6d 73 65 74 28 73 70 68 2d 3e 6e 61 | memset(sph->na|
00002860 6d 65 2c 20 30 2c 20 31 32 29 3b 0a 20 20 73 74 |me, 0, 12);. st|
00002870 72 63 70 79 28 73 70 68 2d 3e 6e 61 6d 65 2c 20 |rcpy(sph->name, |
00002880 69 6d 6e 61 6d 65 29 3b 0a 20 20 73 70 68 2d 3e |imname);. sph->|
00002890 77 69 64 74 68 09 3d 20 78 2d 31 3b 0a 20 20 73 |width.= x-1;. s|
000028a0 70 68 2d 3e 68 65 69 67 68 74 09 3d 20 79 2d 31 |ph->height.= y-1|
000028b0 3b 0a 20 20 73 70 68 2d 3e 6c 62 69 74 09 3d 20 |;. sph->lbit.= |
000028c0 30 3b 0a 20 20 73 70 68 2d 3e 72 62 69 74 09 3d |0;. sph->rbit.=|
000028d0 20 33 31 3b 0a 20 20 73 70 68 2d 3e 69 6d 61 67 | 31;. sph->imag|
000028e0 65 09 3d 20 73 69 7a 65 6f 66 28 73 70 72 69 74 |e.= sizeof(sprit|
000028f0 65 5f 68 65 61 64 65 72 29 3b 0a 20 20 73 70 68 |e_header);. sph|
00002900 2d 3e 6d 61 73 6b 09 3d 20 73 69 7a 65 6f 66 28 |->mask.= sizeof(|
00002910 73 70 72 69 74 65 5f 68 65 61 64 65 72 29 3b 0a |sprite_header);.|
00002920 20 20 73 70 68 2d 3e 6d 6f 64 65 09 3d 20 31 20 | sph->mode.= 1 |
00002930 7c 20 28 39 30 3c 3c 31 29 20 7c 20 28 39 30 3c || (90<<1) | (90<|
00002940 3c 31 34 29 20 7c 20 28 36 3c 3c 32 37 29 3b 0a |<14) | (6<<27);.|
00002950 20 20 6d 65 6d 73 65 74 28 28 63 68 61 72 20 2a | memset((char *|
00002960 29 73 70 68 2b 73 70 68 2d 3e 69 6d 61 67 65 2c |)sph+sph->image,|
00002970 20 30 2c 20 69 6d 61 67 65 73 69 7a 65 29 3b 0a | 0, imagesize);.|
00002980 0a 20 20 72 65 74 75 72 6e 20 54 52 55 45 3b 0a |. return TRUE;.|
00002990 7d 0a 0a 42 4f 4f 4c 20 64 69 73 63 61 72 64 5f |}..BOOL discard_|
000029a0 69 6d 61 67 65 33 32 28 76 6f 69 64 29 0a 7b 0a |image32(void).{.|
000029b0 20 20 69 66 20 28 73 61 33 32 29 20 66 6c 65 78 | if (sa32) flex|
000029c0 5f 66 72 65 65 28 28 66 6c 65 78 5f 70 74 72 29 |_free((flex_ptr)|
000029d0 26 73 61 33 32 29 3b 0a 20 20 72 65 74 75 72 6e |&sa32);. return|
000029e0 20 54 52 55 45 3b 0a 7d 0a 0a 69 6e 74 20 63 61 | TRUE;.}..int ca|
000029f0 6c 63 5f 69 6e 74 65 6e 73 28 69 6e 74 20 2a 64 |lc_intens(int *d|
00002a00 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 2c 20 |ata, int limit, |
00002a10 69 6e 74 20 69 2c 20 69 6e 74 20 6a 29 0a 7b 0a |int i, int j).{.|
00002a20 20 20 69 6e 74 20 6e 3d 31 3c 3c 6c 69 6d 69 74 | int n=1<<limit|
00002a30 3b 0a 20 20 69 6e 74 20 6e 78 2c 6e 79 2c 6e 7a |;. int nx,ny,nz|
00002a40 2c 72 2c 70 3b 0a 20 20 69 6e 74 20 7a 3b 0a 20 |,r,p;. int z;. |
00002a50 20 69 6e 74 20 69 6e 3d 30 3b 0a 0a 20 20 7a 3d | int in=0;.. z=|
00002a60 64 61 74 61 5b 69 2b 6a 2a 6e 5d 3b 0a 20 20 6e |data[i+j*n];. n|
00002a70 78 3d 7a 2d 64 61 74 61 5b 28 28 69 2b 31 29 26 |x=z-data[((i+1)&|
00002a80 28 6e 2d 31 29 29 2b 6a 2a 6e 5d 3b 0a 20 20 6e |(n-1))+j*n];. n|
00002a90 79 3d 64 61 74 61 5b 69 2b 28 28 6a 2d 31 29 26 |y=data[i+((j-1)&|
00002aa0 28 6e 2d 31 29 29 2a 6e 5d 2d 7a 3b 0a 20 20 72 |(n-1))*n]-z;. r|
00002ab0 20 3d 20 73 71 72 74 5f 66 72 61 63 31 36 28 6e | = sqrt_frac16(n|
00002ac0 78 2a 6e 78 2b 6e 79 2a 6e 79 2b 64 7a 2a 64 7a |x*nx+ny*ny+dz*dz|
00002ad0 29 3e 3e 38 3b 0a 20 20 6e 78 20 3d 20 64 69 76 |)>>8;. nx = div|
00002ae0 5f 66 72 61 63 31 36 28 6e 78 2c 20 72 29 3b 0a |_frac16(nx, r);.|
00002af0 20 20 6e 79 20 3d 20 64 69 76 5f 66 72 61 63 31 | ny = div_frac1|
00002b00 36 28 6e 79 2c 20 72 29 3b 0a 20 20 6e 7a 20 3d |6(ny, r);. nz =|
00002b10 20 64 69 76 5f 66 72 61 63 31 36 28 64 7a 2c 20 | div_frac16(dz, |
00002b20 72 29 3b 09 2f 2a 73 75 72 66 61 63 65 20 6e 6f |r);./*surface no|
00002b30 72 6d 61 6c 2a 2f 0a 0a 20 20 72 20 3d 20 6d 75 |rmal*/.. r = mu|
00002b40 6c 5f 66 72 61 63 31 36 63 28 6e 78 2c 76 78 29 |l_frac16c(nx,vx)|
00002b50 2b 6d 75 6c 5f 66 72 61 63 31 36 63 28 6e 79 2c |+mul_frac16c(ny,|
00002b60 76 79 29 2b 6d 75 6c 5f 66 72 61 63 31 36 63 28 |vy)+mul_frac16c(|
00002b70 6e 7a 2c 76 7a 29 3b 0a 20 20 69 66 20 28 72 3c |nz,vz);. if (r<|
00002b80 30 29 20 72 65 74 75 72 6e 20 69 6e 3b 0a 0a 20 |0) return in;.. |
00002b90 20 69 6e 2b 3d 28 69 6e 74 29 28 30 2e 30 35 2a | in+=(int)(0.05*|
00002ba0 6f 6e 65 29 3b 0a 0a 20 20 72 20 3d 20 6d 75 6c |one);.. r = mul|
00002bb0 5f 66 72 61 63 31 36 63 28 6e 78 2c 6c 78 29 2b |_frac16c(nx,lx)+|
00002bc0 6d 75 6c 5f 66 72 61 63 31 36 63 28 6e 79 2c 6c |mul_frac16c(ny,l|
00002bd0 79 29 2b 6d 75 6c 5f 66 72 61 63 31 36 63 28 6e |y)+mul_frac16c(n|
00002be0 7a 2c 6c 7a 29 3b 0a 20 20 69 66 20 28 72 3c 30 |z,lz);. if (r<0|
00002bf0 29 20 72 65 74 75 72 6e 20 69 6e 3b 0a 20 20 69 |) return in;. i|
00002c00 6e 20 2b 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 |n += mul_frac16c|
00002c10 28 72 2c 20 28 69 6e 74 29 28 30 2e 37 35 2a 6f |(r, (int)(0.75*o|
00002c20 6e 65 29 29 3b 0a 0a 20 20 72 20 3d 20 6d 75 6c |ne));.. r = mul|
00002c30 5f 66 72 61 63 31 36 63 28 6e 78 2c 68 78 29 2b |_frac16c(nx,hx)+|
00002c40 6d 75 6c 5f 66 72 61 63 31 36 63 28 6e 79 2c 68 |mul_frac16c(ny,h|
00002c50 79 29 2b 6d 75 6c 5f 66 72 61 63 31 36 63 28 6e |y)+mul_frac16c(n|
00002c60 7a 2c 68 7a 29 3b 0a 20 20 70 20 3d 20 6d 75 6c |z,hz);. p = mul|
00002c70 5f 66 72 61 63 31 36 63 28 72 2c 72 29 3b 0a 20 |_frac16c(r,r);. |
00002c80 20 70 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 | p = mul_frac16c|
00002c90 28 70 2c 70 29 3b 20 2f 2a 5e 34 2c 20 66 6f 72 |(p,p); /*^4, for|
00002ca0 20 5e 33 20 75 73 65 20 72 2c 70 20 69 6e 73 74 | ^3 use r,p inst|
00002cb0 65 61 64 20 6f 66 20 70 2c 70 2a 2f 0a 20 20 69 |ead of p,p*/. i|
00002cc0 6e 20 2b 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 |n += mul_frac16c|
00002cd0 28 70 2c 20 28 69 6e 74 29 28 30 2e 34 35 2a 6f |(p, (int)(0.45*o|
00002ce0 6e 65 29 29 3b 0a 0a 20 20 72 65 74 75 72 6e 20 |ne));.. return |
00002cf0 69 6e 3b 0a 7d 0a 0a 42 4f 4f 4c 20 72 65 6e 64 |in;.}..BOOL rend|
00002d00 65 72 5f 69 6d 61 67 65 33 32 28 69 6e 74 20 2a |er_image32(int *|
00002d10 64 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 2c |data, int limit,|
00002d20 20 69 6e 74 20 79 6c 2c 20 69 6e 74 20 78 6d 2c | int yl, int xm,|
00002d30 20 69 6e 74 20 79 6d 29 0a 7b 0a 20 20 69 6e 74 | int ym).{. int|
00002d40 20 6e 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 | n=1<<limit;. i|
00002d50 6e 74 20 69 2c 6a 3b 0a 20 20 69 6e 74 20 68 6f |nt i,j;. int ho|
00002d60 72 69 7a 6f 6e 3b 0a 20 20 69 6e 74 20 70 2c 70 |rizon;. int p,p|
00002d70 30 2c 70 31 3b 0a 20 20 69 6e 74 20 63 6f 73 20 |0,p1;. int cos |
00002d80 3d 20 63 6f 73 31 36 28 74 69 6c 74 29 3b 0a 20 |= cos16(tilt);. |
00002d90 20 69 6e 74 20 73 69 6e 20 3d 20 73 69 6e 31 36 | int sin = sin16|
00002da0 28 74 69 6c 74 29 3b 0a 20 20 69 6e 74 20 6c 61 |(tilt);. int la|
00002db0 73 74 69 6e 74 65 6e 73 3b 0a 20 20 69 6e 74 20 |stintens;. int |
00002dc0 6c 61 73 74 72 2c 20 6c 61 73 74 67 2c 20 6c 61 |lastr, lastg, la|
00002dd0 73 74 62 2c 20 6e 65 77 72 2c 20 6e 65 77 67 2c |stb, newr, newg,|
00002de0 20 6e 65 77 62 3b 0a 20 20 42 4f 4f 4c 20 6c 61 | newb;. BOOL la|
00002df0 73 74 76 61 6c 69 64 3b 0a 20 20 73 70 72 69 74 |stvalid;. sprit|
00002e00 65 5f 68 65 61 64 65 72 20 2a 73 70 68 20 3d 20 |e_header *sph = |
00002e10 28 73 70 72 69 74 65 5f 68 65 61 64 65 72 20 2a |(sprite_header *|
00002e20 29 28 28 63 68 61 72 20 2a 29 73 61 33 32 2b 73 |)((char *)sa32+s|
00002e30 61 33 32 2d 3e 73 70 72 6f 66 66 29 3b 0a 20 20 |a32->sproff);. |
00002e40 69 6e 74 20 2a 69 6e 74 65 6e 73 20 3d 20 28 69 |int *intens = (i|
00002e50 6e 74 20 2a 29 28 28 63 68 61 72 20 2a 29 73 70 |nt *)((char *)sp|
00002e60 68 2b 73 70 68 2d 3e 69 6d 61 67 65 29 3b 0a 20 |h+sph->image);. |
00002e70 20 69 6e 74 20 69 6e 3b 0a 20 20 69 6e 74 20 64 | int in;. int d|
00002e80 65 6c 74 61 72 2c 20 64 65 6c 74 61 67 2c 20 64 |eltar, deltag, d|
00002e90 65 6c 74 61 62 3b 0a 20 20 69 6e 74 20 72 3b 0a |eltab;. int r;.|
00002ea0 20 20 69 6e 74 20 69 6e 63 6c 2c 20 62 65 61 72 | int incl, bear|
00002eb0 3b 0a 20 20 75 6e 73 69 67 6e 65 64 20 69 6e 74 |;. unsigned int|
00002ec0 20 63 6f 6c 6f 75 72 3b 0a 20 20 69 6e 74 20 63 | colour;. int c|
00002ed0 72 2c 63 67 2c 63 62 3b 0a 0a 20 20 64 7a 20 3d |r,cg,cb;.. dz =|
00002ee0 20 36 35 35 33 36 2f 28 6e 2d 31 29 3b 0a 20 20 | 65536/(n-1);. |
00002ef0 76 78 20 3d 20 30 3b 0a 20 20 76 79 20 3d 20 2d |vx = 0;. vy = -|
00002f00 63 6f 73 3b 0a 20 20 76 7a 20 3d 20 73 69 6e 3b |cos;. vz = sin;|
00002f10 0a 20 20 69 6e 63 6c 20 3d 20 28 69 6e 74 29 28 |. incl = (int)(|
00002f20 30 2e 35 30 2a 6f 6e 65 29 3b 0a 20 20 62 65 61 |0.50*one);. bea|
00002f30 72 20 3d 20 28 69 6e 74 29 28 31 2e 35 30 2a 6f |r = (int)(1.50*o|
00002f40 6e 65 29 3b 0a 20 20 72 20 3d 20 63 6f 73 31 36 |ne);. r = cos16|
00002f50 28 69 6e 63 6c 29 3b 0a 20 20 6c 78 20 3d 20 6d |(incl);. lx = m|
00002f60 75 6c 5f 66 72 61 63 31 36 63 28 63 6f 73 31 36 |ul_frac16c(cos16|
00002f70 28 62 65 61 72 29 2c 20 72 29 3b 0a 20 20 6c 79 |(bear), r);. ly|
00002f80 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 28 73 | = mul_frac16c(s|
00002f90 69 6e 31 36 28 62 65 61 72 29 2c 2d 72 29 3b 0a |in16(bear),-r);.|
00002fa0 20 20 6c 7a 20 3d 20 73 69 6e 31 36 28 69 6e 63 | lz = sin16(inc|
00002fb0 6c 29 3b 0a 20 20 68 78 20 3d 20 28 6c 78 2b 76 |l);. hx = (lx+v|
00002fc0 78 29 2f 32 3b 0a 20 20 68 79 20 3d 20 28 6c 79 |x)/2;. hy = (ly|
00002fd0 2b 76 79 29 2f 32 3b 0a 20 20 68 7a 20 3d 20 28 |+vy)/2;. hz = (|
00002fe0 6c 7a 2b 76 7a 29 2f 32 3b 0a 20 20 72 20 3d 20 |lz+vz)/2;. r = |
00002ff0 73 71 72 74 5f 66 72 61 63 31 36 28 6d 75 6c 5f |sqrt_frac16(mul_|
00003000 66 72 61 63 31 36 63 28 68 78 2c 68 78 29 2b 6d |frac16c(hx,hx)+m|
00003010 75 6c 5f 66 72 61 63 31 36 63 28 68 79 2c 68 79 |ul_frac16c(hy,hy|
00003020 29 2b 6d 75 6c 5f 66 72 61 63 31 36 63 28 68 7a |)+mul_frac16c(hz|
00003030 2c 68 7a 29 29 3b 0a 20 20 68 78 20 3d 20 64 69 |,hz));. hx = di|
00003040 76 5f 66 72 61 63 31 36 28 68 78 2c 20 72 29 3b |v_frac16(hx, r);|
00003050 0a 20 20 68 79 20 3d 20 64 69 76 5f 66 72 61 63 |. hy = div_frac|
00003060 31 36 28 68 79 2c 20 72 29 3b 0a 20 20 68 7a 20 |16(hy, r);. hz |
00003070 3d 20 64 69 76 5f 66 72 61 63 31 36 28 68 7a 2c |= div_frac16(hz,|
00003080 20 72 29 3b 0a 0a 20 20 66 6f 72 20 28 69 3d 30 | r);.. for (i=0|
00003090 3b 20 69 3c 6e 3b 20 69 2b 2b 29 20 7b 0a 20 20 |; i<n; i++) {. |
000030a0 20 20 70 30 20 3d 20 28 20 28 6d 75 6c 5f 66 72 | p0 = ( (mul_fr|
000030b0 61 63 31 36 28 64 61 74 61 5b 69 5d 2c 20 63 6f |ac16(data[i], co|
000030c0 73 29 2b 6d 75 6c 5f 66 72 61 63 31 36 28 30 3c |s)+mul_frac16(0<|
000030d0 3c 28 31 36 2d 6c 69 6d 69 74 29 2c 20 73 69 6e |<(16-limit), sin|
000030e0 29 29 2a 33 2f 32 35 36 20 2d 20 79 6c 29 20 3e |))*3/256 - yl) >|
000030f0 3e 20 28 31 30 2d 6c 69 6d 69 74 29 3b 0a 20 20 |> (10-limit);. |
00003100 20 20 63 6f 6c 6f 75 72 20 3d 20 65 61 72 74 68 | colour = earth|
00003110 5b 64 61 74 61 5b 69 5d 2f 32 35 36 5d 3b 0a 20 |[data[i]/256];. |
00003120 20 20 20 63 72 20 3d 20 63 6f 6c 6f 75 72 20 26 | cr = colour &|
00003130 20 32 35 35 3b 0a 20 20 20 20 63 67 20 3d 20 28 | 255;. cg = (|
00003140 63 6f 6c 6f 75 72 3e 3e 38 29 20 26 20 32 35 35 |colour>>8) & 255|
00003150 3b 0a 20 20 20 20 63 62 20 3d 20 28 63 6f 6c 6f |;. cb = (colo|
00003160 75 72 3e 3e 31 36 29 20 26 20 32 35 35 3b 0a 20 |ur>>16) & 255;. |
00003170 20 20 20 69 6e 20 3d 20 63 61 6c 63 5f 69 6e 74 | in = calc_int|
00003180 65 6e 73 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c |ens(data, limit,|
00003190 20 69 2c 20 30 29 3b 0a 20 20 20 20 63 72 20 3d | i, 0);. cr =|
000031a0 20 28 63 72 2a 69 6e 29 2f 32 35 36 3b 0a 20 20 | (cr*in)/256;. |
000031b0 20 20 63 67 20 3d 20 28 63 67 2a 69 6e 29 2f 32 | cg = (cg*in)/2|
000031c0 35 36 3b 0a 20 20 20 20 63 62 20 3d 20 28 63 62 |56;. cb = (cb|
000031d0 2a 69 6e 29 2f 32 35 36 3b 0a 20 20 20 20 69 66 |*in)/256;. if|
000031e0 20 28 63 72 3c 30 29 20 63 72 3d 30 3b 0a 20 20 | (cr<0) cr=0;. |
000031f0 20 20 69 66 20 28 63 67 3c 30 29 20 63 67 3d 30 | if (cg<0) cg=0|
00003200 3b 0a 20 20 20 20 69 66 20 28 63 62 3c 30 29 20 |;. if (cb<0) |
00003210 63 62 3d 30 3b 0a 20 20 20 20 69 66 20 28 63 72 |cb=0;. if (cr|
00003220 3e 36 35 35 33 35 29 20 63 72 3d 36 35 35 33 35 |>65535) cr=65535|
00003230 3b 0a 20 20 20 20 69 66 20 28 63 67 3e 36 35 35 |;. if (cg>655|
00003240 33 35 29 20 63 67 3d 36 35 35 33 35 3b 0a 20 20 |35) cg=65535;. |
00003250 20 20 69 66 20 28 63 62 3e 36 35 35 33 35 29 20 | if (cb>65535) |
00003260 63 62 3d 36 35 35 33 35 3b 0a 20 20 20 20 6c 61 |cb=65535;. la|
00003270 73 74 72 20 3d 20 63 72 3b 0a 20 20 20 20 6c 61 |str = cr;. la|
00003280 73 74 67 20 3d 20 63 67 3b 0a 20 20 20 20 6c 61 |stg = cg;. la|
00003290 73 74 62 20 3d 20 63 62 3b 0a 20 20 20 20 6c 61 |stb = cb;. la|
000032a0 73 74 76 61 6c 69 64 3d 54 52 55 45 3b 0a 20 20 |stvalid=TRUE;. |
000032b0 20 20 69 66 20 28 70 30 3e 3d 30 20 26 26 20 70 | if (p0>=0 && p|
000032c0 30 3c 79 6d 29 20 69 6e 74 65 6e 73 5b 69 2b 78 |0<ym) intens[i+x|
000032d0 6d 2a 28 79 6d 2d 31 2d 70 30 29 5d 20 3d 20 28 |m*(ym-1-p0)] = (|
000032e0 63 72 3e 3e 38 29 2b 28 63 67 26 28 32 35 35 3c |cr>>8)+(cg&(255<|
000032f0 3c 38 29 29 2b 28 28 63 62 3e 3e 38 29 3c 3c 31 |<8))+((cb>>8)<<1|
00003300 36 29 3b 0a 20 20 20 20 68 6f 72 69 7a 6f 6e 20 |6);. horizon |
00003310 3d 20 70 30 3b 0a 20 20 20 20 66 6f 72 20 28 6a |= p0;. for (j|
00003320 3d 31 3b 20 6a 3c 6e 3b 20 6a 2b 2b 29 20 7b 0a |=1; j<n; j++) {.|
00003330 20 20 20 20 20 20 70 31 20 3d 20 28 20 28 6d 75 | p1 = ( (mu|
00003340 6c 5f 66 72 61 63 31 36 28 64 61 74 61 5b 69 2b |l_frac16(data[i+|
00003350 6a 2a 6e 5d 2c 20 63 6f 73 29 2b 6d 75 6c 5f 66 |j*n], cos)+mul_f|
00003360 72 61 63 31 36 28 6a 3c 3c 28 31 36 2d 6c 69 6d |rac16(j<<(16-lim|
00003370 69 74 29 2c 20 73 69 6e 29 29 2a 33 2f 32 35 36 |it), sin))*3/256|
00003380 20 2d 20 79 6c 29 20 3e 3e 20 28 31 30 2d 6c 69 | - yl) >> (10-li|
00003390 6d 69 74 29 3b 0a 20 20 20 20 20 20 69 66 20 28 |mit);. if (|
000033a0 70 31 3e 68 6f 72 69 7a 6f 6e 29 20 7b 0a 20 20 |p1>horizon) {. |
000033b0 20 20 20 20 20 20 63 6f 6c 6f 75 72 20 3d 20 65 | colour = e|
000033c0 61 72 74 68 5b 64 61 74 61 5b 69 2b 6a 2a 6e 5d |arth[data[i+j*n]|
000033d0 2f 32 35 36 5d 3b 0a 20 20 20 20 20 20 20 20 63 |/256];. c|
000033e0 72 20 3d 20 63 6f 6c 6f 75 72 20 26 20 32 35 35 |r = colour & 255|
000033f0 3b 0a 20 20 20 20 20 20 20 20 63 67 20 3d 20 28 |;. cg = (|
00003400 63 6f 6c 6f 75 72 3e 3e 38 29 20 26 20 32 35 35 |colour>>8) & 255|
00003410 3b 0a 20 20 20 20 20 20 20 20 63 62 20 3d 20 28 |;. cb = (|
00003420 63 6f 6c 6f 75 72 3e 3e 31 36 29 20 26 20 32 35 |colour>>16) & 25|
00003430 35 3b 0a 20 20 20 20 20 20 20 20 69 6e 20 3d 20 |5;. in = |
00003440 63 61 6c 63 5f 69 6e 74 65 6e 73 28 64 61 74 61 |calc_intens(data|
00003450 2c 20 6c 69 6d 69 74 2c 20 69 2c 20 6a 29 3b 0a |, limit, i, j);.|
00003460 20 20 20 20 20 20 20 20 63 72 20 3d 20 28 63 72 | cr = (cr|
00003470 2a 69 6e 29 2f 32 35 36 3b 0a 20 20 20 20 20 20 |*in)/256;. |
00003480 20 20 63 67 20 3d 20 28 63 67 2a 69 6e 29 2f 32 | cg = (cg*in)/2|
00003490 35 36 3b 0a 20 20 20 20 20 20 20 20 63 62 20 3d |56;. cb =|
000034a0 20 28 63 62 2a 69 6e 29 2f 32 35 36 3b 0a 20 20 | (cb*in)/256;. |
000034b0 20 20 20 20 20 20 69 66 20 28 63 72 3c 30 29 20 | if (cr<0) |
000034c0 63 72 3d 30 3b 0a 20 20 20 20 20 20 20 20 69 66 |cr=0;. if|
000034d0 20 28 63 67 3c 30 29 20 63 67 3d 30 3b 0a 20 20 | (cg<0) cg=0;. |
000034e0 20 20 20 20 20 20 69 66 20 28 63 62 3c 30 29 20 | if (cb<0) |
000034f0 63 62 3d 30 3b 0a 20 20 20 20 20 20 20 20 69 66 |cb=0;. if|
00003500 20 28 63 72 3e 36 35 35 33 35 29 20 63 72 3d 36 | (cr>65535) cr=6|
00003510 35 35 33 35 3b 0a 20 20 20 20 20 20 20 20 69 66 |5535;. if|
00003520 20 28 63 67 3e 36 35 35 33 35 29 20 63 67 3d 36 | (cg>65535) cg=6|
00003530 35 35 33 35 3b 0a 20 20 20 20 20 20 20 20 69 66 |5535;. if|
00003540 20 28 63 62 3e 36 35 35 33 35 29 20 63 62 3d 36 | (cb>65535) cb=6|
00003550 35 35 33 35 3b 0a 20 20 20 20 20 20 20 20 6e 65 |5535;. ne|
00003560 77 72 20 3d 20 63 72 3b 0a 20 20 20 20 20 20 20 |wr = cr;. |
00003570 20 6e 65 77 67 20 3d 20 63 67 3b 0a 20 20 20 20 | newg = cg;. |
00003580 20 20 20 20 6e 65 77 62 20 3d 20 63 62 3b 0a 20 | newb = cb;. |
00003590 20 20 20 20 20 20 20 69 66 20 28 70 31 3e 3d 30 | if (p1>=0|
000035a0 20 26 26 20 70 31 3c 79 6d 29 20 69 6e 74 65 6e | && p1<ym) inten|
000035b0 73 5b 69 2b 78 6d 2a 28 79 6d 2d 31 2d 70 31 29 |s[i+xm*(ym-1-p1)|
000035c0 5d 20 3d 20 28 63 72 3e 3e 38 29 2b 28 63 67 26 |] = (cr>>8)+(cg&|
000035d0 28 32 35 35 3c 3c 38 29 29 2b 28 28 63 62 3e 3e |(255<<8))+((cb>>|
000035e0 38 29 3c 3c 31 36 29 3b 0a 20 20 20 20 20 20 20 |8)<<16);. |
000035f0 20 70 20 3d 20 70 31 2d 31 3b 0a 20 20 20 20 20 | p = p1-1;. |
00003600 20 20 20 69 66 20 28 70 3e 68 6f 72 69 7a 6f 6e | if (p>horizon|
00003610 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 69 66 |) {. if|
00003620 20 28 21 6c 61 73 74 76 61 6c 69 64 29 20 7b 0a | (!lastvalid) {.|
00003630 20 20 20 20 20 20 20 20 20 20 20 20 63 6f 6c 6f | colo|
00003640 75 72 20 3d 20 65 61 72 74 68 5b 64 61 74 61 5b |ur = earth[data[|
00003650 69 2b 28 6a 2d 31 29 2a 6e 5d 2f 32 35 36 5d 3b |i+(j-1)*n]/256];|
00003660 0a 20 20 20 20 20 20 20 20 20 20 20 20 6c 61 73 |. las|
00003670 74 72 20 3d 20 63 6f 6c 6f 75 72 20 26 20 32 35 |tr = colour & 25|
00003680 35 3b 0a 20 20 20 20 20 20 20 20 20 20 20 20 6c |5;. l|
00003690 61 73 74 67 20 3d 20 28 63 6f 6c 6f 75 72 3e 3e |astg = (colour>>|
000036a0 38 29 20 26 20 32 35 35 3b 0a 20 20 20 20 20 20 |8) & 255;. |
000036b0 20 20 20 20 20 20 6c 61 73 74 62 20 3d 20 28 63 | lastb = (c|
000036c0 6f 6c 6f 75 72 3e 3e 31 36 29 20 26 20 32 35 35 |olour>>16) & 255|
000036d0 3b 0a 20 20 20 20 20 20 20 20 20 20 20 20 6c 61 |;. la|
000036e0 73 74 69 6e 74 65 6e 73 20 3d 20 63 61 6c 63 5f |stintens = calc_|
000036f0 69 6e 74 65 6e 73 28 64 61 74 61 2c 20 6c 69 6d |intens(data, lim|
00003700 69 74 2c 20 69 2c 20 6a 2d 31 29 3b 0a 20 20 20 |it, i, j-1);. |
00003710 20 20 20 20 20 20 20 20 20 6c 61 73 74 72 20 3d | lastr =|
00003720 20 28 6c 61 73 74 72 2a 6c 61 73 74 69 6e 74 65 | (lastr*lastinte|
00003730 6e 73 29 2f 32 35 36 3b 0a 20 20 20 20 20 20 20 |ns)/256;. |
00003740 20 20 20 20 20 6c 61 73 74 67 20 3d 20 28 6c 61 | lastg = (la|
00003750 73 74 67 2a 6c 61 73 74 69 6e 74 65 6e 73 29 2f |stg*lastintens)/|
00003760 32 35 36 3b 0a 20 20 20 20 20 20 20 20 20 20 20 |256;. |
00003770 20 6c 61 73 74 62 20 3d 20 28 6c 61 73 74 62 2a | lastb = (lastb*|
00003780 6c 61 73 74 69 6e 74 65 6e 73 29 2f 32 35 36 3b |lastintens)/256;|
00003790 0a 20 20 20 20 20 20 20 20 20 20 20 20 69 66 20 |. if |
000037a0 28 6c 61 73 74 72 3c 30 29 20 6c 61 73 74 72 3d |(lastr<0) lastr=|
000037b0 30 3b 0a 20 20 20 20 20 20 20 20 20 20 20 20 69 |0;. i|
000037c0 66 20 28 6c 61 73 74 67 3c 30 29 20 6c 61 73 74 |f (lastg<0) last|
000037d0 67 3d 30 3b 0a 20 20 20 20 20 20 20 20 20 20 20 |g=0;. |
000037e0 20 69 66 20 28 6c 61 73 74 62 3c 30 29 20 6c 61 | if (lastb<0) la|
000037f0 73 74 62 3d 30 3b 0a 20 20 20 20 20 20 20 20 20 |stb=0;. |
00003800 20 20 20 69 66 20 28 6c 61 73 74 72 3e 36 35 35 | if (lastr>655|
00003810 33 35 29 20 6c 61 73 74 72 3d 36 35 35 33 35 3b |35) lastr=65535;|
00003820 0a 20 20 20 20 20 20 20 20 20 20 20 20 69 66 20 |. if |
00003830 28 6c 61 73 74 67 3e 36 35 35 33 35 29 20 6c 61 |(lastg>65535) la|
00003840 73 74 67 3d 36 35 35 33 35 3b 0a 20 20 20 20 20 |stg=65535;. |
00003850 20 20 20 20 20 20 20 69 66 20 28 6c 61 73 74 62 | if (lastb|
00003860 3e 36 35 35 33 35 29 20 6c 61 73 74 62 3d 36 35 |>65535) lastb=65|
00003870 35 33 35 3b 0a 20 20 20 20 20 20 20 20 20 20 7d |535;. }|
00003880 0a 20 20 20 20 20 20 20 20 20 20 64 65 6c 74 61 |. delta|
00003890 72 20 3d 20 28 6c 61 73 74 72 2d 6e 65 77 72 29 |r = (lastr-newr)|
000038a0 2f 28 70 30 2d 70 31 29 3b 0a 20 20 20 20 20 20 |/(p0-p1);. |
000038b0 20 20 20 20 64 65 6c 74 61 67 20 3d 20 28 6c 61 | deltag = (la|
000038c0 73 74 67 2d 6e 65 77 67 29 2f 28 70 30 2d 70 31 |stg-newg)/(p0-p1|
000038d0 29 3b 0a 20 20 20 20 20 20 20 20 20 20 64 65 6c |);. del|
000038e0 74 61 62 20 3d 20 28 6c 61 73 74 62 2d 6e 65 77 |tab = (lastb-new|
000038f0 62 29 2f 28 70 30 2d 70 31 29 3b 0a 20 20 20 20 |b)/(p0-p1);. |
00003900 20 20 20 20 20 20 77 68 69 6c 65 20 28 70 3e 68 | while (p>h|
00003910 6f 72 69 7a 6f 6e 29 20 7b 0a 20 20 20 20 20 20 |orizon) {. |
00003920 20 20 20 20 20 20 63 72 2d 3d 64 65 6c 74 61 72 | cr-=deltar|
00003930 3b 0a 20 20 20 20 20 20 20 20 20 20 20 20 63 67 |;. cg|
00003940 2d 3d 64 65 6c 74 61 67 3b 0a 20 20 20 20 20 20 |-=deltag;. |
00003950 20 20 20 20 20 20 63 62 2d 3d 64 65 6c 74 61 62 | cb-=deltab|
00003960 3b 0a 20 20 20 20 20 20 20 20 20 20 20 20 69 66 |;. if|
00003970 20 28 70 3e 3d 30 20 26 26 20 70 3c 79 6d 29 20 | (p>=0 && p<ym) |
00003980 69 6e 74 65 6e 73 5b 69 2b 78 6d 2a 28 79 6d 2d |intens[i+xm*(ym-|
00003990 31 2d 70 29 5d 20 3d 20 28 63 72 3e 3e 38 29 2b |1-p)] = (cr>>8)+|
000039a0 28 63 67 26 28 32 35 35 3c 3c 38 29 29 2b 28 28 |(cg&(255<<8))+((|
000039b0 63 62 3e 3e 38 29 3c 3c 31 36 29 3b 0a 20 20 20 |cb>>8)<<16);. |
000039c0 20 20 20 20 20 20 20 20 20 70 2d 3d 31 3b 0a 20 | p-=1;. |
000039d0 20 20 20 20 20 20 20 20 20 7d 0a 20 20 20 20 20 | }. |
000039e0 20 20 20 7d 0a 20 20 20 20 20 20 20 20 68 6f 72 | }. hor|
000039f0 69 7a 6f 6e 20 3d 20 70 31 3b 0a 20 20 20 20 20 |izon = p1;. |
00003a00 20 20 20 6c 61 73 74 72 3d 6e 65 77 72 3b 0a 20 | lastr=newr;. |
00003a10 20 20 20 20 20 20 20 6c 61 73 74 67 3d 6e 65 77 | lastg=new|
00003a20 67 3b 0a 20 20 20 20 20 20 20 20 6c 61 73 74 62 |g;. lastb|
00003a30 3d 6e 65 77 62 3b 0a 20 20 20 20 20 20 20 20 6c |=newb;. l|
00003a40 61 73 74 76 61 6c 69 64 20 3d 20 54 52 55 45 3b |astvalid = TRUE;|
00003a50 0a 20 20 20 20 20 20 7d 0a 20 20 20 20 20 20 65 |. }. e|
00003a60 6c 73 65 20 7b 0a 20 20 20 20 20 20 20 20 6c 61 |lse {. la|
00003a70 73 74 76 61 6c 69 64 3d 46 41 4c 53 45 3b 0a 20 |stvalid=FALSE;. |
00003a80 20 20 20 20 20 7d 0a 20 20 20 20 20 20 70 30 3d | }. p0=|
00003a90 70 31 3b 0a 20 20 20 20 7d 0a 20 20 7d 0a 0a 20 |p1;. }. }.. |
00003aa0 20 72 65 74 75 72 6e 20 54 52 55 45 3b 0a 7d 0a | return TRUE;.}.|
00003ab0 0a 42 4f 4f 4c 20 64 75 6d 70 5f 69 6d 61 67 65 |.BOOL dump_image|
00003ac0 33 32 28 69 6e 74 20 2a 64 61 74 61 2c 20 69 6e |32(int *data, in|
00003ad0 74 20 6c 69 6d 69 74 2c 20 69 6e 74 20 79 6c 2c |t limit, int yl,|
00003ae0 20 69 6e 74 20 68 65 69 67 68 74 2c 20 63 68 61 | int height, cha|
00003af0 72 20 2a 66 69 6c 65 6e 61 6d 65 29 0a 7b 0a 20 |r *filename).{. |
00003b00 20 69 6e 74 20 6e 3d 31 3c 3c 6c 69 6d 69 74 3b | int n=1<<limit;|
00003b10 0a 20 20 42 4f 4f 4c 20 72 65 74 63 6f 64 65 3d |. BOOL retcode=|
00003b20 54 52 55 45 3b 0a 20 20 69 6e 74 20 78 3d 6e 3b |TRUE;. int x=n;|
00003b30 0a 20 20 69 6e 74 20 79 3d 68 65 69 67 68 74 3b |. int y=height;|
00003b40 0a 20 20 63 68 61 72 20 62 75 66 5b 32 35 36 5d |. char buf[256]|
00003b50 3b 0a 0a 20 20 69 66 20 28 21 63 72 65 61 74 65 |;.. if (!create|
00003b60 5f 69 6d 61 67 65 33 32 28 78 2c 20 79 29 29 20 |_image32(x, y)) |
00003b70 72 65 74 75 72 6e 20 46 41 4c 53 45 3b 0a 20 20 |return FALSE;. |
00003b80 69 66 20 28 21 72 65 6e 64 65 72 5f 69 6d 61 67 |if (!render_imag|
00003b90 65 33 32 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c |e32(data, limit,|
00003ba0 20 79 6c 2c 20 78 2c 20 79 29 29 20 7b 0a 20 20 | yl, x, y)) {. |
00003bb0 20 20 64 69 73 63 61 72 64 5f 69 6d 61 67 65 33 | discard_image3|
00003bc0 32 28 29 3b 0a 20 20 20 20 72 65 74 75 72 6e 20 |2();. return |
00003bd0 46 41 4c 53 45 3b 0a 20 20 7d 0a 20 20 69 66 20 |FALSE;. }. if |
00003be0 28 73 61 76 65 33 32 29 20 7b 0a 20 20 20 20 69 |(save32) {. i|
00003bf0 66 20 28 73 70 72 69 74 65 5f 61 72 65 61 5f 73 |f (sprite_area_s|
00003c00 61 76 65 28 73 61 33 32 2c 20 66 69 6c 65 6e 61 |ave(sa32, filena|
00003c10 6d 65 29 29 20 72 65 74 63 6f 64 65 3d 46 41 4c |me)) retcode=FAL|
00003c20 53 45 3b 0a 20 20 7d 0a 20 20 69 66 20 28 73 61 |SE;. }. if (sa|
00003c30 76 65 38 29 20 7b 0a 20 20 20 20 69 66 20 28 64 |ve8) {. if (d|
00003c40 69 74 68 61 72 65 61 29 20 66 6c 65 78 5f 66 72 |itharea) flex_fr|
00003c50 65 65 28 28 66 6c 65 78 5f 70 74 72 29 26 64 69 |ee((flex_ptr)&di|
00003c60 74 68 61 72 65 61 29 3b 0a 20 20 20 20 69 66 20 |tharea);. if |
00003c70 28 21 66 6c 65 78 5f 61 6c 6c 6f 63 28 28 66 6c |(!flex_alloc((fl|
00003c80 65 78 5f 70 74 72 29 26 64 69 74 68 61 72 65 61 |ex_ptr)&ditharea|
00003c90 2c 20 73 69 7a 65 6f 66 28 73 70 72 69 74 65 5f |, sizeof(sprite_|
00003ca0 61 72 65 61 29 29 29 20 7b 0a 20 20 20 20 20 20 |area))) {. |
00003cb0 77 65 72 72 28 30 2c 20 22 64 69 74 68 61 72 65 |werr(0, "dithare|
00003cc0 61 20 69 6e 69 74 20 66 61 69 6c 65 64 22 29 3b |a init failed");|
00003cd0 0a 20 20 20 20 20 20 72 65 74 63 6f 64 65 3d 46 |. retcode=F|
00003ce0 41 4c 53 45 3b 0a 20 20 20 20 7d 0a 20 20 20 20 |ALSE;. }. |
00003cf0 65 6c 73 65 20 7b 0a 20 20 20 20 20 20 73 70 72 |else {. spr|
00003d00 69 74 65 5f 61 72 65 61 5f 69 6e 69 74 69 61 6c |ite_area_initial|
00003d10 69 73 65 28 64 69 74 68 61 72 65 61 2c 20 66 6c |ise(ditharea, fl|
00003d20 65 78 5f 73 69 7a 65 28 28 66 6c 65 78 5f 70 74 |ex_size((flex_pt|
00003d30 72 29 26 64 69 74 68 61 72 65 61 29 29 3b 0a 20 |r)&ditharea));. |
00003d40 20 20 20 20 20 69 66 20 28 67 65 6e 5f 32 35 36 | if (gen_256|
00003d50 77 69 74 68 6e 6f 70 61 6c 65 74 74 65 28 28 66 |withnopalette((f|
00003d60 6c 65 78 5f 70 74 72 29 26 73 61 33 32 2c 20 69 |lex_ptr)&sa32, i|
00003d70 6d 6e 61 6d 65 29 29 20 7b 0a 20 20 20 20 20 20 |mname)) {. |
00003d80 20 20 73 74 72 63 61 74 28 73 74 72 63 70 79 28 | strcat(strcpy(|
00003d90 62 75 66 2c 20 66 69 6c 65 6e 61 6d 65 29 2c 20 |buf, filename), |
00003da0 22 64 22 29 3b 0a 20 20 20 20 20 20 20 20 69 66 |"d");. if|
00003db0 20 28 73 70 72 69 74 65 5f 61 72 65 61 5f 73 61 | (sprite_area_sa|
00003dc0 76 65 28 64 69 74 68 61 72 65 61 2c 20 62 75 66 |ve(ditharea, buf|
00003dd0 29 29 20 72 65 74 63 6f 64 65 3d 46 41 4c 53 45 |)) retcode=FALSE|
00003de0 3b 0a 20 20 20 20 20 20 7d 0a 20 20 20 20 20 20 |;. }. |
00003df0 65 6c 73 65 20 72 65 74 63 6f 64 65 20 3d 20 46 |else retcode = F|
00003e00 41 4c 53 45 3b 0a 20 20 20 20 20 20 66 6c 65 78 |ALSE;. flex|
00003e10 5f 66 72 65 65 28 28 66 6c 65 78 5f 70 74 72 29 |_free((flex_ptr)|
00003e20 26 64 69 74 68 61 72 65 61 29 3b 0a 20 20 20 20 |&ditharea);. |
00003e30 7d 0a 20 20 7d 0a 20 20 64 69 73 63 61 72 64 5f |}. }. discard_|
00003e40 69 6d 61 67 65 33 32 28 29 3b 0a 0a 20 20 72 65 |image32();.. re|
00003e50 74 75 72 6e 20 72 65 74 63 6f 64 65 3b 0a 7d 0a |turn retcode;.}.|
00003e60 0a 2f 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |./**************|
00003e70 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00003e80 2a 2a 20 64 69 74 68 65 72 69 6e 67 20 2a 2a 2a |** dithering ***|
00003e90 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00003ea0 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2f 0a 0a |*************/..|
00003eb0 42 4f 4f 4c 20 67 65 6e 5f 32 35 36 77 69 74 68 |BOOL gen_256with|
00003ec0 6e 6f 70 61 6c 65 74 74 65 28 66 6c 65 78 5f 70 |nopalette(flex_p|
00003ed0 74 72 20 73 61 2c 20 63 68 61 72 20 2a 6e 61 6d |tr sa, char *nam|
00003ee0 65 29 0a 7b 0a 20 20 69 6e 74 20 6e 65 77 73 69 |e).{. int newsi|
00003ef0 7a 65 3b 0a 20 20 6f 73 5f 65 72 72 6f 72 20 2a |ze;. os_error *|
00003f00 65 3b 0a 20 20 63 6f 6d 70 61 63 74 5f 64 69 74 |e;. compact_dit|
00003f10 68 65 72 5f 65 72 72 6f 72 20 2a 64 65 3b 0a 20 |her_error *de;. |
00003f20 20 69 6e 74 20 78 2c 20 79 3b 0a 20 20 73 70 72 | int x, y;. spr|
00003f30 69 74 65 5f 69 64 20 73 69 3b 0a 20 20 73 70 72 |ite_id si;. spr|
00003f40 69 74 65 5f 68 65 61 64 65 72 20 2a 73 68 2c 20 |ite_header *sh, |
00003f50 2a 64 68 3b 0a 20 20 63 68 61 72 20 2a 73 70 2c |*dh;. char *sp,|
00003f60 20 2a 64 70 3b 0a 20 20 75 6e 73 69 67 6e 65 64 | *dp;. unsigned|
00003f70 20 69 6e 74 20 6d 6f 64 65 77 6f 72 64 2c 20 74 | int modeword, t|
00003f80 79 70 65 2c 20 78 72 65 73 2c 20 79 72 65 73 3b |ype, xres, yres;|
00003f90 0a 20 20 69 6e 74 20 73 70 72 73 69 7a 65 2c 20 |. int sprsize, |
00003fa0 65 78 74 72 61 3b 0a 0a 20 20 73 69 2e 74 61 67 |extra;.. si.tag|
00003fb0 20 3d 20 73 70 72 69 74 65 5f 69 64 5f 6e 61 6d | = sprite_id_nam|
00003fc0 65 3b 0a 20 20 73 69 2e 73 2e 6e 61 6d 65 20 3d |e;. si.s.name =|
00003fd0 20 6e 61 6d 65 3b 0a 20 20 69 66 20 28 77 69 6d | name;. if (wim|
00003fe0 70 74 5f 63 6f 6d 70 6c 61 69 6e 28 73 70 72 69 |pt_complain(spri|
00003ff0 74 65 5f 73 65 6c 65 63 74 5f 72 70 28 28 73 70 |te_select_rp((sp|
00004000 72 69 74 65 5f 61 72 65 61 20 2a 29 2a 73 61 2c |rite_area *)*sa,|
00004010 20 26 73 69 2c 20 28 73 70 72 69 74 65 5f 70 74 | &si, (sprite_pt|
00004020 72 20 2a 29 26 73 68 29 29 29 20 72 65 74 75 72 |r *)&sh))) retur|
00004030 6e 20 46 41 4c 53 45 3b 0a 20 20 6d 6f 64 65 77 |n FALSE;. modew|
00004040 6f 72 64 20 3d 20 73 68 2d 3e 6d 6f 64 65 3b 0a |ord = sh->mode;.|
00004050 20 20 74 79 70 65 20 3d 20 6d 6f 64 65 77 6f 72 | type = modewor|
00004060 64 3e 3e 32 37 3b 0a 20 20 78 72 65 73 20 3d 20 |d>>27;. xres = |
00004070 28 6d 6f 64 65 77 6f 72 64 3e 3e 20 31 29 26 30 |(modeword>> 1)&0|
00004080 78 31 66 66 66 3b 0a 20 20 79 72 65 73 20 3d 20 |x1fff;. yres = |
00004090 28 6d 6f 64 65 77 6f 72 64 3e 3e 31 34 29 26 30 |(modeword>>14)&0|
000040a0 78 31 66 66 66 3b 0a 0a 20 20 69 66 20 28 74 79 |x1fff;.. if (ty|
000040b0 70 65 21 3d 36 29 20 7b 0a 20 20 20 20 77 65 72 |pe!=6) {. wer|
000040c0 72 28 30 2c 20 22 53 70 72 69 74 65 20 69 73 20 |r(0, "Sprite is |
000040d0 6e 6f 74 20 6e 65 77 20 74 79 70 65 20 36 20 28 |not new type 6 (|
000040e0 33 32 62 70 70 29 20 2d 20 49 20 6f 6e 6c 79 20 |32bpp) - I only |
000040f0 75 6e 64 65 72 73 74 61 6e 64 20 74 68 65 73 65 |understand these|
00004100 21 22 29 3b 0a 20 20 20 20 72 65 74 75 72 6e 20 |!");. return |
00004110 46 41 4c 53 45 3b 0a 20 20 7d 0a 0a 20 20 78 20 |FALSE;. }.. x |
00004120 3d 20 31 2b 73 68 2d 3e 77 69 64 74 68 3b 0a 20 |= 1+sh->width;. |
00004130 20 79 20 3d 20 73 68 2d 3e 68 65 69 67 68 74 20 | y = sh->height |
00004140 2b 20 31 3b 0a 20 20 73 70 72 73 69 7a 65 20 3d |+ 1;. sprsize =|
00004150 20 73 69 7a 65 6f 66 28 73 70 72 69 74 65 5f 68 | sizeof(sprite_h|
00004160 65 61 64 65 72 29 2b 79 2a 28 28 78 2b 33 29 26 |eader)+y*((x+3)&|
00004170 7e 33 29 3b 0a 0a 20 20 65 78 74 72 61 20 3d 20 |~3);.. extra = |
00004180 73 70 72 73 69 7a 65 20 2d 20 28 64 69 74 68 61 |sprsize - (ditha|
00004190 72 65 61 2d 3e 73 69 7a 65 20 2d 20 64 69 74 68 |rea->size - dith|
000041a0 61 72 65 61 2d 3e 66 72 65 65 6f 66 66 29 3b 0a |area->freeoff);.|
000041b0 20 20 69 66 20 28 65 78 74 72 61 3e 30 29 20 7b | if (extra>0) {|
000041c0 0a 20 20 20 20 6e 65 77 73 69 7a 65 20 3d 20 64 |. newsize = d|
000041d0 69 74 68 61 72 65 61 2d 3e 73 69 7a 65 20 2b 20 |itharea->size + |
000041e0 65 78 74 72 61 3b 0a 20 20 20 20 69 66 20 28 21 |extra;. if (!|
000041f0 66 6c 65 78 5f 65 78 74 65 6e 64 28 28 66 6c 65 |flex_extend((fle|
00004200 78 5f 70 74 72 29 26 64 69 74 68 61 72 65 61 2c |x_ptr)&ditharea,|
00004210 20 6e 65 77 73 69 7a 65 29 29 20 7b 0a 20 20 20 | newsize)) {. |
00004220 20 20 20 77 65 72 72 28 30 2c 20 22 4e 6f 20 72 | werr(0, "No r|
00004230 6f 6f 6d 20 66 6f 72 20 64 65 73 74 69 6e 61 74 |oom for destinat|
00004240 69 6f 6e 20 73 70 72 69 74 65 22 29 3b 0a 20 20 |ion sprite");. |
00004250 20 20 20 20 72 65 74 75 72 6e 20 46 41 4c 53 45 | return FALSE|
00004260 3b 0a 20 20 20 20 7d 0a 20 20 20 20 64 69 74 68 |;. }. dith|
00004270 61 72 65 61 2d 3e 73 69 7a 65 20 3d 20 6e 65 77 |area->size = new|
00004280 73 69 7a 65 3b 0a 20 20 20 20 2f 2a 61 62 6f 76 |size;. /*abov|
00004290 65 20 66 6c 65 78 20 77 69 6c 6c 20 6c 69 6b 65 |e flex will like|
000042a0 6c 79 20 68 61 76 65 20 69 6e 76 61 6c 69 64 61 |ly have invalida|
000042b0 74 65 64 20 73 68 2c 20 73 6f 20 72 65 63 61 6c |ted sh, so recal|
000042c0 63 2a 2f 0a 20 20 20 20 69 66 20 28 77 69 6d 70 |c*/. if (wimp|
000042d0 74 5f 63 6f 6d 70 6c 61 69 6e 28 73 70 72 69 74 |t_complain(sprit|
000042e0 65 5f 73 65 6c 65 63 74 5f 72 70 28 28 73 70 72 |e_select_rp((spr|
000042f0 69 74 65 5f 61 72 65 61 20 2a 29 2a 73 61 2c 20 |ite_area *)*sa, |
00004300 26 73 69 2c 20 28 73 70 72 69 74 65 5f 70 74 72 |&si, (sprite_ptr|
00004310 20 2a 29 26 73 68 29 29 29 20 72 65 74 75 72 6e | *)&sh))) return|
00004320 20 46 41 4c 53 45 3b 0a 20 20 7d 0a 0a 20 20 64 | FALSE;. }.. d|
00004330 65 20 3d 20 6d 61 6c 6c 6f 63 28 78 2a 73 69 7a |e = malloc(x*siz|
00004340 65 6f 66 28 63 6f 6d 70 61 63 74 5f 64 69 74 68 |eof(compact_dith|
00004350 65 72 5f 65 72 72 6f 72 29 29 3b 0a 20 20 69 66 |er_error));. if|
00004360 20 28 64 65 3d 3d 30 29 20 7b 0a 20 20 20 20 77 | (de==0) {. w|
00004370 65 72 72 28 30 2c 20 22 4e 6f 20 72 6f 6f 6d 20 |err(0, "No room |
00004380 66 6f 72 20 64 69 74 68 65 72 20 65 72 72 6f 72 |for dither error|
00004390 20 6c 69 6e 65 20 62 75 66 66 65 72 22 29 3b 0a | line buffer");.|
000043a0 20 20 20 20 72 65 74 75 72 6e 20 46 41 4c 53 45 | return FALSE|
000043b0 3b 0a 20 20 7d 0a 0a 20 20 65 3d 73 70 72 69 74 |;. }.. e=sprit|
000043c0 65 5f 63 72 65 61 74 65 5f 72 70 28 64 69 74 68 |e_create_rp(dith|
000043d0 61 72 65 61 2c 20 6e 61 6d 65 2c 20 73 70 72 69 |area, name, spri|
000043e0 74 65 5f 6e 6f 70 61 6c 65 74 74 65 2c 20 78 2c |te_nopalette, x,|
000043f0 20 79 2c 20 78 72 65 73 3d 3d 32 2a 79 72 65 73 | y, xres==2*yres|
00004400 20 3f 20 31 35 20 3a 20 32 38 2c 20 28 73 70 72 | ? 15 : 28, (spr|
00004410 69 74 65 5f 70 74 72 20 2a 29 26 64 68 29 3b 0a |ite_ptr *)&dh);.|
00004420 20 20 69 66 20 28 65 29 20 7b 0a 20 20 20 20 66 | if (e) {. f|
00004430 72 65 65 28 64 65 29 3b 0a 20 20 20 20 77 69 6d |ree(de);. wim|
00004440 70 74 5f 63 6f 6d 70 6c 61 69 6e 28 65 29 3b 0a |pt_complain(e);.|
00004450 20 20 20 20 72 65 74 75 72 6e 20 46 41 4c 53 45 | return FALSE|
00004460 3b 0a 20 20 7d 0a 20 20 64 70 20 3d 20 28 63 68 |;. }. dp = (ch|
00004470 61 72 20 2a 29 64 68 20 2b 20 64 68 2d 3e 69 6d |ar *)dh + dh->im|
00004480 61 67 65 3b 0a 20 20 73 70 20 3d 20 28 63 68 61 |age;. sp = (cha|
00004490 72 20 2a 29 73 68 20 2b 20 73 68 2d 3e 69 6d 61 |r *)sh + sh->ima|
000044a0 67 65 3b 0a 0a 20 20 67 32 35 36 77 6e 70 5f 69 |ge;.. g256wnp_i|
000044b0 6e 74 65 72 6e 61 6c 28 78 2c 20 79 2c 20 64 70 |nternal(x, y, dp|
000044c0 2c 20 73 70 2c 20 64 65 2c 20 6c 75 74 2c 20 70 |, sp, de, lut, p|
000044d0 61 6c 32 35 36 64 65 66 29 3b 0a 0a 20 20 66 72 |al256def);.. fr|
000044e0 65 65 28 64 65 29 3b 0a 20 20 72 65 74 75 72 6e |ee(de);. return|
000044f0 20 54 52 55 45 3b 0a 7d 0a 0a 2f 2a 2a 2a 2a 2a | TRUE;.}../*****|
00004500 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00004510 2a 2a 2a 2a 2a 2a 2a 2a 2a 20 69 6e 69 74 69 61 |********* initia|
00004520 6c 69 73 61 74 69 6f 6e 20 2a 2a 2a 2a 2a 2a 2a |lisation *******|
00004530 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00004540 2a 2a 2a 2a 2a 2a 2a 2f 0a 0a 2f 2a 2d 2d 2d 20 |*******/../*--- |
00004550 49 6e 69 74 69 61 6c 69 73 65 20 74 68 65 20 70 |Initialise the p|
00004560 72 6f 67 72 61 6d 2c 20 72 65 74 75 72 6e 69 6e |rogram, returnin|
00004570 67 20 54 52 55 45 20 69 66 20 69 74 20 77 61 73 |g TRUE if it was|
00004580 20 61 6c 6c 20 4f 4b 2e 20 2d 2d 2d 2a 2f 0a 73 | all OK. ---*/.s|
00004590 74 61 74 69 63 20 42 4f 4f 4c 20 69 6e 69 74 69 |tatic BOOL initi|
000045a0 61 6c 69 73 65 28 69 6e 74 20 61 72 67 63 2c 20 |alise(int argc, |
000045b0 63 68 61 72 20 2a 61 72 67 76 5b 5d 29 0a 7b 0a |char *argv[]).{.|
000045c0 20 20 6f 73 5f 66 69 6c 65 73 74 72 20 66 3b 0a | os_filestr f;.|
000045d0 20 20 69 6e 74 20 69 3b 0a 20 20 63 68 61 72 20 | int i;. char |
000045e0 70 74 5b 31 35 33 36 5d 2c 20 2a 70 3b 0a 0a 20 |pt[1536], *p;.. |
000045f0 20 69 66 20 28 61 72 67 63 3e 3d 32 29 20 7b 0a | if (argc>=2) {.|
00004600 20 20 20 20 69 3d 61 74 6f 69 28 61 72 67 76 5b | i=atoi(argv[|
00004610 31 5d 29 3b 0a 20 20 20 20 73 61 76 65 33 32 20 |1]);. save32 |
00004620 3d 20 69 26 31 20 3f 20 54 52 55 45 20 3a 20 46 |= i&1 ? TRUE : F|
00004630 41 4c 53 45 3b 0a 20 20 20 20 73 61 76 65 38 20 |ALSE;. save8 |
00004640 20 3d 20 69 26 32 20 3f 20 54 52 55 45 20 3a 20 | = i&2 ? TRUE : |
00004650 46 41 4c 53 45 3b 0a 20 20 7d 0a 0a 20 20 66 2e |FALSE;. }.. f.|
00004660 61 63 74 69 6f 6e 3d 32 35 35 3b 0a 20 20 66 2e |action=255;. f.|
00004670 6e 61 6d 65 3d 22 3c 53 63 61 70 65 32 24 44 69 |name="<Scape2$Di|
00004680 72 3e 2e 72 65 73 6f 75 72 63 65 73 2e 6c 75 74 |r>.resources.lut|
00004690 22 3b 0a 20 20 66 2e 6c 6f 61 64 61 64 64 72 3d |";. f.loadaddr=|
000046a0 28 69 6e 74 29 6c 75 74 3b 0a 20 20 66 2e 65 78 |(int)lut;. f.ex|
000046b0 65 63 61 64 64 72 3d 30 3b 0a 20 20 69 66 20 28 |ecaddr=0;. if (|
000046c0 77 69 6d 70 74 5f 63 6f 6d 70 6c 61 69 6e 28 6f |wimpt_complain(o|
000046d0 73 5f 66 69 6c 65 28 26 66 29 29 29 20 72 65 74 |s_file(&f))) ret|
000046e0 75 72 6e 20 46 41 4c 53 45 3b 0a 20 20 69 66 20 |urn FALSE;. if |
000046f0 28 66 2e 73 74 61 72 74 21 3d 33 32 37 36 38 29 |(f.start!=32768)|
00004700 20 7b 20 2f 2a 66 2e 73 74 61 72 74 20 69 73 20 | { /*f.start is |
00004710 72 34 2c 20 68 65 72 65 20 68 6f 6c 64 69 6e 67 |r4, here holding|
00004720 20 66 69 6c 65 20 6c 65 6e 67 74 68 2a 2f 0a 20 | file length*/. |
00004730 20 20 20 77 65 72 72 28 30 2c 20 22 49 6e 74 65 | werr(0, "Inte|
00004740 72 6e 61 6c 3a 20 72 65 73 6f 75 72 63 65 73 2e |rnal: resources.|
00004750 6c 75 74 20 66 69 6c 65 20 63 6f 72 72 75 70 74 |lut file corrupt|
00004760 65 64 20 2d 20 73 68 6f 75 6c 64 20 62 65 20 33 |ed - should be 3|
00004770 32 37 36 38 20 62 79 74 65 73 20 6c 6f 6e 67 22 |2768 bytes long"|
00004780 29 3b 0a 20 20 20 20 72 65 74 75 72 6e 20 46 41 |);. return FA|
00004790 4c 53 45 3b 0a 20 20 7d 0a 0a 20 20 66 2e 61 63 |LSE;. }.. f.ac|
000047a0 74 69 6f 6e 3d 32 35 35 3b 0a 20 20 66 2e 6e 61 |tion=255;. f.na|
000047b0 6d 65 3d 22 3c 53 63 61 70 65 32 24 44 69 72 3e |me="<Scape2$Dir>|
000047c0 2e 72 65 73 6f 75 72 63 65 73 2e 70 61 6c 32 35 |.resources.pal25|
000047d0 36 64 65 66 22 3b 0a 20 20 66 2e 6c 6f 61 64 61 |6def";. f.loada|
000047e0 64 64 72 3d 28 69 6e 74 29 70 61 6c 32 35 36 64 |ddr=(int)pal256d|
000047f0 65 66 3b 0a 20 20 66 2e 65 78 65 63 61 64 64 72 |ef;. f.execaddr|
00004800 3d 30 3b 0a 20 20 69 66 20 28 77 69 6d 70 74 5f |=0;. if (wimpt_|
00004810 63 6f 6d 70 6c 61 69 6e 28 6f 73 5f 66 69 6c 65 |complain(os_file|
00004820 28 26 66 29 29 29 20 72 65 74 75 72 6e 20 46 41 |(&f))) return FA|
00004830 4c 53 45 3b 0a 20 20 69 66 20 28 66 2e 73 74 61 |LSE;. if (f.sta|
00004840 72 74 21 3d 31 30 32 34 29 20 7b 20 2f 2a 66 2e |rt!=1024) { /*f.|
00004850 73 74 61 72 74 20 69 73 20 72 34 2c 20 68 65 72 |start is r4, her|
00004860 65 20 68 6f 6c 64 69 6e 67 20 66 69 6c 65 20 6c |e holding file l|
00004870 65 6e 67 74 68 2a 2f 0a 20 20 20 20 77 65 72 72 |ength*/. werr|
00004880 28 30 2c 20 22 49 6e 74 65 72 6e 61 6c 3a 20 72 |(0, "Internal: r|
00004890 65 73 6f 75 72 63 65 73 2e 70 61 6c 32 35 36 64 |esources.pal256d|
000048a0 65 66 20 66 69 6c 65 20 63 6f 72 72 75 70 74 65 |ef file corrupte|
000048b0 64 20 2d 20 73 68 6f 75 6c 64 20 62 65 20 31 30 |d - should be 10|
000048c0 32 34 20 62 79 74 65 73 20 6c 6f 6e 67 22 29 3b |24 bytes long");|
000048d0 0a 20 20 20 20 72 65 74 75 72 6e 20 46 41 4c 53 |. return FALS|
000048e0 45 3b 0a 20 20 7d 0a 0a 20 20 66 2e 61 63 74 69 |E;. }.. f.acti|
000048f0 6f 6e 3d 32 35 35 3b 0a 20 20 66 2e 6e 61 6d 65 |on=255;. f.name|
00004900 3d 22 3c 53 63 61 70 65 32 24 44 69 72 3e 2e 72 |="<Scape2$Dir>.r|
00004910 65 73 6f 75 72 63 65 73 2e 65 61 72 74 68 22 3b |esources.earth";|
00004920 0a 20 20 66 2e 6c 6f 61 64 61 64 64 72 3d 28 69 |. f.loadaddr=(i|
00004930 6e 74 29 70 74 3b 0a 20 20 66 2e 65 78 65 63 61 |nt)pt;. f.execa|
00004940 64 64 72 3d 30 3b 0a 20 20 69 66 20 28 77 69 6d |ddr=0;. if (wim|
00004950 70 74 5f 63 6f 6d 70 6c 61 69 6e 28 6f 73 5f 66 |pt_complain(os_f|
00004960 69 6c 65 28 26 66 29 29 29 20 72 65 74 75 72 6e |ile(&f))) return|
00004970 20 46 41 4c 53 45 3b 0a 20 20 69 66 20 28 66 2e | FALSE;. if (f.|
00004980 73 74 61 72 74 21 3d 31 35 33 36 29 20 7b 20 2f |start!=1536) { /|
00004990 2a 66 2e 73 74 61 72 74 20 69 73 20 72 34 2c 20 |*f.start is r4, |
000049a0 68 65 72 65 20 68 6f 6c 64 69 6e 67 20 66 69 6c |here holding fil|
000049b0 65 20 6c 65 6e 67 74 68 2a 2f 0a 20 20 20 20 77 |e length*/. w|
000049c0 65 72 72 28 30 2c 20 22 49 6e 74 65 72 6e 61 6c |err(0, "Internal|
000049d0 3a 20 72 65 73 6f 75 72 63 65 73 2e 65 61 72 74 |: resources.eart|
000049e0 68 20 66 69 6c 65 20 63 6f 72 72 75 70 74 65 64 |h file corrupted|
000049f0 20 2d 20 73 68 6f 75 6c 64 20 62 65 20 31 35 33 | - should be 153|
00004a00 36 20 62 79 74 65 73 20 6c 6f 6e 67 22 29 3b 0a |6 bytes long");.|
00004a10 20 20 20 20 72 65 74 75 72 6e 20 46 41 4c 53 45 | return FALSE|
00004a20 3b 0a 20 20 7d 0a 20 20 66 6f 72 20 28 69 3d 30 |;. }. for (i=0|
00004a30 2c 20 70 3d 70 74 3b 20 69 3c 32 35 36 3b 20 69 |, p=pt; i<256; i|
00004a40 2b 2b 2c 20 70 2b 3d 36 29 20 65 61 72 74 68 5b |++, p+=6) earth[|
00004a50 69 5d 3d 70 5b 33 5d 2b 28 70 5b 34 5d 3c 3c 38 |i]=p[3]+(p[4]<<8|
00004a60 29 2b 28 70 5b 35 5d 3c 3c 31 36 29 3b 0a 0a 20 |)+(p[5]<<16);.. |
00004a70 20 2f 2a 20 41 6c 6c 20 77 65 6e 74 20 6f 6b 20 | /* All went ok |
00004a80 2a 2f 0a 20 20 72 65 74 75 72 6e 20 54 52 55 45 |*/. return TRUE|
00004a90 3b 0a 7d 0a 0a 2f 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |;.}../**********|
00004aa0 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
00004ab0 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 20 |*************** |
00004ac0 6d 61 69 6e 20 66 75 6e 63 74 69 6f 6e 73 20 2a |main functions *|
00004ad0 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a 2a |****************|
*
00004b00 2a 2f 0a 0a 42 4f 4f 4c 20 67 72 61 70 68 28 69 |*/..BOOL graph(i|
00004b10 6e 74 20 2a 64 61 74 61 2c 20 69 6e 74 20 6c 69 |nt *data, int li|
00004b20 6d 69 74 29 0a 7b 0a 20 20 69 6e 74 20 6e 3d 31 |mit).{. int n=1|
00004b30 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 20 69 |<<limit;. int i|
00004b40 2c 20 6a 2c 20 6b 3b 0a 20 20 69 6e 74 20 70 6c |, j, k;. int pl|
00004b50 2c 20 70 72 3b 0a 20 20 69 6e 74 20 76 6c 2c 20 |, pr;. int vl, |
00004b60 76 72 2c 20 76 2c 20 76 32 3b 0a 20 20 69 6e 74 |vr, v, v2;. int|
00004b70 20 68 6f 72 69 7a 6f 6e 5b 31 30 32 36 5d 3b 0a | horizon[1026];.|
00004b80 20 20 69 6e 74 20 68 6f 72 69 7a 6f 6e 32 5b 31 | int horizon2[1|
00004b90 30 32 36 5d 3b 0a 20 20 69 6e 74 20 63 6f 73 20 |026];. int cos |
00004ba0 3d 20 63 6f 73 31 36 28 74 69 6c 74 29 3b 0a 20 |= cos16(tilt);. |
00004bb0 20 69 6e 74 20 73 69 6e 20 3d 20 73 69 6e 31 36 | int sin = sin16|
00004bc0 28 74 69 6c 74 29 3b 0a 20 20 69 6e 74 20 69 6d |(tilt);. int im|
00004bd0 61 67 65 5f 79 6c 3d 31 30 32 34 3b 0a 20 20 69 |age_yl=1024;. i|
00004be0 6e 74 20 69 6d 61 67 65 5f 79 68 3d 30 3b 0a 0a |nt image_yh=0;..|
00004bf0 20 20 66 6f 72 20 28 69 3d 30 3b 20 69 3c 31 30 | for (i=0; i<10|
00004c00 32 36 3b 20 69 2b 2b 29 20 7b 0a 20 20 20 20 68 |26; i++) {. h|
00004c10 6f 72 69 7a 6f 6e 5b 69 5d 3d 30 3b 0a 20 20 20 |orizon[i]=0;. |
00004c20 20 68 6f 72 69 7a 6f 6e 32 5b 69 5d 3d 31 30 32 | horizon2[i]=102|
00004c30 34 3b 0a 20 20 7d 0a 20 20 69 6d 33 32 5f 79 6c |4;. }. im32_yl|
00004c40 3d 31 30 32 34 3b 0a 20 20 62 62 63 5f 6f 72 69 |=1024;. bbc_ori|
00004c50 67 69 6e 28 36 34 30 2c 30 29 3b 0a 20 20 66 6f |gin(640,0);. fo|
00004c60 72 20 28 69 3d 30 3b 20 69 3c 6e 3b 20 69 2b 2b |r (i=0; i<n; i++|
00004c70 29 20 7b 0a 20 20 20 20 66 6f 72 20 28 6a 3d 31 |) {. for (j=1|
00004c80 2c 20 70 6c 3d 2d 35 31 32 2c 20 76 6c 20 3d 20 |, pl=-512, vl = |
00004c90 64 61 74 61 5b 30 2b 6e 2a 69 5d 3b 20 6a 3c 6e |data[0+n*i]; j<n|
00004ca0 3b 20 6a 2b 2b 2c 20 70 6c 3d 70 72 2c 20 76 6c |; j++, pl=pr, vl|
00004cb0 3d 76 72 29 20 7b 0a 20 20 20 20 20 20 66 6f 72 |=vr) {. for|
00004cc0 20 28 70 72 3d 70 6c 2b 28 31 30 32 34 3e 3e 6c | (pr=pl+(1024>>l|
00004cd0 69 6d 69 74 29 2c 20 76 72 3d 64 61 74 61 5b 6a |imit), vr=data[j|
00004ce0 2b 6e 2a 69 5d 2c 20 6b 3d 70 6c 2c 20 76 3d 76 |+n*i], k=pl, v=v|
00004cf0 6c 3b 20 6b 3c 70 72 3b 20 6b 2b 3d 32 2c 20 76 |l; k<pr; k+=2, v|
00004d00 2b 3d 28 76 72 2d 76 6c 29 3e 3e 28 39 2d 6c 69 |+=(vr-vl)>>(9-li|
00004d10 6d 69 74 29 29 20 7b 0a 20 20 20 20 20 20 20 20 |mit)) {. |
00004d20 76 32 20 3d 20 28 6d 75 6c 5f 66 72 61 63 31 36 |v2 = (mul_frac16|
00004d30 28 76 2c 20 63 6f 73 29 2b 6d 75 6c 5f 66 72 61 |(v, cos)+mul_fra|
00004d40 63 31 36 28 69 3c 3c 28 31 36 2d 6c 69 6d 69 74 |c16(i<<(16-limit|
00004d50 29 2c 20 73 69 6e 29 29 2a 33 2f 32 35 36 3b 0a |), sin))*3/256;.|
00004d60 20 20 20 20 20 20 20 20 69 66 20 28 69 3d 3d 30 | if (i==0|
00004d70 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 68 6f |) {. ho|
00004d80 72 69 7a 6f 6e 32 5b 6b 2b 35 31 32 5d 3d 76 32 |rizon2[k+512]=v2|
00004d90 3b 0a 20 20 20 20 20 20 20 20 20 20 69 66 20 28 |;. if (|
00004da0 76 32 3c 69 6d 33 32 5f 79 6c 29 20 69 6d 33 32 |v2<im32_yl) im32|
00004db0 5f 79 6c 3d 76 32 3b 0a 20 20 20 20 20 20 20 20 |_yl=v2;. |
00004dc0 7d 0a 20 20 20 20 20 20 20 20 69 66 20 28 68 6f |}. if (ho|
00004dd0 72 69 7a 6f 6e 5b 6b 2b 35 31 32 5d 3c 76 32 29 |rizon[k+512]<v2)|
00004de0 20 7b 0a 20 20 20 20 20 20 20 20 20 20 62 62 63 | {. bbc|
00004df0 5f 67 63 6f 6c 28 30 2c 20 34 2b 36 2a 69 2f 6e |_gcol(0, 4+6*i/n|
00004e00 29 3b 0a 20 20 20 20 20 20 20 20 20 20 62 62 63 |);. bbc|
00004e10 5f 70 6c 6f 74 28 36 39 2c 20 6b 2c 20 76 32 29 |_plot(69, k, v2)|
00004e20 3b 0a 20 20 20 20 20 20 20 20 20 20 68 6f 72 69 |;. hori|
00004e30 7a 6f 6e 5b 6b 2b 35 31 32 5d 3d 76 32 3b 0a 20 |zon[k+512]=v2;. |
00004e40 20 20 20 20 20 20 20 7d 0a 20 20 20 20 20 20 20 | }. |
00004e50 20 65 6c 73 65 20 69 66 20 28 68 6f 72 69 7a 6f | else if (horizo|
00004e60 6e 32 5b 6b 2b 35 31 32 5d 3e 76 32 29 20 7b 0a |n2[k+512]>v2) {.|
00004e70 20 20 20 20 20 20 20 20 20 20 62 62 63 5f 67 63 | bbc_gc|
00004e80 6f 6c 28 30 2c 20 31 30 2b 36 2a 69 2f 6e 29 3b |ol(0, 10+6*i/n);|
00004e90 0a 20 20 20 20 20 20 20 20 20 20 62 62 63 5f 70 |. bbc_p|
00004ea0 6c 6f 74 28 36 39 2c 20 6b 2c 20 76 32 29 3b 0a |lot(69, k, v2);.|
00004eb0 20 20 20 20 20 20 20 20 20 20 68 6f 72 69 7a 6f | horizo|
00004ec0 6e 32 5b 6b 2b 35 31 32 5d 3d 76 32 3b 0a 20 20 |n2[k+512]=v2;. |
00004ed0 20 20 20 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a | }. }.|
00004ee0 20 20 20 20 7d 0a 20 20 20 20 76 32 20 3d 20 28 | }. v2 = (|
00004ef0 6d 75 6c 5f 66 72 61 63 31 36 28 76 72 2c 20 63 |mul_frac16(vr, c|
00004f00 6f 73 29 2b 6d 75 6c 5f 66 72 61 63 31 36 28 69 |os)+mul_frac16(i|
00004f10 3c 3c 28 31 36 2d 6c 69 6d 69 74 29 2c 20 73 69 |<<(16-limit), si|
00004f20 6e 29 29 2a 33 2f 32 35 36 3b 0a 20 20 20 20 69 |n))*3/256;. i|
00004f30 66 20 28 69 3d 3d 30 29 20 68 6f 72 69 7a 6f 6e |f (i==0) horizon|
00004f40 32 5b 6b 2b 35 31 32 5d 3d 76 32 3b 0a 20 20 20 |2[k+512]=v2;. |
00004f50 20 69 66 20 28 68 6f 72 69 7a 6f 6e 5b 6b 2b 35 | if (horizon[k+5|
00004f60 31 32 5d 3c 76 32 29 20 7b 0a 20 20 20 20 20 20 |12]<v2) {. |
00004f70 62 62 63 5f 67 63 6f 6c 28 30 2c 20 34 2b 36 2a |bbc_gcol(0, 4+6*|
00004f80 69 2f 6e 29 3b 0a 20 20 20 20 20 20 62 62 63 5f |i/n);. bbc_|
00004f90 70 6c 6f 74 28 36 39 2c 20 6b 2c 20 76 32 29 3b |plot(69, k, v2);|
00004fa0 0a 20 20 20 20 20 20 68 6f 72 69 7a 6f 6e 5b 6b |. horizon[k|
00004fb0 2b 35 31 32 5d 3d 76 32 3b 0a 20 20 20 20 7d 0a |+512]=v2;. }.|
00004fc0 20 20 20 20 65 6c 73 65 20 69 66 20 28 68 6f 72 | else if (hor|
00004fd0 69 7a 6f 6e 32 5b 6b 2b 35 31 32 5d 3e 76 32 29 |izon2[k+512]>v2)|
00004fe0 20 7b 0a 20 20 20 20 20 20 62 62 63 5f 67 63 6f | {. bbc_gco|
00004ff0 6c 28 30 2c 20 31 30 2b 36 2a 69 2f 6e 29 3b 0a |l(0, 10+6*i/n);.|
00005000 20 20 20 20 20 20 62 62 63 5f 70 6c 6f 74 28 36 | bbc_plot(6|
00005010 39 2c 20 6b 2c 20 76 32 29 3b 0a 20 20 20 20 20 |9, k, v2);. |
00005020 20 68 6f 72 69 7a 6f 6e 32 5b 6b 2b 35 31 32 5d | horizon2[k+512]|
00005030 3d 76 32 3b 0a 20 20 20 20 7d 0a 20 20 7d 0a 20 |=v2;. }. }. |
00005040 20 66 6f 72 20 28 69 3d 30 3b 20 69 3c 31 30 32 | for (i=0; i<102|
00005050 36 3b 20 69 2b 2b 29 20 7b 0a 20 20 20 20 69 66 |6; i++) {. if|
00005060 20 28 68 6f 72 69 7a 6f 6e 5b 69 5d 20 3e 69 6d | (horizon[i] >im|
00005070 61 67 65 5f 79 68 29 20 69 6d 61 67 65 5f 79 68 |age_yh) image_yh|
00005080 3d 68 6f 72 69 7a 6f 6e 5b 69 5d 3b 0a 20 20 20 |=horizon[i];. |
00005090 20 69 66 20 28 68 6f 72 69 7a 6f 6e 32 5b 69 5d | if (horizon2[i]|
000050a0 3c 69 6d 61 67 65 5f 79 6c 29 20 69 6d 61 67 65 |<image_yl) image|
000050b0 5f 79 6c 3d 68 6f 72 69 7a 6f 6e 32 5b 69 5d 3b |_yl=horizon2[i];|
000050c0 0a 20 20 7d 0a 20 20 62 62 63 5f 6d 6f 76 65 28 |. }. bbc_move(|
000050d0 2d 35 31 32 2c 20 69 6d 61 67 65 5f 79 6c 29 3b |-512, image_yl);|
000050e0 0a 20 20 62 62 63 5f 6d 6f 76 65 28 35 31 32 2c |. bbc_move(512,|
000050f0 20 69 6d 61 67 65 5f 79 68 29 3b 0a 20 20 62 62 | image_yh);. bb|
00005100 63 5f 70 6c 6f 74 28 31 38 39 2c 20 2d 35 31 32 |c_plot(189, -512|
00005110 2c 20 35 31 32 2b 28 69 6d 61 67 65 5f 79 6c 2d |, 512+(image_yl-|
00005120 69 6d 61 67 65 5f 79 68 29 2f 32 29 3b 0a 0a 20 |image_yh)/2);.. |
00005130 20 69 6d 33 32 5f 79 6c 09 3d 20 69 6d 33 32 5f | im32_yl.= im32_|
00005140 79 6c 2d 32 3b 0a 20 20 69 6d 33 32 5f 68 65 69 |yl-2;. im32_hei|
00005150 67 68 74 09 3d 20 33 20 2b 20 28 20 28 69 6d 61 |ght.= 3 + ( (ima|
00005160 67 65 5f 79 68 2d 69 6d 33 32 5f 79 6c 2d 32 29 |ge_yh-im32_yl-2)|
00005170 20 3e 3e 20 28 31 30 2d 6c 69 6d 69 74 29 20 29 | >> (10-limit) )|
00005180 3b 0a 0a 20 20 72 65 74 75 72 6e 20 54 52 55 45 |;.. return TRUE|
00005190 3b 0a 7d 0a 0a 6d 6f 6e 69 74 6f 72 20 72 65 61 |;.}..monitor rea|
000051a0 64 5f 6d 6f 6e 69 74 6f 72 5f 74 79 70 65 28 76 |d_monitor_type(v|
000051b0 6f 69 64 29 0a 7b 0a 20 20 69 6e 74 20 72 32 3b |oid).{. int r2;|
000051c0 0a 20 20 6f 73 5f 73 77 69 33 72 28 36 2c 20 31 |. os_swi3r(6, 1|
000051d0 36 31 2c 31 33 33 2c 30 2c 20 30 2c 30 2c 26 72 |61,133,0, 0,0,&r|
000051e0 32 29 3b 0a 20 20 69 66 20 28 28 72 32 2f 34 20 |2);. if ((r2/4 |
000051f0 26 20 33 29 20 3d 3d 20 31 29 20 72 65 74 75 72 |& 3) == 1) retur|
00005200 6e 20 6d 75 6c 74 69 3b 0a 20 20 65 6c 73 65 20 |n multi;. else |
00005210 72 65 74 75 72 6e 20 6e 6f 6e 6d 75 6c 74 69 3b |return nonmulti;|
00005220 0a 7d 0a 0a 69 6e 74 20 6d 61 69 6e 28 69 6e 74 |.}..int main(int|
00005230 20 61 72 67 63 2c 20 63 68 61 72 20 2a 61 72 67 | argc, char *arg|
00005240 76 5b 5d 29 0a 7b 0a 20 20 69 6e 74 20 73 65 65 |v[]).{. int see|
00005250 64 2c 20 69 3b 0a 20 20 63 68 61 72 20 6b 65 79 |d, i;. char key|
00005260 3b 0a 20 20 64 6f 75 62 6c 65 20 64 69 6d 3d 32 |;. double dim=2|
00005270 2e 32 3b 0a 20 20 64 6f 75 62 6c 65 20 64 3b 0a |.2;. double d;.|
00005280 20 20 63 68 61 72 20 62 75 66 5b 31 30 32 34 5d | char buf[1024]|
00005290 3b 0a 20 20 69 6e 74 20 74 6f 72 69 65 6e 74 2c |;. int torient,|
000052a0 20 74 73 63 72 6f 6c 6c 78 2c 20 74 73 63 72 6f | tscrollx, tscro|
000052b0 6c 6c 79 3b 0a 20 20 6d 6f 6e 69 74 6f 72 74 79 |lly;. monitorty|
000052c0 70 65 20 3d 20 72 65 61 64 5f 6d 6f 6e 69 74 6f |pe = read_monito|
000052d0 72 5f 74 79 70 65 28 29 3b 0a 20 20 66 6c 65 78 |r_type();. flex|
000052e0 5f 69 6e 69 74 28 29 3b 0a 0a 20 20 69 66 20 28 |_init();.. if (|
000052f0 21 69 6e 69 74 69 61 6c 69 73 65 28 61 72 67 63 |!initialise(argc|
00005300 2c 20 61 72 67 76 29 29 20 72 65 74 75 72 6e 20 |, argv)) return |
00005310 30 3b 0a 0a 20 20 66 6f 72 20 28 3b 3b 29 20 7b |0;.. for (;;) {|
00005320 0a 0a 20 20 20 20 62 62 63 5f 6d 6f 64 65 28 6d |.. bbc_mode(m|
00005330 6f 6e 69 74 6f 72 74 79 70 65 3d 3d 6d 75 6c 74 |onitortype==mult|
00005340 69 20 3f 20 32 30 20 3a 20 31 32 29 3b 0a 20 20 |i ? 20 : 12);. |
00005350 20 20 70 72 69 6e 74 66 28 22 50 6c 65 61 73 65 | printf("Please|
00005360 20 73 65 6c 65 63 74 20 77 68 69 63 68 20 6c 61 | select which la|
00005370 6e 64 73 63 61 70 65 20 74 6f 20 63 6f 6e 73 74 |ndscape to const|
00005380 72 75 63 74 20 28 61 6e 79 20 69 6e 74 65 67 65 |ruct (any intege|
00005390 72 29 20 22 29 3b 0a 20 20 20 20 73 63 61 6e 66 |r) ");. scanf|
000053a0 28 22 25 69 22 2c 20 26 73 65 65 64 29 3b 0a 20 |("%i", &seed);. |
000053b0 20 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 70 72 | do {. pr|
000053c0 69 6e 74 66 28 22 50 6c 65 61 73 65 20 65 6e 74 |intf("Please ent|
000053d0 65 72 20 73 63 61 70 65 20 72 65 73 6f 6c 75 74 |er scape resolut|
000053e0 69 6f 6e 20 28 73 75 67 67 65 73 74 20 35 2e 2e |ion (suggest 5..|
000053f0 2e 38 2c 20 74 68 6f 75 67 68 20 33 2e 2e 2e 39 |.8, though 3...9|
00005400 20 70 65 72 6d 69 73 73 69 62 6c 65 29 20 22 29 | permissible) ")|
00005410 3b 0a 20 20 20 20 20 20 73 63 61 6e 66 28 22 25 |;. scanf("%|
00005420 69 22 2c 20 26 6c 69 6d 69 74 29 3b 0a 20 20 20 |i", &limit);. |
00005430 20 7d 20 77 68 69 6c 65 20 28 6c 69 6d 69 74 3c | } while (limit<|
00005440 33 20 7c 7c 20 6c 69 6d 69 74 3e 39 29 3b 0a 20 |3 || limit>9);. |
00005450 20 20 20 4e 20 3d 20 31 3c 3c 6c 69 6d 69 74 3b | N = 1<<limit;|
00005460 0a 20 20 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 |. do {. |
00005470 70 72 69 6e 74 66 28 22 50 6c 65 61 73 65 20 65 |printf("Please e|
00005480 6e 74 65 72 20 73 63 61 70 65 20 64 69 6d 65 6e |nter scape dimen|
00005490 73 69 6f 6e 20 28 73 75 67 67 65 73 74 20 32 2e |sion (suggest 2.|
000054a0 31 35 2c 20 32 2d 33 20 70 65 72 6d 69 73 73 69 |15, 2-3 permissi|
000054b0 62 6c 65 29 20 22 29 3b 0a 20 20 20 20 20 20 73 |ble) ");. s|
000054c0 63 61 6e 66 28 22 25 4c 66 22 2c 20 26 64 69 6d |canf("%Lf", &dim|
000054d0 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 20 28 |);. } while (|
000054e0 64 69 6d 3c 32 20 7c 7c 20 64 69 6d 3e 33 29 3b |dim<2 || dim>3);|
000054f0 0a 20 20 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 |. do {. |
00005500 70 72 69 6e 74 66 28 22 50 6c 65 61 73 65 20 65 |printf("Please e|
00005510 6e 74 65 72 20 66 72 65 71 75 65 6e 63 79 20 72 |nter frequency r|
00005520 61 6e 67 65 20 28 73 75 67 67 65 73 74 20 31 2c |ange (suggest 1,|
00005530 20 30 2d 31 20 70 65 72 6d 69 73 73 69 62 6c 65 | 0-1 permissible|
00005540 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 61 6e |) ");. scan|
00005550 66 28 22 25 4c 66 22 2c 20 26 64 29 3b 0a 20 20 |f("%Lf", &d);. |
00005560 20 20 7d 20 77 68 69 6c 65 20 28 64 3c 30 20 7c | } while (d<0 ||
00005570 7c 20 64 3e 31 29 3b 0a 20 20 20 20 72 61 6e 67 || d>1);. rang|
00005580 65 20 3d 20 28 69 6e 74 29 28 30 2e 35 2b 64 2a |e = (int)(0.5+d*|
00005590 28 64 6f 75 62 6c 65 29 28 69 6e 74 29 28 4e 2f |(double)(int)(N/|
000055a0 32 29 29 3b 0a 20 20 20 20 69 66 20 28 72 61 6e |2));. if (ran|
000055b0 67 65 3c 30 29 20 72 61 6e 67 65 3d 30 3b 0a 20 |ge<0) range=0;. |
000055c0 20 20 20 69 66 20 28 72 61 6e 67 65 3e 4e 2f 32 | if (range>N/2|
000055d0 29 20 72 61 6e 67 65 3d 4e 2f 32 3b 0a 20 20 20 |) range=N/2;. |
000055e0 20 64 6f 20 7b 0a 20 20 20 20 20 20 70 72 69 6e | do {. prin|
000055f0 74 66 28 22 50 6c 65 61 73 65 20 65 6e 74 65 72 |tf("Please enter|
00005600 20 73 65 61 2d 6c 65 76 65 6c 20 28 73 75 67 67 | sea-level (sugg|
00005610 65 73 74 20 30 2c 20 30 2d 31 20 70 65 72 6d 69 |est 0, 0-1 permi|
00005620 73 73 69 62 6c 65 29 20 22 29 3b 0a 20 20 20 20 |ssible) ");. |
00005630 20 20 73 63 61 6e 66 28 22 25 4c 66 22 2c 20 26 | scanf("%Lf", &|
00005640 64 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 20 |d);. } while |
00005650 28 64 3c 30 20 7c 7c 20 64 3e 31 29 3b 0a 20 20 |(d<0 || d>1);. |
00005660 20 20 73 65 61 6c 65 76 65 6c 20 3d 20 28 69 6e | sealevel = (in|
00005670 74 29 28 64 2a 6f 6e 65 29 3b 0a 20 20 20 20 64 |t)(d*one);. d|
00005680 6f 20 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 |o {. printf|
00005690 28 22 50 6c 65 61 73 65 20 65 6e 74 65 72 20 70 |("Please enter p|
000056a0 6f 77 65 72 20 73 63 61 6c 69 6e 67 20 6c 61 77 |ower scaling law|
000056b0 20 28 73 75 67 67 65 73 74 20 31 2c 20 6f 72 20 | (suggest 1, or |
000056c0 30 2e 33 2d 33 2e 30 2c 20 74 68 6f 75 67 68 20 |0.3-3.0, though |
000056d0 3e 3d 20 30 20 70 65 72 6d 69 73 73 69 62 6c 65 |>= 0 permissible|
000056e0 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 61 6e |) ");. scan|
000056f0 66 28 22 25 4c 66 22 2c 20 26 64 29 3b 0a 20 20 |f("%Lf", &d);. |
00005700 20 20 7d 20 77 68 69 6c 65 20 28 64 3c 30 29 3b | } while (d<0);|
00005710 0a 20 20 20 20 70 6f 77 73 63 61 6c 65 20 3d 20 |. powscale = |
00005720 28 69 6e 74 29 28 64 2a 6f 6e 65 29 3b 0a 20 20 |(int)(d*one);. |
00005730 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 70 72 69 | do {. pri|
00005740 6e 74 66 28 22 50 6c 65 61 73 65 20 65 6e 74 65 |ntf("Please ente|
00005750 72 20 6c 69 6e 65 61 72 20 73 63 61 6c 69 6e 67 |r linear scaling|
00005760 20 66 61 63 74 6f 72 20 28 73 75 67 67 65 73 74 | factor (suggest|
00005770 20 31 2c 20 30 2d 31 20 70 65 72 6d 69 73 73 69 | 1, 0-1 permissi|
00005780 62 6c 65 29 20 22 29 3b 0a 20 20 20 20 20 20 73 |ble) ");. s|
00005790 63 61 6e 66 28 22 25 4c 66 22 2c 20 26 64 29 3b |canf("%Lf", &d);|
000057a0 0a 20 20 20 20 7d 20 77 68 69 6c 65 20 28 64 3c |. } while (d<|
000057b0 30 20 7c 7c 20 64 3e 31 29 3b 0a 20 20 20 20 6c |0 || d>1);. l|
000057c0 69 6e 73 63 61 6c 65 20 3d 20 28 69 6e 74 29 28 |inscale = (int)(|
000057d0 64 2a 6f 6e 65 29 3b 0a 20 20 20 20 64 6f 20 7b |d*one);. do {|
000057e0 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 50 |. printf("P|
000057f0 6c 65 61 73 65 20 65 6e 74 65 72 20 6f 72 69 65 |lease enter orie|
00005800 6e 74 61 74 69 6f 6e 20 28 73 75 67 67 65 73 74 |ntation (suggest|
00005810 20 30 2c 20 30 2e 2e 2e 33 20 70 65 72 6d 69 73 | 0, 0...3 permis|
00005820 73 69 62 6c 65 29 20 22 29 3b 0a 20 20 20 20 20 |sible) ");. |
00005830 20 73 63 61 6e 66 28 22 25 69 22 2c 20 26 74 6f | scanf("%i", &to|
00005840 72 69 65 6e 74 29 3b 0a 20 20 20 20 7d 20 77 68 |rient);. } wh|
00005850 69 6c 65 20 28 74 6f 72 69 65 6e 74 3c 30 20 7c |ile (torient<0 ||
00005860 7c 20 74 6f 72 69 65 6e 74 3e 33 29 3b 0a 20 20 || torient>3);. |
00005870 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 70 72 69 | do {. pri|
00005880 6e 74 66 28 22 50 6c 65 61 73 65 20 65 6e 74 65 |ntf("Please ente|
00005890 72 20 78 20 73 63 72 6f 6c 6c 20 76 61 6c 75 65 |r x scroll value|
000058a0 20 28 73 75 67 67 65 73 74 20 30 2c 20 30 2e 2e | (suggest 0, 0..|
000058b0 2e 35 31 31 20 70 65 72 6d 69 73 73 69 62 6c 65 |.511 permissible|
000058c0 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 61 6e |) ");. scan|
000058d0 66 28 22 25 69 22 2c 20 26 74 73 63 72 6f 6c 6c |f("%i", &tscroll|
000058e0 78 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 20 |x);. } while |
000058f0 28 74 73 63 72 6f 6c 6c 78 3c 30 20 7c 7c 20 74 |(tscrollx<0 || t|
00005900 73 63 72 6f 6c 6c 78 3e 3d 35 31 32 29 3b 0a 20 |scrollx>=512);. |
00005910 20 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 70 72 | do {. pr|
00005920 69 6e 74 66 28 22 50 6c 65 61 73 65 20 65 6e 74 |intf("Please ent|
00005930 65 72 20 79 20 73 63 72 6f 6c 6c 20 76 61 6c 75 |er y scroll valu|
00005940 65 20 28 73 75 67 67 65 73 74 20 30 2c 20 30 2e |e (suggest 0, 0.|
00005950 2e 2e 35 31 31 20 70 65 72 6d 69 73 73 69 62 6c |..511 permissibl|
00005960 65 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 61 |e) ");. sca|
00005970 6e 66 28 22 25 69 22 2c 20 26 74 73 63 72 6f 6c |nf("%i", &tscrol|
00005980 6c 79 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 |ly);. } while|
00005990 20 28 74 73 63 72 6f 6c 6c 79 3c 30 20 7c 7c 20 | (tscrolly<0 || |
000059a0 74 73 63 72 6f 6c 6c 79 3e 3d 35 31 32 29 3b 0a |tscrolly>=512);.|
000059b0 20 20 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 70 | do {. p|
000059c0 72 69 6e 74 66 28 22 50 6c 65 61 73 65 20 65 6e |rintf("Please en|
000059d0 74 65 72 20 76 69 65 77 69 6e 67 20 61 6e 67 6c |ter viewing angl|
000059e0 65 20 28 73 75 67 67 65 73 74 20 30 2e 33 2c 20 |e (suggest 0.3, |
000059f0 6f 72 20 30 2e 33 2d 30 2e 37 2c 20 74 68 6f 75 |or 0.3-0.7, thou|
00005a00 67 68 20 30 2d 31 20 70 65 72 6d 69 73 73 69 62 |gh 0-1 permissib|
00005a10 6c 65 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 |le) ");. sc|
00005a20 61 6e 66 28 22 25 4c 66 22 2c 20 26 64 29 3b 0a |anf("%Lf", &d);.|
00005a30 20 20 20 20 7d 20 77 68 69 6c 65 20 28 64 3c 30 | } while (d<0|
00005a40 20 7c 7c 20 64 3e 31 29 3b 0a 20 20 20 20 74 69 | || d>1);. ti|
00005a50 6c 74 20 3d 20 28 69 6e 74 29 28 64 2a 6f 6e 65 |lt = (int)(d*one|
00005a60 29 3b 0a 0a 20 20 20 20 69 66 20 28 21 66 6c 65 |);.. if (!fle|
00005a70 78 5f 61 6c 6c 6f 63 28 28 66 6c 65 78 5f 70 74 |x_alloc((flex_pt|
00005a80 72 29 26 64 61 74 61 2c 20 4e 2a 4e 2a 73 69 7a |r)&data, N*N*siz|
00005a90 65 6f 66 28 63 6f 6d 70 6c 65 78 29 29 29 20 77 |eof(complex))) w|
00005aa0 65 72 72 28 31 2c 20 22 49 6e 73 75 66 66 69 63 |err(1, "Insuffic|
00005ab0 69 65 6e 74 20 6d 65 6d 6f 72 79 22 29 3b 0a 20 |ient memory");. |
00005ac0 20 20 20 2f 2a 6e 6f 74 65 2c 20 66 75 72 74 68 | /*note, furth|
00005ad0 65 72 20 66 6c 65 78 20 63 61 6c 6c 73 20 6d 61 |er flex calls ma|
00005ae0 79 20 73 68 69 66 74 20 64 61 74 61 2c 20 73 6f |y shift data, so|
00005af0 20 62 65 20 63 61 72 65 66 75 6c 20 77 69 74 68 | be careful with|
00005b00 20 70 74 72 73 20 69 6e 74 6f 20 64 61 74 61 2a | ptrs into data*|
00005b10 2f 0a 20 20 20 20 62 62 63 5f 6d 6f 64 65 28 6d |/. bbc_mode(m|
00005b20 6f 6e 69 74 6f 72 74 79 70 65 3d 3d 6d 75 6c 74 |onitortype==mult|
00005b30 69 20 3f 20 32 30 20 3a 20 31 32 29 3b 0a 20 20 |i ? 20 : 12);. |
00005b40 20 20 62 62 63 5f 76 64 75 71 28 32 38 2c 37 32 | bbc_vduq(28,72|
00005b50 2c 33 31 2c 37 39 2c 30 29 3b 0a 20 20 20 20 62 |,31,79,0);. b|
00005b60 62 63 5f 70 61 6c 65 74 74 65 28 31 2c 20 31 36 |bc_palette(1, 16|
00005b70 2c 20 35 31 2c 35 31 2c 35 31 29 3b 0a 20 20 20 |, 51,51,51);. |
00005b80 20 62 62 63 5f 70 61 6c 65 74 74 65 28 32 2c 20 | bbc_palette(2, |
00005b90 31 36 2c 20 32 35 35 2c 32 35 35 2c 32 35 35 29 |16, 255,255,255)|
00005ba0 3b 0a 20 20 20 20 62 62 63 5f 70 61 6c 65 74 74 |;. bbc_palett|
00005bb0 65 28 30 2c 20 31 36 2c 20 31 39 32 2c 31 39 32 |e(0, 16, 192,192|
00005bc0 2c 31 39 32 29 3b 0a 20 20 20 20 66 6f 72 20 28 |,192);. for (|
00005bd0 69 3d 30 3b 20 69 3c 36 3b 20 69 2b 2b 29 20 7b |i=0; i<6; i++) {|
00005be0 0a 20 20 20 20 20 20 62 62 63 5f 70 61 6c 65 74 |. bbc_palet|
00005bf0 74 65 28 20 34 2b 69 2c 20 31 36 2c 20 33 32 2b |te( 4+i, 16, 32+|
00005c00 32 30 2a 69 2c 31 33 36 2b 37 2a 69 2c 33 32 2b |20*i,136+7*i,32+|
00005c10 32 30 2a 69 29 3b 0a 20 20 20 20 20 20 62 62 63 |20*i);. bbc|
00005c20 5f 70 61 6c 65 74 74 65 28 31 30 2b 69 2c 20 31 |_palette(10+i, 1|
00005c30 36 2c 20 31 32 38 2b 38 2a 69 2c 39 36 2b 31 32 |6, 128+8*i,96+12|
00005c40 2a 69 2c 36 34 2b 31 36 2a 69 29 3b 0a 20 20 20 |*i,64+16*i);. |
00005c50 20 7d 0a 20 20 20 20 6c 69 6e 65 73 65 70 20 3d | }. linesep =|
00005c60 20 6d 6f 6e 69 74 6f 72 74 79 70 65 3d 3d 6d 75 | monitortype==mu|
00005c70 6c 74 69 20 3f 20 32 30 20 3a 20 34 30 3b 0a 20 |lti ? 20 : 40;. |
00005c80 20 20 20 70 72 69 6e 74 66 28 22 5c 6e 50 6c 65 | printf("\nPle|
00005c90 61 73 65 20 20 77 61 69 74 20 2e 2e 2e 5c 6e 22 |ase wait ...\n"|
00005ca0 29 3b 0a 20 20 20 20 61 73 73 69 67 6e 5f 66 74 |);. assign_ft|
00005cb0 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 64 69 |(data, limit, di|
00005cc0 6d 2c 20 73 65 65 64 29 3b 0a 20 20 20 20 69 66 |m, seed);. if|
00005cd0 20 28 72 61 6e 67 65 3c 4e 2f 32 29 20 63 6c 69 | (range<N/2) cli|
00005ce0 70 5f 66 74 28 64 61 74 61 2c 20 6c 69 6d 69 74 |p_ft(data, limit|
00005cf0 2c 20 72 61 6e 67 65 29 3b 0a 20 20 20 20 69 6e |, range);. in|
00005d00 76 66 66 74 5f 32 64 28 64 61 74 61 2c 20 6c 69 |vfft_2d(data, li|
00005d10 6d 69 74 29 3b 0a 20 20 20 20 69 66 20 28 21 6e |mit);. if (!n|
00005d20 6f 72 6d 61 6c 69 73 65 5f 61 6e 64 5f 63 6f 6d |ormalise_and_com|
00005d30 70 61 63 74 28 26 64 61 74 61 2c 20 6c 69 6d 69 |pact(&data, limi|
00005d40 74 29 29 20 77 65 72 72 28 31 2c 20 22 43 61 6e |t)) werr(1, "Can|
00005d50 27 74 20 63 6f 6d 70 61 63 74 20 64 61 74 61 22 |'t compact data"|
00005d60 29 3b 0a 20 20 20 20 6f 72 69 65 6e 74 3d 73 63 |);. orient=sc|
00005d70 72 6f 6c 6c 78 3d 73 63 72 6f 6c 6c 79 3d 30 3b |rollx=scrolly=0;|
00005d80 0a 20 20 20 20 74 73 63 72 6f 6c 6c 78 3e 3e 3d |. tscrollx>>=|
00005d90 39 2d 6c 69 6d 69 74 3b 0a 20 20 20 20 74 73 63 |9-limit;. tsc|
00005da0 72 6f 6c 6c 79 3e 3e 3d 39 2d 6c 69 6d 69 74 3b |rolly>>=9-limit;|
00005db0 0a 20 20 20 20 69 66 20 28 74 73 63 72 6f 6c 6c |. if (tscroll|
00005dc0 78 3e 30 29 20 73 63 72 6f 6c 6c 28 64 61 74 61 |x>0) scroll(data|
00005dd0 2c 20 6c 69 6d 69 74 2c 20 72 69 67 68 74 2c 20 |, limit, right, |
00005de0 74 73 63 72 6f 6c 6c 78 29 3b 0a 20 20 20 20 69 |tscrollx);. i|
00005df0 66 20 28 74 73 63 72 6f 6c 6c 78 3c 30 29 20 73 |f (tscrollx<0) s|
00005e00 63 72 6f 6c 6c 28 64 61 74 61 2c 20 6c 69 6d 69 |croll(data, limi|
00005e10 74 2c 20 6c 65 66 74 2c 20 2d 74 73 63 72 6f 6c |t, left, -tscrol|
00005e20 6c 78 29 3b 0a 20 20 20 20 69 66 20 28 74 73 63 |lx);. if (tsc|
00005e30 72 6f 6c 6c 79 3e 30 29 20 73 63 72 6f 6c 6c 28 |rolly>0) scroll(|
00005e40 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 75 70 2c |data, limit, up,|
00005e50 20 20 20 20 74 73 63 72 6f 6c 6c 79 29 3b 0a 20 | tscrolly);. |
00005e60 20 20 20 69 66 20 28 74 73 63 72 6f 6c 6c 79 3c | if (tscrolly<|
00005e70 30 29 20 73 63 72 6f 6c 6c 28 64 61 74 61 2c 20 |0) scroll(data, |
00005e80 6c 69 6d 69 74 2c 20 64 6f 77 6e 2c 20 2d 74 73 |limit, down, -ts|
00005e90 63 72 6f 6c 6c 79 29 3b 0a 20 20 20 20 73 77 69 |crolly);. swi|
00005ea0 74 63 68 20 28 74 6f 72 69 65 6e 74 29 20 7b 0a |tch (torient) {.|
00005eb0 20 20 20 20 20 20 63 61 73 65 20 31 3a 0a 20 20 | case 1:. |
00005ec0 20 20 20 20 20 20 72 6f 74 61 74 65 28 64 61 74 | rotate(dat|
00005ed0 61 2c 20 6c 69 6d 69 74 2c 20 54 52 55 45 29 3b |a, limit, TRUE);|
00005ee0 0a 20 20 20 20 20 20 20 20 62 72 65 61 6b 3b 0a |. break;.|
00005ef0 20 20 20 20 20 20 63 61 73 65 20 32 3a 0a 20 20 | case 2:. |
00005f00 20 20 20 20 20 20 72 6f 74 61 74 65 28 64 61 74 | rotate(dat|
00005f10 61 2c 20 6c 69 6d 69 74 2c 20 54 52 55 45 29 3b |a, limit, TRUE);|
00005f20 0a 20 20 20 20 20 20 20 20 72 6f 74 61 74 65 28 |. rotate(|
00005f30 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 54 52 55 |data, limit, TRU|
00005f40 45 29 3b 0a 20 20 20 20 20 20 20 20 62 72 65 61 |E);. brea|
00005f50 6b 3b 0a 20 20 20 20 20 20 63 61 73 65 20 33 3a |k;. case 3:|
00005f60 0a 20 20 20 20 20 20 20 20 72 6f 74 61 74 65 28 |. rotate(|
00005f70 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 46 41 4c |data, limit, FAL|
00005f80 53 45 29 3b 0a 20 20 20 20 20 20 20 20 62 72 65 |SE);. bre|
00005f90 61 6b 3b 0a 20 20 20 20 7d 0a 20 20 20 20 64 6f |ak;. }. do|
00005fa0 20 7b 0a 20 20 20 20 20 20 62 62 63 5f 63 6c 67 | {. bbc_clg|
00005fb0 28 29 3b 0a 20 20 20 20 20 20 67 72 61 70 68 28 |();. graph(|
00005fc0 64 61 74 61 2c 20 6c 69 6d 69 74 29 3b 0a 20 20 |data, limit);. |
00005fd0 20 20 20 20 62 62 63 5f 76 64 75 71 28 32 38 2c | bbc_vduq(28,|
00005fe0 30 2c 33 31 2c 37 2c 30 29 3b 0a 20 20 20 20 20 |0,31,7,0);. |
00005ff0 20 70 72 69 6e 74 66 28 22 20 20 73 63 61 70 65 | printf(" scape|
00006000 5c 6e 25 37 69 5c 6e 5c 6e 20 20 20 20 72 65 73 |\n%7i\n\n res|
00006010 5c 6e 25 37 69 5c 6e 5c 6e 20 20 20 20 64 69 6d |\n%7i\n\n dim|
00006020 5c 6e 25 37 2e 32 66 5c 6e 5c 6e 20 20 72 61 6e |\n%7.2f\n\n ran|
00006030 67 65 5c 6e 25 37 2e 32 66 5c 6e 5c 6e 20 20 20 |ge\n%7.2f\n\n |
00006040 20 73 65 61 5c 6e 25 37 2e 32 66 5c 6e 5c 6e 20 | sea\n%7.2f\n\n |
00006050 70 73 63 61 6c 65 5c 6e 25 37 2e 32 66 5c 6e 5c |pscale\n%7.2f\n\|
00006060 6e 20 6c 73 63 61 6c 65 5c 6e 25 37 2e 32 66 5c |n lscale\n%7.2f\|
00006070 6e 5c 6e 20 6f 72 69 65 6e 74 5c 6e 25 37 69 5c |n\n orient\n%7i\|
00006080 6e 5c 6e 20 20 20 20 20 20 78 5c 6e 25 37 69 5c |n\n x\n%7i\|
00006090 6e 5c 6e 20 20 20 20 20 20 79 5c 6e 25 37 69 5c |n\n y\n%7i\|
000060a0 6e 5c 6e 20 20 20 74 69 6c 74 5c 6e 25 37 2e 32 |n\n tilt\n%7.2|
000060b0 66 22 2c 20 73 65 65 64 2c 20 6c 69 6d 69 74 2c |f", seed, limit,|
000060c0 20 64 69 6d 2c 20 28 69 6e 74 29 72 61 6e 67 65 | dim, (int)range|
000060d0 2f 28 64 6f 75 62 6c 65 29 28 69 6e 74 29 28 4e |/(double)(int)(N|
000060e0 2f 32 29 2c 20 28 69 6e 74 29 73 65 61 6c 65 76 |/2), (int)sealev|
000060f0 65 6c 2f 36 35 35 33 36 2e 30 2c 20 28 69 6e 74 |el/65536.0, (int|
00006100 29 70 6f 77 73 63 61 6c 65 2f 36 35 35 33 36 2e |)powscale/65536.|
00006110 30 2c 20 28 69 6e 74 29 6c 69 6e 73 63 61 6c 65 |0, (int)linscale|
00006120 2f 36 35 35 33 36 2e 30 2c 20 6f 72 69 65 6e 74 |/65536.0, orient|
00006130 2c 20 73 63 72 6f 6c 6c 78 2c 20 73 63 72 6f 6c |, scrollx, scrol|
00006140 6c 79 2c 20 28 69 6e 74 29 74 69 6c 74 2f 36 35 |ly, (int)tilt/65|
00006150 35 33 36 2e 30 29 3b 0a 20 20 20 20 20 20 62 62 |536.0);. bb|
00006160 63 5f 76 64 75 71 28 32 38 2c 37 32 2c 33 31 2c |c_vduq(28,72,31,|
00006170 37 39 2c 30 29 3b 0a 20 20 20 20 20 20 62 62 63 |79,0);. bbc|
00006180 5f 63 6c 73 28 29 3b 0a 20 20 20 20 20 20 70 72 |_cls();. pr|
00006190 69 6e 74 66 28 22 5c 6e 5c 6e 5c 6e 5c 6e 5c 6e |intf("\n\n\n\n\n|
000061a0 5c 6e 5c 6e 5c 6e 50 72 65 73 73 20 61 20 6b 65 |\n\n\nPress a ke|
000061b0 79 20 2e 2e 2e 5c 6e 5c 6e 22 29 3b 0a 20 20 20 |y ...\n\n");. |
000061c0 20 20 20 70 72 69 6e 74 66 28 22 5b 20 27 3c 27 | printf("[ '<'|
000061d0 5c 6e 20 20 6f 72 5c 6e 20 20 27 3e 27 5c 6e 20 |\n or\n '>'\n |
000061e0 20 74 6f 5c 6e 72 6f 74 61 74 65 5c 6e 20 73 63 | to\nrotate\n sc|
000061f0 61 70 65 5c 6e 5c 6e 63 75 72 73 6f 72 73 5c 6e |ape\n\ncursors\n|
00006200 20 28 2b 5c 78 38 62 29 5c 6e 20 20 74 6f 5c 6e | (+\x8b)\n to\n|
00006210 73 63 72 6f 6c 6c 5c 6e 20 73 63 61 70 65 5c 6e |scroll\n scape\n|
00006220 5c 6e 20 63 6f 70 79 5c 6e 20 20 74 6f 5c 6e 20 |\n copy\n to\n |
00006230 73 61 76 65 5c 6e 20 69 6d 61 67 65 5c 6e 5c 6e |save\n image\n\n|
00006240 5c 78 38 62 20 63 6f 70 79 5c 6e 20 20 74 6f 5c |\x8b copy\n to\|
00006250 6e 20 73 61 76 65 5c 6e 63 6f 6c 6f 75 72 5c 6e |n save\ncolour\n|
00006260 20 69 6d 61 67 65 5c 6e 20 20 20 20 20 20 20 5d | image\n ]|
00006270 5c 6e 22 29 3b 0a 20 20 20 20 20 20 64 6f 20 7b |\n");. do {|
00006280 0a 20 20 20 20 20 20 20 20 6f 73 5f 63 6c 69 28 |. os_cli(|
00006290 22 66 78 20 31 35 20 31 22 29 3b 0a 20 20 20 20 |"fx 15 1");. |
000062a0 20 20 20 20 62 62 63 5f 63 75 72 73 6f 72 28 30 | bbc_cursor(0|
000062b0 29 3b 0a 20 20 20 20 20 20 20 20 6f 73 5f 73 77 |);. os_sw|
000062c0 69 32 28 4f 53 5f 42 79 74 65 2c 34 2c 31 29 3b |i2(OS_Byte,4,1);|
000062d0 0a 20 20 20 20 20 20 20 20 6b 65 79 20 3d 20 62 |. key = b|
000062e0 62 63 5f 67 65 74 28 29 3b 0a 20 20 20 20 20 20 |bc_get();. |
000062f0 20 20 6f 73 5f 73 77 69 32 28 4f 53 5f 42 79 74 | os_swi2(OS_Byt|
00006300 65 2c 34 2c 30 29 3b 0a 20 20 20 20 20 20 20 20 |e,4,0);. |
00006310 69 66 20 28 6b 65 79 3d 3d 31 33 35 29 20 7b 0a |if (key==135) {.|
00006320 20 20 20 20 20 20 20 20 20 20 73 63 72 65 65 6e | screen|
00006330 2b 3d 31 3b 0a 20 20 20 20 20 20 20 20 20 20 73 |+=1;. s|
00006340 70 72 69 6e 74 66 28 62 75 66 2c 20 22 3c 53 63 |printf(buf, "<Sc|
00006350 61 70 65 32 24 44 69 72 3e 2e 5e 2e 49 6d 61 67 |ape2$Dir>.^.Imag|
00006360 65 25 30 32 69 22 2c 20 73 63 72 65 65 6e 29 3b |e%02i", screen);|
00006370 0a 20 20 20 20 20 20 20 20 20 20 69 66 20 28 61 |. if (a|
00006380 6b 62 64 5f 70 6f 6c 6c 73 68 28 29 29 20 7b 0a |kbd_pollsh()) {.|
00006390 20 20 20 20 20 20 20 20 20 20 20 20 64 75 6d 70 | dump|
000063a0 5f 69 6d 61 67 65 33 32 28 64 61 74 61 2c 20 6c |_image32(data, l|
000063b0 69 6d 69 74 2c 20 69 6d 33 32 5f 79 6c 2c 20 69 |imit, im32_yl, i|
000063c0 6d 33 32 5f 68 65 69 67 68 74 2c 20 62 75 66 29 |m32_height, buf)|
000063d0 3b 0a 20 20 20 20 20 20 20 20 20 20 7d 0a 20 20 |;. }. |
000063e0 20 20 20 20 20 20 20 20 65 6c 73 65 20 7b 0a 20 | else {. |
000063f0 20 20 20 20 20 20 20 20 20 20 20 77 69 6d 70 74 | wimpt|
00006400 5f 63 6f 6d 70 6c 61 69 6e 28 73 70 72 69 74 65 |_complain(sprite|
00006410 5f 73 63 72 65 65 6e 73 61 76 65 28 62 75 66 2c |_screensave(buf,|
00006420 20 54 52 55 45 29 29 3b 0a 20 20 20 20 20 20 20 | TRUE));. |
00006430 20 20 20 7d 0a 20 20 20 20 20 20 20 20 7d 0a 20 | }. }. |
00006440 20 20 20 20 20 7d 20 77 68 69 6c 65 20 28 6b 65 | } while (ke|
00006450 79 3d 3d 31 33 35 29 3b 0a 20 20 20 20 20 20 69 |y==135);. i|
00006460 66 20 28 6b 65 79 3d 3d 34 34 20 7c 7c 20 6b 65 |f (key==44 || ke|
00006470 79 3d 3d 34 36 29 20 72 6f 74 61 74 65 28 64 61 |y==46) rotate(da|
00006480 74 61 2c 20 6c 69 6d 69 74 2c 20 6b 65 79 3d 3d |ta, limit, key==|
00006490 34 36 29 3b 0a 20 20 20 20 20 20 69 66 20 28 6b |46);. if (k|
000064a0 65 79 3d 3d 31 33 36 20 7c 7c 20 6b 65 79 3d 3d |ey==136 || key==|
000064b0 31 33 37 20 7c 7c 20 6b 65 79 3d 3d 31 33 38 20 |137 || key==138 |
000064c0 7c 7c 20 6b 65 79 3d 3d 31 33 39 29 20 7b 0a 20 ||| key==139) {. |
000064d0 20 20 20 20 20 20 20 69 3d 6c 69 6d 69 74 2d 32 | i=limit-2|
000064e0 3b 0a 20 20 20 20 20 20 20 20 69 66 20 28 61 6b |;. if (ak|
000064f0 62 64 5f 70 6f 6c 6c 73 68 28 29 29 20 69 2d 3d |bd_pollsh()) i-=|
00006500 32 3b 0a 20 20 20 20 20 20 20 20 69 66 20 28 69 |2;. if (i|
00006510 3c 30 29 20 69 3d 30 3b 0a 20 20 20 20 20 20 20 |<0) i=0;. |
00006520 20 73 63 72 6f 6c 6c 28 64 61 74 61 2c 20 6c 69 | scroll(data, li|
00006530 6d 69 74 2c 20 6b 65 79 2c 20 31 3c 3c 69 29 3b |mit, key, 1<<i);|
00006540 0a 20 20 20 20 20 20 7d 0a 20 20 20 20 7d 20 77 |. }. } w|
00006550 68 69 6c 65 20 28 6b 65 79 3d 3d 34 34 20 7c 7c |hile (key==44 |||
00006560 20 6b 65 79 3d 3d 34 36 20 7c 7c 20 6b 65 79 3d | key==46 || key=|
00006570 3d 31 33 36 20 7c 7c 20 6b 65 79 3d 3d 31 33 37 |=136 || key==137|
00006580 20 7c 7c 20 6b 65 79 3d 3d 31 33 38 20 7c 7c 20 | || key==138 || |
00006590 6b 65 79 3d 3d 31 33 39 29 3b 0a 20 20 20 20 66 |key==139);. f|
000065a0 6c 65 78 5f 66 72 65 65 28 28 66 6c 65 78 5f 70 |lex_free((flex_p|
000065b0 74 72 29 26 64 61 74 61 29 3b 0a 20 20 20 20 62 |tr)&data);. b|
000065c0 62 63 5f 63 6c 73 28 29 3b 0a 0a 20 7d 0a 0a 20 |bc_cls();.. }.. |
000065d0 20 72 65 74 75 72 6e 20 30 3b 0a 7d 0a | return 0;.}.|
000065dd 
.