Home » Archimedes archive » Archimedes World » AW_Vol15_Issue7 CD ROM Index Disc.adf » coders » vsn1/c/main

vsn1/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_Vol15_Issue7 CD ROM Index Disc.adf » coders
Filename: vsn1/c/main
Read OK:
File size: 0F89 bytes
Load address: 0000
Exec address: 0000
File contents
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "bbc.h"
#include "os.h"

#include "ass.h"

#define OS_Byte 6
#define PI 3.1415926536
#define lone 16
#define one (1<<lone)
#define fptod(i) ( (int)(0.5+(int)(i)/6.5536)/10000.0 )
#define dtofp(i) ( (int)(0.5+(double)(i)*one) )
typedef enum {multi, nonmulti} monitor;

int maxsx=1280;
int maxsy=1024;

int e[3]={dtofp(1.2F),dtofp(1.3F),dtofp(1.6F)};
int d[3]={0,0,0};
double theta=0, ctheta=1, stheta=0;
int mx[3][3]={one,0,0, 0,one,0, 0,0,one};
int ppd=6*one, xsc=3*one;

int bank=1;
int mode;

void press(int i)
{
  for (;bbc_inkey(i)!=-1;);
  return;
}

void release(int i)
{
  for (;bbc_inkey(i)!=0;);
  return;
}

monitor read_monitor_type(void)
{
  int r1;
  os_swi3r(0x58, 1,0,0, 0,&r1,0);
  if (r1 != 0) return multi;
  else return nonmulti;
}

void tidy(void)
{
  os_swi2(OS_Byte, 112, 1);
  os_swi2(OS_Byte, 113, 1);
}

void swap_banks(void)
{
  os_swi2(OS_Byte, 113, bank);
  os_swi2(OS_Byte, 112, bank=3-bank);
}

void prepare_screen(void)
{
  bank=1;
  os_swi2(OS_Byte, 112, bank);
  os_swi2(OS_Byte, 113, 3-bank);
}

void calc_trans(int *e, int *d, int mx[][3])
{
  double f[3], m[2][3];
  double r2, r, s;
  double a=stheta, b=ctheta;
  f[0]=fptod(d[0]-e[0]);	f[1]=fptod(d[1]-e[1]);	f[2]=fptod(d[2]-e[2]);
  r2=f[0]*f[0]+f[2]*f[2];	r=sqrt(r2);		s=sqrt(r2+f[1]*f[1]);  
  if (!s) {
    mx[0][0]=dtofp(b);				mx[0][1]=dtofp(-a);			mx[0][2]=0;
    mx[1][0]=dtofp(a);				mx[1][1]=dtofp(b);			mx[1][2]=0;
    mx[2][0]=0;					mx[2][1]=0;				mx[2][2]=one;
  }
  else
    if (!r) {
      s=1/s;
      mx[0][0]=dtofp(b);			mx[0][1]=dtofp(-a*r*s);			mx[0][2]=dtofp(a*f[1]*s);
      mx[1][0]=dtofp(a);			mx[1][1]=dtofp(b*r*s);			mx[1][2]=dtofp(-b*f[1]*s);
      mx[2][0]=0;				mx[2][1]=dtofp(f[1]*s);			mx[2][2]=dtofp(r*s);
    }
    else {
      s=1/s; r2=1/r;
      m[0][0]=r2*f[2];		m[0][1]=0;	m[0][2]=-r2*f[0];
      m[1][0]=-r2*s*f[0]*f[1];	m[1][1]=s*r;	m[1][2]=-r2*s*f[1]*f[2];
      mx[2][0]=dtofp(s*f[0]);			mx[2][1]=dtofp(s*f[1]);			mx[2][2]=dtofp(s*f[2]);
      mx[0][0]=dtofp(b*m[0][0]-a*m[1][0]);	mx[0][1]=dtofp(b*m[0][1]-a*m[1][1]);	mx[0][2]=dtofp(b*m[0][2]-a*m[1][2]);
      mx[1][0]=dtofp(a*m[0][0]+b*m[1][0]);	mx[1][1]=dtofp(a*m[0][1]+b*m[1][1]);	mx[1][2]=dtofp(a*m[0][2]+b*m[1][2]);
    }
}

int project_vertex(int *v, int *s, int mx[][3])
{
  int c[3];
  int b[3], m;
  int x, y, flag;
  c[0] = *v++-e[0]; c[1] = *v++-e[1]; c[2] = *v-e[2];
  b[0] = mul_frac16c(c[0],mx[0][0]) + mul_frac16c(c[1],mx[0][1]) + mul_frac16c(c[2],mx[0][2]);
  b[1] = mul_frac16c(c[0],mx[1][0]) + mul_frac16c(c[1],mx[1][1]) + mul_frac16c(c[2],mx[1][2]);
  b[2] = mul_frac16c(c[0],mx[2][0]) + mul_frac16c(c[1],mx[2][1]) + mul_frac16c(c[2],mx[2][2]);
  m = div_frac16(ppd*640, mul_frac16(b[2],xsc));
  x = (mul_frac16(-m,b[0])+one/2)/one+640;
  y = (mul_frac16(m,b[1])+one/2)/one+512;
  flag = b[2]<one/16 ? 16 : 0;
  if (x<0) flag|=1;
  if (x>maxsx) flag|=2;
  if (y<0) flag|=4;
  if (y>maxsy) flag|=8;
  *s++=x; *s=y;      
  return flag;
}

int wav(int x, int y, int t)
{
  int r=sqrt_frac16(mul_frac16(x,x)+mul_frac16(y,y))-t;
  return mul_frac16c(x*5, sin16(r*32)/8);
}

int main(void)
{
  int s[2], p[3];
  int i, j, t, di;

  mode = read_monitor_type()==multi ? 18 : 0;
  atexit(tidy);
  tidy();
  bbc_mode(mode+128);
  bbc_mode(mode);
  printf("Please enter dot separation (suggest 1-8) ");
  scanf("%i", &di);
  prepare_screen();

  e[0]=dtofp(1.2F);
  e[1]=dtofp(0.3F);
  e[2]=dtofp(0.9F);
  d[0]=dtofp(0.0F);
  d[1]=dtofp(0.01F);
  d[2]=dtofp(0.1F);

  calc_trans(e, d, mx);
  
  for (t=8*one/128;;) {
    for (; t<one/8; t+=one/128) {
      bbc_cls(); 
      for (i=0; i<140; i+=di)
        for (j=0; j<210; j+=di) {
          p[0]=wav((210-j)*one/1024, (i-70)*one/1024, t);
          p[1]=(i-70)*one/160;
          p[2]=(j-105)*one/160;
          project_vertex(p, s, mx);
          bbc_plot(69, s[0], s[1]);
        }
      os_swi1(OS_Byte, 19);
      swap_banks();
    }
    t = 0;
  }

  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 6d 61 74 68 2e 68 3e  0a 23 69 6e 63 6c 75 64  |<math.h>.#includ|
00000040  65 20 22 62 62 63 2e 68  22 0a 23 69 6e 63 6c 75  |e "bbc.h".#inclu|
00000050  64 65 20 22 6f 73 2e 68  22 0a 0a 23 69 6e 63 6c  |de "os.h"..#incl|
00000060  75 64 65 20 22 61 73 73  2e 68 22 0a 0a 23 64 65  |ude "ass.h"..#de|
00000070  66 69 6e 65 20 4f 53 5f  42 79 74 65 20 36 0a 23  |fine OS_Byte 6.#|
00000080  64 65 66 69 6e 65 20 50  49 20 33 2e 31 34 31 35  |define PI 3.1415|
00000090  39 32 36 35 33 36 0a 23  64 65 66 69 6e 65 20 6c  |926536.#define l|
000000a0  6f 6e 65 20 31 36 0a 23  64 65 66 69 6e 65 20 6f  |one 16.#define o|
000000b0  6e 65 20 28 31 3c 3c 6c  6f 6e 65 29 0a 23 64 65  |ne (1<<lone).#de|
000000c0  66 69 6e 65 20 66 70 74  6f 64 28 69 29 20 28 20  |fine fptod(i) ( |
000000d0  28 69 6e 74 29 28 30 2e  35 2b 28 69 6e 74 29 28  |(int)(0.5+(int)(|
000000e0  69 29 2f 36 2e 35 35 33  36 29 2f 31 30 30 30 30  |i)/6.5536)/10000|
000000f0  2e 30 20 29 0a 23 64 65  66 69 6e 65 20 64 74 6f  |.0 ).#define dto|
00000100  66 70 28 69 29 20 28 20  28 69 6e 74 29 28 30 2e  |fp(i) ( (int)(0.|
00000110  35 2b 28 64 6f 75 62 6c  65 29 28 69 29 2a 6f 6e  |5+(double)(i)*on|
00000120  65 29 20 29 0a 74 79 70  65 64 65 66 20 65 6e 75  |e) ).typedef enu|
00000130  6d 20 7b 6d 75 6c 74 69  2c 20 6e 6f 6e 6d 75 6c  |m {multi, nonmul|
00000140  74 69 7d 20 6d 6f 6e 69  74 6f 72 3b 0a 0a 69 6e  |ti} monitor;..in|
00000150  74 20 6d 61 78 73 78 3d  31 32 38 30 3b 0a 69 6e  |t maxsx=1280;.in|
00000160  74 20 6d 61 78 73 79 3d  31 30 32 34 3b 0a 0a 69  |t maxsy=1024;..i|
00000170  6e 74 20 65 5b 33 5d 3d  7b 64 74 6f 66 70 28 31  |nt e[3]={dtofp(1|
00000180  2e 32 46 29 2c 64 74 6f  66 70 28 31 2e 33 46 29  |.2F),dtofp(1.3F)|
00000190  2c 64 74 6f 66 70 28 31  2e 36 46 29 7d 3b 0a 69  |,dtofp(1.6F)};.i|
000001a0  6e 74 20 64 5b 33 5d 3d  7b 30 2c 30 2c 30 7d 3b  |nt d[3]={0,0,0};|
000001b0  0a 64 6f 75 62 6c 65 20  74 68 65 74 61 3d 30 2c  |.double theta=0,|
000001c0  20 63 74 68 65 74 61 3d  31 2c 20 73 74 68 65 74  | ctheta=1, sthet|
000001d0  61 3d 30 3b 0a 69 6e 74  20 6d 78 5b 33 5d 5b 33  |a=0;.int mx[3][3|
000001e0  5d 3d 7b 6f 6e 65 2c 30  2c 30 2c 20 30 2c 6f 6e  |]={one,0,0, 0,on|
000001f0  65 2c 30 2c 20 30 2c 30  2c 6f 6e 65 7d 3b 0a 69  |e,0, 0,0,one};.i|
00000200  6e 74 20 70 70 64 3d 36  2a 6f 6e 65 2c 20 78 73  |nt ppd=6*one, xs|
00000210  63 3d 33 2a 6f 6e 65 3b  0a 0a 69 6e 74 20 62 61  |c=3*one;..int ba|
00000220  6e 6b 3d 31 3b 0a 69 6e  74 20 6d 6f 64 65 3b 0a  |nk=1;.int mode;.|
00000230  0a 76 6f 69 64 20 70 72  65 73 73 28 69 6e 74 20  |.void press(int |
00000240  69 29 0a 7b 0a 20 20 66  6f 72 20 28 3b 62 62 63  |i).{.  for (;bbc|
00000250  5f 69 6e 6b 65 79 28 69  29 21 3d 2d 31 3b 29 3b  |_inkey(i)!=-1;);|
00000260  0a 20 20 72 65 74 75 72  6e 3b 0a 7d 0a 0a 76 6f  |.  return;.}..vo|
00000270  69 64 20 72 65 6c 65 61  73 65 28 69 6e 74 20 69  |id release(int i|
00000280  29 0a 7b 0a 20 20 66 6f  72 20 28 3b 62 62 63 5f  |).{.  for (;bbc_|
00000290  69 6e 6b 65 79 28 69 29  21 3d 30 3b 29 3b 0a 20  |inkey(i)!=0;);. |
000002a0  20 72 65 74 75 72 6e 3b  0a 7d 0a 0a 6d 6f 6e 69  | return;.}..moni|
000002b0  74 6f 72 20 72 65 61 64  5f 6d 6f 6e 69 74 6f 72  |tor read_monitor|
000002c0  5f 74 79 70 65 28 76 6f  69 64 29 0a 7b 0a 20 20  |_type(void).{.  |
000002d0  69 6e 74 20 72 31 3b 0a  20 20 6f 73 5f 73 77 69  |int r1;.  os_swi|
000002e0  33 72 28 30 78 35 38 2c  20 31 2c 30 2c 30 2c 20  |3r(0x58, 1,0,0, |
000002f0  30 2c 26 72 31 2c 30 29  3b 0a 20 20 69 66 20 28  |0,&r1,0);.  if (|
00000300  72 31 20 21 3d 20 30 29  20 72 65 74 75 72 6e 20  |r1 != 0) return |
00000310  6d 75 6c 74 69 3b 0a 20  20 65 6c 73 65 20 72 65  |multi;.  else re|
00000320  74 75 72 6e 20 6e 6f 6e  6d 75 6c 74 69 3b 0a 7d  |turn nonmulti;.}|
00000330  0a 0a 76 6f 69 64 20 74  69 64 79 28 76 6f 69 64  |..void tidy(void|
00000340  29 0a 7b 0a 20 20 6f 73  5f 73 77 69 32 28 4f 53  |).{.  os_swi2(OS|
00000350  5f 42 79 74 65 2c 20 31  31 32 2c 20 31 29 3b 0a  |_Byte, 112, 1);.|
00000360  20 20 6f 73 5f 73 77 69  32 28 4f 53 5f 42 79 74  |  os_swi2(OS_Byt|
00000370  65 2c 20 31 31 33 2c 20  31 29 3b 0a 7d 0a 0a 76  |e, 113, 1);.}..v|
00000380  6f 69 64 20 73 77 61 70  5f 62 61 6e 6b 73 28 76  |oid swap_banks(v|
00000390  6f 69 64 29 0a 7b 0a 20  20 6f 73 5f 73 77 69 32  |oid).{.  os_swi2|
000003a0  28 4f 53 5f 42 79 74 65  2c 20 31 31 33 2c 20 62  |(OS_Byte, 113, b|
000003b0  61 6e 6b 29 3b 0a 20 20  6f 73 5f 73 77 69 32 28  |ank);.  os_swi2(|
000003c0  4f 53 5f 42 79 74 65 2c  20 31 31 32 2c 20 62 61  |OS_Byte, 112, ba|
000003d0  6e 6b 3d 33 2d 62 61 6e  6b 29 3b 0a 7d 0a 0a 76  |nk=3-bank);.}..v|
000003e0  6f 69 64 20 70 72 65 70  61 72 65 5f 73 63 72 65  |oid prepare_scre|
000003f0  65 6e 28 76 6f 69 64 29  0a 7b 0a 20 20 62 61 6e  |en(void).{.  ban|
00000400  6b 3d 31 3b 0a 20 20 6f  73 5f 73 77 69 32 28 4f  |k=1;.  os_swi2(O|
00000410  53 5f 42 79 74 65 2c 20  31 31 32 2c 20 62 61 6e  |S_Byte, 112, ban|
00000420  6b 29 3b 0a 20 20 6f 73  5f 73 77 69 32 28 4f 53  |k);.  os_swi2(OS|
00000430  5f 42 79 74 65 2c 20 31  31 33 2c 20 33 2d 62 61  |_Byte, 113, 3-ba|
00000440  6e 6b 29 3b 0a 7d 0a 0a  76 6f 69 64 20 63 61 6c  |nk);.}..void cal|
00000450  63 5f 74 72 61 6e 73 28  69 6e 74 20 2a 65 2c 20  |c_trans(int *e, |
00000460  69 6e 74 20 2a 64 2c 20  69 6e 74 20 6d 78 5b 5d  |int *d, int mx[]|
00000470  5b 33 5d 29 0a 7b 0a 20  20 64 6f 75 62 6c 65 20  |[3]).{.  double |
00000480  66 5b 33 5d 2c 20 6d 5b  32 5d 5b 33 5d 3b 0a 20  |f[3], m[2][3];. |
00000490  20 64 6f 75 62 6c 65 20  72 32 2c 20 72 2c 20 73  | double r2, r, s|
000004a0  3b 0a 20 20 64 6f 75 62  6c 65 20 61 3d 73 74 68  |;.  double a=sth|
000004b0  65 74 61 2c 20 62 3d 63  74 68 65 74 61 3b 0a 20  |eta, b=ctheta;. |
000004c0  20 66 5b 30 5d 3d 66 70  74 6f 64 28 64 5b 30 5d  | f[0]=fptod(d[0]|
000004d0  2d 65 5b 30 5d 29 3b 09  66 5b 31 5d 3d 66 70 74  |-e[0]);.f[1]=fpt|
000004e0  6f 64 28 64 5b 31 5d 2d  65 5b 31 5d 29 3b 09 66  |od(d[1]-e[1]);.f|
000004f0  5b 32 5d 3d 66 70 74 6f  64 28 64 5b 32 5d 2d 65  |[2]=fptod(d[2]-e|
00000500  5b 32 5d 29 3b 0a 20 20  72 32 3d 66 5b 30 5d 2a  |[2]);.  r2=f[0]*|
00000510  66 5b 30 5d 2b 66 5b 32  5d 2a 66 5b 32 5d 3b 09  |f[0]+f[2]*f[2];.|
00000520  72 3d 73 71 72 74 28 72  32 29 3b 09 09 73 3d 73  |r=sqrt(r2);..s=s|
00000530  71 72 74 28 72 32 2b 66  5b 31 5d 2a 66 5b 31 5d  |qrt(r2+f[1]*f[1]|
00000540  29 3b 20 20 0a 20 20 69  66 20 28 21 73 29 20 7b  |);  .  if (!s) {|
00000550  0a 20 20 20 20 6d 78 5b  30 5d 5b 30 5d 3d 64 74  |.    mx[0][0]=dt|
00000560  6f 66 70 28 62 29 3b 09  09 09 09 6d 78 5b 30 5d  |ofp(b);....mx[0]|
00000570  5b 31 5d 3d 64 74 6f 66  70 28 2d 61 29 3b 09 09  |[1]=dtofp(-a);..|
00000580  09 6d 78 5b 30 5d 5b 32  5d 3d 30 3b 0a 20 20 20  |.mx[0][2]=0;.   |
00000590  20 6d 78 5b 31 5d 5b 30  5d 3d 64 74 6f 66 70 28  | mx[1][0]=dtofp(|
000005a0  61 29 3b 09 09 09 09 6d  78 5b 31 5d 5b 31 5d 3d  |a);....mx[1][1]=|
000005b0  64 74 6f 66 70 28 62 29  3b 09 09 09 6d 78 5b 31  |dtofp(b);...mx[1|
000005c0  5d 5b 32 5d 3d 30 3b 0a  20 20 20 20 6d 78 5b 32  |][2]=0;.    mx[2|
000005d0  5d 5b 30 5d 3d 30 3b 09  09 09 09 09 6d 78 5b 32  |][0]=0;.....mx[2|
000005e0  5d 5b 31 5d 3d 30 3b 09  09 09 09 6d 78 5b 32 5d  |][1]=0;....mx[2]|
000005f0  5b 32 5d 3d 6f 6e 65 3b  0a 20 20 7d 0a 20 20 65  |[2]=one;.  }.  e|
00000600  6c 73 65 0a 20 20 20 20  69 66 20 28 21 72 29 20  |lse.    if (!r) |
00000610  7b 0a 20 20 20 20 20 20  73 3d 31 2f 73 3b 0a 20  |{.      s=1/s;. |
00000620  20 20 20 20 20 6d 78 5b  30 5d 5b 30 5d 3d 64 74  |     mx[0][0]=dt|
00000630  6f 66 70 28 62 29 3b 09  09 09 6d 78 5b 30 5d 5b  |ofp(b);...mx[0][|
00000640  31 5d 3d 64 74 6f 66 70  28 2d 61 2a 72 2a 73 29  |1]=dtofp(-a*r*s)|
00000650  3b 09 09 09 6d 78 5b 30  5d 5b 32 5d 3d 64 74 6f  |;...mx[0][2]=dto|
00000660  66 70 28 61 2a 66 5b 31  5d 2a 73 29 3b 0a 20 20  |fp(a*f[1]*s);.  |
00000670  20 20 20 20 6d 78 5b 31  5d 5b 30 5d 3d 64 74 6f  |    mx[1][0]=dto|
00000680  66 70 28 61 29 3b 09 09  09 6d 78 5b 31 5d 5b 31  |fp(a);...mx[1][1|
00000690  5d 3d 64 74 6f 66 70 28  62 2a 72 2a 73 29 3b 09  |]=dtofp(b*r*s);.|
000006a0  09 09 6d 78 5b 31 5d 5b  32 5d 3d 64 74 6f 66 70  |..mx[1][2]=dtofp|
000006b0  28 2d 62 2a 66 5b 31 5d  2a 73 29 3b 0a 20 20 20  |(-b*f[1]*s);.   |
000006c0  20 20 20 6d 78 5b 32 5d  5b 30 5d 3d 30 3b 09 09  |   mx[2][0]=0;..|
000006d0  09 09 6d 78 5b 32 5d 5b  31 5d 3d 64 74 6f 66 70  |..mx[2][1]=dtofp|
000006e0  28 66 5b 31 5d 2a 73 29  3b 09 09 09 6d 78 5b 32  |(f[1]*s);...mx[2|
000006f0  5d 5b 32 5d 3d 64 74 6f  66 70 28 72 2a 73 29 3b  |][2]=dtofp(r*s);|
00000700  0a 20 20 20 20 7d 0a 20  20 20 20 65 6c 73 65 20  |.    }.    else |
00000710  7b 0a 20 20 20 20 20 20  73 3d 31 2f 73 3b 20 72  |{.      s=1/s; r|
00000720  32 3d 31 2f 72 3b 0a 20  20 20 20 20 20 6d 5b 30  |2=1/r;.      m[0|
00000730  5d 5b 30 5d 3d 72 32 2a  66 5b 32 5d 3b 09 09 6d  |][0]=r2*f[2];..m|
00000740  5b 30 5d 5b 31 5d 3d 30  3b 09 6d 5b 30 5d 5b 32  |[0][1]=0;.m[0][2|
00000750  5d 3d 2d 72 32 2a 66 5b  30 5d 3b 0a 20 20 20 20  |]=-r2*f[0];.    |
00000760  20 20 6d 5b 31 5d 5b 30  5d 3d 2d 72 32 2a 73 2a  |  m[1][0]=-r2*s*|
00000770  66 5b 30 5d 2a 66 5b 31  5d 3b 09 6d 5b 31 5d 5b  |f[0]*f[1];.m[1][|
00000780  31 5d 3d 73 2a 72 3b 09  6d 5b 31 5d 5b 32 5d 3d  |1]=s*r;.m[1][2]=|
00000790  2d 72 32 2a 73 2a 66 5b  31 5d 2a 66 5b 32 5d 3b  |-r2*s*f[1]*f[2];|
000007a0  0a 20 20 20 20 20 20 6d  78 5b 32 5d 5b 30 5d 3d  |.      mx[2][0]=|
000007b0  64 74 6f 66 70 28 73 2a  66 5b 30 5d 29 3b 09 09  |dtofp(s*f[0]);..|
000007c0  09 6d 78 5b 32 5d 5b 31  5d 3d 64 74 6f 66 70 28  |.mx[2][1]=dtofp(|
000007d0  73 2a 66 5b 31 5d 29 3b  09 09 09 6d 78 5b 32 5d  |s*f[1]);...mx[2]|
000007e0  5b 32 5d 3d 64 74 6f 66  70 28 73 2a 66 5b 32 5d  |[2]=dtofp(s*f[2]|
000007f0  29 3b 0a 20 20 20 20 20  20 6d 78 5b 30 5d 5b 30  |);.      mx[0][0|
00000800  5d 3d 64 74 6f 66 70 28  62 2a 6d 5b 30 5d 5b 30  |]=dtofp(b*m[0][0|
00000810  5d 2d 61 2a 6d 5b 31 5d  5b 30 5d 29 3b 09 6d 78  |]-a*m[1][0]);.mx|
00000820  5b 30 5d 5b 31 5d 3d 64  74 6f 66 70 28 62 2a 6d  |[0][1]=dtofp(b*m|
00000830  5b 30 5d 5b 31 5d 2d 61  2a 6d 5b 31 5d 5b 31 5d  |[0][1]-a*m[1][1]|
00000840  29 3b 09 6d 78 5b 30 5d  5b 32 5d 3d 64 74 6f 66  |);.mx[0][2]=dtof|
00000850  70 28 62 2a 6d 5b 30 5d  5b 32 5d 2d 61 2a 6d 5b  |p(b*m[0][2]-a*m[|
00000860  31 5d 5b 32 5d 29 3b 0a  20 20 20 20 20 20 6d 78  |1][2]);.      mx|
00000870  5b 31 5d 5b 30 5d 3d 64  74 6f 66 70 28 61 2a 6d  |[1][0]=dtofp(a*m|
00000880  5b 30 5d 5b 30 5d 2b 62  2a 6d 5b 31 5d 5b 30 5d  |[0][0]+b*m[1][0]|
00000890  29 3b 09 6d 78 5b 31 5d  5b 31 5d 3d 64 74 6f 66  |);.mx[1][1]=dtof|
000008a0  70 28 61 2a 6d 5b 30 5d  5b 31 5d 2b 62 2a 6d 5b  |p(a*m[0][1]+b*m[|
000008b0  31 5d 5b 31 5d 29 3b 09  6d 78 5b 31 5d 5b 32 5d  |1][1]);.mx[1][2]|
000008c0  3d 64 74 6f 66 70 28 61  2a 6d 5b 30 5d 5b 32 5d  |=dtofp(a*m[0][2]|
000008d0  2b 62 2a 6d 5b 31 5d 5b  32 5d 29 3b 0a 20 20 20  |+b*m[1][2]);.   |
000008e0  20 7d 0a 7d 0a 0a 69 6e  74 20 70 72 6f 6a 65 63  | }.}..int projec|
000008f0  74 5f 76 65 72 74 65 78  28 69 6e 74 20 2a 76 2c  |t_vertex(int *v,|
00000900  20 69 6e 74 20 2a 73 2c  20 69 6e 74 20 6d 78 5b  | int *s, int mx[|
00000910  5d 5b 33 5d 29 0a 7b 0a  20 20 69 6e 74 20 63 5b  |][3]).{.  int c[|
00000920  33 5d 3b 0a 20 20 69 6e  74 20 62 5b 33 5d 2c 20  |3];.  int b[3], |
00000930  6d 3b 0a 20 20 69 6e 74  20 78 2c 20 79 2c 20 66  |m;.  int x, y, f|
00000940  6c 61 67 3b 0a 20 20 63  5b 30 5d 20 3d 20 2a 76  |lag;.  c[0] = *v|
00000950  2b 2b 2d 65 5b 30 5d 3b  20 63 5b 31 5d 20 3d 20  |++-e[0]; c[1] = |
00000960  2a 76 2b 2b 2d 65 5b 31  5d 3b 20 63 5b 32 5d 20  |*v++-e[1]; c[2] |
00000970  3d 20 2a 76 2d 65 5b 32  5d 3b 0a 20 20 62 5b 30  |= *v-e[2];.  b[0|
00000980  5d 20 3d 20 6d 75 6c 5f  66 72 61 63 31 36 63 28  |] = mul_frac16c(|
00000990  63 5b 30 5d 2c 6d 78 5b  30 5d 5b 30 5d 29 20 2b  |c[0],mx[0][0]) +|
000009a0  20 6d 75 6c 5f 66 72 61  63 31 36 63 28 63 5b 31  | mul_frac16c(c[1|
000009b0  5d 2c 6d 78 5b 30 5d 5b  31 5d 29 20 2b 20 6d 75  |],mx[0][1]) + mu|
000009c0  6c 5f 66 72 61 63 31 36  63 28 63 5b 32 5d 2c 6d  |l_frac16c(c[2],m|
000009d0  78 5b 30 5d 5b 32 5d 29  3b 0a 20 20 62 5b 31 5d  |x[0][2]);.  b[1]|
000009e0  20 3d 20 6d 75 6c 5f 66  72 61 63 31 36 63 28 63  | = mul_frac16c(c|
000009f0  5b 30 5d 2c 6d 78 5b 31  5d 5b 30 5d 29 20 2b 20  |[0],mx[1][0]) + |
00000a00  6d 75 6c 5f 66 72 61 63  31 36 63 28 63 5b 31 5d  |mul_frac16c(c[1]|
00000a10  2c 6d 78 5b 31 5d 5b 31  5d 29 20 2b 20 6d 75 6c  |,mx[1][1]) + mul|
00000a20  5f 66 72 61 63 31 36 63  28 63 5b 32 5d 2c 6d 78  |_frac16c(c[2],mx|
00000a30  5b 31 5d 5b 32 5d 29 3b  0a 20 20 62 5b 32 5d 20  |[1][2]);.  b[2] |
00000a40  3d 20 6d 75 6c 5f 66 72  61 63 31 36 63 28 63 5b  |= mul_frac16c(c[|
00000a50  30 5d 2c 6d 78 5b 32 5d  5b 30 5d 29 20 2b 20 6d  |0],mx[2][0]) + m|
00000a60  75 6c 5f 66 72 61 63 31  36 63 28 63 5b 31 5d 2c  |ul_frac16c(c[1],|
00000a70  6d 78 5b 32 5d 5b 31 5d  29 20 2b 20 6d 75 6c 5f  |mx[2][1]) + mul_|
00000a80  66 72 61 63 31 36 63 28  63 5b 32 5d 2c 6d 78 5b  |frac16c(c[2],mx[|
00000a90  32 5d 5b 32 5d 29 3b 0a  20 20 6d 20 3d 20 64 69  |2][2]);.  m = di|
00000aa0  76 5f 66 72 61 63 31 36  28 70 70 64 2a 36 34 30  |v_frac16(ppd*640|
00000ab0  2c 20 6d 75 6c 5f 66 72  61 63 31 36 28 62 5b 32  |, mul_frac16(b[2|
00000ac0  5d 2c 78 73 63 29 29 3b  0a 20 20 78 20 3d 20 28  |],xsc));.  x = (|
00000ad0  6d 75 6c 5f 66 72 61 63  31 36 28 2d 6d 2c 62 5b  |mul_frac16(-m,b[|
00000ae0  30 5d 29 2b 6f 6e 65 2f  32 29 2f 6f 6e 65 2b 36  |0])+one/2)/one+6|
00000af0  34 30 3b 0a 20 20 79 20  3d 20 28 6d 75 6c 5f 66  |40;.  y = (mul_f|
00000b00  72 61 63 31 36 28 6d 2c  62 5b 31 5d 29 2b 6f 6e  |rac16(m,b[1])+on|
00000b10  65 2f 32 29 2f 6f 6e 65  2b 35 31 32 3b 0a 20 20  |e/2)/one+512;.  |
00000b20  66 6c 61 67 20 3d 20 62  5b 32 5d 3c 6f 6e 65 2f  |flag = b[2]<one/|
00000b30  31 36 20 3f 20 31 36 20  3a 20 30 3b 0a 20 20 69  |16 ? 16 : 0;.  i|
00000b40  66 20 28 78 3c 30 29 20  66 6c 61 67 7c 3d 31 3b  |f (x<0) flag|=1;|
00000b50  0a 20 20 69 66 20 28 78  3e 6d 61 78 73 78 29 20  |.  if (x>maxsx) |
00000b60  66 6c 61 67 7c 3d 32 3b  0a 20 20 69 66 20 28 79  |flag|=2;.  if (y|
00000b70  3c 30 29 20 66 6c 61 67  7c 3d 34 3b 0a 20 20 69  |<0) flag|=4;.  i|
00000b80  66 20 28 79 3e 6d 61 78  73 79 29 20 66 6c 61 67  |f (y>maxsy) flag|
00000b90  7c 3d 38 3b 0a 20 20 2a  73 2b 2b 3d 78 3b 20 2a  ||=8;.  *s++=x; *|
00000ba0  73 3d 79 3b 20 20 20 20  20 20 0a 20 20 72 65 74  |s=y;      .  ret|
00000bb0  75 72 6e 20 66 6c 61 67  3b 0a 7d 0a 0a 69 6e 74  |urn flag;.}..int|
00000bc0  20 77 61 76 28 69 6e 74  20 78 2c 20 69 6e 74 20  | wav(int x, int |
00000bd0  79 2c 20 69 6e 74 20 74  29 0a 7b 0a 20 20 69 6e  |y, int t).{.  in|
00000be0  74 20 72 3d 73 71 72 74  5f 66 72 61 63 31 36 28  |t r=sqrt_frac16(|
00000bf0  6d 75 6c 5f 66 72 61 63  31 36 28 78 2c 78 29 2b  |mul_frac16(x,x)+|
00000c00  6d 75 6c 5f 66 72 61 63  31 36 28 79 2c 79 29 29  |mul_frac16(y,y))|
00000c10  2d 74 3b 0a 20 20 72 65  74 75 72 6e 20 6d 75 6c  |-t;.  return mul|
00000c20  5f 66 72 61 63 31 36 63  28 78 2a 35 2c 20 73 69  |_frac16c(x*5, si|
00000c30  6e 31 36 28 72 2a 33 32  29 2f 38 29 3b 0a 7d 0a  |n16(r*32)/8);.}.|
00000c40  0a 69 6e 74 20 6d 61 69  6e 28 76 6f 69 64 29 0a  |.int main(void).|
00000c50  7b 0a 20 20 69 6e 74 20  73 5b 32 5d 2c 20 70 5b  |{.  int s[2], p[|
00000c60  33 5d 3b 0a 20 20 69 6e  74 20 69 2c 20 6a 2c 20  |3];.  int i, j, |
00000c70  74 2c 20 64 69 3b 0a 0a  20 20 6d 6f 64 65 20 3d  |t, di;..  mode =|
00000c80  20 72 65 61 64 5f 6d 6f  6e 69 74 6f 72 5f 74 79  | read_monitor_ty|
00000c90  70 65 28 29 3d 3d 6d 75  6c 74 69 20 3f 20 31 38  |pe()==multi ? 18|
00000ca0  20 3a 20 30 3b 0a 20 20  61 74 65 78 69 74 28 74  | : 0;.  atexit(t|
00000cb0  69 64 79 29 3b 0a 20 20  74 69 64 79 28 29 3b 0a  |idy);.  tidy();.|
00000cc0  20 20 62 62 63 5f 6d 6f  64 65 28 6d 6f 64 65 2b  |  bbc_mode(mode+|
00000cd0  31 32 38 29 3b 0a 20 20  62 62 63 5f 6d 6f 64 65  |128);.  bbc_mode|
00000ce0  28 6d 6f 64 65 29 3b 0a  20 20 70 72 69 6e 74 66  |(mode);.  printf|
00000cf0  28 22 50 6c 65 61 73 65  20 65 6e 74 65 72 20 64  |("Please enter d|
00000d00  6f 74 20 73 65 70 61 72  61 74 69 6f 6e 20 28 73  |ot separation (s|
00000d10  75 67 67 65 73 74 20 31  2d 38 29 20 22 29 3b 0a  |uggest 1-8) ");.|
00000d20  20 20 73 63 61 6e 66 28  22 25 69 22 2c 20 26 64  |  scanf("%i", &d|
00000d30  69 29 3b 0a 20 20 70 72  65 70 61 72 65 5f 73 63  |i);.  prepare_sc|
00000d40  72 65 65 6e 28 29 3b 0a  0a 20 20 65 5b 30 5d 3d  |reen();..  e[0]=|
00000d50  64 74 6f 66 70 28 31 2e  32 46 29 3b 0a 20 20 65  |dtofp(1.2F);.  e|
00000d60  5b 31 5d 3d 64 74 6f 66  70 28 30 2e 33 46 29 3b  |[1]=dtofp(0.3F);|
00000d70  0a 20 20 65 5b 32 5d 3d  64 74 6f 66 70 28 30 2e  |.  e[2]=dtofp(0.|
00000d80  39 46 29 3b 0a 20 20 64  5b 30 5d 3d 64 74 6f 66  |9F);.  d[0]=dtof|
00000d90  70 28 30 2e 30 46 29 3b  0a 20 20 64 5b 31 5d 3d  |p(0.0F);.  d[1]=|
00000da0  64 74 6f 66 70 28 30 2e  30 31 46 29 3b 0a 20 20  |dtofp(0.01F);.  |
00000db0  64 5b 32 5d 3d 64 74 6f  66 70 28 30 2e 31 46 29  |d[2]=dtofp(0.1F)|
00000dc0  3b 0a 0a 20 20 63 61 6c  63 5f 74 72 61 6e 73 28  |;..  calc_trans(|
00000dd0  65 2c 20 64 2c 20 6d 78  29 3b 0a 20 20 0a 20 20  |e, d, mx);.  .  |
00000de0  66 6f 72 20 28 74 3d 38  2a 6f 6e 65 2f 31 32 38  |for (t=8*one/128|
00000df0  3b 3b 29 20 7b 0a 20 20  20 20 66 6f 72 20 28 3b  |;;) {.    for (;|
00000e00  20 74 3c 6f 6e 65 2f 38  3b 20 74 2b 3d 6f 6e 65  | t<one/8; t+=one|
00000e10  2f 31 32 38 29 20 7b 0a  20 20 20 20 20 20 62 62  |/128) {.      bb|
00000e20  63 5f 63 6c 73 28 29 3b  20 0a 20 20 20 20 20 20  |c_cls(); .      |
00000e30  66 6f 72 20 28 69 3d 30  3b 20 69 3c 31 34 30 3b  |for (i=0; i<140;|
00000e40  20 69 2b 3d 64 69 29 0a  20 20 20 20 20 20 20 20  | i+=di).        |
00000e50  66 6f 72 20 28 6a 3d 30  3b 20 6a 3c 32 31 30 3b  |for (j=0; j<210;|
00000e60  20 6a 2b 3d 64 69 29 20  7b 0a 20 20 20 20 20 20  | j+=di) {.      |
00000e70  20 20 20 20 70 5b 30 5d  3d 77 61 76 28 28 32 31  |    p[0]=wav((21|
00000e80  30 2d 6a 29 2a 6f 6e 65  2f 31 30 32 34 2c 20 28  |0-j)*one/1024, (|
00000e90  69 2d 37 30 29 2a 6f 6e  65 2f 31 30 32 34 2c 20  |i-70)*one/1024, |
00000ea0  74 29 3b 0a 20 20 20 20  20 20 20 20 20 20 70 5b  |t);.          p[|
00000eb0  31 5d 3d 28 69 2d 37 30  29 2a 6f 6e 65 2f 31 36  |1]=(i-70)*one/16|
00000ec0  30 3b 0a 20 20 20 20 20  20 20 20 20 20 70 5b 32  |0;.          p[2|
00000ed0  5d 3d 28 6a 2d 31 30 35  29 2a 6f 6e 65 2f 31 36  |]=(j-105)*one/16|
00000ee0  30 3b 0a 20 20 20 20 20  20 20 20 20 20 70 72 6f  |0;.          pro|
00000ef0  6a 65 63 74 5f 76 65 72  74 65 78 28 70 2c 20 73  |ject_vertex(p, s|
00000f00  2c 20 6d 78 29 3b 0a 20  20 20 20 20 20 20 20 20  |, mx);.         |
00000f10  20 62 62 63 5f 70 6c 6f  74 28 36 39 2c 20 73 5b  | bbc_plot(69, s[|
00000f20  30 5d 2c 20 73 5b 31 5d  29 3b 0a 20 20 20 20 20  |0], s[1]);.     |
00000f30  20 20 20 7d 0a 20 20 20  20 20 20 6f 73 5f 73 77  |   }.      os_sw|
00000f40  69 31 28 4f 53 5f 42 79  74 65 2c 20 31 39 29 3b  |i1(OS_Byte, 19);|
00000f50  0a 20 20 20 20 20 20 73  77 61 70 5f 62 61 6e 6b  |.      swap_bank|
00000f60  73 28 29 3b 0a 20 20 20  20 7d 0a 20 20 20 20 74  |s();.    }.    t|
00000f70  20 3d 20 30 3b 0a 20 20  7d 0a 0a 20 20 72 65 74  | = 0;.  }..  ret|
00000f80  75 72 6e 20 30 3b 0a 7d  0a                       |urn 0;.}.|
00000f89