home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: SHell self-extracting ARchive
(archive/shar).
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| SHell self-extracting ARchive (archive/shar)
| magic
| Supported |
| 1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
| 100%
| file
| ASCII text
| default
| |
| 100%
| checkBytes
| Printable ASCII
| default
| |
| 100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
| 100%
| siegfried
| fmt/329 Shell Archive Format
| default
| |
| 100%
| detectItEasy
| Format: plain text[LF]
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 6d 69 73 63 0a 58 2d |.sources|.misc.X-|
|00000020| 55 4e 49 58 2d 46 72 6f | 6d 3a 20 63 72 61 69 67 |UNIX-Fro|m: craig|
|00000030| 40 77 65 65 64 65 61 74 | 65 72 2e 6d 61 74 68 2e |@weedeat|er.math.|
|00000040| 79 61 6c 65 2e 65 64 75 | 0a 73 75 62 6a 65 63 74 |yale.edu|.subject|
|00000050| 3a 20 76 31 35 69 30 32 | 35 3a 20 47 72 61 70 68 |: v15i02|5: Graph|
|00000060| 69 63 73 20 47 65 6d 73 | 2c 20 50 61 72 74 20 33 |ics Gems|, Part 3|
|00000070| 2f 35 0a 66 72 6f 6d 3a | 20 43 72 61 69 67 20 4b |/5.from:| Craig K|
|00000080| 6f 6c 62 20 3c 63 72 61 | 69 67 40 77 65 65 64 65 |olb <cra|ig@weede|
|00000090| 61 74 65 72 2e 6d 61 74 | 68 2e 79 61 6c 65 2e 65 |ater.mat|h.yale.e|
|000000a0| 64 75 3e 0a 53 65 6e 64 | 65 72 3a 20 61 6c 6c 62 |du>.Send|er: allb|
|000000b0| 65 72 79 40 75 75 6e 65 | 74 2e 55 55 2e 4e 45 54 |ery@uune|t.UU.NET|
|000000c0| 20 28 42 72 61 6e 64 6f | 6e 20 53 2e 20 41 6c 6c | (Brando|n S. All|
|000000d0| 62 65 72 79 20 2d 20 63 | 6f 6d 70 2e 73 6f 75 72 |bery - c|omp.sour|
|000000e0| 63 65 73 2e 6d 69 73 63 | 29 0a 0a 50 6f 73 74 69 |ces.misc|)..Posti|
|000000f0| 6e 67 2d 6e 75 6d 62 65 | 72 3a 20 56 6f 6c 75 6d |ng-numbe|r: Volum|
|00000100| 65 20 31 35 2c 20 49 73 | 73 75 65 20 32 35 0a 53 |e 15, Is|sue 25.S|
|00000110| 75 62 6d 69 74 74 65 64 | 2d 62 79 3a 20 43 72 61 |ubmitted|-by: Cra|
|00000120| 69 67 20 4b 6f 6c 62 20 | 3c 63 72 61 69 67 40 77 |ig Kolb |<craig@w|
|00000130| 65 65 64 65 61 74 65 72 | 2e 6d 61 74 68 2e 79 61 |eedeater|.math.ya|
|00000140| 6c 65 2e 65 64 75 3e 0a | 41 72 63 68 69 76 65 2d |le.edu>.|Archive-|
|00000150| 6e 61 6d 65 3a 20 67 67 | 65 6d 73 2f 70 61 72 74 |name: gg|ems/part|
|00000160| 30 33 0a 0a 23 21 20 2f | 62 69 6e 2f 73 68 0a 23 |03..#! /|bin/sh.#|
|00000170| 20 54 68 69 73 20 69 73 | 20 61 20 73 68 65 6c 6c | This is| a shell|
|00000180| 20 61 72 63 68 69 76 65 | 2e 20 20 52 65 6d 6f 76 | archive|. Remov|
|00000190| 65 20 61 6e 79 74 68 69 | 6e 67 20 62 65 66 6f 72 |e anythi|ng befor|
|000001a0| 65 20 74 68 69 73 20 6c | 69 6e 65 2c 20 74 68 65 |e this l|ine, the|
|000001b0| 6e 20 75 6e 70 61 63 6b | 0a 23 20 69 74 20 62 79 |n unpack|.# it by|
|000001c0| 20 73 61 76 69 6e 67 20 | 69 74 20 69 6e 74 6f 20 | saving |it into |
|000001d0| 61 20 66 69 6c 65 20 61 | 6e 64 20 74 79 70 69 6e |a file a|nd typin|
|000001e0| 67 20 22 73 68 20 66 69 | 6c 65 22 2e 20 20 54 6f |g "sh fi|le". To|
|000001f0| 20 6f 76 65 72 77 72 69 | 74 65 20 65 78 69 73 74 | overwri|te exist|
|00000200| 69 6e 67 0a 23 20 66 69 | 6c 65 73 2c 20 74 79 70 |ing.# fi|les, typ|
|00000210| 65 20 22 73 68 20 66 69 | 6c 65 20 2d 63 22 2e 20 |e "sh fi|le -c". |
|00000220| 20 59 6f 75 20 63 61 6e | 20 61 6c 73 6f 20 66 65 | You can| also fe|
|00000230| 65 64 20 74 68 69 73 20 | 61 73 20 73 74 61 6e 64 |ed this |as stand|
|00000240| 61 72 64 20 69 6e 70 75 | 74 20 76 69 61 0a 23 20 |ard inpu|t via.# |
|00000250| 75 6e 73 68 61 72 2c 20 | 6f 72 20 62 79 20 74 79 |unshar, |or by ty|
|00000260| 70 69 6e 67 20 22 73 68 | 20 3c 66 69 6c 65 22 2c |ping "sh| <file",|
|00000270| 20 65 2e 67 2e 2e 20 20 | 49 66 20 74 68 69 73 20 | e.g.. |If this |
|00000280| 61 72 63 68 69 76 65 20 | 69 73 20 63 6f 6d 70 6c |archive |is compl|
|00000290| 65 74 65 2c 20 79 6f 75 | 0a 23 20 77 69 6c 6c 20 |ete, you|.# will |
|000002a0| 73 65 65 20 74 68 65 20 | 66 6f 6c 6c 6f 77 69 6e |see the |followin|
|000002b0| 67 20 6d 65 73 73 61 67 | 65 20 61 74 20 74 68 65 |g messag|e at the|
|000002c0| 20 65 6e 64 3a 0a 23 09 | 09 22 45 6e 64 20 6f 66 | end:.#.|."End of|
|000002d0| 20 61 72 63 68 69 76 65 | 20 33 20 28 6f 66 20 35 | archive| 3 (of 5|
|000002e0| 29 2e 22 0a 23 20 43 6f | 6e 74 65 6e 74 73 3a 20 |).".# Co|ntents: |
|000002f0| 20 41 41 4c 69 6e 65 73 | 2f 75 74 61 68 2e 63 20 | AALines|/utah.c |
|00000300| 41 6c 62 65 72 73 2e 63 | 20 42 6f 75 6e 64 53 70 |Albers.c| BoundSp|
|00000310| 68 65 72 65 2e 63 20 44 | 69 73 73 6f 6c 76 65 2e |here.c D|issolve.|
|00000320| 63 0a 23 20 20 20 44 6f | 75 62 6c 65 4c 69 6e 65 |c.# Do|ubleLine|
|00000330| 2e 63 20 47 72 61 70 68 | 69 63 73 47 65 6d 73 2e |.c Graph|icsGems.|
|00000340| 68 20 4c 69 6e 65 45 64 | 67 65 2e 63 20 4d 61 74 |h LineEd|ge.c Mat|
|00000350| 72 69 78 49 6e 76 65 72 | 74 2e 63 0a 23 20 20 20 |rixInver|t.c.# |
|00000360| 50 6f 6c 79 53 63 61 6e | 2f 70 6f 6c 79 5f 63 6c |PolyScan|/poly_cl|
|00000370| 69 70 2e 63 20 50 6f 6c | 79 53 63 61 6e 2f 70 6f |ip.c Pol|yScan/po|
|00000380| 6c 79 5f 73 63 61 6e 2e | 63 20 52 6f 6f 74 73 33 |ly_scan.|c Roots3|
|00000390| 41 6e 64 34 2e 63 0a 23 | 20 57 72 61 70 70 65 64 |And4.c.#| Wrapped|
|000003a0| 20 62 79 20 63 72 61 69 | 67 40 77 65 65 64 65 61 | by crai|g@weedea|
|000003b0| 74 65 72 20 6f 6e 20 46 | 72 69 20 4f 63 74 20 31 |ter on F|ri Oct 1|
|000003c0| 32 20 31 35 3a 35 33 3a | 31 33 20 31 39 39 30 0a |2 15:53:|13 1990.|
|000003d0| 50 41 54 48 3d 2f 62 69 | 6e 3a 2f 75 73 72 2f 62 |PATH=/bi|n:/usr/b|
|000003e0| 69 6e 3a 2f 75 73 72 2f | 75 63 62 20 3b 20 65 78 |in:/usr/|ucb ; ex|
|000003f0| 70 6f 72 74 20 50 41 54 | 48 0a 69 66 20 74 65 73 |port PAT|H.if tes|
|00000400| 74 20 2d 66 20 27 41 41 | 4c 69 6e 65 73 2f 75 74 |t -f 'AA|Lines/ut|
|00000410| 61 68 2e 63 27 20 2d 61 | 20 22 24 7b 31 7d 22 20 |ah.c' -a| "${1}" |
|00000420| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|00000430| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c | echo s|har: Wil|
|00000440| 6c 20 6e 6f 74 20 63 6c | 6f 62 62 65 72 20 65 78 |l not cl|obber ex|
|00000450| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 27 41 |isting f|ile \"'A|
|00000460| 41 4c 69 6e 65 73 2f 75 | 74 61 68 2e 63 27 5c 22 |ALines/u|tah.c'\"|
|00000470| 0a 65 6c 73 65 0a 65 63 | 68 6f 20 73 68 61 72 3a |.else.ec|ho shar:|
|00000480| 20 45 78 74 72 61 63 74 | 69 6e 67 20 5c 22 27 41 | Extract|ing \"'A|
|00000490| 41 4c 69 6e 65 73 2f 75 | 74 61 68 2e 63 27 5c 22 |ALines/u|tah.c'\"|
|000004a0| 20 5c 28 34 35 33 32 20 | 63 68 61 72 61 63 74 65 | \(4532 |characte|
|000004b0| 72 73 5c 29 0a 73 65 64 | 20 22 73 2f 5e 58 2f 2f |rs\).sed| "s/^X//|
|000004c0| 22 20 3e 27 41 41 4c 69 | 6e 65 73 2f 75 74 61 68 |" >'AALi|nes/utah|
|000004d0| 2e 63 27 20 3c 3c 27 45 | 4e 44 5f 4f 46 5f 46 49 |.c' <<'E|ND_OF_FI|
|000004e0| 4c 45 27 0a 58 2f 2a 0a | 58 09 66 69 6c 65 3a 09 |LE'.X/*.|X.file:.|
|000004f0| 09 75 74 61 68 2e 63 0a | 58 09 64 65 73 63 72 69 |.utah.c.|X.descri|
|00000500| 70 74 69 6f 6e 3a 09 69 | 6e 74 65 72 66 61 63 65 |ption:.i|nterface|
|00000510| 20 74 6f 20 55 74 61 68 | 20 52 4c 45 20 74 6f 6f | to Utah| RLE too|
|00000520| 6c 6b 69 74 0a 58 09 61 | 75 74 68 6f 72 3a 09 09 |lkit.X.a|uthor:..|
|00000530| 41 2e 20 54 2e 20 43 61 | 6d 70 62 65 6c 6c 0a 58 |A. T. Ca|mpbell.X|
|00000540| 09 64 61 74 65 3a 09 09 | 4f 63 74 6f 62 65 72 20 |.date:..|October |
|00000550| 32 37 2c 20 31 39 38 39 | 0a 58 2a 2f 0a 58 0a 58 |27, 1989|.X*/.X.X|
|00000560| 23 69 66 6e 64 65 66 20 | 6c 69 6e 74 0a 58 73 74 |#ifndef |lint.Xst|
|00000570| 61 74 69 63 20 63 68 61 | 72 09 73 63 63 73 69 64 |atic cha|r.sccsid|
|00000580| 5b 5d 20 3d 20 22 25 57 | 25 20 25 47 25 22 3b 09 |[] = "%W|% %G%";.|
|00000590| 09 2f 2a 20 53 43 43 53 | 20 69 6e 66 6f 20 2a 2f |./* SCCS| info */|
|000005a0| 0a 58 23 65 6e 64 69 66 | 20 6c 69 6e 74 0a 58 09 |.X#endif| lint.X.|
|000005b0| 0a 58 23 69 6e 63 6c 75 | 64 65 20 3c 6d 61 74 68 |.X#inclu|de <math|
|000005c0| 2e 68 3e 0a 58 23 69 6e | 63 6c 75 64 65 20 3c 73 |.h>.X#in|clude <s|
|000005d0| 74 64 69 6f 2e 68 3e 0a | 58 23 69 66 64 65 66 20 |tdio.h>.|X#ifdef |
|000005e0| 73 65 71 75 65 6e 74 0a | 58 23 69 6e 63 6c 75 64 |sequent.|X#includ|
|000005f0| 65 20 3c 73 74 72 69 6e | 67 73 2e 68 3e 0a 58 23 |e <strin|gs.h>.X#|
|00000600| 65 6c 73 65 0a 58 23 69 | 6e 63 6c 75 64 65 20 3c |else.X#i|nclude <|
|00000610| 73 74 72 69 6e 67 2e 68 | 3e 0a 58 23 65 6e 64 69 |string.h|>.X#endi|
|00000620| 66 0a 58 23 69 6e 63 6c | 75 64 65 20 22 75 74 61 |f.X#incl|ude "uta|
|00000630| 68 2e 68 22 0a 58 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a |h.h".X.X|/*******|
|00000640| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000650| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000660| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000670| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000680| 2a 2a 2a 2a 2a 2a 2a 2f | 0a 58 0a 58 2f 2a 20 72 |*******/|.X.X/* r|
|00000690| 65 74 75 72 6e 20 76 61 | 6c 75 65 73 20 2a 2f 0a |eturn va|lues */.|
|000006a0| 58 65 78 74 65 72 6e 20 | 76 6f 69 64 09 66 72 65 |Xextern |void.fre|
|000006b0| 65 28 29 3b 0a 58 65 78 | 74 65 72 6e 20 63 68 61 |e();.Xex|tern cha|
|000006c0| 72 09 2a 6d 61 6c 6c 6f | 63 28 29 3b 0a 58 0a 58 |r.*mallo|c();.X.X|
|000006d0| 2f 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |/*******|********|
|000006e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000006f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000700| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000710| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2f |********|*******/|
|00000720| 0a 58 0a 58 75 74 61 68 | 5f 72 65 61 64 5f 63 6c |.X.Xutah|_read_cl|
|00000730| 6f 73 65 28 75 66 70 29 | 0a 58 55 54 41 48 5f 46 |ose(ufp)|.XUTAH_F|
|00000740| 49 4c 45 09 2a 75 66 70 | 3b 0a 58 7b 0a 58 09 72 |ILE.*ufp|;.X{.X.r|
|00000750| 65 74 75 72 6e 28 30 29 | 3b 0a 58 7d 0a 58 0a 58 |eturn(0)|;.X}.X.X|
|00000760| 2f 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |/*******|********|
|00000770| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000780| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000790| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000007a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2f |********|*******/|
|000007b0| 0a 58 0a 58 55 54 41 48 | 5f 46 49 4c 45 20 2a 0a |.X.XUTAH|_FILE *.|
|000007c0| 58 75 74 61 68 5f 72 65 | 61 64 5f 69 6e 69 74 28 |Xutah_re|ad_init(|
|000007d0| 66 6e 61 6d 65 2c 20 68 | 74 2c 20 77 64 29 0a 58 |fname, h|t, wd).X|
|000007e0| 0a 58 63 68 61 72 09 2a | 66 6e 61 6d 65 3b 0a 58 |.Xchar.*|fname;.X|
|000007f0| 69 6e 74 09 2a 68 74 2c | 20 2a 77 64 3b 0a 58 7b |int.*ht,| *wd;.X{|
|00000800| 0a 58 09 46 49 4c 45 09 | 09 2a 66 70 3b 0a 58 09 |.X.FILE.|.*fp;.X.|
|00000810| 55 54 41 48 5f 46 49 4c | 45 09 2a 75 66 70 3b 0a |UTAH_FIL|E.*ufp;.|
|00000820| 58 0a 58 09 2f 2a 20 6f | 70 65 6e 20 6f 75 74 70 |X.X./* o|pen outp|
|00000830| 75 74 20 73 74 72 65 61 | 6d 20 2a 2f 0a 58 09 69 |ut strea|m */.X.i|
|00000840| 66 20 28 21 73 74 72 63 | 6d 70 28 66 6e 61 6d 65 |f (!strc|mp(fname|
|00000850| 2c 20 22 22 29 29 0a 58 | 09 09 66 70 20 3d 20 73 |, "")).X|..fp = s|
|00000860| 74 64 69 6e 3b 0a 58 09 | 65 6c 73 65 20 7b 0a 58 |tdin;.X.|else {.X|
|00000870| 09 09 69 66 20 28 28 66 | 70 20 3d 20 66 6f 70 65 |..if ((f|p = fope|
|00000880| 6e 28 66 6e 61 6d 65 2c | 20 22 72 22 29 29 20 3d |n(fname,| "r")) =|
|00000890| 3d 20 4e 55 4c 4c 29 0a | 58 09 09 72 65 74 75 72 |= NULL).|X..retur|
|000008a0| 6e 28 4e 55 4c 4c 29 3b | 0a 58 09 7d 0a 58 0a 58 |n(NULL);|.X.}.X.X|
|000008b0| 20 09 2f 2a 20 63 68 61 | 6e 67 65 20 74 68 65 20 | ./* cha|nge the |
|000008c0| 64 65 66 61 75 6c 74 20 | 73 76 5f 67 6c 6f 62 61 |default |sv_globa|
|000008d0| 6c 73 20 73 74 72 75 63 | 74 20 74 6f 20 6d 61 74 |ls struc|t to mat|
|000008e0| 63 68 20 77 68 61 74 20 | 77 65 20 6e 65 65 64 20 |ch what |we need |
|000008f0| 2a 2f 0a 58 09 75 66 70 | 20 3d 20 28 55 54 41 48 |*/.X.ufp| = (UTAH|
|00000900| 5f 46 49 4c 45 20 2a 29 | 20 6d 61 6c 6c 6f 63 28 |_FILE *)| malloc(|
|00000910| 73 69 7a 65 6f 66 28 55 | 54 41 48 5f 46 49 4c 45 |sizeof(U|TAH_FILE|
|00000920| 29 29 3b 0a 58 09 2a 75 | 66 70 20 3d 20 73 76 5f |));.X.*u|fp = sv_|
|00000930| 67 6c 6f 62 61 6c 73 3b | 0a 58 09 75 66 70 2d 3e |globals;|.X.ufp->|
|00000940| 73 76 66 62 5f 66 64 20 | 3d 20 66 70 3b 0a 58 0a |svfb_fd |= fp;.X.|
|00000950| 58 09 2f 2a 20 72 65 61 | 64 20 74 68 65 20 68 65 |X./* rea|d the he|
|00000960| 61 64 65 72 20 69 6e 20 | 74 68 65 20 69 6e 70 75 |ader in |the inpu|
|00000970| 74 20 66 69 6c 65 20 2a | 2f 0a 58 20 20 09 69 66 |t file *|/.X .if|
|00000980| 20 28 72 6c 65 5f 67 65 | 74 5f 73 65 74 75 70 28 | (rle_ge|t_setup(|
|00000990| 75 66 70 29 20 21 3d 20 | 30 29 0a 58 09 09 72 65 |ufp) != |0).X..re|
|000009a0| 74 75 72 6e 28 4e 55 4c | 4c 29 3b 0a 58 0a 58 09 |turn(NUL|L);.X.X.|
|000009b0| 2f 2a 20 67 65 74 20 69 | 6d 61 67 65 20 73 69 7a |/* get i|mage siz|
|000009c0| 65 20 2a 2f 0a 58 09 2a | 77 64 20 3d 20 75 66 70 |e */.X.*|wd = ufp|
|000009d0| 2d 3e 73 76 5f 78 6d 61 | 78 20 2d 20 75 66 70 2d |->sv_xma|x - ufp-|
|000009e0| 3e 73 76 5f 78 6d 69 6e | 20 2b 20 31 3b 0a 58 09 |>sv_xmin| + 1;.X.|
|000009f0| 2a 68 74 20 3d 20 75 66 | 70 2d 3e 73 76 5f 79 6d |*ht = uf|p->sv_ym|
|00000a00| 61 78 20 2d 20 75 66 70 | 2d 3e 73 76 5f 79 6d 69 |ax - ufp|->sv_ymi|
|00000a10| 6e 20 2b 20 31 3b 0a 58 | 0a 58 09 2f 2a 20 6e 6f |n + 1;.X|.X./* no|
|00000a20| 72 6d 61 6c 20 74 65 72 | 6d 69 6e 61 74 69 6f 6e |rmal ter|mination|
|00000a30| 20 2a 2f 0a 58 09 72 65 | 74 75 72 6e 28 75 66 70 | */.X.re|turn(ufp|
|00000a40| 29 3b 0a 58 7d 0a 58 0a | 58 2f 2a 2a 2a 2a 2a 2a |);.X}.X.|X/******|
|00000a50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000a60| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000a70| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000a80| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000a90| 2a 2a 2a 2a 2a 2a 2a 2a | 2f 0a 58 0a 58 75 74 61 |********|/.X.Xuta|
|00000aa0| 68 5f 72 65 61 64 5f 70 | 69 78 65 6c 73 28 75 66 |h_read_p|ixels(uf|
|00000ab0| 70 2c 20 70 69 78 65 6c | 73 29 0a 58 0a 58 55 54 |p, pixel|s).X.XUT|
|00000ac0| 41 48 5f 46 49 4c 45 20 | 09 2a 75 66 70 3b 0a 58 |AH_FILE |.*ufp;.X|
|00000ad0| 75 6e 73 69 67 6e 65 64 | 20 63 68 61 72 09 70 69 |unsigned| char.pi|
|00000ae0| 78 65 6c 73 5b 5d 5b 33 | 5d 3b 0a 58 7b 0a 58 09 |xels[][3|];.X{.X.|
|00000af0| 73 74 61 74 69 63 20 75 | 6e 73 69 67 6e 65 64 09 |static u|nsigned.|
|00000b00| 6e 20 3d 20 30 3b 0a 58 | 09 73 74 61 74 69 63 20 |n = 0;.X|.static |
|00000b10| 75 6e 73 69 67 6e 65 64 | 20 63 68 61 72 09 2a 72 |unsigned| char.*r|
|00000b20| 20 3d 20 4e 55 4c 4c 2c | 20 2a 67 20 3d 20 4e 55 | = NULL,| *g = NU|
|00000b30| 4c 4c 2c 20 2a 62 20 3d | 20 4e 55 4c 4c 3b 0a 58 |LL, *b =| NULL;.X|
|00000b40| 09 69 6e 74 09 09 69 2c | 20 77 69 64 74 68 3b 0a |.int..i,| width;.|
|00000b50| 58 0a 58 09 2f 2a 20 61 | 6c 6c 6f 63 61 74 65 20 |X.X./* a|llocate |
|00000b60| 73 74 6f 72 61 67 65 20 | 2a 2f 0a 58 09 77 69 64 |storage |*/.X.wid|
|00000b70| 74 68 20 3d 20 75 66 70 | 2d 3e 73 76 5f 78 6d 61 |th = ufp|->sv_xma|
|00000b80| 78 20 2b 20 31 3b 0a 58 | 09 69 66 20 28 77 69 64 |x + 1;.X|.if (wid|
|00000b90| 74 68 20 3e 20 6e 29 20 | 7b 0a 58 09 09 69 66 20 |th > n) |{.X..if |
|00000ba0| 28 6e 20 3e 20 30 29 20 | 7b 0a 58 09 09 09 66 72 |(n > 0) |{.X...fr|
|00000bb0| 65 65 28 28 63 68 61 72 | 20 2a 29 72 29 3b 0a 58 |ee((char| *)r);.X|
|00000bc0| 09 09 09 66 72 65 65 28 | 28 63 68 61 72 20 2a 29 |...free(|(char *)|
|00000bd0| 67 29 3b 0a 58 09 09 09 | 66 72 65 65 28 28 63 68 |g);.X...|free((ch|
|00000be0| 61 72 20 2a 29 62 29 3b | 0a 58 09 09 7d 0a 58 09 |ar *)b);|.X..}.X.|
|00000bf0| 09 6e 20 3d 20 77 69 64 | 74 68 3b 0a 58 09 09 72 |.n = wid|th;.X..r|
|00000c00| 20 3d 20 28 75 6e 73 69 | 67 6e 65 64 20 63 68 61 | = (unsi|gned cha|
|00000c10| 72 20 2a 29 20 6d 61 6c | 6c 6f 63 28 6e 20 2a 20 |r *) mal|loc(n * |
|00000c20| 73 69 7a 65 6f 66 28 75 | 6e 73 69 67 6e 65 64 20 |sizeof(u|nsigned |
|00000c30| 63 68 61 72 29 29 3b 0a | 58 09 09 67 20 3d 20 28 |char));.|X..g = (|
|00000c40| 75 6e 73 69 67 6e 65 64 | 20 63 68 61 72 20 2a 29 |unsigned| char *)|
|00000c50| 20 6d 61 6c 6c 6f 63 28 | 6e 20 2a 20 73 69 7a 65 | malloc(|n * size|
|00000c60| 6f 66 28 75 6e 73 69 67 | 6e 65 64 20 63 68 61 72 |of(unsig|ned char|
|00000c70| 29 29 3b 0a 58 09 09 62 | 20 3d 20 28 75 6e 73 69 |));.X..b| = (unsi|
|00000c80| 67 6e 65 64 20 63 68 61 | 72 20 2a 29 20 6d 61 6c |gned cha|r *) mal|
|00000c90| 6c 6f 63 28 6e 20 2a 20 | 73 69 7a 65 6f 66 28 75 |loc(n * |sizeof(u|
|00000ca0| 6e 73 69 67 6e 65 64 20 | 63 68 61 72 29 29 3b 0a |nsigned |char));.|
|00000cb0| 58 09 7d 0a 58 0a 58 09 | 2f 2a 20 72 65 61 64 20 |X.}.X.X.|/* read |
|00000cc0| 74 68 69 73 20 72 6f 77 | 20 2a 2f 0a 58 09 75 74 |this row| */.X.ut|
|00000cd0| 61 68 5f 72 65 61 64 5f | 72 67 62 28 75 66 70 2c |ah_read_|rgb(ufp,|
|00000ce0| 20 72 2c 20 67 2c 20 62 | 29 3b 0a 58 0a 58 09 2f | r, g, b|);.X.X./|
|00000cf0| 2a 20 63 6f 6e 76 65 72 | 74 20 74 6f 20 70 69 78 |* conver|t to pix|
|00000d00| 65 6c 73 20 2a 2f 0a 58 | 09 66 6f 72 20 28 69 20 |els */.X|.for (i |
|00000d10| 3d 20 30 3b 20 69 20 3c | 20 77 69 64 74 68 3b 20 |= 0; i <| width; |
|00000d20| 69 2b 2b 29 20 7b 0a 58 | 09 09 70 69 78 65 6c 73 |i++) {.X|..pixels|
|00000d30| 5b 69 5d 5b 30 5d 20 3d | 20 72 5b 69 5d 3b 0a 58 |[i][0] =| r[i];.X|
|00000d40| 09 09 70 69 78 65 6c 73 | 5b 69 5d 5b 31 5d 20 3d |..pixels|[i][1] =|
|00000d50| 20 67 5b 69 5d 3b 0a 58 | 09 09 70 69 78 65 6c 73 | g[i];.X|..pixels|
|00000d60| 5b 69 5d 5b 32 5d 20 3d | 20 62 5b 69 5d 3b 0a 58 |[i][2] =| b[i];.X|
|00000d70| 09 7d 0a 58 0a 58 09 72 | 65 74 75 72 6e 28 30 29 |.}.X.X.r|eturn(0)|
|00000d80| 3b 0a 58 7d 0a 58 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a |;.X}.X.X|/*******|
|00000d90| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000da0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000db0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000dc0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000dd0| 2a 2a 2a 2a 2a 2a 2a 2f | 0a 58 0a 58 75 74 61 68 |*******/|.X.Xutah|
|00000de0| 5f 72 65 61 64 5f 72 67 | 62 28 75 66 70 2c 20 72 |_read_rg|b(ufp, r|
|00000df0| 2c 20 67 2c 20 62 29 0a | 58 0a 58 55 54 41 48 5f |, g, b).|X.XUTAH_|
|00000e00| 46 49 4c 45 09 2a 75 66 | 70 3b 0a 58 75 6e 73 69 |FILE.*uf|p;.Xunsi|
|00000e10| 67 6e 65 64 20 63 68 61 | 72 09 72 5b 5d 2c 20 67 |gned cha|r.r[], g|
|00000e20| 5b 5d 2c 20 62 5b 5d 3b | 0a 58 7b 0a 58 09 72 6c |[], b[];|.X{.X.rl|
|00000e30| 65 5f 70 69 78 65 6c 09 | 2a 72 6f 77 73 5b 33 5d |e_pixel.|*rows[3]|
|00000e40| 3b 0a 58 0a 58 09 2f 2a | 20 73 65 74 20 63 6f 6c |;.X.X./*| set col|
|00000e50| 6f 72 20 63 68 61 6e 6e | 65 6c 73 20 2a 2f 0a 58 |or chann|els */.X|
|00000e60| 09 72 6f 77 73 5b 30 5d | 20 3d 20 72 3b 0a 58 09 |.rows[0]| = r;.X.|
|00000e70| 72 6f 77 73 5b 31 5d 20 | 3d 20 67 3b 0a 58 09 72 |rows[1] |= g;.X.r|
|00000e80| 6f 77 73 5b 32 5d 20 3d | 20 62 3b 0a 58 0a 58 09 |ows[2] =| b;.X.X.|
|00000e90| 2f 2a 20 72 65 61 64 20 | 74 68 69 73 20 72 6f 77 |/* read |this row|
|00000ea0| 20 2a 2f 0a 58 09 72 6c | 65 5f 67 65 74 72 6f 77 | */.X.rl|e_getrow|
|00000eb0| 28 75 66 70 2c 20 72 6f | 77 73 29 3b 0a 58 09 72 |(ufp, ro|ws);.X.r|
|00000ec0| 65 74 75 72 6e 28 30 29 | 3b 0a 58 7d 0a 58 0a 58 |eturn(0)|;.X}.X.X|
|00000ed0| 2f 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |/*******|********|
|00000ee0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000ef0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000f00| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000f10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2f |********|*******/|
|00000f20| 0a 58 0a 58 75 74 61 68 | 5f 77 72 69 74 65 5f 63 |.X.Xutah|_write_c|
|00000f30| 6c 6f 73 65 28 75 66 70 | 29 0a 58 0a 58 55 54 41 |lose(ufp|).X.XUTA|
|00000f40| 48 5f 46 49 4c 45 09 2a | 75 66 70 3b 0a 58 7b 0a |H_FILE.*|ufp;.X{.|
|00000f50| 58 09 69 66 20 28 21 75 | 66 70 29 20 72 65 74 75 |X.if (!u|fp) retu|
|00000f60| 72 6e 28 2d 31 29 3b 0a | 58 09 73 76 5f 70 75 74 |rn(-1);.|X.sv_put|
|00000f70| 65 6f 66 28 75 66 70 29 | 3b 0a 58 09 72 65 74 75 |eof(ufp)|;.X.retu|
|00000f80| 72 6e 28 30 29 3b 0a 58 | 7d 0a 58 0a 58 2f 2a 2a |rn(0);.X|}.X.X/**|
|00000f90| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000fa0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000fb0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000fc0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000fd0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2f 0a 58 0a |********|****/.X.|
|00000fe0| 58 55 54 41 48 5f 46 49 | 4c 45 20 2a 0a 58 75 74 |XUTAH_FI|LE *.Xut|
|00000ff0| 61 68 5f 77 72 69 74 65 | 5f 69 6e 69 74 28 66 6e |ah_write|_init(fn|
|00001000| 61 6d 65 2c 20 68 74 2c | 20 77 64 29 0a 58 0a 58 |ame, ht,| wd).X.X|
|00001010| 63 68 61 72 09 2a 66 6e | 61 6d 65 3b 0a 58 69 6e |char.*fn|ame;.Xin|
|00001020| 74 09 68 74 2c 20 77 64 | 3b 0a 58 7b 0a 58 09 46 |t.ht, wd|;.X{.X.F|
|00001030| 49 4c 45 09 09 2a 66 70 | 3b 0a 58 09 55 54 41 48 |ILE..*fp|;.X.UTAH|
|00001040| 5f 46 49 4c 45 09 2a 75 | 66 70 3b 0a 58 0a 58 09 |_FILE.*u|fp;.X.X.|
|00001050| 2f 2a 20 6f 70 65 6e 20 | 6f 75 74 70 75 74 20 73 |/* open |output s|
|00001060| 74 72 65 61 6d 20 2a 2f | 0a 58 09 69 66 20 28 21 |tream */|.X.if (!|
|00001070| 73 74 72 63 6d 70 28 66 | 6e 61 6d 65 2c 20 22 22 |strcmp(f|name, ""|
|00001080| 29 29 0a 58 09 09 66 70 | 20 3d 20 73 74 64 6f 75 |)).X..fp| = stdou|
|00001090| 74 3b 0a 58 09 65 6c 73 | 65 20 7b 0a 58 09 09 69 |t;.X.els|e {.X..i|
|000010a0| 66 20 28 28 66 70 20 3d | 20 66 6f 70 65 6e 28 66 |f ((fp =| fopen(f|
|000010b0| 6e 61 6d 65 2c 20 22 77 | 22 29 29 20 3d 3d 20 4e |name, "w|")) == N|
|000010c0| 55 4c 4c 29 0a 58 09 09 | 72 65 74 75 72 6e 28 4e |ULL).X..|return(N|
|000010d0| 55 4c 4c 29 3b 0a 58 09 | 7d 0a 58 0a 58 20 09 2f |ULL);.X.|}.X.X ./|
|000010e0| 2a 20 63 68 61 6e 67 65 | 20 74 68 65 20 64 65 66 |* change| the def|
|000010f0| 61 75 6c 74 20 73 76 5f | 67 6c 6f 62 61 6c 73 20 |ault sv_|globals |
|00001100| 73 74 72 75 63 74 20 74 | 6f 20 6d 61 74 63 68 20 |struct t|o match |
|00001110| 77 68 61 74 20 77 65 20 | 6e 65 65 64 20 2a 2f 0a |what we |need */.|
|00001120| 58 09 75 66 70 20 3d 20 | 28 55 54 41 48 5f 46 49 |X.ufp = |(UTAH_FI|
|00001130| 4c 45 20 2a 29 20 6d 61 | 6c 6c 6f 63 28 73 69 7a |LE *) ma|lloc(siz|
|00001140| 65 6f 66 28 55 54 41 48 | 5f 46 49 4c 45 29 29 3b |eof(UTAH|_FILE));|
|00001150| 0a 58 09 2a 75 66 70 20 | 3d 20 73 76 5f 67 6c 6f |.X.*ufp |= sv_glo|
|00001160| 62 61 6c 73 3b 0a 58 09 | 75 66 70 2d 3e 73 76 66 |bals;.X.|ufp->svf|
|00001170| 62 5f 66 64 20 3d 20 66 | 70 3b 0a 58 09 75 66 70 |b_fd = f|p;.X.ufp|
|00001180| 2d 3e 73 76 5f 78 6d 61 | 78 20 3d 20 77 64 20 2d |->sv_xma|x = wd -|
|00001190| 20 31 3b 0a 58 09 75 66 | 70 2d 3e 73 76 5f 79 6d | 1;.X.uf|p->sv_ym|
|000011a0| 61 78 20 3d 20 68 74 20 | 2d 20 31 3b 0a 58 09 75 |ax = ht |- 1;.X.u|
|000011b0| 66 70 2d 3e 73 76 5f 61 | 6c 70 68 61 20 3d 20 30 |fp->sv_a|lpha = 0|
|000011c0| 3b 09 2f 2a 20 4e 6f 20 | 63 6f 76 65 72 61 67 65 |;./* No |coverage|
|000011d0| 20 28 61 6c 70 68 61 29 | 20 2a 2f 0a 58 0a 58 09 | (alpha)| */.X.X.|
|000011e0| 2f 2a 20 63 72 65 61 74 | 65 20 74 68 65 20 68 65 |/* creat|e the he|
|000011f0| 61 64 65 72 20 69 6e 20 | 74 68 65 20 6f 75 74 70 |ader in |the outp|
|00001200| 75 74 20 66 69 6c 65 20 | 2a 2f 0a 58 20 20 09 73 |ut file |*/.X .s|
|00001210| 76 5f 73 65 74 75 70 28 | 52 55 4e 5f 44 49 53 50 |v_setup(|RUN_DISP|
|00001220| 41 54 43 48 2c 20 75 66 | 70 29 3b 0a 58 0a 58 09 |ATCH, uf|p);.X.X.|
|00001230| 2f 2a 20 6e 6f 72 6d 61 | 6c 20 74 65 72 6d 69 6e |/* norma|l termin|
|00001240| 61 74 69 6f 6e 20 2a 2f | 0a 58 09 72 65 74 75 72 |ation */|.X.retur|
|00001250| 6e 28 75 66 70 29 3b 0a | 58 7d 0a 58 0a 58 2f 2a |n(ufp);.|X}.X.X/*|
|00001260| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001270| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001280| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001290| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000012a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2f 0a 58 |********|*****/.X|
|000012b0| 0a 58 75 74 61 68 5f 77 | 72 69 74 65 5f 70 69 78 |.Xutah_w|rite_pix|
|000012c0| 65 6c 73 28 75 66 70 2c | 20 70 69 78 65 6c 73 29 |els(ufp,| pixels)|
|000012d0| 0a 58 0a 58 55 54 41 48 | 5f 46 49 4c 45 09 2a 75 |.X.XUTAH|_FILE.*u|
|000012e0| 66 70 3b 0a 58 75 6e 73 | 69 67 6e 65 64 20 63 68 |fp;.Xuns|igned ch|
|000012f0| 61 72 09 70 69 78 65 6c | 73 5b 5d 5b 33 5d 3b 0a |ar.pixel|s[][3];.|
|00001300| 58 7b 0a 58 09 73 74 61 | 74 69 63 20 75 6e 73 69 |X{.X.sta|tic unsi|
|00001310| 67 6e 65 64 09 6e 20 3d | 20 30 3b 0a 58 09 73 74 |gned.n =| 0;.X.st|
|00001320| 61 74 69 63 20 75 6e 73 | 69 67 6e 65 64 20 63 68 |atic uns|igned ch|
|00001330| 61 72 09 2a 72 20 3d 20 | 4e 55 4c 4c 2c 20 2a 67 |ar.*r = |NULL, *g|
|00001340| 20 3d 20 4e 55 4c 4c 2c | 20 2a 62 20 3d 20 4e 55 | = NULL,| *b = NU|
|00001350| 4c 4c 3b 0a 58 09 69 6e | 74 09 09 69 2c 20 77 69 |LL;.X.in|t..i, wi|
|00001360| 64 74 68 3b 0a 58 0a 58 | 09 2f 2a 20 61 6c 6c 6f |dth;.X.X|./* allo|
|00001370| 63 61 74 65 20 73 74 6f | 72 61 67 65 20 2a 2f 0a |cate sto|rage */.|
|00001380| 58 09 77 69 64 74 68 20 | 3d 20 75 66 70 2d 3e 73 |X.width |= ufp->s|
|00001390| 76 5f 78 6d 61 78 20 2b | 20 31 3b 0a 58 09 69 66 |v_xmax +| 1;.X.if|
|000013a0| 20 28 77 69 64 74 68 20 | 3e 20 6e 29 20 7b 0a 58 | (width |> n) {.X|
|000013b0| 09 09 69 66 20 28 6e 20 | 3e 20 30 29 20 7b 0a 58 |..if (n |> 0) {.X|
|000013c0| 09 09 09 66 72 65 65 28 | 28 63 68 61 72 20 2a 29 |...free(|(char *)|
|000013d0| 72 29 3b 0a 58 09 09 09 | 66 72 65 65 28 28 63 68 |r);.X...|free((ch|
|000013e0| 61 72 20 2a 29 67 29 3b | 0a 58 09 09 09 66 72 65 |ar *)g);|.X...fre|
|000013f0| 65 28 28 63 68 61 72 20 | 2a 29 62 29 3b 0a 58 09 |e((char |*)b);.X.|
|00001400| 09 7d 0a 58 09 09 6e 20 | 3d 20 77 69 64 74 68 3b |.}.X..n |= width;|
|00001410| 0a 58 09 09 72 20 3d 20 | 28 75 6e 73 69 67 6e 65 |.X..r = |(unsigne|
|00001420| 64 20 63 68 61 72 20 2a | 29 20 6d 61 6c 6c 6f 63 |d char *|) malloc|
|00001430| 28 6e 20 2a 20 73 69 7a | 65 6f 66 28 75 6e 73 69 |(n * siz|eof(unsi|
|00001440| 67 6e 65 64 20 63 68 61 | 72 29 29 3b 0a 58 09 09 |gned cha|r));.X..|
|00001450| 67 20 3d 20 28 75 6e 73 | 69 67 6e 65 64 20 63 68 |g = (uns|igned ch|
|00001460| 61 72 20 2a 29 20 6d 61 | 6c 6c 6f 63 28 6e 20 2a |ar *) ma|lloc(n *|
|00001470| 20 73 69 7a 65 6f 66 28 | 75 6e 73 69 67 6e 65 64 | sizeof(|unsigned|
|00001480| 20 63 68 61 72 29 29 3b | 0a 58 09 09 62 20 3d 20 | char));|.X..b = |
|00001490| 28 75 6e 73 69 67 6e 65 | 64 20 63 68 61 72 20 2a |(unsigne|d char *|
|000014a0| 29 20 6d 61 6c 6c 6f 63 | 28 6e 20 2a 20 73 69 7a |) malloc|(n * siz|
|000014b0| 65 6f 66 28 75 6e 73 69 | 67 6e 65 64 20 63 68 61 |eof(unsi|gned cha|
|000014c0| 72 29 29 3b 0a 58 09 7d | 0a 58 0a 58 09 2f 2a 20 |r));.X.}|.X.X./* |
|000014d0| 63 6f 6e 76 65 72 74 20 | 74 6f 20 63 6f 6c 6f 72 |convert |to color|
|000014e0| 20 63 68 61 6e 6e 65 6c | 73 20 2a 2f 0a 58 09 66 | channel|s */.X.f|
|000014f0| 6f 72 20 28 69 20 3d 20 | 30 3b 20 69 20 3c 20 77 |or (i = |0; i < w|
|00001500| 69 64 74 68 3b 20 69 2b | 2b 29 20 7b 0a 58 09 09 |idth; i+|+) {.X..|
|00001510| 72 5b 69 5d 20 3d 20 70 | 69 78 65 6c 73 5b 69 5d |r[i] = p|ixels[i]|
|00001520| 5b 30 5d 3b 0a 58 09 09 | 67 5b 69 5d 20 3d 20 70 |[0];.X..|g[i] = p|
|00001530| 69 78 65 6c 73 5b 69 5d | 5b 31 5d 3b 0a 58 09 09 |ixels[i]|[1];.X..|
|00001540| 62 5b 69 5d 20 3d 20 70 | 69 78 65 6c 73 5b 69 5d |b[i] = p|ixels[i]|
|00001550| 5b 32 5d 3b 0a 58 09 7d | 0a 58 0a 58 09 2f 2a 20 |[2];.X.}|.X.X./* |
|00001560| 77 72 69 74 65 20 74 68 | 69 73 20 72 6f 77 20 2a |write th|is row *|
|00001570| 2f 0a 58 09 75 74 61 68 | 5f 77 72 69 74 65 5f 72 |/.X.utah|_write_r|
|00001580| 67 62 28 75 66 70 2c 20 | 72 2c 20 67 2c 20 62 29 |gb(ufp, |r, g, b)|
|00001590| 3b 0a 58 09 72 65 74 75 | 72 6e 28 30 29 3b 0a 58 |;.X.retu|rn(0);.X|
|000015a0| 7d 0a 58 0a 58 2f 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |}.X.X/**|********|
|000015b0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000015c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000015d0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000015e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000015f0| 2a 2a 2a 2a 2f 0a 58 0a | 58 75 74 61 68 5f 77 72 |****/.X.|Xutah_wr|
|00001600| 69 74 65 5f 72 67 62 28 | 75 66 70 2c 20 72 2c 20 |ite_rgb(|ufp, r, |
|00001610| 67 2c 20 62 29 0a 58 0a | 58 55 54 41 48 5f 46 49 |g, b).X.|XUTAH_FI|
|00001620| 4c 45 09 2a 75 66 70 3b | 0a 58 75 6e 73 69 67 6e |LE.*ufp;|.Xunsign|
|00001630| 65 64 20 63 68 61 72 09 | 72 5b 5d 2c 20 67 5b 5d |ed char.|r[], g[]|
|00001640| 2c 20 62 5b 5d 3b 0a 58 | 7b 0a 58 09 72 6c 65 5f |, b[];.X|{.X.rle_|
|00001650| 70 69 78 65 6c 09 2a 72 | 6f 77 73 5b 33 5d 3b 0a |pixel.*r|ows[3];.|
|00001660| 58 09 69 6e 74 09 09 77 | 69 64 74 68 3b 0a 58 0a |X.int..w|idth;.X.|
|00001670| 58 09 2f 2a 20 73 65 74 | 20 63 6f 6c 6f 72 20 63 |X./* set| color c|
|00001680| 68 61 6e 6e 65 6c 73 20 | 2a 2f 0a 58 09 72 6f 77 |hannels |*/.X.row|
|00001690| 73 5b 30 5d 20 3d 20 72 | 3b 0a 58 09 72 6f 77 73 |s[0] = r|;.X.rows|
|000016a0| 5b 31 5d 20 3d 20 67 3b | 0a 58 09 72 6f 77 73 5b |[1] = g;|.X.rows[|
|000016b0| 32 5d 20 3d 20 62 3b 0a | 58 0a 58 09 2f 2a 20 77 |2] = b;.|X.X./* w|
|000016c0| 72 69 74 65 20 74 68 69 | 73 20 72 6f 77 20 2a 2f |rite thi|s row */|
|000016d0| 0a 58 09 77 69 64 74 68 | 20 3d 20 75 66 70 2d 3e |.X.width| = ufp->|
|000016e0| 73 76 5f 78 6d 61 78 20 | 2d 20 75 66 70 2d 3e 73 |sv_xmax |- ufp->s|
|000016f0| 76 5f 78 6d 69 6e 20 2b | 20 31 3b 0a 58 09 73 76 |v_xmin +| 1;.X.sv|
|00001700| 5f 70 75 74 72 6f 77 28 | 72 6f 77 73 2c 20 77 69 |_putrow(|rows, wi|
|00001710| 64 74 68 2c 20 75 66 70 | 29 3b 0a 58 09 72 65 74 |dth, ufp|);.X.ret|
|00001720| 75 72 6e 28 30 29 3b 0a | 58 7d 0a 58 0a 58 2f 2a |urn(0);.|X}.X.X/*|
|00001730| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001740| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001750| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001760| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001770| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2f 0a 45 |********|*****/.E|
|00001780| 4e 44 5f 4f 46 5f 46 49 | 4c 45 0a 69 66 20 74 65 |ND_OF_FI|LE.if te|
|00001790| 73 74 20 34 35 33 32 20 | 2d 6e 65 20 60 77 63 20 |st 4532 |-ne `wc |
|000017a0| 2d 63 20 3c 27 41 41 4c | 69 6e 65 73 2f 75 74 61 |-c <'AAL|ines/uta|
|000017b0| 68 2e 63 27 60 3b 20 74 | 68 65 6e 0a 20 20 20 20 |h.c'`; t|hen. |
|000017c0| 65 63 68 6f 20 73 68 61 | 72 3a 20 5c 22 27 41 41 |echo sha|r: \"'AA|
|000017d0| 4c 69 6e 65 73 2f 75 74 | 61 68 2e 63 27 5c 22 20 |Lines/ut|ah.c'\" |
|000017e0| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|000017f0| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|00001800| 6e 64 20 6f 66 20 27 41 | 41 4c 69 6e 65 73 2f 75 |nd of 'A|ALines/u|
|00001810| 74 61 68 2e 63 27 0a 66 | 69 0a 69 66 20 74 65 73 |tah.c'.f|i.if tes|
|00001820| 74 20 2d 66 20 27 41 6c | 62 65 72 73 2e 63 27 20 |t -f 'Al|bers.c' |
|00001830| 2d 61 20 22 24 7b 31 7d | 22 20 21 3d 20 22 2d 63 |-a "${1}|" != "-c|
|00001840| 22 20 3b 20 74 68 65 6e | 20 0a 20 20 65 63 68 6f |" ; then| . echo|
|00001850| 20 73 68 61 72 3a 20 57 | 69 6c 6c 20 6e 6f 74 20 | shar: W|ill not |
|00001860| 63 6c 6f 62 62 65 72 20 | 65 78 69 73 74 69 6e 67 |clobber |existing|
|00001870| 20 66 69 6c 65 20 5c 22 | 27 41 6c 62 65 72 73 2e | file \"|'Albers.|
|00001880| 63 27 5c 22 0a 65 6c 73 | 65 0a 65 63 68 6f 20 73 |c'\".els|e.echo s|
|00001890| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|000018a0| 5c 22 27 41 6c 62 65 72 | 73 2e 63 27 5c 22 20 5c |\"'Alber|s.c'\" \|
|000018b0| 28 33 39 39 34 20 63 68 | 61 72 61 63 74 65 72 73 |(3994 ch|aracters|
|000018c0| 5c 29 0a 73 65 64 20 22 | 73 2f 5e 58 2f 2f 22 20 |\).sed "|s/^X//" |
|000018d0| 3e 27 41 6c 62 65 72 73 | 2e 63 27 20 3c 3c 27 45 |>'Albers|.c' <<'E|
|000018e0| 4e 44 5f 4f 46 5f 46 49 | 4c 45 27 0a 58 2f 2a 0a |ND_OF_FI|LE'.X/*.|
|000018f0| 58 41 6c 62 65 72 73 20 | 45 71 75 61 6c 2d 41 72 |XAlbers |Equal-Ar|
|00001900| 65 61 20 43 6f 6e 69 63 | 20 4d 61 70 20 50 72 6f |ea Conic| Map Pro|
|00001910| 6a 65 63 74 69 6f 6e 0a | 58 62 79 20 50 61 75 6c |jection.|Xby Paul|
|00001920| 20 42 61 6d 65 0a 58 66 | 72 6f 6d 20 22 47 72 61 | Bame.Xf|rom "Gra|
|00001930| 70 68 69 63 73 20 47 65 | 6d 73 22 2c 20 41 63 61 |phics Ge|ms", Aca|
|00001940| 64 65 6d 69 63 20 50 72 | 65 73 73 2c 20 31 39 39 |demic Pr|ess, 199|
|00001950| 30 0a 58 2a 2f 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |0.X*/.X.|X/*.X * |
|00001960| 41 6c 62 65 72 73 20 43 | 6f 6e 69 63 20 45 71 75 |Albers C|onic Equ|
|00001970| 61 6c 2d 41 72 65 61 20 | 50 72 6f 6a 65 63 74 69 |al-Area |Projecti|
|00001980| 6f 6e 0a 58 20 2a 20 46 | 6f 72 6d 75 6c 61 65 20 |on.X * F|ormulae |
|00001990| 74 61 6b 65 6e 20 66 72 | 6f 6d 20 22 4d 61 70 20 |taken fr|om "Map |
|000019a0| 50 72 6f 6a 65 63 74 69 | 6f 6e 73 20 55 73 65 64 |Projecti|ons Used|
|000019b0| 20 62 79 20 74 68 65 20 | 55 2e 53 2e 20 0a 58 20 | by the |U.S. .X |
|000019c0| 2a 20 47 65 6f 6c 6f 67 | 69 63 61 6c 20 53 75 72 |* Geolog|ical Sur|
|000019d0| 76 65 79 22 20 42 75 6c | 6c 65 74 69 6e 20 23 31 |vey" Bul|letin #1|
|000019e0| 35 33 32 0a 58 20 2a 0a | 58 20 2a 20 45 71 75 61 |532.X *.|X * Equa|
|000019f0| 74 69 6f 6e 20 72 65 66 | 65 72 65 6e 63 65 20 6e |tion ref|erence n|
|00001a00| 75 6d 62 65 72 73 20 61 | 6e 64 20 73 6f 6d 65 20 |umbers a|nd some |
|00001a10| 76 61 72 69 61 62 6c 65 | 20 6e 61 6d 65 73 20 74 |variable| names t|
|00001a20| 61 6b 65 6e 0a 58 20 2a | 20 66 72 6f 6d 20 74 68 |aken.X *| from th|
|00001a30| 65 20 72 65 66 65 72 65 | 6e 63 65 2e 0a 58 20 2a |e refere|nce..X *|
|00001a40| 20 54 6f 20 75 73 65 2c | 20 63 61 6c 6c 20 61 6c | To use,| call al|
|00001a50| 62 65 72 73 20 73 65 74 | 75 70 28 29 20 6f 6e 63 |bers set|up() onc|
|00001a60| 65 20 61 6e 64 20 74 68 | 65 6e 20 61 6c 62 65 72 |e and th|en alber|
|00001a70| 73 5f 70 72 6f 6a 65 63 | 74 28 29 0a 58 20 2a 20 |s_projec|t().X * |
|00001a80| 66 6f 72 20 65 61 63 68 | 20 63 6f 6f 72 64 69 6e |for each| coordin|
|00001a90| 61 74 65 20 70 61 69 72 | 20 6f 66 20 69 6e 74 65 |ate pair| of inte|
|00001aa0| 72 65 73 74 2e 0a 58 2a | 2f 0a 58 0a 58 23 69 6e |rest..X*|/.X.X#in|
|00001ab0| 63 6c 75 64 65 20 22 47 | 72 61 70 68 69 63 73 47 |clude "G|raphicsG|
|00001ac0| 65 6d 73 2e 68 22 09 09 | 0a 58 23 69 6e 63 6c 75 |ems.h"..|.X#inclu|
|00001ad0| 64 65 20 3c 73 74 64 69 | 6f 2e 68 3e 0a 58 23 69 |de <stdi|o.h>.X#i|
|00001ae0| 6e 63 6c 75 64 65 20 3c | 6d 61 74 68 2e 68 3e 0a |nclude <|math.h>.|
|00001af0| 58 0a 58 2f 2a 0a 58 20 | 2a 20 54 68 69 73 20 69 |X.X/*.X |* This i|
|00001b00| 73 20 74 68 65 20 43 6c | 61 72 6b 65 20 31 38 36 |s the Cl|arke 186|
|00001b10| 36 20 45 61 72 74 68 20 | 73 70 68 65 72 6f 69 64 |6 Earth |spheroid|
|00001b20| 20 64 61 74 61 20 77 68 | 69 63 68 20 69 73 20 6f | data wh|ich is o|
|00001b30| 66 74 65 6e 0a 58 20 2a | 20 75 73 65 64 20 62 79 |ften.X *| used by|
|00001b40| 20 74 68 65 20 55 53 47 | 53 20 2d 20 74 68 65 72 | the USG|S - ther|
|00001b50| 65 20 61 72 65 20 6f 74 | 68 65 72 20 73 70 68 65 |e are ot|her sphe|
|00001b60| 72 6f 69 64 73 20 68 6f | 77 65 76 65 72 20 2d 20 |roids ho|wever - |
|00001b70| 73 65 65 20 74 68 65 0a | 58 20 2a 20 62 6f 6f 6b |see the.|X * book|
|00001b80| 2e 0a 58 20 2a 2f 0a 58 | 20 0a 58 2f 2a 0a 58 20 |..X */.X| .X/*.X |
|00001b90| 2a 20 45 61 72 74 68 20 | 72 61 64 69 69 20 69 6e |* Earth |radii in|
|00001ba0| 20 64 69 66 66 65 72 65 | 6e 74 20 75 6e 69 74 73 | differe|nt units|
|00001bb0| 20 2a 2f 0a 58 23 64 65 | 66 69 6e 65 20 43 4f 4e | */.X#de|fine CON|
|00001bc0| 53 54 5f 45 72 61 64 69 | 75 73 4b 4d 20 28 36 33 |ST_Eradi|usKM (63|
|00001bd0| 37 38 2e 32 30 36 34 29 | 20 09 2f 2a 20 4b 69 6c |78.2064)| ./* Kil|
|00001be0| 6f 6d 65 74 65 72 73 20 | 2a 2f 0a 58 23 64 65 66 |ometers |*/.X#def|
|00001bf0| 69 6e 65 20 43 4f 4e 53 | 54 5f 45 72 61 64 69 75 |ine CONS|T_Eradiu|
|00001c00| 73 4d 49 20 28 43 4f 4e | 53 54 5f 45 72 61 64 69 |sMI (CON|ST_Eradi|
|00001c10| 75 73 4b 4d 2f 31 2e 36 | 30 39 29 09 2f 2a 20 4d |usKM/1.6|09)./* M|
|00001c20| 69 6c 65 73 20 2a 2f 0a | 58 23 64 65 66 69 6e 65 |iles */.|X#define|
|00001c30| 20 43 4f 4e 53 54 5f 45 | 63 09 09 28 30 2e 30 38 | CONST_E|c..(0.08|
|00001c40| 32 32 37 31 38 35 34 29 | 09 2f 2a 20 45 63 63 65 |2271854)|./* Ecce|
|00001c50| 6e 74 72 69 63 69 74 79 | 20 2a 2f 0a 58 23 64 65 |ntricity| */.X#de|
|00001c60| 66 69 6e 65 20 43 4f 4e | 53 54 5f 45 63 73 71 09 |fine CON|ST_Ecsq.|
|00001c70| 28 30 2e 30 30 36 37 36 | 38 36 35 38 29 09 2f 2a |(0.00676|8658)./*|
|00001c80| 20 45 63 63 65 6e 74 72 | 69 63 69 74 79 20 73 71 | Eccentr|icity sq|
|00001c90| 75 61 72 65 64 20 2a 2f | 0a 58 0a 58 0a 58 0a 58 |uared */|.X.X.X.X|
|00001ca0| 2f 2a 0a 58 20 2a 20 54 | 6f 20 6b 65 65 70 20 74 |/*.X * T|o keep t|
|00001cb0| 68 69 6e 67 73 20 73 69 | 6d 70 6c 65 2c 20 61 73 |hings si|mple, as|
|00001cc0| 73 75 6d 65 20 45 61 72 | 74 68 20 72 61 64 69 75 |sume Ear|th radiu|
|00001cd0| 73 20 69 73 20 31 2e 30 | 2e 20 50 72 6f 6a 65 63 |s is 1.0|. Projec|
|00001ce0| 74 65 64 0a 58 20 2a 20 | 63 6f 6f 72 64 69 6e 61 |ted.X * |coordina|
|00001cf0| 74 65 73 20 28 58 20 61 | 6e 64 20 59 20 6f 62 74 |tes (X a|nd Y obt|
|00001d00| 61 69 6e 65 64 20 66 72 | 6f 6d 20 61 6c 62 65 72 |ained fr|om alber|
|00001d10| 73 20 70 72 6f 6a 65 63 | 74 20 28 29 29 20 61 72 |s projec|t ()) ar|
|00001d20| 65 0a 58 20 2a 20 20 64 | 69 6d 65 6e 73 69 6f 6e |e.X * d|imension|
|00001d30| 6c 65 73 73 20 61 6e 64 | 20 6d 61 79 20 62 65 20 |less and| may be |
|00001d40| 6d 75 6c 74 69 70 6c 69 | 65 64 20 62 79 20 61 6e |multipli|ed by an|
|00001d50| 79 20 64 65 73 69 72 65 | 64 20 72 61 64 69 75 73 |y desire|d radius|
|00001d60| 0a 58 20 2a 20 20 74 6f | 20 63 6f 6e 76 65 72 74 |.X * to| convert|
|00001d70| 20 74 6f 20 64 65 73 69 | 72 65 64 20 75 6e 69 74 | to desi|red unit|
|00001d80| 73 20 28 75 73 75 61 6c | 6c 79 20 4b 69 6c 6f 6d |s (usual|ly Kilom|
|00001d90| 65 74 65 72 73 20 6f 72 | 20 4d 69 6c 65 73 29 2e |eters or| Miles).|
|00001da0| 0a 58 2a 2f 0a 58 23 64 | 65 66 69 6e 65 20 43 4f |.X*/.X#d|efine CO|
|00001db0| 4e 53 54 5f 45 72 61 64 | 69 75 73 20 31 2e 30 0a |NST_Erad|ius 1.0.|
|00001dc0| 58 0a 58 2f 2a 20 50 72 | 65 2d 63 6f 6d 70 75 74 |X.X/* Pr|e-comput|
|00001dd0| 65 64 20 76 61 72 69 61 | 62 6c 65 73 20 2a 2f 0a |ed varia|bles */.|
|00001de0| 58 73 74 61 74 69 63 20 | 64 6f 75 62 6c 65 20 6d |Xstatic |double m|
|00001df0| 69 64 64 6c 65 6c 6f 6e | 3b 09 09 2f 2a 20 6c 6f |iddlelon|;../* lo|
|00001e00| 6e 67 69 74 75 64 65 20 | 61 74 20 63 65 6e 74 65 |ngitude |at cente|
|00001e10| 72 20 6f 66 20 6d 61 70 | 20 2a 2f 0a 58 73 74 61 |r of map| */.Xsta|
|00001e20| 74 69 63 20 64 6f 75 62 | 6c 65 20 62 69 67 43 2c |tic doub|le bigC,|
|00001e30| 20 63 6f 6e 65 5f 63 6f | 6e 73 74 2c 20 72 30 3b | cone_co|nst, r0;|
|00001e40| 09 2f 2a 20 53 65 65 20 | 74 68 65 20 72 65 66 65 |./* See |the refe|
|00001e50| 72 65 6e 63 65 20 2a 2f | 0a 58 0a 58 73 74 61 74 |rence */|.X.Xstat|
|00001e60| 69 63 0a 58 63 61 6c 63 | 5f 71 5f 6d 73 71 28 6c |ic.Xcalc|_q_msq(l|
|00001e70| 61 74 2c 20 71 70 2c 20 | 6d 73 71 70 29 0a 58 64 |at, qp, |msqp).Xd|
|00001e80| 6f 75 62 6c 65 20 6c 61 | 74 3b 0a 58 64 6f 75 62 |ouble la|t;.Xdoub|
|00001e90| 6c 65 20 2a 71 70 3b 0a | 58 64 6f 75 62 6c 65 20 |le *qp;.|Xdouble |
|00001ea0| 2a 6d 73 71 70 3b 0a 58 | 2f 2a 0a 58 20 2a 20 47 |*msqp;.X|/*.X * G|
|00001eb0| 69 76 65 6e 20 6c 61 74 | 69 74 75 64 65 2c 20 63 |iven lat|itude, c|
|00001ec0| 61 6c 63 75 6c 61 74 65 | 20 27 71 27 20 5b 65 71 |alculate| 'q' [eq|
|00001ed0| 20 33 2d 31 32 5d 20 0a | 58 20 2a 20 69 66 20 6d | 3-12] .|X * if m|
|00001ee0| 73 71 70 20 69 73 20 21 | 3d 20 4e 55 4c 4c 2c 20 |sqp is !|= NULL, |
|00001ef0| 6d 5e 32 20 20 5b 65 71 | 2e 20 31 32 2d 31 35 5d |m^2 [eq|. 12-15]|
|00001f00| 2e 0a 58 2a 2f 0a 58 7b | 0a 58 09 64 6f 75 62 6c |..X*/.X{|.X.doubl|
|00001f10| 65 20 73 2c 20 63 2c 20 | 65 73 3b 0a 58 0a 58 09 |e s, c, |es;.X.X.|
|00001f20| 73 20 3d 20 73 69 6e 28 | 6c 61 74 29 3b 0a 58 09 |s = sin(|lat);.X.|
|00001f30| 65 73 20 3d 20 73 20 2a | 20 43 4f 4e 53 54 5f 45 |es = s *| CONST_E|
|00001f40| 63 3b 0a 58 0a 58 09 2a | 71 70 20 3d 20 28 31 2e |c;.X.X.*|qp = (1.|
|00001f50| 30 20 2d 20 43 4f 4e 53 | 54 5f 45 63 73 71 29 20 |0 - CONS|T_Ecsq) |
|00001f60| 2a 20 28 28 73 20 2f 20 | 28 31 20 2d 20 65 73 20 |* ((s / |(1 - es |
|00001f70| 2a 20 65 73 29 29 2d 0a | 58 09 09 28 31 20 2f 20 |* es))-.|X..(1 / |
|00001f80| 28 32 20 2a 20 43 4f 4e | 53 54 5f 45 63 29 29 20 |(2 * CON|ST_Ec)) |
|00001f90| 2a 20 6c 6f 67 28 28 31 | 20 2d 20 65 73 29 20 2f |* log((1| - es) /|
|00001fa0| 20 28 31 20 2b 20 65 73 | 29 29 29 3b 0a 58 0a 58 | (1 + es|)));.X.X|
|00001fb0| 09 69 66 20 28 6d 73 71 | 70 20 21 3d 20 4e 55 4c |.if (msq|p != NUL|
|00001fc0| 4c 29 0a 58 09 7b 0a 58 | 09 09 63 20 3d 20 63 6f |L).X.{.X|..c = co|
|00001fd0| 73 28 6c 61 74 29 3b 0a | 58 09 09 2a 6d 73 71 70 |s(lat);.|X..*msqp|
|00001fe0| 20 3d 20 63 20 2a 20 63 | 2f 20 28 31 20 2d 20 65 | = c * c|/ (1 - e|
|00001ff0| 73 20 2a 20 65 73 29 3b | 0a 58 09 7d 0a 58 7d 0a |s * es);|.X.}.X}.|
|00002000| 58 0a 58 0a 58 0a 58 0a | 58 61 6c 62 65 72 73 5f |X.X.X.X.|Xalbers_|
|00002010| 73 65 74 75 70 28 73 6f | 75 74 68 6c 61 74 2c 20 |setup(so|uthlat, |
|00002020| 6e 6f 72 74 68 6c 61 74 | 2c 20 6f 72 69 67 69 6e |northlat|, origin|
|00002030| 6c 6f 6e 2c 20 6f 72 69 | 67 69 6e 6c 61 74 29 0a |lon, ori|ginlat).|
|00002040| 58 64 6f 75 62 6c 65 20 | 73 6f 75 74 68 6c 61 74 |Xdouble |southlat|
|00002050| 2c 20 6e 6f 72 74 68 6c | 61 74 3b 0a 58 64 6f 75 |, northl|at;.Xdou|
|00002060| 62 6c 65 20 6f 72 69 67 | 69 6e 6c 6f 6e 3b 0a 58 |ble orig|inlon;.X|
|00002070| 64 6f 75 62 6c 65 20 6f | 72 69 67 69 6e 6c 61 74 |double o|riginlat|
|00002080| 3b 0a 58 2f 2a 0a 58 20 | 2a 20 50 72 65 2d 63 6f |;.X/*.X |* Pre-co|
|00002090| 6d 70 75 74 65 20 61 20 | 62 75 6e 63 68 20 6f 66 |mpute a |bunch of|
|000020a0| 20 76 61 72 69 61 62 6c | 65 73 20 77 68 69 63 68 | variabl|es which|
|000020b0| 20 61 72 65 20 75 73 65 | 64 20 62 79 20 74 68 65 | are use|d by the|
|000020c0| 20 0a 58 20 2a 20 61 6c | 62 65 72 73 5f 70 72 6f | .X * al|bers_pro|
|000020d0| 6a 65 63 74 28 29 0a 58 | 20 2a 0a 58 20 2a 20 73 |ject().X| *.X * s|
|000020e0| 6f 75 74 68 6c 61 74 20 | 09 53 6f 75 74 68 65 72 |outhlat |.Souther|
|000020f0| 6e 20 6c 61 74 69 74 75 | 64 65 20 66 6f 72 20 41 |n latitu|de for A|
|00002100| 6c 62 65 72 73 20 70 72 | 6f 6a 65 63 74 69 6f 6e |lbers pr|ojection|
|00002110| 0a 58 20 2a 20 6e 6f 72 | 74 68 6c 61 74 09 4e 6f |.X * nor|thlat.No|
|00002120| 72 74 68 65 72 6e 20 6c | 61 74 69 74 75 74 65 20 |rthern l|atitute |
|00002130| 66 6f 72 20 41 6c 62 65 | 72 73 20 70 72 6f 6a 65 |for Albe|rs proje|
|00002140| 63 74 69 6f 6e 0a 58 20 | 2a 20 6f 72 69 67 69 6e |ction.X |* origin|
|00002150| 6c 6f 6e 09 4c 6f 6e 67 | 69 74 75 64 65 20 66 6f |lon.Long|itude fo|
|00002160| 72 20 6f 72 69 67 69 6e | 20 6f 66 20 70 72 6f 6a |r origin| of proj|
|00002170| 65 63 74 65 64 20 6d 61 | 70 0a 58 20 2a 20 6f 72 |ected ma|p.X * or|
|00002180| 69 67 69 6e 6c 61 74 09 | 4c 61 74 69 74 75 64 65 |iginlat.|Latitude|
|00002190| 20 66 6f 72 20 6f 72 69 | 67 69 6e 20 6f 66 20 70 | for ori|gin of p|
|000021a0| 72 6f 6a 65 63 74 65 64 | 20 6d 61 70 20 2d 0a 58 |rojected| map -.X|
|000021b0| 20 2a 09 09 09 09 6f 66 | 74 65 6e 20 28 6e 6f 72 | *....of|ten (nor|
|000021c0| 74 68 6c 61 74 20 2b 20 | 73 6f 75 74 68 6c 61 74 |thlat + |southlat|
|000021d0| 29 20 2f 20 32 0a 58 2a | 2f 0a 58 7b 0a 58 09 64 |) / 2.X*|/.X{.X.d|
|000021e0| 6f 75 62 6c 65 20 71 31 | 2c 20 71 32 2c 20 71 30 |ouble q1|, q2, q0|
|000021f0| 3b 0a 58 09 64 6f 75 62 | 6c 65 20 6d 31 73 71 2c |;.X.doub|le m1sq,|
|00002200| 20 6d 32 73 71 3b 0a 58 | 0a 58 09 6d 69 64 64 6c | m2sq;.X|.X.middl|
|00002210| 65 6c 6f 6e 20 3d 20 6f | 72 69 67 69 6e 6c 6f 6e |elon = o|riginlon|
|00002220| 3b 0a 58 09 0a 58 09 63 | 61 6c 5f 71 5f 6d 73 71 |;.X..X.c|al_q_msq|
|00002230| 28 73 6f 75 74 68 6c 61 | 74 2c 20 26 71 31 2c 20 |(southla|t, &q1, |
|00002240| 26 6d 31 73 71 29 3b 0a | 58 09 63 61 6c 5f 71 5f |&m1sq);.|X.cal_q_|
|00002250| 6d 73 71 28 6e 6f 72 74 | 68 6c 61 74 2c 20 26 71 |msq(nort|hlat, &q|
|00002260| 32 2c 20 26 6d 32 73 71 | 29 3b 0a 58 09 63 61 6c |2, &m2sq|);.X.cal|
|00002270| 5f 71 5f 6d 73 71 28 6f | 72 69 67 69 6e 6c 61 74 |_q_msq(o|riginlat|
|00002280| 2c 20 26 71 30 2c 20 4e | 55 4c 4c 29 3b 0a 58 0a |, &q0, N|ULL);.X.|
|00002290| 58 09 63 6f 6e 65 5f 63 | 6f 6e 73 74 20 3d 20 28 |X.cone_c|onst = (|
|000022a0| 6d 31 73 71 20 2d 20 6d | 32 73 71 29 20 2f 20 28 |m1sq - m|2sq) / (|
|000022b0| 71 32 20 2d 20 71 31 29 | 3b 0a 58 09 62 69 67 43 |q2 - q1)|;.X.bigC|
|000022c0| 20 3d 20 6d 31 73 71 20 | 2b 20 63 6f 6e 65 5f 63 | = m1sq |+ cone_c|
|000022d0| 6f 6e 73 74 20 2a 20 71 | 31 3b 0a 58 09 72 30 20 |onst * q|1;.X.r0 |
|000022e0| 3d 20 43 4f 4e 53 54 5f | 45 72 61 64 69 75 73 20 |= CONST_|Eradius |
|000022f0| 2a 20 73 71 72 74 28 62 | 69 67 43 20 2d 20 63 6f |* sqrt(b|igC - co|
|00002300| 6e 65 5f 63 6f 6e 73 74 | 20 2a 71 30 29 20 2f 20 |ne_const| *q0) / |
|00002310| 63 6f 6e 65 5f 63 6f 6e | 73 74 3b 0a 58 7d 0a 58 |cone_con|st;.X}.X|
|00002320| 0a 58 2f 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X/*****|********|
|00002330| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002340| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002350| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002360| 2a 2a 2f 0a 58 0a 58 61 | 6c 62 65 72 73 5f 70 72 |**/.X.Xa|lbers_pr|
|00002370| 6f 6a 65 63 74 28 6c 6f | 6e 2c 20 6c 61 74 2c 20 |oject(lo|n, lat, |
|00002380| 78 70 2c 20 79 70 29 0a | 58 64 6f 75 62 6c 65 20 |xp, yp).|Xdouble |
|00002390| 6c 6f 6e 2c 20 6c 61 74 | 3b 0a 58 64 6f 75 62 6c |lon, lat|;.Xdoubl|
|000023a0| 65 20 2a 78 70 2c 20 2a | 79 70 3b 0a 58 2f 2a 0a |e *xp, *|yp;.X/*.|
|000023b0| 58 20 2a 20 50 72 6f 6a | 65 63 74 20 6c 6f 6e 2f |X * Proj|ect lon/|
|000023c0| 6c 61 74 20 28 69 6e 20 | 72 61 64 69 61 6e 73 29 |lat (in |radians)|
|000023d0| 20 61 63 63 6f 72 64 69 | 6e 67 20 74 6f 20 61 6c | accordi|ng to al|
|000023e0| 62 65 72 73 5f 73 65 74 | 75 70 20 61 6e 64 20 0a |bers_set|up and .|
|000023f0| 58 20 2a 20 72 65 74 75 | 72 6e 20 74 68 65 20 72 |X * retu|rn the r|
|00002400| 65 73 75 6c 74 73 20 76 | 69 61 20 78 70 2c 20 79 |esults v|ia xp, y|
|00002410| 70 2e 20 55 6e 69 74 73 | 20 6f 66 20 2a 78 70 20 |p. Units| of *xp |
|00002420| 61 6e 64 20 2a 79 70 20 | 61 72 65 20 73 61 6d 65 |and *yp |are same|
|00002430| 0a 58 20 2a 20 61 73 20 | 74 68 65 20 75 6e 69 74 |.X * as |the unit|
|00002440| 73 20 75 73 65 64 20 62 | 79 20 43 4f 4e 53 54 5f |s used b|y CONST_|
|00002450| 45 72 61 64 69 75 73 2e | 0a 58 2a 2f 0a 58 7b 0a |Eradius.|.X*/.X{.|
|00002460| 58 09 64 6f 75 62 6c 65 | 20 71 2c 20 72 2c 20 74 |X.double| q, r, t|
|00002470| 68 65 74 61 3b 0a 58 09 | 63 61 6c 63 5f 71 5f 6d |heta;.X.|calc_q_m|
|00002480| 73 71 28 6c 61 74 2c 20 | 26 71 2c 20 4e 55 4c 4c |sq(lat, |&q, NULL|
|00002490| 29 3b 0a 58 09 74 68 65 | 74 61 20 3d 20 63 6f 6e |);.X.the|ta = con|
|000024a0| 65 5f 63 6f 6e 73 74 20 | 2a 20 28 6c 6f 6e 20 2d |e_const |* (lon -|
|000024b0| 6d 69 64 64 6c 65 6c 6f | 6e 29 3b 0a 58 09 72 20 |middlelo|n);.X.r |
|000024c0| 3d 20 43 4f 4e 53 54 5f | 45 72 61 64 69 75 73 20 |= CONST_|Eradius |
|000024d0| 2a 20 73 71 72 74 28 62 | 69 67 43 20 2d 20 63 6f |* sqrt(b|igC - co|
|000024e0| 6e 65 5f 63 6f 6e 73 74 | 20 2a 20 71 29 20 2f 20 |ne_const| * q) / |
|000024f0| 63 6f 6e 65 5f 63 6f 6e | 73 74 3b 0a 58 09 2a 78 |cone_con|st;.X.*x|
|00002500| 70 20 3d 20 72 20 2a 20 | 73 69 6e 28 74 68 65 74 |p = r * |sin(thet|
|00002510| 61 29 3b 0a 58 09 2a 79 | 70 20 3d 20 72 30 20 2d |a);.X.*y|p = r0 -|
|00002520| 20 72 20 2a 20 63 6f 73 | 28 74 68 65 74 61 29 3b | r * cos|(theta);|
|00002530| 0a 58 7d 0a 58 0a 58 23 | 69 66 64 65 66 20 54 45 |.X}.X.X#|ifdef TE|
|00002540| 53 54 50 52 4f 47 52 41 | 4d 0a 58 0a 58 2f 2a 0a |STPROGRA|M.X.X/*.|
|00002550| 58 20 2a 20 54 65 73 74 | 20 76 61 6c 75 65 20 66 |X * Test| value f|
|00002560| 72 6f 6d 20 74 68 65 20 | 55 53 47 53 20 62 6f 6f |rom the |USGS boo|
|00002570| 6b 2e 20 42 65 63 61 75 | 73 65 20 6f 66 20 72 6f |k. Becau|se of ro|
|00002580| 75 6e 64 6f 66 66 2c 20 | 74 68 65 20 0a 58 20 2a |undoff, |the .X *|
|00002590| 20 61 63 74 75 61 6c 20 | 76 61 6c 75 65 73 20 61 | actual |values a|
|000025a0| 72 65 20 70 72 69 6e 74 | 65 64 20 66 6f 72 20 76 |re print|ed for v|
|000025b0| 69 73 75 61 6c 20 69 6e | 73 70 65 63 74 69 6f 6e |isual in|spection|
|000025c0| 20 72 61 74 68 65 72 20 | 0a 58 20 2a 20 74 68 61 | rather |.X * tha|
|000025d0| 6e 20 67 75 65 73 73 69 | 6e 67 20 77 68 61 74 20 |n guessi|ng what |
|000025e0| 63 6f 6e 73 74 69 74 75 | 74 65 73 20 22 63 6c 6f |constitu|tes "clo|
|000025f0| 73 65 20 65 6e 6f 75 67 | 68 22 2e 0a 58 2a 2f 0a |se enoug|h"..X*/.|
|00002600| 58 2f 2a 20 43 6f 6e 76 | 65 72 74 20 61 20 64 65 |X/* Conv|ert a de|
|00002610| 67 72 65 73 73 2c 20 6d | 69 6e 75 74 65 73 20 70 |gress, m|inutes p|
|00002620| 61 69 72 20 74 6f 20 72 | 61 64 69 61 6e 73 20 2a |air to r|adians *|
|00002630| 2f 0a 58 23 64 65 66 69 | 6e 65 20 44 45 47 5f 4d |/.X#defi|ne DEG_M|
|00002640| 49 4e 32 52 41 44 28 64 | 65 67 72 65 65 73 2c 20 |IN2RAD(d|egrees, |
|00002650| 6d 69 6e 75 74 65 73 29 | 20 5c 0a 58 09 28 28 64 |minutes)| \.X.((d|
|00002660| 6f 75 62 6c 65 29 20 28 | 28 64 65 67 72 65 65 73 |ouble) (|(degrees|
|00002670| 20 2b 20 6d 69 6e 75 74 | 65 73 20 2f 20 36 30 2e | + minut|es / 60.|
|00002680| 30 29 20 2a 20 4d 5f 50 | 49 20 2f 20 31 38 30 2e |0) * M_P|I / 180.|
|00002690| 30 29 29 0a 58 0a 58 23 | 64 65 66 69 6e 65 20 4e |0)).X.X#|define N|
|000026a0| 6c 61 74 20 44 45 47 5f | 4d 49 4e 32 52 41 44 28 |lat DEG_|MIN2RAD(|
|000026b0| 32 39 2c 20 33 30 29 09 | 2f 2a 20 32 39 20 64 65 |29, 30).|/* 29 de|
|000026c0| 67 72 65 65 73 2c 20 33 | 30 27 20 4e 6f 72 74 68 |grees, 3|0' North|
|000026d0| 20 4c 61 74 20 2a 2f 0a | 58 23 64 65 66 69 6e 65 | Lat */.|X#define|
|000026e0| 20 53 6c 61 74 20 44 45 | 47 5f 4d 49 4e 32 52 41 | Slat DE|G_MIN2RA|
|000026f0| 44 28 34 35 2c 20 33 30 | 29 09 0a 58 23 64 65 66 |D(45, 30|)..X#def|
|00002700| 69 6e 65 20 4f 72 69 67 | 69 6e 6c 61 74 20 44 45 |ine Orig|inlat DE|
|00002710| 47 5f 4d 49 4e 32 52 41 | 44 28 32 33 2c 20 30 29 |G_MIN2RA|D(23, 0)|
|00002720| 09 0a 58 23 64 65 66 69 | 6e 65 20 4f 72 69 67 69 |..X#defi|ne Origi|
|00002730| 6e 6c 6f 6e 20 44 45 47 | 5f 4d 49 4e 32 52 41 44 |nlon DEG|_MIN2RAD|
|00002740| 28 2d 39 36 2c 20 30 29 | 20 2f 2a 20 57 65 73 74 |(-96, 0)| /* West|
|00002750| 20 6c 6f 6e 67 69 74 75 | 64 65 20 69 73 20 6e 65 | longitu|de is ne|
|00002760| 67 61 74 69 76 65 20 2a | 2f 0a 58 0a 58 23 64 65 |gative *|/.X.X#de|
|00002770| 66 69 6e 65 20 54 65 73 | 74 6c 61 74 20 44 45 47 |fine Tes|tlat DEG|
|00002780| 5f 4d 49 4e 32 52 41 44 | 28 33 35 2c 20 30 29 0a |_MIN2RAD|(35, 0).|
|00002790| 58 23 64 65 66 69 6e 65 | 20 54 65 73 74 6c 6f 6e |X#define| Testlon|
|000027a0| 20 44 45 47 5f 4d 49 4e | 32 52 41 44 28 2d 37 35 | DEG_MIN|2RAD(-75|
|000027b0| 2c 20 30 29 0a 58 0a 58 | 23 64 65 66 69 6e 65 20 |, 0).X.X|#define |
|000027c0| 54 65 73 74 58 20 31 38 | 38 35 2e 34 37 32 37 0a |TestX 18|85.4727.|
|000027d0| 58 23 64 65 66 69 6e 65 | 20 54 65 73 74 59 20 31 |X#define| TestY 1|
|000027e0| 35 33 35 2e 39 32 35 30 | 0a 58 0a 58 6d 61 69 6e |535.9250|.X.Xmain|
|000027f0| 28 29 0a 58 7b 0a 58 09 | 69 6e 74 20 69 3b 0a 58 |().X{.X.|int i;.X|
|00002800| 09 64 6f 75 62 6c 65 20 | 78 2c 20 79 3b 0a 58 0a |.double |x, y;.X.|
|00002810| 58 2f 2a 20 53 65 74 75 | 70 20 69 73 20 61 6c 73 |X/* Setu|p is als|
|00002820| 6f 20 66 72 6f 6d 20 55 | 53 47 53 20 62 6f 6f 6b |o from U|SGS book|
|00002830| 20 74 65 73 74 20 73 65 | 74 20 2a 2f 0a 58 09 61 | test se|t */.X.a|
|00002840| 6c 62 65 72 73 5f 73 65 | 74 75 70 28 53 6c 61 74 |lbers_se|tup(Slat|
|00002850| 2c 20 4e 6c 61 74 2c 20 | 4f 72 69 67 69 6e 6c 6f |, Nlat, |Originlo|
|00002860| 6e 2c 20 4f 72 69 67 69 | 6e 6c 61 74 29 3b 0a 58 |n, Origi|nlat);.X|
|00002870| 0a 58 09 61 6c 62 65 72 | 73 5f 70 72 6f 6a 65 63 |.X.alber|s_projec|
|00002880| 74 28 54 65 73 74 6c 6f | 6e 2c 20 54 65 73 74 6c |t(Testlo|n, Testl|
|00002890| 61 74 2c 20 26 78 2c 20 | 26 79 29 3b 0a 58 09 70 |at, &x, |&y);.X.p|
|000028a0| 72 69 6e 74 66 28 22 25 | 2e 34 31 66 2c 20 25 2e |rintf("%|.41f, %.|
|000028b0| 34 31 66 20 3d 3f 3d 20 | 25 2e 34 31 66 2c 20 25 |41f =?= |%.41f, %|
|000028c0| 2e 34 31 66 2f 6e 22 2c | 0a 58 09 09 78 20 2a 20 |.41f/n",|.X..x * |
|000028d0| 43 4f 4e 53 54 5f 45 72 | 61 64 69 75 73 4b 4d 2c |CONST_Er|adiusKM,|
|000028e0| 20 79 20 2a 20 43 4f 4e | 53 54 5f 45 72 61 64 69 | y * CON|ST_Eradi|
|000028f0| 75 73 4b 4d 2c 0a 58 09 | 09 54 65 73 74 58 2c 20 |usKM,.X.|.TestX, |
|00002900| 54 65 73 74 59 29 3b 0a | 58 7d 0a 58 23 65 6e 64 |TestY);.|X}.X#end|
|00002910| 69 66 09 09 2f 2a 20 54 | 45 53 54 50 52 4f 47 52 |if../* T|ESTPROGR|
|00002920| 41 4d 20 2a 2f 0a 45 4e | 44 5f 4f 46 5f 46 49 4c |AM */.EN|D_OF_FIL|
|00002930| 45 0a 69 66 20 74 65 73 | 74 20 33 39 39 34 20 2d |E.if tes|t 3994 -|
|00002940| 6e 65 20 60 77 63 20 2d | 63 20 3c 27 41 6c 62 65 |ne `wc -|c <'Albe|
|00002950| 72 73 2e 63 27 60 3b 20 | 74 68 65 6e 0a 20 20 20 |rs.c'`; |then. |
|00002960| 20 65 63 68 6f 20 73 68 | 61 72 3a 20 5c 22 27 41 | echo sh|ar: \"'A|
|00002970| 6c 62 65 72 73 2e 63 27 | 5c 22 20 75 6e 70 61 63 |lbers.c'|\" unpac|
|00002980| 6b 65 64 20 77 69 74 68 | 20 77 72 6f 6e 67 20 73 |ked with| wrong s|
|00002990| 69 7a 65 21 0a 66 69 0a | 23 20 65 6e 64 20 6f 66 |ize!.fi.|# end of|
|000029a0| 20 27 41 6c 62 65 72 73 | 2e 63 27 0a 66 69 0a 69 | 'Albers|.c'.fi.i|
|000029b0| 66 20 74 65 73 74 20 2d | 66 20 27 42 6f 75 6e 64 |f test -|f 'Bound|
|000029c0| 53 70 68 65 72 65 2e 63 | 27 20 2d 61 20 22 24 7b |Sphere.c|' -a "${|
|000029d0| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|000029e0| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en . ec|ho shar:|
|000029f0| 20 57 69 6c 6c 20 6e 6f | 74 20 63 6c 6f 62 62 65 | Will no|t clobbe|
|00002a00| 72 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |r existi|ng file |
|00002a10| 5c 22 27 42 6f 75 6e 64 | 53 70 68 65 72 65 2e 63 |\"'Bound|Sphere.c|
|00002a20| 27 5c 22 0a 65 6c 73 65 | 0a 65 63 68 6f 20 73 68 |'\".else|.echo sh|
|00002a30| 61 72 3a 20 45 78 74 72 | 61 63 74 69 6e 67 20 5c |ar: Extr|acting \|
|00002a40| 22 27 42 6f 75 6e 64 53 | 70 68 65 72 65 2e 63 27 |"'BoundS|phere.c'|
|00002a50| 5c 22 20 5c 28 33 37 36 | 37 20 63 68 61 72 61 63 |\" \(376|7 charac|
|00002a60| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00002a70| 2f 2f 22 20 3e 27 42 6f | 75 6e 64 53 70 68 65 72 |//" >'Bo|undSpher|
|00002a80| 65 2e 63 27 20 3c 3c 27 | 45 4e 44 5f 4f 46 5f 46 |e.c' <<'|END_OF_F|
|00002a90| 49 4c 45 27 0a 58 2f 2a | 0a 58 41 6e 20 45 66 66 |ILE'.X/*|.XAn Eff|
|00002aa0| 69 63 69 65 6e 74 20 42 | 6f 75 6e 64 69 6e 67 20 |icient B|ounding |
|00002ab0| 53 70 68 65 72 65 0a 58 | 62 79 20 4a 61 63 6b 20 |Sphere.X|by Jack |
|00002ac0| 52 69 74 74 65 72 0a 58 | 66 72 6f 6d 20 22 47 72 |Ritter.X|from "Gr|
|00002ad0| 61 70 68 69 63 73 20 47 | 65 6d 73 22 2c 20 41 63 |aphics G|ems", Ac|
|00002ae0| 61 64 65 6d 69 63 20 50 | 72 65 73 73 2c 20 31 39 |ademic P|ress, 19|
|00002af0| 39 30 0a 58 2a 2f 0a 58 | 0a 58 2f 2a 20 52 6f 75 |90.X*/.X|.X/* Rou|
|00002b00| 74 69 6e 65 20 74 6f 20 | 63 61 6c 63 75 6c 61 74 |tine to |calculat|
|00002b10| 65 20 74 69 67 68 74 20 | 62 6f 75 6e 64 69 6e 67 |e tight |bounding|
|00002b20| 20 73 70 68 65 72 65 20 | 6f 76 65 72 20 20 20 20 | sphere |over |
|00002b30| 2a 2f 0a 58 2f 2a 20 61 | 20 73 65 74 20 6f 66 20 |*/.X/* a| set of |
|00002b40| 70 6f 69 6e 74 73 20 69 | 6e 20 33 44 20 2a 2f 0a |points i|n 3D */.|
|00002b50| 58 2f 2a 20 54 68 69 73 | 20 63 6f 6e 74 61 69 6e |X/* This| contain|
|00002b60| 73 20 74 68 65 20 72 6f | 75 74 69 6e 65 20 66 69 |s the ro|utine fi|
|00002b70| 6e 64 5f 62 6f 75 6e 64 | 69 6e 67 5f 73 70 68 65 |nd_bound|ing_sphe|
|00002b80| 72 65 28 29 2c 20 2a 2f | 0a 58 2f 2a 20 74 68 65 |re(), */|.X/* the|
|00002b90| 20 73 74 72 75 63 74 20 | 64 65 66 69 6e 69 74 69 | struct |definiti|
|00002ba0| 6f 6e 2c 20 61 6e 64 20 | 74 68 65 20 67 6c 6f 62 |on, and |the glob|
|00002bb0| 61 6c 73 20 75 73 65 64 | 20 66 6f 72 20 70 61 72 |als used| for par|
|00002bc0| 61 6d 65 74 65 72 73 2e | 20 2a 2f 0a 58 2f 2a 20 |ameters.| */.X/* |
|00002bd0| 54 68 65 20 61 62 73 28 | 29 20 6f 66 20 61 6c 6c |The abs(|) of all|
|00002be0| 20 63 6f 6f 72 64 69 6e | 61 74 65 73 20 6d 75 73 | coordin|ates mus|
|00002bf0| 74 20 62 65 20 3c 20 42 | 49 47 4e 55 4d 42 45 52 |t be < B|IGNUMBER|
|00002c00| 20 2a 2f 0a 58 2f 2a 20 | 43 6f 64 65 20 77 72 69 | */.X/* |Code wri|
|00002c10| 74 74 65 6e 20 62 79 20 | 4a 61 63 6b 20 52 69 74 |tten by |Jack Rit|
|00002c20| 74 65 72 20 61 6e 64 20 | 4c 79 6c 65 20 52 61 69 |ter and |Lyle Rai|
|00002c30| 6e 73 2e 20 2a 2f 0a 58 | 0a 58 23 69 6e 63 6c 75 |ns. */.X|.X#inclu|
|00002c40| 64 65 20 3c 73 74 64 69 | 6f 2e 68 3e 0a 58 23 69 |de <stdi|o.h>.X#i|
|00002c50| 6e 63 6c 75 64 65 20 3c | 6d 61 74 68 2e 68 3e 0a |nclude <|math.h>.|
|00002c60| 58 23 69 6e 63 6c 75 64 | 65 20 22 47 72 61 70 68 |X#includ|e "Graph|
|00002c70| 69 63 73 47 65 6d 73 2e | 68 22 0a 58 0a 58 23 64 |icsGems.|h".X.X#d|
|00002c80| 65 66 69 6e 65 20 42 49 | 47 4e 55 4d 42 45 52 20 |efine BI|GNUMBER |
|00002c90| 31 30 30 30 30 30 30 30 | 30 2e 30 20 20 09 09 2f |10000000|0.0 ../|
|00002ca0| 2a 20 68 75 6e 64 72 65 | 64 20 6d 69 6c 6c 69 6f |* hundre|d millio|
|00002cb0| 6e 20 2a 2f 0a 58 0a 58 | 2f 2a 20 47 4c 4f 42 41 |n */.X.X|/* GLOBA|
|00002cc0| 4c 53 2e 20 54 68 65 73 | 65 20 61 72 65 20 75 73 |LS. Thes|e are us|
|00002cd0| 65 64 20 61 73 20 69 6e | 70 75 74 20 61 6e 64 20 |ed as in|put and |
|00002ce0| 6f 75 74 70 75 74 20 70 | 61 72 61 6d 65 74 65 72 |output p|arameter|
|00002cf0| 73 2e 20 2a 2f 0a 58 20 | 0a 58 73 74 72 75 63 74 |s. */.X |.Xstruct|
|00002d00| 20 50 6f 69 6e 74 33 53 | 74 72 75 63 74 20 63 61 | Point3S|truct ca|
|00002d10| 6c 6c 65 72 5f 70 2c 63 | 65 6e 3b 0a 58 64 6f 75 |ller_p,c|en;.Xdou|
|00002d20| 62 6c 65 20 72 61 64 3b | 0a 58 69 6e 74 20 4e 55 |ble rad;|.Xint NU|
|00002d30| 4d 5f 50 4f 49 4e 54 53 | 3b 0a 58 0a 58 2f 2a 20 |M_POINTS|;.X.X/* |
|00002d40| 43 61 6c 6c 20 77 69 74 | 68 20 6e 6f 20 70 61 72 |Call wit|h no par|
|00002d50| 61 6d 65 74 65 72 73 2e | 20 43 61 6c 6c 65 72 20 |ameters.| Caller |
|00002d60| 6d 75 73 74 20 73 65 74 | 20 4e 55 4d 5f 50 4f 49 |must set| NUM_POI|
|00002d70| 4e 54 53 20 2a 2f 0a 58 | 2f 2a 20 62 65 66 6f 72 |NTS */.X|/* befor|
|00002d80| 65 20 63 61 6c 6c 69 6e | 67 2e 20 2a 2f 0a 58 2f |e callin|g. */.X/|
|00002d90| 2a 20 43 61 6c 6c 65 72 | 20 6d 75 73 74 20 73 75 |* Caller| must su|
|00002da0| 70 70 6c 79 20 74 68 65 | 20 72 6f 75 74 69 6e 65 |pply the| routine|
|00002db0| 20 47 45 54 5f 69 54 48 | 5f 50 4f 49 4e 54 28 69 | GET_iTH|_POINT(i|
|00002dc0| 29 2c 20 77 68 69 63 68 | 20 6c 6f 61 64 73 20 68 |), which| loads h|
|00002dd0| 69 73 20 2a 2f 0a 58 2f | 2a 20 69 74 68 20 70 6f |is */.X/|* ith po|
|00002de0| 69 6e 74 20 69 6e 74 6f | 20 74 68 65 20 67 6c 6f |int into| the glo|
|00002df0| 62 61 6c 20 73 74 72 75 | 63 74 20 63 61 6c 6c 65 |bal stru|ct calle|
|00002e00| 72 5f 70 2e 20 28 30 20 | 3c 3d 20 69 20 3c 20 4e |r_p. (0 |<= i < N|
|00002e10| 55 4d 5f 50 4f 49 4e 54 | 53 29 2e 20 2a 2f 0a 58 |UM_POINT|S). */.X|
|00002e20| 2f 2a 20 54 68 65 20 63 | 61 6c 6c 69 6e 67 20 6f |/* The c|alling o|
|00002e30| 72 64 65 72 20 6f 66 20 | 74 68 65 20 70 6f 69 6e |rder of |the poin|
|00002e40| 74 73 20 69 73 20 69 72 | 72 65 6c 65 76 61 6e 74 |ts is ir|relevant|
|00002e50| 2e 20 2a 2f 0a 58 2f 2a | 20 54 68 65 20 66 69 6e |. */.X/*| The fin|
|00002e60| 61 6c 20 62 6f 75 6e 64 | 69 6e 67 20 73 70 68 65 |al bound|ing sphe|
|00002e70| 72 65 20 77 69 6c 6c 20 | 62 65 20 70 75 74 20 69 |re will |be put i|
|00002e80| 6e 74 6f 20 74 68 65 20 | 67 6c 6f 62 61 6c 73 20 |nto the |globals |
|00002e90| 2a 2f 0a 58 2f 2a 20 63 | 65 6e 20 61 6e 64 20 72 |*/.X/* c|en and r|
|00002ea0| 61 64 2e 20 2a 2f 0a 58 | 0a 58 0a 58 66 69 6e 64 |ad. */.X|.X.Xfind|
|00002eb0| 5f 62 6f 75 6e 64 69 6e | 67 5f 73 70 68 65 72 65 |_boundin|g_sphere|
|00002ec0| 28 29 0a 58 7b 0a 58 72 | 65 67 69 73 74 65 72 20 |().X{.Xr|egister |
|00002ed0| 69 6e 74 20 69 3b 0a 58 | 72 65 67 69 73 74 65 72 |int i;.X|register|
|00002ee0| 20 64 6f 75 62 6c 65 20 | 64 78 2c 64 79 2c 64 7a | double |dx,dy,dz|
|00002ef0| 3b 0a 58 72 65 67 69 73 | 74 65 72 20 64 6f 75 62 |;.Xregis|ter doub|
|00002f00| 6c 65 20 72 61 64 5f 73 | 71 2c 78 73 70 61 6e 2c |le rad_s|q,xspan,|
|00002f10| 79 73 70 61 6e 2c 7a 73 | 70 61 6e 2c 6d 61 78 73 |yspan,zs|pan,maxs|
|00002f20| 70 61 6e 3b 0a 58 64 6f | 75 62 6c 65 20 6f 6c 64 |pan;.Xdo|uble old|
|00002f30| 5f 74 6f 5f 70 2c 6f 6c | 64 5f 74 6f 5f 70 5f 73 |_to_p,ol|d_to_p_s|
|00002f40| 71 2c 6f 6c 64 5f 74 6f | 5f 6e 65 77 3b 0a 58 73 |q,old_to|_new;.Xs|
|00002f50| 74 72 75 63 74 20 50 6f | 69 6e 74 33 53 74 72 75 |truct Po|int3Stru|
|00002f60| 63 74 20 78 6d 69 6e 2c | 78 6d 61 78 2c 79 6d 69 |ct xmin,|xmax,ymi|
|00002f70| 6e 2c 79 6d 61 78 2c 7a | 6d 69 6e 2c 7a 6d 61 78 |n,ymax,z|min,zmax|
|00002f80| 2c 64 69 61 31 2c 64 69 | 61 32 3b 0a 58 0a 58 0a |,dia1,di|a2;.X.X.|
|00002f90| 58 2f 2a 20 46 49 52 53 | 54 20 50 41 53 53 3a 20 |X/* FIRS|T PASS: |
|00002fa0| 66 69 6e 64 20 36 20 6d | 69 6e 69 6d 61 2f 6d 61 |find 6 m|inima/ma|
|00002fb0| 78 69 6d 61 20 70 6f 69 | 6e 74 73 20 2a 2f 0a 58 |xima poi|nts */.X|
|00002fc0| 78 6d 69 6e 2e 78 3d 79 | 6d 69 6e 2e 79 3d 7a 6d |xmin.x=y|min.y=zm|
|00002fd0| 69 6e 2e 7a 3d 20 42 49 | 47 4e 55 4d 42 45 52 3b |in.z= BI|GNUMBER;|
|00002fe0| 20 2f 2a 20 69 6e 69 74 | 69 61 6c 69 7a 65 20 66 | /* init|ialize f|
|00002ff0| 6f 72 20 6d 69 6e 2f 6d | 61 78 20 63 6f 6d 70 61 |or min/m|ax compa|
|00003000| 72 65 20 2a 2f 0a 58 78 | 6d 61 78 2e 78 3d 79 6d |re */.Xx|max.x=ym|
|00003010| 61 78 2e 79 3d 7a 6d 61 | 78 2e 7a 3d 20 2d 42 49 |ax.y=zma|x.z= -BI|
|00003020| 47 4e 55 4d 42 45 52 3b | 0a 58 66 6f 72 20 28 69 |GNUMBER;|.Xfor (i|
|00003030| 3d 30 3b 69 3c 4e 55 4d | 5f 50 4f 49 4e 54 53 3b |=0;i<NUM|_POINTS;|
|00003040| 69 2b 2b 29 0a 58 09 7b | 0a 58 09 47 45 54 5f 69 |i++).X.{|.X.GET_i|
|00003050| 54 48 5f 50 4f 49 4e 54 | 28 69 29 3b 20 2f 2a 20 |TH_POINT|(i); /* |
|00003060| 6c 6f 61 64 20 67 6c 6f | 62 61 6c 20 73 74 72 75 |load glo|bal stru|
|00003070| 63 74 20 63 61 6c 6c 65 | 72 5f 70 20 77 69 74 68 |ct calle|r_p with|
|00003080| 20 2a 2f 0a 58 09 09 09 | 09 09 2f 2a 20 68 69 73 | */.X...|../* his|
|00003090| 20 69 74 68 20 70 6f 69 | 6e 74 2e 20 2a 2f 0a 58 | ith poi|nt. */.X|
|000030a0| 09 69 66 20 28 63 61 6c | 6c 65 72 5f 70 2e 78 3c |.if (cal|ler_p.x<|
|000030b0| 78 6d 69 6e 2e 78 29 0a | 58 09 09 78 6d 69 6e 20 |xmin.x).|X..xmin |
|000030c0| 3d 20 63 61 6c 6c 65 72 | 5f 70 3b 20 2f 2a 20 4e |= caller|_p; /* N|
|000030d0| 65 77 20 78 6d 69 6e 69 | 6d 75 6d 20 70 6f 69 6e |ew xmini|mum poin|
|000030e0| 74 20 2a 2f 0a 58 09 69 | 66 20 28 63 61 6c 6c 65 |t */.X.i|f (calle|
|000030f0| 72 5f 70 2e 78 3e 78 6d | 61 78 2e 78 29 0a 58 09 |r_p.x>xm|ax.x).X.|
|00003100| 09 78 6d 61 78 20 3d 20 | 63 61 6c 6c 65 72 5f 70 |.xmax = |caller_p|
|00003110| 3b 0a 58 09 69 66 20 28 | 63 61 6c 6c 65 72 5f 70 |;.X.if (|caller_p|
|00003120| 2e 79 3c 79 6d 69 6e 2e | 79 29 0a 58 09 09 79 6d |.y<ymin.|y).X..ym|
|00003130| 69 6e 20 3d 20 63 61 6c | 6c 65 72 5f 70 3b 0a 58 |in = cal|ler_p;.X|
|00003140| 09 69 66 20 28 63 61 6c | 6c 65 72 5f 70 2e 79 3e |.if (cal|ler_p.y>|
|00003150| 79 6d 61 78 2e 79 29 0a | 58 09 09 79 6d 61 78 20 |ymax.y).|X..ymax |
|00003160| 3d 20 63 61 6c 6c 65 72 | 5f 70 3b 0a 58 09 69 66 |= caller|_p;.X.if|
|00003170| 20 28 63 61 6c 6c 65 72 | 5f 70 2e 7a 3c 7a 6d 69 | (caller|_p.z<zmi|
|00003180| 6e 2e 7a 29 0a 58 09 09 | 7a 6d 69 6e 20 3d 20 63 |n.z).X..|zmin = c|
|00003190| 61 6c 6c 65 72 5f 70 3b | 0a 58 09 69 66 20 28 63 |aller_p;|.X.if (c|
|000031a0| 61 6c 6c 65 72 5f 70 2e | 7a 3e 7a 6d 61 78 2e 7a |aller_p.|z>zmax.z|
|000031b0| 29 0a 58 09 09 7a 6d 61 | 78 20 3d 20 63 61 6c 6c |).X..zma|x = call|
|000031c0| 65 72 5f 70 3b 0a 58 09 | 7d 0a 58 2f 2a 20 53 65 |er_p;.X.|}.X/* Se|
|000031d0| 74 20 78 73 70 61 6e 20 | 3d 20 64 69 73 74 61 6e |t xspan |= distan|
|000031e0| 63 65 20 62 65 74 77 65 | 65 6e 20 74 68 65 20 32 |ce betwe|en the 2|
|000031f0| 20 70 6f 69 6e 74 73 20 | 78 6d 69 6e 20 26 20 78 | points |xmin & x|
|00003200| 6d 61 78 20 28 73 71 75 | 61 72 65 64 29 20 2a 2f |max (squ|ared) */|
|00003210| 0a 58 64 78 20 3d 20 78 | 6d 61 78 2e 78 20 2d 20 |.Xdx = x|max.x - |
|00003220| 78 6d 69 6e 2e 78 3b 0a | 58 64 79 20 3d 20 78 6d |xmin.x;.|Xdy = xm|
|00003230| 61 78 2e 79 20 2d 20 78 | 6d 69 6e 2e 79 3b 0a 58 |ax.y - x|min.y;.X|
|00003240| 64 7a 20 3d 20 78 6d 61 | 78 2e 7a 20 2d 20 78 6d |dz = xma|x.z - xm|
|00003250| 69 6e 2e 7a 3b 0a 58 78 | 73 70 61 6e 20 3d 20 64 |in.z;.Xx|span = d|
|00003260| 78 2a 64 78 20 2b 20 64 | 79 2a 64 79 20 2b 20 64 |x*dx + d|y*dy + d|
|00003270| 7a 2a 64 7a 3b 0a 58 0a | 58 2f 2a 20 53 61 6d 65 |z*dz;.X.|X/* Same|
|00003280| 20 66 6f 72 20 79 20 26 | 20 7a 20 73 70 61 6e 73 | for y &| z spans|
|00003290| 20 2a 2f 0a 58 64 78 20 | 3d 20 79 6d 61 78 2e 78 | */.Xdx |= ymax.x|
|000032a0| 20 2d 20 79 6d 69 6e 2e | 78 3b 0a 58 64 79 20 3d | - ymin.|x;.Xdy =|
|000032b0| 20 79 6d 61 78 2e 79 20 | 2d 20 79 6d 69 6e 2e 79 | ymax.y |- ymin.y|
|000032c0| 3b 0a 58 64 7a 20 3d 20 | 79 6d 61 78 2e 7a 20 2d |;.Xdz = |ymax.z -|
|000032d0| 20 79 6d 69 6e 2e 7a 3b | 0a 58 79 73 70 61 6e 20 | ymin.z;|.Xyspan |
|000032e0| 3d 20 64 78 2a 64 78 20 | 2b 20 64 79 2a 64 79 20 |= dx*dx |+ dy*dy |
|000032f0| 2b 20 64 7a 2a 64 7a 3b | 0a 58 0a 58 64 78 20 3d |+ dz*dz;|.X.Xdx =|
|00003300| 20 7a 6d 61 78 2e 78 20 | 2d 20 7a 6d 69 6e 2e 78 | zmax.x |- zmin.x|
|00003310| 3b 0a 58 64 79 20 3d 20 | 7a 6d 61 78 2e 79 20 2d |;.Xdy = |zmax.y -|
|00003320| 20 7a 6d 69 6e 2e 79 3b | 0a 58 64 7a 20 3d 20 7a | zmin.y;|.Xdz = z|
|00003330| 6d 61 78 2e 7a 20 2d 20 | 7a 6d 69 6e 2e 7a 3b 0a |max.z - |zmin.z;.|
|00003340| 58 7a 73 70 61 6e 20 3d | 20 64 78 2a 64 78 20 2b |Xzspan =| dx*dx +|
|00003350| 20 64 79 2a 64 79 20 2b | 20 64 7a 2a 64 7a 3b 0a | dy*dy +| dz*dz;.|
|00003360| 58 0a 58 2f 2a 20 53 65 | 74 20 70 6f 69 6e 74 73 |X.X/* Se|t points|
|00003370| 20 64 69 61 31 20 26 20 | 64 69 61 32 20 74 6f 20 | dia1 & |dia2 to |
|00003380| 74 68 65 20 6d 61 78 69 | 6d 61 6c 6c 79 20 73 65 |the maxi|mally se|
|00003390| 70 65 72 61 74 65 64 20 | 70 61 69 72 20 2a 2f 0a |perated |pair */.|
|000033a0| 58 64 69 61 31 20 3d 20 | 78 6d 69 6e 3b 20 64 69 |Xdia1 = |xmin; di|
|000033b0| 61 32 20 3d 20 78 6d 61 | 78 3b 20 2f 2a 20 61 73 |a2 = xma|x; /* as|
|000033c0| 73 75 6d 65 20 78 73 70 | 61 6e 20 62 69 67 67 65 |sume xsp|an bigge|
|000033d0| 73 74 20 2a 2f 0a 58 6d | 61 78 73 70 61 6e 20 3d |st */.Xm|axspan =|
|000033e0| 20 78 73 70 61 6e 3b 0a | 58 69 66 20 28 79 73 70 | xspan;.|Xif (ysp|
|000033f0| 61 6e 3e 6d 61 78 73 70 | 61 6e 29 0a 58 09 7b 0a |an>maxsp|an).X.{.|
|00003400| 58 09 6d 61 78 73 70 61 | 6e 20 3d 20 79 73 70 61 |X.maxspa|n = yspa|
|00003410| 6e 3b 0a 58 09 64 69 61 | 31 20 3d 20 79 6d 69 6e |n;.X.dia|1 = ymin|
|00003420| 3b 20 64 69 61 32 20 3d | 20 79 6d 61 78 3b 0a 58 |; dia2 =| ymax;.X|
|00003430| 09 7d 0a 58 69 66 20 28 | 7a 73 70 61 6e 3e 6d 61 |.}.Xif (|zspan>ma|
|00003440| 78 73 70 61 6e 29 0a 58 | 09 7b 0a 58 09 64 69 61 |xspan).X|.{.X.dia|
|00003450| 31 20 3d 20 7a 6d 69 6e | 3b 20 64 69 61 32 20 3d |1 = zmin|; dia2 =|
|00003460| 20 7a 6d 61 78 3b 0a 58 | 09 7d 0a 58 0a 58 0a 58 | zmax;.X|.}.X.X.X|
|00003470| 2f 2a 20 64 69 61 31 2c | 64 69 61 32 20 69 73 20 |/* dia1,|dia2 is |
|00003480| 61 20 64 69 61 6d 65 74 | 65 72 20 6f 66 20 69 6e |a diamet|er of in|
|00003490| 69 74 69 61 6c 20 73 70 | 68 65 72 65 20 2a 2f 0a |itial sp|here */.|
|000034a0| 58 2f 2a 20 63 61 6c 63 | 20 69 6e 69 74 69 61 6c |X/* calc| initial|
|000034b0| 20 63 65 6e 74 65 72 20 | 2a 2f 0a 58 63 65 6e 2e | center |*/.Xcen.|
|000034c0| 78 20 3d 20 28 64 69 61 | 31 2e 78 2b 64 69 61 32 |x = (dia|1.x+dia2|
|000034d0| 2e 78 29 2f 32 2e 30 3b | 0a 58 63 65 6e 2e 79 20 |.x)/2.0;|.Xcen.y |
|000034e0| 3d 20 28 64 69 61 31 2e | 79 2b 64 69 61 32 2e 79 |= (dia1.|y+dia2.y|
|000034f0| 29 2f 32 2e 30 3b 0a 58 | 63 65 6e 2e 7a 20 3d 20 |)/2.0;.X|cen.z = |
|00003500| 28 64 69 61 31 2e 7a 2b | 64 69 61 32 2e 7a 29 2f |(dia1.z+|dia2.z)/|
|00003510| 32 2e 30 3b 0a 58 2f 2a | 20 63 61 6c 63 75 6c 61 |2.0;.X/*| calcula|
|00003520| 74 65 20 69 6e 69 74 69 | 61 6c 20 72 61 64 69 75 |te initi|al radiu|
|00003530| 73 2a 2a 32 20 61 6e 64 | 20 72 61 64 69 75 73 20 |s**2 and| radius |
|00003540| 2a 2f 0a 58 64 78 20 3d | 20 64 69 61 32 2e 78 2d |*/.Xdx =| dia2.x-|
|00003550| 63 65 6e 2e 78 3b 20 2f | 2a 20 78 20 63 6f 6d 70 |cen.x; /|* x comp|
|00003560| 6f 6e 61 6e 74 20 6f 66 | 20 72 61 64 69 75 73 20 |onant of| radius |
|00003570| 76 65 63 74 6f 72 20 2a | 2f 0a 58 64 79 20 3d 20 |vector *|/.Xdy = |
|00003580| 64 69 61 32 2e 79 2d 63 | 65 6e 2e 79 3b 20 2f 2a |dia2.y-c|en.y; /*|
|00003590| 20 79 20 63 6f 6d 70 6f | 6e 61 6e 74 20 6f 66 20 | y compo|nant of |
|000035a0| 72 61 64 69 75 73 20 76 | 65 63 74 6f 72 20 2a 2f |radius v|ector */|
|000035b0| 0a 58 64 7a 20 3d 20 64 | 69 61 32 2e 7a 2d 63 65 |.Xdz = d|ia2.z-ce|
|000035c0| 6e 2e 7a 3b 20 2f 2a 20 | 7a 20 63 6f 6d 70 6f 6e |n.z; /* |z compon|
|000035d0| 61 6e 74 20 6f 66 20 72 | 61 64 69 75 73 20 76 65 |ant of r|adius ve|
|000035e0| 63 74 6f 72 20 2a 2f 0a | 58 72 61 64 5f 73 71 20 |ctor */.|Xrad_sq |
|000035f0| 3d 20 64 78 2a 64 78 20 | 2b 20 64 79 2a 64 79 20 |= dx*dx |+ dy*dy |
|00003600| 2b 20 64 7a 2a 64 7a 3b | 0a 58 72 61 64 20 3d 20 |+ dz*dz;|.Xrad = |
|00003610| 73 71 72 74 28 72 61 64 | 5f 73 71 29 3b 0a 58 0a |sqrt(rad|_sq);.X.|
|00003620| 58 2f 2a 20 53 45 43 4f | 4e 44 20 50 41 53 53 3a |X/* SECO|ND PASS:|
|00003630| 20 69 6e 63 72 65 6d 65 | 6e 74 20 63 75 72 72 65 | increme|nt curre|
|00003640| 6e 74 20 73 70 68 65 72 | 65 20 2a 2f 0a 58 0a 58 |nt spher|e */.X.X|
|00003650| 66 6f 72 20 28 69 3d 30 | 3b 69 3c 4e 55 4d 5f 50 |for (i=0|;i<NUM_P|
|00003660| 4f 49 4e 54 53 3b 69 2b | 2b 29 0a 58 09 7b 0a 58 |OINTS;i+|+).X.{.X|
|00003670| 09 47 45 54 5f 69 54 48 | 5f 50 4f 49 4e 54 28 69 |.GET_iTH|_POINT(i|
|00003680| 29 3b 20 2f 2a 20 6c 6f | 61 64 20 67 6c 6f 62 61 |); /* lo|ad globa|
|00003690| 6c 20 73 74 72 75 63 74 | 20 63 61 6c 6c 65 72 5f |l struct| caller_|
|000036a0| 70 20 20 2a 2f 0a 58 09 | 09 09 09 09 2f 2a 20 77 |p */.X.|..../* w|
|000036b0| 69 74 68 20 68 69 73 20 | 69 74 68 20 70 6f 69 6e |ith his |ith poin|
|000036c0| 74 2e 20 2a 2f 0a 58 09 | 64 78 20 3d 20 63 61 6c |t. */.X.|dx = cal|
|000036d0| 6c 65 72 5f 70 2e 78 2d | 63 65 6e 2e 78 3b 0a 58 |ler_p.x-|cen.x;.X|
|000036e0| 09 64 79 20 3d 20 63 61 | 6c 6c 65 72 5f 70 2e 79 |.dy = ca|ller_p.y|
|000036f0| 2d 63 65 6e 2e 79 3b 0a | 58 09 64 7a 20 3d 20 63 |-cen.y;.|X.dz = c|
|00003700| 61 6c 6c 65 72 5f 70 2e | 7a 2d 63 65 6e 2e 7a 3b |aller_p.|z-cen.z;|
|00003710| 0a 58 09 6f 6c 64 5f 74 | 6f 5f 70 5f 73 71 20 3d |.X.old_t|o_p_sq =|
|00003720| 20 64 78 2a 64 78 20 2b | 20 64 79 2a 64 79 20 2b | dx*dx +| dy*dy +|
|00003730| 20 64 7a 2a 64 7a 3b 0a | 58 09 69 66 20 28 6f 6c | dz*dz;.|X.if (ol|
|00003740| 64 5f 74 6f 5f 70 5f 73 | 71 20 3e 20 72 61 64 5f |d_to_p_s|q > rad_|
|00003750| 73 71 29 20 09 2f 2a 20 | 64 6f 20 72 2a 2a 32 20 |sq) ./* |do r**2 |
|00003760| 74 65 73 74 20 66 69 72 | 73 74 20 2a 2f 0a 58 09 |test fir|st */.X.|
|00003770| 09 7b 20 09 2f 2a 20 74 | 68 69 73 20 70 6f 69 6e |.{ ./* t|his poin|
|00003780| 74 20 69 73 20 6f 75 74 | 73 69 64 65 20 6f 66 20 |t is out|side of |
|00003790| 63 75 72 72 65 6e 74 20 | 73 70 68 65 72 65 20 2a |current |sphere *|
|000037a0| 2f 0a 58 09 09 6f 6c 64 | 5f 74 6f 5f 70 20 3d 20 |/.X..old|_to_p = |
|000037b0| 73 71 72 74 28 6f 6c 64 | 5f 74 6f 5f 70 5f 73 71 |sqrt(old|_to_p_sq|
|000037c0| 29 3b 0a 58 09 09 2f 2a | 20 63 61 6c 63 20 72 61 |);.X../*| calc ra|
|000037d0| 64 69 75 73 20 6f 66 20 | 6e 65 77 20 73 70 68 65 |dius of |new sphe|
|000037e0| 72 65 20 2a 2f 0a 58 09 | 09 72 61 64 20 3d 20 28 |re */.X.|.rad = (|
|000037f0| 72 61 64 20 2b 20 6f 6c | 64 5f 74 6f 5f 70 29 20 |rad + ol|d_to_p) |
|00003800| 2f 20 32 2e 30 3b 0a 58 | 09 09 72 61 64 5f 73 71 |/ 2.0;.X|..rad_sq|
|00003810| 20 3d 20 72 61 64 2a 72 | 61 64 3b 20 09 2f 2a 20 | = rad*r|ad; ./* |
|00003820| 66 6f 72 20 6e 65 78 74 | 20 72 2a 2a 32 20 63 6f |for next| r**2 co|
|00003830| 6d 70 61 72 65 20 2a 2f | 0a 58 09 09 6f 6c 64 5f |mpare */|.X..old_|
|00003840| 74 6f 5f 6e 65 77 20 3d | 20 6f 6c 64 5f 74 6f 5f |to_new =| old_to_|
|00003850| 70 20 2d 20 72 61 64 3b | 0a 58 09 09 2f 2a 20 63 |p - rad;|.X../* c|
|00003860| 61 6c 63 20 63 65 6e 74 | 65 72 20 6f 66 20 6e 65 |alc cent|er of ne|
|00003870| 77 20 73 70 68 65 72 65 | 20 2a 2f 0a 58 09 09 63 |w sphere| */.X..c|
|00003880| 65 6e 2e 78 20 3d 20 28 | 72 61 64 2a 63 65 6e 2e |en.x = (|rad*cen.|
|00003890| 78 20 2b 20 6f 6c 64 5f | 74 6f 5f 6e 65 77 2a 63 |x + old_|to_new*c|
|000038a0| 61 6c 6c 65 72 5f 70 2e | 78 29 20 2f 20 6f 6c 64 |aller_p.|x) / old|
|000038b0| 5f 74 6f 5f 70 3b 0a 58 | 09 09 63 65 6e 2e 79 20 |_to_p;.X|..cen.y |
|000038c0| 3d 20 28 72 61 64 2a 63 | 65 6e 2e 79 20 2b 20 6f |= (rad*c|en.y + o|
|000038d0| 6c 64 5f 74 6f 5f 6e 65 | 77 2a 63 61 6c 6c 65 72 |ld_to_ne|w*caller|
|000038e0| 5f 70 2e 79 29 20 2f 20 | 6f 6c 64 5f 74 6f 5f 70 |_p.y) / |old_to_p|
|000038f0| 3b 0a 58 09 09 63 65 6e | 2e 7a 20 3d 20 28 72 61 |;.X..cen|.z = (ra|
|00003900| 64 2a 63 65 6e 2e 7a 20 | 2b 20 6f 6c 64 5f 74 6f |d*cen.z |+ old_to|
|00003910| 5f 6e 65 77 2a 63 61 6c | 6c 65 72 5f 70 2e 7a 29 |_new*cal|ler_p.z)|
|00003920| 20 2f 20 6f 6c 64 5f 74 | 6f 5f 70 3b 0a 58 09 09 | / old_t|o_p;.X..|
|00003930| 2f 2a 20 53 75 70 70 72 | 65 73 73 20 69 66 20 64 |/* Suppr|ess if d|
|00003940| 65 73 69 72 65 64 20 2a | 2f 0a 58 09 09 70 72 69 |esired *|/.X..pri|
|00003950| 6e 74 66 28 22 5c 6e 20 | 4e 65 77 20 73 70 68 65 |ntf("\n |New sphe|
|00003960| 72 65 3a 20 63 65 6e 2c | 72 61 64 20 3d 20 25 66 |re: cen,|rad = %f|
|00003970| 20 25 66 20 25 66 20 20 | 20 25 66 22 2c 0a 58 09 | %f %f | %f",.X.|
|00003980| 09 09 63 65 6e 2e 78 2c | 63 65 6e 2e 79 2c 63 65 |..cen.x,|cen.y,ce|
|00003990| 6e 2e 7a 2c 20 72 61 64 | 29 3b 0a 58 09 09 7d 09 |n.z, rad|);.X..}.|
|000039a0| 0a 58 09 7d 0a 58 7d 20 | 09 09 09 20 2f 2a 20 65 |.X.}.X} |... /* e|
|000039b0| 6e 64 20 6f 66 20 66 69 | 6e 64 5f 62 6f 75 6e 64 |nd of fi|nd_bound|
|000039c0| 69 6e 67 5f 73 70 68 65 | 72 65 28 29 20 20 2a 2f |ing_sphe|re() */|
|000039d0| 0a 45 4e 44 5f 4f 46 5f | 46 49 4c 45 0a 69 66 20 |.END_OF_|FILE.if |
|000039e0| 74 65 73 74 20 33 37 36 | 37 20 2d 6e 65 20 60 77 |test 376|7 -ne `w|
|000039f0| 63 20 2d 63 20 3c 27 42 | 6f 75 6e 64 53 70 68 65 |c -c <'B|oundSphe|
|00003a00| 72 65 2e 63 27 60 3b 20 | 74 68 65 6e 0a 20 20 20 |re.c'`; |then. |
|00003a10| 20 65 63 68 6f 20 73 68 | 61 72 3a 20 5c 22 27 42 | echo sh|ar: \"'B|
|00003a20| 6f 75 6e 64 53 70 68 65 | 72 65 2e 63 27 5c 22 20 |oundSphe|re.c'\" |
|00003a30| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|00003a40| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|00003a50| 6e 64 20 6f 66 20 27 42 | 6f 75 6e 64 53 70 68 65 |nd of 'B|oundSphe|
|00003a60| 72 65 2e 63 27 0a 66 69 | 0a 69 66 20 74 65 73 74 |re.c'.fi|.if test|
|00003a70| 20 2d 66 20 27 44 69 73 | 73 6f 6c 76 65 2e 63 27 | -f 'Dis|solve.c'|
|00003a80| 20 2d 61 20 22 24 7b 31 | 7d 22 20 21 3d 20 22 2d | -a "${1|}" != "-|
|00003a90| 63 22 20 3b 20 74 68 65 | 6e 20 0a 20 20 65 63 68 |c" ; the|n . ech|
|00003aa0| 6f 20 73 68 61 72 3a 20 | 57 69 6c 6c 20 6e 6f 74 |o shar: |Will not|
|00003ab0| 20 63 6c 6f 62 62 65 72 | 20 65 78 69 73 74 69 6e | clobber| existin|
|00003ac0| 67 20 66 69 6c 65 20 5c | 22 27 44 69 73 73 6f 6c |g file \|"'Dissol|
|00003ad0| 76 65 2e 63 27 5c 22 0a | 65 6c 73 65 0a 65 63 68 |ve.c'\".|else.ech|
|00003ae0| 6f 20 73 68 61 72 3a 20 | 45 78 74 72 61 63 74 69 |o shar: |Extracti|
|00003af0| 6e 67 20 5c 22 27 44 69 | 73 73 6f 6c 76 65 2e 63 |ng \"'Di|ssolve.c|
|00003b00| 27 5c 22 20 5c 28 34 35 | 39 36 20 63 68 61 72 61 |'\" \(45|96 chara|
|00003b10| 63 74 65 72 73 5c 29 0a | 73 65 64 20 22 73 2f 5e |cters\).|sed "s/^|
|00003b20| 58 2f 2f 22 20 3e 27 44 | 69 73 73 6f 6c 76 65 2e |X//" >'D|issolve.|
|00003b30| 63 27 20 3c 3c 27 45 4e | 44 5f 4f 46 5f 46 49 4c |c' <<'EN|D_OF_FIL|
|00003b40| 45 27 0a 58 2f 2a 20 41 | 20 44 69 67 69 74 61 6c |E'.X/* A| Digital|
|00003b50| 20 44 69 73 73 6f 6c 76 | 65 20 45 66 66 65 63 74 | Dissolv|e Effect|
|00003b60| 0a 58 62 79 20 4d 69 6b | 65 20 4d 6f 72 74 6f 6e |.Xby Mik|e Morton|
|00003b70| 0a 58 66 72 6f 6d 20 22 | 47 72 61 70 68 69 63 73 |.Xfrom "|Graphics|
|00003b80| 20 47 65 6d 73 22 2c 20 | 41 63 61 64 65 6d 69 63 | Gems", |Academic|
|00003b90| 20 50 72 65 73 73 2c 20 | 31 39 39 30 0a 58 2a 2f | Press, |1990.X*/|
|00003ba0| 0a 58 0a 58 2f 2a 0a 58 | 20 2a 20 43 6f 64 65 20 |.X.X/*.X| * Code |
|00003bb0| 66 72 61 67 6d 65 6e 74 | 20 74 6f 20 61 64 76 61 |fragment| to adva|
|00003bc0| 6e 63 65 20 66 72 6f 6d | 20 6f 6e 65 20 65 6c 65 |nce from| one ele|
|00003bd0| 6d 65 6e 74 20 74 6f 20 | 74 68 65 20 6e 65 78 74 |ment to |the next|
|00003be0| 2e 0a 58 20 2a 0a 58 20 | 2a 20 69 6e 74 20 72 65 |..X *.X |* int re|
|00003bf0| 67 3b 09 09 09 09 2f 2a | 20 63 75 72 72 65 6e 74 |g;..../*| current|
|00003c00| 20 73 65 71 75 65 6e 63 | 65 20 65 6c 65 6d 65 6e | sequenc|e elemen|
|00003c10| 74 0a 58 20 2a 20 72 65 | 67 20 3d 20 31 3b 09 09 |t.X * re|g = 1;..|
|00003c20| 09 09 2f 2a 20 73 74 61 | 72 74 20 69 6e 20 61 6e |../* sta|rt in an|
|00003c30| 79 20 6e 6f 6e 2d 7a 65 | 72 6f 20 73 74 61 74 65 |y non-ze|ro state|
|00003c40| 0a 58 20 2a 20 69 66 20 | 28 72 65 67 20 26 20 31 |.X * if |(reg & 1|
|00003c50| 29 09 09 09 09 2f 2a 20 | 69 73 20 74 68 65 20 62 |)..../* |is the b|
|00003c60| 6f 74 74 6f 6d 20 62 69 | 74 20 73 65 74 3f 0a 58 |ottom bi|t set?.X|
|00003c70| 20 2a 20 09 72 65 67 20 | 3d 20 28 72 65 67 20 3e | * .reg |= (reg >|
|00003c80| 3e 31 29 20 5e 20 4d 41 | 53 4b 3b 09 09 2f 2a 20 |>1) ^ MA|SK;../* |
|00003c90| 79 65 73 3a 20 74 6f 73 | 73 20 6f 75 74 20 31 20 |yes: tos|s out 1 |
|00003ca0| 62 69 74 3b 20 58 4f 52 | 20 69 6e 20 6d 61 73 6b |bit; XOR| in mask|
|00003cb0| 0a 58 20 2a 20 65 6c 73 | 65 20 72 65 67 20 3d 20 |.X * els|e reg = |
|00003cc0| 72 65 67 20 3e 3e 31 3b | 09 09 09 2f 2a 20 6e 6f |reg >>1;|.../* no|
|00003cd0| 3a 20 74 6f 73 73 20 6f | 75 74 20 30 20 62 69 74 |: toss o|ut 0 bit|
|00003ce0| 20 0a 58 20 2a 2f 0a 58 | 0a 58 69 6e 74 20 72 61 | .X */.X|.Xint ra|
|00003cf0| 6e 64 6d 61 73 6b 73 5b | 33 32 5d 3b 09 2f 2a 20 |ndmasks[|32];./* |
|00003d00| 47 6f 74 74 61 20 66 69 | 6c 6c 20 74 68 69 73 20 |Gotta fi|ll this |
|00003d10| 69 6e 20 79 6f 75 72 73 | 65 6c 66 2e 20 2a 2f 0a |in yours|elf. */.|
|00003d20| 58 0a 58 64 69 73 73 6f | 6c 76 65 31 20 28 68 65 |X.Xdisso|lve1 (he|
|00003d30| 69 67 68 74 2c 20 77 69 | 64 74 68 29 09 2f 2a 20 |ight, wi|dth)./* |
|00003d40| 66 69 72 73 74 20 76 65 | 72 73 69 6f 6e 20 6f 66 |first ve|rsion of|
|00003d50| 20 74 68 65 20 64 69 73 | 73 6f 6c 76 65 20 09 09 | the dis|solve ..|
|00003d60| 09 09 09 09 09 09 2f 2a | 20 61 6c 67 6f 72 69 74 |....../*| algorit|
|00003d70| 68 6d 20 2a 2f 0a 58 09 | 69 6e 74 20 68 65 69 67 |hm */.X.|int heig|
|00003d80| 68 74 2c 20 77 69 64 74 | 68 3b 09 2f 2a 20 6e 75 |ht, widt|h;./* nu|
|00003d90| 6d 62 65 72 20 6f 66 20 | 72 6f 77 73 2c 20 63 6f |mber of |rows, co|
|00003da0| 6c 75 6d 6e 73 20 2a 2f | 0a 58 7b 0a 58 09 69 6e |lumns */|.X{.X.in|
|00003db0| 74 20 70 69 78 65 6c 73 | 2c 20 6c 61 73 74 6e 75 |t pixels|, lastnu|
|00003dc0| 6d 3b 09 2f 2a 20 6e 75 | 6d 62 65 72 20 6f 66 20 |m;./* nu|mber of |
|00003dd0| 70 69 78 65 6c 73 3b 20 | 2a 2f 0a 58 09 09 09 09 |pixels; |*/.X....|
|00003de0| 2f 2a 20 6c 61 73 74 20 | 70 69 78 65 6c 27 73 20 |/* last |pixel's |
|00003df0| 6e 75 6d 62 65 72 20 2a | 2f 0a 58 09 69 6e 74 20 |number *|/.X.int |
|00003e00| 72 65 67 77 69 64 74 68 | 3b 09 09 2f 2a 20 22 77 |regwidth|;../* "w|
|00003e10| 69 64 74 68 22 20 6f 66 | 20 73 65 71 75 65 6e 63 |idth" of| sequenc|
|00003e20| 65 20 67 65 6e 65 72 61 | 74 6f 72 20 2a 2f 0a 58 |e genera|tor */.X|
|00003e30| 09 72 65 67 69 73 74 65 | 72 20 6c 6f 6e 67 20 6d |.registe|r long m|
|00003e40| 61 73 6b 3b 09 2f 2a 20 | 6d 61 73 6b 20 74 6f 20 |ask;./* |mask to |
|00003e50| 58 4f 52 20 77 69 74 68 | 20 74 6f 2a 2f 0a 58 09 |XOR with| to*/.X.|
|00003e60| 09 09 09 09 2f 2a 20 63 | 72 65 61 74 65 20 73 65 |..../* c|reate se|
|00003e70| 71 75 65 6e 63 65 20 2a | 2f 0a 58 09 72 65 67 69 |quence *|/.X.regi|
|00003e80| 73 74 65 72 20 75 6e 73 | 69 67 6e 65 64 20 6c 6f |ster uns|igned lo|
|00003e90| 6e 67 20 65 6c 65 6d 65 | 6e 74 3b 20 0a 58 09 09 |ng eleme|nt; .X..|
|00003ea0| 09 09 09 2f 2a 20 6f 6e | 65 20 65 6c 65 6d 65 6e |.../* on|e elemen|
|00003eb0| 74 20 6f 66 20 72 61 6e | 64 6f 6d 20 73 65 71 75 |t of ran|dom sequ|
|00003ec0| 65 6e 63 65 20 2a 2f 0a | 58 09 72 65 67 69 73 74 |ence */.|X.regist|
|00003ed0| 65 72 20 69 6e 74 20 72 | 6f 77 2c 20 63 6f 6c 75 |er int r|ow, colu|
|00003ee0| 6d 6e 3b 0a 58 09 09 09 | 09 09 2f 2a 20 72 6f 77 |mn;.X...|../* row|
|00003ef0| 20 61 6e 64 20 63 6f 6c | 75 6d 6e 20 6e 75 6d 62 | and col|umn numb|
|00003f00| 65 72 73 20 66 6f 72 20 | 61 20 70 69 78 65 6c 20 |ers for |a pixel |
|00003f10| 2a 2f 0a 58 0a 58 09 20 | 20 2f 2a 20 46 69 6e 64 |*/.X.X. | /* Find|
|00003f20| 20 73 6d 61 6c 6c 65 73 | 74 20 72 65 67 69 73 74 | smalles|t regist|
|00003f30| 65 72 20 77 68 69 63 68 | 20 70 72 6f 64 75 63 65 |er which| produce|
|00003f40| 73 20 65 6e 6f 75 67 68 | 20 70 69 78 65 6c 20 6e |s enough| pixel n|
|00003f50| 75 6d 62 65 72 73 20 2a | 2f 0a 58 09 20 70 69 78 |umbers *|/.X. pix|
|00003f60| 65 6c 73 20 3d 20 68 65 | 69 67 68 74 20 2a 20 77 |els = he|ight * w|
|00003f70| 69 64 74 68 3b 20 2f 2a | 20 63 6f 6d 70 75 74 65 |idth; /*| compute|
|00003f80| 20 6e 75 6d 62 65 72 20 | 6f 66 20 70 69 78 65 6c | number |of pixel|
|00003f90| 73 20 2a 2f 0a 58 09 09 | 09 09 09 09 09 2f 2a 20 |s */.X..|...../* |
|00003fa0| 74 6f 20 64 69 73 73 6f | 6c 76 65 20 2a 2f 0a 58 |to disso|lve */.X|
|00003fb0| 09 20 6c 61 73 74 6e 75 | 6d 20 3d 20 70 69 78 65 |. lastnu|m = pixe|
|00003fc0| 6c 73 2d 31 3b 09 2f 2a | 20 66 69 6e 64 20 6c 61 |ls-1;./*| find la|
|00003fd0| 73 74 20 65 6c 65 6d 65 | 6e 74 20 28 74 68 65 79 |st eleme|nt (they|
|00003fe0| 20 67 6f 20 30 2e 2e 6c | 61 73 74 6e 75 6d 29 20 | go 0..l|astnum) |
|00003ff0| 2a 2f 0a 58 09 20 72 65 | 67 77 69 64 74 68 20 3d |*/.X. re|gwidth =|
|00004000| 20 62 69 74 77 69 64 74 | 68 20 28 6c 61 73 74 6e | bitwidt|h (lastn|
|00004010| 75 6d 29 3b 20 2f 2a 20 | 68 6f 77 20 77 69 64 65 |um); /* |how wide|
|00004020| 20 6d 75 73 74 20 74 68 | 65 20 2a 2f 0a 58 09 09 | must th|e */.X..|
|00004030| 09 09 09 2f 2a 20 72 65 | 67 69 73 74 65 72 20 62 |.../* re|gister b|
|00004040| 65 3f 20 2a 2f 0a 58 09 | 20 6d 61 73 6b 20 3d 20 |e? */.X.| mask = |
|00004050| 72 61 6e 64 6d 61 73 6b | 73 20 5b 72 65 67 77 69 |randmask|s [regwi|
|00004060| 64 74 68 5d 3b 09 2f 2a | 20 77 68 69 63 68 20 6d |dth];./*| which m|
|00004070| 61 73 6b 20 69 73 20 66 | 6f 72 20 74 68 61 74 20 |ask is f|or that |
|00004080| 77 69 64 74 68 3f 20 2a | 2f 0a 58 0a 58 09 20 2f |width? *|/.X.X. /|
|00004090| 2a 20 4e 6f 77 20 63 79 | 63 6c 65 20 74 68 72 6f |* Now cy|cle thro|
|000040a0| 75 67 68 20 61 6c 6c 20 | 73 65 71 75 65 6e 63 65 |ugh all |sequence|
|000040b0| 20 65 6c 65 6d 65 6e 74 | 73 2e 20 2a 2f 0a 58 0a | element|s. */.X.|
|000040c0| 58 09 20 20 65 6c 65 6d | 65 6e 74 20 3d 20 31 3b |X. elem|ent = 1;|
|000040d0| 09 2f 2a 20 31 73 74 20 | 65 6c 65 6d 65 6e 74 20 |./* 1st |element |
|000040e0| 28 63 6f 75 6c 64 20 62 | 65 20 61 6e 79 20 6e 6f |(could b|e any no|
|000040f0| 6e 7a 65 72 6f 29 20 2a | 2f 0a 58 0a 58 0a 58 09 |nzero) *|/.X.X.X.|
|00004100| 20 20 64 6f 20 7b 0a 58 | 09 20 20 20 20 72 6f 77 | do {.X|. row|
|00004110| 20 3d 20 65 6c 65 6d 65 | 6e 74 20 2f 20 77 69 64 | = eleme|nt / wid|
|00004120| 74 68 3b 09 2f 2a 20 68 | 6f 77 20 6d 61 6e 79 20 |th;./* h|ow many |
|00004130| 72 6f 77 73 20 64 6f 77 | 6e 20 69 73 20 74 68 69 |rows dow|n is thi|
|00004140| 73 20 70 69 78 65 6c 3f | 20 2a 2f 0a 58 09 20 20 |s pixel?| */.X. |
|00004150| 20 20 63 6f 6c 75 6d 6e | 20 3d 20 65 6c 65 6d 65 | column| = eleme|
|00004160| 6e 74 20 25 20 77 69 64 | 74 68 3b 09 2f 2a 20 61 |nt % wid|th;./* a|
|00004170| 6e 64 20 68 6f 77 20 6d | 61 6e 79 20 63 6f 6c 75 |nd how m|any colu|
|00004180| 6d 6e 73 20 61 63 72 6f | 73 73 3f 20 2a 2f 0a 58 |mns acro|ss? */.X|
|00004190| 09 20 20 20 20 69 66 20 | 28 72 6f 77 20 3c 20 68 |. if |(row < h|
|000041a0| 65 69 67 68 74 29 09 2f | 2a 20 69 73 20 74 68 69 |eight)./|* is thi|
|000041b0| 73 20 73 65 71 20 65 6c | 65 6d 65 6e 74 20 69 6e |s seq el|ement in|
|000041c0| 20 74 68 65 20 61 72 72 | 61 79 3f 20 2a 2f 0a 58 | the arr|ay? */.X|
|000041d0| 09 20 20 20 20 20 20 63 | 6f 70 79 20 28 72 6f 77 |. c|opy (row|
|000041e0| 2c 20 63 6f 6c 75 6d 6e | 29 3b 09 2f 2a 20 79 65 |, column|);./* ye|
|000041f0| 73 3a 20 63 6f 70 79 20 | 74 68 65 20 28 72 2c 63 |s: copy |the (r,c|
|00004200| 29 27 74 68 20 70 69 78 | 65 6c 20 2a 2f 0a 58 0a |)'th pix|el */.X.|
|00004210| 58 09 20 20 20 20 2f 2a | 20 43 6f 6d 70 75 74 65 |X. /*| Compute|
|00004220| 20 74 68 65 20 6e 65 78 | 74 20 73 65 71 75 65 6e | the nex|t sequen|
|00004230| 63 65 20 65 6c 65 6d 65 | 6e 74 20 2a 2f 0a 58 09 |ce eleme|nt */.X.|
|00004240| 20 20 20 20 69 66 20 28 | 65 6c 65 6d 65 6e 74 20 | if (|element |
|00004250| 26 20 31 29 09 09 2f 2a | 20 69 73 20 74 68 65 20 |& 1)../*| is the |
|00004260| 6c 6f 77 20 62 69 74 20 | 73 65 74 3f 20 2a 2f 0a |low bit |set? */.|
|00004270| 58 09 20 20 20 20 20 20 | 65 6c 65 6d 65 6e 74 20 |X. |element |
|00004280| 3d 20 28 65 6c 65 6d 65 | 6e 74 20 3e 3e 31 29 5e |= (eleme|nt >>1)^|
|00004290| 6d 61 73 6b 3b 09 2f 2a | 20 79 65 73 3a 20 73 68 |mask;./*| yes: sh|
|000042a0| 69 66 74 20 76 61 6c 75 | 65 2c 20 2a 2f 0a 58 09 |ift valu|e, */.X.|
|000042b0| 09 09 09 09 09 2f 2a 20 | 58 4f 52 20 69 6e 20 6d |...../* |XOR in m|
|000042c0| 61 73 6b 20 2a 2f 0a 58 | 09 20 20 20 20 65 6c 73 |ask */.X|. els|
|000042d0| 65 20 65 6c 65 6d 65 6e | 74 20 3d 20 28 65 6c 65 |e elemen|t = (ele|
|000042e0| 6d 65 6e 74 20 3e 3e 31 | 29 3b 09 2f 2a 20 6e 6f |ment >>1|);./* no|
|000042f0| 3a 20 6a 75 73 74 20 73 | 68 69 66 74 20 74 68 65 |: just s|hift the|
|00004300| 20 76 61 6c 75 65 20 2a | 2f 0a 58 09 20 7d 20 77 | value *|/.X. } w|
|00004310| 68 69 6c 65 20 28 65 6c | 65 6d 65 6e 74 20 21 3d |hile (el|ement !=|
|00004320| 20 31 29 3b 09 09 2f 2a | 20 6c 6f 6f 70 20 75 6e | 1);../*| loop un|
|00004330| 74 69 6c 20 77 65 20 72 | 65 74 75 72 6e 20 20 2a |til we r|eturn *|
|00004340| 2f 0a 58 09 09 09 09 09 | 09 2f 2a 20 74 6f 20 6f |/.X.....|./* to o|
|00004350| 72 69 67 69 6e 61 6c 20 | 65 6c 65 6d 65 6e 74 20 |riginal |element |
|00004360| 2a 2f 0a 58 09 20 63 6f | 70 79 20 28 30 2c 20 30 |*/.X. co|py (0, 0|
|00004370| 29 3b 09 09 2f 2a 20 6b | 6c 75 64 67 65 3a 20 74 |);../* k|ludge: t|
|00004380| 68 65 20 6c 6f 6f 70 20 | 64 6f 65 73 6e 27 74 20 |he loop |doesn't |
|00004390| 70 72 6f 64 75 63 65 20 | 28 30 2c 30 29 20 2a 2f |produce |(0,0) */|
|000043a0| 0a 58 7d 09 09 09 09 09 | 09 2f 2a 20 65 6e 64 20 |.X}.....|./* end |
|000043b0| 6f 66 20 64 69 73 73 6f | 6c 76 65 31 28 29 20 2a |of disso|lve1() *|
|000043c0| 2f 0a 58 0a 58 0a 58 0a | 58 69 6e 74 20 62 69 74 |/.X.X.X.|Xint bit|
|000043d0| 77 69 64 74 68 20 28 4e | 29 09 2f 2a 20 66 69 6e |width (N|)./* fin|
|000043e0| 64 20 22 62 69 74 2d 77 | 69 64 74 68 22 20 6e 65 |d "bit-w|idth" ne|
|000043f0| 65 64 65 64 20 74 6f 20 | 72 65 70 72 65 73 65 6e |eded to |represen|
|00004400| 74 20 4e 20 2a 2f 0a 58 | 09 75 6e 73 69 67 6e 65 |t N */.X|.unsigne|
|00004410| 64 20 69 6e 74 20 4e 3b | 09 2f 2a 20 6e 75 6d 62 |d int N;|./* numb|
|00004420| 65 72 20 74 6f 20 63 6f | 6d 70 75 74 65 20 74 68 |er to co|mpute th|
|00004430| 65 20 77 69 64 74 68 20 | 6f 66 20 2a 2f 0a 58 7b |e width |of */.X{|
|00004440| 0a 58 09 20 69 6e 74 20 | 77 69 64 74 68 20 3d 20 |.X. int |width = |
|00004450| 30 3b 09 2f 2a 20 69 6e | 69 74 69 61 6c 6c 79 2c |0;./* in|itially,|
|00004460| 20 6e 6f 20 62 69 74 73 | 20 6e 65 65 64 65 64 20 | no bits| needed |
|00004470| 74 6f 20 72 65 70 72 65 | 73 65 6e 74 20 4e 20 2a |to repre|sent N *|
|00004480| 2f 0a 58 09 20 77 68 69 | 6c 65 20 28 4e 20 21 3d |/.X. whi|le (N !=|
|00004490| 20 30 29 20 7b 09 2f 2a | 20 6c 6f 6f 70 20 27 74 | 0) {./*| loop 't|
|000044a0| 69 6c 20 4e 20 68 61 73 | 20 62 65 65 6e 20 77 68 |il N has| been wh|
|000044b0| 69 74 74 6c 65 64 20 64 | 6f 77 6e 20 74 6f 20 30 |ittled d|own to 0|
|000044c0| 20 2a 2f 0a 58 09 20 20 | 20 20 4e 20 3e 3e 3d 20 | */.X. | N >>= |
|000044d0| 31 3b 09 09 2f 2a 20 73 | 68 69 66 74 20 4e 20 72 |1;../* s|hift N r|
|000044e0| 69 67 68 74 20 31 20 62 | 69 74 20 28 4e 42 3a 20 |ight 1 b|it (NB: |
|000044f0| 4e 20 69 73 20 75 6e 73 | 69 67 6e 65 64 29 20 2a |N is uns|igned) *|
|00004500| 2f 0a 58 09 20 20 20 20 | 77 69 64 74 68 2b 2b 3b |/.X. |width++;|
|00004510| 09 09 2f 2a 20 61 6e 64 | 20 72 65 6d 65 6d 62 65 |../* and| remembe|
|00004520| 72 20 68 6f 77 20 77 69 | 64 65 20 4e 20 69 73 20 |r how wi|de N is |
|00004530| 2a 2f 0a 58 09 20 20 7d | 09 09 09 2f 2a 20 65 6e |*/.X. }|.../* en|
|00004540| 64 20 6f 66 20 6c 6f 6f | 70 20 73 68 72 69 6e 6b |d of loo|p shrink|
|00004550| 69 6e 67 20 4e 20 64 6f | 77 6e 20 74 6f 20 6e 6f |ing N do|wn to no|
|00004560| 74 68 69 6e 67 20 2a 0a | 58 09 20 20 72 65 74 75 |thing *.|X. retu|
|00004570| 72 6e 20 28 77 69 64 74 | 68 29 3b 09 2f 2a 20 72 |rn (widt|h);./* r|
|00004580| 65 74 75 72 6e 20 62 69 | 74 20 70 6f 73 69 74 69 |eturn bi|t positi|
|00004590| 6f 6e 73 20 63 6f 75 6e | 74 65 64 20 2a 2f 0a 58 |ons coun|ted */.X|
|000045a0| 0a 58 7d 09 09 09 09 09 | 09 2f 2a 20 65 6e 64 20 |.X}.....|./* end |
|000045b0| 6f 66 20 62 69 74 77 69 | 64 74 68 28 29 20 2a 2f |of bitwi|dth() */|
|000045c0| 0a 58 0a 58 0a 58 0a 58 | 64 69 73 73 6f 6c 76 65 |.X.X.X.X|dissolve|
|000045d0| 32 20 28 68 65 69 67 68 | 74 2c 20 77 69 64 74 68 |2 (heigh|t, width|
|000045e0| 29 09 2f 2a 20 66 61 73 | 74 20 76 65 72 73 69 6f |)./* fas|t versio|
|000045f0| 6e 20 6f 66 20 74 68 65 | 20 64 69 73 73 6f 6c 76 |n of the| dissolv|
|00004600| 65 20 61 6c 67 6f 72 69 | 74 68 6d 20 2a 2f 0a 58 |e algori|thm */.X|
|00004610| 09 69 6e 74 20 68 65 69 | 67 68 74 2c 20 77 69 64 |.int hei|ght, wid|
|00004620| 74 68 3b 09 2f 2a 20 6e | 75 6d 62 65 72 20 6f 66 |th;./* n|umber of|
|00004630| 20 72 6f 77 73 2c 20 63 | 6f 6c 75 6d 6e 73 20 2a | rows, c|olumns *|
|00004640| 2f 0a 58 7b 0a 58 09 69 | 6e 74 20 72 77 69 64 74 |/.X{.X.i|nt rwidt|
|00004650| 68 2c 20 63 77 69 64 74 | 68 3b 09 2f 2a 20 62 69 |h, cwidt|h;./* bi|
|00004660| 74 20 77 69 64 74 68 20 | 66 6f 72 20 72 6f 77 73 |t width |for rows|
|00004670| 2c 20 66 6f 72 20 63 6f | 6c 75 6d 6e 73 20 2a 2f |, for co|lumns */|
|00004680| 0a 58 09 69 6e 74 20 72 | 65 67 77 69 64 74 68 3b |.X.int r|egwidth;|
|00004690| 09 09 2f 2a 20 22 77 69 | 64 74 68 22 20 6f 66 20 |../* "wi|dth" of |
|000046a0| 73 65 71 75 65 6e 63 65 | 20 67 65 6e 65 72 61 74 |sequence| generat|
|000046b0| 6f 72 20 2a 2f 0a 58 09 | 72 65 67 69 73 74 65 72 |or */.X.|register|
|000046c0| 20 6c 6f 6e 67 20 6d 61 | 73 6b 3b 09 2f 2a 20 6d | long ma|sk;./* m|
|000046d0| 61 73 6b 20 74 6f 20 58 | 4f 52 20 77 69 74 68 20 |ask to X|OR with |
|000046e0| 74 6f 20 63 72 65 61 74 | 65 20 73 65 71 75 65 6e |to creat|e sequen|
|000046f0| 63 65 20 2a 2f 0a 58 09 | 72 65 67 69 73 74 65 72 |ce */.X.|register|
|00004700| 20 69 6e 74 20 72 6f 77 | 73 68 69 66 74 3b 09 2f | int row|shift;./|
|00004710| 2a 20 73 68 69 66 74 20 | 64 69 73 74 61 6e 63 65 |* shift |distance|
|00004720| 20 74 6f 20 67 65 74 20 | 72 6f 77 20 20 2a 2f 0a | to get |row */.|
|00004730| 58 09 09 09 09 09 09 09 | 2f 2a 20 66 72 6f 6d 20 |X.......|/* from |
|00004740| 65 6c 65 6d 65 6e 74 20 | 2a 2f 0a 58 09 72 65 67 |element |*/.X.reg|
|00004750| 69 73 74 65 72 20 69 6e | 74 20 63 6f 6c 6d 61 73 |ister in|t colmas|
|00004760| 6b 3b 20 2f 2a 20 6d 61 | 73 6b 20 74 6f 20 65 78 |k; /* ma|sk to ex|
|00004770| 74 72 61 63 74 20 63 6f | 6c 75 6d 6e 20 66 72 6f |tract co|lumn fro|
|00004780| 6d 20 65 6c 65 6d 65 6e | 74 20 2a 2f 0a 58 09 72 |m elemen|t */.X.r|
|00004790| 65 67 69 73 74 65 72 20 | 75 6e 73 69 67 6e 65 64 |egister |unsigned|
|000047a0| 20 6c 6f 6e 67 20 65 6c | 65 6d 65 6e 74 3b 20 2f | long el|ement; /|
|000047b0| 2a 20 6f 6e 65 20 65 6c | 65 6d 65 6e 74 20 6f 66 |* one el|ement of|
|000047c0| 20 72 61 6e 64 6f 6d 20 | 2a 2f 20 09 09 09 09 09 | random |*/ .....|
|000047d0| 09 09 09 20 20 20 20 2f | 2a 20 73 65 71 75 65 6e |... /|* sequen|
|000047e0| 63 65 20 2a 2f 0a 58 09 | 72 65 67 69 73 74 65 72 |ce */.X.|register|
|000047f0| 20 69 6e 74 20 72 6f 77 | 2c 20 63 6f 6c 75 6d 6e | int row|, column|
|00004800| 3b 20 20 20 20 2f 2a 20 | 72 6f 77 20 61 6e 64 20 |; /* |row and |
|00004810| 63 6f 6c 75 6d 6e 20 66 | 6f 72 20 6f 6e 65 20 70 |column f|or one p|
|00004820| 69 78 65 6c 20 2a 2f 0a | 58 0a 58 0a 58 09 20 20 |ixel */.|X.X.X. |
|00004830| 2f 2a 20 46 69 6e 64 20 | 74 68 65 20 6d 61 73 6b |/* Find |the mask|
|00004840| 20 74 6f 20 70 72 6f 64 | 75 63 65 20 61 6c 6c 20 | to prod|uce all |
|00004850| 72 6f 77 73 20 61 6e 64 | 20 63 6f 6c 75 6d 6e 73 |rows and| columns|
|00004860| 2e 20 2a 2f 0a 58 0a 58 | 09 72 77 69 64 74 68 20 |. */.X.X|.rwidth |
|00004870| 3d 20 62 69 74 77 69 64 | 74 68 20 28 68 65 69 67 |= bitwid|th (heig|
|00004880| 68 74 29 3b 20 2f 2a 20 | 68 6f 77 20 6d 61 6e 79 |ht); /* |how many|
|00004890| 20 62 69 74 73 20 6e 65 | 65 64 65 64 20 66 6f 72 | bits ne|eded for|
|000048a0| 20 68 65 69 67 68 74 3f | 20 2a 2f 0a 58 09 63 77 | height?| */.X.cw|
|000048b0| 69 64 74 68 20 3d 20 62 | 69 74 77 69 64 74 68 20 |idth = b|itwidth |
|000048c0| 28 77 69 64 74 68 29 3b | 20 20 2f 2a 20 68 6f 77 |(width);| /* how|
|000048d0| 20 6d 61 6e 79 20 62 69 | 74 73 20 6e 65 65 64 65 | many bi|ts neede|
|000048e0| 64 20 66 6f 72 20 77 69 | 64 74 68 3f 20 2a 2f 0a |d for wi|dth? */.|
|000048f0| 58 09 72 65 67 77 69 64 | 74 68 20 3d 20 72 77 69 |X.regwid|th = rwi|
|00004900| 64 74 68 20 2b 20 63 77 | 69 64 74 68 3b 20 2f 2a |dth + cw|idth; /*|
|00004910| 20 68 6f 77 20 77 69 64 | 65 20 6d 75 73 74 20 74 | how wid|e must t|
|00004920| 68 65 20 72 65 67 69 73 | 74 65 72 20 62 65 3f 20 |he regis|ter be? |
|00004930| 2a 2f 0a 58 09 6d 61 73 | 6b 20 3d 20 72 61 6e 64 |*/.X.mas|k = rand|
|00004940| 6d 61 73 6b 73 20 5b 72 | 65 67 77 69 64 74 68 5d |masks [r|egwidth]|
|00004950| 3b 20 2f 2a 20 77 68 69 | 63 68 20 6d 61 73 6b 20 |; /* whi|ch mask |
|00004960| 69 73 20 66 6f 72 20 74 | 68 61 74 20 77 69 64 74 |is for t|hat widt|
|00004970| 68 3f 20 2a 2f 0a 58 0a | 58 20 2f 2a 20 46 69 6e |h? */.X.|X /* Fin|
|00004980| 64 20 76 61 6c 75 65 73 | 20 74 6f 20 65 78 74 72 |d values| to extr|
|00004990| 61 63 74 20 72 6f 77 20 | 61 6e 64 20 63 6f 6c 20 |act row |and col |
|000049a0| 6e 75 6d 62 65 72 73 20 | 66 72 6f 6d 20 65 61 63 |numbers |from eac|
|000049b0| 68 20 65 6c 65 6d 65 6e | 74 2e 20 2a 2f 0a 58 09 |h elemen|t. */.X.|
|000049c0| 72 6f 77 73 68 69 66 74 | 20 3d 20 63 77 69 64 74 |rowshift| = cwidt|
|000049d0| 68 3b 20 2f 2a 20 66 69 | 6e 64 20 64 69 73 74 20 |h; /* fi|nd dist |
|000049e0| 74 6f 20 73 68 69 66 74 | 20 74 6f 20 67 65 74 20 |to shift| to get |
|000049f0| 74 6f 70 20 62 69 74 73 | 20 28 72 6f 77 29 20 2a |top bits| (row) *|
|00004a00| 2f 0a 58 09 63 6f 6c 6d | 61 73 6b 20 3d 20 28 31 |/.X.colm|ask = (1|
|00004a10| 3c 3c 63 77 69 64 74 68 | 29 2d 31 3b 09 2f 2a 20 |<<cwidth|)-1;./* |
|00004a20| 66 69 6e 64 20 6d 61 73 | 6b 20 74 6f 20 65 78 74 |find mas|k to ext|
|00004a30| 72 61 63 74 20 20 2a 2f | 0a 58 09 09 09 09 09 09 |ract */|.X......|
|00004a40| 2f 2a 20 62 6f 74 74 6f | 6d 20 62 69 74 73 20 28 |/* botto|m bits (|
|00004a50| 63 6f 6c 29 20 2a 2f 0a | 58 0a 58 09 20 20 2f 2a |col) */.|X.X. /*|
|00004a60| 20 4e 6f 77 20 63 79 63 | 6c 65 20 74 68 72 6f 75 | Now cyc|le throu|
|00004a70| 67 68 20 61 6c 6c 20 73 | 65 71 75 65 6e 63 65 20 |gh all s|equence |
|00004a80| 65 6c 65 6d 65 6e 74 73 | 2e 20 2a 2f 0a 58 0a 58 |elements|. */.X.X|
|00004a90| 09 65 6c 65 6d 65 6e 74 | 20 3d 20 31 3b 09 2f 2a |.element| = 1;./*|
|00004aa0| 20 31 73 74 20 65 6c 65 | 6d 65 6e 74 20 28 63 6f | 1st ele|ment (co|
|00004ab0| 75 6c 64 20 62 65 20 61 | 6e 79 20 6e 6f 6e 7a 65 |uld be a|ny nonze|
|00004ac0| 72 6f 29 20 2a 2f 0a 58 | 09 64 6f 20 7b 0a 58 09 |ro) */.X|.do {.X.|
|00004ad0| 09 72 6f 77 20 3d 20 65 | 6c 65 6d 65 6e 74 20 3e |.row = e|lement >|
|00004ae0| 3e 20 72 6f 77 73 68 69 | 66 74 3b 20 2f 2a 20 66 |> rowshi|ft; /* f|
|00004af0| 69 6e 64 20 72 6f 77 20 | 6e 75 6d 62 65 72 20 66 |ind row |number f|
|00004b00| 6f 72 20 74 68 69 73 20 | 70 69 78 65 6c 20 2a 2f |or this |pixel */|
|00004b10| 0a 58 09 09 63 6f 6c 75 | 6d 6e 20 3d 20 65 6c 65 |.X..colu|mn = ele|
|00004b20| 6d 65 6e 74 20 26 20 63 | 6f 6c 6d 61 73 6b 3b 20 |ment & c|olmask; |
|00004b30| 2f 2a 20 61 6e 64 20 68 | 6f 77 20 6d 61 6e 79 20 |/* and h|ow many |
|00004b40| 63 6f 6c 75 6d 6e 73 20 | 61 63 72 6f 73 73 3f 20 |columns |across? |
|00004b50| 2a 2f 0a 58 09 09 69 66 | 20 28 28 72 6f 77 20 3c |*/.X..if| ((row <|
|00004b60| 20 68 65 69 67 68 74 29 | 09 2f 2a 20 64 6f 65 73 | height)|./* does|
|00004b70| 20 65 6c 65 6d 65 6e 74 | 20 66 61 6c 6c 20 69 6e | element| fall in|
|00004b80| 20 74 68 65 20 61 72 72 | 61 79 3f 20 2a 2f 0a 58 | the arr|ay? */.X|
|00004b90| 09 09 09 26 26 20 28 63 | 6f 6c 75 6d 6e 20 3c 20 |...&& (c|olumn < |
|00004ba0| 77 69 64 74 68 29 29 20 | 2f 2a 20 2e 2e 2e 6d 75 |width)) |/* ...mu|
|00004bb0| 73 74 20 63 68 65 63 6b | 20 72 6f 77 20 41 4e 44 |st check| row AND|
|00004bc0| 20 63 6f 6c 75 6d 6e 20 | 2a 2f 0a 58 09 09 63 6f | column |*/.X..co|
|00004bd0| 70 79 20 28 72 6f 77 2c | 20 63 6f 6c 75 6d 6e 29 |py (row,| column)|
|00004be0| 3b 20 2f 2a 20 69 6e 20 | 62 6f 75 6e 64 73 3a 20 |; /* in |bounds: |
|00004bf0| 63 6f 70 79 20 74 68 65 | 20 28 72 2c 63 29 27 74 |copy the| (r,c)'t|
|00004c00| 68 20 70 69 78 65 6c 20 | 2a 2f 0a 58 0a 58 09 20 |h pixel |*/.X.X. |
|00004c10| 20 20 20 2f 2a 20 43 6f | 6d 70 75 74 65 20 74 68 | /* Co|mpute th|
|00004c20| 65 20 6e 65 78 74 20 73 | 65 71 75 65 6e 63 65 20 |e next s|equence |
|00004c30| 65 6c 65 6d 65 6e 74 20 | 2a 2f 0a 58 09 09 69 66 |element |*/.X..if|
|00004c40| 20 28 65 6c 65 6d 65 6e | 74 20 26 20 31 29 09 09 | (elemen|t & 1)..|
|00004c50| 2f 2a 20 69 73 20 74 68 | 65 20 6c 6f 77 20 62 69 |/* is th|e low bi|
|00004c60| 74 20 73 65 74 3f 20 2a | 2f 0a 58 09 09 65 6c 65 |t set? *|/.X..ele|
|00004c70| 6d 65 6e 74 20 3d 20 28 | 65 6c 65 6d 65 6e 74 20 |ment = (|element |
|00004c80| 3e 3e 31 29 5e 6d 61 73 | 6b 3b 20 2f 2a 20 79 65 |>>1)^mas|k; /* ye|
|00004c90| 73 3a 20 73 68 69 66 74 | 20 76 61 6c 75 65 2c 20 |s: shift| value, |
|00004ca0| 2f 2a 0a 58 09 09 09 09 | 09 09 2f 2a 20 58 4f 52 |/*.X....|../* XOR|
|00004cb0| 20 69 6e 20 6d 61 73 6b | 20 2a 2f 0a 58 09 09 65 | in mask| */.X..e|
|00004cc0| 6c 73 65 20 65 6c 65 6d | 65 6e 74 20 3d 20 28 65 |lse elem|ent = (e|
|00004cd0| 6c 65 6d 65 6e 74 20 3e | 3e 31 29 3b 20 2f 2a 20 |lement >|>1); /* |
|00004ce0| 6e 6f 3a 20 6a 75 73 74 | 20 73 68 69 66 74 20 74 |no: just| shift t|
|00004cf0| 68 65 20 76 61 6c 75 65 | 20 2a 2f 0a 58 09 7d 20 |he value| */.X.} |
|00004d00| 77 68 69 6c 65 20 28 65 | 6c 65 6d 65 6e 74 20 21 |while (e|lement !|
|00004d10| 3d 20 31 29 3b 20 09 2f | 2a 20 6c 6f 6f 70 20 75 |= 1); ./|* loop u|
|00004d20| 6e 74 69 6c 20 77 65 20 | 72 65 74 75 72 6e 20 74 |ntil we |return t|
|00004d30| 6f 20 2a 2f 0a 58 09 09 | 09 09 09 2f 2a 20 20 6f |o */.X..|.../* o|
|00004d40| 72 69 67 69 6e 61 6c 20 | 65 6c 65 6d 65 6e 74 20 |riginal |element |
|00004d50| 2a 2f 0a 58 0a 58 09 63 | 6f 70 79 20 28 30 2c 20 |*/.X.X.c|opy (0, |
|00004d60| 30 29 3b 09 09 2f 2a 20 | 6b 6c 75 64 67 65 3a 20 |0);../* |kludge: |
|00004d70| 65 6c 65 6d 65 6e 74 20 | 6e 65 76 65 72 20 63 6f |element |never co|
|00004d80| 6d 65 73 20 75 70 20 7a | 65 72 6f 20 2a 2f 0a 58 |mes up z|ero */.X|
|00004d90| 7d 09 09 09 09 09 2f 2a | 20 65 6e 64 20 6f 66 20 |}...../*| end of |
|00004da0| 64 69 73 73 6f 6c 76 65 | 32 28 29 20 2a 2f 0a 45 |dissolve|2() */.E|
|00004db0| 4e 44 5f 4f 46 5f 46 49 | 4c 45 0a 69 66 20 74 65 |ND_OF_FI|LE.if te|
|00004dc0| 73 74 20 34 35 39 36 20 | 2d 6e 65 20 60 77 63 20 |st 4596 |-ne `wc |
|00004dd0| 2d 63 20 3c 27 44 69 73 | 73 6f 6c 76 65 2e 63 27 |-c <'Dis|solve.c'|
|00004de0| 60 3b 20 74 68 65 6e 0a | 20 20 20 20 65 63 68 6f |`; then.| echo|
|00004df0| 20 73 68 61 72 3a 20 5c | 22 27 44 69 73 73 6f 6c | shar: \|"'Dissol|
|00004e00| 76 65 2e 63 27 5c 22 20 | 75 6e 70 61 63 6b 65 64 |ve.c'\" |unpacked|
|00004e10| 20 77 69 74 68 20 77 72 | 6f 6e 67 20 73 69 7a 65 | with wr|ong size|
|00004e20| 21 0a 66 69 0a 23 20 65 | 6e 64 20 6f 66 20 27 44 |!.fi.# e|nd of 'D|
|00004e30| 69 73 73 6f 6c 76 65 2e | 63 27 0a 66 69 0a 69 66 |issolve.|c'.fi.if|
|00004e40| 20 74 65 73 74 20 2d 66 | 20 27 44 6f 75 62 6c 65 | test -f| 'Double|
|00004e50| 4c 69 6e 65 2e 63 27 20 | 2d 61 20 22 24 7b 31 7d |Line.c' |-a "${1}|
|00004e60| 22 20 21 3d 20 22 2d 63 | 22 20 3b 20 74 68 65 6e |" != "-c|" ; then|
|00004e70| 20 0a 20 20 65 63 68 6f | 20 73 68 61 72 3a 20 57 | . echo| shar: W|
|00004e80| 69 6c 6c 20 6e 6f 74 20 | 63 6c 6f 62 62 65 72 20 |ill not |clobber |
|00004e90| 65 78 69 73 74 69 6e 67 | 20 66 69 6c 65 20 5c 22 |existing| file \"|
|00004ea0| 27 44 6f 75 62 6c 65 4c | 69 6e 65 2e 63 27 5c 22 |'DoubleL|ine.c'\"|
|00004eb0| 0a 65 6c 73 65 0a 65 63 | 68 6f 20 73 68 61 72 3a |.else.ec|ho shar:|
|00004ec0| 20 45 78 74 72 61 63 74 | 69 6e 67 20 5c 22 27 44 | Extract|ing \"'D|
|00004ed0| 6f 75 62 6c 65 4c 69 6e | 65 2e 63 27 5c 22 20 5c |oubleLin|e.c'\" \|
|00004ee0| 28 34 36 34 37 20 63 68 | 61 72 61 63 74 65 72 73 |(4647 ch|aracters|
|00004ef0| 5c 29 0a 73 65 64 20 22 | 73 2f 5e 58 2f 2f 22 20 |\).sed "|s/^X//" |
|00004f00| 3e 27 44 6f 75 62 6c 65 | 4c 69 6e 65 2e 63 27 20 |>'Double|Line.c' |
|00004f10| 3c 3c 27 45 4e 44 5f 4f | 46 5f 46 49 4c 45 27 0a |<<'END_O|F_FILE'.|
|00004f20| 58 2f 2a 0a 58 53 79 6d | 6d 65 74 72 69 63 20 44 |X/*.XSym|metric D|
|00004f30| 6f 75 62 6c 65 20 53 74 | 65 70 20 4c 69 6e 65 20 |ouble St|ep Line |
|00004f40| 41 6c 67 6f 72 69 74 68 | 6d 0a 58 62 79 20 42 72 |Algorith|m.Xby Br|
|00004f50| 69 61 6e 20 57 79 76 69 | 6c 6c 0a 58 66 72 6f 6d |ian Wyvi|ll.Xfrom|
|00004f60| 20 22 47 72 61 70 68 69 | 63 73 20 47 65 6d 73 22 | "Graphi|cs Gems"|
|00004f70| 2c 20 41 63 61 64 65 6d | 69 63 20 50 72 65 73 73 |, Academ|ic Press|
|00004f80| 2c 20 31 39 39 30 0a 58 | 2a 2f 0a 58 0a 58 23 64 |, 1990.X|*/.X.X#d|
|00004f90| 65 66 69 6e 65 20 73 77 | 61 70 28 61 2c 62 29 20 |efine sw|ap(a,b) |
|00004fa0| 20 20 20 20 20 20 20 20 | 20 20 7b 61 5e 3d 62 3b | | {a^=b;|
|00004fb0| 20 62 5e 3d 61 3b 20 61 | 5e 3d 62 3b 7d 0a 58 23 | b^=a; a|^=b;}.X#|
|00004fc0| 64 65 66 69 6e 65 20 61 | 62 73 6f 6c 75 74 65 28 |define a|bsolute(|
|00004fd0| 69 2c 6a 2c 6b 29 20 20 | 20 20 20 28 20 28 69 2d |i,j,k) | ( (i-|
|00004fe0| 6a 29 2a 28 6b 20 3d 20 | 28 20 28 69 2d 6a 29 3c |j)*(k = |( (i-j)<|
|00004ff0| 30 20 3f 20 2d 31 20 3a | 20 31 29 29 29 0a 58 69 |0 ? -1 :| 1))).Xi|
|00005000| 6e 74 0a 58 73 79 6d 77 | 75 6c 69 6e 65 28 61 31 |nt.Xsymw|uline(a1|
|00005010| 2c 20 62 31 2c 20 61 32 | 2c 20 62 32 29 20 69 6e |, b1, a2|, b2) in|
|00005020| 74 20 61 31 2c 20 62 31 | 2c 20 61 32 2c 20 62 32 |t a1, b1|, a2, b2|
|00005030| 3b 0a 58 7b 0a 58 09 69 | 6e 74 20 20 20 20 20 20 |;.X{.X.i|nt |
|00005040| 20 20 20 20 20 64 78 2c | 20 64 79 2c 20 69 6e 63 | dx,| dy, inc|
|00005050| 72 31 2c 20 69 6e 63 72 | 32 2c 20 44 2c 20 78 2c |r1, incr|2, D, x,|
|00005060| 20 79 2c 20 78 65 6e 64 | 2c 20 63 2c 20 70 69 78 | y, xend|, c, pix|
|00005070| 65 6c 73 5f 6c 65 66 74 | 3b 0a 58 09 69 6e 74 20 |els_left|;.X.int |
|00005080| 20 20 20 20 20 20 20 20 | 20 20 78 31 2c 20 79 31 | | x1, y1|
|00005090| 3b 0a 58 09 69 6e 74 20 | 20 20 20 20 20 20 20 20 |;.X.int | |
|000050a0| 20 20 73 69 67 6e 5f 78 | 2c 20 73 69 67 6e 5f 79 | sign_x|, sign_y|
|000050b0| 2c 20 73 74 65 70 2c 20 | 72 65 76 65 72 73 65 2c |, step, |reverse,|
|000050c0| 20 69 3b 0a 58 0a 58 09 | 64 78 20 3d 20 61 62 73 | i;.X.X.|dx = abs|
|000050d0| 6f 6c 75 74 65 28 61 32 | 2c 20 61 31 2c 20 73 69 |olute(a2|, a1, si|
|000050e0| 67 6e 5f 78 29 3b 0a 58 | 09 64 79 20 3d 20 61 62 |gn_x);.X|.dy = ab|
|000050f0| 73 6f 6c 75 74 65 28 62 | 32 2c 20 62 31 2c 20 73 |solute(b|2, b1, s|
|00005100| 69 67 6e 5f 79 29 3b 0a | 58 09 2f 2a 20 64 65 63 |ign_y);.|X./* dec|
|00005110| 69 64 65 20 69 6e 63 72 | 65 6d 65 6e 74 20 73 69 |ide incr|ement si|
|00005120| 67 6e 20 62 79 20 74 68 | 65 20 73 6c 6f 70 65 20 |gn by th|e slope |
|00005130| 73 69 67 6e 20 2a 2f 0a | 58 09 69 66 20 28 73 69 |sign */.|X.if (si|
|00005140| 67 6e 5f 78 20 3d 3d 20 | 73 69 67 6e 5f 79 29 0a |gn_x == |sign_y).|
|00005150| 58 09 09 73 74 65 70 20 | 3d 20 31 3b 0a 58 09 65 |X..step |= 1;.X.e|
|00005160| 6c 73 65 0a 58 09 09 73 | 74 65 70 20 3d 20 2d 31 |lse.X..s|tep = -1|
|00005170| 3b 0a 58 0a 58 09 69 66 | 20 28 64 79 20 3e 20 64 |;.X.X.if| (dy > d|
|00005180| 78 29 20 7b 09 09 2f 2a | 20 63 68 6f 6f 73 65 73 |x) {../*| chooses|
|00005190| 20 61 78 69 73 20 6f 66 | 20 67 72 65 61 74 65 73 | axis of| greates|
|000051a0| 74 20 6d 6f 76 65 6d 65 | 6e 74 20 28 6d 61 6b 65 |t moveme|nt (make|
|000051b0| 0a 58 09 09 09 09 20 09 | 09 20 2a 20 64 78 29 20 |.X.... .|. * dx) |
|000051c0| 2a 2f 0a 58 09 09 73 77 | 61 70 28 61 31 2c 20 62 |*/.X..sw|ap(a1, b|
|000051d0| 31 29 3b 0a 58 09 09 73 | 77 61 70 28 61 32 2c 20 |1);.X..s|wap(a2, |
|000051e0| 62 32 29 3b 0a 58 09 09 | 73 77 61 70 28 64 78 2c |b2);.X..|swap(dx,|
|000051f0| 20 64 79 29 3b 0a 58 09 | 09 72 65 76 65 72 73 65 | dy);.X.|.reverse|
|00005200| 20 3d 20 31 3b 0a 58 09 | 7d 20 65 6c 73 65 0a 58 | = 1;.X.|} else.X|
|00005210| 09 09 72 65 76 65 72 73 | 65 20 3d 20 30 3b 0a 58 |..revers|e = 0;.X|
|00005220| 09 2f 2a 20 6e 6f 74 65 | 20 65 72 72 6f 72 20 63 |./* note| error c|
|00005230| 68 65 63 6b 20 66 6f 72 | 20 64 78 3d 3d 30 20 73 |heck for| dx==0 s|
|00005240| 68 6f 75 6c 64 20 62 65 | 20 69 6e 63 6c 75 64 65 |hould be| include|
|00005250| 64 20 68 65 72 65 20 2a | 2f 0a 58 09 69 66 20 28 |d here *|/.X.if (|
|00005260| 61 31 20 3e 20 61 32 29 | 20 7b 09 09 2f 2a 20 73 |a1 > a2)| {../* s|
|00005270| 74 61 72 74 20 66 72 6f | 6d 20 74 68 65 20 73 6d |tart fro|m the sm|
|00005280| 61 6c 6c 65 72 20 63 6f | 6f 72 64 69 6e 61 74 65 |aller co|ordinate|
|00005290| 20 2a 2f 0a 58 09 09 78 | 20 3d 20 61 32 3b 0a 58 | */.X..x| = a2;.X|
|000052a0| 09 09 79 20 3d 20 62 32 | 3b 0a 58 09 09 78 31 20 |..y = b2|;.X..x1 |
|000052b0| 3d 20 61 31 3b 0a 58 09 | 09 79 31 20 3d 20 62 31 |= a1;.X.|.y1 = b1|
|000052c0| 3b 0a 58 09 7d 20 65 6c | 73 65 20 7b 0a 58 09 09 |;.X.} el|se {.X..|
|000052d0| 78 20 3d 20 61 31 3b 0a | 58 09 09 79 20 3d 20 62 |x = a1;.|X..y = b|
|000052e0| 31 3b 0a 58 09 09 78 31 | 20 3d 20 61 32 3b 0a 58 |1;.X..x1| = a2;.X|
|000052f0| 09 09 79 31 20 3d 20 62 | 32 3b 0a 58 09 7d 0a 58 |..y1 = b|2;.X.}.X|
|00005300| 0a 58 0a 58 09 2f 2a 20 | 4e 6f 74 65 20 64 78 3d |.X.X./* |Note dx=|
|00005310| 6e 20 69 6d 70 6c 69 65 | 73 20 30 20 2d 20 6e 20 |n implie|s 0 - n |
|00005320| 6f 72 20 28 64 78 2b 31 | 29 20 70 69 78 65 6c 73 |or (dx+1|) pixels|
|00005330| 20 74 6f 20 62 65 20 73 | 65 74 20 2a 2f 0a 58 09 | to be s|et */.X.|
|00005340| 2f 2a 20 47 6f 20 72 6f | 75 6e 64 20 6c 6f 6f 70 |/* Go ro|und loop|
|00005350| 20 64 78 2f 34 20 74 69 | 6d 65 73 20 74 68 65 6e | dx/4 ti|mes then|
|00005360| 20 70 6c 6f 74 20 6c 61 | 73 74 20 30 2c 31 2c 32 | plot la|st 0,1,2|
|00005370| 20 6f 72 20 33 20 70 69 | 78 65 6c 73 20 2a 2f 0a | or 3 pi|xels */.|
|00005380| 58 09 2f 2a 20 49 6e 20 | 66 61 63 74 20 28 64 78 |X./* In |fact (dx|
|00005390| 2d 31 29 2f 34 20 61 73 | 20 32 20 70 69 78 65 6c |-1)/4 as| 2 pixel|
|000053a0| 73 20 61 72 65 20 61 6c | 72 65 61 64 79 20 70 6c |s are al|ready pl|
|000053b0| 6f 74 74 74 65 64 20 2a | 2f 0a 58 09 78 65 6e 64 |ottted *|/.X.xend|
|000053c0| 20 3d 20 28 64 78 20 2d | 20 31 29 20 2f 20 34 3b | = (dx -| 1) / 4;|
|000053d0| 0a 58 09 70 69 78 65 6c | 73 5f 6c 65 66 74 20 3d |.X.pixel|s_left =|
|000053e0| 20 28 64 78 20 2d 20 31 | 29 20 25 20 34 3b 09 2f | (dx - 1|) % 4;./|
|000053f0| 2a 20 6e 75 6d 62 65 72 | 20 6f 66 20 70 69 78 65 |* number| of pixe|
|00005400| 6c 73 20 6c 65 66 74 20 | 6f 76 65 72 20 61 74 20 |ls left |over at |
|00005410| 74 68 65 0a 58 09 09 09 | 09 09 20 09 09 09 20 2a |the.X...|.. ... *|
|00005420| 20 65 6e 64 20 2a 2f 0a | 58 09 70 6c 6f 74 28 78 | end */.|X.plot(x|
|00005430| 2c 20 79 2c 20 72 65 76 | 65 72 73 65 29 3b 0a 58 |, y, rev|erse);.X|
|00005440| 09 70 6c 6f 74 28 78 31 | 2c 20 79 31 2c 20 72 65 |.plot(x1|, y1, re|
|00005450| 76 65 72 73 65 29 3b 09 | 2f 2a 20 70 6c 6f 74 20 |verse);.|/* plot |
|00005460| 66 69 72 73 74 20 74 77 | 6f 20 70 6f 69 6e 74 73 |first tw|o points|
|00005470| 20 2a 2f 0a 58 09 69 6e | 63 72 32 20 3d 20 34 20 | */.X.in|cr2 = 4 |
|00005480| 2a 20 64 79 20 2d 20 32 | 20 2a 20 64 78 3b 0a 58 |* dy - 2| * dx;.X|
|00005490| 09 69 66 20 28 69 6e 63 | 72 32 20 3c 20 30 29 20 |.if (inc|r2 < 0) |
|000054a0| 7b 09 2f 2a 20 73 6c 6f | 70 65 20 6c 65 73 73 20 |{./* slo|pe less |
|000054b0| 74 68 61 6e 20 31 2f 32 | 20 2a 2f 0a 58 09 09 63 |than 1/2| */.X..c|
|000054c0| 20 3d 20 32 20 2a 20 64 | 79 3b 0a 58 09 09 69 6e | = 2 * d|y;.X..in|
|000054d0| 63 72 31 20 3d 20 32 20 | 2a 20 63 3b 0a 58 09 09 |cr1 = 2 |* c;.X..|
|000054e0| 44 20 3d 20 69 6e 63 72 | 31 20 2d 20 64 78 3b 0a |D = incr|1 - dx;.|
|000054f0| 58 0a 58 09 09 66 6f 72 | 20 28 69 20 3d 20 30 3b |X.X..for| (i = 0;|
|00005500| 20 69 20 3c 20 78 65 6e | 64 3b 20 69 2b 2b 29 20 | i < xen|d; i++) |
|00005510| 7b 09 2f 2a 20 70 6c 6f | 74 74 69 6e 67 20 6c 6f |{./* plo|tting lo|
|00005520| 6f 70 20 2a 2f 0a 58 09 | 09 09 2b 2b 78 3b 0a 58 |op */.X.|..++x;.X|
|00005530| 09 09 09 2d 2d 78 31 3b | 0a 58 09 09 09 69 66 20 |...--x1;|.X...if |
|00005540| 28 44 20 3c 20 30 29 20 | 7b 0a 58 20 20 20 20 20 |(D < 0) |{.X |
|00005550| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 09 09 09 | | ...|
|00005560| 2f 2a 20 70 61 74 74 65 | 72 6e 20 31 20 66 6f 72 |/* patte|rn 1 for|
|00005570| 77 61 72 64 73 20 2a 2f | 0a 58 09 09 09 09 70 6c |wards */|.X....pl|
|00005580| 6f 74 28 78 2c 20 79 2c | 20 72 65 76 65 72 73 65 |ot(x, y,| reverse|
|00005590| 29 3b 0a 58 09 09 09 09 | 70 6c 6f 74 28 2b 2b 78 |);.X....|plot(++x|
|000055a0| 2c 20 79 2c 20 72 65 76 | 65 72 73 65 29 3b 0a 58 |, y, rev|erse);.X|
|000055b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000055c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000055d0| 2f 2a 20 70 61 74 74 65 | 72 6e 20 31 20 62 61 63 |/* patte|rn 1 bac|
|000055e0| 6b 77 61 72 64 73 20 2a | 2f 0a 58 09 09 09 09 70 |kwards *|/.X....p|
|000055f0| 6c 6f 74 28 78 31 2c 20 | 79 31 2c 20 72 65 76 65 |lot(x1, |y1, reve|
|00005600| 72 73 65 29 3b 0a 58 09 | 09 09 09 70 6c 6f 74 28 |rse);.X.|...plot(|
|00005610| 2d 2d 78 31 2c 20 79 31 | 2c 20 72 65 76 65 72 73 |--x1, y1|, revers|
|00005620| 65 29 3b 0a 58 09 09 09 | 09 44 20 2b 3d 20 69 6e |e);.X...|.D += in|
|00005630| 63 72 31 3b 0a 58 09 09 | 09 7d 20 65 6c 73 65 20 |cr1;.X..|.} else |
|00005640| 7b 0a 58 09 09 09 09 69 | 66 20 28 44 20 3c 20 63 |{.X....i|f (D < c|
|00005650| 29 20 7b 0a 58 09 09 09 | 09 09 2f 2a 20 70 61 74 |) {.X...|../* pat|
|00005660| 74 65 72 6e 20 32 20 66 | 6f 72 77 61 72 64 73 20 |tern 2 f|orwards |
|00005670| 2a 2f 0a 58 09 09 09 09 | 09 70 6c 6f 74 28 78 2c |*/.X....|.plot(x,|
|00005680| 20 79 2c 20 72 65 76 65 | 72 73 65 29 3b 0a 58 09 | y, reve|rse);.X.|
|00005690| 09 09 09 09 70 6c 6f 74 | 28 2b 2b 78 2c 20 79 20 |....plot|(++x, y |
|000056a0| 2b 3d 20 73 74 65 70 2c | 20 72 65 76 65 72 73 65 |+= step,| reverse|
|000056b0| 29 3b 0a 58 09 09 09 09 | 09 2f 2a 20 70 61 74 74 |);.X....|./* patt|
|000056c0| 65 72 6e 20 32 20 62 61 | 63 6b 77 61 72 64 73 20 |ern 2 ba|ckwards |
|000056d0| 2a 2f 0a 58 09 09 09 09 | 09 70 6c 6f 74 28 78 31 |*/.X....|.plot(x1|
|000056e0| 2c 20 79 31 2c 20 72 65 | 76 65 72 73 65 29 3b 0a |, y1, re|verse);.|
|000056f0| 58 09 09 09 09 09 70 6c | 6f 74 28 2d 2d 78 31 2c |X.....pl|ot(--x1,|
|00005700| 20 79 31 20 2d 3d 20 73 | 74 65 70 2c 20 72 65 76 | y1 -= s|tep, rev|
|00005710| 65 72 73 65 29 3b 09 0a | 58 09 09 09 09 7d 20 65 |erse);..|X....} e|
|00005720| 6c 73 65 20 7b 0a 58 09 | 09 09 09 20 20 20 20 20 |lse {.X.|... |
|00005730| 20 20 20 2f 2a 20 70 61 | 74 74 65 72 6e 20 33 20 | /* pa|ttern 3 |
|00005740| 66 6f 72 77 61 72 64 73 | 20 2a 2f 0a 58 09 09 09 |forwards| */.X...|
|00005750| 09 09 70 6c 6f 74 28 78 | 2c 20 79 20 2b 3d 20 73 |..plot(x|, y += s|
|00005760| 74 65 70 2c 20 72 65 76 | 65 72 73 65 29 3b 0a 58 |tep, rev|erse);.X|
|00005770| 09 09 09 09 09 70 6c 6f | 74 28 2b 2b 78 2c 20 79 |.....plo|t(++x, y|
|00005780| 2c 20 72 65 76 65 72 73 | 65 29 3b 0a 58 09 09 09 |, revers|e);.X...|
|00005790| 09 09 2f 2a 20 70 61 74 | 74 65 72 6e 20 33 20 62 |../* pat|tern 3 b|
|000057a0| 61 63 6b 77 61 72 64 73 | 20 2a 2f 0a 58 09 09 09 |ackwards| */.X...|
|000057b0| 09 09 70 6c 6f 74 28 78 | 31 2c 20 79 31 20 2d 3d |..plot(x|1, y1 -=|
|000057c0| 20 73 74 65 70 2c 20 72 | 65 76 65 72 73 65 29 3b | step, r|everse);|
|000057d0| 0a 58 09 09 09 09 09 70 | 6c 6f 74 28 2d 2d 78 31 |.X.....p|lot(--x1|
|000057e0| 2c 20 79 31 2c 20 72 65 | 76 65 72 73 65 29 3b 0a |, y1, re|verse);.|
|000057f0| 58 09 09 09 09 7d 0a 58 | 09 09 09 09 44 20 2b 20 |X....}.X|....D + |
|00005800| 3d 20 69 6e 63 72 32 3b | 0a 58 09 09 09 7d 0a 58 |= incr2;|.X...}.X|
|00005810| 09 09 7d 09 09 2f 2a 20 | 65 6e 64 20 66 6f 72 20 |..}../* |end for |
|00005820| 2a 2f 0a 58 0a 58 09 09 | 2f 2a 20 70 6c 6f 74 20 |*/.X.X..|/* plot |
|00005830| 6c 61 73 74 20 70 61 74 | 74 65 72 6e 20 2a 2f 0a |last pat|tern */.|
|00005840| 58 09 09 69 66 20 28 70 | 69 78 65 6c 73 5f 6c 65 |X..if (p|ixels_le|
|00005850| 66 74 29 20 7b 0a 58 09 | 09 09 69 66 20 28 44 20 |ft) {.X.|..if (D |
|00005860| 3c 20 30 29 20 7b 0a 58 | 09 09 09 09 70 6c 6f 74 |< 0) {.X|....plot|
|00005870| 28 2b 2b 78 2c 20 79 2c | 20 72 65 76 65 72 73 65 |(++x, y,| reverse|
|00005880| 29 3b 09 2f 2a 20 70 61 | 74 74 65 72 6e 20 31 20 |);./* pa|ttern 1 |
|00005890| 2a 2f 0a 58 09 09 09 09 | 69 66 20 28 70 69 78 65 |*/.X....|if (pixe|
|000058a0| 6c 73 5f 6c 65 66 74 20 | 3e 20 31 29 0a 58 09 09 |ls_left |> 1).X..|
|000058b0| 09 09 09 70 6c 6f 74 28 | 2b 2b 78 2c 20 79 2c 20 |...plot(|++x, y, |
|000058c0| 72 65 76 65 72 73 65 29 | 3b 0a 58 09 09 09 09 69 |reverse)|;.X....i|
|000058d0| 66 20 28 70 69 78 65 6c | 73 5f 6c 65 66 74 20 3e |f (pixel|s_left >|
|000058e0| 20 32 29 0a 58 09 09 09 | 09 09 70 6c 6f 74 28 2d | 2).X...|..plot(-|
|000058f0| 2d 78 31 2c 20 79 31 2c | 20 72 65 76 65 72 73 65 |-x1, y1,| reverse|
|00005900| 29 3b 0a 58 09 09 09 7d | 20 65 6c 73 65 20 7b 0a |);.X...}| else {.|
|00005910| 58 09 09 09 09 69 66 20 | 28 44 20 3c 20 63 29 20 |X....if |(D < c) |
|00005920| 7b 0a 58 09 09 09 09 09 | 70 6c 6f 74 28 2b 2b 78 |{.X.....|plot(++x|
|00005930| 2c 20 79 2c 20 72 65 76 | 65 72 73 65 29 3b 09 2f |, y, rev|erse);./|
|00005940| 2a 20 70 61 74 74 65 72 | 6e 20 32 20 20 2a 2f 0a |* patter|n 2 */.|
|00005950| 58 09 09 09 09 09 69 66 | 20 28 70 69 78 65 6c 73 |X.....if| (pixels|
|00005960| 5f 6c 65 66 74 20 3e 20 | 31 29 0a 58 09 09 09 09 |_left > |1).X....|
|00005970| 09 09 70 6c 6f 74 28 2b | 2b 78 2c 20 79 20 2b 3d |..plot(+|+x, y +=|
|00005980| 20 73 74 65 70 2c 20 72 | 65 76 65 72 73 65 29 3b | step, r|everse);|
|00005990| 0a 58 09 09 09 09 09 69 | 66 20 28 70 69 78 65 6c |.X.....i|f (pixel|
|000059a0| 73 5f 6c 65 66 74 20 3e | 20 32 29 0a 58 09 09 09 |s_left >| 2).X...|
|000059b0| 09 09 09 70 6c 6f 74 28 | 2d 2d 78 31 2c 20 79 31 |...plot(|--x1, y1|
|000059c0| 2c 20 72 65 76 65 72 73 | 65 29 3b 0a 58 09 09 09 |, revers|e);.X...|
|000059d0| 09 7d 20 65 6c 73 65 20 | 7b 0a 58 09 09 09 09 20 |.} else |{.X.... |
|000059e0| 20 2f 2a 20 70 61 74 74 | 65 72 6e 20 33 20 2a 2f | /* patt|ern 3 */|
|000059f0| 0a 58 09 09 09 09 09 70 | 6c 6f 74 28 2b 2b 78 2c |.X.....p|lot(++x,|
|00005a00| 20 79 20 2b 3d 20 73 74 | 65 70 2c 20 72 65 76 65 | y += st|ep, reve|
|00005a10| 72 73 65 29 3b 0a 58 09 | 09 09 09 09 69 66 20 28 |rse);.X.|....if (|
|00005a20| 70 69 78 65 6c 73 5f 6c | 65 66 74 20 3e 20 31 29 |pixels_l|eft > 1)|
|00005a30| 0a 58 09 09 09 09 09 09 | 70 6c 6f 74 28 2b 2b 78 |.X......|plot(++x|
|00005a40| 2c 20 79 2c 20 72 65 76 | 65 72 73 65 29 3b 0a 58 |, y, rev|erse);.X|
|00005a50| 09 09 09 09 09 69 66 20 | 28 70 69 78 65 6c 73 5f |.....if |(pixels_|
|00005a60| 6c 65 66 74 20 3e 20 32 | 29 0a 58 09 09 09 09 09 |left > 2|).X.....|
|00005a70| 09 70 6c 6f 74 28 2d 2d | 78 31 2c 20 79 31 20 2d |.plot(--|x1, y1 -|
|00005a80| 3d 20 73 74 65 70 2c 20 | 72 65 76 65 72 73 65 29 |= step, |reverse)|
|00005a90| 3b 0a 58 09 09 09 09 7d | 0a 58 09 09 09 7d 0a 58 |;.X....}|.X...}.X|
|00005aa0| 09 09 7d 09 09 2f 2a 20 | 65 6e 64 20 69 66 20 70 |..}../* |end if p|
|00005ab0| 69 78 65 6c 73 5f 6c 65 | 66 74 20 2a 2f 0a 58 09 |ixels_le|ft */.X.|
|00005ac0| 7d 0a 58 09 2f 2a 20 65 | 6e 64 20 73 6c 6f 70 65 |}.X./* e|nd slope|
|00005ad0| 20 3c 20 31 2f 32 20 2a | 2f 0a 58 09 65 6c 73 65 | < 1/2 *|/.X.else|
|00005ae0| 20 7b 09 09 09 2f 2a 20 | 73 6c 6f 70 65 20 67 72 | {.../* |slope gr|
|00005af0| 65 61 74 65 72 20 74 68 | 61 6e 20 31 2f 32 20 2a |eater th|an 1/2 *|
|00005b00| 2f 0a 58 09 09 63 20 3d | 20 32 20 2a 20 28 64 79 |/.X..c =| 2 * (dy|
|00005b10| 20 2d 20 64 78 29 3b 0a | 58 09 09 69 6e 63 72 31 | - dx);.|X..incr1|
|00005b20| 20 3d 20 32 20 2a 20 63 | 3b 0a 58 09 09 44 20 3d | = 2 * c|;.X..D =|
|00005b30| 20 69 6e 63 72 31 20 2b | 20 64 78 3b 0a 58 09 09 | incr1 +| dx;.X..|
|00005b40| 66 6f 72 20 28 69 20 3d | 20 30 3b 20 69 20 3c 20 |for (i =| 0; i < |
|00005b50| 78 65 6e 64 3b 20 69 2b | 2b 29 20 7b 0a 58 09 09 |xend; i+|+) {.X..|
|00005b60| 09 2b 2b 78 3b 0a 58 09 | 09 09 2d 2d 78 31 3b 0a |.++x;.X.|..--x1;.|
|00005b70| 58 09 09 09 69 66 20 28 | 44 20 3e 20 30 29 20 7b |X...if (|D > 0) {|
|00005b80| 0a 58 09 09 09 20 20 2f | 2a 20 70 61 74 74 65 72 |.X... /|* patter|
|00005b90| 6e 20 34 20 66 6f 72 77 | 61 72 64 73 20 2a 2f 0a |n 4 forw|ards */.|
|00005ba0| 58 09 09 09 09 70 6c 6f | 74 28 78 2c 20 79 20 2b |X....plo|t(x, y +|
|00005bb0| 3d 20 73 74 65 70 2c 20 | 72 65 76 65 72 73 65 29 |= step, |reverse)|
|00005bc0| 3b 0a 58 09 09 09 09 70 | 6c 6f 74 28 2b 2b 78 2c |;.X....p|lot(++x,|
|00005bd0| 20 79 20 2b 3d 20 73 74 | 65 70 2c 20 72 65 76 65 | y += st|ep, reve|
|00005be0| 72 73 65 29 3b 0a 58 09 | 09 09 20 20 2f 2a 20 70 |rse);.X.|.. /* p|
|00005bf0| 61 74 74 65 72 6e 20 34 | 20 62 61 63 6b 77 61 72 |attern 4| backwar|
|00005c00| 64 73 20 2a 2f 0a 58 09 | 09 09 09 70 6c 6f 74 28 |ds */.X.|...plot(|
|00005c10| 78 31 2c 20 79 31 20 2d | 3d 20 73 74 65 70 2c 20 |x1, y1 -|= step, |
|00005c20| 72 65 76 65 72 73 65 29 | 3b 0a 58 09 09 09 09 70 |reverse)|;.X....p|
|00005c30| 6c 6f 74 28 2d 2d 78 31 | 2c 20 79 31 20 2d 3d 20 |lot(--x1|, y1 -= |
|00005c40| 73 74 65 70 2c 20 72 65 | 76 65 72 73 65 29 3b 0a |step, re|verse);.|
|00005c50| 58 09 09 09 09 44 20 2b | 3d 20 69 6e 63 72 31 3b |X....D +|= incr1;|
|00005c60| 0a 58 09 09 09 7d 20 65 | 6c 73 65 20 7b 0a 58 09 |.X...} e|lse {.X.|
|00005c70| 09 09 09 69 66 20 28 44 | 20 3c 20 63 29 20 7b 0a |...if (D| < c) {.|
|00005c80| 58 09 09 09 09 20 20 2f | 2a 20 70 61 74 74 65 72 |X.... /|* patter|
|00005c90| 6e 20 32 20 66 6f 72 77 | 61 72 64 73 20 2a 2f 0a |n 2 forw|ards */.|
|00005ca0| 58 09 09 09 09 09 70 6c | 6f 74 28 78 2c 20 79 2c |X.....pl|ot(x, y,|
|00005cb0| 20 72 65 76 65 72 73 65 | 29 3b 0a 58 09 09 09 09 | reverse|);.X....|
|00005cc0| 09 70 6c 6f 74 28 2b 2b | 78 2c 20 79 20 2b 3d 20 |.plot(++|x, y += |
|00005cd0| 73 74 65 70 2c 20 72 65 | 76 65 72 73 65 29 3b 0a |step, re|verse);.|
|00005ce0| 58 0a 58 20 09 09 09 09 | 20 20 2f 2a 20 70 61 74 |X.X ....| /* pat|
|00005cf0| 74 65 72 6e 20 32 20 62 | 61 63 6b 77 61 72 64 73 |tern 2 b|ackwards|
|00005d00| 20 2a 2f 0a 58 09 09 09 | 09 09 70 6c 6f 74 28 78 | */.X...|..plot(x|
|00005d10| 31 2c 20 79 31 2c 20 72 | 65 76 65 72 73 65 29 3b |1, y1, r|everse);|
|00005d20| 0a 58 09 09 09 09 09 70 | 6c 6f 74 28 2d 2d 78 31 |.X.....p|lot(--x1|
|00005d30| 2c 20 79 31 20 2d 3d 20 | 73 74 65 70 2c 20 72 65 |, y1 -= |step, re|
|00005d40| 76 65 72 73 65 29 3b 0a | 58 09 09 09 09 7d 20 65 |verse);.|X....} e|
|00005d50| 6c 73 65 20 7b 0a 58 09 | 09 09 09 20 20 2f 2a 20 |lse {.X.|... /* |
|00005d60| 70 61 74 74 65 72 6e 20 | 33 20 66 6f 72 77 61 72 |pattern |3 forwar|
|00005d70| 64 73 20 2a 2f 0a 58 09 | 09 09 09 09 70 6c 6f 74 |ds */.X.|....plot|
|00005d80| 28 78 2c 20 79 20 2b 3d | 20 73 74 65 70 2c 20 72 |(x, y +=| step, r|
|00005d90| 65 76 65 72 73 65 29 3b | 0a 58 09 09 09 09 09 70 |everse);|.X.....p|
|00005da0| 6c 6f 74 28 2b 2b 78 2c | 20 79 2c 20 72 65 76 65 |lot(++x,| y, reve|
|00005db0| 72 73 65 29 3b 0a 58 09 | 09 09 09 20 20 2f 2a 20 |rse);.X.|... /* |
|00005dc0| 70 61 74 74 65 72 6e 20 | 33 20 62 61 63 6b 77 61 |pattern |3 backwa|
|00005dd0| 72 64 73 20 2a 2f 0a 58 | 09 09 09 09 09 70 6c 6f |rds */.X|.....plo|
|00005de0| 74 28 78 31 2c 20 79 31 | 20 2d 3d 20 73 74 65 70 |t(x1, y1| -= step|
|00005df0| 2c 20 72 65 76 65 72 73 | 65 29 3b 0a 58 09 09 09 |, revers|e);.X...|
|00005e00| 09 09 70 6c 6f 74 28 2d | 2d 78 31 2c 20 79 31 2c |..plot(-|-x1, y1,|
|00005e10| 20 72 65 76 65 72 73 65 | 29 3b 0a 58 09 09 09 09 | reverse|);.X....|
|00005e20| 7d 0a 58 09 09 09 09 44 | 20 2b 3d 20 69 6e 63 72 |}.X....D| += incr|
|00005e30| 32 3b 0a 58 09 09 09 7d | 0a 58 09 09 7d 09 09 2f |2;.X...}|.X..}../|
|00005e40| 2a 20 65 6e 64 20 66 6f | 72 20 2a 2f 0a 58 09 09 |* end fo|r */.X..|
|00005e50| 2f 2a 20 70 6c 6f 74 20 | 6c 61 73 74 20 70 61 74 |/* plot |last pat|
|00005e60| 74 65 72 6e 20 2a 2f 0a | 58 09 09 69 66 20 28 70 |tern */.|X..if (p|
|00005e70| 69 78 65 6c 73 5f 6c 65 | 66 74 29 20 7b 0a 58 09 |ixels_le|ft) {.X.|
|00005e80| 09 09 69 66 20 28 44 20 | 3e 20 30 29 20 7b 0a 58 |..if (D |> 0) {.X|
|00005e90| 09 09 09 09 70 6c 6f 74 | 28 2b 2b 78 2c 20 79 20 |....plot|(++x, y |
|00005ea0| 2b 3d 20 73 74 65 70 2c | 20 72 65 76 65 72 73 65 |+= step,| reverse|
|00005eb0| 29 3b 09 2f 2a 20 70 61 | 74 74 65 72 6e 20 34 20 |);./* pa|ttern 4 |
|00005ec0| 2a 2f 0a 58 09 09 09 09 | 69 66 20 28 70 69 78 65 |*/.X....|if (pixe|
|00005ed0| 6c 73 5f 6c 65 66 74 20 | 3e 20 31 29 0a 58 09 09 |ls_left |> 1).X..|
|00005ee0| 09 09 09 70 6c 6f 74 28 | 2b 2b 78 2c 20 79 20 2b |...plot(|++x, y +|
|00005ef0| 3d 20 73 74 65 70 2c 20 | 72 65 76 65 72 73 65 29 |= step, |reverse)|
|00005f00| 3b 0a 58 09 09 09 09 69 | 66 20 28 70 69 78 65 6c |;.X....i|f (pixel|
|00005f10| 73 5f 6c 65 66 74 20 3e | 20 32 29 0a 58 09 09 09 |s_left >| 2).X...|
|00005f20| 09 09 70 6c 6f 74 28 2d | 2d 78 31 2c 20 79 31 20 |..plot(-|-x1, y1 |
|00005f30| 2d 3d 20 73 74 65 70 2c | 20 72 65 76 65 72 73 65 |-= step,| reverse|
|00005f40| 29 3b 0a 58 09 09 09 7d | 20 65 6c 73 65 20 7b 0a |);.X...}| else {.|
|00005f50| 58 09 09 09 09 69 66 20 | 28 44 20 3c 20 63 29 20 |X....if |(D < c) |
|00005f60| 7b 0a 58 09 09 09 09 09 | 70 6c 6f 74 28 2b 2b 78 |{.X.....|plot(++x|
|00005f70| 2c 20 79 2c 20 72 65 76 | 65 72 73 65 29 3b 09 2f |, y, rev|erse);./|
|00005f80| 2a 20 70 61 74 74 65 72 | 6e 20 32 20 20 2a 2f 0a |* patter|n 2 */.|
|00005f90| 58 09 09 09 09 09 69 66 | 20 28 70 69 78 65 6c 73 |X.....if| (pixels|
|00005fa0| 5f 6c 65 66 74 20 3e 20 | 31 29 0a 58 09 09 09 09 |_left > |1).X....|
|00005fb0| 09 09 70 6c 6f 74 28 2b | 2b 78 2c 20 79 20 2b 3d |..plot(+|+x, y +=|
|00005fc0| 20 73 74 65 70 2c 20 72 | 65 76 65 72 73 65 29 3b | step, r|everse);|
|00005fd0| 0a 58 09 09 09 09 09 69 | 66 20 28 70 69 78 65 6c |.X.....i|f (pixel|
|00005fe0| 73 5f 6c 65 66 74 20 3e | 20 32 29 0a 58 09 09 09 |s_left >| 2).X...|
|00005ff0| 09 09 09 70 6c 6f 74 28 | 2d 2d 78 31 2c 20 79 31 |...plot(|--x1, y1|
|00006000| 2c 20 72 65 76 65 72 73 | 65 29 3b 0a 58 09 09 09 |, revers|e);.X...|
|00006010| 09 7d 20 65 6c 73 65 20 | 7b 0a 58 09 09 09 09 20 |.} else |{.X.... |
|00006020| 20 2f 2a 20 70 61 74 74 | 65 72 6e 20 33 20 2a 2f | /* patt|ern 3 */|
|00006030| 0a 58 09 09 09 09 09 70 | 6c 6f 74 28 2b 2b 78 2c |.X.....p|lot(++x,|
|00006040| 20 79 20 2b 3d 20 73 74 | 65 70 2c 20 72 65 76 65 | y += st|ep, reve|
|00006050| 72 73 65 29 3b 0a 58 09 | 09 09 09 09 69 66 20 28 |rse);.X.|....if (|
|00006060| 70 69 78 65 6c 73 5f 6c | 65 66 74 20 3e 20 31 29 |pixels_l|eft > 1)|
|00006070| 0a 58 09 09 09 09 09 09 | 70 6c 6f 74 28 2b 2b 78 |.X......|plot(++x|
|00006080| 2c 20 79 2c 20 72 65 76 | 65 72 73 65 29 3b 0a 58 |, y, rev|erse);.X|
|00006090| 09 09 09 09 09 69 66 20 | 28 70 69 78 65 6c 73 5f |.....if |(pixels_|
|000060a0| 6c 65 66 74 20 3e 20 32 | 29 20 7b 0a 58 09 09 09 |left > 2|) {.X...|
|000060b0| 09 09 09 69 66 20 28 44 | 20 3e 20 63 29 20 2f 2a |...if (D| > c) /*|
|000060c0| 20 73 74 65 70 20 33 20 | 2a 2f 0a 58 09 09 09 09 | step 3 |*/.X....|
|000060d0| 09 09 20 20 20 70 6c 6f | 74 28 2d 2d 78 31 2c 20 |.. plo|t(--x1, |
|000060e0| 79 31 20 2d 3d 20 73 74 | 65 70 2c 20 72 65 76 65 |y1 -= st|ep, reve|
|000060f0| 72 73 65 29 3b 0a 58 09 | 09 09 09 09 09 65 6c 73 |rse);.X.|.....els|
|00006100| 65 20 2f 2a 20 73 74 65 | 70 20 32 20 2a 2f 0a 58 |e /* ste|p 2 */.X|
|00006110| 09 09 09 09 09 09 09 70 | 6c 6f 74 28 2d 2d 78 31 |.......p|lot(--x1|
|00006120| 2c 20 79 31 2c 20 72 65 | 76 65 72 73 65 29 3b 0a |, y1, re|verse);.|
|00006130| 58 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X | |
|00006140| 20 20 20 20 20 20 20 20 | 20 20 09 09 7d 0a 58 09 | | ..}.X.|
|00006150| 09 09 09 7d 0a 58 09 09 | 09 7d 0a 58 09 09 7d 0a |...}.X..|.}.X..}.|
|00006160| 58 09 7d 0a 58 7d 0a 58 | 2f 2a 20 6e 6f 6e 2d 7a |X.}.X}.X|/* non-z|
|00006170| 65 72 6f 20 66 6c 61 67 | 20 69 6e 64 69 63 61 74 |ero flag| indicat|
|00006180| 65 73 20 74 68 65 20 70 | 69 78 65 6c 73 20 6e 65 |es the p|ixels ne|
|00006190| 65 64 69 6e 67 20 73 77 | 61 70 20 62 61 63 6b 2e |eding sw|ap back.|
|000061a0| 20 2a 2f 0a 58 70 6c 6f | 74 28 78 2c 20 79 2c 20 | */.Xplo|t(x, y, |
|000061b0| 66 6c 61 67 29 20 69 6e | 74 20 78 2c 20 79 2c 20 |flag) in|t x, y, |
|000061c0| 66 6c 61 67 3b 0a 58 7b | 0a 58 09 69 66 20 28 66 |flag;.X{|.X.if (f|
|000061d0| 6c 61 67 29 0a 58 09 09 | 73 65 74 70 69 78 65 6c |lag).X..|setpixel|
|000061e0| 28 79 2c 20 78 29 3b 0a | 58 09 65 6c 73 65 0a 58 |(y, x);.|X.else.X|
|000061f0| 09 09 73 65 74 70 69 78 | 65 6c 28 78 2c 20 79 29 |..setpix|el(x, y)|
|00006200| 3b 0a 58 7d 0a 58 0a 45 | 4e 44 5f 4f 46 5f 46 49 |;.X}.X.E|ND_OF_FI|
|00006210| 4c 45 0a 69 66 20 74 65 | 73 74 20 34 36 34 37 20 |LE.if te|st 4647 |
|00006220| 2d 6e 65 20 60 77 63 20 | 2d 63 20 3c 27 44 6f 75 |-ne `wc |-c <'Dou|
|00006230| 62 6c 65 4c 69 6e 65 2e | 63 27 60 3b 20 74 68 65 |bleLine.|c'`; the|
|00006240| 6e 0a 20 20 20 20 65 63 | 68 6f 20 73 68 61 72 3a |n. ec|ho shar:|
|00006250| 20 5c 22 27 44 6f 75 62 | 6c 65 4c 69 6e 65 2e 63 | \"'Doub|leLine.c|
|00006260| 27 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |'\" unpa|cked wit|
|00006270| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 66 69 |h wrong |size!.fi|
|00006280| 0a 23 20 65 6e 64 20 6f | 66 20 27 44 6f 75 62 6c |.# end o|f 'Doubl|
|00006290| 65 4c 69 6e 65 2e 63 27 | 0a 66 69 0a 69 66 20 74 |eLine.c'|.fi.if t|
|000062a0| 65 73 74 20 2d 66 20 27 | 47 72 61 70 68 69 63 73 |est -f '|Graphics|
|000062b0| 47 65 6d 73 2e 68 27 20 | 2d 61 20 22 24 7b 31 7d |Gems.h' |-a "${1}|
|000062c0| 22 20 21 3d 20 22 2d 63 | 22 20 3b 20 74 68 65 6e |" != "-c|" ; then|
|000062d0| 20 0a 20 20 65 63 68 6f | 20 73 68 61 72 3a 20 57 | . echo| shar: W|
|000062e0| 69 6c 6c 20 6e 6f 74 20 | 63 6c 6f 62 62 65 72 20 |ill not |clobber |
|000062f0| 65 78 69 73 74 69 6e 67 | 20 66 69 6c 65 20 5c 22 |existing| file \"|
|00006300| 27 47 72 61 70 68 69 63 | 73 47 65 6d 73 2e 68 27 |'Graphic|sGems.h'|
|00006310| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|00006320| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|00006330| 27 47 72 61 70 68 69 63 | 73 47 65 6d 73 2e 68 27 |'Graphic|sGems.h'|
|00006340| 5c 22 20 5c 28 34 30 34 | 39 20 63 68 61 72 61 63 |\" \(404|9 charac|
|00006350| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00006360| 2f 2f 22 20 3e 27 47 72 | 61 70 68 69 63 73 47 65 |//" >'Gr|aphicsGe|
|00006370| 6d 73 2e 68 27 20 3c 3c | 27 45 4e 44 5f 4f 46 5f |ms.h' <<|'END_OF_|
|00006380| 46 49 4c 45 27 0a 58 2f | 2a 20 0a 58 20 2a 20 47 |FILE'.X/|* .X * G|
|00006390| 72 61 70 68 69 63 73 47 | 65 6d 73 2e 68 20 20 0a |raphicsG|ems.h .|
|000063a0| 58 20 2a 20 56 65 72 73 | 69 6f 6e 20 31 2e 30 20 |X * Vers|ion 1.0 |
|000063b0| 2d 20 41 6e 64 72 65 77 | 20 47 6c 61 73 73 6e 65 |- Andrew| Glassne|
|000063c0| 72 0a 58 20 2a 20 66 72 | 6f 6d 20 22 47 72 61 70 |r.X * fr|om "Grap|
|000063d0| 68 69 63 73 20 47 65 6d | 73 22 2c 20 41 63 61 64 |hics Gem|s", Acad|
|000063e0| 65 6d 69 63 20 50 72 65 | 73 73 2c 20 31 39 39 30 |emic Pre|ss, 1990|
|000063f0| 0a 58 20 2a 2f 0a 58 0a | 58 23 69 66 6e 64 65 66 |.X */.X.|X#ifndef|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
