home *** CD-ROM | disk | FTP | other *** search
/ Source Code 1992 March / Source_Code_CD-ROM_Walnut_Creek_March_1992.iso / usenet / compsrcs / unix / volume10 / complex.lib < prev    next >
SHell self-extracting ARchive  |  1988-09-11  |  24.2 KB
open in: MacOS 8.1     |     Win98     |     DOS

view JSON data     |     view as text

This file was processed as: SHell self-extracting ARchive (archive/shar).

You can browse this item here: complex.lib

ConfidenceProgramDetectionMatch TypeSupport
100% dexvert Newsgroup Content (archive/news) magic Supported
100% dexvert SHell self-extracting ARchive (archive/shar) magic Supported
1% dexvert Sibelius Music / Score (music/sibeliusMusicScore) ext Unsupported
1% dexvert IntelliFont Font (font/intelliFont) ext Unsupported
1% dexvert Amiga Hunk Library/Object (other/amigaLibrary) ext Unsupported
1% dexvert Amiga Shared Library (other/amigaSharedLibrary) ext Unsupported
1% dexvert Chess Assistant File (other/chessAssistant) ext Unsupported
1% dexvert COFF Library (other/coffLibrary) ext Unsupported
1% dexvert Maple Library (other/mapleLibrary) ext Unsupported
1% dexvert TopSpeed Modula-2 Library (other/topSpeedModula2Library) ext Unsupported
1% dexvert Text File (text/txt) fallback Supported
100% file news text default
99% file C source, ASCII text default
100% checkBytes Printable ASCII default
100% dexmagic PrintFox/Pagefox WEAK default
100% perlTextCheck Likely Text (Perl) default
100% siegfried fmt/329 Shell Archive Format default
100% detectItEasy Format: plain text[LF] default (weak)
100% xdgMime message/news default


hex view
+--------+-------------------------+-------------------------+--------+--------+
|00000000| 50 61 74 68 3a 20 75 75 | 6e 65 74 21 72 73 0a 46 |Path: uu|net!rs.F|
|00000010| 72 6f 6d 3a 20 72 73 40 | 75 75 6e 65 74 2e 55 55 |rom: rs@|uunet.UU|
|00000020| 2e 4e 45 54 20 28 52 69 | 63 68 20 53 61 6c 7a 29 |.NET (Ri|ch Salz)|
|00000030| 0a 4e 65 77 73 67 72 6f | 75 70 73 3a 20 63 6f 6d |.Newsgro|ups: com|
|00000040| 70 2e 73 6f 75 72 63 65 | 73 2e 75 6e 69 78 0a 53 |p.source|s.unix.S|
|00000050| 75 62 6a 65 63 74 3a 20 | 76 31 30 69 30 35 36 3a |ubject: |v10i056:|
|00000060| 20 20 43 6f 6d 70 6c 65 | 78 20 61 72 69 74 68 6d |  Comple|x arithm|
|00000070| 65 74 69 63 20 6c 69 62 | 72 61 72 79 0a 4d 65 73 |etic lib|rary.Mes|
|00000080| 73 61 67 65 2d 49 44 3a | 20 3c 36 34 37 40 75 75 |sage-ID:| <647@uu|
|00000090| 6e 65 74 2e 55 55 2e 4e | 45 54 3e 0a 44 61 74 65 |net.UU.N|ET>.Date|
|000000a0| 3a 20 31 36 20 4a 75 6c | 20 38 37 20 30 30 3a 31 |: 16 Jul| 87 00:1|
|000000b0| 32 3a 31 35 20 47 4d 54 | 0a 4f 72 67 61 6e 69 7a |2:15 GMT|.Organiz|
|000000c0| 61 74 69 6f 6e 3a 20 55 | 55 4e 45 54 20 43 6f 6d |ation: U|UNET Com|
|000000d0| 6d 75 6e 69 63 61 74 69 | 6f 6e 73 20 53 65 72 76 |municati|ons Serv|
|000000e0| 69 63 65 73 2c 20 41 72 | 6c 69 6e 67 74 6f 6e 2c |ices, Ar|lington,|
|000000f0| 20 56 41 0a 4c 69 6e 65 | 73 3a 20 31 30 39 35 0a | VA.Line|s: 1095.|
|00000100| 41 70 70 72 6f 76 65 64 | 3a 20 72 73 40 75 75 6e |Approved|: rs@uun|
|00000110| 65 74 2e 55 55 2e 4e 45 | 54 0a 0a 53 75 62 6d 69 |et.UU.NE|T..Submi|
|00000120| 74 74 65 64 2d 62 79 3a | 20 67 77 79 6e 40 62 72 |tted-by:| gwyn@br|
|00000130| 6c 2e 61 72 70 61 20 28 | 44 6f 75 67 20 47 77 79 |l.arpa (|Doug Gwy|
|00000140| 6e 29 0a 50 6f 73 74 69 | 6e 67 2d 4e 75 6d 62 65 |n).Posti|ng-Numbe|
|00000150| 72 3a 20 56 6f 6c 75 6d | 65 20 31 30 2c 20 49 73 |r: Volum|e 10, Is|
|00000160| 73 75 65 20 35 36 0a 41 | 72 63 68 69 76 65 2d 6e |sue 56.A|rchive-n|
|00000170| 61 6d 65 3a 20 63 6f 6d | 70 6c 65 78 2d 6c 69 62 |ame: com|plex-lib|
|00000180| 0a 0a 5b 20 20 48 65 72 | 65 20 69 73 20 61 20 70 |..[  Her|e is a p|
|00000190| 61 63 6b 61 67 65 20 74 | 6f 20 64 6f 20 63 6f 6d |ackage t|o do com|
|000001a0| 70 6c 65 78 20 61 72 69 | 74 68 6d 65 74 69 63 3b |plex ari|thmetic;|
|000001b0| 20 62 6c 61 6d 65 20 6d | 65 20 66 6f 72 20 74 68 | blame m|e for th|
|000001c0| 65 20 4d 61 6b 65 66 69 | 6c 65 2e 0a 09 2d 2d 72 |e Makefi|le...--r|
|000001d0| 24 20 20 5d 0a 0a 23 21 | 20 2f 62 69 6e 2f 73 68 |$  ]..#!| /bin/sh|
|000001e0| 0a 23 20 54 68 69 73 20 | 69 73 20 61 20 73 68 65 |.# This |is a she|
|000001f0| 6c 6c 20 61 72 63 68 69 | 76 65 2e 20 20 52 65 6d |ll archi|ve.  Rem|
|00000200| 6f 76 65 20 61 6e 79 74 | 68 69 6e 67 20 62 65 66 |ove anyt|hing bef|
|00000210| 6f 72 65 20 74 68 69 73 | 20 6c 69 6e 65 2c 20 74 |ore this| line, t|
|00000220| 68 65 6e 20 75 6e 70 61 | 63 6b 0a 23 20 69 74 20 |hen unpa|ck.# it |
|00000230| 62 79 20 73 61 76 69 6e | 67 20 69 74 20 69 6e 74 |by savin|g it int|
|00000240| 6f 20 61 20 66 69 6c 65 | 20 61 6e 64 20 74 79 70 |o a file| and typ|
|00000250| 69 6e 67 20 22 73 68 20 | 66 69 6c 65 22 2e 20 20 |ing "sh |file".  |
|00000260| 54 6f 20 6f 76 65 72 77 | 72 69 74 65 20 65 78 69 |To overw|rite exi|
|00000270| 73 74 69 6e 67 0a 23 20 | 66 69 6c 65 73 2c 20 74 |sting.# |files, t|
|00000280| 79 70 65 20 22 73 68 20 | 66 69 6c 65 20 2d 63 22 |ype "sh |file -c"|
|00000290| 2e 20 20 59 6f 75 20 63 | 61 6e 20 61 6c 73 6f 20 |.  You c|an also |
|000002a0| 66 65 65 64 20 74 68 69 | 73 20 61 73 20 73 74 61 |feed thi|s as sta|
|000002b0| 6e 64 61 72 64 20 69 6e | 70 75 74 20 76 69 61 0a |ndard in|put via.|
|000002c0| 23 20 75 6e 73 68 61 72 | 2c 20 6f 72 20 62 79 20 |# unshar|, or by |
|000002d0| 74 79 70 69 6e 67 20 22 | 73 68 20 3c 66 69 6c 65 |typing "|sh <file|
|000002e0| 22 2c 20 65 2e 67 2e 2e | 20 20 49 66 20 74 68 69 |", e.g..|  If thi|
|000002f0| 73 20 61 72 63 68 69 76 | 65 20 69 73 20 63 6f 6d |s archiv|e is com|
|00000300| 70 6c 65 74 65 2c 20 79 | 6f 75 0a 23 20 77 69 6c |plete, y|ou.# wil|
|00000310| 6c 20 73 65 65 20 74 68 | 65 20 66 6f 6c 6c 6f 77 |l see th|e follow|
|00000320| 69 6e 67 20 6d 65 73 73 | 61 67 65 20 61 74 20 74 |ing mess|age at t|
|00000330| 68 65 20 65 6e 64 3a 0a | 23 09 09 22 45 6e 64 20 |he end:.|#.."End |
|00000340| 6f 66 20 73 68 65 6c 6c | 20 61 72 63 68 69 76 65 |of shell| archive|
|00000350| 2e 22 0a 23 20 43 6f 6e | 74 65 6e 74 73 3a 20 20 |.".# Con|tents:  |
|00000360| 4d 61 6b 65 66 69 6c 65 | 20 63 6f 6d 70 6c 65 78 |Makefile| complex|
|00000370| 2e 33 20 63 6f 6d 70 6c | 65 78 2e 68 20 63 78 5f |.3 compl|ex.h cx_|
|00000380| 74 65 73 74 2e 63 20 63 | 78 61 64 64 2e 63 20 63 |test.c c|xadd.c c|
|00000390| 78 61 6d 70 6c 2e 63 0a | 23 20 20 20 63 78 63 6f |xampl.c.|#   cxco|
|000003a0| 6e 6a 2e 63 20 63 78 63 | 6f 6e 73 2e 63 20 63 78 |nj.c cxc|ons.c cx|
|000003b0| 63 6f 70 79 2e 63 20 63 | 78 64 69 76 2e 63 20 63 |copy.c c|xdiv.c c|
|000003c0| 78 6d 75 6c 2e 63 20 63 | 78 70 68 61 73 2e 63 20 |xmul.c c|xphas.c |
|000003d0| 63 78 70 68 73 72 2e 63 | 0a 23 20 20 20 63 78 73 |cxphsr.c|.#   cxs|
|000003e0| 63 61 6c 2e 63 20 63 78 | 73 71 72 74 2e 63 20 63 |cal.c cx|sqrt.c c|
|000003f0| 78 73 75 62 2e 63 0a 50 | 41 54 48 3d 2f 62 69 6e |xsub.c.P|ATH=/bin|
|00000400| 3a 2f 75 73 72 2f 62 69 | 6e 3a 2f 75 73 72 2f 75 |:/usr/bi|n:/usr/u|
|00000410| 63 62 20 3b 20 65 78 70 | 6f 72 74 20 50 41 54 48 |cb ; exp|ort PATH|
|00000420| 0a 69 66 20 74 65 73 74 | 20 2d 66 20 4d 61 6b 65 |.if test| -f Make|
|00000430| 66 69 6c 65 20 2d 61 20 | 22 24 7b 31 7d 22 20 21 |file -a |"${1}" !|
|00000440| 3d 20 22 2d 63 22 20 3b | 20 74 68 65 6e 20 0a 20 |= "-c" ;| then . |
|00000450| 20 65 63 68 6f 20 73 68 | 61 72 3a 20 57 69 6c 6c | echo sh|ar: Will|
|00000460| 20 6e 6f 74 20 6f 76 65 | 72 2d 77 72 69 74 65 20 | not ove|r-write |
|00000470| 65 78 69 73 74 69 6e 67 | 20 66 69 6c 65 20 5c 22 |existing| file \"|
|00000480| 4d 61 6b 65 66 69 6c 65 | 5c 22 0a 65 6c 73 65 0a |Makefile|\".else.|
|00000490| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|000004a0| 63 74 69 6e 67 20 5c 22 | 4d 61 6b 65 66 69 6c 65 |cting \"|Makefile|
|000004b0| 5c 22 20 5c 28 33 39 35 | 20 63 68 61 72 61 63 74 |\" \(395| charact|
|000004c0| 65 72 73 5c 29 0a 73 65 | 64 20 22 73 2f 5e 58 2f |ers\).se|d "s/^X/|
|000004d0| 2f 22 20 3e 4d 61 6b 65 | 66 69 6c 65 20 3c 3c 27 |/" >Make|file <<'|
|000004e0| 45 4e 44 5f 4f 46 5f 4d | 61 6b 65 66 69 6c 65 27 |END_OF_M|akefile'|
|000004f0| 0a 58 41 4c 4c 3d 63 6f | 6d 70 6c 65 78 2e 33 20 |.XALL=co|mplex.3 |
|00000500| 63 6f 6d 70 6c 65 78 2e | 68 20 6c 69 62 63 6f 6d |complex.|h libcom|
|00000510| 70 6c 65 78 2e 61 0a 58 | 4f 42 4a 53 3d 5c 0a 58 |plex.a.X|OBJS=\.X|
|00000520| 20 20 20 20 63 78 61 64 | 64 2e 6f 20 63 78 61 6d |    cxad|d.o cxam|
|00000530| 70 6c 2e 6f 20 63 78 63 | 6f 6e 6a 2e 6f 20 63 78 |pl.o cxc|onj.o cx|
|00000540| 63 6f 6e 73 2e 6f 20 63 | 78 63 6f 70 79 2e 6f 20 |cons.o c|xcopy.o |
|00000550| 63 78 64 69 76 2e 6f 20 | 63 78 6d 75 6c 2e 6f 20 |cxdiv.o |cxmul.o |
|00000560| 5c 0a 58 20 20 20 20 63 | 78 70 68 61 73 2e 6f 20 |\.X    c|xphas.o |
|00000570| 63 78 70 68 73 72 2e 6f | 20 63 78 73 63 61 6c 2e |cxphsr.o| cxscal.|
|00000580| 6f 20 63 78 73 71 72 74 | 2e 6f 20 63 78 73 75 62 |o cxsqrt|.o cxsub|
|00000590| 2e 6f 0a 58 0a 58 61 6c | 6c 3a 09 24 28 41 4c 4c |.o.X.Xal|l:.$(ALL|
|000005a0| 29 0a 58 0a 58 69 6e 73 | 74 61 6c 6c 3a 09 24 28 |).X.Xins|tall:.$(|
|000005b0| 41 4c 4c 29 0a 58 09 40 | 65 63 68 6f 20 69 6e 73 |ALL).X.@|echo ins|
|000005c0| 74 61 6c 6c 20 24 28 41 | 4c 4c 29 20 61 63 63 6f |tall $(A|LL) acco|
|000005d0| 72 64 69 6e 67 20 74 6f | 20 6c 6f 63 61 6c 20 63 |rding to| local c|
|000005e0| 6f 6e 76 65 6e 74 69 6f | 6e 2e 0a 58 0a 58 63 78 |onventio|n..X.Xcx|
|000005f0| 5f 74 65 73 74 3a 09 63 | 78 5f 74 65 73 74 2e 63 |_test:.c|x_test.c|
|00000600| 20 6c 69 62 63 6f 6d 70 | 6c 65 78 2e 61 0a 58 09 | libcomp|lex.a.X.|
|00000610| 24 28 43 43 29 20 24 28 | 43 46 4c 41 47 53 29 20 |$(CC) $(|CFLAGS) |
|00000620| 2d 6f 20 63 78 5f 74 65 | 73 74 20 63 78 5f 74 65 |-o cx_te|st cx_te|
|00000630| 73 74 2e 63 20 6c 69 62 | 63 6f 6d 70 6c 65 78 2e |st.c lib|complex.|
|00000640| 61 0a 58 0a 58 6c 69 62 | 63 6f 6d 70 6c 65 78 2e |a.X.Xlib|complex.|
|00000650| 61 3a 09 24 28 4f 42 4a | 53 29 0a 58 09 61 72 20 |a:.$(OBJ|S).X.ar |
|00000660| 72 20 6c 69 62 63 6f 6d | 70 6c 65 78 2e 61 20 24 |r libcom|plex.a $|
|00000670| 28 4f 42 4a 53 29 0a 58 | 0a 58 24 28 4f 42 4a 53 |(OBJS).X|.X$(OBJS|
|00000680| 29 3a 09 63 6f 6d 70 6c | 65 78 2e 68 0a 45 4e 44 |):.compl|ex.h.END|
|00000690| 5f 4f 46 5f 4d 61 6b 65 | 66 69 6c 65 0a 69 66 20 |_OF_Make|file.if |
|000006a0| 74 65 73 74 20 33 39 35 | 20 2d 6e 65 20 60 77 63 |test 395| -ne `wc|
|000006b0| 20 2d 63 20 3c 4d 61 6b | 65 66 69 6c 65 60 3b 20 | -c <Mak|efile`; |
|000006c0| 74 68 65 6e 0a 20 20 20 | 20 65 63 68 6f 20 73 68 |then.   | echo sh|
|000006d0| 61 72 3a 20 5c 22 4d 61 | 6b 65 66 69 6c 65 5c 22 |ar: \"Ma|kefile\"|
|000006e0| 20 75 6e 70 61 63 6b 65 | 64 20 77 69 74 68 20 77 | unpacke|d with w|
|000006f0| 72 6f 6e 67 20 73 69 7a | 65 21 0a 66 69 0a 23 20 |rong siz|e!.fi.# |
|00000700| 65 6e 64 20 6f 66 20 6f | 76 65 72 77 72 69 74 69 |end of o|verwriti|
|00000710| 6e 67 20 63 68 65 63 6b | 0a 66 69 0a 69 66 20 74 |ng check|.fi.if t|
|00000720| 65 73 74 20 2d 66 20 63 | 6f 6d 70 6c 65 78 2e 33 |est -f c|omplex.3|
|00000730| 20 2d 61 20 22 24 7b 31 | 7d 22 20 21 3d 20 22 2d | -a "${1|}" != "-|
|00000740| 63 22 20 3b 20 74 68 65 | 6e 20 0a 20 20 65 63 68 |c" ; the|n .  ech|
|00000750| 6f 20 73 68 61 72 3a 20 | 57 69 6c 6c 20 6e 6f 74 |o shar: |Will not|
|00000760| 20 6f 76 65 72 2d 77 72 | 69 74 65 20 65 78 69 73 | over-wr|ite exis|
|00000770| 74 69 6e 67 20 66 69 6c | 65 20 5c 22 63 6f 6d 70 |ting fil|e \"comp|
|00000780| 6c 65 78 2e 33 5c 22 0a | 65 6c 73 65 0a 65 63 68 |lex.3\".|else.ech|
|00000790| 6f 20 73 68 61 72 3a 20 | 45 78 74 72 61 63 74 69 |o shar: |Extracti|
|000007a0| 6e 67 20 5c 22 63 6f 6d | 70 6c 65 78 2e 33 5c 22 |ng \"com|plex.3\"|
|000007b0| 20 5c 28 35 38 33 30 20 | 63 68 61 72 61 63 74 65 | \(5830 |characte|
|000007c0| 72 73 5c 29 0a 73 65 64 | 20 22 73 2f 5e 58 2f 2f |rs\).sed| "s/^X//|
|000007d0| 22 20 3e 63 6f 6d 70 6c | 65 78 2e 33 20 3c 3c 27 |" >compl|ex.3 <<'|
|000007e0| 45 4e 44 5f 4f 46 5f 63 | 6f 6d 70 6c 65 78 2e 33 |END_OF_c|omplex.3|
|000007f0| 27 0a 58 27 5c 22 20 65 | 0a 58 2e 54 48 20 43 4f |'.X'\" e|.X.TH CO|
|00000800| 4d 50 4c 45 58 20 33 56 | 20 4c 4f 43 41 4c 0a 58 |MPLEX 3V| LOCAL.X|
|00000810| 27 5c 22 09 6c 61 73 74 | 20 65 64 69 74 3a 09 38 |'\".last| edit:.8|
|00000820| 36 2f 30 32 2f 30 33 09 | 44 20 41 20 47 77 79 6e |6/02/03.|D A Gwyn|
|00000830| 0a 58 27 5c 22 09 53 43 | 43 53 20 49 44 3a 09 40 |.X'\".SC|CS ID:.@|
|00000840| 28 23 29 63 6f 6d 70 6c | 65 78 2e 33 09 31 2e 32 |(#)compl|ex.3.1.2|
|00000850| 20 28 6d 6f 64 69 66 69 | 65 64 20 66 6f 72 20 70 | (modifi|ed for p|
|00000860| 75 62 6c 69 63 20 76 65 | 72 73 69 6f 6e 29 0a 58 |ublic ve|rsion).X|
|00000870| 2e 45 51 0a 58 64 65 6c | 69 6d 20 40 40 0a 58 2e |.EQ.Xdel|im @@.X.|
|00000880| 45 4e 0a 58 2e 53 48 20 | 4e 41 4d 45 0a 58 63 6f |EN.X.SH |NAME.Xco|
|00000890| 6d 70 6c 65 78 20 5c 2d | 20 63 6f 6d 70 6c 65 78 |mplex \-| complex|
|000008a0| 20 61 72 69 74 68 6d 65 | 74 69 63 20 6f 70 65 72 | arithme|tic oper|
|000008b0| 61 74 69 6f 6e 73 0a 58 | 2e 53 48 20 53 59 4e 4f |ations.X|.SH SYNO|
|000008c0| 50 53 49 53 0a 58 2e 42 | 0a 58 23 69 6e 63 6c 75 |PSIS.X.B|.X#inclu|
|000008d0| 64 65 20 3c 63 6f 6d 70 | 6c 65 78 2e 68 3e 09 2f |de <comp|lex.h>./|
|000008e0| 2a 20 61 73 73 75 6d 69 | 6e 67 20 61 70 70 72 6f |* assumi|ng appro|
|000008f0| 70 72 69 61 74 65 20 63 | 63 20 5c 2d 49 20 6f 70 |priate c|c \-I op|
|00000900| 74 69 6f 6e 20 2a 2f 0a | 58 2e 62 72 0a 58 2f 2a |tion */.|X.br.X/*|
|00000910| 20 41 6c 6c 20 74 68 65 | 20 66 6f 6c 6c 6f 77 69 | All the| followi|
|00000920| 6e 67 20 66 75 6e 63 74 | 69 6f 6e 73 20 61 72 65 |ng funct|ions are|
|00000930| 20 64 65 63 6c 61 72 65 | 64 20 69 6e 20 74 68 69 | declare|d in thi|
|00000940| 73 20 68 65 61 64 65 72 | 20 66 69 6c 65 2e 20 2a |s header| file. *|
|00000950| 2f 0a 58 2e 50 0a 58 2e | 42 20 63 6f 6d 70 6c 65 |/.X.P.X.|B comple|
|00000960| 78 20 2a 43 78 41 64 64 | 28 61 70 2c 62 70 29 3b |x *CxAdd|(ap,bp);|
|00000970| 0a 58 2e 62 72 0a 58 2e | 42 20 63 6f 6d 70 6c 65 |.X.br.X.|B comple|
|00000980| 78 20 2a 61 70 2c 20 2a | 62 70 3b 0a 58 2e 50 0a |x *ap, *|bp;.X.P.|
|00000990| 58 2e 42 20 63 6f 6d 70 | 6c 65 78 20 2a 43 78 53 |X.B comp|lex *CxS|
|000009a0| 75 62 28 61 70 2c 62 70 | 29 3b 0a 58 2e 62 72 0a |ub(ap,bp|);.X.br.|
|000009b0| 58 2e 42 20 63 6f 6d 70 | 6c 65 78 20 2a 61 70 2c |X.B comp|lex *ap,|
|000009c0| 20 2a 62 70 3b 0a 58 2e | 50 0a 58 2e 42 20 63 6f | *bp;.X.|P.X.B co|
|000009d0| 6d 70 6c 65 78 20 2a 43 | 78 4d 75 6c 28 61 70 2c |mplex *C|xMul(ap,|
|000009e0| 62 70 29 3b 0a 58 2e 62 | 72 0a 58 2e 42 20 63 6f |bp);.X.b|r.X.B co|
|000009f0| 6d 70 6c 65 78 20 2a 61 | 70 2c 20 2a 62 70 3b 0a |mplex *a|p, *bp;.|
|00000a00| 58 2e 50 0a 58 2e 42 20 | 63 6f 6d 70 6c 65 78 20 |X.P.X.B |complex |
|00000a10| 2a 43 78 44 69 76 28 61 | 70 2c 62 70 29 3b 0a 58 |*CxDiv(a|p,bp);.X|
|00000a20| 2e 62 72 0a 58 2e 42 20 | 63 6f 6d 70 6c 65 78 20 |.br.X.B |complex |
|00000a30| 2a 61 70 2c 20 2a 62 70 | 3b 0a 58 2e 50 0a 58 2e |*ap, *bp|;.X.P.X.|
|00000a40| 42 20 63 6f 6d 70 6c 65 | 78 20 2a 43 78 53 71 72 |B comple|x *CxSqr|
|00000a50| 74 28 63 70 29 3b 0a 58 | 2e 62 72 0a 58 2e 42 20 |t(cp);.X|.br.X.B |
|00000a60| 63 6f 6d 70 6c 65 78 20 | 2a 63 70 3b 0a 58 2e 50 |complex |*cp;.X.P|
|00000a70| 0a 58 2e 42 20 63 6f 6d | 70 6c 65 78 20 2a 43 78 |.X.B com|plex *Cx|
|00000a80| 53 63 61 6c 28 63 70 2c | 5c 5e 73 29 3b 0a 58 2e |Scal(cp,|\^s);.X.|
|00000a90| 62 72 0a 58 2e 42 20 63 | 6f 6d 70 6c 65 78 20 2a |br.X.B c|omplex *|
|00000aa0| 63 70 3b 0a 58 2e 62 72 | 0a 58 2e 42 20 64 6f 75 |cp;.X.br|.X.B dou|
|00000ab0| 62 6c 65 20 73 3b 0a 58 | 2e 50 0a 58 2e 42 20 63 |ble s;.X|.P.X.B c|
|00000ac0| 6f 6d 70 6c 65 78 20 2a | 43 78 4e 65 67 28 63 70 |omplex *|CxNeg(cp|
|00000ad0| 29 3b 0a 58 2e 62 72 0a | 58 2e 42 20 63 6f 6d 70 |);.X.br.|X.B comp|
|00000ae0| 6c 65 78 20 2a 63 70 3b | 0a 58 2e 50 0a 58 2e 42 |lex *cp;|.X.P.X.B|
|00000af0| 20 63 6f 6d 70 6c 65 78 | 20 2a 43 78 43 6f 6e 6a | complex| *CxConj|
|00000b00| 28 63 70 29 3b 0a 58 2e | 62 72 0a 58 2e 42 20 63 |(cp);.X.|br.X.B c|
|00000b10| 6f 6d 70 6c 65 78 20 2a | 63 70 3b 0a 58 2e 50 0a |omplex *|cp;.X.P.|
|00000b20| 58 2e 42 20 63 6f 6d 70 | 6c 65 78 20 2a 43 78 43 |X.B comp|lex *CxC|
|00000b30| 6f 70 79 28 61 70 2c 62 | 70 29 3b 0a 58 2e 62 72 |opy(ap,b|p);.X.br|
|00000b40| 0a 58 2e 42 20 63 6f 6d | 70 6c 65 78 20 2a 61 70 |.X.B com|plex *ap|
|00000b50| 2c 20 2a 62 70 3b 0a 58 | 2e 50 0a 58 2e 42 20 63 |, *bp;.X|.P.X.B c|
|00000b60| 6f 6d 70 6c 65 78 20 2a | 43 78 43 6f 6e 73 28 63 |omplex *|CxCons(c|
|00000b70| 70 2c 5c 5e 72 2c 5c 5e | 69 29 3b 0a 58 2e 62 72 |p,\^r,\^|i);.X.br|
|00000b80| 0a 58 2e 42 20 63 6f 6d | 70 6c 65 78 20 2a 63 70 |.X.B com|plex *cp|
|00000b90| 3b 0a 58 2e 62 72 0a 58 | 2e 42 20 64 6f 75 62 6c |;.X.br.X|.B doubl|
|00000ba0| 65 20 72 2c 20 69 3b 0a | 58 2e 50 0a 58 2e 42 20 |e r, i;.|X.P.X.B |
|00000bb0| 63 6f 6d 70 6c 65 78 20 | 2a 43 78 50 68 73 72 28 |complex |*CxPhsr(|
|00000bc0| 63 70 2c 6d 2c 70 29 3b | 0a 58 2e 62 72 0a 58 2e |cp,m,p);|.X.br.X.|
|00000bd0| 42 20 63 6f 6d 70 6c 65 | 78 20 2a 63 70 3b 0a 58 |B comple|x *cp;.X|
|00000be0| 2e 62 72 0a 58 2e 42 20 | 64 6f 75 62 6c 65 20 6d |.br.X.B |double m|
|00000bf0| 2c 20 70 3b 0a 58 2e 50 | 0a 58 2e 42 20 64 6f 75 |, p;.X.P|.X.B dou|
|00000c00| 62 6c 65 20 43 78 52 65 | 61 6c 28 63 70 29 3b 0a |ble CxRe|al(cp);.|
|00000c10| 58 2e 62 72 0a 58 2e 42 | 20 63 6f 6d 70 6c 65 78 |X.br.X.B| complex|
|00000c20| 20 2a 63 70 3b 0a 58 2e | 50 0a 58 2e 42 20 64 6f | *cp;.X.|P.X.B do|
|00000c30| 75 62 6c 65 20 43 78 49 | 6d 61 67 28 63 70 29 3b |uble CxI|mag(cp);|
|00000c40| 0a 58 2e 62 72 0a 58 2e | 42 20 63 6f 6d 70 6c 65 |.X.br.X.|B comple|
|00000c50| 78 20 2a 63 70 3b 0a 58 | 2e 50 0a 58 2e 42 20 64 |x *cp;.X|.P.X.B d|
|00000c60| 6f 75 62 6c 65 20 43 78 | 41 6d 70 6c 28 63 70 29 |ouble Cx|Ampl(cp)|
|00000c70| 3b 0a 58 2e 62 72 0a 58 | 2e 42 20 63 6f 6d 70 6c |;.X.br.X|.B compl|
|00000c80| 65 78 20 2a 63 70 3b 0a | 58 2e 50 0a 58 2e 42 20 |ex *cp;.|X.P.X.B |
|00000c90| 64 6f 75 62 6c 65 20 43 | 78 50 68 61 73 28 63 70 |double C|xPhas(cp|
|00000ca0| 29 3b 0a 58 2e 62 72 0a | 58 2e 42 20 63 6f 6d 70 |);.X.br.|X.B comp|
|00000cb0| 6c 65 78 20 2a 63 70 3b | 0a 58 2e 50 0a 58 2e 42 |lex *cp;|.X.P.X.B|
|00000cc0| 20 63 6f 6d 70 6c 65 78 | 20 2a 43 78 41 6c 6c 6f | complex| *CxAllo|
|00000cd0| 28 5c 20 29 3b 0a 58 2e | 50 0a 58 2e 42 20 76 6f |(\ );.X.|P.X.B vo|
|00000ce0| 69 64 20 43 78 46 72 65 | 65 28 63 70 29 3b 0a 58 |id CxFre|e(cp);.X|
|00000cf0| 2e 62 72 0a 58 2e 42 20 | 63 6f 6d 70 6c 65 78 20 |.br.X.B |complex |
|00000d00| 2a 63 70 3b 0a 58 2e 53 | 48 20 44 45 53 43 52 49 |*cp;.X.S|H DESCRI|
|00000d10| 50 54 49 4f 4e 0a 58 54 | 68 65 73 65 20 72 6f 75 |PTION.XT|hese rou|
|00000d20| 74 69 6e 65 73 20 70 65 | 72 66 6f 72 6d 20 61 72 |tines pe|rform ar|
|00000d30| 69 74 68 6d 65 74 69 63 | 0a 58 61 6e 64 20 6f 74 |ithmetic|.Xand ot|
|00000d40| 68 65 72 20 75 73 65 66 | 75 6c 20 6f 70 65 72 61 |her usef|ul opera|
|00000d50| 74 69 6f 6e 73 20 6f 6e | 20 63 6f 6d 70 6c 65 78 |tions on| complex|
|00000d60| 20 6e 75 6d 62 65 72 73 | 2e 0a 58 41 6e 20 61 70 | numbers|..XAn ap|
|00000d70| 70 72 6f 70 72 69 61 74 | 65 20 64 61 74 61 20 73 |propriat|e data s|
|00000d80| 74 72 75 63 74 75 72 65 | 0a 58 2e 42 20 63 6f 6d |tructure|.X.B com|
|00000d90| 70 6c 65 78 0a 58 69 73 | 20 64 65 66 69 6e 65 64 |plex.Xis| defined|
|00000da0| 20 69 6e 20 74 68 65 20 | 68 65 61 64 65 72 20 66 | in the |header f|
|00000db0| 69 6c 65 3b 0a 58 61 6c | 6c 20 61 63 63 65 73 73 |ile;.Xal|l access|
|00000dc0| 20 74 6f 0a 58 2e 42 20 | 63 6f 6d 70 6c 65 78 0a | to.X.B |complex.|
|00000dd0| 58 64 61 74 61 20 73 68 | 6f 75 6c 64 20 62 65 0a |Xdata sh|ould be.|
|00000de0| 58 2e 49 20 76 69 61 0a | 58 74 68 65 73 65 20 70 |X.I via.|Xthese p|
|00000df0| 72 65 64 65 66 69 6e 65 | 64 20 66 75 6e 63 74 69 |redefine|d functi|
|00000e00| 6f 6e 73 2e 0a 58 28 53 | 65 65 0a 58 2e 53 4d 20 |ons..X(S|ee.X.SM |
|00000e10| 48 49 4e 54 53 0a 58 66 | 6f 72 20 66 75 72 74 68 |HINTS.Xf|or furth|
|00000e20| 65 72 20 69 6e 66 6f 72 | 6d 61 74 69 6f 6e 2e 29 |er infor|mation.)|
|00000e30| 0a 58 2e 50 0a 58 49 6e | 20 74 68 65 20 66 6f 6c |.X.P.XIn| the fol|
|00000e40| 6c 6f 77 69 6e 67 20 64 | 65 73 63 72 69 70 74 69 |lowing d|escripti|
|00000e50| 6f 6e 73 2c 0a 58 74 68 | 65 20 6e 61 6d 65 73 0a |ons,.Xth|e names.|
|00000e60| 58 2e 49 52 20 61 20 2c | 0a 58 2e 49 52 20 62 20 |X.IR a ,|.X.IR b |
|00000e70| 2c 0a 58 61 6e 64 0a 58 | 2e 49 20 63 0a 58 72 65 |,.Xand.X|.I c.Xre|
|00000e80| 70 72 65 73 65 6e 74 20 | 74 68 65 0a 58 2e 42 20 |present |the.X.B |
|00000e90| 63 6f 6d 70 6c 65 78 0a | 58 64 61 74 61 20 61 64 |complex.|Xdata ad|
|00000ea0| 64 72 65 73 73 65 64 20 | 62 79 20 74 68 65 20 63 |dressed |by the c|
|00000eb0| 6f 72 72 65 73 70 6f 6e | 64 69 6e 67 20 70 6f 69 |orrespon|ding poi|
|00000ec0| 6e 74 65 72 73 0a 58 2e | 49 52 20 61 70 20 2c 0a |nters.X.|IR ap ,.|
|00000ed0| 58 2e 49 52 20 62 70 20 | 2c 0a 58 61 6e 64 0a 58 |X.IR bp |,.Xand.X|
|00000ee0| 2e 49 52 20 63 70 20 2e | 0a 58 2e 50 0a 58 2e 49 |.IR cp .|.X.P.X.I|
|00000ef0| 20 43 78 41 64 64 5c 5e | 0a 58 61 64 64 73 0a 58 | CxAdd\^|.Xadds.X|
|00000f00| 2e 49 20 62 0a 58 74 6f | 0a 58 2e 49 20 61 0a 58 |.I b.Xto|.X.I a.X|
|00000f10| 61 6e 64 20 72 65 74 75 | 72 6e 73 20 61 20 70 6f |and retu|rns a po|
|00000f20| 69 6e 74 65 72 20 74 6f | 20 74 68 65 20 72 65 73 |inter to| the res|
|00000f30| 75 6c 74 2e 0a 58 2e 50 | 0a 58 2e 49 20 43 78 53 |ult..X.P|.X.I CxS|
|00000f40| 75 62 0a 58 73 75 62 74 | 72 61 63 74 73 0a 58 2e |ub.Xsubt|racts.X.|
|00000f50| 49 20 62 0a 58 66 72 6f | 6d 0a 58 2e 49 20 61 0a |I b.Xfro|m.X.I a.|
|00000f60| 58 61 6e 64 20 72 65 74 | 75 72 6e 73 20 61 20 70 |Xand ret|urns a p|
|00000f70| 6f 69 6e 74 65 72 20 74 | 6f 20 74 68 65 20 72 65 |ointer t|o the re|
|00000f80| 73 75 6c 74 2e 0a 58 2e | 50 0a 58 2e 49 20 43 78 |sult..X.|P.X.I Cx|
|00000f90| 4d 75 6c 5c 5e 0a 58 6d | 75 6c 74 69 70 6c 69 65 |Mul\^.Xm|ultiplie|
|00000fa0| 73 0a 58 2e 49 20 61 0a | 58 62 79 0a 58 2e 49 20 |s.X.I a.|Xby.X.I |
|00000fb0| 62 0a 58 61 6e 64 20 72 | 65 74 75 72 6e 73 20 61 |b.Xand r|eturns a|
|00000fc0| 20 70 6f 69 6e 74 65 72 | 20 74 6f 20 74 68 65 20 | pointer| to the |
|00000fd0| 72 65 73 75 6c 74 2e 0a | 58 2e 50 0a 58 2e 49 20 |result..|X.P.X.I |
|00000fe0| 43 78 44 69 76 0a 58 64 | 69 76 69 64 65 73 0a 58 |CxDiv.Xd|ivides.X|
|00000ff0| 2e 49 20 61 0a 58 62 79 | 0a 58 2e 49 20 62 0a 58 |.I a.Xby|.X.I b.X|
|00001000| 61 6e 64 20 72 65 74 75 | 72 6e 73 20 61 20 70 6f |and retu|rns a po|
|00001010| 69 6e 74 65 72 20 74 6f | 20 74 68 65 20 72 65 73 |inter to| the res|
|00001020| 75 6c 74 2e 0a 58 54 68 | 65 20 64 69 76 69 73 6f |ult..XTh|e diviso|
|00001030| 72 20 6d 75 73 74 20 6e | 6f 74 20 62 65 20 70 72 |r must n|ot be pr|
|00001040| 65 63 69 73 65 6c 79 20 | 7a 65 72 6f 2e 0a 58 2e |ecisely |zero..X.|
|00001050| 50 0a 58 2e 49 20 43 78 | 53 71 72 74 0a 58 72 65 |P.X.I Cx|Sqrt.Xre|
|00001060| 70 6c 61 63 65 73 0a 58 | 2e 49 20 63 0a 58 62 79 |places.X|.I c.Xby|
|00001070| 20 74 68 65 20 60 60 70 | 72 69 6e 63 69 70 61 6c | the ``p|rincipal|
|00001080| 20 76 61 6c 75 65 27 27 | 20 6f 66 20 69 74 73 20 | value''| of its |
|00001090| 73 71 75 61 72 65 20 72 | 6f 6f 74 0a 58 28 6f 6e |square r|oot.X(on|
|000010a0| 65 20 68 61 76 69 6e 67 | 20 61 20 6e 6f 6e 2d 6e |e having| a non-n|
|000010b0| 65 67 61 74 69 76 65 20 | 69 6d 61 67 69 6e 61 72 |egative |imaginar|
|000010c0| 79 20 70 61 72 74 29 0a | 58 61 6e 64 20 72 65 74 |y part).|Xand ret|
|000010d0| 75 72 6e 73 20 61 20 70 | 6f 69 6e 74 65 72 20 74 |urns a p|ointer t|
|000010e0| 6f 20 74 68 65 20 72 65 | 73 75 6c 74 2e 0a 58 2e |o the re|sult..X.|
|000010f0| 50 0a 58 2e 49 20 43 78 | 53 63 61 6c 5c 5e 0a 58 |P.X.I Cx|Scal\^.X|
|00001100| 6d 75 6c 74 69 70 6c 69 | 65 73 0a 58 2e 49 20 63 |multipli|es.X.I c|
|00001110| 0a 58 62 79 20 74 68 65 | 20 73 63 61 6c 61 72 0a |.Xby the| scalar.|
|00001120| 58 2e 49 20 73 0a 58 61 | 6e 64 20 72 65 74 75 72 |X.I s.Xa|nd retur|
|00001130| 6e 73 20 61 20 70 6f 69 | 6e 74 65 72 20 74 6f 20 |ns a poi|nter to |
|00001140| 74 68 65 20 72 65 73 75 | 6c 74 2e 0a 58 2e 50 0a |the resu|lt..X.P.|
|00001150| 58 2e 49 20 43 78 4e 65 | 67 0a 58 6e 65 67 61 74 |X.I CxNe|g.Xnegat|
|00001160| 65 73 0a 58 2e 49 20 63 | 0a 58 61 6e 64 20 72 65 |es.X.I c|.Xand re|
|00001170| 74 75 72 6e 73 20 61 20 | 70 6f 69 6e 74 65 72 20 |turns a |pointer |
|00001180| 74 6f 20 74 68 65 20 72 | 65 73 75 6c 74 2e 0a 58 |to the r|esult..X|
|00001190| 2e 50 0a 58 2e 49 20 43 | 78 43 6f 6e 6a 0a 58 63 |.P.X.I C|xConj.Xc|
|000011a0| 6f 6e 6a 75 67 61 74 65 | 73 0a 58 2e 49 20 63 0a |onjugate|s.X.I c.|
|000011b0| 58 61 6e 64 20 72 65 74 | 75 72 6e 73 20 61 20 70 |Xand ret|urns a p|
|000011c0| 6f 69 6e 74 65 72 20 74 | 6f 20 74 68 65 20 72 65 |ointer t|o the re|
|000011d0| 73 75 6c 74 2e 0a 58 2e | 50 0a 58 2e 49 20 43 78 |sult..X.|P.X.I Cx|
|000011e0| 43 6f 70 79 0a 58 61 73 | 73 69 67 6e 73 20 74 68 |Copy.Xas|signs th|
|000011f0| 65 20 76 61 6c 75 65 20 | 6f 66 0a 58 2e 49 20 62 |e value |of.X.I b|
|00001200| 0a 58 74 6f 0a 58 2e 49 | 20 61 0a 58 61 6e 64 20 |.Xto.X.I| a.Xand |
|00001210| 72 65 74 75 72 6e 73 20 | 61 20 70 6f 69 6e 74 65 |returns |a pointe|
|00001220| 72 20 74 6f 20 74 68 65 | 20 72 65 73 75 6c 74 2e |r to the| result.|
|00001230| 0a 58 2e 50 0a 58 2e 49 | 20 43 78 43 6f 6e 73 0a |.X.P.X.I| CxCons.|
|00001240| 58 63 6f 6e 73 74 72 75 | 63 74 73 20 74 68 65 20 |Xconstru|cts the |
|00001250| 63 6f 6d 70 6c 65 78 20 | 6e 75 6d 62 65 72 0a 58 |complex |number.X|
|00001260| 2e 49 20 63 0a 58 66 72 | 6f 6d 20 69 74 73 20 72 |.I c.Xfr|om its r|
|00001270| 65 61 6c 20 61 6e 64 20 | 69 6d 61 67 69 6e 61 72 |eal and |imaginar|
|00001280| 79 20 70 61 72 74 73 0a | 58 2e 49 20 72 0a 58 61 |y parts.|X.I r.Xa|
|00001290| 6e 64 0a 58 2e 49 52 20 | 69 20 2c 0a 58 72 65 73 |nd.X.IR |i ,.Xres|
|000012a0| 70 65 63 74 69 76 65 6c | 79 2c 0a 58 61 6e 64 20 |pectivel|y,.Xand |
|000012b0| 72 65 74 75 72 6e 73 20 | 61 20 70 6f 69 6e 74 65 |returns |a pointe|
|000012c0| 72 20 74 6f 20 74 68 65 | 20 72 65 73 75 6c 74 2e |r to the| result.|
|000012d0| 0a 58 2e 50 0a 58 2e 49 | 20 43 78 50 68 73 72 0a |.X.P.X.I| CxPhsr.|
|000012e0| 58 63 6f 6e 73 74 72 75 | 63 74 73 20 74 68 65 20 |Xconstru|cts the |
|000012f0| 63 6f 6d 70 6c 65 78 20 | 6e 75 6d 62 65 72 0a 58 |complex |number.X|
|00001300| 2e 49 20 63 0a 58 66 72 | 6f 6d 20 69 74 73 20 60 |.I c.Xfr|om its `|
|00001310| 60 70 68 61 73 6f 72 27 | 27 20 61 6d 70 6c 69 74 |`phasor'|' amplit|
|00001320| 75 64 65 20 61 6e 64 20 | 70 68 61 73 65 20 28 67 |ude and |phase (g|
|00001330| 69 76 65 6e 20 69 6e 20 | 72 61 64 69 61 6e 73 29 |iven in |radians)|
|00001340| 0a 58 2e 49 20 6d 0a 58 | 61 6e 64 0a 58 2e 49 52 |.X.I m.X|and.X.IR|
|00001350| 20 70 20 2c 0a 58 72 65 | 73 70 65 63 74 69 76 65 | p ,.Xre|spective|
|00001360| 6c 79 2c 0a 58 61 6e 64 | 20 72 65 74 75 72 6e 73 |ly,.Xand| returns|
|00001370| 20 61 20 70 6f 69 6e 74 | 65 72 20 74 6f 20 74 68 | a point|er to th|
|00001380| 65 20 72 65 73 75 6c 74 | 2e 0a 58 2e 50 0a 58 2e |e result|..X.P.X.|
|00001390| 49 20 43 78 52 65 61 6c | 5c 5e 0a 58 72 65 74 75 |I CxReal|\^.Xretu|
|000013a0| 72 6e 73 20 74 68 65 20 | 72 65 61 6c 20 70 61 72 |rns the |real par|
|000013b0| 74 20 6f 66 20 74 68 65 | 20 63 6f 6d 70 6c 65 78 |t of the| complex|
|000013c0| 20 6e 75 6d 62 65 72 0a | 58 2e 49 52 20 63 20 2e | number.|X.IR c .|
|000013d0| 0a 58 2e 50 0a 58 2e 49 | 20 43 78 49 6d 61 67 0a |.X.P.X.I| CxImag.|
|000013e0| 58 72 65 74 75 72 6e 73 | 20 74 68 65 20 69 6d 61 |Xreturns| the ima|
|000013f0| 67 69 6e 61 72 79 20 70 | 61 72 74 20 6f 66 20 74 |ginary p|art of t|
|00001400| 68 65 20 63 6f 6d 70 6c | 65 78 20 6e 75 6d 62 65 |he compl|ex numbe|
|00001410| 72 0a 58 2e 49 52 20 63 | 20 2e 0a 58 2e 50 0a 58 |r.X.IR c| ..X.P.X|
|00001420| 2e 49 20 43 78 41 6d 70 | 6c 5c 5e 0a 58 72 65 74 |.I CxAmp|l\^.Xret|
|00001430| 75 72 6e 73 20 74 68 65 | 20 61 6d 70 6c 69 74 75 |urns the| amplitu|
|00001440| 64 65 20 6f 66 20 74 68 | 65 20 63 6f 6d 70 6c 65 |de of th|e comple|
|00001450| 78 20 6e 75 6d 62 65 72 | 0a 58 2e 49 52 20 63 20 |x number|.X.IR c |
|00001460| 2e 0a 58 2e 50 0a 58 2e | 49 20 43 78 50 68 61 73 |..X.P.X.|I CxPhas|
|00001470| 0a 58 72 65 74 75 72 6e | 73 20 74 68 65 20 70 68 |.Xreturn|s the ph|
|00001480| 61 73 65 20 6f 66 20 74 | 68 65 20 63 6f 6d 70 6c |ase of t|he compl|
|00001490| 65 78 20 6e 75 6d 62 65 | 72 0a 58 2e 49 52 20 63 |ex numbe|r.X.IR c|
|000014a0| 20 2c 0a 58 61 73 20 72 | 61 64 69 61 6e 73 20 69 | ,.Xas r|adians i|
|000014b0| 6e 20 74 68 65 20 72 61 | 6e 67 65 20 40 28 2d 20 |n the ra|nge @(- |
|000014c0| 70 69 20 2c 20 70 69 20 | 5d 40 2e 0a 58 2e 50 0a |pi , pi |]@..X.P.|
|000014d0| 58 2e 49 20 43 78 41 6c | 6c 6f 0a 58 61 6c 6c 6f |X.I CxAl|lo.Xallo|
|000014e0| 63 61 74 65 73 20 73 74 | 6f 72 61 67 65 20 66 6f |cates st|orage fo|
|000014f0| 72 20 61 0a 58 2e 42 20 | 63 6f 6d 70 6c 65 78 0a |r a.X.B |complex.|
|00001500| 58 64 61 74 75 6d 3b 20 | 69 74 20 72 65 74 75 72 |Xdatum; |it retur|
|00001510| 6e 73 0a 58 2e 53 4d 0a | 58 2e 42 20 4e 55 4c 4c |ns.X.SM.|X.B NULL|
|00001520| 0a 58 28 64 65 66 69 6e | 65 64 20 61 73 20 30 20 |.X(defin|ed as 0 |
|00001530| 69 6e 0a 58 2e 42 52 20 | 3c 73 74 64 69 6f 2e 68 |in.X.BR |<stdio.h|
|00001540| 3e 20 29 0a 58 69 66 20 | 6e 6f 74 20 65 6e 6f 75 |> ).Xif |not enou|
|00001550| 67 68 20 73 74 6f 72 61 | 67 65 20 69 73 20 61 76 |gh stora|ge is av|
|00001560| 61 69 6c 61 62 6c 65 2e | 0a 58 2e 50 0a 58 2e 49 |ailable.|.X.P.X.I|
|00001570| 20 43 78 46 72 65 65 0a | 58 72 65 6c 65 61 73 65 | CxFree.|Xrelease|
|00001580| 73 20 73 74 6f 72 61 67 | 65 20 70 72 65 76 69 6f |s storag|e previo|
|00001590| 75 73 6c 79 20 61 6c 6c | 6f 63 61 74 65 64 20 62 |usly all|ocated b|
|000015a0| 79 0a 58 2e 49 52 20 43 | 78 41 6c 6c 6f 20 2e 0a |y.X.IR C|xAllo ..|
|000015b0| 58 54 68 65 20 63 6f 6e | 74 65 6e 74 73 20 6f 66 |XThe con|tents of|
|000015c0| 20 73 75 63 68 20 73 74 | 6f 72 61 67 65 20 6d 75 | such st|orage mu|
|000015d0| 73 74 20 6e 6f 74 20 62 | 65 20 75 73 65 64 20 61 |st not b|e used a|
|000015e0| 66 74 65 72 77 61 72 64 | 2e 0a 58 2e 53 48 20 48 |fterward|..X.SH H|
|000015f0| 49 4e 54 53 0a 58 54 68 | 65 0a 58 2e 42 20 63 6f |INTS.XTh|e.X.B co|
|00001600| 6d 70 6c 65 78 0a 58 64 | 61 74 61 20 74 79 70 65 |mplex.Xd|ata type|
|00001610| 20 63 6f 6e 73 69 73 74 | 73 20 6f 66 20 72 65 61 | consist|s of rea|
|00001620| 6c 20 61 6e 64 20 69 6d | 61 67 69 6e 61 72 79 20 |l and im|aginary |
|00001630| 63 6f 6d 70 6f 6e 65 6e | 74 73 3b 0a 58 2e 49 20 |componen|ts;.X.I |
|00001640| 43 78 52 65 61 6c 5c 5e | 0a 58 61 6e 64 0a 58 2e |CxReal\^|.Xand.X.|
|00001650| 49 20 43 78 49 6d 61 67 | 0a 58 61 72 65 20 61 63 |I CxImag|.Xare ac|
|00001660| 74 75 61 6c 6c 79 20 6d | 61 63 72 6f 73 20 74 68 |tually m|acros th|
|00001670| 61 74 20 61 63 63 65 73 | 73 20 74 68 65 73 65 20 |at acces|s these |
|00001680| 63 6f 6d 70 6f 6e 65 6e | 74 73 20 64 69 72 65 63 |componen|ts direc|
|00001690| 74 6c 79 2e 0a 58 54 68 | 69 73 20 61 6c 6c 6f 77 |tly..XTh|is allow|
|000016a0| 73 20 61 64 64 72 65 73 | 73 65 73 20 6f 66 20 74 |s addres|ses of t|
|000016b0| 68 65 20 63 6f 6d 70 6f | 6e 65 6e 74 73 20 74 6f |he compo|nents to|
|000016c0| 20 62 65 20 74 61 6b 65 | 6e 2c 0a 58 61 73 20 69 | be take|n,.Xas i|
|000016d0| 6e 20 74 68 65 20 66 6f | 6c 6c 6f 77 69 6e 67 20 |n the fo|llowing |
|000016e0| 5c 73 2d 31 45 58 41 4d | 50 4c 45 5c 73 30 2e 0a |\s-1EXAM|PLE\s0..|
|000016f0| 58 2e 50 0a 58 54 68 65 | 20 63 6f 6d 70 6c 65 78 |X.P.XThe| complex|
|00001700| 20 66 75 6e 63 74 69 6f | 6e 73 20 61 72 65 20 64 | functio|ns are d|
|00001710| 65 73 69 67 6e 65 64 20 | 74 6f 20 62 65 20 6e 65 |esigned |to be ne|
|00001720| 73 74 65 64 3b 0a 58 73 | 65 65 20 74 68 65 20 66 |sted;.Xs|ee the f|
|00001730| 6f 6c 6c 6f 77 69 6e 67 | 20 5c 73 2d 31 45 58 41 |ollowing| \s-1EXA|
|00001740| 4d 50 4c 45 5c 73 30 2e | 0a 58 46 6f 72 20 74 68 |MPLE\s0.|.XFor th|
|00001750| 69 73 20 72 65 61 73 6f | 6e 2c 0a 58 6d 61 6e 79 |is reaso|n,.Xmany|
|00001760| 20 6f 66 20 74 68 65 6d | 20 6d 6f 64 69 66 79 20 | of them| modify |
|00001770| 74 68 65 20 63 6f 6e 74 | 65 6e 74 73 20 6f 66 20 |the cont|ents of |
|00001780| 74 68 65 69 72 20 66 69 | 72 73 74 20 70 61 72 61 |their fi|rst para|
|00001790| 6d 65 74 65 72 2e 0a 58 | 2e 49 20 43 78 43 6f 70 |meter..X|.I CxCop|
|000017a0| 79 0a 58 63 61 6e 20 62 | 65 20 75 73 65 64 20 74 |y.Xcan b|e used t|
|000017b0| 6f 20 63 72 65 61 74 65 | 20 61 20 60 60 77 6f 72 |o create| a ``wor|
|000017c0| 6b 69 6e 67 20 63 6f 70 | 79 27 27 20 6f 66 0a 58 |king cop|y'' of.X|
|000017d0| 2e 42 20 63 6f 6d 70 6c | 65 78 0a 58 64 61 74 61 |.B compl|ex.Xdata|
|000017e0| 20 74 68 61 74 20 77 6f | 75 6c 64 20 6f 74 68 65 | that wo|uld othe|
|000017f0| 72 77 69 73 65 20 62 65 | 20 6d 6f 64 69 66 69 65 |rwise be| modifie|
|00001800| 64 2e 0a 58 2e 50 0a 58 | 54 68 65 20 73 71 75 61 |d..X.P.X|The squa|
|00001810| 72 65 2d 72 6f 6f 74 20 | 66 75 6e 63 74 69 6f 6e |re-root |function|
|00001820| 20 69 73 20 69 6e 68 65 | 72 65 6e 74 6c 79 20 64 | is inhe|rently d|
|00001830| 6f 75 62 6c 65 2d 76 61 | 6c 75 65 64 3b 0a 58 69 |ouble-va|lued;.Xi|
|00001840| 6e 20 6d 6f 73 74 20 61 | 70 70 6c 69 63 61 74 69 |n most a|pplicati|
|00001850| 6f 6e 73 2c 20 62 6f 74 | 68 20 72 6f 6f 74 73 20 |ons, bot|h roots |
|00001860| 73 68 6f 75 6c 64 20 72 | 65 63 65 69 76 65 20 65 |should r|eceive e|
|00001870| 71 75 61 6c 20 63 6f 6e | 73 69 64 65 72 61 74 69 |qual con|siderati|
|00001880| 6f 6e 2e 0a 58 54 68 65 | 20 73 65 63 6f 6e 64 20 |on..XThe| second |
|00001890| 72 6f 6f 74 20 69 73 20 | 74 68 65 20 6e 65 67 61 |root is |the nega|
|000018a0| 74 69 76 65 20 6f 66 20 | 74 68 65 20 60 60 70 72 |tive of |the ``pr|
|000018b0| 69 6e 63 69 70 61 6c 20 | 76 61 6c 75 65 27 27 2e |incipal |value''.|
|000018c0| 0a 58 2e 62 70 0a 58 2e | 53 48 20 45 58 41 4d 50 |.X.bp.X.|SH EXAMP|
|000018d0| 4c 45 0a 58 54 68 65 20 | 66 6f 6c 6c 6f 77 69 6e |LE.XThe |followin|
|000018e0| 67 20 70 72 6f 67 72 61 | 6d 20 69 73 20 63 6f 6d |g progra|m is com|
|000018f0| 70 69 6c 65 64 20 62 79 | 20 74 68 65 20 63 6f 6d |piled by| the com|
|00001900| 6d 61 6e 64 0a 58 2e 62 | 72 0a 58 09 24 20 5c 66 |mand.X.b|r.X.$ \f|
|00001910| 49 63 63 20 5c 7c 5c 2d | 49 2f 75 73 72 2f 6c 6f |Icc \|\-|I/usr/lo|
|00001920| 63 61 6c 2f 69 6e 63 6c | 75 64 65 20 5c 7c 65 78 |cal/incl|ude \|ex|
|00001930| 61 6d 70 6c 65 2e 63 20 | 5c 7c 2f 75 73 72 2f 6c |ample.c |\|/usr/l|
|00001940| 6f 63 61 6c 2f 6c 69 62 | 2f 6c 69 62 63 6f 6d 70 |ocal/lib|/libcomp|
|00001950| 6c 65 78 2e 61 20 5c 7c | 5c 2d 6c 6d 5c 66 50 0a |lex.a \||\-lm\fP.|
|00001960| 58 2e 62 72 0a 58 49 74 | 20 72 65 61 64 73 20 69 |X.br.XIt| reads i|
|00001970| 6e 20 74 77 6f 20 63 6f | 6d 70 6c 65 78 20 76 65 |n two co|mplex ve|
|00001980| 63 74 6f 72 73 2c 0a 58 | 74 68 65 6e 20 63 6f 6d |ctors,.X|then com|
|00001990| 70 75 74 65 73 20 61 6e | 64 20 70 72 69 6e 74 73 |putes an|d prints|
|000019a0| 20 74 68 65 69 72 20 69 | 6e 6e 65 72 20 70 72 6f | their i|nner pro|
|000019b0| 64 75 63 74 2e 0a 58 2e | 73 70 0a 58 2e 50 0a 58 |duct..X.|sp.X.P.X|
|000019c0| 09 23 69 6e 63 6c 75 64 | 65 09 3c 73 74 64 69 6f |.#includ|e.<stdio|
|000019d0| 2e 68 3e 0a 58 2e 62 72 | 0a 58 09 23 69 6e 63 6c |.h>.X.br|.X.#incl|
|000019e0| 75 64 65 09 3c 63 6f 6d | 70 6c 65 78 2e 68 3e 0a |ude.<com|plex.h>.|
|000019f0| 58 2e 73 70 0a 58 09 6d | 61 69 6e 28 20 61 72 67 |X.sp.X.m|ain( arg|
|00001a00| 63 2c 20 61 72 67 76 20 | 29 0a 58 2e 62 72 0a 58 |c, argv |).X.br.X|
|00001a10| 09 09 69 6e 74 09 09 61 | 72 67 63 3b 0a 58 2e 62 |..int..a|rgc;.X.b|
|00001a20| 72 0a 58 09 09 63 68 61 | 72 09 09 2a 61 72 67 76 |r.X..cha|r..*argv|
|00001a30| 5b 5c 7c 5d 3b 0a 58 2e | 62 72 0a 58 09 09 7b 0a |[\|];.X.|br.X..{.|
|00001a40| 58 2e 62 72 0a 58 09 09 | 69 6e 74 09 09 6e 3b 09 |X.br.X..|int..n;.|
|00001a50| 09 2f 2a 20 23 20 65 6c | 65 6d 65 6e 74 73 20 69 |./* # el|ements i|
|00001a60| 6e 20 65 61 63 68 20 61 | 72 72 61 79 20 2a 2f 0a |n each a|rray */.|
|00001a70| 58 2e 62 72 0a 58 09 09 | 69 6e 74 09 09 69 3b 09 |X.br.X..|int..i;.|
|00001a80| 09 2f 2a 20 69 6e 64 65 | 78 65 73 20 61 72 72 61 |./* inde|xes arra|
|00001a90| 79 73 20 2a 2f 0a 58 2e | 62 72 0a 58 09 09 63 6f |ys */.X.|br.X..co|
|00001aa0| 6d 70 6c 65 78 09 09 61 | 5b 31 30 5d 2c 20 62 5b |mplex..a|[10], b[|
|00001ab0| 31 30 5d 3b 09 2f 2a 20 | 69 6e 70 75 74 20 76 65 |10];./* |input ve|
|00001ac0| 63 74 6f 72 73 20 2a 2f | 0a 58 2e 62 72 0a 58 09 |ctors */|.X.br.X.|
|00001ad0| 09 63 6f 6d 70 6c 65 78 | 09 09 73 3b 09 09 2f 2a |.complex|..s;../*|
|00001ae0| 20 61 63 63 75 6d 75 6c | 61 74 65 73 20 73 63 61 | accumul|ates sca|
|00001af0| 6c 61 72 20 70 72 6f 64 | 75 63 74 20 2a 2f 0a 58 |lar prod|uct */.X|
|00001b00| 2e 62 72 0a 58 09 09 63 | 6f 6d 70 6c 65 78 09 09 |.br.X..c|omplex..|
|00001b10| 2a 63 20 3d 20 43 78 41 | 6c 6c 6f 28 5c 7c 29 3b |*c = CxA|llo(\|);|
|00001b20| 09 2f 2a 20 68 6f 6c 64 | 73 20 63 72 6f 73 73 2d |./* hold|s cross-|
|00001b30| 74 65 72 6d 20 2a 2f 0a | 58 2e 73 70 0a 58 09 09 |term */.|X.sp.X..|
|00001b40| 69 66 20 28 20 63 20 3d | 3d 20 4e 55 4c 4c 20 29 |if ( c =|= NULL )|
|00001b50| 0a 58 2e 62 72 0a 58 09 | 09 09 7b 0a 58 2e 62 72 |.X.br.X.|..{.X.br|
|00001b60| 0a 58 09 09 09 28 76 6f | 69 64 29 66 70 72 69 6e |.X...(vo|id)fprin|
|00001b70| 74 66 28 20 73 74 64 65 | 72 72 2c 20 60 60 6e 6f |tf( stde|rr, ``no|
|00001b80| 74 20 65 6e 6f 75 67 68 | 20 6d 65 6d 6f 72 79 5c |t enough| memory\|
|00001b90| 65 6e 27 27 20 29 3b 0a | 58 2e 62 72 0a 58 09 09 |en'' );.|X.br.X..|
|00001ba0| 09 72 65 74 75 72 6e 20 | 31 3b 0a 58 2e 62 72 0a |.return |1;.X.br.|
|00001bb0| 58 09 09 09 7d 0a 58 2e | 62 72 0a 58 09 09 28 76 |X...}.X.|br.X..(v|
|00001bc0| 6f 69 64 29 70 72 69 6e | 74 66 28 20 60 60 5c 65 |oid)prin|tf( ``\e|
|00001bd0| 6e 65 6e 74 65 72 20 6e | 75 6d 62 65 72 20 6f 66 |nenter n|umber of|
|00001be0| 20 65 6c 65 6d 65 6e 74 | 73 3a 20 27 27 20 29 3b | element|s: '' );|
|00001bf0| 0a 58 2e 62 72 0a 58 09 | 09 28 76 6f 69 64 29 73 |.X.br.X.|.(void)s|
|00001c00| 63 61 6e 66 28 20 60 60 | 20 25 64 27 27 2c 20 26 |canf( ``| %d'', &|
|00001c10| 6e 20 29 3b 0a 58 2e 62 | 72 0a 58 09 09 2f 2a 20 |n );.X.b|r.X../* |
|00001c20| 28 54 68 65 72 65 20 72 | 65 61 6c 6c 79 20 73 68 |(There r|eally sh|
|00001c30| 6f 75 6c 64 20 62 65 20 | 73 6f 6d 65 20 69 6e 70 |ould be |some inp|
|00001c40| 75 74 20 76 61 6c 69 64 | 61 74 69 6f 6e 20 68 65 |ut valid|ation he|
|00001c50| 72 65 2e 29 20 2a 2f 0a | 58 2e 62 72 0a 58 09 09 |re.) */.|X.br.X..|
|00001c60| 28 76 6f 69 64 29 20 70 | 72 69 6e 74 66 28 20 60 |(void) p|rintf( `|
|00001c70| 60 5c 65 6e 65 6e 74 65 | 72 20 72 65 61 6c 2c 20 |`\enente|r real, |
|00001c80| 69 6d 61 67 69 6e 61 72 | 79 20 70 61 69 72 73 20 |imaginar|y pairs |
|00001c90| 66 6f 72 20 66 69 72 73 | 74 20 61 72 72 61 79 3a |for firs|t array:|
|00001ca0| 5c 65 6e 27 27 20 29 3b | 0a 58 2e 62 72 0a 58 09 |\en'' );|.X.br.X.|
|00001cb0| 09 66 6f 72 20 28 20 69 | 20 3d 20 30 3b 20 69 20 |.for ( i| = 0; i |
|00001cc0| 3c 20 6e 3b 20 2b 2b 69 | 20 29 0a 58 2e 62 72 0a |< n; ++i| ).X.br.|
|00001cd0| 58 09 09 09 28 76 6f 69 | 64 29 73 63 61 6e 66 28 |X...(voi|d)scanf(|
|00001ce0| 20 60 60 20 25 6c 67 20 | 25 6c 67 27 27 2c 20 26 | `` %lg |%lg'', &|
|00001cf0| 43 78 52 65 61 6c 28 20 | 26 61 5b 69 5d 20 29 2c |CxReal( |&a[i] ),|
|00001d00| 20 26 43 78 49 6d 61 67 | 28 20 26 61 5b 69 5d 20 | &CxImag|( &a[i] |
|00001d10| 29 20 29 3b 0a 58 2e 62 | 72 0a 58 09 09 28 76 6f |) );.X.b|r.X..(vo|
|00001d20| 69 64 29 70 72 69 6e 74 | 66 28 20 60 60 5c 65 6e |id)print|f( ``\en|
|00001d30| 65 6e 74 65 72 20 72 65 | 61 6c 2c 20 69 6d 61 67 |enter re|al, imag|
|00001d40| 69 6e 61 72 79 20 70 61 | 69 72 73 20 66 6f 72 20 |inary pa|irs for |
|00001d50| 73 65 63 6f 6e 64 20 61 | 72 72 61 79 3a 5c 65 6e |second a|rray:\en|
|00001d60| 27 27 20 29 3b 0a 58 2e | 62 72 0a 58 09 09 66 6f |'' );.X.|br.X..fo|
|00001d70| 72 20 28 20 69 20 3d 20 | 30 3b 20 69 20 3c 20 6e |r ( i = |0; i < n|
|00001d80| 3b 20 2b 2b 69 20 29 0a | 58 2e 62 72 0a 58 09 09 |; ++i ).|X.br.X..|
|00001d90| 09 28 76 6f 69 64 29 73 | 63 61 6e 66 28 20 60 60 |.(void)s|canf( ``|
|00001da0| 20 25 6c 67 20 25 6c 67 | 27 27 2c 20 26 43 78 52 | %lg %lg|'', &CxR|
|00001db0| 65 61 6c 28 20 26 62 5b | 69 5d 20 29 2c 20 26 43 |eal( &b[|i] ), &C|
|00001dc0| 78 49 6d 61 67 28 20 26 | 62 5b 69 5d 20 29 20 29 |xImag( &|b[i] ) )|
|00001dd0| 3b 0a 58 2e 62 72 0a 58 | 09 09 28 76 6f 69 64 29 |;.X.br.X|..(void)|
|00001de0| 43 78 43 6f 6e 73 28 20 | 26 73 2c 20 30 2e 30 2c |CxCons( |&s, 0.0,|
|00001df0| 20 30 2e 30 20 29 3b 09 | 2f 2a 20 69 6e 69 74 69 | 0.0 );.|/* initi|
|00001e00| 61 6c 69 7a 65 20 61 63 | 63 75 6d 75 6c 61 74 6f |alize ac|cumulato|
|00001e10| 72 20 2a 2f 0a 58 2e 62 | 72 0a 58 09 09 66 6f 72 |r */.X.b|r.X..for|
|00001e20| 20 28 20 69 20 3d 20 30 | 3b 20 69 20 3c 20 6e 3b | ( i = 0|; i < n;|
|00001e30| 20 2b 2b 69 20 29 0a 58 | 2e 62 72 0a 58 09 09 09 | ++i ).X|.br.X...|
|00001e40| 28 76 6f 69 64 29 43 78 | 41 64 64 28 20 26 73 2c |(void)Cx|Add( &s,|
|00001e50| 20 43 78 4d 75 6c 28 20 | 26 61 5b 69 5d 2c 20 43 | CxMul( |&a[i], C|
|00001e60| 78 43 6f 6e 6a 28 20 43 | 78 43 6f 70 79 28 20 63 |xConj( C|xCopy( c|
|00001e70| 2c 20 26 62 5b 69 5d 20 | 29 20 29 20 29 20 29 3b |, &b[i] |) ) ) );|
|00001e80| 0a 58 2e 62 72 0a 58 09 | 09 28 76 6f 69 64 29 70 |.X.br.X.|.(void)p|
|00001e90| 72 69 6e 74 66 28 20 60 | 60 5c 65 6e 70 72 6f 64 |rintf( `|`\enprod|
|00001ea0| 75 63 74 20 69 73 20 28 | 25 67 2c 25 67 29 5c 65 |uct is (|%g,%g)\e|
|00001eb0| 6e 27 27 2c 20 43 78 52 | 65 61 6c 28 20 26 73 20 |n'', CxR|eal( &s |
|00001ec0| 29 2c 20 43 78 49 6d 61 | 67 28 20 26 73 20 29 20 |), CxIma|g( &s ) |
|00001ed0| 29 3b 0a 58 2e 62 72 0a | 58 09 09 43 78 46 72 65 |);.X.br.|X..CxFre|
|00001ee0| 65 28 20 63 20 29 3b 0a | 58 2e 62 72 0a 58 09 09 |e( c );.|X.br.X..|
|00001ef0| 72 65 74 75 72 6e 20 30 | 3b 0a 58 2e 62 72 0a 58 |return 0|;.X.br.X|
|00001f00| 09 09 7d 0a 58 2e 53 48 | 20 46 49 4c 45 53 0a 58 |..}.X.SH| FILES.X|
|00001f10| 2f 75 73 72 2f 6c 6f 63 | 61 6c 2f 69 6e 63 6c 75 |/usr/loc|al/inclu|
|00001f20| 64 65 2f 63 6f 6d 70 6c | 65 78 2e 68 09 09 68 65 |de/compl|ex.h..he|
|00001f30| 61 64 65 72 20 66 69 6c | 65 20 63 6f 6e 74 61 69 |ader fil|e contai|
|00001f40| 6e 69 6e 67 20 64 65 66 | 69 6e 69 74 69 6f 6e 73 |ning def|initions|
|00001f50| 0a 58 2e 62 72 0a 58 2f | 75 73 72 2f 6c 6f 63 61 |.X.br.X/|usr/loca|
|00001f60| 6c 2f 6c 69 62 2f 6c 69 | 62 63 6f 6d 70 6c 65 78 |l/lib/li|bcomplex|
|00001f70| 2e 61 09 09 63 6f 6d 70 | 6c 65 78 20 72 75 6e 2d |.a..comp|lex run-|
|00001f80| 74 69 6d 65 20 73 75 70 | 70 6f 72 74 20 6c 69 62 |time sup|port lib|
|00001f90| 72 61 72 79 0a 58 2e 53 | 48 20 41 55 54 48 4f 52 |rary.X.S|H AUTHOR|
|00001fa0| 53 0a 58 44 6f 75 67 6c | 61 73 20 41 2e 20 47 77 |S.XDougl|as A. Gw|
|00001fb0| 79 6e 2c 20 42 52 4c 2f | 56 4c 44 2d 56 4d 42 0a |yn, BRL/|VLD-VMB.|
|00001fc0| 58 2e 62 72 0a 58 4a 65 | 66 66 20 48 61 6e 65 73 |X.br.XJe|ff Hanes|
|00001fd0| 2c 20 42 52 4c 2f 56 4c | 44 2d 56 4d 42 20 28 6f |, BRL/VL|D-VMB (o|
|00001fe0| 72 69 67 69 6e 61 6c 20 | 76 65 72 73 69 6f 6e 20 |riginal |version |
|00001ff0| 6f 66 0a 58 2e 49 52 20 | 43 78 53 71 72 74 5c 5e |of.X.IR |CxSqrt\^|
|00002000| 20 29 0a 45 4e 44 5f 4f | 46 5f 63 6f 6d 70 6c 65 | ).END_O|F_comple|
|00002010| 78 2e 33 0a 69 66 20 74 | 65 73 74 20 35 38 33 30 |x.3.if t|est 5830|
|00002020| 20 2d 6e 65 20 60 77 63 | 20 2d 63 20 3c 63 6f 6d | -ne `wc| -c <com|
|00002030| 70 6c 65 78 2e 33 60 3b | 20 74 68 65 6e 0a 20 20 |plex.3`;| then.  |
|00002040| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 63 |  echo s|har: \"c|
|00002050| 6f 6d 70 6c 65 78 2e 33 | 5c 22 20 75 6e 70 61 63 |omplex.3|\" unpac|
|00002060| 6b 65 64 20 77 69 74 68 | 20 77 72 6f 6e 67 20 73 |ked with| wrong s|
|00002070| 69 7a 65 21 0a 66 69 0a | 23 20 65 6e 64 20 6f 66 |ize!.fi.|# end of|
|00002080| 20 6f 76 65 72 77 72 69 | 74 69 6e 67 20 63 68 65 | overwri|ting che|
|00002090| 63 6b 0a 66 69 0a 69 66 | 20 74 65 73 74 20 2d 66 |ck.fi.if| test -f|
|000020a0| 20 63 6f 6d 70 6c 65 78 | 2e 68 20 2d 61 20 22 24 | complex|.h -a "$|
|000020b0| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|000020c0| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen .  e|cho shar|
|000020d0| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 6f 76 65 72 2d |: Will n|ot over-|
|000020e0| 77 72 69 74 65 20 65 78 | 69 73 74 69 6e 67 20 66 |write ex|isting f|
|000020f0| 69 6c 65 20 5c 22 63 6f | 6d 70 6c 65 78 2e 68 5c |ile \"co|mplex.h\|
|00002100| 22 0a 65 6c 73 65 0a 65 | 63 68 6f 20 73 68 61 72 |".else.e|cho shar|
|00002110| 3a 20 45 78 74 72 61 63 | 74 69 6e 67 20 5c 22 63 |: Extrac|ting \"c|
|00002120| 6f 6d 70 6c 65 78 2e 68 | 5c 22 20 5c 28 39 36 36 |omplex.h|\" \(966|
|00002130| 20 63 68 61 72 61 63 74 | 65 72 73 5c 29 0a 73 65 | charact|ers\).se|
|00002140| 64 20 22 73 2f 5e 58 2f | 2f 22 20 3e 63 6f 6d 70 |d "s/^X/|/" >comp|
|00002150| 6c 65 78 2e 68 20 3c 3c | 27 45 4e 44 5f 4f 46 5f |lex.h <<|'END_OF_|
|00002160| 63 6f 6d 70 6c 65 78 2e | 68 27 0a 58 2f 2a 0a 58 |complex.|h'.X/*.X|
|00002170| 09 3c 63 6f 6d 70 6c 65 | 78 2e 68 3e 20 2d 2d 20 |.<comple|x.h> -- |
|00002180| 64 65 66 69 6e 69 74 69 | 6f 6e 73 20 66 6f 72 20 |definiti|ons for |
|00002190| 63 6f 6d 70 6c 65 78 20 | 61 72 69 74 68 6d 65 74 |complex |arithmet|
|000021a0| 69 63 20 72 6f 75 74 69 | 6e 65 73 0a 58 0a 58 09 |ic routi|nes.X.X.|
|000021b0| 6c 61 73 74 20 65 64 69 | 74 3a 09 38 36 2f 30 31 |last edi|t:.86/01|
|000021c0| 2f 30 34 09 44 20 41 20 | 47 77 79 6e 0a 58 0a 58 |/04.D A |Gwyn.X.X|
|000021d0| 09 53 43 43 53 20 49 44 | 3a 09 40 28 23 29 63 6f |.SCCS ID|:.@(#)co|
|000021e0| 6d 70 6c 65 78 2e 68 09 | 31 2e 31 20 28 6d 6f 64 |mplex.h.|1.1 (mod|
|000021f0| 69 66 69 65 64 20 66 6f | 72 20 70 75 62 6c 69 63 |ified fo|r public|
|00002200| 20 76 65 72 73 69 6f 6e | 29 0a 58 2a 2f 0a 58 0a | version|).X*/.X.|
|00002210| 58 2f 2a 20 22 63 6f 6d | 70 6c 65 78 20 6e 75 6d |X/* "com|plex num|
|00002220| 62 65 72 22 20 64 61 74 | 61 20 74 79 70 65 3a 20 |ber" dat|a type: |
|00002230| 2a 2f 0a 58 0a 58 74 79 | 70 65 64 65 66 20 73 74 |*/.X.Xty|pedef st|
|00002240| 72 75 63 74 0a 58 09 7b | 0a 58 09 64 6f 75 62 6c |ruct.X.{|.X.doubl|
|00002250| 65 09 09 72 65 3b 09 09 | 2f 2a 20 72 65 61 6c 20 |e..re;..|/* real |
|00002260| 70 61 72 74 20 2a 2f 0a | 58 09 64 6f 75 62 6c 65 |part */.|X.double|
|00002270| 09 09 69 6d 3b 09 09 2f | 2a 20 69 6d 61 67 69 6e |..im;../|* imagin|
|00002280| 61 72 79 20 70 61 72 74 | 20 2a 2f 0a 58 09 7d 09 |ary part| */.X.}.|
|00002290| 63 6f 6d 70 6c 65 78 3b | 0a 58 0a 58 2f 2a 20 22 |complex;|.X.X/* "|
|000022a0| 54 68 65 20 66 75 74 75 | 72 65 20 69 73 20 6e 6f |The futu|re is no|
|000022b0| 77 22 3a 20 2a 2f 0a 58 | 0a 58 23 69 66 64 65 66 |w": */.X|.X#ifdef|
|000022c0| 20 5f 5f 53 54 44 43 5f | 5f 09 2f 2a 20 58 33 4a | __STDC_|_./* X3J|
|000022d0| 31 31 20 2a 2f 0a 58 23 | 64 65 66 69 6e 65 09 5f |11 */.X#|define._|
|000022e0| 43 78 47 65 6e 50 74 72 | 09 76 6f 69 64 20 2a 09 |CxGenPtr|.void *.|
|000022f0| 09 2f 2a 20 67 65 6e 65 | 72 69 63 20 70 6f 69 6e |./* gene|ric poin|
|00002300| 74 65 72 20 74 79 70 65 | 20 2a 2f 0a 58 23 65 6c |ter type| */.X#el|
|00002310| 73 65 09 09 2f 2a 20 4b | 26 52 20 2a 2f 0a 58 23 |se../* K|&R */.X#|
|00002320| 64 65 66 69 6e 65 09 5f | 43 78 47 65 6e 50 74 72 |define._|CxGenPtr|
|00002330| 09 63 68 61 72 20 2a 09 | 09 2f 2a 20 67 65 6e 65 |.char *.|./* gene|
|00002340| 72 69 63 20 70 6f 69 6e | 74 65 72 20 74 79 70 65 |ric poin|ter type|
|00002350| 20 2a 2f 0a 58 23 65 6e | 64 69 66 0a 58 0a 58 2f | */.X#en|dif.X.X/|
|00002360| 2a 20 66 75 6e 63 74 69 | 6f 6e 73 20 74 68 61 74 |* functi|ons that|
|00002370| 20 61 72 65 20 63 6f 72 | 72 65 63 74 6c 79 20 64 | are cor|rectly d|
|00002380| 6f 6e 65 20 61 73 20 6d | 61 63 72 6f 73 3a 20 2a |one as m|acros: *|
|00002390| 2f 0a 58 0a 58 23 64 65 | 66 69 6e 65 09 43 78 41 |/.X.X#de|fine.CxA|
|000023a0| 6c 6c 6f 28 29 09 09 28 | 28 63 6f 6d 70 6c 65 78 |llo()..(|(complex|
|000023b0| 20 2a 29 6d 61 6c 6c 6f | 63 28 20 73 69 7a 65 6f | *)mallo|c( sizeo|
|000023c0| 66 20 28 63 6f 6d 70 6c | 65 78 29 20 29 29 0a 58 |f (compl|ex) )).X|
|000023d0| 23 64 65 66 69 6e 65 09 | 43 78 46 72 65 65 28 20 |#define.|CxFree( |
|000023e0| 63 70 20 29 09 09 66 72 | 65 65 28 20 28 5f 43 78 |cp )..fr|ee( (_Cx|
|000023f0| 47 65 6e 50 74 72 29 28 | 63 70 29 20 29 0a 58 23 |GenPtr)(|cp) ).X#|
|00002400| 64 65 66 69 6e 65 09 43 | 78 4e 65 67 28 20 63 70 |define.C|xNeg( cp|
|00002410| 20 29 09 09 43 78 53 63 | 61 6c 28 20 63 70 2c 20 | )..CxSc|al( cp, |
|00002420| 2d 31 2e 30 20 29 0a 58 | 23 64 65 66 69 6e 65 09 |-1.0 ).X|#define.|
|00002430| 43 78 52 65 61 6c 28 20 | 63 70 20 29 09 09 28 63 |CxReal( |cp )..(c|
|00002440| 70 29 2d 3e 72 65 0a 58 | 23 64 65 66 69 6e 65 09 |p)->re.X|#define.|
|00002450| 43 78 49 6d 61 67 28 20 | 63 70 20 29 09 09 28 63 |CxImag( |cp )..(c|
|00002460| 70 29 2d 3e 69 6d 0a 58 | 0a 58 65 78 74 65 72 6e |p)->im.X|.Xextern|
|00002470| 20 76 6f 69 64 09 09 66 | 72 65 65 28 29 3b 0a 58 | void..f|ree();.X|
|00002480| 65 78 74 65 72 6e 20 5f | 43 78 47 65 6e 50 74 72 |extern _|CxGenPtr|
|00002490| 09 6d 61 6c 6c 6f 63 28 | 29 3b 0a 58 0a 58 2f 2a |.malloc(|);.X.X/*|
|000024a0| 20 6c 69 62 72 61 72 79 | 20 66 75 6e 63 74 69 6f | library| functio|
|000024b0| 6e 73 3a 20 2a 2f 0a 58 | 0a 58 65 78 74 65 72 6e |ns: */.X|.Xextern|
|000024c0| 20 64 6f 75 62 6c 65 09 | 43 78 41 6d 70 6c 28 29 | double.|CxAmpl()|
|000024d0| 2c 20 43 78 50 68 61 73 | 28 29 3b 0a 58 65 78 74 |, CxPhas|();.Xext|
|000024e0| 65 72 6e 20 63 6f 6d 70 | 6c 65 78 09 2a 43 78 41 |ern comp|lex.*CxA|
|000024f0| 64 64 28 29 2c 20 2a 43 | 78 43 6f 6e 6a 28 29 2c |dd(), *C|xConj(),|
|00002500| 20 2a 43 78 43 6f 6e 73 | 28 29 2c 20 2a 43 78 43 | *CxCons|(), *CxC|
|00002510| 6f 70 79 28 29 2c 20 2a | 43 78 44 69 76 28 29 2c |opy(), *|CxDiv(),|
|00002520| 0a 58 09 09 2a 43 78 4d | 75 6c 28 29 2c 20 2a 43 |.X..*CxM|ul(), *C|
|00002530| 78 50 68 73 72 28 29 2c | 20 2a 43 78 53 63 61 6c |xPhsr(),| *CxScal|
|00002540| 28 29 2c 20 2a 43 78 53 | 71 72 74 28 29 2c 20 2a |(), *CxS|qrt(), *|
|00002550| 43 78 53 75 62 28 29 3b | 0a 45 4e 44 5f 4f 46 5f |CxSub();|.END_OF_|
|00002560| 63 6f 6d 70 6c 65 78 2e | 68 0a 69 66 20 74 65 73 |complex.|h.if tes|
|00002570| 74 20 39 36 36 20 2d 6e | 65 20 60 77 63 20 2d 63 |t 966 -n|e `wc -c|
|00002580| 20 3c 63 6f 6d 70 6c 65 | 78 2e 68 60 3b 20 74 68 | <comple|x.h`; th|
|00002590| 65 6e 0a 20 20 20 20 65 | 63 68 6f 20 73 68 61 72 |en.    e|cho shar|
|000025a0| 3a 20 5c 22 63 6f 6d 70 | 6c 65 78 2e 68 5c 22 20 |: \"comp|lex.h\" |
|000025b0| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|000025c0| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|000025d0| 6e 64 20 6f 66 20 6f 76 | 65 72 77 72 69 74 69 6e |nd of ov|erwritin|
|000025e0| 67 20 63 68 65 63 6b 0a | 66 69 0a 69 66 20 74 65 |g check.|fi.if te|
|000025f0| 73 74 20 2d 66 20 63 78 | 5f 74 65 73 74 2e 63 20 |st -f cx|_test.c |
|00002600| 2d 61 20 22 24 7b 31 7d | 22 20 21 3d 20 22 2d 63 |-a "${1}|" != "-c|
|00002610| 22 20 3b 20 74 68 65 6e | 20 0a 20 20 65 63 68 6f |" ; then| .  echo|
|00002620| 20 73 68 61 72 3a 20 57 | 69 6c 6c 20 6e 6f 74 20 | shar: W|ill not |
|00002630| 6f 76 65 72 2d 77 72 69 | 74 65 20 65 78 69 73 74 |over-wri|te exist|
|00002640| 69 6e 67 20 66 69 6c 65 | 20 5c 22 63 78 5f 74 65 |ing file| \"cx_te|
|00002650| 73 74 2e 63 5c 22 0a 65 | 6c 73 65 0a 65 63 68 6f |st.c\".e|lse.echo|
|00002660| 20 73 68 61 72 3a 20 45 | 78 74 72 61 63 74 69 6e | shar: E|xtractin|
|00002670| 67 20 5c 22 63 78 5f 74 | 65 73 74 2e 63 5c 22 20 |g \"cx_t|est.c\" |
|00002680| 5c 28 34 32 35 30 20 63 | 68 61 72 61 63 74 65 72 |\(4250 c|haracter|
|00002690| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|000026a0| 20 3e 63 78 5f 74 65 73 | 74 2e 63 20 3c 3c 27 45 | >cx_tes|t.c <<'E|
|000026b0| 4e 44 5f 4f 46 5f 63 78 | 5f 74 65 73 74 2e 63 27 |ND_OF_cx|_test.c'|
|000026c0| 0a 58 2f 2a 0a 58 09 63 | 74 65 73 74 20 2d 2d 20 |.X/*.X.c|test -- |
|000026d0| 63 6f 6d 70 6c 65 78 20 | 61 72 69 74 68 6d 65 74 |complex |arithmet|
|000026e0| 69 63 20 74 65 73 74 0a | 58 0a 58 09 6c 61 73 74 |ic test.|X.X.last|
|000026f0| 20 65 64 69 74 3a 09 38 | 36 2f 30 31 2f 30 34 09 | edit:.8|6/01/04.|
|00002700| 44 20 41 20 47 77 79 6e | 0a 58 0a 58 09 53 43 43 |D A Gwyn|.X.X.SCC|
|00002710| 53 20 49 44 3a 09 40 28 | 23 29 63 78 5f 74 65 73 |S ID:.@(|#)cx_tes|
|00002720| 74 2e 63 09 31 2e 31 20 | 28 6d 6f 64 69 66 69 65 |t.c.1.1 |(modifie|
|00002730| 64 20 66 6f 72 20 70 75 | 62 6c 69 63 20 76 65 72 |d for pu|blic ver|
|00002740| 73 69 6f 6e 29 0a 58 2a | 2f 0a 58 0a 58 23 69 6e |sion).X*|/.X.X#in|
|00002750| 63 6c 75 64 65 09 3c 73 | 74 64 69 6f 2e 68 3e 0a |clude.<s|tdio.h>.|
|00002760| 58 23 69 6e 63 6c 75 64 | 65 09 3c 6d 61 74 68 2e |X#includ|e.<math.|
|00002770| 68 3e 0a 58 0a 58 23 69 | 6e 63 6c 75 64 65 09 3c |h>.X.X#i|nclude.<|
|00002780| 63 6f 6d 70 6c 65 78 2e | 68 3e 0a 58 0a 58 23 64 |complex.|h>.X.X#d|
|00002790| 65 66 69 6e 65 20 44 45 | 47 52 41 44 09 35 37 2e |efine DE|GRAD.57.|
|000027a0| 32 39 35 37 37 39 35 31 | 33 30 38 32 33 32 30 38 |29577951|30823208|
|000027b0| 37 36 37 39 38 31 35 34 | 38 31 34 31 30 35 31 37 |76798154|81410517|
|000027c0| 30 33 33 32 34 30 35 34 | 37 32 34 36 36 35 36 34 |03324054|72466564|
|000027d0| 32 0a 58 09 09 09 09 09 | 2f 2a 20 64 65 67 72 65 |2.X.....|/* degre|
|000027e0| 65 73 20 70 65 72 20 72 | 61 64 69 61 6e 20 2a 2f |es per r|adian */|
|000027f0| 0a 58 23 64 65 66 69 6e | 65 20 41 62 73 28 20 78 |.X#defin|e Abs( x|
|00002800| 20 29 09 28 28 78 29 20 | 3c 20 30 20 3f 20 2d 28 | ).((x) |< 0 ? -(|
|00002810| 78 29 20 3a 20 28 78 29 | 29 0a 58 23 64 65 66 69 |x) : (x)|).X#defi|
|00002820| 6e 65 20 4d 61 78 28 20 | 61 2c 20 62 20 29 09 28 |ne Max( |a, b ).(|
|00002830| 28 61 29 20 3e 20 28 62 | 29 20 3f 20 28 61 29 20 |(a) > (b|) ? (a) |
|00002840| 3a 20 28 62 29 29 0a 58 | 0a 58 65 78 74 65 72 6e |: (b)).X|.Xextern|
|00002850| 20 76 6f 69 64 09 65 78 | 69 74 28 29 3b 0a 58 0a | void.ex|it();.X.|
|00002860| 58 23 64 65 66 69 6e 65 | 09 50 72 69 6e 74 66 09 |X#define|.Printf.|
|00002870| 28 76 6f 69 64 29 70 72 | 69 6e 74 66 0a 58 0a 58 |(void)pr|intf.X.X|
|00002880| 23 64 65 66 69 6e 65 09 | 54 4f 4c 09 31 2e 30 65 |#define.|TOL.1.0e|
|00002890| 2d 31 30 09 09 09 2f 2a | 20 74 6f 6c 65 72 61 6e |-10.../*| toleran|
|000028a0| 63 65 20 66 6f 72 20 63 | 68 65 63 6b 73 20 2a 2f |ce for c|hecks */|
|000028b0| 0a 58 0a 58 73 74 61 74 | 69 63 20 69 6e 74 09 65 |.X.Xstat|ic int.e|
|000028c0| 72 72 73 20 3d 20 30 3b | 09 09 2f 2a 20 74 61 6c |rrs = 0;|../* tal|
|000028d0| 6c 79 20 65 72 72 6f 72 | 73 20 2a 2f 0a 58 0a 58 |ly error|s */.X.X|
|000028e0| 73 74 61 74 69 63 20 76 | 6f 69 64 09 43 43 68 65 |static v|oid.CChe|
|000028f0| 63 6b 28 29 2c 20 52 43 | 68 65 63 6b 28 29 3b 0a |ck(), RC|heck();.|
|00002900| 58 73 74 61 74 69 63 20 | 64 6f 75 62 6c 65 09 52 |Xstatic |double.R|
|00002910| 65 6c 44 69 66 28 29 3b | 0a 58 0a 58 0a 58 2f 2a |elDif();|.X.X.X/*|
|00002920| 41 52 47 53 55 53 45 44 | 2a 2f 0a 58 6d 61 69 6e |ARGSUSED|*/.Xmain|
|00002930| 28 20 61 72 67 63 2c 20 | 61 72 67 76 20 29 0a 58 |( argc, |argv ).X|
|00002940| 09 69 6e 74 09 61 72 67 | 63 3b 0a 58 09 63 68 61 |.int.arg|c;.X.cha|
|00002950| 72 09 2a 61 72 67 76 5b | 5d 3b 0a 58 09 7b 0a 58 |r.*argv[|];.X.{.X|
|00002960| 09 63 6f 6d 70 6c 65 78 | 20 61 2c 20 2a 62 70 2c |.complex| a, *bp,|
|00002970| 20 2a 63 70 3b 0a 58 0a | 58 09 2f 2a 20 43 78 41 | *cp;.X.|X./* CxA|
|00002980| 6c 6c 6f 20 74 65 73 74 | 20 2a 2f 0a 58 09 62 70 |llo test| */.X.bp|
|00002990| 20 3d 20 43 78 41 6c 6c | 6f 28 29 3b 0a 58 09 69 | = CxAll|o();.X.i|
|000029a0| 66 20 28 20 62 70 20 3d | 3d 20 4e 55 4c 4c 20 29 |f ( bp =|= NULL )|
|000029b0| 0a 58 09 09 7b 0a 58 09 | 09 50 72 69 6e 74 66 28 |.X..{.X.|.Printf(|
|000029c0| 20 22 43 78 41 6c 6c 6f | 20 66 61 69 6c 65 64 5c | "CxAllo| failed\|
|000029d0| 6e 22 20 29 3b 0a 58 09 | 09 65 78 69 74 28 20 31 |n" );.X.|.exit( 1|
|000029e0| 20 29 3b 0a 58 09 09 7d | 0a 58 0a 58 09 2f 2a 20 | );.X..}|.X.X./* |
|000029f0| 43 78 52 65 61 6c 2c 20 | 43 78 49 6d 61 67 20 74 |CxReal, |CxImag t|
|00002a00| 65 73 74 20 2a 2f 0a 58 | 09 43 78 52 65 61 6c 28 |est */.X|.CxReal(|
|00002a10| 20 62 70 20 29 20 3d 20 | 31 2e 30 3b 0a 58 09 43 | bp ) = |1.0;.X.C|
|00002a20| 78 49 6d 61 67 28 20 62 | 70 20 29 20 3d 20 32 2e |xImag( b|p ) = 2.|
|00002a30| 30 3b 0a 58 09 52 43 68 | 65 63 6b 28 20 22 43 78 |0;.X.RCh|eck( "Cx|
|00002a40| 52 65 61 6c 22 2c 20 43 | 78 52 65 61 6c 28 20 62 |Real", C|xReal( b|
|00002a50| 70 20 29 2c 20 31 2e 30 | 20 29 3b 0a 58 09 52 43 |p ), 1.0| );.X.RC|
|00002a60| 68 65 63 6b 28 20 22 43 | 78 49 6d 61 67 22 2c 20 |heck( "C|xImag", |
|00002a70| 43 78 49 6d 61 67 28 20 | 62 70 20 29 2c 20 32 2e |CxImag( |bp ), 2.|
|00002a80| 30 20 29 3b 0a 58 0a 58 | 09 2f 2a 20 43 78 43 6f |0 );.X.X|./* CxCo|
|00002a90| 6e 73 20 74 65 73 74 20 | 2a 2f 0a 58 09 63 70 20 |ns test |*/.X.cp |
|00002aa0| 3d 20 43 78 43 6f 6e 73 | 28 20 26 61 2c 20 2d 33 |= CxCons|( &a, -3|
|00002ab0| 2e 30 2c 20 2d 34 2e 30 | 29 3b 0a 58 09 43 43 68 |.0, -4.0|);.X.CCh|
|00002ac0| 65 63 6b 28 20 22 43 78 | 43 6f 6e 73 20 31 22 2c |eck( "Cx|Cons 1",|
|00002ad0| 20 61 2c 20 2d 33 2e 30 | 2c 20 2d 34 2e 30 20 29 | a, -3.0|, -4.0 )|
|00002ae0| 3b 0a 58 09 43 43 68 65 | 63 6b 28 20 22 43 78 43 |;.X.CChe|ck( "CxC|
|00002af0| 6f 6e 73 20 32 22 2c 20 | 2a 63 70 2c 20 2d 33 2e |ons 2", |*cp, -3.|
|00002b00| 30 2c 20 2d 34 2e 30 20 | 29 3b 0a 58 0a 58 09 2f |0, -4.0 |);.X.X./|
|00002b10| 2a 20 43 78 4e 65 67 20 | 74 65 73 74 20 2a 2f 0a |* CxNeg |test */.|
|00002b20| 58 09 63 70 20 3d 20 43 | 78 4e 65 67 28 20 26 61 |X.cp = C|xNeg( &a|
|00002b30| 20 29 3b 0a 58 09 43 43 | 68 65 63 6b 28 20 22 43 | );.X.CC|heck( "C|
|00002b40| 78 4e 65 67 20 31 22 2c | 20 61 2c 20 33 2e 30 2c |xNeg 1",| a, 3.0,|
|00002b50| 20 34 2e 30 20 29 3b 0a | 58 09 43 43 68 65 63 6b | 4.0 );.|X.CCheck|
|00002b60| 28 20 22 43 78 4e 65 67 | 20 32 22 2c 20 2a 63 70 |( "CxNeg| 2", *cp|
|00002b70| 2c 20 33 2e 30 2c 20 34 | 2e 30 20 29 3b 0a 58 0a |, 3.0, 4|.0 );.X.|
|00002b80| 58 09 2f 2a 20 43 78 43 | 6f 70 79 20 74 65 73 74 |X./* CxC|opy test|
|00002b90| 20 2a 2f 0a 58 09 63 70 | 20 3d 20 43 78 43 6f 70 | */.X.cp| = CxCop|
|00002ba0| 79 28 20 62 70 2c 20 26 | 61 20 29 3b 0a 58 09 28 |y( bp, &|a );.X.(|
|00002bb0| 76 6f 69 64 29 43 78 43 | 6f 6e 73 28 20 26 61 2c |void)CxC|ons( &a,|
|00002bc0| 20 31 2e 30 2c 20 73 71 | 72 74 28 20 33 2e 30 20 | 1.0, sq|rt( 3.0 |
|00002bd0| 29 20 29 3b 0a 58 09 43 | 43 68 65 63 6b 28 20 22 |) );.X.C|Check( "|
|00002be0| 43 78 43 6f 70 79 20 31 | 22 2c 20 2a 62 70 2c 20 |CxCopy 1|", *bp, |
|00002bf0| 33 2e 30 2c 20 34 2e 30 | 20 29 3b 0a 58 09 43 43 |3.0, 4.0| );.X.CC|
|00002c00| 68 65 63 6b 28 20 22 43 | 78 43 6f 70 79 20 32 22 |heck( "C|xCopy 2"|
|00002c10| 2c 20 2a 63 70 2c 20 33 | 2e 30 2c 20 34 2e 30 20 |, *cp, 3|.0, 4.0 |
|00002c20| 29 3b 0a 58 0a 58 09 2f | 2a 20 43 78 41 6d 70 6c |);.X.X./|* CxAmpl|
|00002c30| 2c 20 43 78 50 68 61 73 | 20 74 65 73 74 20 2a 2f |, CxPhas| test */|
|00002c40| 0a 58 09 52 43 68 65 63 | 6b 28 20 22 43 78 41 6d |.X.RChec|k( "CxAm|
|00002c50| 70 6c 20 31 22 2c 20 43 | 78 41 6d 70 6c 28 20 26 |pl 1", C|xAmpl( &|
|00002c60| 61 20 29 2c 20 32 2e 30 | 20 29 3b 0a 58 09 52 43 |a ), 2.0| );.X.RC|
|00002c70| 68 65 63 6b 28 20 22 43 | 78 50 68 61 73 20 31 22 |heck( "C|xPhas 1"|
|00002c80| 2c 20 43 78 50 68 61 73 | 28 20 26 61 20 29 20 2a |, CxPhas|( &a ) *|
|00002c90| 20 44 45 47 52 41 44 2c | 20 36 30 2e 30 20 29 3b | DEGRAD,| 60.0 );|
|00002ca0| 0a 58 09 2f 2a 20 74 72 | 79 20 6f 74 68 65 72 20 |.X./* tr|y other |
|00002cb0| 71 75 61 64 72 61 6e 74 | 73 20 2a 2f 0a 58 09 61 |quadrant|s */.X.a|
|00002cc0| 2e 72 65 20 3d 20 2d 61 | 2e 72 65 3b 0a 58 09 52 |.re = -a|.re;.X.R|
|00002cd0| 43 68 65 63 6b 28 20 22 | 43 78 41 6d 70 6c 20 32 |Check( "|CxAmpl 2|
|00002ce0| 22 2c 20 43 78 41 6d 70 | 6c 28 20 26 61 20 29 2c |", CxAmp|l( &a ),|
|00002cf0| 20 32 2e 30 20 29 3b 0a | 58 09 52 43 68 65 63 6b | 2.0 );.|X.RCheck|
|00002d00| 28 20 22 43 78 50 68 61 | 73 20 32 22 2c 20 43 78 |( "CxPha|s 2", Cx|
|00002d10| 50 68 61 73 28 20 26 61 | 20 29 20 2a 20 44 45 47 |Phas( &a| ) * DEG|
|00002d20| 52 41 44 2c 20 31 32 30 | 2e 30 20 29 3b 0a 58 09 |RAD, 120|.0 );.X.|
|00002d30| 61 2e 69 6d 20 3d 20 2d | 61 2e 69 6d 3b 0a 58 09 |a.im = -|a.im;.X.|
|00002d40| 52 43 68 65 63 6b 28 20 | 22 43 78 41 6d 70 6c 20 |RCheck( |"CxAmpl |
|00002d50| 33 22 2c 20 43 78 41 6d | 70 6c 28 20 26 61 20 29 |3", CxAm|pl( &a )|
|00002d60| 2c 20 32 2e 30 20 29 3b | 0a 58 09 52 43 68 65 63 |, 2.0 );|.X.RChec|
|00002d70| 6b 28 20 22 43 78 50 68 | 61 73 20 33 22 2c 20 43 |k( "CxPh|as 3", C|
|00002d80| 78 50 68 61 73 28 20 26 | 61 20 29 20 2a 20 44 45 |xPhas( &|a ) * DE|
|00002d90| 47 52 41 44 2c 20 2d 31 | 32 30 2e 30 20 29 3b 0a |GRAD, -1|20.0 );.|
|00002da0| 58 09 61 2e 72 65 20 3d | 20 2d 61 2e 72 65 3b 0a |X.a.re =| -a.re;.|
|00002db0| 58 09 52 43 68 65 63 6b | 28 20 22 43 78 41 6d 70 |X.RCheck|( "CxAmp|
|00002dc0| 6c 20 34 22 2c 20 43 78 | 41 6d 70 6c 28 20 26 61 |l 4", Cx|Ampl( &a|
|00002dd0| 20 29 2c 20 32 2e 30 20 | 29 3b 0a 58 09 52 43 68 | ), 2.0 |);.X.RCh|
|00002de0| 65 63 6b 28 20 22 43 78 | 50 68 61 73 20 34 22 2c |eck( "Cx|Phas 4",|
|00002df0| 20 43 78 50 68 61 73 28 | 20 26 61 20 29 20 2a 20 | CxPhas(| &a ) * |
|00002e00| 44 45 47 52 41 44 2c 20 | 2d 36 30 2e 30 20 29 3b |DEGRAD, |-60.0 );|
|00002e10| 0a 58 09 2f 2a 20 6f 6e | 65 20 6d 6f 72 65 20 66 |.X./* on|e more f|
|00002e20| 6f 72 20 67 6f 6f 64 20 | 6d 65 61 73 75 72 65 20 |or good |measure |
|00002e30| 2a 2f 0a 58 09 52 43 68 | 65 63 6b 28 20 22 43 78 |*/.X.RCh|eck( "Cx|
|00002e40| 41 6d 70 6c 20 35 22 2c | 20 43 78 41 6d 70 6c 28 |Ampl 5",| CxAmpl(|
|00002e50| 20 62 70 20 29 2c 20 35 | 2e 30 20 29 3b 0a 58 0a | bp ), 5|.0 );.X.|
|00002e60| 58 09 2f 2a 20 43 78 50 | 68 73 72 20 74 65 73 74 |X./* CxP|hsr test|
|00002e70| 20 2a 2f 0a 58 09 63 70 | 20 3d 20 43 78 50 68 73 | */.X.cp| = CxPhs|
|00002e80| 72 28 20 26 61 2c 20 31 | 30 30 2e 30 2c 20 2d 32 |r( &a, 1|00.0, -2|
|00002e90| 30 2e 30 20 2f 20 44 45 | 47 52 41 44 20 29 3b 0a |0.0 / DE|GRAD );.|
|00002ea0| 58 09 52 43 68 65 63 6b | 28 20 22 43 78 50 68 73 |X.RCheck|( "CxPhs|
|00002eb0| 72 20 31 22 2c 20 43 78 | 41 6d 70 6c 28 20 26 61 |r 1", Cx|Ampl( &a|
|00002ec0| 20 29 2c 20 31 30 30 2e | 30 20 29 3b 0a 58 09 52 | ), 100.|0 );.X.R|
|00002ed0| 43 68 65 63 6b 28 20 22 | 43 78 50 68 73 72 20 32 |Check( "|CxPhsr 2|
|00002ee0| 22 2c 20 43 78 50 68 61 | 73 28 20 26 61 20 29 20 |", CxPha|s( &a ) |
|00002ef0| 2a 20 44 45 47 52 41 44 | 2c 20 2d 32 30 2e 30 20 |* DEGRAD|, -20.0 |
|00002f00| 29 3b 0a 58 09 52 43 68 | 65 63 6b 28 20 22 43 78 |);.X.RCh|eck( "Cx|
|00002f10| 50 68 73 72 20 33 22 2c | 20 43 78 41 6d 70 6c 28 |Phsr 3",| CxAmpl(|
|00002f20| 20 63 70 20 29 2c 20 31 | 30 30 2e 30 20 29 3b 0a | cp ), 1|00.0 );.|
|00002f30| 58 09 52 43 68 65 63 6b | 28 20 22 43 78 50 68 73 |X.RCheck|( "CxPhs|
|00002f40| 72 20 34 22 2c 20 43 78 | 50 68 61 73 28 20 63 70 |r 4", Cx|Phas( cp|
|00002f50| 20 29 20 2a 20 44 45 47 | 52 41 44 2c 20 2d 32 30 | ) * DEG|RAD, -20|
|00002f60| 2e 30 20 29 3b 0a 58 0a | 58 09 2f 2a 20 43 78 43 |.0 );.X.|X./* CxC|
|00002f70| 6f 6e 6a 20 74 65 73 74 | 20 2a 2f 0a 58 09 63 70 |onj test| */.X.cp|
|00002f80| 20 3d 20 43 78 43 6f 6e | 6a 28 20 62 70 20 29 3b | = CxCon|j( bp );|
|00002f90| 0a 58 09 43 43 68 65 63 | 6b 28 20 22 43 78 43 6f |.X.CChec|k( "CxCo|
|00002fa0| 6e 6a 20 31 22 2c 20 2a | 62 70 2c 20 33 2e 30 2c |nj 1", *|bp, 3.0,|
|00002fb0| 20 2d 34 2e 30 20 29 3b | 0a 58 09 43 43 68 65 63 | -4.0 );|.X.CChec|
|00002fc0| 6b 28 20 22 43 78 43 6f | 6e 6a 20 32 22 2c 20 2a |k( "CxCo|nj 2", *|
|00002fd0| 63 70 2c 20 33 2e 30 2c | 20 2d 34 2e 30 20 29 3b |cp, 3.0,| -4.0 );|
|00002fe0| 0a 58 0a 58 09 2f 2a 20 | 43 78 53 63 61 6c 20 74 |.X.X./* |CxScal t|
|00002ff0| 65 73 74 20 2a 2f 0a 58 | 09 63 70 20 3d 20 43 78 |est */.X|.cp = Cx|
|00003000| 53 63 61 6c 28 20 62 70 | 2c 20 32 2e 30 20 29 3b |Scal( bp|, 2.0 );|
|00003010| 0a 58 09 43 43 68 65 63 | 6b 28 20 22 43 78 53 63 |.X.CChec|k( "CxSc|
|00003020| 61 6c 20 31 22 2c 20 2a | 62 70 2c 20 36 2e 30 2c |al 1", *|bp, 6.0,|
|00003030| 20 2d 38 2e 30 20 29 3b | 0a 58 09 43 43 68 65 63 | -8.0 );|.X.CChec|
|00003040| 6b 28 20 22 43 78 53 63 | 61 6c 20 32 22 2c 20 2a |k( "CxSc|al 2", *|
|00003050| 63 70 2c 20 36 2e 30 2c | 20 2d 38 2e 30 20 29 3b |cp, 6.0,| -8.0 );|
|00003060| 0a 58 0a 58 09 2f 2a 20 | 43 78 41 64 64 20 74 65 |.X.X./* |CxAdd te|
|00003070| 73 74 20 2a 2f 0a 58 09 | 63 70 20 3d 20 43 78 41 |st */.X.|cp = CxA|
|00003080| 64 64 28 20 43 78 43 6f | 6e 73 28 20 26 61 2c 20 |dd( CxCo|ns( &a, |
|00003090| 2d 34 2e 30 2c 20 31 31 | 2e 30 20 29 2c 20 62 70 |-4.0, 11|.0 ), bp|
|000030a0| 20 29 3b 0a 58 09 43 43 | 68 65 63 6b 28 20 22 43 | );.X.CC|heck( "C|
|000030b0| 78 41 64 64 20 31 22 2c | 20 61 2c 20 32 2e 30 2c |xAdd 1",| a, 2.0,|
|000030c0| 20 33 2e 30 20 29 3b 0a | 58 09 43 43 68 65 63 6b | 3.0 );.|X.CCheck|
|000030d0| 28 20 22 43 78 41 64 64 | 20 32 22 2c 20 2a 63 70 |( "CxAdd| 2", *cp|
|000030e0| 2c 20 32 2e 30 2c 20 33 | 2e 30 20 29 3b 0a 58 0a |, 2.0, 3|.0 );.X.|
|000030f0| 58 09 2f 2a 20 43 78 53 | 75 62 20 74 65 73 74 20 |X./* CxS|ub test |
|00003100| 2a 2f 0a 58 09 63 70 20 | 3d 20 43 78 53 75 62 28 |*/.X.cp |= CxSub(|
|00003110| 20 43 78 43 6f 6e 73 28 | 20 26 61 2c 20 34 2e 30 | CxCons(| &a, 4.0|
|00003120| 2c 20 37 2e 30 20 29 2c | 20 62 70 20 29 3b 0a 58 |, 7.0 ),| bp );.X|
|00003130| 09 43 43 68 65 63 6b 28 | 20 22 43 78 53 75 62 20 |.CCheck(| "CxSub |
|00003140| 31 22 2c 20 61 2c 20 2d | 32 2e 30 2c 20 31 35 2e |1", a, -|2.0, 15.|
|00003150| 30 20 29 3b 0a 58 09 43 | 43 68 65 63 6b 28 20 22 |0 );.X.C|Check( "|
|00003160| 43 78 53 75 62 20 32 22 | 2c 20 2a 63 70 2c 20 2d |CxSub 2"|, *cp, -|
|00003170| 32 2e 30 2c 20 31 35 2e | 30 20 29 3b 0a 58 0a 58 |2.0, 15.|0 );.X.X|
|00003180| 09 2f 2a 20 43 78 4d 75 | 6c 20 74 65 73 74 20 2a |./* CxMu|l test *|
|00003190| 2f 0a 58 09 63 70 20 3d | 20 43 78 4d 75 6c 28 20 |/.X.cp =| CxMul( |
|000031a0| 43 78 43 6f 6e 73 28 20 | 62 70 2c 20 2d 31 2e 30 |CxCons( |bp, -1.0|
|000031b0| 2c 20 33 2e 30 20 29 2c | 20 43 78 43 6f 6e 73 28 |, 3.0 ),| CxCons(|
|000031c0| 20 26 61 2c 20 31 2e 30 | 2c 20 32 2e 30 20 29 20 | &a, 1.0|, 2.0 ) |
|000031d0| 29 3b 0a 58 09 43 43 68 | 65 63 6b 28 20 22 43 78 |);.X.CCh|eck( "Cx|
|000031e0| 4d 75 6c 20 31 22 2c 20 | 2a 62 70 2c 20 2d 37 2e |Mul 1", |*bp, -7.|
|000031f0| 30 2c 20 31 2e 30 20 29 | 3b 0a 58 09 43 43 68 65 |0, 1.0 )|;.X.CChe|
|00003200| 63 6b 28 20 22 43 78 4d | 75 6c 20 32 22 2c 20 2a |ck( "CxM|ul 2", *|
|00003210| 63 70 2c 20 2d 37 2e 30 | 2c 20 31 2e 30 20 29 3b |cp, -7.0|, 1.0 );|
|00003220| 0a 58 0a 58 09 2f 2a 20 | 43 78 44 69 76 20 74 65 |.X.X./* |CxDiv te|
|00003230| 73 74 20 2a 2f 0a 58 09 | 63 70 20 3d 20 43 78 44 |st */.X.|cp = CxD|
|00003240| 69 76 28 20 62 70 2c 20 | 26 61 20 29 3b 0a 58 09 |iv( bp, |&a );.X.|
|00003250| 43 43 68 65 63 6b 28 20 | 22 43 78 44 69 76 20 31 |CCheck( |"CxDiv 1|
|00003260| 22 2c 20 2a 62 70 2c 20 | 2d 31 2e 30 2c 20 33 2e |", *bp, |-1.0, 3.|
|00003270| 30 20 29 3b 0a 58 09 43 | 43 68 65 63 6b 28 20 22 |0 );.X.C|Check( "|
|00003280| 43 78 44 69 76 20 32 22 | 2c 20 2a 63 70 2c 20 2d |CxDiv 2"|, *cp, -|
|00003290| 31 2e 30 2c 20 33 2e 30 | 20 29 3b 0a 58 0a 58 09 |1.0, 3.0| );.X.X.|
|000032a0| 2f 2a 20 43 78 53 71 72 | 74 20 61 6e 64 20 6f 76 |/* CxSqr|t and ov|
|000032b0| 65 72 6c 61 70 70 69 6e | 67 20 43 78 4d 75 6c 20 |erlappin|g CxMul |
|000032c0| 74 65 73 74 73 20 2a 2f | 0a 58 09 28 76 6f 69 64 |tests */|.X.(void|
|000032d0| 29 43 78 43 6f 6e 73 28 | 20 26 61 2c 20 2d 31 2e |)CxCons(| &a, -1.|
|000032e0| 30 2c 20 32 2e 30 20 29 | 3b 0a 58 09 63 70 20 3d |0, 2.0 )|;.X.cp =|
|000032f0| 20 43 78 53 71 72 74 28 | 20 43 78 4d 75 6c 28 20 | CxSqrt(| CxMul( |
|00003300| 26 61 2c 20 26 61 20 29 | 20 29 3b 0a 58 09 43 43 |&a, &a )| );.X.CC|
|00003310| 68 65 63 6b 28 20 22 43 | 78 53 71 72 74 20 31 22 |heck( "C|xSqrt 1"|
|00003320| 2c 20 61 2c 20 2d 31 2e | 30 2c 20 32 2e 30 20 29 |, a, -1.|0, 2.0 )|
|00003330| 3b 0a 58 09 43 43 68 65 | 63 6b 28 20 22 43 78 53 |;.X.CChe|ck( "CxS|
|00003340| 71 72 74 20 32 22 2c 20 | 2a 63 70 2c 20 2d 31 2e |qrt 2", |*cp, -1.|
|00003350| 30 2c 20 32 2e 30 20 29 | 3b 0a 58 09 28 76 6f 69 |0, 2.0 )|;.X.(voi|
|00003360| 64 29 43 78 43 6f 6e 73 | 28 20 26 61 2c 20 33 2e |d)CxCons|( &a, 3.|
|00003370| 30 2c 20 32 2e 30 20 29 | 3b 0a 58 09 63 70 20 3d |0, 2.0 )|;.X.cp =|
|00003380| 20 43 78 53 71 72 74 28 | 20 43 78 4d 75 6c 28 20 | CxSqrt(| CxMul( |
|00003390| 26 61 2c 20 26 61 20 29 | 20 29 3b 0a 58 09 43 43 |&a, &a )| );.X.CC|
|000033a0| 68 65 63 6b 28 20 22 43 | 78 53 71 72 74 20 33 22 |heck( "C|xSqrt 3"|
|000033b0| 2c 20 61 2c 20 33 2e 30 | 2c 20 32 2e 30 20 29 3b |, a, 3.0|, 2.0 );|
|000033c0| 0a 58 09 43 43 68 65 63 | 6b 28 20 22 43 78 53 71 |.X.CChec|k( "CxSq|
|000033d0| 72 74 20 34 22 2c 20 2a | 63 70 2c 20 33 2e 30 2c |rt 4", *|cp, 3.0,|
|000033e0| 20 32 2e 30 20 29 3b 0a | 58 0a 58 09 2f 2a 20 43 | 2.0 );.|X.X./* C|
|000033f0| 78 46 72 65 65 20 22 74 | 65 73 74 22 20 2a 2f 0a |xFree "t|est" */.|
|00003400| 58 09 43 78 46 72 65 65 | 28 20 62 70 20 29 3b 0a |X.CxFree|( bp );.|
|00003410| 58 0a 58 09 72 65 74 75 | 72 6e 20 65 72 72 73 3b |X.X.retu|rn errs;|
|00003420| 0a 58 09 7d 0a 58 0a 58 | 0a 58 73 74 61 74 69 63 |.X.}.X.X|.Xstatic|
|00003430| 20 76 6f 69 64 0a 58 43 | 43 68 65 63 6b 28 20 73 | void.XC|Check( s|
|00003440| 2c 20 63 2c 20 72 2c 20 | 69 20 29 09 09 09 2f 2a |, c, r, |i ).../*|
|00003450| 20 63 68 65 63 6b 20 63 | 6f 6d 70 6c 65 78 20 6e | check c|omplex n|
|00003460| 75 6d 62 65 72 20 2a 2f | 0a 58 09 63 68 61 72 09 |umber */|.X.char.|
|00003470| 2a 73 3b 09 09 09 2f 2a | 20 6d 65 73 73 61 67 65 |*s;.../*| message|
|00003480| 20 73 74 72 69 6e 67 20 | 66 6f 72 20 66 61 69 6c | string |for fail|
|00003490| 75 72 65 20 2a 2f 0a 58 | 09 63 6f 6d 70 6c 65 78 |ure */.X|.complex|
|000034a0| 09 63 3b 09 09 09 2f 2a | 20 63 6f 6d 70 6c 65 78 |.c;.../*| complex|
|000034b0| 20 74 6f 20 62 65 20 63 | 68 65 63 6b 65 64 20 2a | to be c|hecked *|
|000034c0| 2f 0a 58 09 64 6f 75 62 | 6c 65 09 72 2c 20 69 3b |/.X.doub|le.r, i;|
|000034d0| 09 09 09 2f 2a 20 65 78 | 70 65 63 74 65 64 20 72 |.../* ex|pected r|
|000034e0| 65 61 6c 2c 20 69 6d 61 | 67 69 6e 61 72 79 20 70 |eal, ima|ginary p|
|000034f0| 61 72 74 73 20 2a 2f 0a | 58 09 7b 0a 58 09 69 66 |arts */.|X.{.X.if|
|00003500| 20 28 20 52 65 6c 44 69 | 66 28 20 43 78 52 65 61 | ( RelDi|f( CxRea|
|00003510| 6c 28 20 26 63 20 29 2c | 20 72 20 29 20 3e 20 54 |l( &c ),| r ) > T|
|00003520| 4f 4c 0a 58 09 20 20 7c | 7c 20 52 65 6c 44 69 66 |OL.X.  ||| RelDif|
|00003530| 28 20 43 78 49 6d 61 67 | 28 20 26 63 20 29 2c 20 |( CxImag|( &c ), |
|00003540| 69 20 29 20 3e 20 54 4f | 4c 0a 58 09 20 20 20 29 |i ) > TO|L.X.   )|
|00003550| 09 7b 0a 58 09 09 2b 2b | 65 72 72 73 3b 0a 58 09 |.{.X..++|errs;.X.|
|00003560| 09 50 72 69 6e 74 66 28 | 20 22 25 73 3b 20 73 2e |.Printf(| "%s; s.|
|00003570| 62 2e 20 28 25 66 2c 25 | 66 29 2c 20 77 61 73 20 |b. (%f,%|f), was |
|00003580| 28 25 67 2c 25 67 29 5c | 6e 22 2c 0a 58 09 09 09 |(%g,%g)\|n",.X...|
|00003590| 73 2c 20 72 2c 20 69 2c | 20 63 2e 72 65 2c 20 63 |s, r, i,| c.re, c|
|000035a0| 2e 69 6d 0a 58 09 09 20 | 20 20 20 20 20 29 3b 0a |.im.X.. |     );.|
|000035b0| 58 09 09 7d 0a 58 09 7d | 0a 58 0a 58 0a 58 73 74 |X..}.X.}|.X.X.Xst|
|000035c0| 61 74 69 63 20 76 6f 69 | 64 0a 58 52 43 68 65 63 |atic voi|d.XRChec|
|000035d0| 6b 28 20 73 2c 20 64 2c | 20 72 20 29 09 09 09 2f |k( s, d,| r ).../|
|000035e0| 2a 20 63 68 65 63 6b 20 | 72 65 61 6c 20 6e 75 6d |* check |real num|
|000035f0| 62 65 72 20 2a 2f 0a 58 | 09 63 68 61 72 09 2a 73 |ber */.X|.char.*s|
|00003600| 3b 09 09 09 2f 2a 20 6d | 65 73 73 61 67 65 20 73 |;.../* m|essage s|
|00003610| 74 72 69 6e 67 20 66 6f | 72 20 66 61 69 6c 75 72 |tring fo|r failur|
|00003620| 65 20 2a 2f 0a 58 09 64 | 6f 75 62 6c 65 09 64 3b |e */.X.d|ouble.d;|
|00003630| 09 09 09 2f 2a 20 72 65 | 61 6c 20 74 6f 20 62 65 |.../* re|al to be|
|00003640| 20 63 68 65 63 6b 65 64 | 20 2a 2f 0a 58 09 64 6f | checked| */.X.do|
|00003650| 75 62 6c 65 09 72 3b 09 | 09 09 2f 2a 20 65 78 70 |uble.r;.|../* exp|
|00003660| 65 63 74 65 64 20 76 61 | 6c 75 65 20 2a 2f 0a 58 |ected va|lue */.X|
|00003670| 09 7b 0a 58 09 69 66 20 | 28 20 52 65 6c 44 69 66 |.{.X.if |( RelDif|
|00003680| 28 20 64 2c 20 72 20 29 | 20 3e 20 54 4f 4c 20 29 |( d, r )| > TOL )|
|00003690| 0a 58 09 09 7b 0a 58 09 | 09 2b 2b 65 72 72 73 3b |.X..{.X.|.++errs;|
|000036a0| 0a 58 09 09 50 72 69 6e | 74 66 28 20 22 25 73 3b |.X..Prin|tf( "%s;|
|000036b0| 20 73 2e 62 2e 20 25 66 | 2c 20 77 61 73 20 25 67 | s.b. %f|, was %g|
|000036c0| 5c 6e 22 2c 20 73 2c 20 | 72 2c 20 64 20 29 3b 0a |\n", s, |r, d );.|
|000036d0| 58 09 09 7d 0a 58 09 7d | 0a 58 0a 58 0a 58 73 74 |X..}.X.}|.X.X.Xst|
|000036e0| 61 74 69 63 20 64 6f 75 | 62 6c 65 0a 58 52 65 6c |atic dou|ble.XRel|
|000036f0| 44 69 66 28 20 61 2c 20 | 62 20 29 09 09 09 2f 2a |Dif( a, |b ).../*|
|00003700| 20 72 65 74 75 72 6e 73 | 20 72 65 6c 61 74 69 76 | returns| relativ|
|00003710| 65 20 64 69 66 66 65 72 | 65 6e 63 65 3a 09 2a 2f |e differ|ence:.*/|
|00003720| 0a 58 09 64 6f 75 62 6c | 65 09 61 2c 20 62 3b 09 |.X.doubl|e.a, b;.|
|00003730| 09 2f 2a 20 30 2e 30 20 | 69 66 20 65 78 61 63 74 |./* 0.0 |if exact|
|00003740| 6c 79 20 74 68 65 20 73 | 61 6d 65 2c 0a 58 09 09 |ly the s|ame,.X..|
|00003750| 09 09 20 20 20 6f 74 68 | 65 72 77 69 73 65 20 72 |..   oth|erwise r|
|00003760| 61 74 69 6f 20 6f 66 20 | 64 69 66 66 65 72 65 6e |atio of |differen|
|00003770| 63 65 0a 58 09 09 09 09 | 20 20 20 74 6f 20 74 68 |ce.X....|   to th|
|00003780| 65 20 6c 61 72 67 65 72 | 20 6f 66 20 74 68 65 20 |e larger| of the |
|00003790| 74 77 6f 09 2a 2f 0a 58 | 09 7b 0a 58 09 64 6f 75 |two.*/.X|.{.X.dou|
|000037a0| 62 6c 65 09 63 20 3d 20 | 41 62 73 28 20 61 20 29 |ble.c = |Abs( a )|
|000037b0| 3b 0a 58 09 64 6f 75 62 | 6c 65 09 64 20 3d 20 41 |;.X.doub|le.d = A|
|000037c0| 62 73 28 20 62 20 29 3b | 0a 58 0a 58 09 64 20 3d |bs( b );|.X.X.d =|
|000037d0| 20 4d 61 78 28 20 63 2c | 20 64 20 29 3b 0a 58 0a | Max( c,| d );.X.|
|000037e0| 58 09 72 65 74 75 72 6e | 20 64 20 3d 3d 20 30 2e |X.return| d == 0.|
|000037f0| 30 20 3f 20 30 2e 30 20 | 3a 20 41 62 73 28 20 61 |0 ? 0.0 |: Abs( a|
|00003800| 20 2d 20 62 20 29 20 2f | 20 64 3b 0a 58 09 7d 0a | - b ) /| d;.X.}.|
|00003810| 45 4e 44 5f 4f 46 5f 63 | 78 5f 74 65 73 74 2e 63 |END_OF_c|x_test.c|
|00003820| 0a 69 66 20 74 65 73 74 | 20 34 32 35 30 20 2d 6e |.if test| 4250 -n|
|00003830| 65 20 60 77 63 20 2d 63 | 20 3c 63 78 5f 74 65 73 |e `wc -c| <cx_tes|
|00003840| 74 2e 63 60 3b 20 74 68 | 65 6e 0a 20 20 20 20 65 |t.c`; th|en.    e|
|00003850| 63 68 6f 20 73 68 61 72 | 3a 20 5c 22 63 78 5f 74 |cho shar|: \"cx_t|
|00003860| 65 73 74 2e 63 5c 22 20 | 75 6e 70 61 63 6b 65 64 |est.c\" |unpacked|
|00003870| 20 77 69 74 68 20 77 72 | 6f 6e 67 20 73 69 7a 65 | with wr|ong size|
|00003880| 21 0a 66 69 0a 23 20 65 | 6e 64 20 6f 66 20 6f 76 |!.fi.# e|nd of ov|
|00003890| 65 72 77 72 69 74 69 6e | 67 20 63 68 65 63 6b 0a |erwritin|g check.|
|000038a0| 66 69 0a 69 66 20 74 65 | 73 74 20 2d 66 20 63 78 |fi.if te|st -f cx|
|000038b0| 61 64 64 2e 63 20 2d 61 | 20 22 24 7b 31 7d 22 20 |add.c -a| "${1}" |
|000038c0| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|000038d0| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c |  echo s|har: Wil|
|000038e0| 6c 20 6e 6f 74 20 6f 76 | 65 72 2d 77 72 69 74 65 |l not ov|er-write|
|000038f0| 20 65 78 69 73 74 69 6e | 67 20 66 69 6c 65 20 5c | existin|g file \|
|00003900| 22 63 78 61 64 64 2e 63 | 5c 22 0a 65 6c 73 65 0a |"cxadd.c|\".else.|
|00003910| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|00003920| 63 74 69 6e 67 20 5c 22 | 63 78 61 64 64 2e 63 5c |cting \"|cxadd.c\|
|00003930| 22 20 5c 28 33 30 34 20 | 63 68 61 72 61 63 74 65 |" \(304 |characte|
|00003940| 72 73 5c 29 0a 73 65 64 | 20 22 73 2f 5e 58 2f 2f |rs\).sed| "s/^X//|
|00003950| 22 20 3e 63 78 61 64 64 | 2e 63 20 3c 3c 27 45 4e |" >cxadd|.c <<'EN|
|00003960| 44 5f 4f 46 5f 63 78 61 | 64 64 2e 63 27 0a 58 2f |D_OF_cxa|dd.c'.X/|
|00003970| 2a 0a 58 09 43 78 41 64 | 64 20 2d 2d 20 61 64 64 |*.X.CxAd|d -- add|
|00003980| 20 6f 6e 65 20 63 6f 6d | 70 6c 65 78 20 74 6f 20 | one com|plex to |
|00003990| 61 6e 6f 74 68 65 72 0a | 58 0a 58 09 6c 61 73 74 |another.|X.X.last|
|000039a0| 20 65 64 69 74 3a 09 38 | 36 2f 30 31 2f 30 34 09 | edit:.8|6/01/04.|
|000039b0| 44 20 41 20 47 77 79 6e | 0a 58 0a 58 09 53 43 43 |D A Gwyn|.X.X.SCC|
|000039c0| 53 20 49 44 3a 09 40 28 | 23 29 63 78 61 64 64 2e |S ID:.@(|#)cxadd.|
|000039d0| 63 09 31 2e 31 0a 58 0a | 58 09 43 78 41 64 64 28 |c.1.1.X.|X.CxAdd(|
|000039e0| 20 26 61 2c 20 26 62 20 | 29 09 61 64 64 73 20 20 | &a, &b |).adds  |
|000039f0| 62 20 20 74 6f 20 20 61 | 20 20 61 6e 64 20 72 65 |b  to  a|  and re|
|00003a00| 74 75 72 6e 73 20 20 26 | 61 0a 58 2a 2f 0a 58 0a |turns  &|a.X*/.X.|
|00003a10| 58 23 69 6e 63 6c 75 64 | 65 09 3c 63 6f 6d 70 6c |X#includ|e.<compl|
|00003a20| 65 78 2e 68 3e 0a 58 0a | 58 63 6f 6d 70 6c 65 78 |ex.h>.X.|Xcomplex|
|00003a30| 20 2a 0a 58 43 78 41 64 | 64 28 20 61 70 2c 20 62 | *.XCxAd|d( ap, b|
|00003a40| 70 20 29 0a 58 09 72 65 | 67 69 73 74 65 72 20 63 |p ).X.re|gister c|
|00003a50| 6f 6d 70 6c 65 78 09 2a | 61 70 2c 20 2a 62 70 3b |omplex.*|ap, *bp;|
|00003a60| 09 2f 2a 20 6d 61 79 20 | 63 6f 69 6e 63 69 64 65 |./* may |coincide|
|00003a70| 20 2a 2f 0a 58 09 7b 0a | 58 09 61 70 2d 3e 72 65 | */.X.{.|X.ap->re|
|00003a80| 20 2b 3d 20 62 70 2d 3e | 72 65 3b 0a 58 09 61 70 | += bp->|re;.X.ap|
|00003a90| 2d 3e 69 6d 20 2b 3d 20 | 62 70 2d 3e 69 6d 3b 0a |->im += |bp->im;.|
|00003aa0| 58 0a 58 09 72 65 74 75 | 72 6e 20 61 70 3b 0a 58 |X.X.retu|rn ap;.X|
|00003ab0| 09 7d 0a 45 4e 44 5f 4f | 46 5f 63 78 61 64 64 2e |.}.END_O|F_cxadd.|
|00003ac0| 63 0a 69 66 20 74 65 73 | 74 20 33 30 34 20 2d 6e |c.if tes|t 304 -n|
|00003ad0| 65 20 60 77 63 20 2d 63 | 20 3c 63 78 61 64 64 2e |e `wc -c| <cxadd.|
|00003ae0| 63 60 3b 20 74 68 65 6e | 0a 20 20 20 20 65 63 68 |c`; then|.    ech|
|00003af0| 6f 20 73 68 61 72 3a 20 | 5c 22 63 78 61 64 64 2e |o shar: |\"cxadd.|
|00003b00| 63 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |c\" unpa|cked wit|
|00003b10| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 66 69 |h wrong |size!.fi|
|00003b20| 0a 23 20 65 6e 64 20 6f | 66 20 6f 76 65 72 77 72 |.# end o|f overwr|
|00003b30| 69 74 69 6e 67 20 63 68 | 65 63 6b 0a 66 69 0a 69 |iting ch|eck.fi.i|
|00003b40| 66 20 74 65 73 74 20 2d | 66 20 63 78 61 6d 70 6c |f test -|f cxampl|
|00003b50| 2e 63 20 2d 61 20 22 24 | 7b 31 7d 22 20 21 3d 20 |.c -a "$|{1}" != |
|00003b60| 22 2d 63 22 20 3b 20 74 | 68 65 6e 20 0a 20 20 65 |"-c" ; t|hen .  e|
|00003b70| 63 68 6f 20 73 68 61 72 | 3a 20 57 69 6c 6c 20 6e |cho shar|: Will n|
|00003b80| 6f 74 20 6f 76 65 72 2d | 77 72 69 74 65 20 65 78 |ot over-|write ex|
|00003b90| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 63 78 |isting f|ile \"cx|
|00003ba0| 61 6d 70 6c 2e 63 5c 22 | 0a 65 6c 73 65 0a 65 63 |ampl.c\"|.else.ec|
|00003bb0| 68 6f 20 73 68 61 72 3a | 20 45 78 74 72 61 63 74 |ho shar:| Extract|
|00003bc0| 69 6e 67 20 5c 22 63 78 | 61 6d 70 6c 2e 63 5c 22 |ing \"cx|ampl.c\"|
|00003bd0| 20 5c 28 32 37 38 20 63 | 68 61 72 61 63 74 65 72 | \(278 c|haracter|
|00003be0| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|00003bf0| 20 3e 63 78 61 6d 70 6c | 2e 63 20 3c 3c 27 45 4e | >cxampl|.c <<'EN|
|00003c00| 44 5f 4f 46 5f 63 78 61 | 6d 70 6c 2e 63 27 0a 58 |D_OF_cxa|mpl.c'.X|
|00003c10| 2f 2a 0a 58 09 43 78 41 | 6d 70 6c 20 2d 2d 20 61 |/*.X.CxA|mpl -- a|
|00003c20| 6d 70 6c 69 74 75 64 65 | 20 28 6d 61 67 6e 69 74 |mplitude| (magnit|
|00003c30| 75 64 65 2c 20 6d 6f 64 | 75 6c 75 73 2c 20 6e 6f |ude, mod|ulus, no|
|00003c40| 72 6d 29 20 6f 66 20 61 | 20 63 6f 6d 70 6c 65 78 |rm) of a| complex|
|00003c50| 0a 58 0a 58 09 43 78 41 | 6d 70 6c 28 20 26 63 20 |.X.X.CxA|mpl( &c |
|00003c60| 29 09 72 65 74 75 72 6e | 73 20 20 7c 63 7c 0a 58 |).return|s  |c|.X|
|00003c70| 0a 58 09 6c 61 73 74 20 | 65 64 69 74 3a 09 38 36 |.X.last |edit:.86|
|00003c80| 2f 30 31 2f 30 34 09 44 | 20 41 20 47 77 79 6e 0a |/01/04.D| A Gwyn.|
|00003c90| 58 0a 58 09 53 43 43 53 | 20 49 44 3a 09 40 28 23 |X.X.SCCS| ID:.@(#|
|00003ca0| 29 63 78 61 6d 70 6c 2e | 63 09 31 2e 31 0a 58 2a |)cxampl.|c.1.1.X*|
|00003cb0| 2f 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 09 3c 6d |/.X.X#in|clude.<m|
|00003cc0| 61 74 68 2e 68 3e 0a 58 | 0a 58 23 69 6e 63 6c 75 |ath.h>.X|.X#inclu|
|00003cd0| 64 65 09 3c 63 6f 6d 70 | 6c 65 78 2e 68 3e 0a 58 |de.<comp|lex.h>.X|
|00003ce0| 0a 58 64 6f 75 62 6c 65 | 0a 58 43 78 41 6d 70 6c |.Xdouble|.XCxAmpl|
|00003cf0| 28 20 63 70 20 29 0a 58 | 09 72 65 67 69 73 74 65 |( cp ).X|.registe|
|00003d00| 72 20 63 6f 6d 70 6c 65 | 78 09 2a 63 70 3b 0a 58 |r comple|x.*cp;.X|
|00003d10| 09 7b 0a 58 09 72 65 74 | 75 72 6e 20 68 79 70 6f |.{.X.ret|urn hypo|
|00003d20| 74 28 20 63 70 2d 3e 72 | 65 2c 20 63 70 2d 3e 69 |t( cp->r|e, cp->i|
|00003d30| 6d 20 29 3b 0a 58 09 7d | 0a 45 4e 44 5f 4f 46 5f |m );.X.}|.END_OF_|
|00003d40| 63 78 61 6d 70 6c 2e 63 | 0a 69 66 20 74 65 73 74 |cxampl.c|.if test|
|00003d50| 20 32 37 38 20 2d 6e 65 | 20 60 77 63 20 2d 63 20 | 278 -ne| `wc -c |
|00003d60| 3c 63 78 61 6d 70 6c 2e | 63 60 3b 20 74 68 65 6e |<cxampl.|c`; then|
|00003d70| 0a 20 20 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |.    ech|o shar: |
|00003d80| 5c 22 63 78 61 6d 70 6c | 2e 63 5c 22 20 75 6e 70 |\"cxampl|.c\" unp|
|00003d90| 61 63 6b 65 64 20 77 69 | 74 68 20 77 72 6f 6e 67 |acked wi|th wrong|
|00003da0| 20 73 69 7a 65 21 0a 66 | 69 0a 23 20 65 6e 64 20 | size!.f|i.# end |
|00003db0| 6f 66 20 6f 76 65 72 77 | 72 69 74 69 6e 67 20 63 |of overw|riting c|
|00003dc0| 68 65 63 6b 0a 66 69 0a | 69 66 20 74 65 73 74 20 |heck.fi.|if test |
|00003dd0| 2d 66 20 63 78 63 6f 6e | 6a 2e 63 20 2d 61 20 22 |-f cxcon|j.c -a "|
|00003de0| 24 7b 31 7d 22 20 21 3d | 20 22 2d 63 22 20 3b 20 |${1}" !=| "-c" ; |
|00003df0| 74 68 65 6e 20 0a 20 20 | 65 63 68 6f 20 73 68 61 |then .  |echo sha|
|00003e00| 72 3a 20 57 69 6c 6c 20 | 6e 6f 74 20 6f 76 65 72 |r: Will |not over|
|00003e10| 2d 77 72 69 74 65 20 65 | 78 69 73 74 69 6e 67 20 |-write e|xisting |
|00003e20| 66 69 6c 65 20 5c 22 63 | 78 63 6f 6e 6a 2e 63 5c |file \"c|xconj.c\|
|00003e30| 22 0a 65 6c 73 65 0a 65 | 63 68 6f 20 73 68 61 72 |".else.e|cho shar|
|00003e40| 3a 20 45 78 74 72 61 63 | 74 69 6e 67 20 5c 22 63 |: Extrac|ting \"c|
|00003e50| 78 63 6f 6e 6a 2e 63 5c | 22 20 5c 28 32 37 38 20 |xconj.c\|" \(278 |
|00003e60| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00003e70| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 63 78 63 6f 6e | "s/^X//|" >cxcon|
|00003e80| 6a 2e 63 20 3c 3c 27 45 | 4e 44 5f 4f 46 5f 63 78 |j.c <<'E|ND_OF_cx|
|00003e90| 63 6f 6e 6a 2e 63 27 0a | 58 2f 2a 0a 58 09 43 78 |conj.c'.|X/*.X.Cx|
|00003ea0| 43 6f 6e 6a 20 2d 2d 20 | 63 6f 6e 6a 75 67 61 74 |Conj -- |conjugat|
|00003eb0| 65 20 61 20 63 6f 6d 70 | 6c 65 78 0a 58 0a 58 09 |e a comp|lex.X.X.|
|00003ec0| 43 78 43 6f 6e 6a 28 20 | 26 63 20 29 09 63 6f 6e |CxConj( |&c ).con|
|00003ed0| 6a 75 67 61 74 65 73 20 | 20 63 20 20 61 6e 64 20 |jugates | c  and |
|00003ee0| 72 65 74 75 72 6e 73 20 | 20 26 63 0a 58 0a 58 09 |returns | &c.X.X.|
|00003ef0| 6c 61 73 74 20 65 64 69 | 74 3a 09 38 36 2f 30 31 |last edi|t:.86/01|
|00003f00| 2f 30 34 09 44 20 41 20 | 47 77 79 6e 0a 58 0a 58 |/04.D A |Gwyn.X.X|
|00003f10| 09 53 43 43 53 20 49 44 | 3a 09 40 28 23 29 63 78 |.SCCS ID|:.@(#)cx|
|00003f20| 63 6f 6e 6a 2e 63 09 31 | 2e 31 0a 58 2a 2f 0a 58 |conj.c.1|.1.X*/.X|
|00003f30| 0a 58 23 69 6e 63 6c 75 | 64 65 09 3c 63 6f 6d 70 |.X#inclu|de.<comp|
|00003f40| 6c 65 78 2e 68 3e 0a 58 | 0a 58 63 6f 6d 70 6c 65 |lex.h>.X|.Xcomple|
|00003f50| 78 20 2a 0a 58 43 78 43 | 6f 6e 6a 28 20 63 70 20 |x *.XCxC|onj( cp |
|00003f60| 29 0a 58 09 72 65 67 69 | 73 74 65 72 20 63 6f 6d |).X.regi|ster com|
|00003f70| 70 6c 65 78 09 2a 63 70 | 3b 0a 58 09 7b 0a 58 09 |plex.*cp|;.X.{.X.|
|00003f80| 2f 2a 20 28 72 65 61 6c | 20 70 61 72 74 20 75 6e |/* (real| part un|
|00003f90| 63 68 61 6e 67 65 64 29 | 20 2a 2f 0a 58 09 63 70 |changed)| */.X.cp|
|00003fa0| 2d 3e 69 6d 20 3d 20 2d | 63 70 2d 3e 69 6d 3b 0a |->im = -|cp->im;.|
|00003fb0| 58 0a 58 09 72 65 74 75 | 72 6e 20 63 70 3b 0a 58 |X.X.retu|rn cp;.X|
|00003fc0| 09 7d 0a 45 4e 44 5f 4f | 46 5f 63 78 63 6f 6e 6a |.}.END_O|F_cxconj|
|00003fd0| 2e 63 0a 69 66 20 74 65 | 73 74 20 32 37 38 20 2d |.c.if te|st 278 -|
|00003fe0| 6e 65 20 60 77 63 20 2d | 63 20 3c 63 78 63 6f 6e |ne `wc -|c <cxcon|
|00003ff0| 6a 2e 63 60 3b 20 74 68 | 65 6e 0a 20 20 20 20 65 |j.c`; th|en.    e|
|00004000| 63 68 6f 20 73 68 61 72 | 3a 20 5c 22 63 78 63 6f |cho shar|: \"cxco|
|00004010| 6e 6a 2e 63 5c 22 20 75 | 6e 70 61 63 6b 65 64 20 |nj.c\" u|npacked |
|00004020| 77 69 74 68 20 77 72 6f | 6e 67 20 73 69 7a 65 21 |with wro|ng size!|
|00004030| 0a 66 69 0a 23 20 65 6e | 64 20 6f 66 20 6f 76 65 |.fi.# en|d of ove|
|00004040| 72 77 72 69 74 69 6e 67 | 20 63 68 65 63 6b 0a 66 |rwriting| check.f|
|00004050| 69 0a 69 66 20 74 65 73 | 74 20 2d 66 20 63 78 63 |i.if tes|t -f cxc|
|00004060| 6f 6e 73 2e 63 20 2d 61 | 20 22 24 7b 31 7d 22 20 |ons.c -a| "${1}" |
|00004070| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|00004080| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c |  echo s|har: Wil|
|00004090| 6c 20 6e 6f 74 20 6f 76 | 65 72 2d 77 72 69 74 65 |l not ov|er-write|
|000040a0| 20 65 78 69 73 74 69 6e | 67 20 66 69 6c 65 20 5c | existin|g file \|
|000040b0| 22 63 78 63 6f 6e 73 2e | 63 5c 22 0a 65 6c 73 65 |"cxcons.|c\".else|
|000040c0| 0a 65 63 68 6f 20 73 68 | 61 72 3a 20 45 78 74 72 |.echo sh|ar: Extr|
|000040d0| 61 63 74 69 6e 67 20 5c | 22 63 78 63 6f 6e 73 2e |acting \|"cxcons.|
|000040e0| 63 5c 22 20 5c 28 33 32 | 39 20 63 68 61 72 61 63 |c\" \(32|9 charac|
|000040f0| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00004100| 2f 2f 22 20 3e 63 78 63 | 6f 6e 73 2e 63 20 3c 3c |//" >cxc|ons.c <<|
|00004110| 27 45 4e 44 5f 4f 46 5f | 63 78 63 6f 6e 73 2e 63 |'END_OF_|cxcons.c|
|00004120| 27 0a 58 2f 2a 0a 58 09 | 43 78 43 6f 6e 73 20 2d |'.X/*.X.|CxCons -|
|00004130| 2d 20 63 6f 6e 73 74 72 | 75 63 74 20 61 20 63 6f |- constr|uct a co|
|00004140| 6d 70 6c 65 78 20 66 72 | 6f 6d 20 72 65 61 6c 20 |mplex fr|om real |
|00004150| 61 6e 64 20 69 6d 61 67 | 69 6e 61 72 79 20 70 61 |and imag|inary pa|
|00004160| 72 74 73 0a 58 0a 58 09 | 43 78 43 6f 6e 73 28 20 |rts.X.X.|CxCons( |
|00004170| 26 63 2c 20 72 65 2c 20 | 69 6d 20 29 09 6d 61 6b |&c, re, |im ).mak|
|00004180| 65 73 20 20 63 20 3d 20 | 72 65 20 2b 20 69 20 69 |es  c = |re + i i|
|00004190| 6d 20 20 61 6e 64 20 72 | 65 74 75 72 6e 73 20 20 |m  and r|eturns  |
|000041a0| 26 63 0a 58 0a 58 09 6c | 61 73 74 20 65 64 69 74 |&c.X.X.l|ast edit|
|000041b0| 3a 09 38 36 2f 30 31 2f | 30 34 09 44 20 41 20 47 |:.86/01/|04.D A G|
|000041c0| 77 79 6e 0a 58 0a 58 09 | 53 43 43 53 20 49 44 3a |wyn.X.X.|SCCS ID:|
|000041d0| 09 40 28 23 29 63 78 63 | 6f 6e 73 2e 63 09 31 2e |.@(#)cxc|ons.c.1.|
|000041e0| 31 0a 58 2a 2f 0a 58 0a | 58 23 69 6e 63 6c 75 64 |1.X*/.X.|X#includ|
|000041f0| 65 09 3c 63 6f 6d 70 6c | 65 78 2e 68 3e 0a 58 0a |e.<compl|ex.h>.X.|
|00004200| 58 63 6f 6d 70 6c 65 78 | 20 2a 0a 58 43 78 43 6f |Xcomplex| *.XCxCo|
|00004210| 6e 73 28 20 63 70 2c 20 | 72 65 2c 20 69 6d 20 29 |ns( cp, |re, im )|
|00004220| 0a 58 09 72 65 67 69 73 | 74 65 72 20 63 6f 6d 70 |.X.regis|ter comp|
|00004230| 6c 65 78 09 2a 63 70 3b | 0a 58 09 64 6f 75 62 6c |lex.*cp;|.X.doubl|
|00004240| 65 09 09 09 72 65 2c 20 | 69 6d 3b 0a 58 09 7b 0a |e...re, |im;.X.{.|
|00004250| 58 09 63 70 2d 3e 72 65 | 20 3d 20 72 65 3b 0a 58 |X.cp->re| = re;.X|
|00004260| 09 63 70 2d 3e 69 6d 20 | 3d 20 69 6d 3b 0a 58 0a |.cp->im |= im;.X.|
|00004270| 58 09 72 65 74 75 72 6e | 20 63 70 3b 0a 58 09 7d |X.return| cp;.X.}|
|00004280| 0a 45 4e 44 5f 4f 46 5f | 63 78 63 6f 6e 73 2e 63 |.END_OF_|cxcons.c|
|00004290| 0a 69 66 20 74 65 73 74 | 20 33 32 39 20 2d 6e 65 |.if test| 329 -ne|
|000042a0| 20 60 77 63 20 2d 63 20 | 3c 63 78 63 6f 6e 73 2e | `wc -c |<cxcons.|
|000042b0| 63 60 3b 20 74 68 65 6e | 0a 20 20 20 20 65 63 68 |c`; then|.    ech|
|000042c0| 6f 20 73 68 61 72 3a 20 | 5c 22 63 78 63 6f 6e 73 |o shar: |\"cxcons|
|000042d0| 2e 63 5c 22 20 75 6e 70 | 61 63 6b 65 64 20 77 69 |.c\" unp|acked wi|
|000042e0| 74 68 20 77 72 6f 6e 67 | 20 73 69 7a 65 21 0a 66 |th wrong| size!.f|
|000042f0| 69 0a 23 20 65 6e 64 20 | 6f 66 20 6f 76 65 72 77 |i.# end |of overw|
|00004300| 72 69 74 69 6e 67 20 63 | 68 65 63 6b 0a 66 69 0a |riting c|heck.fi.|
|00004310| 69 66 20 74 65 73 74 20 | 2d 66 20 63 78 63 6f 70 |if test |-f cxcop|
|00004320| 79 2e 63 20 2d 61 20 22 | 24 7b 31 7d 22 20 21 3d |y.c -a "|${1}" !=|
|00004330| 20 22 2d 63 22 20 3b 20 | 74 68 65 6e 20 0a 20 20 | "-c" ; |then .  |
|00004340| 65 63 68 6f 20 73 68 61 | 72 3a 20 57 69 6c 6c 20 |echo sha|r: Will |
|00004350| 6e 6f 74 20 6f 76 65 72 | 2d 77 72 69 74 65 20 65 |not over|-write e|
|00004360| 78 69 73 74 69 6e 67 20 | 66 69 6c 65 20 5c 22 63 |xisting |file \"c|
|00004370| 78 63 6f 70 79 2e 63 5c | 22 0a 65 6c 73 65 0a 65 |xcopy.c\|".else.e|
|00004380| 63 68 6f 20 73 68 61 72 | 3a 20 45 78 74 72 61 63 |cho shar|: Extrac|
|00004390| 74 69 6e 67 20 5c 22 63 | 78 63 6f 70 79 2e 63 5c |ting \"c|xcopy.c\|
|000043a0| 22 20 5c 28 32 36 34 20 | 63 68 61 72 61 63 74 65 |" \(264 |characte|
|000043b0| 72 73 5c 29 0a 73 65 64 | 20 22 73 2f 5e 58 2f 2f |rs\).sed| "s/^X//|
|000043c0| 22 20 3e 63 78 63 6f 70 | 79 2e 63 20 3c 3c 27 45 |" >cxcop|y.c <<'E|
|000043d0| 4e 44 5f 4f 46 5f 63 78 | 63 6f 70 79 2e 63 27 0a |ND_OF_cx|copy.c'.|
|000043e0| 58 2f 2a 0a 58 09 43 78 | 43 6f 70 79 20 2d 2d 20 |X/*.X.Cx|Copy -- |
|000043f0| 63 6f 70 79 20 61 20 63 | 6f 6d 70 6c 65 78 0a 58 |copy a c|omplex.X|
|00004400| 0a 58 09 6c 61 73 74 20 | 65 64 69 74 3a 09 38 36 |.X.last |edit:.86|
|00004410| 2f 30 31 2f 30 34 09 44 | 20 41 20 47 77 79 6e 0a |/01/04.D| A Gwyn.|
|00004420| 58 0a 58 09 53 43 43 53 | 20 49 44 3a 09 40 28 23 |X.X.SCCS| ID:.@(#|
|00004430| 29 63 78 63 6f 70 79 2e | 63 09 31 2e 31 0a 58 0a |)cxcopy.|c.1.1.X.|
|00004440| 58 09 43 78 43 6f 70 79 | 28 20 26 61 2c 20 26 62 |X.CxCopy|( &a, &b|
|00004450| 20 29 09 63 6f 70 69 65 | 73 20 20 62 20 20 74 6f | ).copie|s  b  to|
|00004460| 20 20 61 20 20 61 6e 64 | 20 72 65 74 75 72 6e 73 |  a  and| returns|
|00004470| 20 20 26 61 0a 58 2a 2f | 0a 58 0a 58 23 69 6e 63 |  &a.X*/|.X.X#inc|
|00004480| 6c 75 64 65 09 3c 63 6f | 6d 70 6c 65 78 2e 68 3e |lude.<co|mplex.h>|
|00004490| 0a 58 0a 58 63 6f 6d 70 | 6c 65 78 20 2a 0a 58 43 |.X.Xcomp|lex *.XC|
|000044a0| 78 43 6f 70 79 28 20 61 | 70 2c 20 62 70 20 29 0a |xCopy( a|p, bp ).|
|000044b0| 58 09 63 6f 6d 70 6c 65 | 78 09 2a 61 70 2c 20 2a |X.comple|x.*ap, *|
|000044c0| 62 70 3b 09 09 2f 2a 20 | 6d 61 79 20 63 6f 69 6e |bp;../* |may coin|
|000044d0| 63 69 64 65 20 2a 2f 0a | 58 09 7b 0a 58 09 2a 61 |cide */.|X.{.X.*a|
|000044e0| 70 20 3d 20 2a 62 70 3b | 0a 58 0a 58 09 72 65 74 |p = *bp;|.X.X.ret|
|000044f0| 75 72 6e 20 61 70 3b 0a | 58 09 7d 0a 45 4e 44 5f |urn ap;.|X.}.END_|
|00004500| 4f 46 5f 63 78 63 6f 70 | 79 2e 63 0a 69 66 20 74 |OF_cxcop|y.c.if t|
|00004510| 65 73 74 20 32 36 34 20 | 2d 6e 65 20 60 77 63 20 |est 264 |-ne `wc |
|00004520| 2d 63 20 3c 63 78 63 6f | 70 79 2e 63 60 3b 20 74 |-c <cxco|py.c`; t|
|00004530| 68 65 6e 0a 20 20 20 20 | 65 63 68 6f 20 73 68 61 |hen.    |echo sha|
|00004540| 72 3a 20 5c 22 63 78 63 | 6f 70 79 2e 63 5c 22 20 |r: \"cxc|opy.c\" |
|00004550| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|00004560| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|00004570| 6e 64 20 6f 66 20 6f 76 | 65 72 77 72 69 74 69 6e |nd of ov|erwritin|
|00004580| 67 20 63 68 65 63 6b 0a | 66 69 0a 69 66 20 74 65 |g check.|fi.if te|
|00004590| 73 74 20 2d 66 20 63 78 | 64 69 76 2e 63 20 2d 61 |st -f cx|div.c -a|
|000045a0| 20 22 24 7b 31 7d 22 20 | 21 3d 20 22 2d 63 22 20 | "${1}" |!= "-c" |
|000045b0| 3b 20 74 68 65 6e 20 0a | 20 20 65 63 68 6f 20 73 |; then .|  echo s|
|000045c0| 68 61 72 3a 20 57 69 6c | 6c 20 6e 6f 74 20 6f 76 |har: Wil|l not ov|
|000045d0| 65 72 2d 77 72 69 74 65 | 20 65 78 69 73 74 69 6e |er-write| existin|
|000045e0| 67 20 66 69 6c 65 20 5c | 22 63 78 64 69 76 2e 63 |g file \|"cxdiv.c|
|000045f0| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|00004600| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|00004610| 63 78 64 69 76 2e 63 5c | 22 20 5c 28 38 32 30 20 |cxdiv.c\|" \(820 |
|00004620| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00004630| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 63 78 64 69 76 | "s/^X//|" >cxdiv|
|00004640| 2e 63 20 3c 3c 27 45 4e | 44 5f 4f 46 5f 63 78 64 |.c <<'EN|D_OF_cxd|
|00004650| 69 76 2e 63 27 0a 58 2f | 2a 0a 58 09 43 78 44 69 |iv.c'.X/|*.X.CxDi|
|00004660| 76 20 2d 2d 20 64 69 76 | 69 64 65 20 6f 6e 65 20 |v -- div|ide one |
|00004670| 63 6f 6d 70 6c 65 78 20 | 62 79 20 61 6e 6f 74 68 |complex |by anoth|
|00004680| 65 72 0a 58 0a 58 09 43 | 78 44 69 76 28 20 26 61 |er.X.X.C|xDiv( &a|
|00004690| 2c 20 26 62 20 29 09 64 | 69 76 69 64 65 73 20 20 |, &b ).d|ivides  |
|000046a0| 61 20 20 62 79 20 20 62 | 20 20 61 6e 64 20 72 65 |a  by  b|  and re|
|000046b0| 74 75 72 6e 73 20 20 26 | 61 3b 0a 58 09 09 09 7a |turns  &|a;.X...z|
|000046c0| 65 72 6f 20 64 69 76 69 | 73 6f 72 20 66 61 69 6c |ero divi|sor fail|
|000046d0| 73 0a 58 0a 58 09 6c 61 | 73 74 20 65 64 69 74 3a |s.X.X.la|st edit:|
|000046e0| 09 38 36 2f 30 31 2f 30 | 34 09 44 20 41 20 47 77 |.86/01/0|4.D A Gw|
|000046f0| 79 6e 0a 58 0a 58 09 53 | 43 43 53 20 49 44 3a 09 |yn.X.X.S|CCS ID:.|
|00004700| 40 28 23 29 63 78 64 69 | 76 2e 63 09 31 2e 31 20 |@(#)cxdi|v.c.1.1 |
|00004710| 28 6d 6f 64 69 66 69 65 | 64 20 66 6f 72 20 70 75 |(modifie|d for pu|
|00004720| 62 6c 69 63 20 76 65 72 | 73 69 6f 6e 29 0a 58 2a |blic ver|sion).X*|
|00004730| 2f 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 09 3c 63 |/.X.X#in|clude.<c|
|00004740| 6f 6d 70 6c 65 78 2e 68 | 3e 0a 58 0a 58 23 64 65 |omplex.h|>.X.X#de|
|00004750| 66 69 6e 65 20 41 62 73 | 28 20 78 20 29 09 28 28 |fine Abs|( x ).((|
|00004760| 78 29 20 3c 20 30 20 3f | 20 2d 28 78 29 20 3a 20 |x) < 0 ?| -(x) : |
|00004770| 28 78 29 29 0a 58 0a 58 | 63 6f 6d 70 6c 65 78 20 |(x)).X.X|complex |
|00004780| 2a 0a 58 43 78 44 69 76 | 28 20 61 70 2c 20 62 70 |*.XCxDiv|( ap, bp|
|00004790| 20 29 0a 58 09 72 65 67 | 69 73 74 65 72 20 63 6f | ).X.reg|ister co|
|000047a0| 6d 70 6c 65 78 09 2a 61 | 70 2c 20 2a 62 70 3b 09 |mplex.*a|p, *bp;.|
|000047b0| 2f 2a 20 6d 61 79 20 63 | 6f 69 6e 63 69 64 65 20 |/* may c|oincide |
|000047c0| 28 3f 29 20 2a 2f 0a 58 | 09 7b 0a 58 09 64 6f 75 |(?) */.X|.{.X.dou|
|000047d0| 62 6c 65 09 72 2c 20 73 | 3b 0a 58 09 64 6f 75 62 |ble.r, s|;.X.doub|
|000047e0| 6c 65 09 61 70 5f 5f 72 | 65 20 3d 20 61 70 2d 3e |le.ap__r|e = ap->|
|000047f0| 72 65 3b 0a 58 0a 58 09 | 2f 2a 20 4e 6f 74 65 3a |re;.X.X.|/* Note:|
|00004800| 20 63 6c 61 73 73 69 63 | 61 6c 20 66 6f 72 6d 75 | classic|al formu|
|00004810| 6c 61 20 6d 61 79 20 63 | 61 75 73 65 20 75 6e 6e |la may c|ause unn|
|00004820| 65 63 65 73 73 61 72 79 | 20 6f 76 65 72 66 6c 6f |ecessary| overflo|
|00004830| 77 20 2a 2f 0a 58 09 72 | 20 3d 20 62 70 2d 3e 72 |w */.X.r| = bp->r|
|00004840| 65 3b 0a 58 09 73 20 3d | 20 62 70 2d 3e 69 6d 3b |e;.X.s =| bp->im;|
|00004850| 0a 58 09 69 66 20 28 20 | 41 62 73 28 20 72 20 29 |.X.if ( |Abs( r )|
|00004860| 20 3e 3d 20 41 62 73 28 | 20 73 20 29 20 29 0a 58 | >= Abs(| s ) ).X|
|00004870| 09 09 7b 0a 58 09 09 72 | 20 3d 20 73 20 2f 20 72 |..{.X..r| = s / r|
|00004880| 3b 09 09 2f 2a 20 3c 3d | 20 31 20 2a 2f 0a 58 09 |;../* <=| 1 */.X.|
|00004890| 09 73 20 3d 20 62 70 2d | 3e 72 65 20 2b 20 72 20 |.s = bp-|>re + r |
|000048a0| 2a 20 73 3b 0a 58 09 09 | 61 70 2d 3e 72 65 20 3d |* s;.X..|ap->re =|
|000048b0| 20 28 61 70 2d 3e 72 65 | 20 2b 20 61 70 2d 3e 69 | (ap->re| + ap->i|
|000048c0| 6d 20 2a 20 72 29 20 2f | 20 73 3b 0a 58 09 09 61 |m * r) /| s;.X..a|
|000048d0| 70 2d 3e 69 6d 20 3d 20 | 28 61 70 2d 3e 69 6d 20 |p->im = |(ap->im |
|000048e0| 2d 20 61 70 5f 5f 72 65 | 20 2a 20 72 29 20 2f 20 |- ap__re| * r) / |
|000048f0| 73 3b 0a 58 09 09 7d 0a | 58 09 65 6c 73 65 20 2f |s;.X..}.|X.else /|
|00004900| 2a 20 41 62 73 28 20 73 | 20 29 20 3e 20 41 62 73 |* Abs( s| ) > Abs|
|00004910| 28 20 72 20 29 20 2a 2f | 0a 58 09 09 7b 0a 58 09 |( r ) */|.X..{.X.|
|00004920| 09 72 20 3d 20 72 20 2f | 20 73 3b 09 09 2f 2a 20 |.r = r /| s;../* |
|00004930| 3c 20 31 20 2a 2f 0a 58 | 09 09 73 20 3d 20 73 20 |< 1 */.X|..s = s |
|00004940| 2b 20 72 20 2a 20 62 70 | 2d 3e 72 65 3b 0a 58 09 |+ r * bp|->re;.X.|
|00004950| 09 61 70 2d 3e 72 65 20 | 3d 20 28 61 70 2d 3e 72 |.ap->re |= (ap->r|
|00004960| 65 20 2a 20 72 20 2b 20 | 61 70 2d 3e 69 6d 29 20 |e * r + |ap->im) |
|00004970| 2f 20 73 3b 0a 58 09 09 | 61 70 2d 3e 69 6d 20 3d |/ s;.X..|ap->im =|
|00004980| 20 28 61 70 2d 3e 69 6d | 20 2a 20 72 20 2d 20 61 | (ap->im| * r - a|
|00004990| 70 5f 5f 72 65 29 20 2f | 20 73 3b 0a 58 09 09 7d |p__re) /| s;.X..}|
|000049a0| 0a 58 0a 58 09 72 65 74 | 75 72 6e 20 61 70 3b 0a |.X.X.ret|urn ap;.|
|000049b0| 58 09 7d 0a 45 4e 44 5f | 4f 46 5f 63 78 64 69 76 |X.}.END_|OF_cxdiv|
|000049c0| 2e 63 0a 69 66 20 74 65 | 73 74 20 38 32 30 20 2d |.c.if te|st 820 -|
|000049d0| 6e 65 20 60 77 63 20 2d | 63 20 3c 63 78 64 69 76 |ne `wc -|c <cxdiv|
|000049e0| 2e 63 60 3b 20 74 68 65 | 6e 0a 20 20 20 20 65 63 |.c`; the|n.    ec|
|000049f0| 68 6f 20 73 68 61 72 3a | 20 5c 22 63 78 64 69 76 |ho shar:| \"cxdiv|
|00004a00| 2e 63 5c 22 20 75 6e 70 | 61 63 6b 65 64 20 77 69 |.c\" unp|acked wi|
|00004a10| 74 68 20 77 72 6f 6e 67 | 20 73 69 7a 65 21 0a 66 |th wrong| size!.f|
|00004a20| 69 0a 23 20 65 6e 64 20 | 6f 66 20 6f 76 65 72 77 |i.# end |of overw|
|00004a30| 72 69 74 69 6e 67 20 63 | 68 65 63 6b 0a 66 69 0a |riting c|heck.fi.|
|00004a40| 69 66 20 74 65 73 74 20 | 2d 66 20 63 78 6d 75 6c |if test |-f cxmul|
|00004a50| 2e 63 20 2d 61 20 22 24 | 7b 31 7d 22 20 21 3d 20 |.c -a "$|{1}" != |
|00004a60| 22 2d 63 22 20 3b 20 74 | 68 65 6e 20 0a 20 20 65 |"-c" ; t|hen .  e|
|00004a70| 63 68 6f 20 73 68 61 72 | 3a 20 57 69 6c 6c 20 6e |cho shar|: Will n|
|00004a80| 6f 74 20 6f 76 65 72 2d | 77 72 69 74 65 20 65 78 |ot over-|write ex|
|00004a90| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 63 78 |isting f|ile \"cx|
|00004aa0| 6d 75 6c 2e 63 5c 22 0a | 65 6c 73 65 0a 65 63 68 |mul.c\".|else.ech|
|00004ab0| 6f 20 73 68 61 72 3a 20 | 45 78 74 72 61 63 74 69 |o shar: |Extracti|
|00004ac0| 6e 67 20 5c 22 63 78 6d | 75 6c 2e 63 5c 22 20 5c |ng \"cxm|ul.c\" \|
|00004ad0| 28 34 32 34 20 63 68 61 | 72 61 63 74 65 72 73 5c |(424 cha|racters\|
|00004ae0| 29 0a 73 65 64 20 22 73 | 2f 5e 58 2f 2f 22 20 3e |).sed "s|/^X//" >|
|00004af0| 63 78 6d 75 6c 2e 63 20 | 3c 3c 27 45 4e 44 5f 4f |cxmul.c |<<'END_O|
|00004b00| 46 5f 63 78 6d 75 6c 2e | 63 27 0a 58 2f 2a 0a 58 |F_cxmul.|c'.X/*.X|
|00004b10| 09 43 78 4d 75 6c 20 2d | 2d 20 6d 75 6c 74 69 70 |.CxMul -|- multip|
|00004b20| 6c 79 20 6f 6e 65 20 63 | 6f 6d 70 6c 65 78 20 62 |ly one c|omplex b|
|00004b30| 79 20 61 6e 6f 74 68 65 | 72 0a 58 0a 58 09 43 78 |y anothe|r.X.X.Cx|
|00004b40| 4d 75 6c 28 20 26 61 2c | 20 26 62 20 29 09 6d 75 |Mul( &a,| &b ).mu|
|00004b50| 6c 74 69 70 6c 69 65 73 | 20 20 61 20 20 62 79 20 |ltiplies|  a  by |
|00004b60| 20 62 20 20 61 6e 64 20 | 72 65 74 75 72 6e 73 20 | b  and |returns |
|00004b70| 20 26 61 0a 58 0a 58 09 | 6c 61 73 74 20 65 64 69 | &a.X.X.|last edi|
|00004b80| 74 3a 09 38 36 2f 30 31 | 2f 30 34 09 44 20 41 20 |t:.86/01|/04.D A |
|00004b90| 47 77 79 6e 0a 58 0a 58 | 09 53 43 43 53 20 49 44 |Gwyn.X.X|.SCCS ID|
|00004ba0| 3a 09 40 28 23 29 63 78 | 6d 75 6c 2e 63 09 31 2e |:.@(#)cx|mul.c.1.|
|00004bb0| 31 0a 58 2a 2f 0a 58 0a | 58 23 69 6e 63 6c 75 64 |1.X*/.X.|X#includ|
|00004bc0| 65 09 3c 63 6f 6d 70 6c | 65 78 2e 68 3e 0a 58 0a |e.<compl|ex.h>.X.|
|00004bd0| 58 63 6f 6d 70 6c 65 78 | 20 2a 0a 58 43 78 4d 75 |Xcomplex| *.XCxMu|
|00004be0| 6c 28 20 61 70 2c 20 62 | 70 20 29 0a 58 09 72 65 |l( ap, b|p ).X.re|
|00004bf0| 67 69 73 74 65 72 20 63 | 6f 6d 70 6c 65 78 09 2a |gister c|omplex.*|
|00004c00| 61 70 2c 20 2a 62 70 3b | 09 2f 2a 20 28 6d 61 79 |ap, *bp;|./* (may|
|00004c10| 20 63 6f 69 6e 63 69 64 | 65 29 20 2a 2f 0a 58 09 | coincid|e) */.X.|
|00004c20| 7b 0a 58 09 64 6f 75 62 | 6c 65 09 09 09 61 70 5f |{.X.doub|le...ap_|
|00004c30| 5f 72 65 20 3d 20 61 70 | 2d 3e 72 65 3b 0a 58 09 |_re = ap|->re;.X.|
|00004c40| 64 6f 75 62 6c 65 09 09 | 09 62 70 5f 5f 72 65 20 |double..|.bp__re |
|00004c50| 3d 20 62 70 2d 3e 72 65 | 3b 0a 58 0a 58 09 61 70 |= bp->re|;.X.X.ap|
|00004c60| 2d 3e 72 65 20 3d 20 61 | 70 5f 5f 72 65 20 2a 20 |->re = a|p__re * |
|00004c70| 62 70 5f 5f 72 65 20 2d | 20 61 70 2d 3e 69 6d 20 |bp__re -| ap->im |
|00004c80| 2a 20 62 70 2d 3e 69 6d | 3b 0a 58 09 61 70 2d 3e |* bp->im|;.X.ap->|
|00004c90| 69 6d 20 3d 20 61 70 5f | 5f 72 65 20 2a 20 62 70 |im = ap_|_re * bp|
|00004ca0| 2d 3e 69 6d 20 2b 20 61 | 70 2d 3e 69 6d 20 2a 20 |->im + a|p->im * |
|00004cb0| 62 70 5f 5f 72 65 3b 0a | 58 0a 58 09 72 65 74 75 |bp__re;.|X.X.retu|
|00004cc0| 72 6e 20 61 70 3b 0a 58 | 09 7d 0a 45 4e 44 5f 4f |rn ap;.X|.}.END_O|
|00004cd0| 46 5f 63 78 6d 75 6c 2e | 63 0a 69 66 20 74 65 73 |F_cxmul.|c.if tes|
|00004ce0| 74 20 34 32 34 20 2d 6e | 65 20 60 77 63 20 2d 63 |t 424 -n|e `wc -c|
|00004cf0| 20 3c 63 78 6d 75 6c 2e | 63 60 3b 20 74 68 65 6e | <cxmul.|c`; then|
|00004d00| 0a 20 20 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |.    ech|o shar: |
|00004d10| 5c 22 63 78 6d 75 6c 2e | 63 5c 22 20 75 6e 70 61 |\"cxmul.|c\" unpa|
|00004d20| 63 6b 65 64 20 77 69 74 | 68 20 77 72 6f 6e 67 20 |cked wit|h wrong |
|00004d30| 73 69 7a 65 21 0a 66 69 | 0a 23 20 65 6e 64 20 6f |size!.fi|.# end o|
|00004d40| 66 20 6f 76 65 72 77 72 | 69 74 69 6e 67 20 63 68 |f overwr|iting ch|
|00004d50| 65 63 6b 0a 66 69 0a 69 | 66 20 74 65 73 74 20 2d |eck.fi.i|f test -|
|00004d60| 66 20 63 78 70 68 61 73 | 2e 63 20 2d 61 20 22 24 |f cxphas|.c -a "$|
|00004d70| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|00004d80| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen .  e|cho shar|
|00004d90| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 6f 76 65 72 2d |: Will n|ot over-|
|00004da0| 77 72 69 74 65 20 65 78 | 69 73 74 69 6e 67 20 66 |write ex|isting f|
|00004db0| 69 6c 65 20 5c 22 63 78 | 70 68 61 73 2e 63 5c 22 |ile \"cx|phas.c\"|
|00004dc0| 0a 65 6c 73 65 0a 65 63 | 68 6f 20 73 68 61 72 3a |.else.ec|ho shar:|
|00004dd0| 20 45 78 74 72 61 63 74 | 69 6e 67 20 5c 22 63 78 | Extract|ing \"cx|
|00004de0| 70 68 61 73 2e 63 5c 22 | 20 5c 28 34 30 36 20 63 |phas.c\"| \(406 c|
|00004df0| 68 61 72 61 63 74 65 72 | 73 5c 29 0a 73 65 64 20 |haracter|s\).sed |
|00004e00| 22 73 2f 5e 58 2f 2f 22 | 20 3e 63 78 70 68 61 73 |"s/^X//"| >cxphas|
|00004e10| 2e 63 20 3c 3c 27 45 4e | 44 5f 4f 46 5f 63 78 70 |.c <<'EN|D_OF_cxp|
|00004e20| 68 61 73 2e 63 27 0a 58 | 2f 2a 0a 58 09 43 78 50 |has.c'.X|/*.X.CxP|
|00004e30| 68 61 73 20 2d 2d 20 70 | 68 61 73 65 20 28 61 6e |has -- p|hase (an|
|00004e40| 67 6c 65 2c 20 61 72 67 | 75 6d 65 6e 74 29 20 6f |gle, arg|ument) o|
|00004e50| 66 20 61 20 63 6f 6d 70 | 6c 65 78 0a 58 0a 58 09 |f a comp|lex.X.X.|
|00004e60| 43 78 50 68 61 73 28 20 | 26 63 20 29 09 72 65 74 |CxPhas( |&c ).ret|
|00004e70| 75 72 6e 73 20 20 61 72 | 67 28 63 29 20 20 69 6e |urns  ar|g(c)  in|
|00004e80| 20 72 61 64 69 61 6e 73 | 20 28 2d 50 69 2c 50 69 | radians| (-Pi,Pi|
|00004e90| 5d 0a 58 0a 58 09 6c 61 | 73 74 20 65 64 69 74 3a |].X.X.la|st edit:|
|00004ea0| 09 38 36 2f 30 31 2f 30 | 34 09 44 20 41 20 47 77 |.86/01/0|4.D A Gw|
|00004eb0| 79 6e 0a 58 0a 58 09 53 | 43 43 53 20 49 44 3a 09 |yn.X.X.S|CCS ID:.|
|00004ec0| 40 28 23 29 63 78 70 68 | 61 73 2e 63 09 31 2e 31 |@(#)cxph|as.c.1.1|
|00004ed0| 20 28 6d 6f 64 69 66 69 | 65 64 20 66 6f 72 20 70 | (modifi|ed for p|
|00004ee0| 75 62 6c 69 63 20 76 65 | 72 73 69 6f 6e 29 0a 58 |ublic ve|rsion).X|
|00004ef0| 2a 2f 0a 58 0a 58 23 69 | 6e 63 6c 75 64 65 09 3c |*/.X.X#i|nclude.<|
|00004f00| 6d 61 74 68 2e 68 3e 0a | 58 0a 58 23 69 6e 63 6c |math.h>.|X.X#incl|
|00004f10| 75 64 65 09 3c 63 6f 6d | 70 6c 65 78 2e 68 3e 0a |ude.<com|plex.h>.|
|00004f20| 58 0a 58 64 6f 75 62 6c | 65 0a 58 43 78 50 68 61 |X.Xdoubl|e.XCxPha|
|00004f30| 73 28 20 63 70 20 29 0a | 58 09 72 65 67 69 73 74 |s( cp ).|X.regist|
|00004f40| 65 72 20 63 6f 6d 70 6c | 65 78 09 2a 63 70 3b 0a |er compl|ex.*cp;.|
|00004f50| 58 09 7b 0a 58 09 69 66 | 20 28 20 63 70 2d 3e 72 |X.{.X.if| ( cp->r|
|00004f60| 65 20 3d 3d 20 30 2e 30 | 20 26 26 20 63 70 2d 3e |e == 0.0| && cp->|
|00004f70| 69 6d 20 3d 3d 20 30 2e | 30 20 29 0a 58 09 09 72 |im == 0.|0 ).X..r|
|00004f80| 65 74 75 72 6e 20 30 2e | 30 3b 09 09 2f 2a 20 63 |eturn 0.|0;../* c|
|00004f90| 61 6e 27 74 20 74 72 75 | 73 74 20 61 74 61 6e 32 |an't tru|st atan2|
|00004fa0| 28 29 20 2a 2f 0a 58 09 | 65 6c 73 65 0a 58 09 09 |() */.X.|else.X..|
|00004fb0| 72 65 74 75 72 6e 20 61 | 74 61 6e 32 28 20 63 70 |return a|tan2( cp|
|00004fc0| 2d 3e 69 6d 2c 20 63 70 | 2d 3e 72 65 20 29 3b 0a |->im, cp|->re );.|
|00004fd0| 58 09 7d 0a 45 4e 44 5f | 4f 46 5f 63 78 70 68 61 |X.}.END_|OF_cxpha|
|00004fe0| 73 2e 63 0a 69 66 20 74 | 65 73 74 20 34 30 36 20 |s.c.if t|est 406 |
|00004ff0| 2d 6e 65 20 60 77 63 20 | 2d 63 20 3c 63 78 70 68 |-ne `wc |-c <cxph|
|00005000| 61 73 2e 63 60 3b 20 74 | 68 65 6e 0a 20 20 20 20 |as.c`; t|hen.    |
|00005010| 65 63 68 6f 20 73 68 61 | 72 3a 20 5c 22 63 78 70 |echo sha|r: \"cxp|
|00005020| 68 61 73 2e 63 5c 22 20 | 75 6e 70 61 63 6b 65 64 |has.c\" |unpacked|
|00005030| 20 77 69 74 68 20 77 72 | 6f 6e 67 20 73 69 7a 65 | with wr|ong size|
|00005040| 21 0a 66 69 0a 23 20 65 | 6e 64 20 6f 66 20 6f 76 |!.fi.# e|nd of ov|
|00005050| 65 72 77 72 69 74 69 6e | 67 20 63 68 65 63 6b 0a |erwritin|g check.|
|00005060| 66 69 0a 69 66 20 74 65 | 73 74 20 2d 66 20 63 78 |fi.if te|st -f cx|
|00005070| 70 68 73 72 2e 63 20 2d | 61 20 22 24 7b 31 7d 22 |phsr.c -|a "${1}"|
|00005080| 20 21 3d 20 22 2d 63 22 | 20 3b 20 74 68 65 6e 20 | != "-c"| ; then |
|00005090| 0a 20 20 65 63 68 6f 20 | 73 68 61 72 3a 20 57 69 |.  echo |shar: Wi|
|000050a0| 6c 6c 20 6e 6f 74 20 6f | 76 65 72 2d 77 72 69 74 |ll not o|ver-writ|
|000050b0| 65 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |e existi|ng file |
|000050c0| 5c 22 63 78 70 68 73 72 | 2e 63 5c 22 0a 65 6c 73 |\"cxphsr|.c\".els|
|000050d0| 65 0a 65 63 68 6f 20 73 | 68 61 72 3a 20 45 78 74 |e.echo s|har: Ext|
|000050e0| 72 61 63 74 69 6e 67 20 | 5c 22 63 78 70 68 73 72 |racting |\"cxphsr|
|000050f0| 2e 63 5c 22 20 5c 28 33 | 39 35 20 63 68 61 72 61 |.c\" \(3|95 chara|
|00005100| 63 74 65 72 73 5c 29 0a | 73 65 64 20 22 73 2f 5e |cters\).|sed "s/^|
|00005110| 58 2f 2f 22 20 3e 63 78 | 70 68 73 72 2e 63 20 3c |X//" >cx|phsr.c <|
|00005120| 3c 27 45 4e 44 5f 4f 46 | 5f 63 78 70 68 73 72 2e |<'END_OF|_cxphsr.|
|00005130| 63 27 0a 58 2f 2a 0a 58 | 09 43 78 50 68 73 72 20 |c'.X/*.X|.CxPhsr |
|00005140| 2d 2d 20 63 6f 6e 73 74 | 72 75 63 74 20 61 20 63 |-- const|ruct a c|
|00005150| 6f 6d 70 6c 65 78 20 22 | 70 68 61 73 6f 72 22 20 |omplex "|phasor" |
|00005160| 66 72 6f 6d 20 61 6d 70 | 6c 69 74 75 64 65 20 61 |from amp|litude a|
|00005170| 6e 64 20 70 68 61 73 65 | 0a 58 0a 58 09 43 78 50 |nd phase|.X.X.CxP|
|00005180| 68 73 72 28 20 26 63 2c | 20 61 6d 70 2c 20 70 68 |hsr( &c,| amp, ph|
|00005190| 73 20 29 09 6d 61 6b 65 | 73 20 20 63 20 3d 20 61 |s ).make|s  c = a|
|000051a0| 6d 70 20 65 78 70 28 69 | 20 70 68 73 29 0a 58 09 |mp exp(i| phs).X.|
|000051b0| 09 09 09 09 61 6e 64 20 | 72 65 74 75 72 6e 73 20 |....and |returns |
|000051c0| 20 26 63 0a 58 0a 58 09 | 6c 61 73 74 20 65 64 69 | &c.X.X.|last edi|
|000051d0| 74 3a 09 38 36 2f 30 31 | 2f 30 34 09 44 20 41 20 |t:.86/01|/04.D A |
|000051e0| 47 77 79 6e 0a 58 0a 58 | 09 53 43 43 53 20 49 44 |Gwyn.X.X|.SCCS ID|
|000051f0| 3a 09 40 28 23 29 63 78 | 70 68 73 72 2e 63 09 31 |:.@(#)cx|phsr.c.1|
|00005200| 2e 31 0a 58 2a 2f 0a 58 | 0a 58 23 69 6e 63 6c 75 |.1.X*/.X|.X#inclu|
|00005210| 64 65 09 3c 6d 61 74 68 | 2e 68 3e 0a 58 0a 58 23 |de.<math|.h>.X.X#|
|00005220| 69 6e 63 6c 75 64 65 09 | 3c 63 6f 6d 70 6c 65 78 |include.|<complex|
|00005230| 2e 68 3e 0a 58 0a 58 63 | 6f 6d 70 6c 65 78 20 2a |.h>.X.Xc|omplex *|
|00005240| 0a 58 43 78 50 68 73 72 | 28 20 63 70 2c 20 61 6d |.XCxPhsr|( cp, am|
|00005250| 70 2c 20 70 68 73 20 29 | 0a 58 09 72 65 67 69 73 |p, phs )|.X.regis|
|00005260| 74 65 72 20 63 6f 6d 70 | 6c 65 78 09 2a 63 70 3b |ter comp|lex.*cp;|
|00005270| 0a 58 09 64 6f 75 62 6c | 65 09 09 09 61 6d 70 2c |.X.doubl|e...amp,|
|00005280| 20 70 68 73 3b 0a 58 09 | 7b 0a 58 09 63 70 2d 3e | phs;.X.|{.X.cp->|
|00005290| 72 65 20 3d 20 61 6d 70 | 20 2a 20 63 6f 73 28 20 |re = amp| * cos( |
|000052a0| 70 68 73 20 29 3b 0a 58 | 09 63 70 2d 3e 69 6d 20 |phs );.X|.cp->im |
|000052b0| 3d 20 61 6d 70 20 2a 20 | 73 69 6e 28 20 70 68 73 |= amp * |sin( phs|
|000052c0| 20 29 3b 0a 58 0a 58 09 | 72 65 74 75 72 6e 20 63 | );.X.X.|return c|
|000052d0| 70 3b 0a 58 09 7d 0a 45 | 4e 44 5f 4f 46 5f 63 78 |p;.X.}.E|ND_OF_cx|
|000052e0| 70 68 73 72 2e 63 0a 69 | 66 20 74 65 73 74 20 33 |phsr.c.i|f test 3|
|000052f0| 39 35 20 2d 6e 65 20 60 | 77 63 20 2d 63 20 3c 63 |95 -ne `|wc -c <c|
|00005300| 78 70 68 73 72 2e 63 60 | 3b 20 74 68 65 6e 0a 20 |xphsr.c`|; then. |
|00005310| 20 20 20 65 63 68 6f 20 | 73 68 61 72 3a 20 5c 22 |   echo |shar: \"|
|00005320| 63 78 70 68 73 72 2e 63 | 5c 22 20 75 6e 70 61 63 |cxphsr.c|\" unpac|
|00005330| 6b 65 64 20 77 69 74 68 | 20 77 72 6f 6e 67 20 73 |ked with| wrong s|
|00005340| 69 7a 65 21 0a 66 69 0a | 23 20 65 6e 64 20 6f 66 |ize!.fi.|# end of|
|00005350| 20 6f 76 65 72 77 72 69 | 74 69 6e 67 20 63 68 65 | overwri|ting che|
|00005360| 63 6b 0a 66 69 0a 69 66 | 20 74 65 73 74 20 2d 66 |ck.fi.if| test -f|
|00005370| 20 63 78 73 63 61 6c 2e | 63 20 2d 61 20 22 24 7b | cxscal.|c -a "${|
|00005380| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|00005390| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en .  ec|ho shar:|
|000053a0| 20 57 69 6c 6c 20 6e 6f | 74 20 6f 76 65 72 2d 77 | Will no|t over-w|
|000053b0| 72 69 74 65 20 65 78 69 | 73 74 69 6e 67 20 66 69 |rite exi|sting fi|
|000053c0| 6c 65 20 5c 22 63 78 73 | 63 61 6c 2e 63 5c 22 0a |le \"cxs|cal.c\".|
|000053d0| 65 6c 73 65 0a 65 63 68 | 6f 20 73 68 61 72 3a 20 |else.ech|o shar: |
|000053e0| 45 78 74 72 61 63 74 69 | 6e 67 20 5c 22 63 78 73 |Extracti|ng \"cxs|
|000053f0| 63 61 6c 2e 63 5c 22 20 | 5c 28 32 39 31 20 63 68 |cal.c\" |\(291 ch|
|00005400| 61 72 61 63 74 65 72 73 | 5c 29 0a 73 65 64 20 22 |aracters|\).sed "|
|00005410| 73 2f 5e 58 2f 2f 22 20 | 3e 63 78 73 63 61 6c 2e |s/^X//" |>cxscal.|
|00005420| 63 20 3c 3c 27 45 4e 44 | 5f 4f 46 5f 63 78 73 63 |c <<'END|_OF_cxsc|
|00005430| 61 6c 2e 63 27 0a 58 2f | 2a 0a 58 09 43 78 53 63 |al.c'.X/|*.X.CxSc|
|00005440| 61 6c 20 2d 2d 20 6d 75 | 6c 74 69 70 6c 79 20 61 |al -- mu|ltiply a|
|00005450| 20 63 6f 6d 70 6c 65 78 | 20 62 79 20 61 20 73 63 | complex| by a sc|
|00005460| 61 6c 61 72 0a 58 0a 58 | 09 43 78 53 63 61 6c 28 |alar.X.X|.CxScal(|
|00005470| 20 26 63 2c 20 73 20 29 | 09 73 63 61 6c 65 73 20 | &c, s )|.scales |
|00005480| 20 63 20 20 62 79 20 20 | 73 20 20 61 6e 64 20 72 | c  by  |s  and r|
|00005490| 65 74 75 72 6e 73 20 20 | 26 63 0a 58 0a 58 09 6c |eturns  |&c.X.X.l|
|000054a0| 61 73 74 20 65 64 69 74 | 3a 09 38 36 2f 30 31 2f |ast edit|:.86/01/|
|000054b0| 30 34 09 44 20 41 20 47 | 77 79 6e 0a 58 0a 58 09 |04.D A G|wyn.X.X.|
|000054c0| 53 43 43 53 20 49 44 3a | 09 40 28 23 29 63 78 73 |SCCS ID:|.@(#)cxs|
|000054d0| 63 61 6c 2e 63 09 31 2e | 31 0a 58 2a 2f 0a 58 0a |cal.c.1.|1.X*/.X.|
|000054e0| 58 23 69 6e 63 6c 75 64 | 65 09 3c 63 6f 6d 70 6c |X#includ|e.<compl|
|000054f0| 65 78 2e 68 3e 0a 58 0a | 58 63 6f 6d 70 6c 65 78 |ex.h>.X.|Xcomplex|
|00005500| 20 2a 0a 58 43 78 53 63 | 61 6c 28 20 63 70 2c 20 | *.XCxSc|al( cp, |
|00005510| 73 20 29 0a 58 09 72 65 | 67 69 73 74 65 72 20 63 |s ).X.re|gister c|
|00005520| 6f 6d 70 6c 65 78 09 2a | 63 70 3b 0a 58 09 64 6f |omplex.*|cp;.X.do|
|00005530| 75 62 6c 65 09 09 09 73 | 3b 0a 58 09 7b 0a 58 09 |uble...s|;.X.{.X.|
|00005540| 63 70 2d 3e 72 65 20 2a | 3d 20 73 3b 0a 58 09 63 |cp->re *|= s;.X.c|
|00005550| 70 2d 3e 69 6d 20 2a 3d | 20 73 3b 0a 58 0a 58 09 |p->im *=| s;.X.X.|
|00005560| 72 65 74 75 72 6e 20 63 | 70 3b 0a 58 09 7d 0a 45 |return c|p;.X.}.E|
|00005570| 4e 44 5f 4f 46 5f 63 78 | 73 63 61 6c 2e 63 0a 69 |ND_OF_cx|scal.c.i|
|00005580| 66 20 74 65 73 74 20 32 | 39 31 20 2d 6e 65 20 60 |f test 2|91 -ne `|
|00005590| 77 63 20 2d 63 20 3c 63 | 78 73 63 61 6c 2e 63 60 |wc -c <c|xscal.c`|
|000055a0| 3b 20 74 68 65 6e 0a 20 | 20 20 20 65 63 68 6f 20 |; then. |   echo |
|000055b0| 73 68 61 72 3a 20 5c 22 | 63 78 73 63 61 6c 2e 63 |shar: \"|cxscal.c|
|000055c0| 5c 22 20 75 6e 70 61 63 | 6b 65 64 20 77 69 74 68 |\" unpac|ked with|
|000055d0| 20 77 72 6f 6e 67 20 73 | 69 7a 65 21 0a 66 69 0a | wrong s|ize!.fi.|
|000055e0| 23 20 65 6e 64 20 6f 66 | 20 6f 76 65 72 77 72 69 |# end of| overwri|
|000055f0| 74 69 6e 67 20 63 68 65 | 63 6b 0a 66 69 0a 69 66 |ting che|ck.fi.if|
|00005600| 20 74 65 73 74 20 2d 66 | 20 63 78 73 71 72 74 2e | test -f| cxsqrt.|
|00005610| 63 20 2d 61 20 22 24 7b | 31 7d 22 20 21 3d 20 22 |c -a "${|1}" != "|
|00005620| 2d 63 22 20 3b 20 74 68 | 65 6e 20 0a 20 20 65 63 |-c" ; th|en .  ec|
|00005630| 68 6f 20 73 68 61 72 3a | 20 57 69 6c 6c 20 6e 6f |ho shar:| Will no|
|00005640| 74 20 6f 76 65 72 2d 77 | 72 69 74 65 20 65 78 69 |t over-w|rite exi|
|00005650| 73 74 69 6e 67 20 66 69 | 6c 65 20 5c 22 63 78 73 |sting fi|le \"cxs|
|00005660| 71 72 74 2e 63 5c 22 0a | 65 6c 73 65 0a 65 63 68 |qrt.c\".|else.ech|
|00005670| 6f 20 73 68 61 72 3a 20 | 45 78 74 72 61 63 74 69 |o shar: |Extracti|
|00005680| 6e 67 20 5c 22 63 78 73 | 71 72 74 2e 63 5c 22 20 |ng \"cxs|qrt.c\" |
|00005690| 5c 28 31 33 33 39 20 63 | 68 61 72 61 63 74 65 72 |\(1339 c|haracter|
|000056a0| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|000056b0| 20 3e 63 78 73 71 72 74 | 2e 63 20 3c 3c 27 45 4e | >cxsqrt|.c <<'EN|
|000056c0| 44 5f 4f 46 5f 63 78 73 | 71 72 74 2e 63 27 0a 58 |D_OF_cxs|qrt.c'.X|
|000056d0| 2f 2a 0a 58 09 43 78 53 | 71 72 74 20 2d 2d 20 63 |/*.X.CxS|qrt -- c|
|000056e0| 6f 6d 70 75 74 65 20 73 | 71 75 61 72 65 20 72 6f |ompute s|quare ro|
|000056f0| 6f 74 20 6f 66 20 63 6f | 6d 70 6c 65 78 20 6e 75 |ot of co|mplex nu|
|00005700| 6d 62 65 72 0a 58 0a 58 | 09 43 78 53 71 72 74 28 |mber.X.X|.CxSqrt(|
|00005710| 20 26 63 20 29 09 72 65 | 70 6c 61 63 65 73 20 20 | &c ).re|places  |
|00005720| 63 20 20 62 79 20 20 73 | 71 72 74 28 63 29 20 20 |c  by  s|qrt(c)  |
|00005730| 61 6e 64 20 72 65 74 75 | 72 6e 73 20 20 26 63 0a |and retu|rns  &c.|
|00005740| 58 0a 58 09 4e 6f 74 65 | 3a 09 54 68 69 73 20 69 |X.X.Note|:.This i|
|00005750| 73 20 61 20 64 6f 75 62 | 6c 65 2d 76 61 6c 75 65 |s a doub|le-value|
|00005760| 64 20 66 75 6e 63 74 69 | 6f 6e 3b 20 74 68 65 20 |d functi|on; the |
|00005770| 72 65 73 75 6c 74 20 6f | 66 0a 58 09 09 43 78 53 |result o|f.X..CxS|
|00005780| 71 72 74 28 29 20 61 6c | 77 61 79 73 20 68 61 73 |qrt() al|ways has|
|00005790| 20 6e 6f 6e 6e 65 67 61 | 74 69 76 65 20 69 6d 61 | nonnega|tive ima|
|000057a0| 67 69 6e 61 72 79 20 70 | 61 72 74 2e 0a 58 0a 58 |ginary p|art..X.X|
|000057b0| 09 69 6e 73 70 69 72 65 | 64 20 62 79 20 4a 65 66 |.inspire|d by Jef|
|000057c0| 66 20 48 61 6e 65 73 27 | 20 76 65 72 73 69 6f 6e |f Hanes'| version|
|000057d0| 0a 58 0a 58 09 6c 61 73 | 74 20 65 64 69 74 3a 09 |.X.X.las|t edit:.|
|000057e0| 38 36 2f 30 31 2f 30 34 | 09 44 20 41 20 47 77 79 |86/01/04|.D A Gwy|
|000057f0| 6e 0a 58 0a 58 09 53 43 | 43 53 20 49 44 3a 09 40 |n.X.X.SC|CS ID:.@|
|00005800| 28 23 29 63 78 73 71 72 | 74 2e 63 09 31 2e 31 20 |(#)cxsqr|t.c.1.1 |
|00005810| 28 6d 6f 64 69 66 69 65 | 64 20 66 6f 72 20 70 75 |(modifie|d for pu|
|00005820| 62 6c 69 63 20 76 65 72 | 73 69 6f 6e 29 0a 58 2a |blic ver|sion).X*|
|00005830| 2f 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 09 3c 6d |/.X.X#in|clude.<m|
|00005840| 61 74 68 2e 68 3e 0a 58 | 0a 58 23 69 6e 63 6c 75 |ath.h>.X|.X#inclu|
|00005850| 64 65 09 3c 63 6f 6d 70 | 6c 65 78 2e 68 3e 0a 58 |de.<comp|lex.h>.X|
|00005860| 0a 58 23 64 65 66 69 6e | 65 09 53 67 6e 28 20 78 |.X#defin|e.Sgn( x|
|00005870| 20 29 09 28 28 78 29 20 | 3d 3d 20 30 20 3f 20 30 | ).((x) |== 0 ? 0|
|00005880| 20 3a 20 28 78 29 20 3e | 20 30 20 3f 20 31 20 3a | : (x) >| 0 ? 1 :|
|00005890| 20 2d 31 29 0a 58 0a 58 | 63 6f 6d 70 6c 65 78 20 | -1).X.X|complex |
|000058a0| 2a 0a 58 43 78 53 71 72 | 74 28 20 63 70 20 29 0a |*.XCxSqr|t( cp ).|
|000058b0| 58 09 72 65 67 69 73 74 | 65 72 20 63 6f 6d 70 6c |X.regist|er compl|
|000058c0| 65 78 09 2a 63 70 3b 0a | 58 09 7b 0a 58 09 2f 2a |ex.*cp;.|X.{.X./*|
|000058d0| 20 72 65 63 6f 72 64 20 | 73 69 67 6e 73 20 6f 66 | record |signs of|
|000058e0| 20 6f 72 69 67 69 6e 61 | 6c 20 72 65 61 6c 20 26 | origina|l real &|
|000058f0| 20 69 6d 61 67 69 6e 61 | 72 79 20 70 61 72 74 73 | imagina|ry parts|
|00005900| 20 2a 2f 0a 58 09 69 6e | 74 09 09 09 72 65 5f 73 | */.X.in|t...re_s|
|00005910| 69 67 6e 20 3d 20 53 67 | 6e 28 20 63 70 2d 3e 72 |ign = Sg|n( cp->r|
|00005920| 65 20 29 3b 0a 58 09 69 | 6e 74 09 09 09 69 6d 5f |e );.X.i|nt...im_|
|00005930| 73 69 67 6e 20 3d 20 53 | 67 6e 28 20 63 70 2d 3e |sign = S|gn( cp->|
|00005940| 69 6d 20 29 3b 0a 58 0a | 58 09 2f 2a 20 73 70 65 |im );.X.|X./* spe|
|00005950| 63 69 61 6c 20 63 61 73 | 65 73 20 61 72 65 20 6e |cial cas|es are n|
|00005960| 6f 74 20 6e 65 63 65 73 | 73 61 72 79 3b 20 74 68 |ot neces|sary; th|
|00005970| 65 79 20 61 72 65 20 68 | 65 72 65 20 66 6f 72 20 |ey are h|ere for |
|00005980| 73 70 65 65 64 20 2a 2f | 0a 58 0a 58 09 69 66 20 |speed */|.X.X.if |
|00005990| 28 20 72 65 5f 73 69 67 | 6e 20 3d 3d 20 30 20 29 |( re_sig|n == 0 )|
|000059a0| 0a 58 09 09 69 66 20 28 | 20 69 6d 5f 73 69 67 6e |.X..if (| im_sign|
|000059b0| 20 3d 3d 20 30 20 29 0a | 58 09 09 09 3b 09 09 2f | == 0 ).|X...;../|
|000059c0| 2a 20 28 30 2c 30 29 20 | 61 6c 72 65 61 64 79 20 |* (0,0) |already |
|000059d0| 74 68 65 72 65 20 2a 2f | 0a 58 09 09 65 6c 73 65 |there */|.X..else|
|000059e0| 20 69 66 20 28 20 69 6d | 5f 73 69 67 6e 20 3e 20 | if ( im|_sign > |
|000059f0| 30 20 29 0a 58 09 09 09 | 63 70 2d 3e 72 65 20 3d |0 ).X...|cp->re =|
|00005a00| 20 63 70 2d 3e 69 6d 20 | 3d 20 73 71 72 74 28 20 | cp->im |= sqrt( |
|00005a10| 63 70 2d 3e 69 6d 20 2f | 20 32 2e 30 20 29 3b 0a |cp->im /| 2.0 );.|
|00005a20| 58 09 09 65 6c 73 65 09 | 09 09 2f 2a 20 69 6d 5f |X..else.|../* im_|
|00005a30| 73 69 67 6e 20 3c 20 30 | 20 2a 2f 0a 58 09 09 09 |sign < 0| */.X...|
|00005a40| 63 70 2d 3e 72 65 20 3d | 20 2d 28 63 70 2d 3e 69 |cp->re =| -(cp->i|
|00005a50| 6d 20 3d 20 73 71 72 74 | 28 20 2d 63 70 2d 3e 69 |m = sqrt|( -cp->i|
|00005a60| 6d 20 2f 20 32 2e 30 20 | 29 29 3b 0a 58 09 65 6c |m / 2.0 |));.X.el|
|00005a70| 73 65 20 69 66 20 28 20 | 69 6d 5f 73 69 67 6e 20 |se if ( |im_sign |
|00005a80| 3d 3d 20 30 20 29 0a 58 | 09 09 69 66 20 28 20 72 |== 0 ).X|..if ( r|
|00005a90| 65 5f 73 69 67 6e 20 3e | 20 30 20 29 0a 58 09 09 |e_sign >| 0 ).X..|
|00005aa0| 09 63 70 2d 3e 72 65 20 | 3d 20 73 71 72 74 28 20 |.cp->re |= sqrt( |
|00005ab0| 63 70 2d 3e 72 65 20 29 | 3b 0a 58 2f 2a 09 09 09 |cp->re )|;.X/*...|
|00005ac0| 63 70 2d 3e 69 6d 20 3d | 20 30 2e 30 3b 09 2f 2a |cp->im =| 0.0;./*|
|00005ad0| 20 30 20 61 6c 72 65 61 | 64 79 20 74 68 65 72 65 | 0 alrea|dy there|
|00005ae0| 20 2a 2f 0a 58 09 09 65 | 6c 73 65 09 7b 09 09 2f | */.X..e|lse.{../|
|00005af0| 2a 20 72 65 5f 73 69 67 | 6e 20 3c 20 30 20 2a 2f |* re_sig|n < 0 */|
|00005b00| 0a 58 09 09 09 63 70 2d | 3e 69 6d 20 3d 20 73 71 |.X...cp-|>im = sq|
|00005b10| 72 74 28 20 2d 63 70 2d | 3e 72 65 20 29 3b 0a 58 |rt( -cp-|>re );.X|
|00005b20| 09 09 09 63 70 2d 3e 72 | 65 20 3d 20 30 2e 30 3b |...cp->r|e = 0.0;|
|00005b30| 0a 58 09 09 09 7d 0a 58 | 09 65 6c 73 65 09 7b 09 |.X...}.X|.else.{.|
|00005b40| 09 09 2f 2a 20 6e 6f 20 | 73 68 6f 72 74 63 75 74 |../* no |shortcut|
|00005b50| 73 20 2a 2f 0a 58 09 09 | 64 6f 75 62 6c 65 09 61 |s */.X..|double.a|
|00005b60| 6d 70 6c 20 3d 20 43 78 | 41 6d 70 6c 28 20 63 70 |mpl = Cx|Ampl( cp|
|00005b70| 20 29 3b 0a 58 0a 58 09 | 09 63 70 2d 3e 69 6d 20 | );.X.X.|.cp->im |
|00005b80| 3d 20 73 71 72 74 28 20 | 28 61 6d 70 6c 20 2d 20 |= sqrt( |(ampl - |
|00005b90| 63 70 2d 3e 72 65 29 20 | 2f 32 2e 30 20 29 3b 0a |cp->re) |/2.0 );.|
|00005ba0| 58 0a 58 09 09 69 66 20 | 28 20 69 6d 5f 73 69 67 |X.X..if |( im_sig|
|00005bb0| 6e 20 3e 20 30 20 29 0a | 58 09 09 09 63 70 2d 3e |n > 0 ).|X...cp->|
|00005bc0| 72 65 20 3d 20 73 71 72 | 74 28 20 28 61 6d 70 6c |re = sqr|t( (ampl|
|00005bd0| 20 2b 20 63 70 2d 3e 72 | 65 29 20 2f 20 32 2e 30 | + cp->r|e) / 2.0|
|00005be0| 20 29 3b 0a 58 09 09 65 | 6c 73 65 09 09 09 2f 2a | );.X..e|lse.../*|
|00005bf0| 20 69 6d 5f 73 69 67 6e | 20 3c 20 30 20 2a 2f 0a | im_sign| < 0 */.|
|00005c00| 58 09 09 09 63 70 2d 3e | 72 65 20 3d 20 2d 73 71 |X...cp->|re = -sq|
|00005c10| 72 74 28 20 28 61 6d 70 | 6c 20 2b 20 63 70 2d 3e |rt( (amp|l + cp->|
|00005c20| 72 65 29 20 2f 20 32 2e | 30 20 29 3b 0a 58 09 09 |re) / 2.|0 );.X..|
|00005c30| 7d 0a 58 0a 58 09 72 65 | 74 75 72 6e 20 63 70 3b |}.X.X.re|turn cp;|
|00005c40| 0a 58 09 7d 0a 45 4e 44 | 5f 4f 46 5f 63 78 73 71 |.X.}.END|_OF_cxsq|
|00005c50| 72 74 2e 63 0a 69 66 20 | 74 65 73 74 20 31 33 33 |rt.c.if |test 133|
|00005c60| 39 20 2d 6e 65 20 60 77 | 63 20 2d 63 20 3c 63 78 |9 -ne `w|c -c <cx|
|00005c70| 73 71 72 74 2e 63 60 3b | 20 74 68 65 6e 0a 20 20 |sqrt.c`;| then.  |
|00005c80| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 63 |  echo s|har: \"c|
|00005c90| 78 73 71 72 74 2e 63 5c | 22 20 75 6e 70 61 63 6b |xsqrt.c\|" unpack|
|00005ca0| 65 64 20 77 69 74 68 20 | 77 72 6f 6e 67 20 73 69 |ed with |wrong si|
|00005cb0| 7a 65 21 0a 66 69 0a 23 | 20 65 6e 64 20 6f 66 20 |ze!.fi.#| end of |
|00005cc0| 6f 76 65 72 77 72 69 74 | 69 6e 67 20 63 68 65 63 |overwrit|ing chec|
|00005cd0| 6b 0a 66 69 0a 69 66 20 | 74 65 73 74 20 2d 66 20 |k.fi.if |test -f |
|00005ce0| 63 78 73 75 62 2e 63 20 | 2d 61 20 22 24 7b 31 7d |cxsub.c |-a "${1}|
|00005cf0| 22 20 21 3d 20 22 2d 63 | 22 20 3b 20 74 68 65 6e |" != "-c|" ; then|
|00005d00| 20 0a 20 20 65 63 68 6f | 20 73 68 61 72 3a 20 57 | .  echo| shar: W|
|00005d10| 69 6c 6c 20 6e 6f 74 20 | 6f 76 65 72 2d 77 72 69 |ill not |over-wri|
|00005d20| 74 65 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |te exist|ing file|
|00005d30| 20 5c 22 63 78 73 75 62 | 2e 63 5c 22 0a 65 6c 73 | \"cxsub|.c\".els|
|00005d40| 65 0a 65 63 68 6f 20 73 | 68 61 72 3a 20 45 78 74 |e.echo s|har: Ext|
|00005d50| 72 61 63 74 69 6e 67 20 | 5c 22 63 78 73 75 62 2e |racting |\"cxsub.|
|00005d60| 63 5c 22 20 5c 28 33 32 | 30 20 63 68 61 72 61 63 |c\" \(32|0 charac|
|00005d70| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00005d80| 2f 2f 22 20 3e 63 78 73 | 75 62 2e 63 20 3c 3c 27 |//" >cxs|ub.c <<'|
|00005d90| 45 4e 44 5f 4f 46 5f 63 | 78 73 75 62 2e 63 27 0a |END_OF_c|xsub.c'.|
|00005da0| 58 2f 2a 0a 58 09 43 78 | 53 75 62 20 2d 2d 20 73 |X/*.X.Cx|Sub -- s|
|00005db0| 75 62 74 72 61 63 74 20 | 6f 6e 65 20 63 6f 6d 70 |ubtract |one comp|
|00005dc0| 6c 65 78 20 66 72 6f 6d | 20 61 6e 6f 74 68 65 72 |lex from| another|
|00005dd0| 0a 58 0a 58 09 43 78 53 | 75 62 28 20 26 61 2c 20 |.X.X.CxS|ub( &a, |
|00005de0| 26 62 20 29 09 73 75 62 | 74 72 61 63 74 73 20 20 |&b ).sub|tracts  |
|00005df0| 62 20 20 66 72 6f 6d 20 | 20 61 20 20 61 6e 64 20 |b  from | a  and |
|00005e00| 72 65 74 75 72 6e 73 20 | 20 26 61 0a 58 0a 58 09 |returns | &a.X.X.|
|00005e10| 6c 61 73 74 20 65 64 69 | 74 3a 09 38 36 2f 30 31 |last edi|t:.86/01|
|00005e20| 2f 30 34 09 44 20 41 20 | 47 77 79 6e 0a 58 0a 58 |/04.D A |Gwyn.X.X|
|00005e30| 09 53 43 43 53 20 49 44 | 3a 09 40 28 23 29 63 78 |.SCCS ID|:.@(#)cx|
|00005e40| 73 75 62 2e 63 09 31 2e | 31 0a 58 2a 2f 0a 58 0a |sub.c.1.|1.X*/.X.|
|00005e50| 58 23 69 6e 63 6c 75 64 | 65 09 3c 63 6f 6d 70 6c |X#includ|e.<compl|
|00005e60| 65 78 2e 68 3e 0a 58 0a | 58 63 6f 6d 70 6c 65 78 |ex.h>.X.|Xcomplex|
|00005e70| 20 2a 0a 58 43 78 53 75 | 62 28 20 61 70 2c 20 62 | *.XCxSu|b( ap, b|
|00005e80| 70 20 29 0a 58 09 72 65 | 67 69 73 74 65 72 20 63 |p ).X.re|gister c|
|00005e90| 6f 6d 70 6c 65 78 09 2a | 61 70 2c 20 2a 62 70 3b |omplex.*|ap, *bp;|
|00005ea0| 09 2f 2a 20 28 6d 61 79 | 20 63 6f 69 6e 63 69 64 |./* (may| coincid|
|00005eb0| 65 29 20 2a 2f 0a 58 09 | 7b 0a 58 09 61 70 2d 3e |e) */.X.|{.X.ap->|
|00005ec0| 72 65 20 2d 3d 20 62 70 | 2d 3e 72 65 3b 0a 58 09 |re -= bp|->re;.X.|
|00005ed0| 61 70 2d 3e 69 6d 20 2d | 3d 20 62 70 2d 3e 69 6d |ap->im -|= bp->im|
|00005ee0| 3b 0a 58 0a 58 09 72 65 | 74 75 72 6e 20 61 70 3b |;.X.X.re|turn ap;|
|00005ef0| 0a 58 09 7d 0a 45 4e 44 | 5f 4f 46 5f 63 78 73 75 |.X.}.END|_OF_cxsu|
|00005f00| 62 2e 63 0a 69 66 20 74 | 65 73 74 20 33 32 30 20 |b.c.if t|est 320 |
|00005f10| 2d 6e 65 20 60 77 63 20 | 2d 63 20 3c 63 78 73 75 |-ne `wc |-c <cxsu|
|00005f20| 62 2e 63 60 3b 20 74 68 | 65 6e 0a 20 20 20 20 65 |b.c`; th|en.    e|
|00005f30| 63 68 6f 20 73 68 61 72 | 3a 20 5c 22 63 78 73 75 |cho shar|: \"cxsu|
|00005f40| 62 2e 63 5c 22 20 75 6e | 70 61 63 6b 65 64 20 77 |b.c\" un|packed w|
|00005f50| 69 74 68 20 77 72 6f 6e | 67 20 73 69 7a 65 21 0a |ith wron|g size!.|
|00005f60| 66 69 0a 23 20 65 6e 64 | 20 6f 66 20 6f 76 65 72 |fi.# end| of over|
|00005f70| 77 72 69 74 69 6e 67 20 | 63 68 65 63 6b 0a 66 69 |writing |check.fi|
|00005f80| 0a 65 63 68 6f 20 73 68 | 61 72 3a 20 45 6e 64 20 |.echo sh|ar: End |
|00005f90| 6f 66 20 73 68 65 6c 6c | 20 61 72 63 68 69 76 65 |of shell| archive|
|00005fa0| 2e 0a 65 78 69 74 20 30 | 0a 2d 2d 20 0a 0a 52 69 |..exit 0|.-- ..Ri|
|00005fb0| 63 68 20 24 61 6c 7a 09 | 09 09 22 41 6e 67 65 72 |ch $alz.|.."Anger|
|00005fc0| 20 69 73 20 61 6e 20 65 | 6e 65 72 67 79 22 0a 43 | is an e|nergy".C|
|00005fd0| 72 6f 6e 75 73 20 50 72 | 6f 6a 65 63 74 2c 20 42 |ronus Pr|oject, B|
|00005fe0| 42 4e 20 4c 61 62 73 09 | 72 73 61 6c 7a 40 62 62 |BN Labs.|rsalz@bb|
|00005ff0| 6e 2e 63 6f 6d 0a 4d 6f | 64 65 72 61 74 6f 72 2c |n.com.Mo|derator,|
|00006000| 20 63 6f 6d 70 2e 73 6f | 75 72 63 65 73 2e 75 6e | comp.so|urces.un|
|00006010| 69 78 09 73 6f 75 72 63 | 65 73 40 75 75 6e 65 74 |ix.sourc|es@uunet|
|00006020| 2e 75 75 2e 6e 65 74 0a | 2d 2d 20 0a 0a 52 69 63 |.uu.net.|-- ..Ric|
|00006030| 68 20 24 61 6c 7a 09 09 | 09 22 41 6e 67 65 72 20 |h $alz..|."Anger |
|00006040| 69 73 20 61 6e 20 65 6e | 65 72 67 79 22 0a 43 72 |is an en|ergy".Cr|
|00006050| 6f 6e 75 73 20 50 72 6f | 6a 65 63 74 2c 20 42 42 |onus Pro|ject, BB|
|00006060| 4e 20 4c 61 62 73 09 72 | 73 61 6c 7a 40 62 62 6e |N Labs.r|salz@bbn|
|00006070| 2e 63 6f 6d 0a 4d 6f 64 | 65 72 61 74 6f 72 2c 20 |.com.Mod|erator, |
|00006080| 63 6f 6d 70 2e 73 6f 75 | 72 63 65 73 2e 75 6e 69 |comp.sou|rces.uni|
|00006090| 78 09 73 6f 75 72 63 65 | 73 40 75 75 6e 65 74 2e |x.source|s@uunet.|
|000060a0| 75 75 2e 6e 65 74 0a    |                         |uu.net. |        |
+--------+-------------------------+-------------------------+--------+--------+