Home » Archimedes archive » Archimedes World » AW-1996-03-Disc 2.adf » !AcornAns_AcornAns » !Trans2d/c/main
!Trans2d/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: | !Trans2d/c/main |
| Read OK: | ✔ |
| File size: | 2317 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 "swis.h"
#define one (1<<16)
#define dimsn(f) ((int)((f)-5)/5.0+2.0)
typedef enum {multi, nonmulti} monitor;
monitor monitortype;
int linesep;
typedef struct {
int re;
int im;
} complex;
int N;
int limit;
double T;
complex *data = 0;
double F;
char *fn_desc;
char *ft_desc;
char it_desc[256];
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 fft_2d(complex *data, int limit, int sign)
{
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 = sign<0 ? -2*one : 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 (l=sign>=0?limit-3:3, d1=data, d1e=d1+n*n; d1<d1e; d1++) {
d1->re >>= l;
d1->im >>= l;
}
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 = sign<0 ? -2*one : 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;
}
}
}
}
for (l=sign>=0?limit-3:3, d1=data, d1e=d1+n*n; d1<d1e; d1++) {
d1->re >>= l;
d1->im >>= l;
}
return TRUE;
}
BOOL filter(complex *data, int limit, double dimension)
{
int n=1<<limit;
int i, j, k2, l, m;
int mbpobt = (int)(-(4-dimension)*one);
int scale;
complex *d;
for (i=0, d=data; i<n; i++) {
k2 = i<n/2 ? i : n-i;
k2 = k2*k2;
for (j=0; j<n; j++, d++) {
if (i==0 && j==0) scale=0;
else {
l = j<n/2 ? j : n-j;
/*want m^(mbpobt/2), but range on m too great for pow16, so tweak m to allow application of
sqrt_frac16, then raise to power mbpobt*/
m = k2+l*l;
scale = pow16(sqrt_frac16(m*(one/4))*2, mbpobt);
}
d->re = mul_frac16(d->re, scale);
d->im = mul_frac16(d->im, scale);
}
}
return TRUE;
}
double assign(complex *data, int limit, int function)
{
int n=1<<limit;
int i, j;
int v;
double T;
double dim;
complex *d1;
switch (function)
{
case 0:
T = 1;
fn_desc = "sine wave packet";
ft_desc = "";
*it_desc = 0;
for (i=0, d1=data; i<n; i++) {
v = sin16(div_frac16(4*i, n));
for (j=0; j<n; j++, d1++) {
d1->re = 16*mul_frac16c(sin16(div_frac16(4*j, n)), v);
d1->im = 0;
}
}
break;
case 1:
T = 5;
fn_desc = "damped sine wave";
ft_desc = "";
*it_desc = 0;
for (i=0, d1=data; i<n; i++) {
v = i==0 ? one : mul_frac16c(div_frac16(n, i), sin16(div_frac16(4*i, n)));
for (j=0; j<n; j++, d1++) {
d1->re = mul_frac16c(j==0 ? one : mul_frac16c(div_frac16(n, j), sin16(div_frac16(4*j, n))), v);
d1->im = 0;
}
}
break;
case 2:
T = 1;
fn_desc = "square wave packet";
ft_desc = "";
*it_desc = 0;
for (i=0, d1=data; i<n; i++) {
for (j=0; j<n; j++, d1++) {
d1->re = i<n/2 && j<n/2 ? 16*one : 0;
d1->im = 0;
}
}
break;
case 3:
T = 1;
fn_desc = "white noise";
ft_desc = "";
*it_desc = 0;
for (i=0, d1=data; i<n; i++) {
for (j=0; j<n; j++, d1++) {
d1->re = 16*((rand()>>15)-one/2);
d1->im = 0;
}
}
break;
case 4:
T = 1;
fn_desc = "gaussian white noise";
ft_desc = "";
*it_desc = 0;
for (i=0, d1=data; i<n; i++) {
for (j=0; j<n; j++, d1++) {
d1->re = 16*gauss16();
d1->im = 0;
}
}
break;
default: /*for 6 cases, 5-10*/
T = 1;
fn_desc = "gaussian white noise";
ft_desc = "Filtered by 1/(f^k)";
dim = dimsn(function);
sprintf(it_desc, "k = %.3f, Fractal dim. = %.3f", 4.0-dim, dim);
for (i=0, d1=data; i<n; i++) {
for (j=0; j<n; j++, d1++) {
d1->re = 16*gauss16();
d1->im = 0;
}
}
break;
}
return T;
}
#define reorder(i, n, flag) ((flag) ? ((i)<(n)/2 ? (i)+(n)/2 : (i)-(n)/2) : (i))
#define yscale 124
BOOL graph(complex *data, int limit, int yo, BOOL repos, char *text, ...)
{
int n=1<<limit;
int i, j;
int min=1<<30;
int max=-min;
int scale;
complex *d1, *d1e;
va_list argp;
for (d1=data, d1e=d1+n*n; d1<d1e; d1++) {
if (d1->re<min) min=d1->re;
if (d1->im<min) min=d1->im;
if (d1->re>max) max=d1->re;
if (d1->im>max) max=d1->im;
}
max = max>0 ? max : -max;
min = min>0 ? min : -min;
scale = max>min ? max : min;
if (scale==0) scale=1;
bbc_origin(640,yo);
bbc_gcol(0, 1);
bbc_rectanglefill(-517,-4,1031,9);
bbc_rectanglefill(repos?-5:-517,-yscale-4,9,2*yscale+7);
bbc_vdu(5);
bbc_move((repos?0:-512)-44,yscale+40);
printf("%.3f", (double)scale/one);
bbc_move(-512+96,yscale+40);
va_start(argp, text);
vprintf(text, argp);
va_end(argp);
bbc_move((repos?0:-512)-60,-yscale-20);
printf("%.3f", -(double)scale/one);
bbc_move(-512-128,linesep/4);
printf("%.3f", repos?-F/2:0.0);
bbc_move(520,linesep/4);
printf("%.3f", repos?F/2:T);
bbc_vdu(4);
for (i=0; i<n; i++) {
bbc_gcol(0, 4+12*i/n);
bbc_move(-512,data[reorder(0, n, repos)+n*reorder(i, n, repos)].re*yscale/scale);
for (j=1; j<n; j++) {
bbc_plot(13, j*1024/n-512,data[reorder(j, n, repos)+n*reorder(i, n, repos)].re*yscale/scale);
}
bbc_gcol(0, 3);
bbc_move(-512,data[reorder(0, n, repos)+n*reorder(i, n, repos)].im*yscale/scale);
for (j=1; j<n; j++) {
bbc_plot(13, j*1024/n-512,data[reorder(j, n, repos)+n*reorder(i, n, repos)].im*yscale/scale);
}
}
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 function, seed, i;
char dummy;
monitortype = read_monitor_type();
for (;;) {
seed = -clock();
do {
printf("Please enter power (suggest 5-8, though 3-9 permissible) ");
scanf("%i", &limit);
} while (limit<3 || limit>9);
N = 1<<limit;
data = calloc(N*N, sizeof(complex));
if (data==0) werr(1, "Insufficient memory");
for (function=0; function<11; function++) {
sgauss16(seed);
srand(seed);
bbc_mode(monitortype==multi ? 20 : 12);
bbc_palette(1, 16, 51,51,51);
bbc_palette(2, 16, 255,255,255);
for (i=0; i<12; i++) bbc_palette(4+i, 16, 0,68+17*i,0);
linesep = monitortype==multi ? 20 : 40;
T=assign(data, limit, function);
F=N/T;
graph(data, limit, 852, FALSE, "Source function '%s' (Sample size N=%i)", fn_desc, N);
fft_2d(data, limit, 1);
if (function>4) filter(data, limit, dimsn(function));
graph(data, limit, 512, TRUE, "Fourier Transform %s", ft_desc);
fft_2d(data, limit, -1);
graph(data, limit, 172, FALSE, "Inverted Transform %s", it_desc);
bbc_vdu(5);
bbc_origin(0,0);
bbc_move(1152, 131);
printf("Press a");
bbc_move(1152, 131-linesep);
printf("key ...");
bbc_vdu(4);
os_cli("fx 15 1");
bbc_cursor(0);
dummy = bbc_get();
}
free(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 0a 23 69 6e 63 6c |de "os.h"..#incl|
000000c0 75 64 65 20 22 73 77 69 73 2e 68 22 0a 23 64 65 |ude "swis.h".#de|
000000d0 66 69 6e 65 20 6f 6e 65 20 28 31 3c 3c 31 36 29 |fine one (1<<16)|
000000e0 0a 0a 23 64 65 66 69 6e 65 20 64 69 6d 73 6e 28 |..#define dimsn(|
000000f0 66 29 20 28 28 69 6e 74 29 28 28 66 29 2d 35 29 |f) ((int)((f)-5)|
00000100 2f 35 2e 30 2b 32 2e 30 29 0a 0a 74 79 70 65 64 |/5.0+2.0)..typed|
00000110 65 66 20 65 6e 75 6d 20 7b 6d 75 6c 74 69 2c 20 |ef enum {multi, |
00000120 6e 6f 6e 6d 75 6c 74 69 7d 20 6d 6f 6e 69 74 6f |nonmulti} monito|
00000130 72 3b 0a 6d 6f 6e 69 74 6f 72 20 6d 6f 6e 69 74 |r;.monitor monit|
00000140 6f 72 74 79 70 65 3b 0a 69 6e 74 20 6c 69 6e 65 |ortype;.int line|
00000150 73 65 70 3b 0a 0a 74 79 70 65 64 65 66 20 73 74 |sep;..typedef st|
00000160 72 75 63 74 20 7b 0a 20 20 69 6e 74 20 72 65 3b |ruct {. int re;|
00000170 0a 20 20 69 6e 74 20 69 6d 3b 0a 7d 20 63 6f 6d |. int im;.} com|
00000180 70 6c 65 78 3b 0a 0a 69 6e 74 20 4e 3b 0a 69 6e |plex;..int N;.in|
00000190 74 20 6c 69 6d 69 74 3b 0a 64 6f 75 62 6c 65 20 |t limit;.double |
000001a0 54 3b 0a 63 6f 6d 70 6c 65 78 20 2a 64 61 74 61 |T;.complex *data|
000001b0 20 3d 20 30 3b 0a 64 6f 75 62 6c 65 20 46 3b 0a | = 0;.double F;.|
000001c0 0a 63 68 61 72 20 2a 66 6e 5f 64 65 73 63 3b 0a |.char *fn_desc;.|
000001d0 63 68 61 72 20 2a 66 74 5f 64 65 73 63 3b 0a 63 |char *ft_desc;.c|
000001e0 68 61 72 20 69 74 5f 64 65 73 63 5b 32 35 36 5d |har it_desc[256]|
000001f0 3b 0a 0a 65 78 74 65 72 6e 20 69 6e 74 20 72 62 |;..extern int rb|
00000200 62 63 69 6e 63 28 69 6e 74 20 72 2c 20 69 6e 74 |bcinc(int r, int|
00000210 20 6b 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 | k);.extern int |
00000220 64 69 76 5f 66 72 61 63 31 36 28 69 6e 74 20 6e |div_frac16(int n|
00000230 75 6d 62 65 72 2c 20 69 6e 74 20 64 69 76 69 73 |umber, int divis|
00000240 6f 72 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 |or);.extern int |
00000250 6d 75 6c 5f 66 72 61 63 31 36 28 69 6e 74 20 78 |mul_frac16(int x|
00000260 2c 20 69 6e 74 20 61 29 3b 0a 65 78 74 65 72 6e |, int a);.extern|
00000270 20 69 6e 74 20 6d 75 6c 5f 66 72 61 63 31 36 63 | int mul_frac16c|
00000280 28 69 6e 74 20 78 2c 20 69 6e 74 20 61 29 3b 0a |(int x, int a);.|
00000290 65 78 74 65 72 6e 20 69 6e 74 20 73 71 72 74 5f |extern int sqrt_|
000002a0 66 72 61 63 31 36 28 75 6e 73 69 67 6e 65 64 20 |frac16(unsigned |
000002b0 69 6e 74 20 78 29 3b 0a 65 78 74 65 72 6e 20 69 |int x);.extern i|
000002c0 6e 74 20 67 61 75 73 73 31 36 28 76 6f 69 64 29 |nt gauss16(void)|
000002d0 3b 0a 65 78 74 65 72 6e 20 76 6f 69 64 20 73 67 |;.extern void sg|
000002e0 61 75 73 73 31 36 28 69 6e 74 20 73 65 65 64 29 |auss16(int seed)|
000002f0 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 63 6f 73 |;.extern int cos|
00000300 31 36 28 69 6e 74 20 61 29 3b 0a 65 78 74 65 72 |16(int a);.exter|
00000310 6e 20 69 6e 74 20 73 69 6e 31 36 28 69 6e 74 20 |n int sin16(int |
00000320 61 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 65 |a);.extern int e|
00000330 78 70 31 36 28 69 6e 74 20 61 29 3b 0a 65 78 74 |xp16(int a);.ext|
00000340 65 72 6e 20 69 6e 74 20 6c 6e 31 36 28 69 6e 74 |ern int ln16(int|
00000350 20 61 29 3b 0a 65 78 74 65 72 6e 20 69 6e 74 20 | a);.extern int |
00000360 70 6f 77 31 36 28 69 6e 74 20 61 2c 20 69 6e 74 |pow16(int a, int|
00000370 20 62 29 3b 0a 0a 42 4f 4f 4c 20 66 66 74 5f 32 | b);..BOOL fft_2|
00000380 64 28 63 6f 6d 70 6c 65 78 20 2a 64 61 74 61 2c |d(complex *data,|
00000390 20 69 6e 74 20 6c 69 6d 69 74 2c 20 69 6e 74 20 | int limit, int |
000003a0 73 69 67 6e 29 0a 7b 0a 20 20 69 6e 74 20 69 2c |sign).{. int i,|
000003b0 20 6a 3b 0a 20 20 69 6e 74 20 6e 20 3d 20 31 3c | j;. int n = 1<|
000003c0 3c 6c 69 6d 69 74 3b 0a 20 20 63 6f 6d 70 6c 65 |<limit;. comple|
000003d0 78 20 74 31 2c 20 2a 64 31 2c 20 2a 64 31 65 2c |x t1, *d1, *d1e,|
000003e0 20 2a 64 32 3b 0a 20 20 69 6e 74 20 6b 2c 20 6c | *d2;. int k, l|
000003f0 3b 0a 20 20 69 6e 74 20 6c 6f 6f 70 3b 0a 20 20 |;. int loop;. |
00000400 69 6e 74 20 72 6f 74 61 6e 67 6c 65 2c 20 61 6e |int rotangle, an|
00000410 67 6c 65 3b 0a 20 20 69 6e 74 20 63 6f 73 2c 20 |gle;. int cos, |
00000420 73 69 6e 3b 0a 20 20 69 66 20 28 6c 69 6d 69 74 |sin;. if (limit|
00000430 3c 33 20 7c 7c 20 6c 69 6d 69 74 3e 33 32 29 20 |<3 || limit>32) |
00000440 72 65 74 75 72 6e 20 46 41 4c 53 45 3b 0a 0a 20 |return FALSE;.. |
00000450 20 66 6f 72 20 28 69 3d 6a 3d 30 3b 20 6a 3c 6e | for (i=j=0; j<n|
00000460 3b 20 69 3d 72 62 62 63 69 6e 63 28 69 2c 20 6c |; i=rbbcinc(i, l|
00000470 69 6d 69 74 29 2c 20 6a 2b 2b 29 20 69 66 20 28 |imit), j++) if (|
00000480 69 3e 6a 29 20 7b 0a 20 20 20 20 66 6f 72 20 28 |i>j) {. for (|
00000490 64 31 3d 64 61 74 61 2b 28 69 3c 3c 6c 69 6d 69 |d1=data+(i<<limi|
000004a0 74 29 2c 20 64 32 3d 64 61 74 61 2b 28 6a 3c 3c |t), d2=data+(j<<|
000004b0 6c 69 6d 69 74 29 2c 20 64 31 65 3d 64 31 2b 6e |limit), d1e=d1+n|
000004c0 3b 20 64 31 3c 64 31 65 3b 29 20 7b 0a 20 20 20 |; d1<d1e;) {. |
000004d0 20 20 20 74 31 20 3d 20 2a 64 31 3b 0a 20 20 20 | t1 = *d1;. |
000004e0 20 20 20 2a 64 31 2b 2b 20 3d 20 2a 64 32 3b 0a | *d1++ = *d2;.|
000004f0 20 20 20 20 20 20 2a 64 32 2b 2b 20 3d 20 74 31 | *d2++ = t1|
00000500 3b 0a 20 20 20 20 7d 0a 20 20 7d 0a 20 20 72 6f |;. }. }. ro|
00000510 74 61 6e 67 6c 65 20 3d 20 73 69 67 6e 3c 30 20 |tangle = sign<0 |
00000520 3f 20 2d 32 2a 6f 6e 65 20 3a 20 32 2a 6f 6e 65 |? -2*one : 2*one|
00000530 3b 0a 20 20 66 6f 72 20 28 6b 3d 6c 6f 6f 70 3d |;. for (k=loop=|
00000540 31 3b 20 6c 6f 6f 70 3c 3d 6c 69 6d 69 74 3b 20 |1; loop<=limit; |
00000550 6c 6f 6f 70 2b 2b 2c 20 6b 3c 3c 3d 31 2c 20 72 |loop++, k<<=1, r|
00000560 6f 74 61 6e 67 6c 65 2f 3d 32 29 20 7b 0a 20 20 |otangle/=2) {. |
00000570 20 20 66 6f 72 20 28 61 6e 67 6c 65 3d 30 2c 20 | for (angle=0, |
00000580 73 69 6e 3d 30 2c 20 63 6f 73 3d 6f 6e 65 2c 20 |sin=0, cos=one, |
00000590 69 3d 31 3b 20 69 3c 3d 6b 3b 20 69 2b 2b 2c 20 |i=1; i<=k; i++, |
000005a0 61 6e 67 6c 65 2b 3d 72 6f 74 61 6e 67 6c 65 2c |angle+=rotangle,|
000005b0 20 73 69 6e 3d 73 69 6e 31 36 28 61 6e 67 6c 65 | sin=sin16(angle|
000005c0 29 2c 20 63 6f 73 3d 63 6f 73 31 36 28 61 6e 67 |), cos=cos16(ang|
000005d0 6c 65 29 29 20 7b 0a 20 20 20 20 20 20 66 6f 72 |le)) {. for|
000005e0 20 28 6a 3d 69 2d 31 3b 20 6a 3c 6e 3b 20 6a 2b | (j=i-1; j<n; j+|
000005f0 3d 6b 3c 3c 31 29 20 7b 0a 20 20 20 20 20 20 20 |=k<<1) {. |
00000600 20 6c 20 3d 20 6a 2b 6b 3b 0a 20 20 20 20 20 20 | l = j+k;. |
00000610 20 20 66 6f 72 20 28 64 31 3d 64 61 74 61 2b 28 | for (d1=data+(|
00000620 6c 3c 3c 6c 69 6d 69 74 29 2c 20 64 32 3d 64 61 |l<<limit), d2=da|
00000630 74 61 2b 28 6a 3c 3c 6c 69 6d 69 74 29 2c 20 64 |ta+(j<<limit), d|
00000640 31 65 3d 64 31 2b 6e 3b 20 64 31 3c 64 31 65 3b |1e=d1+n; d1<d1e;|
00000650 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 74 31 |) {. t1|
00000660 2e 72 65 20 3d 20 28 6d 75 6c 5f 66 72 61 63 31 |.re = (mul_frac1|
00000670 36 63 28 64 31 2d 3e 72 65 2c 20 63 6f 73 29 2d |6c(d1->re, cos)-|
00000680 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e |mul_frac16c(d1->|
00000690 69 6d 2c 20 73 69 6e 29 29 3b 0a 20 20 20 20 20 |im, sin));. |
000006a0 20 20 20 20 20 74 31 2e 69 6d 20 3d 20 28 6d 75 | t1.im = (mu|
000006b0 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e 72 65 |l_frac16c(d1->re|
000006c0 2c 20 73 69 6e 29 2b 6d 75 6c 5f 66 72 61 63 31 |, sin)+mul_frac1|
000006d0 36 63 28 64 31 2d 3e 69 6d 2c 20 63 6f 73 29 29 |6c(d1->im, cos))|
000006e0 3b 0a 20 20 20 20 20 20 20 20 20 20 64 31 2d 3e |;. d1->|
000006f0 72 65 20 20 20 20 3d 20 64 32 2d 3e 72 65 2d 74 |re = d2->re-t|
00000700 31 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 20 20 |1.re;. |
00000710 64 31 2b 2b 2d 3e 69 6d 20 20 3d 20 64 32 2d 3e |d1++->im = d2->|
00000720 69 6d 2d 74 31 2e 69 6d 3b 0a 20 20 20 20 20 20 |im-t1.im;. |
00000730 20 20 20 20 64 32 2d 3e 72 65 20 20 20 2b 3d 20 | d2->re += |
00000740 74 31 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 20 |t1.re;. |
00000750 20 64 32 2b 2b 2d 3e 69 6d 20 2b 3d 20 74 31 2e | d2++->im += t1.|
00000760 69 6d 3b 0a 20 20 20 20 20 20 20 20 7d 0a 20 20 |im;. }. |
00000770 20 20 20 20 7d 0a 20 20 20 20 7d 0a 20 20 7d 0a | }. }. }.|
00000780 0a 20 20 66 6f 72 20 28 6c 3d 73 69 67 6e 3e 3d |. for (l=sign>=|
00000790 30 3f 6c 69 6d 69 74 2d 33 3a 33 2c 20 64 31 3d |0?limit-3:3, d1=|
000007a0 64 61 74 61 2c 20 64 31 65 3d 64 31 2b 6e 2a 6e |data, d1e=d1+n*n|
000007b0 3b 20 64 31 3c 64 31 65 3b 20 64 31 2b 2b 29 20 |; d1<d1e; d1++) |
000007c0 7b 0a 20 20 20 20 64 31 2d 3e 72 65 20 3e 3e 3d |{. d1->re >>=|
000007d0 20 6c 3b 0a 20 20 20 20 64 31 2d 3e 69 6d 20 3e | l;. d1->im >|
000007e0 3e 3d 20 6c 3b 0a 20 20 7d 0a 0a 20 20 66 6f 72 |>= l;. }.. for|
000007f0 20 28 69 3d 6a 3d 30 3b 20 6a 3c 6e 3b 20 69 3d | (i=j=0; j<n; i=|
00000800 72 62 62 63 69 6e 63 28 69 2c 20 6c 69 6d 69 74 |rbbcinc(i, limit|
00000810 29 2c 20 6a 2b 2b 29 20 69 66 20 28 69 3e 6a 29 |), j++) if (i>j)|
00000820 20 7b 0a 20 20 20 20 66 6f 72 20 28 64 31 3d 64 | {. for (d1=d|
00000830 61 74 61 2b 69 2c 20 64 32 3d 64 61 74 61 2b 6a |ata+i, d2=data+j|
00000840 2c 20 64 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 |, d1e=d1+n*n; d1|
00000850 3c 64 31 65 3b 20 64 31 2b 3d 6e 2c 20 64 32 2b |<d1e; d1+=n, d2+|
00000860 3d 6e 29 20 7b 0a 20 20 20 20 20 20 74 31 20 3d |=n) {. t1 =|
00000870 20 2a 64 31 3b 0a 20 20 20 20 20 20 2a 64 31 20 | *d1;. *d1 |
00000880 3d 20 2a 64 32 3b 0a 20 20 20 20 20 20 2a 64 32 |= *d2;. *d2|
00000890 20 3d 20 74 31 3b 0a 20 20 20 20 7d 0a 20 20 7d | = t1;. }. }|
000008a0 0a 20 20 72 6f 74 61 6e 67 6c 65 20 3d 20 73 69 |. rotangle = si|
000008b0 67 6e 3c 30 20 3f 20 2d 32 2a 6f 6e 65 20 3a 20 |gn<0 ? -2*one : |
000008c0 32 2a 6f 6e 65 3b 0a 20 20 66 6f 72 20 28 6b 3d |2*one;. for (k=|
000008d0 6c 6f 6f 70 3d 31 3b 20 6c 6f 6f 70 3c 3d 6c 69 |loop=1; loop<=li|
000008e0 6d 69 74 3b 20 6c 6f 6f 70 2b 2b 2c 20 6b 3c 3c |mit; loop++, k<<|
000008f0 3d 31 2c 20 72 6f 74 61 6e 67 6c 65 2f 3d 32 29 |=1, rotangle/=2)|
00000900 20 7b 0a 20 20 20 20 66 6f 72 20 28 61 6e 67 6c | {. for (angl|
00000910 65 3d 30 2c 20 73 69 6e 3d 30 2c 20 63 6f 73 3d |e=0, sin=0, cos=|
00000920 6f 6e 65 2c 20 69 3d 31 3b 20 69 3c 3d 6b 3b 20 |one, i=1; i<=k; |
00000930 69 2b 2b 2c 20 61 6e 67 6c 65 2b 3d 72 6f 74 61 |i++, angle+=rota|
00000940 6e 67 6c 65 2c 20 73 69 6e 3d 73 69 6e 31 36 28 |ngle, sin=sin16(|
00000950 61 6e 67 6c 65 29 2c 20 63 6f 73 3d 63 6f 73 31 |angle), cos=cos1|
00000960 36 28 61 6e 67 6c 65 29 29 20 7b 0a 20 20 20 20 |6(angle)) {. |
00000970 20 20 66 6f 72 20 28 6a 3d 69 2d 31 3b 20 6a 3c | for (j=i-1; j<|
00000980 6e 3b 20 6a 2b 3d 6b 3c 3c 31 29 20 7b 0a 20 20 |n; j+=k<<1) {. |
00000990 20 20 20 20 20 20 6c 20 3d 20 6a 2b 6b 3b 0a 20 | l = j+k;. |
000009a0 20 20 20 20 20 20 20 66 6f 72 20 28 64 31 3d 64 | for (d1=d|
000009b0 61 74 61 2b 6c 2c 20 64 32 3d 64 61 74 61 2b 6a |ata+l, d2=data+j|
000009c0 2c 20 64 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 |, d1e=d1+n*n; d1|
000009d0 3c 64 31 65 3b 20 64 31 2b 3d 6e 2c 20 64 32 2b |<d1e; d1+=n, d2+|
000009e0 3d 6e 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 |=n) {. |
000009f0 74 31 2e 72 65 20 3d 20 28 6d 75 6c 5f 66 72 61 |t1.re = (mul_fra|
00000a00 63 31 36 63 28 64 31 2d 3e 72 65 2c 20 63 6f 73 |c16c(d1->re, cos|
00000a10 29 2d 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 |)-mul_frac16c(d1|
00000a20 2d 3e 69 6d 2c 20 73 69 6e 29 29 3b 0a 20 20 20 |->im, sin));. |
00000a30 20 20 20 20 20 20 20 74 31 2e 69 6d 20 3d 20 28 | t1.im = (|
00000a40 6d 75 6c 5f 66 72 61 63 31 36 63 28 64 31 2d 3e |mul_frac16c(d1->|
00000a50 72 65 2c 20 73 69 6e 29 2b 6d 75 6c 5f 66 72 61 |re, sin)+mul_fra|
00000a60 63 31 36 63 28 64 31 2d 3e 69 6d 2c 20 63 6f 73 |c16c(d1->im, cos|
00000a70 29 29 3b 0a 20 20 20 20 20 20 20 20 20 20 64 31 |));. d1|
00000a80 2d 3e 72 65 20 20 3d 20 64 32 2d 3e 72 65 2d 74 |->re = d2->re-t|
00000a90 31 2e 72 65 3b 0a 20 20 20 20 20 20 20 20 20 20 |1.re;. |
00000aa0 64 31 2d 3e 69 6d 20 20 3d 20 64 32 2d 3e 69 6d |d1->im = d2->im|
00000ab0 2d 74 31 2e 69 6d 3b 0a 20 20 20 20 20 20 20 20 |-t1.im;. |
00000ac0 20 20 64 32 2d 3e 72 65 20 2b 3d 20 74 31 2e 72 | d2->re += t1.r|
00000ad0 65 3b 0a 20 20 20 20 20 20 20 20 20 20 64 32 2d |e;. d2-|
00000ae0 3e 69 6d 20 2b 3d 20 74 31 2e 69 6d 3b 0a 20 20 |>im += t1.im;. |
00000af0 20 20 20 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a | }. }.|
00000b00 20 20 20 20 7d 0a 20 20 7d 0a 0a 20 20 66 6f 72 | }. }.. for|
00000b10 20 28 6c 3d 73 69 67 6e 3e 3d 30 3f 6c 69 6d 69 | (l=sign>=0?limi|
00000b20 74 2d 33 3a 33 2c 20 64 31 3d 64 61 74 61 2c 20 |t-3:3, d1=data, |
00000b30 64 31 65 3d 64 31 2b 6e 2a 6e 3b 20 64 31 3c 64 |d1e=d1+n*n; d1<d|
00000b40 31 65 3b 20 64 31 2b 2b 29 20 7b 0a 20 20 20 20 |1e; d1++) {. |
00000b50 64 31 2d 3e 72 65 20 3e 3e 3d 20 6c 3b 0a 20 20 |d1->re >>= l;. |
00000b60 20 20 64 31 2d 3e 69 6d 20 3e 3e 3d 20 6c 3b 0a | d1->im >>= l;.|
00000b70 20 20 7d 0a 0a 20 20 72 65 74 75 72 6e 20 54 52 | }.. return TR|
00000b80 55 45 3b 0a 7d 0a 0a 42 4f 4f 4c 20 66 69 6c 74 |UE;.}..BOOL filt|
00000b90 65 72 28 63 6f 6d 70 6c 65 78 20 2a 64 61 74 61 |er(complex *data|
00000ba0 2c 20 69 6e 74 20 6c 69 6d 69 74 2c 20 64 6f 75 |, int limit, dou|
00000bb0 62 6c 65 20 64 69 6d 65 6e 73 69 6f 6e 29 0a 7b |ble dimension).{|
00000bc0 0a 20 20 69 6e 74 20 6e 3d 31 3c 3c 6c 69 6d 69 |. int n=1<<limi|
00000bd0 74 3b 0a 20 20 69 6e 74 20 69 2c 20 6a 2c 20 6b |t;. int i, j, k|
00000be0 32 2c 20 6c 2c 20 6d 3b 0a 20 20 69 6e 74 20 6d |2, l, m;. int m|
00000bf0 62 70 6f 62 74 20 3d 20 28 69 6e 74 29 28 2d 28 |bpobt = (int)(-(|
00000c00 34 2d 64 69 6d 65 6e 73 69 6f 6e 29 2a 6f 6e 65 |4-dimension)*one|
00000c10 29 3b 0a 20 20 69 6e 74 20 73 63 61 6c 65 3b 0a |);. int scale;.|
00000c20 20 20 63 6f 6d 70 6c 65 78 20 2a 64 3b 0a 0a 20 | complex *d;.. |
00000c30 20 66 6f 72 20 28 69 3d 30 2c 20 64 3d 64 61 74 | for (i=0, d=dat|
00000c40 61 3b 20 69 3c 6e 3b 20 69 2b 2b 29 20 7b 0a 20 |a; i<n; i++) {. |
00000c50 20 20 20 6b 32 20 3d 20 69 3c 6e 2f 32 20 3f 20 | k2 = i<n/2 ? |
00000c60 69 20 3a 20 6e 2d 69 3b 0a 20 20 20 20 6b 32 20 |i : n-i;. k2 |
00000c70 3d 20 6b 32 2a 6b 32 3b 0a 20 20 20 20 66 6f 72 |= k2*k2;. for|
00000c80 20 28 6a 3d 30 3b 20 6a 3c 6e 3b 20 6a 2b 2b 2c | (j=0; j<n; j++,|
00000c90 20 64 2b 2b 29 20 7b 0a 20 20 20 20 20 20 69 66 | d++) {. if|
00000ca0 20 28 69 3d 3d 30 20 26 26 20 6a 3d 3d 30 29 20 | (i==0 && j==0) |
00000cb0 73 63 61 6c 65 3d 30 3b 0a 20 20 20 20 20 20 65 |scale=0;. e|
00000cc0 6c 73 65 20 7b 0a 20 20 20 20 20 20 20 20 6c 20 |lse {. l |
00000cd0 3d 20 6a 3c 6e 2f 32 20 3f 20 6a 20 3a 20 6e 2d |= j<n/2 ? j : n-|
00000ce0 6a 3b 0a 09 2f 2a 77 61 6e 74 20 6d 5e 28 6d 62 |j;../*want m^(mb|
00000cf0 70 6f 62 74 2f 32 29 2c 20 62 75 74 20 72 61 6e |pobt/2), but ran|
00000d00 67 65 20 6f 6e 20 6d 20 74 6f 6f 20 67 72 65 61 |ge on m too grea|
00000d10 74 20 66 6f 72 20 70 6f 77 31 36 2c 20 73 6f 20 |t for pow16, so |
00000d20 74 77 65 61 6b 20 6d 20 74 6f 20 61 6c 6c 6f 77 |tweak m to allow|
00000d30 20 61 70 70 6c 69 63 61 74 69 6f 6e 20 6f 66 0a | application of.|
00000d40 09 20 20 73 71 72 74 5f 66 72 61 63 31 36 2c 20 |. sqrt_frac16, |
00000d50 74 68 65 6e 20 72 61 69 73 65 20 74 6f 20 70 6f |then raise to po|
00000d60 77 65 72 20 6d 62 70 6f 62 74 2a 2f 0a 20 20 20 |wer mbpobt*/. |
00000d70 20 20 20 20 20 6d 20 3d 20 6b 32 2b 6c 2a 6c 3b | m = k2+l*l;|
00000d80 0a 20 20 20 20 20 20 20 20 73 63 61 6c 65 20 3d |. scale =|
00000d90 20 70 6f 77 31 36 28 73 71 72 74 5f 66 72 61 63 | pow16(sqrt_frac|
00000da0 31 36 28 6d 2a 28 6f 6e 65 2f 34 29 29 2a 32 2c |16(m*(one/4))*2,|
00000db0 20 6d 62 70 6f 62 74 29 3b 0a 20 20 20 20 20 20 | mbpobt);. |
00000dc0 7d 0a 20 20 20 20 20 20 64 2d 3e 72 65 20 3d 20 |}. d->re = |
00000dd0 6d 75 6c 5f 66 72 61 63 31 36 28 64 2d 3e 72 65 |mul_frac16(d->re|
00000de0 2c 20 73 63 61 6c 65 29 3b 0a 20 20 20 20 20 20 |, scale);. |
00000df0 64 2d 3e 69 6d 20 3d 20 6d 75 6c 5f 66 72 61 63 |d->im = mul_frac|
00000e00 31 36 28 64 2d 3e 69 6d 2c 20 73 63 61 6c 65 29 |16(d->im, scale)|
00000e10 3b 0a 20 20 20 20 7d 0a 20 20 7d 0a 0a 20 20 72 |;. }. }.. r|
00000e20 65 74 75 72 6e 20 54 52 55 45 3b 0a 7d 0a 0a 64 |eturn TRUE;.}..d|
00000e30 6f 75 62 6c 65 20 61 73 73 69 67 6e 28 63 6f 6d |ouble assign(com|
00000e40 70 6c 65 78 20 2a 64 61 74 61 2c 20 69 6e 74 20 |plex *data, int |
00000e50 6c 69 6d 69 74 2c 20 69 6e 74 20 66 75 6e 63 74 |limit, int funct|
00000e60 69 6f 6e 29 0a 7b 0a 20 20 69 6e 74 20 6e 3d 31 |ion).{. int n=1|
00000e70 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 20 69 |<<limit;. int i|
00000e80 2c 20 6a 3b 0a 20 20 69 6e 74 20 76 3b 0a 20 20 |, j;. int v;. |
00000e90 64 6f 75 62 6c 65 20 54 3b 0a 20 20 64 6f 75 62 |double T;. doub|
00000ea0 6c 65 20 64 69 6d 3b 0a 20 20 63 6f 6d 70 6c 65 |le dim;. comple|
00000eb0 78 20 2a 64 31 3b 0a 0a 20 20 73 77 69 74 63 68 |x *d1;.. switch|
00000ec0 20 28 66 75 6e 63 74 69 6f 6e 29 0a 20 20 7b 0a | (function). {.|
00000ed0 20 20 20 20 63 61 73 65 20 30 3a 0a 20 20 20 20 | case 0:. |
00000ee0 20 20 54 20 3d 20 31 3b 0a 20 20 20 20 20 20 66 | T = 1;. f|
00000ef0 6e 5f 64 65 73 63 20 3d 20 22 73 69 6e 65 20 77 |n_desc = "sine w|
00000f00 61 76 65 20 70 61 63 6b 65 74 22 3b 0a 20 20 20 |ave packet";. |
00000f10 20 20 20 66 74 5f 64 65 73 63 20 3d 20 22 22 3b | ft_desc = "";|
00000f20 0a 20 20 20 20 20 20 2a 69 74 5f 64 65 73 63 20 |. *it_desc |
00000f30 3d 20 30 3b 0a 20 20 20 20 20 20 66 6f 72 20 28 |= 0;. for (|
00000f40 69 3d 30 2c 20 64 31 3d 64 61 74 61 3b 20 69 3c |i=0, d1=data; i<|
00000f50 6e 3b 20 69 2b 2b 29 20 7b 0a 20 20 20 20 20 20 |n; i++) {. |
00000f60 20 20 76 20 3d 20 73 69 6e 31 36 28 64 69 76 5f | v = sin16(div_|
00000f70 66 72 61 63 31 36 28 34 2a 69 2c 20 6e 29 29 3b |frac16(4*i, n));|
00000f80 0a 20 20 20 20 20 20 20 20 66 6f 72 20 28 6a 3d |. for (j=|
00000f90 30 3b 20 6a 3c 6e 3b 20 6a 2b 2b 2c 20 64 31 2b |0; j<n; j++, d1+|
00000fa0 2b 29 20 7b 0a 20 20 20 20 20 20 20 20 20 20 64 |+) {. d|
00000fb0 31 2d 3e 72 65 20 3d 20 31 36 2a 6d 75 6c 5f 66 |1->re = 16*mul_f|
00000fc0 72 61 63 31 36 63 28 73 69 6e 31 36 28 64 69 76 |rac16c(sin16(div|
00000fd0 5f 66 72 61 63 31 36 28 34 2a 6a 2c 20 6e 29 29 |_frac16(4*j, n))|
00000fe0 2c 20 76 29 3b 0a 20 20 20 20 20 20 20 20 20 20 |, v);. |
00000ff0 64 31 2d 3e 69 6d 20 3d 20 30 3b 0a 20 20 20 20 |d1->im = 0;. |
00001000 20 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a 20 20 | }. }. |
00001010 20 20 20 20 62 72 65 61 6b 3b 0a 20 20 20 20 63 | break;. c|
00001020 61 73 65 20 31 3a 0a 20 20 20 20 20 20 54 20 3d |ase 1:. T =|
00001030 20 35 3b 0a 20 20 20 20 20 20 66 6e 5f 64 65 73 | 5;. fn_des|
00001040 63 20 3d 20 22 64 61 6d 70 65 64 20 73 69 6e 65 |c = "damped sine|
00001050 20 77 61 76 65 22 3b 0a 20 20 20 20 20 20 66 74 | wave";. ft|
00001060 5f 64 65 73 63 20 3d 20 22 22 3b 0a 20 20 20 20 |_desc = "";. |
00001070 20 20 2a 69 74 5f 64 65 73 63 20 3d 20 30 3b 0a | *it_desc = 0;.|
00001080 20 20 20 20 20 20 66 6f 72 20 28 69 3d 30 2c 20 | for (i=0, |
00001090 64 31 3d 64 61 74 61 3b 20 69 3c 6e 3b 20 69 2b |d1=data; i<n; i+|
000010a0 2b 29 20 7b 0a 20 20 20 20 20 20 20 20 76 20 3d |+) {. v =|
000010b0 20 69 3d 3d 30 20 3f 20 6f 6e 65 20 3a 20 6d 75 | i==0 ? one : mu|
000010c0 6c 5f 66 72 61 63 31 36 63 28 64 69 76 5f 66 72 |l_frac16c(div_fr|
000010d0 61 63 31 36 28 6e 2c 20 69 29 2c 20 73 69 6e 31 |ac16(n, i), sin1|
000010e0 36 28 64 69 76 5f 66 72 61 63 31 36 28 34 2a 69 |6(div_frac16(4*i|
000010f0 2c 20 6e 29 29 29 3b 0a 20 20 20 20 20 20 20 20 |, n)));. |
00001100 66 6f 72 20 28 6a 3d 30 3b 20 6a 3c 6e 3b 20 6a |for (j=0; j<n; j|
00001110 2b 2b 2c 20 64 31 2b 2b 29 20 7b 0a 20 20 20 20 |++, d1++) {. |
00001120 20 20 20 20 20 20 64 31 2d 3e 72 65 20 3d 20 6d | d1->re = m|
00001130 75 6c 5f 66 72 61 63 31 36 63 28 6a 3d 3d 30 20 |ul_frac16c(j==0 |
00001140 3f 20 6f 6e 65 20 3a 20 6d 75 6c 5f 66 72 61 63 |? one : mul_frac|
00001150 31 36 63 28 64 69 76 5f 66 72 61 63 31 36 28 6e |16c(div_frac16(n|
00001160 2c 20 6a 29 2c 20 73 69 6e 31 36 28 64 69 76 5f |, j), sin16(div_|
00001170 66 72 61 63 31 36 28 34 2a 6a 2c 20 6e 29 29 29 |frac16(4*j, n)))|
00001180 2c 20 76 29 3b 0a 20 20 20 20 20 20 20 20 20 20 |, v);. |
00001190 64 31 2d 3e 69 6d 20 3d 20 30 3b 0a 20 20 20 20 |d1->im = 0;. |
000011a0 20 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a 20 20 | }. }. |
000011b0 20 20 20 20 62 72 65 61 6b 3b 0a 20 20 20 20 63 | break;. c|
000011c0 61 73 65 20 32 3a 0a 20 20 20 20 20 20 54 20 3d |ase 2:. T =|
000011d0 20 31 3b 0a 20 20 20 20 20 20 66 6e 5f 64 65 73 | 1;. fn_des|
000011e0 63 20 3d 20 22 73 71 75 61 72 65 20 77 61 76 65 |c = "square wave|
000011f0 20 70 61 63 6b 65 74 22 3b 0a 20 20 20 20 20 20 | packet";. |
00001200 66 74 5f 64 65 73 63 20 3d 20 22 22 3b 0a 20 20 |ft_desc = "";. |
00001210 20 20 20 20 2a 69 74 5f 64 65 73 63 20 3d 20 30 | *it_desc = 0|
00001220 3b 0a 20 20 20 20 20 20 66 6f 72 20 28 69 3d 30 |;. for (i=0|
00001230 2c 20 64 31 3d 64 61 74 61 3b 20 69 3c 6e 3b 20 |, d1=data; i<n; |
00001240 69 2b 2b 29 20 7b 0a 20 20 20 20 20 20 20 20 66 |i++) {. f|
00001250 6f 72 20 28 6a 3d 30 3b 20 6a 3c 6e 3b 20 6a 2b |or (j=0; j<n; j+|
00001260 2b 2c 20 64 31 2b 2b 29 20 7b 0a 20 20 20 20 20 |+, d1++) {. |
00001270 20 20 20 20 20 64 31 2d 3e 72 65 20 3d 20 69 3c | d1->re = i<|
00001280 6e 2f 32 20 26 26 20 6a 3c 6e 2f 32 20 3f 20 31 |n/2 && j<n/2 ? 1|
00001290 36 2a 6f 6e 65 20 3a 20 30 3b 0a 20 20 20 20 20 |6*one : 0;. |
000012a0 20 20 20 20 20 64 31 2d 3e 69 6d 20 3d 20 30 3b | d1->im = 0;|
000012b0 0a 20 20 20 20 20 20 20 20 7d 0a 20 20 20 20 20 |. }. |
000012c0 20 7d 0a 20 20 20 20 20 20 62 72 65 61 6b 3b 0a | }. break;.|
000012d0 20 20 20 20 63 61 73 65 20 33 3a 0a 20 20 20 20 | case 3:. |
000012e0 20 20 54 20 3d 20 31 3b 0a 20 20 20 20 20 20 66 | T = 1;. f|
000012f0 6e 5f 64 65 73 63 20 3d 20 22 77 68 69 74 65 20 |n_desc = "white |
00001300 6e 6f 69 73 65 22 3b 0a 20 20 20 20 20 20 66 74 |noise";. ft|
00001310 5f 64 65 73 63 20 3d 20 22 22 3b 0a 20 20 20 20 |_desc = "";. |
00001320 20 20 2a 69 74 5f 64 65 73 63 20 3d 20 30 3b 0a | *it_desc = 0;.|
00001330 20 20 20 20 20 20 66 6f 72 20 28 69 3d 30 2c 20 | for (i=0, |
00001340 64 31 3d 64 61 74 61 3b 20 69 3c 6e 3b 20 69 2b |d1=data; i<n; i+|
00001350 2b 29 20 7b 0a 20 20 20 20 20 20 20 20 66 6f 72 |+) {. for|
00001360 20 28 6a 3d 30 3b 20 6a 3c 6e 3b 20 6a 2b 2b 2c | (j=0; j<n; j++,|
00001370 20 64 31 2b 2b 29 20 7b 0a 20 20 20 20 20 20 20 | d1++) {. |
00001380 20 20 20 64 31 2d 3e 72 65 20 3d 20 31 36 2a 28 | d1->re = 16*(|
00001390 28 72 61 6e 64 28 29 3e 3e 31 35 29 2d 6f 6e 65 |(rand()>>15)-one|
000013a0 2f 32 29 3b 0a 20 20 20 20 20 20 20 20 20 20 64 |/2);. d|
000013b0 31 2d 3e 69 6d 20 3d 20 30 3b 0a 20 20 20 20 20 |1->im = 0;. |
000013c0 20 20 20 7d 0a 20 20 20 20 20 20 7d 0a 20 20 20 | }. }. |
000013d0 20 20 20 62 72 65 61 6b 3b 0a 20 20 20 20 63 61 | break;. ca|
000013e0 73 65 20 34 3a 0a 20 20 20 20 20 20 54 20 3d 20 |se 4:. T = |
000013f0 31 3b 0a 20 20 20 20 20 20 66 6e 5f 64 65 73 63 |1;. fn_desc|
00001400 20 3d 20 22 67 61 75 73 73 69 61 6e 20 77 68 69 | = "gaussian whi|
00001410 74 65 20 6e 6f 69 73 65 22 3b 0a 20 20 20 20 20 |te noise";. |
00001420 20 66 74 5f 64 65 73 63 20 3d 20 22 22 3b 0a 20 | ft_desc = "";. |
00001430 20 20 20 20 20 2a 69 74 5f 64 65 73 63 20 3d 20 | *it_desc = |
00001440 30 3b 0a 20 20 20 20 20 20 66 6f 72 20 28 69 3d |0;. for (i=|
00001450 30 2c 20 64 31 3d 64 61 74 61 3b 20 69 3c 6e 3b |0, d1=data; i<n;|
00001460 20 69 2b 2b 29 20 7b 0a 20 20 20 20 20 20 20 20 | i++) {. |
00001470 66 6f 72 20 28 6a 3d 30 3b 20 6a 3c 6e 3b 20 6a |for (j=0; j<n; j|
00001480 2b 2b 2c 20 64 31 2b 2b 29 20 7b 0a 20 20 20 20 |++, d1++) {. |
00001490 20 20 20 20 20 20 64 31 2d 3e 72 65 20 3d 20 31 | d1->re = 1|
000014a0 36 2a 67 61 75 73 73 31 36 28 29 3b 0a 20 20 20 |6*gauss16();. |
000014b0 20 20 20 20 20 20 20 64 31 2d 3e 69 6d 20 3d 20 | d1->im = |
000014c0 30 3b 0a 20 20 20 20 20 20 20 20 7d 0a 20 20 20 |0;. }. |
000014d0 20 20 20 7d 0a 20 20 20 20 20 20 62 72 65 61 6b | }. break|
000014e0 3b 0a 20 20 20 20 64 65 66 61 75 6c 74 3a 20 2f |;. default: /|
000014f0 2a 66 6f 72 20 36 20 63 61 73 65 73 2c 20 35 2d |*for 6 cases, 5-|
00001500 31 30 2a 2f 0a 20 20 20 20 20 20 54 20 3d 20 31 |10*/. T = 1|
00001510 3b 0a 20 20 20 20 20 20 66 6e 5f 64 65 73 63 20 |;. fn_desc |
00001520 3d 20 22 67 61 75 73 73 69 61 6e 20 77 68 69 74 |= "gaussian whit|
00001530 65 20 6e 6f 69 73 65 22 3b 0a 20 20 20 20 20 20 |e noise";. |
00001540 66 74 5f 64 65 73 63 20 3d 20 22 46 69 6c 74 65 |ft_desc = "Filte|
00001550 72 65 64 20 62 79 20 31 2f 28 66 5e 6b 29 22 3b |red by 1/(f^k)";|
00001560 0a 20 20 20 20 20 20 64 69 6d 20 3d 20 64 69 6d |. dim = dim|
00001570 73 6e 28 66 75 6e 63 74 69 6f 6e 29 3b 0a 20 20 |sn(function);. |
00001580 20 20 20 20 73 70 72 69 6e 74 66 28 69 74 5f 64 | sprintf(it_d|
00001590 65 73 63 2c 20 22 6b 20 3d 20 25 2e 33 66 2c 20 |esc, "k = %.3f, |
000015a0 46 72 61 63 74 61 6c 20 64 69 6d 2e 20 3d 20 25 |Fractal dim. = %|
000015b0 2e 33 66 22 2c 20 34 2e 30 2d 64 69 6d 2c 20 64 |.3f", 4.0-dim, d|
000015c0 69 6d 29 3b 0a 20 20 20 20 20 20 66 6f 72 20 28 |im);. for (|
000015d0 69 3d 30 2c 20 64 31 3d 64 61 74 61 3b 20 69 3c |i=0, d1=data; i<|
000015e0 6e 3b 20 69 2b 2b 29 20 7b 0a 20 20 20 20 20 20 |n; i++) {. |
000015f0 20 20 66 6f 72 20 28 6a 3d 30 3b 20 6a 3c 6e 3b | for (j=0; j<n;|
00001600 20 6a 2b 2b 2c 20 64 31 2b 2b 29 20 7b 0a 20 20 | j++, d1++) {. |
00001610 20 20 20 20 20 20 20 20 64 31 2d 3e 72 65 20 3d | d1->re =|
00001620 20 31 36 2a 67 61 75 73 73 31 36 28 29 3b 0a 20 | 16*gauss16();. |
00001630 20 20 20 20 20 20 20 20 20 64 31 2d 3e 69 6d 20 | d1->im |
00001640 3d 20 30 3b 0a 20 20 20 20 20 20 20 20 7d 0a 20 |= 0;. }. |
00001650 20 20 20 20 20 7d 0a 20 20 20 20 20 20 62 72 65 | }. bre|
00001660 61 6b 3b 0a 20 20 7d 0a 0a 20 20 72 65 74 75 72 |ak;. }.. retur|
00001670 6e 20 54 3b 0a 7d 0a 0a 23 64 65 66 69 6e 65 20 |n T;.}..#define |
00001680 72 65 6f 72 64 65 72 28 69 2c 20 6e 2c 20 66 6c |reorder(i, n, fl|
00001690 61 67 29 20 28 28 66 6c 61 67 29 20 3f 20 28 28 |ag) ((flag) ? ((|
000016a0 69 29 3c 28 6e 29 2f 32 20 3f 20 28 69 29 2b 28 |i)<(n)/2 ? (i)+(|
000016b0 6e 29 2f 32 20 3a 20 28 69 29 2d 28 6e 29 2f 32 |n)/2 : (i)-(n)/2|
000016c0 29 20 3a 20 28 69 29 29 0a 23 64 65 66 69 6e 65 |) : (i)).#define|
000016d0 20 79 73 63 61 6c 65 20 31 32 34 0a 0a 42 4f 4f | yscale 124..BOO|
000016e0 4c 20 67 72 61 70 68 28 63 6f 6d 70 6c 65 78 20 |L graph(complex |
000016f0 2a 64 61 74 61 2c 20 69 6e 74 20 6c 69 6d 69 74 |*data, int limit|
00001700 2c 20 69 6e 74 20 79 6f 2c 20 42 4f 4f 4c 20 72 |, int yo, BOOL r|
00001710 65 70 6f 73 2c 20 63 68 61 72 20 2a 74 65 78 74 |epos, char *text|
00001720 2c 20 2e 2e 2e 29 0a 7b 0a 20 20 69 6e 74 20 6e |, ...).{. int n|
00001730 3d 31 3c 3c 6c 69 6d 69 74 3b 0a 20 20 69 6e 74 |=1<<limit;. int|
00001740 20 69 2c 20 6a 3b 0a 20 20 69 6e 74 20 6d 69 6e | i, j;. int min|
00001750 3d 31 3c 3c 33 30 3b 0a 20 20 69 6e 74 20 6d 61 |=1<<30;. int ma|
00001760 78 3d 2d 6d 69 6e 3b 0a 20 20 69 6e 74 20 73 63 |x=-min;. int sc|
00001770 61 6c 65 3b 0a 20 20 63 6f 6d 70 6c 65 78 20 2a |ale;. complex *|
00001780 64 31 2c 20 2a 64 31 65 3b 0a 20 20 76 61 5f 6c |d1, *d1e;. va_l|
00001790 69 73 74 20 61 72 67 70 3b 0a 0a 20 20 66 6f 72 |ist argp;.. for|
000017a0 20 28 64 31 3d 64 61 74 61 2c 20 64 31 65 3d 64 | (d1=data, d1e=d|
000017b0 31 2b 6e 2a 6e 3b 20 64 31 3c 64 31 65 3b 20 64 |1+n*n; d1<d1e; d|
000017c0 31 2b 2b 29 20 7b 0a 20 20 20 20 69 66 20 28 64 |1++) {. if (d|
000017d0 31 2d 3e 72 65 3c 6d 69 6e 29 20 6d 69 6e 3d 64 |1->re<min) min=d|
000017e0 31 2d 3e 72 65 3b 0a 20 20 20 20 69 66 20 28 64 |1->re;. if (d|
000017f0 31 2d 3e 69 6d 3c 6d 69 6e 29 20 6d 69 6e 3d 64 |1->im<min) min=d|
00001800 31 2d 3e 69 6d 3b 0a 20 20 20 20 69 66 20 28 64 |1->im;. if (d|
00001810 31 2d 3e 72 65 3e 6d 61 78 29 20 6d 61 78 3d 64 |1->re>max) max=d|
00001820 31 2d 3e 72 65 3b 0a 20 20 20 20 69 66 20 28 64 |1->re;. if (d|
00001830 31 2d 3e 69 6d 3e 6d 61 78 29 20 6d 61 78 3d 64 |1->im>max) max=d|
00001840 31 2d 3e 69 6d 3b 0a 20 20 7d 0a 20 20 6d 61 78 |1->im;. }. max|
00001850 20 3d 20 6d 61 78 3e 30 20 3f 20 6d 61 78 20 3a | = max>0 ? max :|
00001860 20 2d 6d 61 78 3b 0a 20 20 6d 69 6e 20 3d 20 6d | -max;. min = m|
00001870 69 6e 3e 30 20 3f 20 6d 69 6e 20 3a 20 2d 6d 69 |in>0 ? min : -mi|
00001880 6e 3b 0a 20 20 73 63 61 6c 65 20 3d 20 6d 61 78 |n;. scale = max|
00001890 3e 6d 69 6e 20 3f 20 6d 61 78 20 3a 20 6d 69 6e |>min ? max : min|
000018a0 3b 0a 20 20 69 66 20 28 73 63 61 6c 65 3d 3d 30 |;. if (scale==0|
000018b0 29 20 73 63 61 6c 65 3d 31 3b 0a 0a 20 20 62 62 |) scale=1;.. bb|
000018c0 63 5f 6f 72 69 67 69 6e 28 36 34 30 2c 79 6f 29 |c_origin(640,yo)|
000018d0 3b 0a 20 20 62 62 63 5f 67 63 6f 6c 28 30 2c 20 |;. bbc_gcol(0, |
000018e0 31 29 3b 0a 20 20 62 62 63 5f 72 65 63 74 61 6e |1);. bbc_rectan|
000018f0 67 6c 65 66 69 6c 6c 28 2d 35 31 37 2c 2d 34 2c |glefill(-517,-4,|
00001900 31 30 33 31 2c 39 29 3b 0a 20 20 62 62 63 5f 72 |1031,9);. bbc_r|
00001910 65 63 74 61 6e 67 6c 65 66 69 6c 6c 28 72 65 70 |ectanglefill(rep|
00001920 6f 73 3f 2d 35 3a 2d 35 31 37 2c 2d 79 73 63 61 |os?-5:-517,-ysca|
00001930 6c 65 2d 34 2c 39 2c 32 2a 79 73 63 61 6c 65 2b |le-4,9,2*yscale+|
00001940 37 29 3b 0a 20 20 62 62 63 5f 76 64 75 28 35 29 |7);. bbc_vdu(5)|
00001950 3b 0a 20 20 62 62 63 5f 6d 6f 76 65 28 28 72 65 |;. bbc_move((re|
00001960 70 6f 73 3f 30 3a 2d 35 31 32 29 2d 34 34 2c 79 |pos?0:-512)-44,y|
00001970 73 63 61 6c 65 2b 34 30 29 3b 0a 20 20 70 72 69 |scale+40);. pri|
00001980 6e 74 66 28 22 25 2e 33 66 22 2c 20 28 64 6f 75 |ntf("%.3f", (dou|
00001990 62 6c 65 29 73 63 61 6c 65 2f 6f 6e 65 29 3b 0a |ble)scale/one);.|
000019a0 20 20 62 62 63 5f 6d 6f 76 65 28 2d 35 31 32 2b | bbc_move(-512+|
000019b0 39 36 2c 79 73 63 61 6c 65 2b 34 30 29 3b 0a 20 |96,yscale+40);. |
000019c0 20 76 61 5f 73 74 61 72 74 28 61 72 67 70 2c 20 | va_start(argp, |
000019d0 74 65 78 74 29 3b 0a 20 20 76 70 72 69 6e 74 66 |text);. vprintf|
000019e0 28 74 65 78 74 2c 20 61 72 67 70 29 3b 0a 20 20 |(text, argp);. |
000019f0 76 61 5f 65 6e 64 28 61 72 67 70 29 3b 0a 20 20 |va_end(argp);. |
00001a00 62 62 63 5f 6d 6f 76 65 28 28 72 65 70 6f 73 3f |bbc_move((repos?|
00001a10 30 3a 2d 35 31 32 29 2d 36 30 2c 2d 79 73 63 61 |0:-512)-60,-ysca|
00001a20 6c 65 2d 32 30 29 3b 0a 20 20 70 72 69 6e 74 66 |le-20);. printf|
00001a30 28 22 25 2e 33 66 22 2c 20 2d 28 64 6f 75 62 6c |("%.3f", -(doubl|
00001a40 65 29 73 63 61 6c 65 2f 6f 6e 65 29 3b 0a 20 20 |e)scale/one);. |
00001a50 62 62 63 5f 6d 6f 76 65 28 2d 35 31 32 2d 31 32 |bbc_move(-512-12|
00001a60 38 2c 6c 69 6e 65 73 65 70 2f 34 29 3b 0a 20 20 |8,linesep/4);. |
00001a70 70 72 69 6e 74 66 28 22 25 2e 33 66 22 2c 20 72 |printf("%.3f", r|
00001a80 65 70 6f 73 3f 2d 46 2f 32 3a 30 2e 30 29 3b 0a |epos?-F/2:0.0);.|
00001a90 20 20 62 62 63 5f 6d 6f 76 65 28 35 32 30 2c 6c | bbc_move(520,l|
00001aa0 69 6e 65 73 65 70 2f 34 29 3b 0a 20 20 70 72 69 |inesep/4);. pri|
00001ab0 6e 74 66 28 22 25 2e 33 66 22 2c 20 72 65 70 6f |ntf("%.3f", repo|
00001ac0 73 3f 46 2f 32 3a 54 29 3b 0a 20 20 62 62 63 5f |s?F/2:T);. bbc_|
00001ad0 76 64 75 28 34 29 3b 0a 20 20 66 6f 72 20 28 69 |vdu(4);. for (i|
00001ae0 3d 30 3b 20 69 3c 6e 3b 20 69 2b 2b 29 20 7b 0a |=0; i<n; i++) {.|
00001af0 20 20 20 20 62 62 63 5f 67 63 6f 6c 28 30 2c 20 | bbc_gcol(0, |
00001b00 34 2b 31 32 2a 69 2f 6e 29 3b 0a 20 20 20 20 62 |4+12*i/n);. b|
00001b10 62 63 5f 6d 6f 76 65 28 2d 35 31 32 2c 64 61 74 |bc_move(-512,dat|
00001b20 61 5b 72 65 6f 72 64 65 72 28 30 2c 20 6e 2c 20 |a[reorder(0, n, |
00001b30 72 65 70 6f 73 29 2b 6e 2a 72 65 6f 72 64 65 72 |repos)+n*reorder|
00001b40 28 69 2c 20 6e 2c 20 72 65 70 6f 73 29 5d 2e 72 |(i, n, repos)].r|
00001b50 65 2a 79 73 63 61 6c 65 2f 73 63 61 6c 65 29 3b |e*yscale/scale);|
00001b60 0a 20 20 20 20 66 6f 72 20 28 6a 3d 31 3b 20 6a |. for (j=1; j|
00001b70 3c 6e 3b 20 6a 2b 2b 29 20 7b 0a 20 20 20 20 20 |<n; j++) {. |
00001b80 20 62 62 63 5f 70 6c 6f 74 28 31 33 2c 20 6a 2a | bbc_plot(13, j*|
00001b90 31 30 32 34 2f 6e 2d 35 31 32 2c 64 61 74 61 5b |1024/n-512,data[|
00001ba0 72 65 6f 72 64 65 72 28 6a 2c 20 6e 2c 20 72 65 |reorder(j, n, re|
00001bb0 70 6f 73 29 2b 6e 2a 72 65 6f 72 64 65 72 28 69 |pos)+n*reorder(i|
00001bc0 2c 20 6e 2c 20 72 65 70 6f 73 29 5d 2e 72 65 2a |, n, repos)].re*|
00001bd0 79 73 63 61 6c 65 2f 73 63 61 6c 65 29 3b 0a 20 |yscale/scale);. |
00001be0 20 20 20 7d 0a 20 20 20 20 62 62 63 5f 67 63 6f | }. bbc_gco|
00001bf0 6c 28 30 2c 20 33 29 3b 0a 20 20 20 20 62 62 63 |l(0, 3);. bbc|
00001c00 5f 6d 6f 76 65 28 2d 35 31 32 2c 64 61 74 61 5b |_move(-512,data[|
00001c10 72 65 6f 72 64 65 72 28 30 2c 20 6e 2c 20 72 65 |reorder(0, n, re|
00001c20 70 6f 73 29 2b 6e 2a 72 65 6f 72 64 65 72 28 69 |pos)+n*reorder(i|
00001c30 2c 20 6e 2c 20 72 65 70 6f 73 29 5d 2e 69 6d 2a |, n, repos)].im*|
00001c40 79 73 63 61 6c 65 2f 73 63 61 6c 65 29 3b 0a 20 |yscale/scale);. |
00001c50 20 20 20 66 6f 72 20 28 6a 3d 31 3b 20 6a 3c 6e | for (j=1; j<n|
00001c60 3b 20 6a 2b 2b 29 20 7b 0a 20 20 20 20 20 20 62 |; j++) {. b|
00001c70 62 63 5f 70 6c 6f 74 28 31 33 2c 20 6a 2a 31 30 |bc_plot(13, j*10|
00001c80 32 34 2f 6e 2d 35 31 32 2c 64 61 74 61 5b 72 65 |24/n-512,data[re|
00001c90 6f 72 64 65 72 28 6a 2c 20 6e 2c 20 72 65 70 6f |order(j, n, repo|
00001ca0 73 29 2b 6e 2a 72 65 6f 72 64 65 72 28 69 2c 20 |s)+n*reorder(i, |
00001cb0 6e 2c 20 72 65 70 6f 73 29 5d 2e 69 6d 2a 79 73 |n, repos)].im*ys|
00001cc0 63 61 6c 65 2f 73 63 61 6c 65 29 3b 0a 20 20 20 |cale/scale);. |
00001cd0 20 7d 0a 20 20 7d 0a 20 20 72 65 74 75 72 6e 20 | }. }. return |
00001ce0 54 52 55 45 3b 0a 7d 0a 0a 6d 6f 6e 69 74 6f 72 |TRUE;.}..monitor|
00001cf0 20 72 65 61 64 5f 6d 6f 6e 69 74 6f 72 5f 74 79 | read_monitor_ty|
00001d00 70 65 28 76 6f 69 64 29 0a 7b 0a 20 20 69 6e 74 |pe(void).{. int|
00001d10 20 72 32 3b 0a 20 20 6f 73 5f 73 77 69 33 72 28 | r2;. os_swi3r(|
00001d20 36 2c 20 31 36 31 2c 31 33 33 2c 30 2c 20 30 2c |6, 161,133,0, 0,|
00001d30 30 2c 26 72 32 29 3b 0a 20 20 69 66 20 28 28 72 |0,&r2);. if ((r|
00001d40 32 2f 34 20 26 20 33 29 20 3d 3d 20 31 29 20 72 |2/4 & 3) == 1) r|
00001d50 65 74 75 72 6e 20 6d 75 6c 74 69 3b 0a 20 20 65 |eturn multi;. e|
00001d60 6c 73 65 20 72 65 74 75 72 6e 20 6e 6f 6e 6d 75 |lse return nonmu|
00001d70 6c 74 69 3b 0a 7d 0a 0a 69 6e 74 20 6d 61 69 6e |lti;.}..int main|
00001d80 28 76 6f 69 64 29 0a 7b 0a 20 20 69 6e 74 20 66 |(void).{. int f|
00001d90 75 6e 63 74 69 6f 6e 2c 20 73 65 65 64 2c 20 69 |unction, seed, i|
00001da0 3b 0a 20 20 63 68 61 72 20 64 75 6d 6d 79 3b 0a |;. char dummy;.|
00001db0 20 20 6d 6f 6e 69 74 6f 72 74 79 70 65 20 3d 20 | monitortype = |
00001dc0 72 65 61 64 5f 6d 6f 6e 69 74 6f 72 5f 74 79 70 |read_monitor_typ|
00001dd0 65 28 29 3b 0a 0a 20 20 66 6f 72 20 28 3b 3b 29 |e();.. for (;;)|
00001de0 20 7b 0a 0a 20 20 20 20 73 65 65 64 20 3d 20 2d | {.. seed = -|
00001df0 63 6c 6f 63 6b 28 29 3b 0a 20 20 20 20 64 6f 20 |clock();. do |
00001e00 7b 0a 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 |{. printf("|
00001e10 50 6c 65 61 73 65 20 65 6e 74 65 72 20 70 6f 77 |Please enter pow|
00001e20 65 72 20 28 73 75 67 67 65 73 74 20 35 2d 38 2c |er (suggest 5-8,|
00001e30 20 74 68 6f 75 67 68 20 33 2d 39 20 70 65 72 6d | though 3-9 perm|
00001e40 69 73 73 69 62 6c 65 29 20 22 29 3b 0a 20 20 20 |issible) ");. |
00001e50 20 20 20 73 63 61 6e 66 28 22 25 69 22 2c 20 26 | scanf("%i", &|
00001e60 6c 69 6d 69 74 29 3b 0a 20 20 20 20 7d 20 77 68 |limit);. } wh|
00001e70 69 6c 65 20 28 6c 69 6d 69 74 3c 33 20 7c 7c 20 |ile (limit<3 || |
00001e80 6c 69 6d 69 74 3e 39 29 3b 0a 20 20 20 20 4e 20 |limit>9);. N |
00001e90 3d 20 31 3c 3c 6c 69 6d 69 74 3b 0a 0a 20 20 20 |= 1<<limit;.. |
00001ea0 20 64 61 74 61 20 3d 20 63 61 6c 6c 6f 63 28 4e | data = calloc(N|
00001eb0 2a 4e 2c 20 73 69 7a 65 6f 66 28 63 6f 6d 70 6c |*N, sizeof(compl|
00001ec0 65 78 29 29 3b 0a 20 20 20 20 69 66 20 28 64 61 |ex));. if (da|
00001ed0 74 61 3d 3d 30 29 20 77 65 72 72 28 31 2c 20 22 |ta==0) werr(1, "|
00001ee0 49 6e 73 75 66 66 69 63 69 65 6e 74 20 6d 65 6d |Insufficient mem|
00001ef0 6f 72 79 22 29 3b 0a 0a 20 20 20 20 66 6f 72 20 |ory");.. for |
00001f00 28 66 75 6e 63 74 69 6f 6e 3d 30 3b 20 66 75 6e |(function=0; fun|
00001f10 63 74 69 6f 6e 3c 31 31 3b 20 66 75 6e 63 74 69 |ction<11; functi|
00001f20 6f 6e 2b 2b 29 20 7b 0a 20 20 20 20 20 20 73 67 |on++) {. sg|
00001f30 61 75 73 73 31 36 28 73 65 65 64 29 3b 0a 20 20 |auss16(seed);. |
00001f40 20 20 20 20 73 72 61 6e 64 28 73 65 65 64 29 3b | srand(seed);|
00001f50 0a 20 20 20 20 20 20 62 62 63 5f 6d 6f 64 65 28 |. bbc_mode(|
00001f60 6d 6f 6e 69 74 6f 72 74 79 70 65 3d 3d 6d 75 6c |monitortype==mul|
00001f70 74 69 20 3f 20 32 30 20 3a 20 31 32 29 3b 0a 20 |ti ? 20 : 12);. |
00001f80 20 20 20 20 20 62 62 63 5f 70 61 6c 65 74 74 65 | bbc_palette|
00001f90 28 31 2c 20 31 36 2c 20 35 31 2c 35 31 2c 35 31 |(1, 16, 51,51,51|
00001fa0 29 3b 0a 20 20 20 20 20 20 62 62 63 5f 70 61 6c |);. bbc_pal|
00001fb0 65 74 74 65 28 32 2c 20 31 36 2c 20 32 35 35 2c |ette(2, 16, 255,|
00001fc0 32 35 35 2c 32 35 35 29 3b 0a 20 20 20 20 20 20 |255,255);. |
00001fd0 66 6f 72 20 28 69 3d 30 3b 20 69 3c 31 32 3b 20 |for (i=0; i<12; |
00001fe0 69 2b 2b 29 20 62 62 63 5f 70 61 6c 65 74 74 65 |i++) bbc_palette|
00001ff0 28 34 2b 69 2c 20 31 36 2c 20 30 2c 36 38 2b 31 |(4+i, 16, 0,68+1|
00002000 37 2a 69 2c 30 29 3b 0a 20 20 20 20 20 20 6c 69 |7*i,0);. li|
00002010 6e 65 73 65 70 20 3d 20 6d 6f 6e 69 74 6f 72 74 |nesep = monitort|
00002020 79 70 65 3d 3d 6d 75 6c 74 69 20 3f 20 32 30 20 |ype==multi ? 20 |
00002030 3a 20 34 30 3b 0a 20 20 20 20 20 20 54 3d 61 73 |: 40;. T=as|
00002040 73 69 67 6e 28 64 61 74 61 2c 20 6c 69 6d 69 74 |sign(data, limit|
00002050 2c 20 66 75 6e 63 74 69 6f 6e 29 3b 0a 20 20 20 |, function);. |
00002060 20 20 20 46 3d 4e 2f 54 3b 0a 20 20 20 20 20 20 | F=N/T;. |
00002070 67 72 61 70 68 28 64 61 74 61 2c 20 6c 69 6d 69 |graph(data, limi|
00002080 74 2c 20 38 35 32 2c 20 46 41 4c 53 45 2c 20 22 |t, 852, FALSE, "|
00002090 53 6f 75 72 63 65 20 66 75 6e 63 74 69 6f 6e 20 |Source function |
000020a0 27 25 73 27 20 20 28 53 61 6d 70 6c 65 20 73 69 |'%s' (Sample si|
000020b0 7a 65 20 4e 3d 25 69 29 22 2c 20 66 6e 5f 64 65 |ze N=%i)", fn_de|
000020c0 73 63 2c 20 4e 29 3b 0a 20 20 20 20 20 20 66 66 |sc, N);. ff|
000020d0 74 5f 32 64 28 64 61 74 61 2c 20 6c 69 6d 69 74 |t_2d(data, limit|
000020e0 2c 20 31 29 3b 0a 20 20 20 20 20 20 69 66 20 28 |, 1);. if (|
000020f0 66 75 6e 63 74 69 6f 6e 3e 34 29 20 66 69 6c 74 |function>4) filt|
00002100 65 72 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 |er(data, limit, |
00002110 64 69 6d 73 6e 28 66 75 6e 63 74 69 6f 6e 29 29 |dimsn(function))|
00002120 3b 0a 20 20 20 20 20 20 67 72 61 70 68 28 64 61 |;. graph(da|
00002130 74 61 2c 20 6c 69 6d 69 74 2c 20 35 31 32 2c 20 |ta, limit, 512, |
00002140 54 52 55 45 2c 20 22 46 6f 75 72 69 65 72 20 54 |TRUE, "Fourier T|
00002150 72 61 6e 73 66 6f 72 6d 20 20 20 20 20 20 20 20 |ransform |
00002160 20 20 20 20 20 20 20 20 20 25 73 22 2c 20 66 74 | %s", ft|
00002170 5f 64 65 73 63 29 3b 0a 20 20 20 20 20 20 66 66 |_desc);. ff|
00002180 74 5f 32 64 28 64 61 74 61 2c 20 6c 69 6d 69 74 |t_2d(data, limit|
00002190 2c 20 2d 31 29 3b 0a 20 20 20 20 20 20 67 72 61 |, -1);. gra|
000021a0 70 68 28 64 61 74 61 2c 20 6c 69 6d 69 74 2c 20 |ph(data, limit, |
000021b0 31 37 32 2c 20 46 41 4c 53 45 2c 20 22 49 6e 76 |172, FALSE, "Inv|
000021c0 65 72 74 65 64 20 54 72 61 6e 73 66 6f 72 6d 20 |erted Transform |
000021d0 20 20 20 20 20 20 20 20 20 20 20 20 20 25 73 22 | %s"|
000021e0 2c 20 69 74 5f 64 65 73 63 29 3b 0a 20 20 20 20 |, it_desc);. |
000021f0 20 20 62 62 63 5f 76 64 75 28 35 29 3b 0a 20 20 | bbc_vdu(5);. |
00002200 20 20 20 20 62 62 63 5f 6f 72 69 67 69 6e 28 30 | bbc_origin(0|
00002210 2c 30 29 3b 0a 20 20 20 20 20 20 62 62 63 5f 6d |,0);. bbc_m|
00002220 6f 76 65 28 31 31 35 32 2c 20 31 33 31 29 3b 0a |ove(1152, 131);.|
00002230 20 20 20 20 20 20 70 72 69 6e 74 66 28 22 50 72 | printf("Pr|
00002240 65 73 73 20 61 22 29 3b 0a 20 20 20 20 20 20 62 |ess a");. b|
00002250 62 63 5f 6d 6f 76 65 28 31 31 35 32 2c 20 31 33 |bc_move(1152, 13|
00002260 31 2d 6c 69 6e 65 73 65 70 29 3b 0a 20 20 20 20 |1-linesep);. |
00002270 20 20 70 72 69 6e 74 66 28 22 6b 65 79 20 2e 2e | printf("key ..|
00002280 2e 22 29 3b 0a 20 20 20 20 20 20 62 62 63 5f 76 |.");. bbc_v|
00002290 64 75 28 34 29 3b 0a 20 20 20 20 20 20 6f 73 5f |du(4);. os_|
000022a0 63 6c 69 28 22 66 78 20 31 35 20 31 22 29 3b 0a |cli("fx 15 1");.|
000022b0 20 20 20 20 20 20 62 62 63 5f 63 75 72 73 6f 72 | bbc_cursor|
000022c0 28 30 29 3b 0a 20 20 20 20 20 20 64 75 6d 6d 79 |(0);. dummy|
000022d0 20 3d 20 62 62 63 5f 67 65 74 28 29 3b 0a 20 20 | = bbc_get();. |
000022e0 20 20 7d 0a 0a 20 20 20 20 66 72 65 65 28 64 61 | }.. free(da|
000022f0 74 61 29 3b 0a 20 20 20 20 62 62 63 5f 63 6c 73 |ta);. bbc_cls|
00002300 28 29 3b 0a 0a 20 7d 0a 0a 20 20 72 65 74 75 72 |();.. }.. retur|
00002310 6e 20 30 3b 0a 7d 0a |n 0;.}.|
00002317 
.