Home » Archimedes archive » Archimedes World » AW-1996-03-Disc 2.adf » !AcornAns_AcornAns » !Scape1/c/main
!Scape1/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-03-Disc 2.adf » !AcornAns_AcornAns |
| Filename: | !Scape1/c/main |
| Read OK: | ✔ |
| File size: | 3AF0 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;
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);
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;
}
/*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;
}
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;
}
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 (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);
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(void)
{
int seed, i;
char key;
double dim=2.2;
double d;
char buf[1024];
int torient, tscrollx, tscrolly;
monitortype = read_monitor_type();
flex_init();
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 %i, 0...%i permissible) ",N/2,N/2);
scanf("%i", &range);
} while (range<0 || 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);
tilt = 20000; /*0-65536 permissible*/
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("\n\n scape\n%7i\n\n res\n%7i\n\n dim\n%7.2f\n\n range\n%7i\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", seed, limit, dim, range, (int)sealevel/65536.0, (int)powscale/65536.0, (int)linscale/65536.0, orient, scrollx, scrolly);
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");
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, "<Scape1$Dir>.^.Image%03i", screen);
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 65 78 74 65 72 6e 20 69 6e 74 20 |=0;..extern int |
000002b0 72 62 62 63 69 6e 63 28 69 6e 74 20 72 2c 20 69 |rbbcinc(int r, i|
000002c0 6e 74 20 6b 29 3b 0a 65 78 74 65 72 6e 20 69 6e |nt k);.extern in|
000002d0 74 20 64 69 76 5f 66 72 61 63 31 36 28 69 6e 74 |t div_frac16(int|
000002e0 20 6e 75 6d 62 65 72 2c 20 69 6e 74 20 64 69 76 | number, int div|
000002f0 69 73 6f 72 29 3b 0a 65 78 74 65 72 6e 20 69 6e |isor);.extern in|
00000300 74 20 6d 75 6c 5f 66 72 61 63 31 36 28 69 6e 74 |t mul_frac16(int|
00000310 20 78 2c 20 69 6e 74 20 61 29 3b 0a 65 78 74 65 | x, int a);.exte|
00000320 72 6e 20 69 6e 74 20 6d 75 6c 5f 66 72 61 63 31 |rn int mul_frac1|
00000330 36 63 28 69 6e 74 20 78 2c 20 69 6e 74 20 61 29 |6c(int x, int a)|
00000340 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 73 71 72 |;.extern int sqr|
00000350 74 5f 66 72 61 63 31 36 28 75 6e 73 69 67 6e 65 |t_frac16(unsigne|
00000360 64 20 69 6e 74 20 78 29 3b 0a 65 78 74 65 72 6e |d int x);.extern|
00000370 20 69 6e 74 20 67 61 75 73 73 31 36 28 76 6f 69 | int gauss16(voi|
00000380 64 29 3b 0a 65 78 74 65 72 6e 20 76 6f 69 64 20 |d);.extern void |
00000390 73 67 61 75 73 73 31 36 28 69 6e 74 20 73 65 65 |sgauss16(int see|
000003a0 64 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 63 |d);.extern int c|
000003b0 6f 73 31 36 28 69 6e 74 20 61 29 3b 0a 65 78 74 |os16(int a);.ext|
000003c0 65 72 6e 20 69 6e 74 20 73 69 6e 31 36 28 69 6e |ern int sin16(in|
000003d0 74 20 61 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 |t a);.extern int|
000003e0 20 65 78 70 31 36 28 69 6e 74 20 61 29 3b 0a 65 | exp16(int a);.e|
000003f0 78 74 65 72 6e 20 69 6e 74 20 6c 6e 31 36 28 69 |xtern int ln16(i|
00000400 6e 74 20 61 29 3b 0a 65 78 74 65 72 6e 20 69 6e |nt a);.extern in|
00000410 74 20 70 6f 77 31 36 28 69 6e 74 20 61 2c 20 69 |t pow16(int a, i|
00000420 6e 74 20 62 29 3b 0a 0a 42 4f 4f 4c 20 69 6e 76 |nt b);..BOOL inv|
00000430 66 66 74 5f 32 64 28 63 6f 6d 70 6c 65 78 20 2a |fft_2d(complex *|
00000440 64 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 29 |data, int limit)|
00000450 0a 7b 0a 20 20 69 6e 74 20 69 2c 20 6a 3b 0a 20 |.{. int i, j;. |
00000460 20 69 6e 74 20 6e 20 3d 20 31 3c 3c 6c 69 6d 69 | int n = 1<<limi|
00000470 74 3b 0a 20 20 63 6f 6d 70 6c 65 78 20 74 31 2c |t;. complex t1,|
00000480 20 2a 64 31 2c 20 2a 64 31 65 2c 20 2a 64 32 3b | *d1, *d1e, *d2;|
00000490 0a 20 20 69 6e 74 20 6b 2c 20 6c 3b 0a 20 20 69 |. int k, l;. i|
000004a0 6e 74 20 6c 6f 6f 70 3b 0a 20 20 69 6e 74 20 72 |nt loop;. int r|
000004b0 6f 74 61 6e 67 6c 65 2c 20 61 6e 67 6c 65 3b 0a |otangle, angle;.|
000004c0 20 20 69 6e 74 20 63 6f 73 2c 20 73 69 6e 3b 0a | int cos, sin;.|
000004d0 20 20 69 66 20 28 6c 69 6d 69 74 3c 33 20 7c 7c | if (limit<3 |||
000004e0 20 6c 69 6d 69 74 3e 33 32 29 20 72 65 74 75 72 | limit>32) retur|
000004f0 6e 20 46 41 4c 53 45 3b 0a 0a 20 20 66 6f 72 20 |n FALSE;.. for |
00000500 28 69 3d 6a 3d 30 3b 20 6a 3c 6e 3b 20 69 3d 72 |(i=j=0; j<n; i=r|
00000510 62 62 63 69 6e 63 28 69 2c 20 6c 69 6d 69 74 29 |bbcinc(i, limit)|
00000520 2c 20 6a 2b 2b 29 20 69 66 20 28 69 3e 6a 29 20 |, j++) if (i>j) |
00000530 7b 0a 20 20 20 20 66 6f 72 20 28 64 31 3d 64 61 |{. for (d1=da|
00000540 74 61 2b 28 69 3c 3c 6c 69 6d 69 74 29 2c 20 64 |ta+(i<<limit), d|
00000550 32 3d 64 61 74 61 2b 28 6a 3c 3c 6c 69 6d 69 74 |2=data+(j<<limit|
00000560 29 2c 20 64 31 65 3d 64 31 2b 6e 3b 20 64 31 3c |), d1e=d1+n; d1<|
00000570 64 31 65 3b 29 20 7b 0a 20 20 20 20 20 20 74 31 |d1e;) {. t1|
00000580 20 3d 20 2a 64 31 3b 0a 20 20 20 20 20 20 2a 64 | = *d1;. *d|
00000590 31 2b 2b 20 3d 20 2a 64 32 3b 0a 20 20 20 20 20 |1++ = *d2;. |
000005a0 20 2a 64 32 2b 2b 20 3d 20 74 31 3b 0a 20 20 20 | *d2++ = t1;. |
000005b0 20 7d 0a 20 20 7d 0a 20 20 72 6f 74 61 6e 67 6c | }. }. rotangl|
000005c0 65 20 3d 20 2d 32 2a 6f 6e 65 3b 0a 20 20 66 6f |e = -2*one;. fo|
000005d0 72 20 28 6b 3d 6c 6f 6f 70 3d 31 3b 20 6c 6f 6f |r (k=loop=1; loo|
000005e0 70 3c 3d 6c 69 6d 69 74 3b 20 6c 6f 6f 70 2b 2b |p<=limit; loop++|
000005f0 2c 20 6b 3c 3c 3d 31 2c 20 72 6f 74 61 6e 67 6c |, k<<=1, rotangl|
00000600 65 2f 3d 32 29 20 7b 0a 20 20 20 20 66 6f 72 20 |e/=2) {. for |
00000610 28 61 6e 67 6c 65 3d 30 2c 20 73 69 6e 3d 30 2c |(angle=0, sin=0,|
00000620 20 63 6f 73 3d 6f 6e 65 2c 20 69 3d 31 3b 20 69 | cos=one, i=1; i|
00000630 3c 3d 6b 3b 20 69 2b 2b 2c 20 61 6e 67 6c 65 2b |<=k; i++, angle+|
00000640 3d 72 6f 74 61 6e 67 6c 65 2c 20 73 69 6e 3d 73 |=rotangle, sin=s|
00000650 69 6e 31 36 28 61 6e 67 6c 65 29 2c 20 63 6f 73 |in16(angle), cos|
00000660 3d 63 6f 73 31 36 28 61 6e 67 6c 65 29 29 20 7b |=cos16(angle)) {|
00000670 0a 20 20 20 20 20 20 66 6f 72 20 28 6a 3d 69 2d |. for (j=i-|
00000680 31 3b 20 6a 3c 6e 3b 20 6a 2b 3d 6b 3c 3c 31 29 |1; j<n; j+=k<<1)|
00000690 20 7b 0a 20 20 20 20 20 20 20 20 6c 20 3d 20 6a | {. l = j|
000006a0 2b 6b 3b 0a 20 20 20 20 20 20 20 20 66 6f 72 20 |+k;. for |
000006b0 28 64 31 3d 64 61 74 61 2b 28 6c 3c 3c 6c 69 6d |(d1=data+(l<<lim|
000006c0 69 74 29 2c 20 64 32 3d 64 61 74 61 2b 28 6a 3c |it), d2=data+(j<|
000006d0 3c 6c 69 6d 69 74 29 2c 20 64 31 65 3d 64 31 2b |<limit), d1e=d1+|
000006e0 6e 3b 20 64 31 3c 64 31 65 3b 29 20 7b 0a 20 20 |n; d1<d1e;) {. |
000006f0 20 20 20 20 20 20 20 20 74 31 2e 72 65 20 3d 20 | t1.re = |
00000700 28 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d |(mul_frac16c(d1-|
00000710 3e 72 65 2c 20 63 6f 73 29 2d 6d 75 6c 5f 66 72 |>re, cos)-mul_fr|
00000720 61 63 31 36 63 28 64 31 2d 3e 69 6d 2c 20 73 69 |ac16c(d1->im, si|
00000730 6e 29 29 3b 0a 20 20 20 20 20 20 20 20 20 20 74 |n));. t|
00000740 31 2e 69 6d 20 3d 20 28 6d 75 6c 5f 66 72 61 63 |1.im = (mul_frac|
00000750 31 36 63 28 64 31 2d 3e 72 65 2c 20 73 69 6e 29 |16c(d1->re, sin)|
00000760 2b 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d |+mul_frac16c(d1-|
00000770 3e 69 6d 2c 20 63 6f 73 29 29 3b 0a 20 20 20 20 |>im, cos));. |
00000780 20 20 20 20 20 20 64 31 2d 3e 72 65 20 20 20 20 | d1->re |
00000790 3d 20 64 32 2d 3e 72 65 2d 74 31 2e 72 65 3b 0a |= d2->re-t1.re;.|
000007a0 20 20 20 20 20 20 20 20 20 20 64 31 2b 2b 2d 3e | d1++->|
000007b0 69 6d 20 20 3d 20 64 32 2d 3e 69 6d 2d 74 31 2e |im = d2->im-t1.|
000007c0 69 6d 3b 0a 20 20 20 20 20 20 20 20 20 20 64 32 |im;. d2|
000007d0 2d 3e 72 65 20 20 20 2b 3d 20 74 31 2e 72 65 3b |->re += t1.re;|
000007e0 0a 20 20 20 20 20 20 20 20 20 20 64 32 2b 2b 2d |. d2++-|
000007f0 3e 69 6d 20 2b 3d 20 74 31 2e 69 6d 3b 0a 20 20 |>im += t1.im;. |
00000800 20 20 20 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a | }. }.|
00000810 20 20 20 20 7d 0a 20 20 7d 0a 0a 20 20 66 6f 72 | }. }.. for|
00000820 20 28 69 3d 6a 3d 30 3b 20 6a 3c 6e 3b 20 69 3d | (i=j=0; j<n; i=|
00000830 72 62 62 63 69 6e 63 28 69 2c 20 6c 69 6d 69 74 |rbbcinc(i, limit|
00000840 29 2c 20 6a 2b 2b 29 20 69 66 20 28 69 3e 6a 29 |), j++) if (i>j)|
00000850 20 7b 0a 20 20 20 20 66 6f 72 20 28 64 31 3d 64 | {. for (d1=d|
00000860 61 74 61 2b 69 2c 20 64 32 3d 64 61 74 61 2b 6a |ata+i, d2=data+j|
00000870 2c 20 64 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 |, d1e=d1+n*n; d1|
00000880 3c 64 31 65 3b 20 64 31 2b 3d 6e 2c 20 64 32 2b |<d1e; d1+=n, d2+|
00000890 3d 6e 29 20 7b 0a 20 20 20 20 20 20 74 31 20 3d |=n) {. t1 =|
000008a0 20 2a 64 31 3b 0a 20 20 20 20 20 20 2a 64 31 20 | *d1;. *d1 |
000008b0 3d 20 2a 64 32 3b 0a 20 20 20 20 20 20 2a 64 32 |= *d2;. *d2|
000008c0 20 3d 20 74 31 3b 0a 20 20 20 20 7d 0a 20 20 7d | = t1;. }. }|
000008d0 0a 20 20 72 6f 74 61 6e 67 6c 65 20 3d 20 2d 32 |. rotangle = -2|
000008e0 2a 6f 6e 65 3b 0a 20 20 66 6f 72 20 28 6b 3d 6c |*one;. for (k=l|
000008f0 6f 6f 70 3d 31 3b 20 6c 6f 6f 70 3c 3d 6c 69 6d |oop=1; loop<=lim|
00000900 69 74 3b 20 6c 6f 6f 70 2b 2b 2c 20 6b 3c 3c 3d |it; loop++, k<<=|
00000910 31 2c 20 72 6f 74 61 6e 67 6c 65 2f 3d 32 29 20 |1, rotangle/=2) |
00000920 7b 0a 20 20 20 20 66 6f 72 20 28 61 6e 67 6c 65 |{. for (angle|
00000930 3d 30 2c 20 73 69 6e 3d 30 2c 20 63 6f 73 3d 6f |=0, sin=0, cos=o|
00000940 6e 65 2c 20 69 3d 31 3b 20 69 3c 3d 6b 3b 20 69 |ne, i=1; i<=k; i|
00000950 2b 2b 2c 20 61 6e 67 6c 65 2b 3d 72 6f 74 61 6e |++, angle+=rotan|
00000960 67 6c 65 2c 20 73 69 6e 3d 73 69 6e 31 36 28 61 |gle, sin=sin16(a|
00000970 6e 67 6c 65 29 2c 20 63 6f 73 3d 63 6f 73 31 36 |ngle), cos=cos16|
00000980 28 61 6e 67 6c 65 29 29 20 7b 0a 20 20 20 20 20 |(angle)) {. |
00000990 20 66 6f 72 20 28 6a 3d 69 2d 31 3b 20 6a 3c 6e | for (j=i-1; j<n|
000009a0 3b 20 6a 2b 3d 6b 3c 3c 31 29 20 7b 0a 20 20 20 |; j+=k<<1) {. |
000009b0 20 20 20 20 20 6c 20 3d 20 6a 2b 6b 3b 0a 20 20 | l = j+k;. |
000009c0 20 20 20 20 20 20 66 6f 72 20 28 64 31 3d 64 61 | for (d1=da|
000009d0 74 61 2b 6c 2c 20 64 32 3d 64 61 74 61 2b 6a 2c |ta+l, d2=data+j,|
000009e0 20 64 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 3c | d1e=d1+n*n; d1<|
000009f0 64 31 65 3b 20 64 31 2b 3d 6e 2c 20 64 32 2b 3d |d1e; d1+=n, d2+=|
00000a00 6e 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 74 |n) {. t|
00000a10 31 2e 72 65 20 3d 20 28 6d 75 6c 5f 66 72 61 63 |1.re = (mul_frac|
00000a20 31 36 63 28 64 31 2d 3e 72 65 2c 20 63 6f 73 29 |16c(d1->re, cos)|
00000a30 2d 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d |-mul_frac16c(d1-|
00000a40 3e 69 6d 2c 20 73 69 6e 29 29 3b 0a 20 20 20 20 |>im, sin));. |
00000a50 20 20 20 20 20 20 74 31 2e 69 6d 20 3d 20 28 6d | t1.im = (m|
00000a60 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e 72 |ul_frac16c(d1->r|
00000a70 65 2c 20 73 69 6e 29 2b 6d 75 6c 5f 66 72 61 63 |e, sin)+mul_frac|
00000a80 31 36 63 28 64 31 2d 3e 69 6d 2c 20 63 6f 73 29 |16c(d1->im, cos)|
00000a90 29 3b 0a 20 20 20 20 20 20 20 20 20 20 64 31 2d |);. d1-|
00000aa0 3e 72 65 20 20 3d 20 64 32 2d 3e 72 65 2d 74 31 |>re = d2->re-t1|
00000ab0 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 20 20 64 |.re;. d|
00000ac0 31 2d 3e 69 6d 20 20 3d 20 64 32 2d 3e 69 6d 2d |1->im = d2->im-|
00000ad0 74 31 2e 69 6d 3b 0a 20 20 20 20 20 20 20 20 20 |t1.im;. |
00000ae0 20 64 32 2d 3e 72 65 20 2b 3d 20 74 31 2e 72 65 | d2->re += t1.re|
00000af0 3b 0a 20 20 20 20 20 20 20 20 20 20 64 32 2d 3e |;. d2->|
00000b00 69 6d 20 2b 3d 20 74 31 2e 69 6d 3b 0a 20 20 20 |im += t1.im;. |
00000b10 20 20 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a 20 | }. }. |
00000b20 20 20 20 7d 0a 20 20 7d 0a 0a 20 20 72 65 74 75 | }. }.. retu|
00000b30 72 6e 20 54 52 55 45 3b 0a 7d 0a 0a 2f 2a 54 68 |rn TRUE;.}../*Th|
00000b40 65 20 66 6f 6c 6c 6f 77 69 6e 67 20 61 73 73 69 |e following assi|
00000b50 67 6e 73 20 76 61 6c 75 65 73 20 74 6f 20 32 64 |gns values to 2d|
00000b60 20 64 61 74 61 20 61 72 72 61 79 20 77 69 74 68 | data array with|
00000b70 20 73 75 69 74 61 62 6c 65 20 70 72 6f 70 65 72 | suitable proper|
00000b80 74 69 65 73 2e 2a 2f 0a 2f 2a 49 74 20 64 6f 65 |ties.*/./*It doe|
00000b90 73 20 68 6f 77 65 76 65 72 20 64 6f 20 73 6f 20 |s however do so |
00000ba0 62 79 20 73 63 61 6e 6e 69 6e 67 20 74 68 65 20 |by scanning the |
00000bb0 61 72 72 61 79 20 69 6e 20 61 20 6e 6f 6e 2d 6f |array in a non-o|
00000bc0 62 76 69 6f 75 73 20 6f 72 64 65 72 2e 0a 20 20 |bvious order.. |
00000bd0 54 68 69 73 20 68 61 73 20 62 65 65 6e 20 63 68 |This has been ch|
00000be0 6f 73 65 6e 20 73 6f 20 74 68 61 74 20 69 66 20 |osen so that if |
00000bf0 6e 20 69 73 20 69 6e 63 72 65 61 73 65 64 20 62 |n is increased b|
00000c00 79 20 64 6f 75 62 6c 69 6e 67 2c 20 74 68 65 6e |y doubling, then|
00000c10 20 74 68 6f 73 65 20 72 61 6e 64 6f 6d 20 76 61 | those random va|
00000c20 6c 75 65 73 20 75 73 65 64 20 66 6f 72 20 74 68 |lues used for th|
00000c30 65 0a 20 20 6c 61 73 74 20 67 72 69 64 20 77 69 |e. last grid wi|
00000c40 6c 6c 20 62 65 20 6e 6f 77 20 77 72 69 74 74 65 |ll be now writte|
00000c50 6e 20 74 6f 20 74 68 65 20 6c 6f 77 20 66 72 65 |n to the low fre|
00000c60 71 75 65 6e 63 79 20 71 75 61 64 72 61 6e 74 20 |quency quadrant |
00000c70 6f 66 20 74 68 65 20 6e 65 77 20 67 72 69 64 20 |of the new grid |
00000c80 28 69 66 20 73 61 6d 65 20 73 65 65 64 20 75 73 |(if same seed us|
00000c90 65 64 29 2c 0a 20 20 73 75 69 74 61 62 6c 79 20 |ed),. suitably |
00000ca0 72 65 20 66 69 6c 74 65 72 65 64 2e 20 54 68 69 |re filtered. Thi|
00000cb0 73 20 73 68 6f 75 6c 64 20 68 61 76 65 20 65 66 |s should have ef|
00000cc0 66 65 63 74 20 6f 66 20 67 69 76 69 6e 67 20 61 |fect of giving a|
00000cd0 20 64 65 67 72 65 65 20 6f 66 20 63 6f 6e 74 69 | degree of conti|
00000ce0 6e 75 69 74 79 20 74 6f 20 74 68 65 20 63 68 61 |nuity to the cha|
00000cf0 6e 67 65 20 69 6e 20 6e 2c 0a 20 20 73 69 6e 63 |nge in n,. sinc|
00000d00 65 20 66 6f 72 6d 65 72 20 66 72 65 71 75 65 6e |e former frequen|
00000d10 63 69 65 73 20 61 72 65 20 63 61 72 72 69 65 64 |cies are carried|
00000d20 20 6f 76 65 72 20 74 6f 20 74 68 65 20 6c 6f 77 | over to the low|
00000d30 20 66 72 65 71 75 65 6e 63 69 65 73 20 6f 66 20 | frequencies of |
00000d40 74 68 65 20 68 69 67 68 65 72 20 72 65 73 20 67 |the higher res g|
00000d50 72 69 64 2e 0a 2a 2f 0a 42 4f 4f 4c 20 61 73 73 |rid..*/.BOOL ass|
00000d60 69 67 6e 5f 66 74 28 63 6f 6d 70 6c 65 78 20 2a |ign_ft(complex *|
00000d70 64 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 2c |data, int limit,|
00000d80 20 64 6f 75 62 6c 65 20 64 69 6d 2c 20 69 6e 74 | double dim, int|
00000d90 20 73 65 65 64 29 0a 7b 0a 20 20 69 6e 74 20 6e | seed).{. int n|
00000da0 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 |=1<<limit;. int|
00000db0 20 69 2c 20 69 63 2c 20 6a 2c 20 69 32 2c 20 6a | i, ic, j, i2, j|
00000dc0 32 2c 20 6b 2c 20 6c 3b 0a 20 20 63 6f 6d 70 6c |2, k, l;. compl|
00000dd0 65 78 20 2a 64 31 2c 20 2a 64 32 2c 20 2a 64 3b |ex *d1, *d2, *d;|
00000de0 0a 20 20 69 6e 74 20 70 68 61 73 65 2c 20 72 61 |. int phase, ra|
00000df0 64 3b 0a 20 20 69 6e 74 20 69 6d 2c 20 72 65 3b |d;. int im, re;|
00000e00 0a 20 20 69 6e 74 20 6d 62 70 6f 62 74 3b 0a 0a |. int mbpobt;..|
00000e10 20 20 73 72 61 6e 64 28 2d 73 65 65 64 29 3b 0a | srand(-seed);.|
00000e20 20 20 73 67 61 75 73 73 31 36 28 2d 73 65 65 64 | sgauss16(-seed|
00000e30 29 3b 0a 0a 20 20 6d 62 70 6f 62 74 20 3d 20 28 |);.. mbpobt = (|
00000e40 69 6e 74 29 28 2d 28 34 2d 64 69 6d 29 2a 6f 6e |int)(-(4-dim)*on|
00000e50 65 29 3b 0a 20 20 64 61 74 61 2d 3e 72 65 20 3d |e);. data->re =|
00000e60 20 64 61 74 61 2d 3e 69 6d 20 3d 20 30 3b 0a 20 | data->im = 0;. |
00000e70 20 66 6f 72 20 28 69 3d 31 3b 20 69 3c 3d 6e 2f | for (i=1; i<=n/|
00000e80 32 3b 20 69 2b 2b 29 20 7b 0a 20 20 20 20 6a 32 |2; i++) {. j2|
00000e90 20 3d 20 69 32 20 3d 20 69 2a 69 3b 0a 20 20 20 | = i2 = i*i;. |
00000ea0 20 6b 20 20 3d 20 6e 2d 69 3b 0a 20 20 20 20 64 | k = n-i;. d|
00000eb0 31 20 3d 20 64 61 74 61 2b 69 2a 6e 3b 0a 20 20 |1 = data+i*n;. |
00000ec0 20 20 64 32 20 3d 20 64 61 74 61 2b 6b 2a 6e 3b | d2 = data+k*n;|
00000ed0 0a 20 20 20 20 66 6f 72 20 28 6a 3d 30 3b 20 6a |. for (j=0; j|
00000ee0 3c 3d 69 3b 20 6a 2b 2b 29 20 7b 0a 20 20 20 20 |<=i; j++) {. |
00000ef0 20 20 6c 20 3d 20 6a 3d 3d 30 20 3f 20 30 20 3a | l = j==0 ? 0 :|
00000f00 20 6e 2d 6a 3b 0a 20 20 20 20 20 20 70 68 61 73 | n-j;. phas|
00000f10 65 20 3d 20 72 61 6e 64 28 29 3e 3e 31 33 3b 0a |e = rand()>>13;.|
00000f20 20 20 20 20 20 20 72 61 64 20 3d 20 6d 75 6c 5f | rad = mul_|
00000f30 66 72 61 63 31 36 28 70 6f 77 31 36 28 73 71 72 |frac16(pow16(sqr|
00000f40 74 5f 66 72 61 63 31 36 28 28 69 32 2b 6a 2a 6a |t_frac16((i2+j*j|
00000f50 29 2a 28 6f 6e 65 2f 34 29 29 2a 32 2c 20 6d 62 |)*(one/4))*2, mb|
00000f60 70 6f 62 74 29 2c 20 67 61 75 73 73 31 36 28 29 |pobt), gauss16()|
00000f70 29 3b 0a 20 20 20 20 20 20 72 65 20 3d 20 6d 75 |);. re = mu|
00000f80 6c 5f 66 72 61 63 31 36 63 28 72 61 64 2c 20 63 |l_frac16c(rad, c|
00000f90 6f 73 31 36 28 70 68 61 73 65 29 29 3b 0a 20 20 |os16(phase));. |
00000fa0 20 20 20 20 69 6d 20 3d 20 6d 75 6c 5f 66 72 61 | im = mul_fra|
00000fb0 63 31 36 63 28 72 61 64 2c 20 73 69 6e 31 36 28 |c16c(rad, sin16(|
00000fc0 70 68 61 73 65 29 29 3b 0a 20 20 20 20 20 20 64 |phase));. d|
00000fd0 20 3d 20 64 31 2b 6a 3b 0a 20 20 20 20 20 20 64 | = d1+j;. d|
00000fe0 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 20 20 20 20 |->re = re;. |
00000ff0 20 64 2d 3e 69 6d 20 3d 20 69 6d 3b 0a 20 20 20 | d->im = im;. |
00001000 20 20 20 64 20 3d 20 64 32 2b 6c 3b 0a 20 20 20 | d = d2+l;. |
00001010 20 20 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 | d->re = re;. |
00001020 20 20 20 20 20 64 2d 3e 69 6d 20 3d 20 2d 69 6d | d->im = -im|
00001030 3b 0a 20 20 20 20 7d 0a 20 20 20 20 6a 20 3d 20 |;. }. j = |
00001040 69 3b 0a 20 20 20 20 6c 20 3d 20 6e 2d 6a 3b 0a |i;. l = n-j;.|
00001050 20 20 20 20 64 31 20 3d 20 64 61 74 61 2b 6a 3b | d1 = data+j;|
00001060 0a 20 20 20 20 64 32 20 3d 20 64 61 74 61 2b 6c |. d2 = data+l|
00001070 3b 0a 20 20 20 20 66 6f 72 20 28 69 63 3d 69 2d |;. for (ic=i-|
00001080 31 3b 20 69 63 3e 3d 30 3b 20 69 63 2d 2d 29 20 |1; ic>=0; ic--) |
00001090 7b 0a 20 20 20 20 20 20 6b 20 3d 20 69 63 3d 3d |{. k = ic==|
000010a0 30 20 3f 20 30 20 3a 20 6e 2d 69 63 3b 0a 20 20 |0 ? 0 : n-ic;. |
000010b0 20 20 20 20 70 68 61 73 65 20 3d 20 72 61 6e 64 | phase = rand|
000010c0 28 29 3e 3e 31 33 3b 0a 20 20 20 20 20 20 72 61 |()>>13;. ra|
000010d0 64 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 28 70 |d = mul_frac16(p|
000010e0 6f 77 31 36 28 73 71 72 74 5f 66 72 61 63 31 36 |ow16(sqrt_frac16|
000010f0 28 28 69 63 2a 69 63 2b 6a 32 29 2a 28 6f 6e 65 |((ic*ic+j2)*(one|
00001100 2f 34 29 29 2a 32 2c 20 6d 62 70 6f 62 74 29 2c |/4))*2, mbpobt),|
00001110 20 67 61 75 73 73 31 36 28 29 29 3b 0a 20 20 20 | gauss16());. |
00001120 20 20 20 72 65 20 3d 20 6d 75 6c 5f 66 72 61 63 | re = mul_frac|
00001130 31 36 63 28 72 61 64 2c 20 63 6f 73 31 36 28 70 |16c(rad, cos16(p|
00001140 68 61 73 65 29 29 3b 0a 20 20 20 20 20 20 69 6d |hase));. im|
00001150 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 28 72 | = mul_frac16c(r|
00001160 61 64 2c 20 73 69 6e 31 36 28 70 68 61 73 65 29 |ad, sin16(phase)|
00001170 29 3b 0a 20 20 20 20 20 20 64 20 3d 20 64 31 2b |);. d = d1+|
00001180 69 63 2a 6e 3b 0a 20 20 20 20 20 20 64 2d 3e 72 |ic*n;. d->r|
00001190 65 20 3d 20 72 65 3b 0a 20 20 20 20 20 20 64 2d |e = re;. d-|
000011a0 3e 69 6d 20 3d 20 69 6d 3b 0a 20 20 20 20 20 20 |>im = im;. |
000011b0 64 20 3d 20 64 32 2b 6b 2a 6e 3b 0a 20 20 20 20 |d = d2+k*n;. |
000011c0 20 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 20 | d->re = re;. |
000011d0 20 20 20 20 64 2d 3e 69 6d 20 3d 20 2d 69 6d 3b | d->im = -im;|
000011e0 0a 20 20 20 20 7d 0a 20 20 20 20 69 66 20 28 69 |. }. if (i|
000011f0 3c 6e 2f 32 29 20 7b 0a 20 20 20 20 20 20 69 63 |<n/2) {. ic|
00001200 20 3d 20 6e 2d 69 3b 0a 20 20 20 20 20 20 6b 20 | = n-i;. k |
00001210 3d 20 69 3b 0a 20 20 20 20 20 20 64 31 20 3d 20 |= i;. d1 = |
00001220 64 61 74 61 2b 69 63 2a 6e 3b 0a 20 20 20 20 20 |data+ic*n;. |
00001230 20 64 32 20 3d 20 64 61 74 61 2b 6b 2a 6e 3b 0a | d2 = data+k*n;.|
00001240 20 20 20 20 20 20 66 6f 72 20 28 6a 3d 31 3b 20 | for (j=1; |
00001250 6a 3c 3d 69 3b 20 6a 2b 2b 29 20 7b 0a 20 20 20 |j<=i; j++) {. |
00001260 20 20 20 20 20 6c 20 3d 20 6e 2d 6a 3b 0a 20 20 | l = n-j;. |
00001270 20 20 20 20 20 20 70 68 61 73 65 20 3d 20 72 61 | phase = ra|
00001280 6e 64 28 29 3e 3e 31 33 3b 0a 20 20 20 20 20 20 |nd()>>13;. |
00001290 20 20 72 61 64 20 3d 20 6d 75 6c 5f 66 72 61 63 | rad = mul_frac|
000012a0 31 36 28 70 6f 77 31 36 28 73 71 72 74 5f 66 72 |16(pow16(sqrt_fr|
000012b0 61 63 31 36 28 28 69 32 2b 6a 2a 6a 29 2a 28 6f |ac16((i2+j*j)*(o|
000012c0 6e 65 2f 34 29 29 2a 32 2c 20 6d 62 70 6f 62 74 |ne/4))*2, mbpobt|
000012d0 29 2c 20 67 61 75 73 73 31 36 28 29 29 3b 0a 20 |), gauss16());. |
000012e0 20 20 20 20 20 20 20 72 65 20 3d 20 6d 75 6c 5f | re = mul_|
000012f0 66 72 61 63 31 36 63 28 72 61 64 2c 20 63 6f 73 |frac16c(rad, cos|
00001300 31 36 28 70 68 61 73 65 29 29 3b 0a 20 20 20 20 |16(phase));. |
00001310 20 20 20 20 69 6d 20 3d 20 6d 75 6c 5f 66 72 61 | im = mul_fra|
00001320 63 31 36 63 28 72 61 64 2c 20 73 69 6e 31 36 28 |c16c(rad, sin16(|
00001330 70 68 61 73 65 29 29 3b 0a 20 20 20 20 20 20 20 |phase));. |
00001340 20 64 20 3d 20 64 31 2b 6a 3b 0a 20 20 20 20 20 | d = d1+j;. |
00001350 20 20 20 64 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 | d->re = re;. |
00001360 20 20 20 20 20 20 20 64 2d 3e 69 6d 20 3d 20 69 | d->im = i|
00001370 6d 3b 0a 20 20 20 20 20 20 20 20 64 20 3d 20 64 |m;. d = d|
00001380 32 2b 6c 3b 0a 20 20 20 20 20 20 20 20 64 2d 3e |2+l;. d->|
00001390 72 65 20 3d 20 72 65 3b 0a 20 20 20 20 20 20 20 |re = re;. |
000013a0 20 64 2d 3e 69 6d 20 3d 20 2d 69 6d 3b 0a 20 20 | d->im = -im;. |
000013b0 20 20 20 20 7d 0a 20 20 20 20 20 20 6a 3d 69 3b | }. j=i;|
000013c0 0a 20 20 20 20 20 20 6c 20 3d 20 6e 2d 6a 3b 0a |. l = n-j;.|
000013d0 20 20 20 20 20 20 64 31 20 3d 20 64 61 74 61 2b | d1 = data+|
000013e0 6a 3b 0a 20 20 20 20 20 20 64 32 20 3d 20 64 61 |j;. d2 = da|
000013f0 74 61 2b 6c 3b 0a 20 20 20 20 20 20 66 6f 72 20 |ta+l;. for |
00001400 28 69 63 2b 2b 3b 20 69 63 3c 6e 3b 20 69 63 2b |(ic++; ic<n; ic+|
00001410 2b 29 20 7b 0a 20 20 20 20 20 20 20 20 6b 20 3d |+) {. k =|
00001420 20 6e 2d 69 63 3b 0a 20 20 20 20 20 20 20 20 70 | n-ic;. p|
00001430 68 61 73 65 20 3d 20 72 61 6e 64 28 29 3e 3e 31 |hase = rand()>>1|
00001440 33 3b 0a 20 20 20 20 20 20 20 20 72 61 64 20 3d |3;. rad =|
00001450 20 6d 75 6c 5f 66 72 61 63 31 36 28 70 6f 77 31 | mul_frac16(pow1|
00001460 36 28 73 71 72 74 5f 66 72 61 63 31 36 28 28 6b |6(sqrt_frac16((k|
00001470 2a 6b 2b 6a 32 29 2a 28 6f 6e 65 2f 34 29 29 2a |*k+j2)*(one/4))*|
00001480 32 2c 20 6d 62 70 6f 62 74 29 2c 20 67 61 75 73 |2, mbpobt), gaus|
00001490 73 31 36 28 29 29 3b 0a 20 20 20 20 20 20 20 20 |s16());. |
000014a0 72 65 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 |re = mul_frac16c|
000014b0 28 72 61 64 2c 20 63 6f 73 31 36 28 70 68 61 73 |(rad, cos16(phas|
000014c0 65 29 29 3b 0a 20 20 20 20 20 20 20 20 69 6d 20 |e));. im |
000014d0 3d 20 6d 75 6c 5f 66 72 61 63 31 36 63 28 72 61 |= mul_frac16c(ra|
000014e0 64 2c 20 73 69 6e 31 36 28 70 68 61 73 65 29 29 |d, sin16(phase))|
000014f0 3b 0a 20 20 20 20 20 20 20 20 64 20 3d 20 64 31 |;. d = d1|
00001500 2b 69 63 2a 6e 3b 0a 20 20 20 20 20 20 20 20 64 |+ic*n;. d|
00001510 2d 3e 72 65 20 3d 20 72 65 3b 0a 20 20 20 20 20 |->re = re;. |
00001520 20 20 20 64 2d 3e 69 6d 20 3d 20 69 6d 3b 0a 20 | d->im = im;. |
00001530 20 20 20 20 20 20 20 64 20 3d 20 64 32 2b 6b 2a | d = d2+k*|
00001540 6e 3b 0a 20 20 20 20 20 20 20 20 64 2d 3e 72 65 |n;. d->re|
00001550 20 3d 20 72 65 3b 0a 20 20 20 20 20 20 20 20 64 | = re;. d|
00001560 2d 3e 69 6d 20 3d 20 2d 69 6d 3b 0a 20 20 20 20 |->im = -im;. |
00001570 20 20 7d 0a 20 20 20 20 7d 0a 20 20 7d 0a 20 20 | }. }. }. |
00001580 28 64 61 74 61 2b 6e 2a 6e 2f 32 29 2d 3e 69 6d |(data+n*n/2)->im|
00001590 20 3d 20 28 64 61 74 61 2b 6e 2f 32 29 2d 3e 69 | = (data+n/2)->i|
000015a0 6d 20 3d 20 28 64 61 74 61 2b 6e 2a 28 6e 2b 31 |m = (data+n*(n+1|
000015b0 29 2f 32 29 2d 3e 69 6d 20 3d 20 30 3b 0a 0a 20 |)/2)->im = 0;.. |
000015c0 20 72 65 74 75 72 6e 20 54 52 55 45 3b 0a 7d 0a | return TRUE;.}.|
000015d0 0a 42 4f 4f 4c 20 6e 6f 72 6d 61 6c 69 73 65 5f |.BOOL normalise_|
000015e0 61 6e 64 5f 63 6f 6d 70 61 63 74 28 76 6f 69 64 |and_compact(void|
000015f0 20 2a 2a 64 61 74 61 2c 20 69 6e 74 20 6c 69 6d | **data, int lim|
00001600 69 74 29 0a 7b 0a 20 20 69 6e 74 20 6e 3d 31 3c |it).{. int n=1<|
00001610 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 20 6d 69 |<limit;. int mi|
00001620 6e 3d 31 3c 3c 33 30 3b 0a 20 20 69 6e 74 20 6d |n=1<<30;. int m|
00001630 61 78 3d 2d 6d 69 6e 3b 0a 20 20 69 6e 74 20 73 |ax=-min;. int s|
00001640 63 61 6c 65 2c 20 76 3b 0a 20 20 63 6f 6d 70 6c |cale, v;. compl|
00001650 65 78 20 2a 64 31 2c 20 2a 64 31 65 3b 0a 20 20 |ex *d1, *d1e;. |
00001660 69 6e 74 20 2a 64 32 3b 0a 0a 20 20 66 6f 72 20 |int *d2;.. for |
00001670 28 64 31 3d 2a 64 61 74 61 2c 20 64 31 65 3d 64 |(d1=*data, d1e=d|
00001680 31 2b 6e 2a 6e 3b 20 64 31 3c 64 31 65 3b 20 64 |1+n*n; d1<d1e; d|
00001690 31 2b 2b 29 20 7b 0a 20 20 20 20 69 66 20 28 64 |1++) {. if (d|
000016a0 31 2d 3e 72 65 3c 6d 69 6e 29 20 6d 69 6e 3d 64 |1->re<min) min=d|
000016b0 31 2d 3e 72 65 3b 0a 20 20 20 20 69 66 20 28 64 |1->re;. if (d|
000016c0 31 2d 3e 72 65 3e 6d 61 78 29 20 6d 61 78 3d 64 |1->re>max) max=d|
000016d0 31 2d 3e 72 65 3b 0a 20 20 7d 0a 20 20 73 63 61 |1->re;. }. sca|
000016e0 6c 65 20 3d 20 64 69 76 5f 66 72 61 63 31 36 28 |le = div_frac16(|
000016f0 6f 6e 65 2c 20 6d 61 78 2d 6d 69 6e 29 3b 0a 20 |one, max-min);. |
00001700 20 66 6f 72 20 28 64 32 3d 2a 64 61 74 61 2c 20 | for (d2=*data, |
00001710 64 31 3d 2a 64 61 74 61 2c 20 64 31 65 3d 64 31 |d1=*data, d1e=d1|
00001720 2b 6e 2a 6e 3b 20 64 31 3c 64 31 65 3b 20 64 31 |+n*n; d1<d1e; d1|
00001730 2b 2b 2c 20 64 32 2b 2b 29 20 7b 0a 20 20 20 20 |++, d2++) {. |
00001740 76 20 3d 20 6d 75 6c 5f 66 72 61 63 31 36 28 64 |v = mul_frac16(d|
00001750 31 2d 3e 72 65 2d 6d 69 6e 2c 20 73 63 61 6c 65 |1->re-min, scale|
00001760 29 3b 0a 20 20 20 20 69 66 20 28 76 3e 3d 6f 6e |);. if (v>=on|
00001770 65 29 20 76 3d 6f 6e 65 2d 31 3b 0a 20 20 20 20 |e) v=one-1;. |
00001780 69 66 20 28 76 3c 30 29 20 20 20 20 76 3d 30 3b |if (v<0) v=0;|
00001790 0a 20 20 20 20 76 20 3d 20 70 6f 77 31 36 28 76 |. v = pow16(v|
000017a0 2c 20 70 6f 77 73 63 61 6c 65 29 3b 0a 20 20 20 |, powscale);. |
000017b0 20 69 66 20 28 76 3c 73 65 61 6c 65 76 65 6c 29 | if (v<sealevel)|
000017c0 20 76 3d 73 65 61 6c 65 76 65 6c 3b 0a 20 20 20 | v=sealevel;. |
000017d0 20 69 66 20 28 6c 69 6e 73 63 61 6c 65 21 3d 6f | if (linscale!=o|
000017e0 6e 65 29 20 76 3d 6d 75 6c 5f 66 72 61 63 31 36 |ne) v=mul_frac16|
000017f0 63 28 76 2c 20 6c 69 6e 73 63 61 6c 65 29 3b 0a |c(v, linscale);.|
00001800 20 20 20 20 69 66 20 28 76 3e 3d 6f 6e 65 29 20 | if (v>=one) |
00001810 76 3d 6f 6e 65 2d 31 3b 0a 20 20 20 20 69 66 20 |v=one-1;. if |
00001820 28 76 3c 30 29 20 20 20 20 76 3d 30 3b 0a 20 20 |(v<0) v=0;. |
00001830 20 20 2a 64 32 20 3d 20 76 3b 0a 20 20 7d 0a 20 | *d2 = v;. }. |
00001840 20 2f 2a 4e 42 20 67 69 76 65 6e 20 72 61 6e 67 | /*NB given rang|
00001850 65 20 6f 6e 20 64 61 74 61 2c 20 63 6f 75 6c 64 |e on data, could|
00001860 20 70 61 63 6b 20 65 61 63 68 20 76 61 6c 75 65 | pack each value|
00001870 20 69 6e 74 6f 20 32 20 63 68 61 72 73 2c 20 62 | into 2 chars, b|
00001880 75 74 20 69 73 20 73 6c 6f 77 65 72 20 74 6f 20 |ut is slower to |
00001890 72 65 61 64 20 26 20 77 72 69 74 65 0a 20 20 20 |read & write. |
000018a0 20 68 61 6c 66 20 77 6f 72 64 73 20 74 68 61 6e | half words than|
000018b0 20 77 68 6f 6c 65 20 77 6f 72 64 73 2c 20 73 6f | whole words, so|
000018c0 20 6a 75 73 74 20 72 65 64 75 63 65 20 66 72 6f | just reduce fro|
000018d0 6d 20 32 20 69 6e 74 73 20 64 6f 77 6e 20 74 6f |m 2 ints down to|
000018e0 20 31 20 69 6e 74 2e 2a 2f 0a 20 20 72 65 74 75 | 1 int.*/. retu|
000018f0 72 6e 20 66 6c 65 78 5f 65 78 74 65 6e 64 28 28 |rn flex_extend((|
00001900 66 6c 65 78 5f 70 74 72 29 64 61 74 61 2c 20 6e |flex_ptr)data, n|
00001910 2a 6e 2a 73 69 7a 65 6f 66 28 69 6e 74 29 29 3b |*n*sizeof(int));|
00001920 0a 7d 0a 0a 42 4f 4f 4c 20 63 6c 69 70 5f 66 74 |.}..BOOL clip_ft|
00001930 28 63 6f 6d 70 6c 65 78 20 2a 64 61 74 61 2c 20 |(complex *data, |
00001940 69 6e 74 20 6c 69 6d 69 74 2c 20 69 6e 74 20 72 |int limit, int r|
00001950 6c 29 0a 7b 0a 20 20 69 6e 74 20 6e 3d 31 3c 3c |l).{. int n=1<<|
00001960 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 20 69 2c 20 |limit;. int i, |
00001970 6a 3b 0a 20 20 69 6e 74 20 72 68 3d 6e 2d 72 6c |j;. int rh=n-rl|
00001980 3b 0a 0a 20 20 66 6f 72 20 28 69 3d 30 3b 20 69 |;.. for (i=0; i|
00001990 3c 6e 3b 20 69 2b 2b 29 0a 20 20 20 20 66 6f 72 |<n; i++). for|
000019a0 20 28 6a 3d 30 3b 20 6a 3c 6e 3b 20 6a 2b 2b 2c | (j=0; j<n; j++,|
000019b0 20 64 61 74 61 2b 2b 29 0a 20 20 20 20 20 20 69 | data++). i|
000019c0 66 20 28 69 3e 72 6c 20 26 26 20 69 3c 72 68 20 |f (i>rl && i<rh |
000019d0 7c 7c 20 6a 3e 72 6c 20 26 26 20 6a 3c 72 68 29 ||| j>rl && j<rh)|
000019e0 20 64 61 74 61 2d 3e 72 65 3d 64 61 74 61 2d 3e | data->re=data->|
000019f0 69 6d 3d 30 3b 0a 0a 20 20 72 65 74 75 72 6e 20 |im=0;.. return |
00001a00 54 52 55 45 3b 0a 7d 0a 0a 42 4f 4f 4c 20 73 63 |TRUE;.}..BOOL sc|
00001a10 72 6f 6c 6c 28 69 6e 74 20 2a 64 61 74 61 2c 20 |roll(int *data, |
00001a20 69 6e 74 20 6c 69 6d 69 74 2c 20 73 63 72 6f 6c |int limit, scrol|
00001a30 6c 64 69 72 20 6b 65 79 2c 20 69 6e 74 20 6f 66 |ldir key, int of|
00001a40 66 73 65 74 29 0a 7b 0a 20 20 69 6e 74 20 6e 3d |fset).{. int n=|
00001a50 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 20 |1<<limit;. int |
00001a60 69 2c 20 6a 3b 0a 20 20 69 6e 74 20 2a 70 2c 20 |i, j;. int *p, |
00001a70 2a 70 31 2c 20 2a 70 32 3b 0a 20 20 69 6e 74 20 |*p1, *p2;. int |
00001a80 62 75 66 5b 35 31 32 5d 3b 0a 20 20 69 6e 74 20 |buf[512];. int |
00001a90 73 69 3d 6f 66 66 73 65 74 3c 3c 28 39 2d 6c 69 |si=offset<<(9-li|
00001aa0 6d 69 74 29 3b 0a 0a 20 20 73 77 69 74 63 68 20 |mit);.. switch |
00001ab0 28 6b 65 79 29 20 7b 0a 20 20 20 20 63 61 73 65 |(key) {. case|
00001ac0 20 6c 65 66 74 3a 0a 20 20 20 20 20 20 73 63 72 | left:. scr|
00001ad0 6f 6c 6c 78 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d |ollx += orient==|
00001ae0 30 20 3f 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e |0 ? -si : (orien|
00001af0 74 3d 3d 32 20 3f 20 73 69 20 3a 20 30 29 3b 0a |t==2 ? si : 0);.|
00001b00 20 20 20 20 20 20 73 63 72 6f 6c 6c 79 20 2b 3d | scrolly +=|
00001b10 20 6f 72 69 65 6e 74 3d 3d 31 20 3f 20 2d 73 69 | orient==1 ? -si|
00001b20 20 3a 20 28 6f 72 69 65 6e 74 3d 3d 33 20 3f 20 | : (orient==3 ? |
00001b30 73 69 20 3a 20 30 29 3b 0a 20 20 20 20 20 20 62 |si : 0);. b|
00001b40 72 65 61 6b 3b 0a 20 20 20 20 63 61 73 65 20 72 |reak;. case r|
00001b50 69 67 68 74 3a 0a 20 20 20 20 20 20 73 63 72 6f |ight:. scro|
00001b60 6c 6c 78 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d 32 |llx += orient==2|
00001b70 20 3f 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e 74 | ? -si : (orient|
00001b80 3d 3d 30 20 3f 20 73 69 20 3a 20 30 29 3b 0a 20 |==0 ? si : 0);. |
00001b90 20 20 20 20 20 73 63 72 6f 6c 6c 79 20 2b 3d 20 | scrolly += |
00001ba0 6f 72 69 65 6e 74 3d 3d 33 20 3f 20 2d 73 69 20 |orient==3 ? -si |
00001bb0 3a 20 28 6f 72 69 65 6e 74 3d 3d 31 20 3f 20 73 |: (orient==1 ? s|
00001bc0 69 20 3a 20 30 29 3b 0a 20 20 20 20 20 20 6f 66 |i : 0);. of|
00001bd0 66 73 65 74 3d 6e 2d 6f 66 66 73 65 74 3b 0a 20 |fset=n-offset;. |
00001be0 20 20 20 20 20 6b 65 79 3d 6c 65 66 74 3b 0a 20 | key=left;. |
00001bf0 20 20 20 20 20 62 72 65 61 6b 3b 0a 20 20 20 20 | break;. |
00001c00 63 61 73 65 20 64 6f 77 6e 3a 0a 20 20 20 20 20 |case down:. |
00001c10 20 73 63 72 6f 6c 6c 78 20 2b 3d 20 6f 72 69 65 | scrollx += orie|
00001c20 6e 74 3d 3d 33 20 3f 20 2d 73 69 20 3a 20 28 6f |nt==3 ? -si : (o|
00001c30 72 69 65 6e 74 3d 3d 31 20 3f 20 73 69 20 3a 20 |rient==1 ? si : |
00001c40 30 29 3b 0a 20 20 20 20 20 20 73 63 72 6f 6c 6c |0);. scroll|
00001c50 79 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d 30 20 3f |y += orient==0 ?|
00001c60 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e 74 3d 3d | -si : (orient==|
00001c70 32 20 3f 20 73 69 20 3a 20 30 29 3b 0a 20 20 20 |2 ? si : 0);. |
00001c80 20 20 20 62 72 65 61 6b 3b 0a 20 20 20 20 63 61 | break;. ca|
00001c90 73 65 20 75 70 3a 0a 20 20 20 20 20 20 73 63 72 |se up:. scr|
00001ca0 6f 6c 6c 78 20 2b 3d 20 6f 72 69 65 6e 74 3d 3d |ollx += orient==|
00001cb0 31 20 3f 20 2d 73 69 20 3a 20 28 6f 72 69 65 6e |1 ? -si : (orien|
00001cc0 74 3d 3d 33 20 3f 20 73 69 20 3a 20 30 29 3b 0a |t==3 ? si : 0);.|
00001cd0 20 20 20 20 20 20 73 63 72 6f 6c 6c 79 20 2b 3d | scrolly +=|
00001ce0 20 6f 72 69 65 6e 74 3d 3d 32 20 3f 20 2d 73 69 | orient==2 ? -si|
00001cf0 20 3a 20 28 6f 72 69 65 6e 74 3d 3d 30 20 3f 20 | : (orient==0 ? |
00001d00 73 69 20 3a 20 30 29 3b 0a 20 20 20 20 20 20 6f |si : 0);. o|
00001d10 66 66 73 65 74 3d 6e 2d 6f 66 66 73 65 74 3b 0a |ffset=n-offset;.|
00001d20 20 20 20 20 20 20 6b 65 79 3d 64 6f 77 6e 3b 0a | key=down;.|
00001d30 20 20 20 20 20 20 62 72 65 61 6b 3b 0a 20 20 7d | break;. }|
00001d40 0a 20 20 73 63 72 6f 6c 6c 78 20 26 3d 20 35 31 |. scrollx &= 51|
00001d50 31 3b 0a 20 20 73 63 72 6f 6c 6c 79 20 26 3d 20 |1;. scrolly &= |
00001d60 35 31 31 3b 0a 20 20 6f 66 66 73 65 74 20 20 26 |511;. offset &|
00001d70 3d 20 6e 2d 31 3b 0a 0a 20 20 69 66 20 28 6b 65 |= n-1;.. if (ke|
00001d80 79 3d 3d 6c 65 66 74 29 20 7b 0a 20 20 20 20 66 |y==left) {. f|
00001d90 6f 72 20 28 69 3d 30 2c 20 70 3d 64 61 74 61 3b |or (i=0, p=data;|
00001da0 20 69 3c 6e 3b 20 69 2b 2b 2c 20 70 2b 3d 6e 29 | i<n; i++, p+=n)|
00001db0 20 7b 0a 20 20 20 20 20 20 66 6f 72 20 28 6a 3d | {. for (j=|
00001dc0 6e 2d 6f 66 66 73 65 74 2c 20 70 31 3d 70 2b 6a |n-offset, p1=p+j|
00001dd0 2c 20 70 32 3d 62 75 66 3b 20 6a 3c 6e 3b 20 6a |, p2=buf; j<n; j|
00001de0 2b 2b 2c 20 70 31 2b 2b 2c 20 70 32 2b 2b 29 20 |++, p1++, p2++) |
00001df0 2a 70 32 3d 2a 70 31 3b 0a 20 20 20 20 20 20 66 |*p2=*p1;. f|
00001e00 6f 72 20 28 6a 3d 6e 2d 6f 66 66 73 65 74 2d 31 |or (j=n-offset-1|
00001e10 2c 20 70 31 3d 70 2b 6a 2c 20 70 32 3d 70 2b 6a |, p1=p+j, p2=p+j|
00001e20 2b 6f 66 66 73 65 74 3b 20 6a 3e 3d 30 3b 20 6a |+offset; j>=0; j|
00001e30 2d 2d 2c 20 70 31 2d 2d 2c 20 70 32 2d 2d 29 20 |--, p1--, p2--) |
00001e40 2a 70 32 3d 2a 70 31 3b 0a 20 20 20 20 20 20 66 |*p2=*p1;. f|
00001e50 6f 72 20 28 6a 3d 30 2c 20 70 31 3d 62 75 66 2c |or (j=0, p1=buf,|
00001e60 20 70 32 3d 70 3b 20 6a 3c 6f 66 66 73 65 74 3b | p2=p; j<offset;|
00001e70 20 6a 2b 2b 2c 20 70 31 2b 2b 2c 20 70 32 2b 2b | j++, p1++, p2++|
00001e80 29 20 2a 70 32 3d 2a 70 31 3b 0a 20 20 20 20 7d |) *p2=*p1;. }|
00001e90 0a 20 20 7d 0a 20 20 65 6c 73 65 20 7b 0a 20 20 |. }. else {. |
00001ea0 20 20 66 6f 72 20 28 6a 3d 30 2c 20 70 3d 64 61 | for (j=0, p=da|
00001eb0 74 61 3b 20 6a 3c 6e 3b 20 6a 2b 2b 2c 20 70 2b |ta; j<n; j++, p+|
00001ec0 2b 29 20 7b 0a 20 20 20 20 20 20 66 6f 72 20 28 |+) {. for (|
00001ed0 69 3d 6e 2d 6f 66 66 73 65 74 2c 20 70 31 3d 70 |i=n-offset, p1=p|
00001ee0 2b 69 2a 6e 2c 20 70 32 3d 62 75 66 3b 20 69 3c |+i*n, p2=buf; i<|
00001ef0 6e 3b 20 69 2b 2b 2c 20 70 31 2b 3d 6e 2c 20 70 |n; i++, p1+=n, p|
00001f00 32 2b 2b 29 20 2a 70 32 3d 2a 70 31 3b 0a 20 20 |2++) *p2=*p1;. |
00001f10 20 20 20 20 66 6f 72 20 28 69 3d 6e 2d 6f 66 66 | for (i=n-off|
00001f20 73 65 74 2d 31 2c 20 70 31 3d 70 2b 69 2a 6e 2c |set-1, p1=p+i*n,|
00001f30 20 70 32 3d 70 2b 28 69 2b 6f 66 66 73 65 74 29 | p2=p+(i+offset)|
00001f40 2a 6e 3b 20 69 3e 3d 30 3b 20 69 2d 2d 2c 20 70 |*n; i>=0; i--, p|
00001f50 31 2d 3d 6e 2c 20 70 32 2d 3d 6e 29 20 2a 70 32 |1-=n, p2-=n) *p2|
00001f60 3d 2a 70 31 3b 0a 20 20 20 20 20 20 66 6f 72 20 |=*p1;. for |
00001f70 28 69 3d 30 2c 20 70 31 3d 62 75 66 2c 20 70 32 |(i=0, p1=buf, p2|
00001f80 3d 70 3b 20 69 3c 6f 66 66 73 65 74 3b 20 69 2b |=p; i<offset; i+|
00001f90 2b 2c 20 70 31 2b 2b 2c 20 70 32 2b 3d 6e 29 20 |+, p1++, p2+=n) |
00001fa0 2a 70 32 3d 2a 70 31 3b 0a 20 20 20 20 7d 0a 20 |*p2=*p1;. }. |
00001fb0 20 7d 0a 0a 20 20 72 65 74 75 72 6e 20 54 52 55 | }.. return TRU|
00001fc0 45 3b 0a 7d 0a 0a 42 4f 4f 4c 20 72 6f 74 61 74 |E;.}..BOOL rotat|
00001fd0 65 28 69 6e 74 20 2a 64 61 74 61 2c 20 69 6e 74 |e(int *data, int|
00001fe0 20 6c 69 6d 69 74 2c 20 42 4f 4f 4c 20 63 77 29 | limit, BOOL cw)|
00001ff0 0a 7b 0a 20 20 69 6e 74 20 6e 3d 31 3c 3c 6c 69 |.{. int n=1<<li|
00002000 6d 69 74 3b 0a 20 20 69 6e 74 20 69 2c 20 6a 3b |mit;. int i, j;|
00002010 0a 20 20 69 6e 74 20 6b 20 3d 20 6e 2d 31 3b 0a |. int k = n-1;.|
00002020 20 20 69 6e 74 20 6c 20 3d 20 6e 2f 32 3b 0a 20 | int l = n/2;. |
00002030 20 69 6e 74 20 2a 71 31 2c 20 2a 71 32 2c 20 2a | int *q1, *q2, *|
00002040 71 33 2c 20 2a 71 34 3b 0a 20 20 69 6e 74 20 74 |q3, *q4;. int t|
00002050 3b 0a 0a 20 20 69 66 20 28 63 77 29 20 7b 0a 20 |;.. if (cw) {. |
00002060 20 20 20 6f 72 69 65 6e 74 2b 2b 3b 0a 20 20 20 | orient++;. |
00002070 20 66 6f 72 20 28 69 3d 30 2c 20 71 31 3d 64 61 | for (i=0, q1=da|
00002080 74 61 2c 20 71 32 3d 64 61 74 61 2b 6e 2a 6b 2c |ta, q2=data+n*k,|
00002090 20 71 33 3d 64 61 74 61 2b 28 6e 2b 31 29 2a 6b | q3=data+(n+1)*k|
000020a0 2c 20 71 34 3d 64 61 74 61 2b 6b 3b 20 69 3c 6c |, q4=data+k; i<l|
000020b0 3b 0a 20 20 09 09 09 20 20 20 69 2b 2b 2c 20 71 |;. ... i++, q|
000020c0 31 2b 3d 6c 2c 20 71 32 2b 3d 6c 2a 6e 2b 31 2c |1+=l, q2+=l*n+1,|
000020d0 20 71 33 2d 3d 6c 2c 20 71 34 2d 3d 6c 2a 6e 2b | q3-=l, q4-=l*n+|
000020e0 31 29 0a 20 20 20 20 20 20 66 6f 72 20 28 6a 3d |1). for (j=|
000020f0 30 3b 20 6a 3c 6c 3b 20 6a 2b 2b 2c 20 71 31 2b |0; j<l; j++, q1+|
00002100 2b 2c 20 71 32 2d 3d 6e 2c 20 71 33 2d 2d 2c 20 |+, q2-=n, q3--, |
00002110 71 34 2b 3d 6e 29 20 7b 0a 20 20 20 20 20 20 20 |q4+=n) {. |
00002120 20 74 3d 2a 71 31 3b 0a 20 20 20 20 20 20 20 20 | t=*q1;. |
00002130 2a 71 31 3d 2a 71 34 3b 0a 20 20 20 20 20 20 20 |*q1=*q4;. |
00002140 20 2a 71 34 3d 2a 71 33 3b 0a 20 20 20 20 20 20 | *q4=*q3;. |
00002150 20 20 2a 71 33 3d 2a 71 32 3b 0a 20 20 20 20 20 | *q3=*q2;. |
00002160 20 20 20 2a 71 32 3d 74 3b 0a 20 20 20 20 20 20 | *q2=t;. |
00002170 7d 0a 20 20 7d 0a 20 20 65 6c 73 65 20 7b 0a 20 |}. }. else {. |
00002180 20 20 20 6f 72 69 65 6e 74 2d 2d 3b 0a 20 20 20 | orient--;. |
00002190 20 66 6f 72 20 28 69 3d 30 2c 20 71 31 3d 64 61 | for (i=0, q1=da|
000021a0 74 61 2c 20 71 32 3d 64 61 74 61 2b 6e 2a 6b 2c |ta, q2=data+n*k,|
000021b0 20 71 33 3d 64 61 74 61 2b 28 6e 2b 31 29 2a 6b | q3=data+(n+1)*k|
000021c0 2c 20 71 34 3d 64 61 74 61 2b 6b 3b 20 69 3c 6c |, q4=data+k; i<l|
000021d0 3b 0a 20 20 09 09 09 20 20 20 69 2b 2b 2c 20 71 |;. ... i++, q|
000021e0 31 2b 3d 6c 2c 20 71 32 2b 3d 6c 2a 6e 2b 31 2c |1+=l, q2+=l*n+1,|
000021f0 20 71 33 2d 3d 6c 2c 20 71 34 2d 3d 6c 2a 6e 2b | q3-=l, q4-=l*n+|
00002200 31 29 0a 20 20 20 20 20 20 66 6f 72 20 28 6a 3d |1). for (j=|
00002210 30 3b 20 6a 3c 6c 3b 20 6a 2b 2b 2c 20 71 31 2b |0; j<l; j++, q1+|
00002220 2b 2c 20 71 32 2d 3d 6e 2c 20 71 33 2d 2d 2c 20 |+, q2-=n, q3--, |
00002230 71 34 2b 3d 6e 29 20 7b 0a 20 20 20 20 20 20 20 |q4+=n) {. |
00002240 20 74 3d 2a 71 31 3b 0a 20 20 20 20 20 20 20 20 | t=*q1;. |
00002250 2a 71 31 3d 2a 71 32 3b 0a 20 20 20 20 20 20 20 |*q1=*q2;. |
00002260 20 2a 71 32 3d 2a 71 33 3b 0a 20 20 20 20 20 20 | *q2=*q3;. |
00002270 20 20 2a 71 33 3d 2a 71 34 3b 0a 20 20 20 20 20 | *q3=*q4;. |
00002280 20 20 20 2a 71 34 3d 74 3b 0a 20 20 20 20 20 20 | *q4=t;. |
00002290 7d 0a 20 20 7d 0a 20 20 6f 72 69 65 6e 74 20 26 |}. }. orient &|
000022a0 3d 20 33 3b 0a 0a 20 20 72 65 74 75 72 6e 20 54 |= 3;.. return T|
000022b0 52 55 45 3b 0a 7d 0a 0a 42 4f 4f 4c 20 67 72 61 |RUE;.}..BOOL gra|
000022c0 70 68 28 69 6e 74 20 2a 64 61 74 61 2c 20 69 6e |ph(int *data, in|
000022d0 74 20 6c 69 6d 69 74 29 0a 7b 0a 20 20 69 6e 74 |t limit).{. int|
000022e0 20 6e 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 | n=1<<limit;. i|
000022f0 6e 74 20 69 2c 20 6a 2c 20 6b 3b 0a 20 20 69 6e |nt i, j, k;. in|
00002300 74 20 70 6c 2c 20 70 72 3b 0a 20 20 69 6e 74 20 |t pl, pr;. int |
00002310 76 6c 2c 20 76 72 2c 20 76 2c 20 76 32 3b 0a 20 |vl, vr, v, v2;. |
00002320 20 69 6e 74 20 68 6f 72 69 7a 6f 6e 5b 31 30 32 | int horizon[102|
00002330 36 5d 3b 0a 20 20 69 6e 74 20 68 6f 72 69 7a 6f |6];. int horizo|
00002340 6e 32 5b 31 30 32 36 5d 3b 0a 20 20 69 6e 74 20 |n2[1026];. int |
00002350 63 6f 73 20 3d 20 63 6f 73 31 36 28 74 69 6c 74 |cos = cos16(tilt|
00002360 29 3b 0a 20 20 69 6e 74 20 73 69 6e 20 3d 20 73 |);. int sin = s|
00002370 69 6e 31 36 28 74 69 6c 74 29 3b 0a 20 20 69 6e |in16(tilt);. in|
00002380 74 20 69 6d 61 67 65 5f 79 6c 3d 31 30 32 34 3b |t image_yl=1024;|
00002390 0a 20 20 69 6e 74 20 69 6d 61 67 65 5f 79 68 3d |. int image_yh=|
000023a0 30 3b 0a 0a 20 20 66 6f 72 20 28 69 3d 30 3b 20 |0;.. for (i=0; |
000023b0 69 3c 31 30 32 36 3b 20 69 2b 2b 29 20 7b 0a 20 |i<1026; i++) {. |
000023c0 20 20 20 68 6f 72 69 7a 6f 6e 5b 69 5d 3d 30 3b | horizon[i]=0;|
000023d0 0a 20 20 20 20 68 6f 72 69 7a 6f 6e 32 5b 69 5d |. horizon2[i]|
000023e0 3d 31 30 32 34 3b 0a 20 20 7d 0a 20 20 62 62 63 |=1024;. }. bbc|
000023f0 5f 6f 72 69 67 69 6e 28 36 34 30 2c 30 29 3b 0a |_origin(640,0);.|
00002400 20 20 66 6f 72 20 28 69 3d 30 3b 20 69 3c 6e 3b | for (i=0; i<n;|
00002410 20 69 2b 2b 29 20 7b 0a 20 20 20 20 66 6f 72 20 | i++) {. for |
00002420 28 6a 3d 31 2c 20 70 6c 3d 2d 35 31 32 2c 20 76 |(j=1, pl=-512, v|
00002430 6c 20 3d 20 64 61 74 61 5b 30 2b 6e 2a 69 5d 3b |l = data[0+n*i];|
00002440 20 6a 3c 6e 3b 20 6a 2b 2b 2c 20 70 6c 3d 70 72 | j<n; j++, pl=pr|
00002450 2c 20 76 6c 3d 76 72 29 20 7b 0a 20 20 20 20 20 |, vl=vr) {. |
00002460 20 66 6f 72 20 28 70 72 3d 70 6c 2b 28 31 30 32 | for (pr=pl+(102|
00002470 34 3e 3e 6c 69 6d 69 74 29 2c 20 76 72 3d 64 61 |4>>limit), vr=da|
00002480 74 61 5b 6a 2b 6e 2a 69 5d 2c 20 6b 3d 70 6c 2c |ta[j+n*i], k=pl,|
00002490 20 76 3d 76 6c 3b 20 6b 3c 70 72 3b 20 6b 2b 3d | v=vl; k<pr; k+=|
000024a0 32 2c 20 76 2b 3d 28 76 72 2d 76 6c 29 3e 3e 28 |2, v+=(vr-vl)>>(|
000024b0 39 2d 6c 69 6d 69 74 29 29 20 7b 0a 20 20 20 20 |9-limit)) {. |
000024c0 20 20 20 20 76 32 20 3d 20 28 6d 75 6c 5f 66 72 | v2 = (mul_fr|
000024d0 61 63 31 36 28 76 2c 20 63 6f 73 29 2b 6d 75 6c |ac16(v, cos)+mul|
000024e0 5f 66 72 61 63 31 36 28 69 3c 3c 28 31 36 2d 6c |_frac16(i<<(16-l|
000024f0 69 6d 69 74 29 2c 20 73 69 6e 29 29 2a 33 2f 32 |imit), sin))*3/2|
00002500 35 36 3b 0a 20 20 20 20 20 20 20 20 69 66 20 28 |56;. if (|
00002510 69 3d 3d 30 29 20 68 6f 72 69 7a 6f 6e 32 5b 6b |i==0) horizon2[k|
00002520 2b 35 31 32 5d 3d 76 32 3b 0a 20 20 20 20 20 20 |+512]=v2;. |
00002530 20 20 69 66 20 28 68 6f 72 69 7a 6f 6e 5b 6b 2b | if (horizon[k+|
00002540 35 31 32 5d 3c 76 32 29 20 7b 0a 20 20 20 20 20 |512]<v2) {. |
00002550 20 20 20 20 20 62 62 63 5f 67 63 6f 6c 28 30 2c | bbc_gcol(0,|
00002560 20 34 2b 36 2a 69 2f 6e 29 3b 0a 20 20 20 20 20 | 4+6*i/n);. |
00002570 20 20 20 20 20 62 62 63 5f 70 6c 6f 74 28 36 39 | bbc_plot(69|
00002580 2c 20 6b 2c 20 76 32 29 3b 0a 20 20 20 20 20 20 |, k, v2);. |
00002590 20 20 20 20 68 6f 72 69 7a 6f 6e 5b 6b 2b 35 31 | horizon[k+51|
000025a0 32 5d 3d 76 32 3b 0a 20 20 20 20 20 20 20 20 7d |2]=v2;. }|
000025b0 0a 20 20 20 20 20 20 20 20 65 6c 73 65 20 69 66 |. else if|
000025c0 20 28 68 6f 72 69 7a 6f 6e 32 5b 6b 2b 35 31 32 | (horizon2[k+512|
000025d0 5d 3e 76 32 29 20 7b 0a 20 20 20 20 20 20 20 20 |]>v2) {. |
000025e0 20 20 62 62 63 5f 67 63 6f 6c 28 30 2c 20 31 30 | bbc_gcol(0, 10|
000025f0 2b 36 2a 69 2f 6e 29 3b 0a 20 20 20 20 20 20 20 |+6*i/n);. |
00002600 20 20 20 62 62 63 5f 70 6c 6f 74 28 36 39 2c 20 | bbc_plot(69, |
00002610 6b 2c 20 76 32 29 3b 0a 20 20 20 20 20 20 20 20 |k, v2);. |
00002620 20 20 68 6f 72 69 7a 6f 6e 32 5b 6b 2b 35 31 32 | horizon2[k+512|
00002630 5d 3d 76 32 3b 0a 20 20 20 20 20 20 20 20 7d 0a |]=v2;. }.|
00002640 20 20 20 20 20 20 7d 0a 20 20 20 20 7d 0a 20 20 | }. }. |
00002650 20 20 76 32 20 3d 20 28 6d 75 6c 5f 66 72 61 63 | v2 = (mul_frac|
00002660 31 36 28 76 72 2c 20 63 6f 73 29 2b 6d 75 6c 5f |16(vr, cos)+mul_|
00002670 66 72 61 63 31 36 28 69 3c 3c 28 31 36 2d 6c 69 |frac16(i<<(16-li|
00002680 6d 69 74 29 2c 20 73 69 6e 29 29 2a 33 2f 32 35 |mit), sin))*3/25|
00002690 36 3b 0a 20 20 20 20 69 66 20 28 69 3d 3d 30 29 |6;. if (i==0)|
000026a0 20 68 6f 72 69 7a 6f 6e 32 5b 6b 2b 35 31 32 5d | horizon2[k+512]|
000026b0 3d 76 32 3b 0a 20 20 20 20 69 66 20 28 68 6f 72 |=v2;. if (hor|
000026c0 69 7a 6f 6e 5b 6b 2b 35 31 32 5d 3c 76 32 29 20 |izon[k+512]<v2) |
000026d0 7b 0a 20 20 20 20 20 20 62 62 63 5f 67 63 6f 6c |{. bbc_gcol|
000026e0 28 30 2c 20 34 2b 36 2a 69 2f 6e 29 3b 0a 20 20 |(0, 4+6*i/n);. |
000026f0 20 20 20 20 62 62 63 5f 70 6c 6f 74 28 36 39 2c | bbc_plot(69,|
00002700 20 6b 2c 20 76 32 29 3b 0a 20 20 20 20 20 20 68 | k, v2);. h|
00002710 6f 72 69 7a 6f 6e 5b 6b 2b 35 31 32 5d 3d 76 32 |orizon[k+512]=v2|
00002720 3b 0a 20 20 20 20 7d 0a 20 20 20 20 65 6c 73 65 |;. }. else|
00002730 20 69 66 20 28 68 6f 72 69 7a 6f 6e 32 5b 6b 2b | if (horizon2[k+|
00002740 35 31 32 5d 3e 76 32 29 20 7b 0a 20 20 20 20 20 |512]>v2) {. |
00002750 20 62 62 63 5f 67 63 6f 6c 28 30 2c 20 31 30 2b | bbc_gcol(0, 10+|
00002760 36 2a 69 2f 6e 29 3b 0a 20 20 20 20 20 20 62 62 |6*i/n);. bb|
00002770 63 5f 70 6c 6f 74 28 36 39 2c 20 6b 2c 20 76 32 |c_plot(69, k, v2|
00002780 29 3b 0a 20 20 20 20 20 20 68 6f 72 69 7a 6f 6e |);. horizon|
00002790 32 5b 6b 2b 35 31 32 5d 3d 76 32 3b 0a 20 20 20 |2[k+512]=v2;. |
000027a0 20 7d 0a 20 20 7d 0a 20 20 66 6f 72 20 28 69 3d | }. }. for (i=|
000027b0 30 3b 20 69 3c 31 30 32 36 3b 20 69 2b 2b 29 20 |0; i<1026; i++) |
000027c0 7b 0a 20 20 20 20 69 66 20 28 68 6f 72 69 7a 6f |{. if (horizo|
000027d0 6e 5b 69 5d 20 3e 69 6d 61 67 65 5f 79 68 29 20 |n[i] >image_yh) |
000027e0 69 6d 61 67 65 5f 79 68 3d 68 6f 72 69 7a 6f 6e |image_yh=horizon|
000027f0 5b 69 5d 3b 0a 20 20 20 20 69 66 20 28 68 6f 72 |[i];. if (hor|
00002800 69 7a 6f 6e 32 5b 69 5d 3c 69 6d 61 67 65 5f 79 |izon2[i]<image_y|
00002810 6c 29 20 69 6d 61 67 65 5f 79 6c 3d 68 6f 72 69 |l) image_yl=hori|
00002820 7a 6f 6e 32 5b 69 5d 3b 0a 20 20 7d 0a 20 20 62 |zon2[i];. }. b|
00002830 62 63 5f 6d 6f 76 65 28 2d 35 31 32 2c 20 69 6d |bc_move(-512, im|
00002840 61 67 65 5f 79 6c 29 3b 0a 20 20 62 62 63 5f 6d |age_yl);. bbc_m|
00002850 6f 76 65 28 35 31 32 2c 20 69 6d 61 67 65 5f 79 |ove(512, image_y|
00002860 68 29 3b 0a 20 20 62 62 63 5f 70 6c 6f 74 28 31 |h);. bbc_plot(1|
00002870 38 39 2c 20 2d 35 31 32 2c 20 35 31 32 2b 28 69 |89, -512, 512+(i|
00002880 6d 61 67 65 5f 79 6c 2d 69 6d 61 67 65 5f 79 68 |mage_yl-image_yh|
00002890 29 2f 32 29 3b 0a 0a 20 20 72 65 74 75 72 6e 20 |)/2);.. return |
000028a0 54 52 55 45 3b 0a 7d 0a 0a 6d 6f 6e 69 74 6f 72 |TRUE;.}..monitor|
000028b0 20 72 65 61 64 5f 6d 6f 6e 69 74 6f 72 5f 74 79 | read_monitor_ty|
000028c0 70 65 28 76 6f 69 64 29 0a 7b 0a 20 20 69 6e 74 |pe(void).{. int|
000028d0 20 72 32 3b 0a 20 20 6f 73 5f 73 77 69 33 72 28 | r2;. os_swi3r(|
000028e0 36 2c 20 31 36 31 2c 31 33 33 2c 30 2c 20 30 2c |6, 161,133,0, 0,|
000028f0 30 2c 26 72 32 29 3b 0a 20 20 69 66 20 28 28 72 |0,&r2);. if ((r|
00002900 32 2f 34 20 26 20 33 29 20 3d 3d 20 31 29 20 72 |2/4 & 3) == 1) r|
00002910 65 74 75 72 6e 20 6d 75 6c 74 69 3b 0a 20 20 65 |eturn multi;. e|
00002920 6c 73 65 20 72 65 74 75 72 6e 20 6e 6f 6e 6d 75 |lse return nonmu|
00002930 6c 74 69 3b 0a 7d 0a 0a 69 6e 74 20 6d 61 69 6e |lti;.}..int main|
00002940 28 76 6f 69 64 29 0a 7b 0a 20 20 69 6e 74 20 73 |(void).{. int s|
00002950 65 65 64 2c 20 69 3b 0a 20 20 63 68 61 72 20 6b |eed, i;. char k|
00002960 65 79 3b 0a 20 20 64 6f 75 62 6c 65 20 64 69 6d |ey;. double dim|
00002970 3d 32 2e 32 3b 0a 20 20 64 6f 75 62 6c 65 20 64 |=2.2;. double d|
00002980 3b 0a 20 20 63 68 61 72 20 62 75 66 5b 31 30 32 |;. char buf[102|
00002990 34 5d 3b 0a 20 20 69 6e 74 20 74 6f 72 69 65 6e |4];. int torien|
000029a0 74 2c 20 74 73 63 72 6f 6c 6c 78 2c 20 74 73 63 |t, tscrollx, tsc|
000029b0 72 6f 6c 6c 79 3b 0a 20 20 6d 6f 6e 69 74 6f 72 |rolly;. monitor|
000029c0 74 79 70 65 20 3d 20 72 65 61 64 5f 6d 6f 6e 69 |type = read_moni|
000029d0 74 6f 72 5f 74 79 70 65 28 29 3b 0a 20 20 66 6c |tor_type();. fl|
000029e0 65 78 5f 69 6e 69 74 28 29 3b 0a 0a 20 20 66 6f |ex_init();.. fo|
000029f0 72 20 28 3b 3b 29 20 7b 0a 0a 20 20 20 20 62 62 |r (;;) {.. bb|
00002a00 63 5f 6d 6f 64 65 28 6d 6f 6e 69 74 6f 72 74 79 |c_mode(monitorty|
00002a10 70 65 3d 3d 6d 75 6c 74 69 20 3f 20 32 30 20 3a |pe==multi ? 20 :|
00002a20 20 31 32 29 3b 0a 20 20 20 20 70 72 69 6e 74 66 | 12);. printf|
00002a30 28 22 50 6c 65 61 73 65 20 73 65 6c 65 63 74 20 |("Please select |
00002a40 77 68 69 63 68 20 6c 61 6e 64 73 63 61 70 65 20 |which landscape |
00002a50 74 6f 20 63 6f 6e 73 74 72 75 63 74 20 28 61 6e |to construct (an|
00002a60 79 20 69 6e 74 65 67 65 72 29 20 22 29 3b 0a 20 |y integer) ");. |
00002a70 20 20 20 73 63 61 6e 66 28 22 25 69 22 2c 20 26 | scanf("%i", &|
00002a80 73 65 65 64 29 3b 0a 20 20 20 20 64 6f 20 7b 0a |seed);. do {.|
00002a90 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 50 6c | printf("Pl|
00002aa0 65 61 73 65 20 65 6e 74 65 72 20 73 63 61 70 65 |ease enter scape|
00002ab0 20 72 65 73 6f 6c 75 74 69 6f 6e 20 28 73 75 67 | resolution (sug|
00002ac0 67 65 73 74 20 35 2e 2e 2e 38 2c 20 74 68 6f 75 |gest 5...8, thou|
00002ad0 67 68 20 33 2e 2e 2e 39 20 70 65 72 6d 69 73 73 |gh 3...9 permiss|
00002ae0 69 62 6c 65 29 20 22 29 3b 0a 20 20 20 20 20 20 |ible) ");. |
00002af0 73 63 61 6e 66 28 22 25 69 22 2c 20 26 6c 69 6d |scanf("%i", &lim|
00002b00 69 74 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 |it);. } while|
00002b10 20 28 6c 69 6d 69 74 3c 33 20 7c 7c 20 6c 69 6d | (limit<3 || lim|
00002b20 69 74 3e 39 29 3b 0a 20 20 20 20 4e 20 3d 20 31 |it>9);. N = 1|
00002b30 3c 3c 6c 69 6d 69 74 3b 0a 20 20 20 20 64 6f 20 |<<limit;. do |
00002b40 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 |{. printf("|
00002b50 50 6c 65 61 73 65 20 65 6e 74 65 72 20 73 63 61 |Please enter sca|
00002b60 70 65 20 64 69 6d 65 6e 73 69 6f 6e 20 28 73 75 |pe dimension (su|
00002b70 67 67 65 73 74 20 32 2e 31 35 2c 20 32 2d 33 20 |ggest 2.15, 2-3 |
00002b80 70 65 72 6d 69 73 73 69 62 6c 65 29 20 22 29 3b |permissible) ");|
00002b90 0a 20 20 20 20 20 20 73 63 61 6e 66 28 22 25 4c |. scanf("%L|
00002ba0 66 22 2c 20 26 64 69 6d 29 3b 0a 20 20 20 20 7d |f", &dim);. }|
00002bb0 20 77 68 69 6c 65 20 28 64 69 6d 3c 32 20 7c 7c | while (dim<2 |||
00002bc0 20 64 69 6d 3e 33 29 3b 0a 20 20 20 20 64 6f 20 | dim>3);. do |
00002bd0 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 |{. printf("|
00002be0 50 6c 65 61 73 65 20 65 6e 74 65 72 20 66 72 65 |Please enter fre|
00002bf0 71 75 65 6e 63 79 20 72 61 6e 67 65 20 28 73 75 |quency range (su|
00002c00 67 67 65 73 74 20 25 69 2c 20 30 2e 2e 2e 25 69 |ggest %i, 0...%i|
00002c10 20 70 65 72 6d 69 73 73 69 62 6c 65 29 20 22 2c | permissible) ",|
00002c20 4e 2f 32 2c 4e 2f 32 29 3b 0a 20 20 20 20 20 20 |N/2,N/2);. |
00002c30 73 63 61 6e 66 28 22 25 69 22 2c 20 26 72 61 6e |scanf("%i", &ran|
00002c40 67 65 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 |ge);. } while|
00002c50 20 28 72 61 6e 67 65 3c 30 20 7c 7c 20 72 61 6e | (range<0 || ran|
00002c60 67 65 3e 4e 2f 32 29 3b 0a 20 20 20 20 64 6f 20 |ge>N/2);. do |
00002c70 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 |{. printf("|
00002c80 50 6c 65 61 73 65 20 65 6e 74 65 72 20 73 65 61 |Please enter sea|
00002c90 2d 6c 65 76 65 6c 20 28 73 75 67 67 65 73 74 20 |-level (suggest |
00002ca0 30 2c 20 30 2d 31 20 70 65 72 6d 69 73 73 69 62 |0, 0-1 permissib|
00002cb0 6c 65 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 |le) ");. sc|
00002cc0 61 6e 66 28 22 25 4c 66 22 2c 20 26 64 29 3b 0a |anf("%Lf", &d);.|
00002cd0 20 20 20 20 7d 20 77 68 69 6c 65 20 28 64 3c 30 | } while (d<0|
00002ce0 20 7c 7c 20 64 3e 31 29 3b 0a 20 20 20 20 73 65 | || d>1);. se|
00002cf0 61 6c 65 76 65 6c 20 3d 20 28 69 6e 74 29 28 64 |alevel = (int)(d|
00002d00 2a 6f 6e 65 29 3b 0a 20 20 20 20 64 6f 20 7b 0a |*one);. do {.|
00002d10 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 50 6c | printf("Pl|
00002d20 65 61 73 65 20 65 6e 74 65 72 20 70 6f 77 65 72 |ease enter power|
00002d30 20 73 63 61 6c 69 6e 67 20 6c 61 77 20 28 73 75 | scaling law (su|
00002d40 67 67 65 73 74 20 31 2c 20 6f 72 20 30 2e 33 2d |ggest 1, or 0.3-|
00002d50 33 2e 30 2c 20 74 68 6f 75 67 68 20 3e 3d 20 30 |3.0, though >= 0|
00002d60 20 70 65 72 6d 69 73 73 69 62 6c 65 29 20 22 29 | permissible) ")|
00002d70 3b 0a 20 20 20 20 20 20 73 63 61 6e 66 28 22 25 |;. scanf("%|
00002d80 4c 66 22 2c 20 26 64 29 3b 0a 20 20 20 20 7d 20 |Lf", &d);. } |
00002d90 77 68 69 6c 65 20 28 64 3c 30 29 3b 0a 20 20 20 |while (d<0);. |
00002da0 20 70 6f 77 73 63 61 6c 65 20 3d 20 28 69 6e 74 | powscale = (int|
00002db0 29 28 64 2a 6f 6e 65 29 3b 0a 20 20 20 20 64 6f |)(d*one);. do|
00002dc0 20 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 | {. printf(|
00002dd0 22 50 6c 65 61 73 65 20 65 6e 74 65 72 20 6c 69 |"Please enter li|
00002de0 6e 65 61 72 20 73 63 61 6c 69 6e 67 20 66 61 63 |near scaling fac|
00002df0 74 6f 72 20 28 73 75 67 67 65 73 74 20 31 2c 20 |tor (suggest 1, |
00002e00 30 2d 31 20 70 65 72 6d 69 73 73 69 62 6c 65 29 |0-1 permissible)|
00002e10 20 22 29 3b 0a 20 20 20 20 20 20 73 63 61 6e 66 | ");. scanf|
00002e20 28 22 25 4c 66 22 2c 20 26 64 29 3b 0a 20 20 20 |("%Lf", &d);. |
00002e30 20 7d 20 77 68 69 6c 65 20 28 64 3c 30 20 7c 7c | } while (d<0 |||
00002e40 20 64 3e 31 29 3b 0a 20 20 20 20 6c 69 6e 73 63 | d>1);. linsc|
00002e50 61 6c 65 20 3d 20 28 69 6e 74 29 28 64 2a 6f 6e |ale = (int)(d*on|
00002e60 65 29 3b 0a 20 20 20 20 64 6f 20 7b 0a 20 20 20 |e);. do {. |
00002e70 20 20 20 70 72 69 6e 74 66 28 22 50 6c 65 61 73 | printf("Pleas|
00002e80 65 20 65 6e 74 65 72 20 6f 72 69 65 6e 74 61 74 |e enter orientat|
00002e90 69 6f 6e 20 28 73 75 67 67 65 73 74 20 30 2c 20 |ion (suggest 0, |
00002ea0 30 2e 2e 2e 33 20 70 65 72 6d 69 73 73 69 62 6c |0...3 permissibl|
00002eb0 65 29 20 22 29 3b 0a 20 20 20 20 20 20 73 63 61 |e) ");. sca|
00002ec0 6e 66 28 22 25 69 22 2c 20 26 74 6f 72 69 65 6e |nf("%i", &torien|
00002ed0 74 29 3b 0a 20 20 20 20 7d 20 77 68 69 6c 65 20 |t);. } while |
00002ee0 28 74 6f 72 69 65 6e 74 3c 30 20 7c 7c 20 74 6f |(torient<0 || to|
00002ef0 72 69 65 6e 74 3e 33 29 3b 0a 20 20 20 20 64 6f |rient>3);. do|
00002f00 20 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 | {. printf(|
00002f10 22 50 6c 65 61 73 65 20 65 6e 74 65 72 20 78 20 |"Please enter x |
00002f20 73 63 72 6f 6c 6c 20 76 61 6c 75 65 20 28 73 75 |scroll value (su|
00002f30 67 67 65 73 74 20 30 2c 20 30 2e 2e 2e 35 31 31 |ggest 0, 0...511|
00002f40 20 70 65 72 6d 69 73 73 69 62 6c 65 29 20 22 29 | permissible) ")|
00002f50 3b 0a 20 20 20 20 20 20 73 63 61 6e 66 28 22 25 |;. scanf("%|
00002f60 69 22 2c 20 26 74 73 63 72 6f 6c 6c 78 29 3b 0a |i", &tscrollx);.|
00002f70 20 20 20 20 7d 20 77 68 69 6c 65 20 28 74 73 63 | } while (tsc|
00002f80 72 6f 6c 6c 78 3c 30 20 7c 7c 20 74 73 63 72 6f |rollx<0 || tscro|
00002f90 6c 6c 78 3e 3d 35 31 32 29 3b 0a 20 20 20 20 64 |llx>=512);. d|
00002fa0 6f 20 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 |o {. printf|
00002fb0 28 22 50 6c 65 61 73 65 20 65 6e 74 65 72 20 79 |("Please enter y|
00002fc0 20 73 63 72 6f 6c 6c 20 76 61 6c 75 65 20 28 73 | scroll value (s|
00002fd0 75 67 67 65 73 74 20 30 2c 20 30 2e 2e 2e 35 31 |uggest 0, 0...51|
00002fe0 31 20 70 65 72 6d 69 73 73 69 62 6c 65 29 20 22 |1 permissible) "|
00002ff0 29 3b 0a 20 20 20 20 20 20 73 63 61 6e 66 28 22 |);. scanf("|
00003000 25 69 22 2c 20 26 74 73 63 72 6f 6c 6c 79 29 3b |%i", &tscrolly);|
00003010 0a 20 20 20 20 7d 20 77 68 69 6c 65 20 28 74 73 |. } while (ts|
00003020 63 72 6f 6c 6c 79 3c 30 20 7c 7c 20 74 73 63 72 |crolly<0 || tscr|
00003030 6f 6c 6c 79 3e 3d 35 31 32 29 3b 0a 0a 20 20 20 |olly>=512);.. |
00003040 20 74 69 6c 74 20 3d 20 32 30 30 30 30 3b 09 2f | tilt = 20000;./|
00003050 2a 30 2d 36 35 35 33 36 20 70 65 72 6d 69 73 73 |*0-65536 permiss|
00003060 69 62 6c 65 2a 2f 0a 0a 20 20 20 20 69 66 20 28 |ible*/.. if (|
00003070 21 66 6c 65 78 5f 61 6c 6c 6f 63 28 28 66 6c 65 |!flex_alloc((fle|
00003080 78 5f 70 74 72 29 26 64 61 74 61 2c 20 4e 2a 4e |x_ptr)&data, N*N|
00003090 2a 73 69 7a 65 6f 66 28 63 6f 6d 70 6c 65 78 29 |*sizeof(complex)|
000030a0 29 29 20 77 65 72 72 28 31 2c 20 22 49 6e 73 75 |)) werr(1, "Insu|
000030b0 66 66 69 63 69 65 6e 74 20 6d 65 6d 6f 72 79 22 |fficient memory"|
000030c0 29 3b 0a 20 20 20 20 2f 2a 6e 6f 74 65 2c 20 66 |);. /*note, f|
000030d0 75 72 74 68 65 72 20 66 6c 65 78 20 63 61 6c 6c |urther flex call|
000030e0 73 20 6d 61 79 20 73 68 69 66 74 20 64 61 74 61 |s may shift data|
000030f0 2c 20 73 6f 20 62 65 20 63 61 72 65 66 75 6c 20 |, so be careful |
00003100 77 69 74 68 20 70 74 72 73 20 69 6e 74 6f 20 64 |with ptrs into d|
00003110 61 74 61 2a 2f 0a 20 20 20 20 62 62 63 5f 6d 6f |ata*/. bbc_mo|
00003120 64 65 28 6d 6f 6e 69 74 6f 72 74 79 70 65 3d 3d |de(monitortype==|
00003130 6d 75 6c 74 69 20 3f 20 32 30 20 3a 20 31 32 29 |multi ? 20 : 12)|
00003140 3b 0a 20 20 20 20 62 62 63 5f 76 64 75 71 28 32 |;. bbc_vduq(2|
00003150 38 2c 37 32 2c 33 31 2c 37 39 2c 30 29 3b 0a 20 |8,72,31,79,0);. |
00003160 20 20 20 62 62 63 5f 70 61 6c 65 74 74 65 28 31 | bbc_palette(1|
00003170 2c 20 31 36 2c 20 35 31 2c 35 31 2c 35 31 29 3b |, 16, 51,51,51);|
00003180 0a 20 20 20 20 62 62 63 5f 70 61 6c 65 74 74 65 |. bbc_palette|
00003190 28 32 2c 20 31 36 2c 20 32 35 35 2c 32 35 35 2c |(2, 16, 255,255,|
000031a0 32 35 35 29 3b 0a 20 20 20 20 62 62 63 5f 70 61 |255);. bbc_pa|
000031b0 6c 65 74 74 65 28 30 2c 20 31 36 2c 20 31 39 32 |lette(0, 16, 192|
000031c0 2c 31 39 32 2c 31 39 32 29 3b 0a 20 20 20 20 66 |,192,192);. f|
000031d0 6f 72 20 28 69 3d 30 3b 20 69 3c 36 3b 20 69 2b |or (i=0; i<6; i+|
000031e0 2b 29 20 7b 0a 20 20 20 20 20 20 62 62 63 5f 70 |+) {. bbc_p|
000031f0 61 6c 65 74 74 65 28 20 34 2b 69 2c 20 31 36 2c |alette( 4+i, 16,|
00003200 20 33 32 2b 32 30 2a 69 2c 31 33 36 2b 37 2a 69 | 32+20*i,136+7*i|
00003210 2c 33 32 2b 32 30 2a 69 29 3b 0a 20 20 20 20 20 |,32+20*i);. |
00003220 20 62 62 63 5f 70 61 6c 65 74 74 65 28 31 30 2b | bbc_palette(10+|
00003230 69 2c 20 31 36 2c 20 31 32 38 2b 38 2a 69 2c 39 |i, 16, 128+8*i,9|
00003240 36 2b 31 32 2a 69 2c 36 34 2b 31 36 2a 69 29 3b |6+12*i,64+16*i);|
00003250 0a 20 20 20 20 7d 0a 20 20 20 20 6c 69 6e 65 73 |. }. lines|
00003260 65 70 20 3d 20 6d 6f 6e 69 74 6f 72 74 79 70 65 |ep = monitortype|
00003270 3d 3d 6d 75 6c 74 69 20 3f 20 32 30 20 3a 20 34 |==multi ? 20 : 4|
00003280 30 3b 0a 20 20 20 20 70 72 69 6e 74 66 28 22 5c |0;. printf("\|
00003290 6e 50 6c 65 61 73 65 20 20 77 61 69 74 20 2e 2e |nPlease wait ..|
000032a0 2e 5c 6e 22 29 3b 0a 20 20 20 20 61 73 73 69 67 |.\n");. assig|
000032b0 6e 5f 66 74 28 64 61 74 61 2c 20 6c 69 6d 69 74 |n_ft(data, limit|
000032c0 2c 20 64 69 6d 2c 20 73 65 65 64 29 3b 0a 20 20 |, dim, seed);. |
000032d0 20 20 69 66 20 28 72 61 6e 67 65 3c 4e 2f 32 29 | if (range<N/2)|
000032e0 20 63 6c 69 70 5f 66 74 28 64 61 74 61 2c 20 6c | clip_ft(data, l|
000032f0 69 6d 69 74 2c 20 72 61 6e 67 65 29 3b 0a 20 20 |imit, range);. |
00003300 20 20 69 6e 76 66 66 74 5f 32 64 28 64 61 74 61 | invfft_2d(data|
00003310 2c 20 6c 69 6d 69 74 29 3b 0a 20 20 20 20 69 66 |, limit);. if|
00003320 20 28 21 6e 6f 72 6d 61 6c 69 73 65 5f 61 6e 64 | (!normalise_and|
00003330 5f 63 6f 6d 70 61 63 74 28 26 64 61 74 61 2c 20 |_compact(&data, |
00003340 6c 69 6d 69 74 29 29 20 77 65 72 72 28 31 2c 20 |limit)) werr(1, |
00003350 22 43 61 6e 27 74 20 63 6f 6d 70 61 63 74 20 64 |"Can't compact d|
00003360 61 74 61 22 29 3b 0a 20 20 20 20 6f 72 69 65 6e |ata");. orien|
00003370 74 3d 73 63 72 6f 6c 6c 78 3d 73 63 72 6f 6c 6c |t=scrollx=scroll|
00003380 79 3d 30 3b 0a 20 20 20 20 74 73 63 72 6f 6c 6c |y=0;. tscroll|
00003390 78 3e 3e 3d 39 2d 6c 69 6d 69 74 3b 0a 20 20 20 |x>>=9-limit;. |
000033a0 20 74 73 63 72 6f 6c 6c 79 3e 3e 3d 39 2d 6c 69 | tscrolly>>=9-li|
000033b0 6d 69 74 3b 0a 20 20 20 20 69 66 20 28 74 73 63 |mit;. if (tsc|
000033c0 72 6f 6c 6c 78 3e 30 29 20 73 63 72 6f 6c 6c 28 |rollx>0) scroll(|
000033d0 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 72 69 67 |data, limit, rig|
000033e0 68 74 2c 20 74 73 63 72 6f 6c 6c 78 29 3b 0a 20 |ht, tscrollx);. |
000033f0 20 20 20 69 66 20 28 74 73 63 72 6f 6c 6c 78 3c | if (tscrollx<|
00003400 30 29 20 73 63 72 6f 6c 6c 28 64 61 74 61 2c 20 |0) scroll(data, |
00003410 6c 69 6d 69 74 2c 20 6c 65 66 74 2c 20 2d 74 73 |limit, left, -ts|
00003420 63 72 6f 6c 6c 78 29 3b 0a 20 20 20 20 69 66 20 |crollx);. if |
00003430 28 74 73 63 72 6f 6c 6c 79 3e 30 29 20 73 63 72 |(tscrolly>0) scr|
00003440 6f 6c 6c 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c |oll(data, limit,|
00003450 20 75 70 2c 20 20 20 20 74 73 63 72 6f 6c 6c 79 | up, tscrolly|
00003460 29 3b 0a 20 20 20 20 69 66 20 28 74 73 63 72 6f |);. if (tscro|
00003470 6c 6c 79 3c 30 29 20 73 63 72 6f 6c 6c 28 64 61 |lly<0) scroll(da|
00003480 74 61 2c 20 6c 69 6d 69 74 2c 20 64 6f 77 6e 2c |ta, limit, down,|
00003490 20 2d 74 73 63 72 6f 6c 6c 79 29 3b 0a 20 20 20 | -tscrolly);. |
000034a0 20 73 77 69 74 63 68 20 28 74 6f 72 69 65 6e 74 | switch (torient|
000034b0 29 20 7b 0a 20 20 20 20 20 20 63 61 73 65 20 31 |) {. case 1|
000034c0 3a 0a 20 20 20 20 20 20 20 20 72 6f 74 61 74 65 |:. rotate|
000034d0 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 54 52 |(data, limit, TR|
000034e0 55 45 29 3b 0a 20 20 20 20 20 20 20 20 62 72 65 |UE);. bre|
000034f0 61 6b 3b 0a 20 20 20 20 20 20 63 61 73 65 20 32 |ak;. case 2|
00003500 3a 0a 20 20 20 20 20 20 20 20 72 6f 74 61 74 65 |:. rotate|
00003510 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 54 52 |(data, limit, TR|
00003520 55 45 29 3b 0a 20 20 20 20 20 20 20 20 72 6f 74 |UE);. rot|
00003530 61 74 65 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c |ate(data, limit,|
00003540 20 54 52 55 45 29 3b 0a 20 20 20 20 20 20 20 20 | TRUE);. |
00003550 62 72 65 61 6b 3b 0a 20 20 20 20 20 20 63 61 73 |break;. cas|
00003560 65 20 33 3a 0a 20 20 20 20 20 20 20 20 72 6f 74 |e 3:. rot|
00003570 61 74 65 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c |ate(data, limit,|
00003580 20 46 41 4c 53 45 29 3b 0a 20 20 20 20 20 20 20 | FALSE);. |
00003590 20 62 72 65 61 6b 3b 0a 20 20 20 20 7d 0a 20 20 | break;. }. |
000035a0 20 20 64 6f 20 7b 0a 20 20 20 20 20 20 62 62 63 | do {. bbc|
000035b0 5f 63 6c 67 28 29 3b 0a 20 20 20 20 20 20 67 72 |_clg();. gr|
000035c0 61 70 68 28 64 61 74 61 2c 20 6c 69 6d 69 74 29 |aph(data, limit)|
000035d0 3b 0a 20 20 20 20 20 20 62 62 63 5f 76 64 75 71 |;. bbc_vduq|
000035e0 28 32 38 2c 30 2c 33 31 2c 37 2c 30 29 3b 0a 20 |(28,0,31,7,0);. |
000035f0 20 20 20 20 20 70 72 69 6e 74 66 28 22 5c 6e 5c | printf("\n\|
00003600 6e 20 20 73 63 61 70 65 5c 6e 25 37 69 5c 6e 5c |n scape\n%7i\n\|
00003610 6e 20 20 20 20 72 65 73 5c 6e 25 37 69 5c 6e 5c |n res\n%7i\n\|
00003620 6e 20 20 20 20 64 69 6d 5c 6e 25 37 2e 32 66 5c |n dim\n%7.2f\|
00003630 6e 5c 6e 20 20 72 61 6e 67 65 5c 6e 25 37 69 5c |n\n range\n%7i\|
00003640 6e 5c 6e 20 20 20 20 73 65 61 5c 6e 25 37 2e 32 |n\n sea\n%7.2|
00003650 66 5c 6e 5c 6e 20 70 73 63 61 6c 65 5c 6e 25 37 |f\n\n pscale\n%7|
00003660 2e 32 66 5c 6e 5c 6e 20 6c 73 63 61 6c 65 5c 6e |.2f\n\n lscale\n|
00003670 25 37 2e 32 66 5c 6e 5c 6e 20 6f 72 69 65 6e 74 |%7.2f\n\n orient|
00003680 5c 6e 25 37 69 5c 6e 5c 6e 20 20 20 20 20 20 78 |\n%7i\n\n x|
00003690 5c 6e 25 37 69 5c 6e 5c 6e 20 20 20 20 20 20 79 |\n%7i\n\n y|
000036a0 5c 6e 25 37 69 22 2c 20 73 65 65 64 2c 20 6c 69 |\n%7i", seed, li|
000036b0 6d 69 74 2c 20 64 69 6d 2c 20 72 61 6e 67 65 2c |mit, dim, range,|
000036c0 20 28 69 6e 74 29 73 65 61 6c 65 76 65 6c 2f 36 | (int)sealevel/6|
000036d0 35 35 33 36 2e 30 2c 20 28 69 6e 74 29 70 6f 77 |5536.0, (int)pow|
000036e0 73 63 61 6c 65 2f 36 35 35 33 36 2e 30 2c 20 28 |scale/65536.0, (|
000036f0 69 6e 74 29 6c 69 6e 73 63 61 6c 65 2f 36 35 35 |int)linscale/655|
00003700 33 36 2e 30 2c 20 6f 72 69 65 6e 74 2c 20 73 63 |36.0, orient, sc|
00003710 72 6f 6c 6c 78 2c 20 73 63 72 6f 6c 6c 79 29 3b |rollx, scrolly);|
00003720 0a 20 20 20 20 20 20 62 62 63 5f 76 64 75 71 28 |. bbc_vduq(|
00003730 32 38 2c 37 32 2c 33 31 2c 37 39 2c 30 29 3b 0a |28,72,31,79,0);.|
00003740 20 20 20 20 20 20 62 62 63 5f 63 6c 73 28 29 3b | bbc_cls();|
00003750 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 5c |. printf("\|
00003760 6e 5c 6e 5c 6e 5c 6e 5c 6e 5c 6e 5c 6e 5c 6e 50 |n\n\n\n\n\n\n\nP|
00003770 72 65 73 73 20 61 20 6b 65 79 20 2e 2e 2e 5c 6e |ress a key ...\n|
00003780 5c 6e 22 29 3b 0a 20 20 20 20 20 20 70 72 69 6e |\n");. prin|
00003790 74 66 28 22 5b 20 27 3c 27 5c 6e 20 20 6f 72 5c |tf("[ '<'\n or\|
000037a0 6e 20 20 27 3e 27 5c 6e 20 20 74 6f 5c 6e 72 6f |n '>'\n to\nro|
000037b0 74 61 74 65 5c 6e 20 73 63 61 70 65 5c 6e 5c 6e |tate\n scape\n\n|
000037c0 63 75 72 73 6f 72 73 5c 6e 20 28 2b 5c 78 38 62 |cursors\n (+\x8b|
000037d0 29 5c 6e 20 20 74 6f 5c 6e 73 63 72 6f 6c 6c 5c |)\n to\nscroll\|
000037e0 6e 20 73 63 61 70 65 5c 6e 5c 6e 20 63 6f 70 79 |n scape\n\n copy|
000037f0 5c 6e 20 20 74 6f 5c 6e 20 73 61 76 65 5c 6e 20 |\n to\n save\n |
00003800 69 6d 61 67 65 5c 6e 20 20 20 20 20 20 20 5d 5c |image\n ]\|
00003810 6e 22 29 3b 0a 20 20 20 20 20 20 64 6f 20 7b 0a |n");. do {.|
00003820 20 20 20 20 20 20 20 20 6f 73 5f 63 6c 69 28 22 | os_cli("|
00003830 66 78 20 31 35 20 31 22 29 3b 0a 20 20 20 20 20 |fx 15 1");. |
00003840 20 20 20 62 62 63 5f 63 75 72 73 6f 72 28 30 29 | bbc_cursor(0)|
00003850 3b 0a 20 20 20 20 20 20 20 20 6f 73 5f 73 77 69 |;. os_swi|
00003860 32 28 4f 53 5f 42 79 74 65 2c 34 2c 31 29 3b 0a |2(OS_Byte,4,1);.|
00003870 20 20 20 20 20 20 20 20 6b 65 79 20 3d 20 62 62 | key = bb|
00003880 63 5f 67 65 74 28 29 3b 0a 20 20 20 20 20 20 20 |c_get();. |
00003890 20 6f 73 5f 73 77 69 32 28 4f 53 5f 42 79 74 65 | os_swi2(OS_Byte|
000038a0 2c 34 2c 30 29 3b 0a 20 20 20 20 20 20 20 20 69 |,4,0);. i|
000038b0 66 20 28 6b 65 79 3d 3d 31 33 35 29 20 7b 0a 20 |f (key==135) {. |
000038c0 20 20 20 20 20 20 20 20 20 73 63 72 65 65 6e 2b | screen+|
000038d0 3d 31 3b 0a 20 20 20 20 20 20 20 20 20 20 73 70 |=1;. sp|
000038e0 72 69 6e 74 66 28 62 75 66 2c 20 22 3c 53 63 61 |rintf(buf, "<Sca|
000038f0 70 65 31 24 44 69 72 3e 2e 5e 2e 49 6d 61 67 65 |pe1$Dir>.^.Image|
00003900 25 30 33 69 22 2c 20 73 63 72 65 65 6e 29 3b 0a |%03i", screen);.|
00003910 20 20 20 20 20 20 20 20 20 20 77 69 6d 70 74 5f | wimpt_|
00003920 63 6f 6d 70 6c 61 69 6e 28 73 70 72 69 74 65 5f |complain(sprite_|
00003930 73 63 72 65 65 6e 73 61 76 65 28 62 75 66 2c 20 |screensave(buf, |
00003940 54 52 55 45 29 29 3b 0a 20 20 20 20 20 20 20 20 |TRUE));. |
00003950 7d 0a 20 20 20 20 20 20 7d 20 77 68 69 6c 65 20 |}. } while |
00003960 28 6b 65 79 3d 3d 31 33 35 29 3b 0a 20 20 20 20 |(key==135);. |
00003970 20 20 69 66 20 28 6b 65 79 3d 3d 34 34 20 7c 7c | if (key==44 |||
00003980 20 6b 65 79 3d 3d 34 36 29 20 72 6f 74 61 74 65 | key==46) rotate|
00003990 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 6b 65 |(data, limit, ke|
000039a0 79 3d 3d 34 36 29 3b 0a 20 20 20 20 20 20 69 66 |y==46);. if|
000039b0 20 28 6b 65 79 3d 3d 31 33 36 20 7c 7c 20 6b 65 | (key==136 || ke|
000039c0 79 3d 3d 31 33 37 20 7c 7c 20 6b 65 79 3d 3d 31 |y==137 || key==1|
000039d0 33 38 20 7c 7c 20 6b 65 79 3d 3d 31 33 39 29 20 |38 || key==139) |
000039e0 7b 0a 20 20 20 20 20 20 20 20 69 3d 6c 69 6d 69 |{. i=limi|
000039f0 74 2d 32 3b 0a 20 20 20 20 20 20 20 20 69 66 20 |t-2;. if |
00003a00 28 61 6b 62 64 5f 70 6f 6c 6c 73 68 28 29 29 20 |(akbd_pollsh()) |
00003a10 69 2d 3d 32 3b 0a 20 20 20 20 20 20 20 20 69 66 |i-=2;. if|
00003a20 20 28 69 3c 30 29 20 69 3d 30 3b 0a 20 20 20 20 | (i<0) i=0;. |
00003a30 20 20 20 20 73 63 72 6f 6c 6c 28 64 61 74 61 2c | scroll(data,|
00003a40 20 6c 69 6d 69 74 2c 20 6b 65 79 2c 20 31 3c 3c | limit, key, 1<<|
00003a50 69 29 3b 0a 20 20 20 20 20 20 7d 0a 20 20 20 20 |i);. }. |
00003a60 7d 20 77 68 69 6c 65 20 28 6b 65 79 3d 3d 34 34 |} while (key==44|
00003a70 20 7c 7c 20 6b 65 79 3d 3d 34 36 20 7c 7c 20 6b | || key==46 || k|
00003a80 65 79 3d 3d 31 33 36 20 7c 7c 20 6b 65 79 3d 3d |ey==136 || key==|
00003a90 31 33 37 20 7c 7c 20 6b 65 79 3d 3d 31 33 38 20 |137 || key==138 |
00003aa0 7c 7c 20 6b 65 79 3d 3d 31 33 39 29 3b 0a 20 20 ||| key==139);. |
00003ab0 20 20 66 6c 65 78 5f 66 72 65 65 28 28 66 6c 65 | flex_free((fle|
00003ac0 78 5f 70 74 72 29 26 64 61 74 61 29 3b 0a 20 20 |x_ptr)&data);. |
00003ad0 20 20 62 62 63 5f 63 6c 73 28 29 3b 0a 0a 20 7d | bbc_cls();.. }|
00003ae0 0a 0a 20 20 72 65 74 75 72 6e 20 30 3b 0a 7d 0a |.. return 0;.}.|
00003af0 
.