/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / unix / volume26 / calc / part05

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 75 6e 69 78 0a 46 72 |.sources|.unix.Fr|
|00000020| 6f 6d 3a 20 64 62 65 6c | 6c 40 70 64 61 63 74 2e |om: dbel|l@pdact.|
|00000030| 70 64 2e 6e 65 63 69 73 | 61 2e 6f 7a 2e 61 75 20 |pd.necis|a.oz.au |
|00000040| 28 44 61 76 69 64 20 49 | 2e 20 42 65 6c 6c 29 0a |(David I|. Bell).|
|00000050| 53 75 62 6a 65 63 74 3a | 20 76 32 36 69 30 33 31 |Subject:| v26i031|
|00000060| 3a 20 43 41 4c 43 20 2d | 20 41 6e 20 61 72 62 69 |: CALC -| An arbi|
|00000070| 74 72 61 72 79 20 70 72 | 65 63 69 73 69 6f 6e 20 |trary pr|ecision |
|00000080| 43 2d 6c 69 6b 65 20 63 | 61 6c 63 75 6c 61 74 6f |C-like c|alculato|
|00000090| 72 2c 20 50 61 72 74 30 | 35 2f 32 31 0a 53 65 6e |r, Part0|5/21.Sen|
|000000a0| 64 65 72 3a 20 75 6e 69 | 78 2d 73 6f 75 72 63 65 |der: uni|x-source|
|000000b0| 73 2d 6d 6f 64 65 72 61 | 74 6f 72 40 70 61 2e 64 |s-modera|tor@pa.d|
|000000c0| 65 63 2e 63 6f 6d 0a 41 | 70 70 72 6f 76 65 64 3a |ec.com.A|pproved:|
|000000d0| 20 76 69 78 69 65 40 70 | 61 2e 64 65 63 2e 63 6f | vixie@p|a.dec.co|
|000000e0| 6d 0a 0a 53 75 62 6d 69 | 74 74 65 64 2d 42 79 3a |m..Submi|tted-By:|
|000000f0| 20 64 62 65 6c 6c 40 70 | 64 61 63 74 2e 70 64 2e | dbell@p|dact.pd.|
|00000100| 6e 65 63 69 73 61 2e 6f | 7a 2e 61 75 20 28 44 61 |necisa.o|z.au (Da|
|00000110| 76 69 64 20 49 2e 20 42 | 65 6c 6c 29 0a 50 6f 73 |vid I. B|ell).Pos|
|00000120| 74 69 6e 67 2d 4e 75 6d | 62 65 72 3a 20 56 6f 6c |ting-Num|ber: Vol|
|00000130| 75 6d 65 20 32 36 2c 20 | 49 73 73 75 65 20 33 31 |ume 26, |Issue 31|
|00000140| 0a 41 72 63 68 69 76 65 | 2d 4e 61 6d 65 3a 20 63 |.Archive|-Name: c|
|00000150| 61 6c 63 2f 70 61 72 74 | 30 35 0a 0a 23 21 20 2f |alc/part|05..#! /|
|00000160| 62 69 6e 2f 73 68 0a 23 | 20 54 68 69 73 20 69 73 |bin/sh.#| This is|
|00000170| 20 61 20 73 68 65 6c 6c | 20 61 72 63 68 69 76 65 | a shell| archive|
|00000180| 2e 20 20 52 65 6d 6f 76 | 65 20 61 6e 79 74 68 69 |.  Remov|e anythi|
|00000190| 6e 67 20 62 65 66 6f 72 | 65 20 74 68 69 73 20 6c |ng befor|e this l|
|000001a0| 69 6e 65 2c 20 74 68 65 | 6e 20 75 6e 70 61 63 6b |ine, the|n unpack|
|000001b0| 0a 23 20 69 74 20 62 79 | 20 73 61 76 69 6e 67 20 |.# it by| saving |
|000001c0| 69 74 20 69 6e 74 6f 20 | 61 20 66 69 6c 65 20 61 |it into |a file a|
|000001d0| 6e 64 20 74 79 70 69 6e | 67 20 22 73 68 20 66 69 |nd typin|g "sh fi|
|000001e0| 6c 65 22 2e 20 20 54 6f | 20 6f 76 65 72 77 72 69 |le".  To| overwri|
|000001f0| 74 65 20 65 78 69 73 74 | 69 6e 67 0a 23 20 66 69 |te exist|ing.# fi|
|00000200| 6c 65 73 2c 20 74 79 70 | 65 20 22 73 68 20 66 69 |les, typ|e "sh fi|
|00000210| 6c 65 20 2d 63 22 2e 20 | 20 59 6f 75 20 63 61 6e |le -c". | You can|
|00000220| 20 61 6c 73 6f 20 66 65 | 65 64 20 74 68 69 73 20 | also fe|ed this |
|00000230| 61 73 20 73 74 61 6e 64 | 61 72 64 20 69 6e 70 75 |as stand|ard inpu|
|00000240| 74 20 76 69 61 0a 23 20 | 75 6e 73 68 61 72 2c 20 |t via.# |unshar, |
|00000250| 6f 72 20 62 79 20 74 79 | 70 69 6e 67 20 22 73 68 |or by ty|ping "sh|
|00000260| 20 3c 66 69 6c 65 22 2c | 20 65 2e 67 2e 2e 20 20 | <file",| e.g..  |
|00000270| 49 66 20 74 68 69 73 20 | 61 72 63 68 69 76 65 20 |If this |archive |
|00000280| 69 73 20 63 6f 6d 70 6c | 65 74 65 2c 20 79 6f 75 |is compl|ete, you|
|00000290| 0a 23 20 77 69 6c 6c 20 | 73 65 65 20 74 68 65 20 |.# will |see the |
|000002a0| 66 6f 6c 6c 6f 77 69 6e | 67 20 6d 65 73 73 61 67 |followin|g messag|
|000002b0| 65 20 61 74 20 74 68 65 | 20 65 6e 64 3a 0a 23 09 |e at the| end:.#.|
|000002c0| 09 22 45 6e 64 20 6f 66 | 20 61 72 63 68 69 76 65 |."End of| archive|
|000002d0| 20 35 20 28 6f 66 20 32 | 31 29 2e 22 0a 23 20 43 | 5 (of 2|1).".# C|
|000002e0| 6f 6e 74 65 6e 74 73 3a | 20 20 63 6f 6d 66 75 6e |ontents:|  comfun|
|000002f0| 63 2e 63 20 63 6f 6d 6d | 61 74 68 2e 63 20 66 69 |c.c comm|ath.c fi|
|00000300| 6c 65 2e 63 20 6c 69 73 | 74 66 75 6e 63 2e 63 20 |le.c lis|tfunc.c |
|00000310| 71 6d 6f 64 2e 63 0a 23 | 20 57 72 61 70 70 65 64 |qmod.c.#| Wrapped|
|00000320| 20 62 79 20 64 62 65 6c | 6c 40 65 6c 6d 20 6f 6e | by dbel|l@elm on|
|00000330| 20 54 75 65 20 46 65 62 | 20 32 35 20 31 35 3a 32 | Tue Feb| 25 15:2|
|00000340| 30 3a 35 39 20 31 39 39 | 32 0a 50 41 54 48 3d 2f |0:59 199|2.PATH=/|
|00000350| 62 69 6e 3a 2f 75 73 72 | 2f 62 69 6e 3a 2f 75 73 |bin:/usr|/bin:/us|
|00000360| 72 2f 75 63 62 20 3b 20 | 65 78 70 6f 72 74 20 50 |r/ucb ; |export P|
|00000370| 41 54 48 0a 69 66 20 74 | 65 73 74 20 2d 66 20 27 |ATH.if t|est -f '|
|00000380| 63 6f 6d 66 75 6e 63 2e | 63 27 20 2d 61 20 22 24 |comfunc.|c' -a "$|
|00000390| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|000003a0| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen .  e|cho shar|
|000003b0| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 63 6c 6f 62 62 |: Will n|ot clobb|
|000003c0| 65 72 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |er exist|ing file|
|000003d0| 20 5c 22 27 63 6f 6d 66 | 75 6e 63 2e 63 27 5c 22 | \"'comf|unc.c'\"|
|000003e0| 0a 65 6c 73 65 0a 65 63 | 68 6f 20 73 68 61 72 3a |.else.ec|ho shar:|
|000003f0| 20 45 78 74 72 61 63 74 | 69 6e 67 20 5c 22 27 63 | Extract|ing \"'c|
|00000400| 6f 6d 66 75 6e 63 2e 63 | 27 5c 22 20 5c 28 31 30 |omfunc.c|'\" \(10|
|00000410| 33 33 32 20 63 68 61 72 | 61 63 74 65 72 73 5c 29 |332 char|acters\)|
|00000420| 0a 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 27 |.sed "s/|^X//" >'|
|00000430| 63 6f 6d 66 75 6e 63 2e | 63 27 20 3c 3c 27 45 4e |comfunc.|c' <<'EN|
|00000440| 44 5f 4f 46 5f 46 49 4c | 45 27 0a 58 2f 2a 0a 58 |D_OF_FIL|E'.X/*.X|
|00000450| 20 2a 20 43 6f 70 79 72 | 69 67 68 74 20 28 63 29 | * Copyr|ight (c)|
|00000460| 20 31 39 39 32 20 44 61 | 76 69 64 20 49 2e 20 42 | 1992 Da|vid I. B|
|00000470| 65 6c 6c 0a 58 20 2a 20 | 50 65 72 6d 69 73 73 69 |ell.X * |Permissi|
|00000480| 6f 6e 20 69 73 20 67 72 | 61 6e 74 65 64 20 74 6f |on is gr|anted to|
|00000490| 20 75 73 65 2c 20 64 69 | 73 74 72 69 62 75 74 65 | use, di|stribute|
|000004a0| 2c 20 6f 72 20 6d 6f 64 | 69 66 79 20 74 68 69 73 |, or mod|ify this|
|000004b0| 20 73 6f 75 72 63 65 2c | 0a 58 20 2a 20 70 72 6f | source,|.X * pro|
|000004c0| 76 69 64 65 64 20 74 68 | 61 74 20 74 68 69 73 20 |vided th|at this |
|000004d0| 63 6f 70 79 72 69 67 68 | 74 20 6e 6f 74 69 63 65 |copyrigh|t notice|
|000004e0| 20 72 65 6d 61 69 6e 73 | 20 69 6e 74 61 63 74 2e | remains| intact.|
|000004f0| 0a 58 20 2a 0a 58 20 2a | 20 45 78 74 65 6e 64 65 |.X *.X *| Extende|
|00000500| 64 20 70 72 65 63 69 73 | 69 6f 6e 20 63 6f 6d 70 |d precis|ion comp|
|00000510| 6c 65 78 20 61 72 69 74 | 68 6d 65 74 69 63 20 6e |lex arit|hmetic n|
|00000520| 6f 6e 2d 70 72 69 6d 69 | 74 69 76 65 20 72 6f 75 |on-primi|tive rou|
|00000530| 74 69 6e 65 73 0a 58 20 | 2a 2f 0a 58 0a 58 23 69 |tines.X |*/.X.X#i|
|00000540| 6e 63 6c 75 64 65 20 22 | 63 61 6c 63 2e 68 22 0a |nclude "|calc.h".|
|00000550| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 52 6f 75 6e |X.X.X/*.|X * Roun|
|00000560| 64 20 61 20 63 6f 6d 70 | 6c 65 78 20 6e 75 6d 62 |d a comp|lex numb|
|00000570| 65 72 20 74 6f 20 74 68 | 65 20 73 70 65 63 69 66 |er to th|e specif|
|00000580| 69 65 64 20 6e 75 6d 62 | 65 72 20 6f 66 20 64 65 |ied numb|er of de|
|00000590| 63 69 6d 61 6c 20 70 6c | 61 63 65 73 2e 0a 58 20 |cimal pl|aces..X |
|000005a0| 2a 20 54 68 69 73 20 73 | 69 6d 70 6c 79 20 6d 65 |* This s|imply me|
|000005b0| 61 6e 73 20 74 6f 20 72 | 6f 75 6e 64 20 65 61 63 |ans to r|ound eac|
|000005c0| 68 20 6f 66 20 74 68 65 | 20 63 6f 6d 70 6f 6e 65 |h of the| compone|
|000005d0| 6e 74 73 20 6f 66 20 74 | 68 65 20 6e 75 6d 62 65 |nts of t|he numbe|
|000005e0| 72 2e 0a 58 20 2a 20 5a | 65 72 6f 20 64 65 63 69 |r..X * Z|ero deci|
|000005f0| 6d 61 6c 20 70 6c 61 63 | 65 73 20 6d 65 61 6e 73 |mal plac|es means|
|00000600| 20 72 6f 75 6e 64 20 74 | 6f 20 74 68 65 20 6e 65 | round t|o the ne|
|00000610| 61 72 65 73 74 20 63 6f | 6d 70 6c 65 78 20 69 6e |arest co|mplex in|
|00000620| 74 65 67 65 72 2e 0a 58 | 20 2a 2f 0a 58 43 4f 4d |teger..X| */.XCOM|
|00000630| 50 4c 45 58 20 2a 0a 58 | 63 72 6f 75 6e 64 28 63 |PLEX *.X|cround(c|
|00000640| 2c 20 70 6c 61 63 65 73 | 29 0a 58 09 43 4f 4d 50 |, places|).X.COMP|
|00000650| 4c 45 58 20 2a 63 3b 0a | 58 09 6c 6f 6e 67 20 70 |LEX *c;.|X.long p|
|00000660| 6c 61 63 65 73 3b 0a 58 | 7b 0a 58 09 43 4f 4d 50 |laces;.X|{.X.COMP|
|00000670| 4c 45 58 20 2a 72 65 73 | 3b 09 09 2f 2a 20 72 65 |LEX *res|;../* re|
|00000680| 73 75 6c 74 20 2a 2f 0a | 58 0a 58 09 72 65 73 20 |sult */.|X.X.res |
|00000690| 3d 20 63 6f 6d 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |= comall|oc();.X.|
|000006a0| 72 65 73 2d 3e 72 65 61 | 6c 20 3d 20 71 72 6f 75 |res->rea|l = qrou|
|000006b0| 6e 64 28 63 2d 3e 72 65 | 61 6c 2c 20 70 6c 61 63 |nd(c->re|al, plac|
|000006c0| 65 73 29 3b 0a 58 09 72 | 65 73 2d 3e 69 6d 61 67 |es);.X.r|es->imag|
|000006d0| 20 3d 20 71 72 6f 75 6e | 64 28 63 2d 3e 69 6d 61 | = qroun|d(c->ima|
|000006e0| 67 2c 20 70 6c 61 63 65 | 73 29 3b 0a 58 09 72 65 |g, place|s);.X.re|
|000006f0| 74 75 72 6e 20 72 65 73 | 3b 0a 58 7d 0a 58 0a 58 |turn res|;.X}.X.X|
|00000700| 0a 58 2f 2a 0a 58 20 2a | 20 52 6f 75 6e 64 20 61 |.X/*.X *| Round a|
|00000710| 20 63 6f 6d 70 6c 65 78 | 20 6e 75 6d 62 65 72 20 | complex| number |
|00000720| 74 6f 20 74 68 65 20 73 | 70 65 63 69 66 69 65 64 |to the s|pecified|
|00000730| 20 6e 75 6d 62 65 72 20 | 6f 66 20 62 69 6e 61 72 | number |of binar|
|00000740| 79 20 64 65 63 69 6d 61 | 6c 20 70 6c 61 63 65 73 |y decima|l places|
|00000750| 2e 0a 58 20 2a 20 54 68 | 69 73 20 73 69 6d 70 6c |..X * Th|is simpl|
|00000760| 79 20 6d 65 61 6e 73 20 | 74 6f 20 72 6f 75 6e 64 |y means |to round|
|00000770| 20 65 61 63 68 20 6f 66 | 20 74 68 65 20 63 6f 6d | each of| the com|
|00000780| 70 6f 6e 65 6e 74 73 20 | 6f 66 20 74 68 65 20 6e |ponents |of the n|
|00000790| 75 6d 62 65 72 2e 0a 58 | 20 2a 20 5a 65 72 6f 20 |umber..X| * Zero |
|000007a0| 62 69 6e 61 72 79 20 70 | 6c 61 63 65 73 20 6d 65 |binary p|laces me|
|000007b0| 61 6e 73 20 72 6f 75 6e | 64 20 74 6f 20 74 68 65 |ans roun|d to the|
|000007c0| 20 6e 65 61 72 65 73 74 | 20 63 6f 6d 70 6c 65 78 | nearest| complex|
|000007d0| 20 69 6e 74 65 67 65 72 | 2e 0a 58 20 2a 2f 0a 58 | integer|..X */.X|
|000007e0| 43 4f 4d 50 4c 45 58 20 | 2a 0a 58 63 62 72 6f 75 |COMPLEX |*.Xcbrou|
|000007f0| 6e 64 28 63 2c 20 70 6c | 61 63 65 73 29 0a 58 09 |nd(c, pl|aces).X.|
|00000800| 43 4f 4d 50 4c 45 58 20 | 2a 63 3b 0a 58 09 6c 6f |COMPLEX |*c;.X.lo|
|00000810| 6e 67 20 70 6c 61 63 65 | 73 3b 0a 58 7b 0a 58 09 |ng place|s;.X{.X.|
|00000820| 43 4f 4d 50 4c 45 58 20 | 2a 72 65 73 3b 09 09 2f |COMPLEX |*res;../|
|00000830| 2a 20 72 65 73 75 6c 74 | 20 2a 2f 0a 58 0a 58 09 |* result| */.X.X.|
|00000840| 72 65 73 20 3d 20 63 6f | 6d 61 6c 6c 6f 63 28 29 |res = co|malloc()|
|00000850| 3b 0a 58 09 72 65 73 2d | 3e 72 65 61 6c 20 3d 20 |;.X.res-|>real = |
|00000860| 71 62 72 6f 75 6e 64 28 | 63 2d 3e 72 65 61 6c 2c |qbround(|c->real,|
|00000870| 20 70 6c 61 63 65 73 29 | 3b 0a 58 09 72 65 73 2d | places)|;.X.res-|
|00000880| 3e 69 6d 61 67 20 3d 20 | 71 62 72 6f 75 6e 64 28 |>imag = |qbround(|
|00000890| 63 2d 3e 69 6d 61 67 2c | 20 70 6c 61 63 65 73 29 |c->imag,| places)|
|000008a0| 3b 0a 58 09 72 65 74 75 | 72 6e 20 72 65 73 3b 0a |;.X.retu|rn res;.|
|000008b0| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 43 |X}.X.X.X|/*.X * C|
|000008c0| 6f 6d 70 75 74 65 20 74 | 68 65 20 72 65 73 75 6c |ompute t|he resul|
|000008d0| 74 20 6f 66 20 72 61 69 | 73 69 6e 67 20 61 20 63 |t of rai|sing a c|
|000008e0| 6f 6d 70 6c 65 78 20 6e | 75 6d 62 65 72 20 74 6f |omplex n|umber to|
|000008f0| 20 61 6e 20 69 6e 74 65 | 67 65 72 20 70 6f 77 65 | an inte|ger powe|
|00000900| 72 2e 0a 58 20 2a 2f 0a | 58 43 4f 4d 50 4c 45 58 |r..X */.|XCOMPLEX|
|00000910| 20 2a 0a 58 63 70 6f 77 | 69 28 63 2c 20 71 29 0a | *.Xcpow|i(c, q).|
|00000920| 58 09 43 4f 4d 50 4c 45 | 58 20 2a 63 3b 09 09 2f |X.COMPLE|X *c;../|
|00000930| 2a 20 63 6f 6d 70 6c 65 | 78 20 6e 75 6d 62 65 72 |* comple|x number|
|00000940| 20 74 6f 20 62 65 20 72 | 61 69 73 65 64 20 2a 2f | to be r|aised */|
|00000950| 0a 58 09 4e 55 4d 42 45 | 52 20 2a 71 3b 09 09 2f |.X.NUMBE|R *q;../|
|00000960| 2a 20 70 6f 77 65 72 20 | 74 6f 20 72 61 69 73 65 |* power |to raise|
|00000970| 20 69 74 20 74 6f 20 2a | 2f 0a 58 7b 0a 58 09 43 | it to *|/.X{.X.C|
|00000980| 4f 4d 50 4c 45 58 20 2a | 74 6d 70 2c 20 2a 72 65 |OMPLEX *|tmp, *re|
|00000990| 73 3b 09 2f 2a 20 74 65 | 6d 70 6f 72 61 72 79 20 |s;./* te|mporary |
|000009a0| 76 61 6c 75 65 73 20 2a | 2f 0a 58 09 6c 6f 6e 67 |values *|/.X.long|
|000009b0| 20 70 6f 77 65 72 3b 09 | 09 2f 2a 20 70 6f 77 65 | power;.|./* powe|
|000009c0| 72 20 74 6f 20 72 61 69 | 73 65 20 74 6f 20 2a 2f |r to rai|se to */|
|000009d0| 0a 58 09 75 6e 73 69 67 | 6e 65 64 20 6c 6f 6e 67 |.X.unsig|ned long|
|000009e0| 20 62 69 74 3b 09 2f 2a | 20 63 75 72 72 65 6e 74 | bit;./*| current|
|000009f0| 20 62 69 74 20 76 61 6c | 75 65 20 2a 2f 0a 58 09 | bit val|ue */.X.|
|00000a00| 69 6e 74 20 73 69 67 6e | 3b 0a 58 0a 58 09 69 66 |int sign|;.X.X.if|
|00000a10| 20 28 71 69 73 66 72 61 | 63 28 71 29 29 0a 58 09 | (qisfra|c(q)).X.|
|00000a20| 09 65 72 72 6f 72 28 22 | 52 61 69 73 69 6e 67 20 |.error("|Raising |
|00000a30| 6e 75 6d 62 65 72 20 74 | 6f 20 6e 6f 6e 2d 69 6e |number t|o non-in|
|00000a40| 74 65 67 72 61 6c 20 70 | 6f 77 65 72 22 29 3b 0a |tegral p|ower");.|
|00000a50| 58 09 69 66 20 28 69 73 | 62 69 67 28 71 2d 3e 6e |X.if (is|big(q->n|
|00000a60| 75 6d 29 29 0a 58 09 09 | 65 72 72 6f 72 28 22 52 |um)).X..|error("R|
|00000a70| 61 69 73 69 6e 67 20 6e | 75 6d 62 65 72 20 74 6f |aising n|umber to|
|00000a80| 20 76 65 72 79 20 6c 61 | 72 67 65 20 70 6f 77 65 | very la|rge powe|
|00000a90| 72 22 29 3b 0a 58 09 70 | 6f 77 65 72 20 3d 20 28 |r");.X.p|ower = (|
|00000aa0| 69 73 74 69 6e 79 28 71 | 2d 3e 6e 75 6d 29 20 3f |istiny(q|->num) ?|
|00000ab0| 20 7a 31 74 6f 6c 28 71 | 2d 3e 6e 75 6d 29 20 3a | z1tol(q|->num) :|
|00000ac0| 20 7a 32 74 6f 6c 28 71 | 2d 3e 6e 75 6d 29 29 3b | z2tol(q|->num));|
|00000ad0| 0a 58 09 69 66 20 28 63 | 69 73 7a 65 72 6f 28 63 |.X.if (c|iszero(c|
|00000ae0| 29 20 26 26 20 28 70 6f | 77 65 72 20 3d 3d 20 30 |) && (po|wer == 0|
|00000af0| 29 29 0a 58 09 09 65 72 | 72 6f 72 28 22 52 61 69 |)).X..er|ror("Rai|
|00000b00| 73 69 6e 67 20 7a 65 72 | 6f 20 74 6f 20 7a 65 72 |sing zer|o to zer|
|00000b10| 6f 74 68 20 70 6f 77 65 | 72 22 29 3b 0a 58 09 73 |oth powe|r");.X.s|
|00000b20| 69 67 6e 20 3d 20 31 3b | 0a 58 09 69 66 20 28 71 |ign = 1;|.X.if (q|
|00000b30| 69 73 6e 65 67 28 71 29 | 29 0a 58 09 09 73 69 67 |isneg(q)|).X..sig|
|00000b40| 6e 20 3d 20 2d 31 3b 0a | 58 09 2f 2a 0a 58 09 20 |n = -1;.|X./*.X. |
|00000b50| 2a 20 48 61 6e 64 6c 65 | 20 73 6f 6d 65 20 6c 6f |* Handle| some lo|
|00000b60| 77 20 70 6f 77 65 72 73 | 20 73 70 65 63 69 61 6c |w powers| special|
|00000b70| 6c 79 0a 58 09 20 2a 2f | 0a 58 09 69 66 20 28 70 |ly.X. */|.X.if (p|
|00000b80| 6f 77 65 72 20 3c 3d 20 | 34 29 20 7b 0a 58 09 09 |ower <= |4) {.X..|
|00000b90| 73 77 69 74 63 68 20 28 | 28 69 6e 74 29 20 28 70 |switch (|(int) (p|
|00000ba0| 6f 77 65 72 20 2a 20 73 | 69 67 6e 29 29 20 7b 0a |ower * s|ign)) {.|
|00000bb0| 58 09 09 09 63 61 73 65 | 20 30 3a 0a 58 09 09 09 |X...case| 0:.X...|
|00000bc0| 09 72 65 74 75 72 6e 20 | 63 6c 69 6e 6b 28 26 5f |.return |clink(&_|
|00000bd0| 63 6f 6e 65 5f 29 3b 0a | 58 09 09 09 63 61 73 65 |cone_);.|X...case|
|00000be0| 20 31 3a 0a 58 09 09 09 | 09 72 65 74 75 72 6e 20 | 1:.X...|.return |
|00000bf0| 63 6c 69 6e 6b 28 63 29 | 3b 0a 58 09 09 09 63 61 |clink(c)|;.X...ca|
|00000c00| 73 65 20 2d 31 3a 0a 58 | 09 09 09 09 72 65 74 75 |se -1:.X|....retu|
|00000c10| 72 6e 20 63 69 6e 76 28 | 63 29 3b 0a 58 09 09 09 |rn cinv(|c);.X...|
|00000c20| 63 61 73 65 20 32 3a 0a | 58 09 09 09 09 72 65 74 |case 2:.|X....ret|
|00000c30| 75 72 6e 20 63 73 71 75 | 61 72 65 28 63 29 3b 0a |urn csqu|are(c);.|
|00000c40| 58 09 09 09 63 61 73 65 | 20 2d 32 3a 0a 58 09 09 |X...case| -2:.X..|
|00000c50| 09 09 74 6d 70 20 3d 20 | 63 73 71 75 61 72 65 28 |..tmp = |csquare(|
|00000c60| 63 29 3b 0a 58 09 09 09 | 09 72 65 73 20 3d 20 63 |c);.X...|.res = c|
|00000c70| 69 6e 76 28 74 6d 70 29 | 3b 0a 58 09 09 09 09 63 |inv(tmp)|;.X....c|
|00000c80| 6f 6d 66 72 65 65 28 74 | 6d 70 29 3b 0a 58 09 09 |omfree(t|mp);.X..|
|00000c90| 09 09 72 65 74 75 72 6e | 20 72 65 73 3b 0a 58 09 |..return| res;.X.|
|00000ca0| 09 09 63 61 73 65 20 33 | 3a 0a 58 09 09 09 09 74 |..case 3|:.X....t|
|00000cb0| 6d 70 20 3d 20 63 73 71 | 75 61 72 65 28 63 29 3b |mp = csq|uare(c);|
|00000cc0| 0a 58 09 09 09 09 72 65 | 73 20 3d 20 63 6d 75 6c |.X....re|s = cmul|
|00000cd0| 28 63 2c 20 74 6d 70 29 | 3b 0a 58 09 09 09 09 63 |(c, tmp)|;.X....c|
|00000ce0| 6f 6d 66 72 65 65 28 74 | 6d 70 29 3b 0a 58 09 09 |omfree(t|mp);.X..|
|00000cf0| 09 09 72 65 74 75 72 6e | 20 72 65 73 3b 0a 58 09 |..return| res;.X.|
|00000d00| 09 09 63 61 73 65 20 34 | 3a 0a 58 09 09 09 09 74 |..case 4|:.X....t|
|00000d10| 6d 70 20 3d 20 63 73 71 | 75 61 72 65 28 63 29 3b |mp = csq|uare(c);|
|00000d20| 0a 58 09 09 09 09 72 65 | 73 20 3d 20 63 73 71 75 |.X....re|s = csqu|
|00000d30| 61 72 65 28 74 6d 70 29 | 3b 0a 58 09 09 09 09 63 |are(tmp)|;.X....c|
|00000d40| 6f 6d 66 72 65 65 28 74 | 6d 70 29 3b 0a 58 09 09 |omfree(t|mp);.X..|
|00000d50| 09 09 72 65 74 75 72 6e | 20 72 65 73 3b 0a 58 09 |..return| res;.X.|
|00000d60| 09 7d 0a 58 09 7d 0a 58 | 09 2f 2a 0a 58 09 20 2a |.}.X.}.X|./*.X. *|
|00000d70| 20 43 6f 6d 70 75 74 65 | 20 74 68 65 20 70 6f 77 | Compute| the pow|
|00000d80| 65 72 20 62 79 20 73 71 | 75 61 72 69 6e 67 20 61 |er by sq|uaring a|
|00000d90| 6e 64 20 6d 75 6c 74 69 | 70 6c 79 69 6e 67 2e 0a |nd multi|plying..|
|00000da0| 58 09 20 2a 20 54 68 69 | 73 20 75 73 65 73 20 74 |X. * Thi|s uses t|
|00000db0| 68 65 20 6c 65 66 74 20 | 74 6f 20 72 69 67 68 74 |he left |to right|
|00000dc0| 20 6d 65 74 68 6f 64 20 | 6f 66 20 70 6f 77 65 72 | method |of power|
|00000dd0| 20 72 61 69 73 69 6e 67 | 2e 0a 58 09 20 2a 2f 0a | raising|..X. */.|
|00000de0| 58 09 62 69 74 20 3d 20 | 54 4f 50 46 55 4c 4c 3b |X.bit = |TOPFULL;|
|00000df0| 0a 58 09 77 68 69 6c 65 | 20 28 28 62 69 74 20 26 |.X.while| ((bit &|
|00000e00| 20 70 6f 77 65 72 29 20 | 3d 3d 20 30 29 0a 58 09 | power) |== 0).X.|
|00000e10| 09 62 69 74 20 3e 3e 3d | 20 31 4c 3b 0a 58 09 62 |.bit >>=| 1L;.X.b|
|00000e20| 69 74 20 3e 3e 3d 20 31 | 4c 3b 0a 58 09 72 65 73 |it >>= 1|L;.X.res|
|00000e30| 20 3d 20 63 73 71 75 61 | 72 65 28 63 29 3b 0a 58 | = csqua|re(c);.X|
|00000e40| 09 69 66 20 28 62 69 74 | 20 26 20 70 6f 77 65 72 |.if (bit| & power|
|00000e50| 29 20 7b 0a 58 09 09 74 | 6d 70 20 3d 20 63 6d 75 |) {.X..t|mp = cmu|
|00000e60| 6c 28 72 65 73 2c 20 63 | 29 3b 0a 58 09 09 63 6f |l(res, c|);.X..co|
|00000e70| 6d 66 72 65 65 28 72 65 | 73 29 3b 0a 58 09 09 72 |mfree(re|s);.X..r|
|00000e80| 65 73 20 3d 20 74 6d 70 | 3b 0a 58 09 7d 0a 58 09 |es = tmp|;.X.}.X.|
|00000e90| 62 69 74 20 3e 3e 3d 20 | 31 4c 3b 0a 58 09 77 68 |bit >>= |1L;.X.wh|
|00000ea0| 69 6c 65 20 28 62 69 74 | 29 20 7b 0a 58 09 09 74 |ile (bit|) {.X..t|
|00000eb0| 6d 70 20 3d 20 63 73 71 | 75 61 72 65 28 72 65 73 |mp = csq|uare(res|
|00000ec0| 29 3b 0a 58 09 09 63 6f | 6d 66 72 65 65 28 72 65 |);.X..co|mfree(re|
|00000ed0| 73 29 3b 0a 58 09 09 72 | 65 73 20 3d 20 74 6d 70 |s);.X..r|es = tmp|
|00000ee0| 3b 0a 58 09 09 69 66 20 | 28 62 69 74 20 26 20 70 |;.X..if |(bit & p|
|00000ef0| 6f 77 65 72 29 20 7b 0a | 58 09 09 09 74 6d 70 20 |ower) {.|X...tmp |
|00000f00| 3d 20 63 6d 75 6c 28 72 | 65 73 2c 20 63 29 3b 0a |= cmul(r|es, c);.|
|00000f10| 58 09 09 09 63 6f 6d 66 | 72 65 65 28 72 65 73 29 |X...comf|ree(res)|
|00000f20| 3b 0a 58 09 09 09 72 65 | 73 20 3d 20 74 6d 70 3b |;.X...re|s = tmp;|
|00000f30| 0a 58 09 09 7d 0a 58 09 | 09 62 69 74 20 3e 3e 3d |.X..}.X.|.bit >>=|
|00000f40| 20 31 4c 3b 0a 58 09 7d | 0a 58 09 69 66 20 28 73 | 1L;.X.}|.X.if (s|
|00000f50| 69 67 6e 20 3c 20 30 29 | 20 7b 0a 58 09 09 74 6d |ign < 0)| {.X..tm|
|00000f60| 70 20 3d 20 63 69 6e 76 | 28 72 65 73 29 3b 0a 58 |p = cinv|(res);.X|
|00000f70| 09 09 63 6f 6d 66 72 65 | 65 28 72 65 73 29 3b 0a |..comfre|e(res);.|
|00000f80| 58 09 09 72 65 73 20 3d | 20 74 6d 70 3b 0a 58 09 |X..res =| tmp;.X.|
|00000f90| 7d 0a 58 09 72 65 74 75 | 72 6e 20 72 65 73 3b 0a |}.X.retu|rn res;.|
|00000fa0| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 43 |X}.X.X.X|/*.X * C|
|00000fb0| 61 6c 63 75 6c 61 74 65 | 20 74 68 65 20 73 71 75 |alculate| the squ|
|00000fc0| 61 72 65 20 72 6f 6f 74 | 20 6f 66 20 61 20 63 6f |are root| of a co|
|00000fd0| 6d 70 6c 65 78 20 6e 75 | 6d 62 65 72 2c 20 77 69 |mplex nu|mber, wi|
|00000fe0| 74 68 20 65 61 63 68 20 | 63 6f 6d 70 6f 6e 65 6e |th each |componen|
|00000ff0| 74 0a 58 20 2a 20 77 69 | 74 68 69 6e 20 74 68 65 |t.X * wi|thin the|
|00001000| 20 73 70 65 63 69 66 69 | 65 64 20 65 72 72 6f 72 | specifi|ed error|
|00001010| 2e 20 54 68 69 73 20 75 | 73 65 73 20 74 68 65 20 |. This u|ses the |
|00001020| 66 6f 72 6d 75 6c 61 3a | 0a 58 20 2a 09 73 71 72 |formula:|.X *.sqr|
|00001030| 74 28 61 2b 62 69 29 20 | 3d 20 73 71 72 74 28 28 |t(a+bi) |= sqrt((|
|00001040| 61 2b 73 71 72 74 28 61 | 5e 32 2b 62 5e 32 29 29 |a+sqrt(a|^2+b^2))|
|00001050| 2f 32 29 20 2b 20 73 71 | 72 74 28 28 2d 61 2b 73 |/2) + sq|rt((-a+s|
|00001060| 71 72 74 28 61 5e 32 2b | 62 5e 32 29 29 2f 32 29 |qrt(a^2+|b^2))/2)|
|00001070| 69 2e 0a 58 20 2a 2f 0a | 58 43 4f 4d 50 4c 45 58 |i..X */.|XCOMPLEX|
|00001080| 20 2a 0a 58 63 73 71 72 | 74 28 63 2c 20 65 70 73 | *.Xcsqr|t(c, eps|
|00001090| 69 6c 6f 6e 29 0a 58 09 | 43 4f 4d 50 4c 45 58 20 |ilon).X.|COMPLEX |
|000010a0| 2a 63 3b 0a 58 09 4e 55 | 4d 42 45 52 20 2a 65 70 |*c;.X.NU|MBER *ep|
|000010b0| 73 69 6c 6f 6e 3b 0a 58 | 7b 0a 58 09 43 4f 4d 50 |silon;.X|{.X.COMP|
|000010c0| 4c 45 58 20 2a 72 3b 0a | 58 09 4e 55 4d 42 45 52 |LEX *r;.|X.NUMBER|
|000010d0| 20 2a 68 79 70 6f 74 2c | 20 2a 74 6d 70 31 2c 20 | *hypot,| *tmp1, |
|000010e0| 2a 74 6d 70 32 3b 0a 58 | 0a 58 09 69 66 20 28 63 |*tmp2;.X|.X.if (c|
|000010f0| 69 73 7a 65 72 6f 28 63 | 29 20 7c 7c 20 63 69 73 |iszero(c|) || cis|
|00001100| 6f 6e 65 28 63 29 29 0a | 58 09 09 72 65 74 75 72 |one(c)).|X..retur|
|00001110| 6e 20 63 6c 69 6e 6b 28 | 63 29 3b 0a 58 09 72 20 |n clink(|c);.X.r |
|00001120| 3d 20 63 6f 6d 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |= comall|oc();.X.|
|00001130| 69 66 20 28 63 69 73 72 | 65 61 6c 28 63 29 29 20 |if (cisr|eal(c)) |
|00001140| 7b 0a 58 09 09 69 66 20 | 28 21 71 69 73 6e 65 67 |{.X..if |(!qisneg|
|00001150| 28 63 2d 3e 72 65 61 6c | 29 29 20 7b 0a 58 09 09 |(c->real|)) {.X..|
|00001160| 09 72 2d 3e 72 65 61 6c | 20 3d 20 71 73 71 72 74 |.r->real| = qsqrt|
|00001170| 28 63 2d 3e 72 65 61 6c | 2c 20 65 70 73 69 6c 6f |(c->real|, epsilo|
|00001180| 6e 29 3b 0a 58 09 09 09 | 72 65 74 75 72 6e 20 72 |n);.X...|return r|
|00001190| 3b 0a 58 09 09 7d 0a 58 | 09 09 74 6d 70 31 20 3d |;.X..}.X|..tmp1 =|
|000011a0| 20 71 6e 65 67 28 63 2d | 3e 72 65 61 6c 29 3b 0a | qneg(c-|>real);.|
|000011b0| 58 09 09 72 2d 3e 69 6d | 61 67 20 3d 20 71 73 71 |X..r->im|ag = qsq|
|000011c0| 72 74 28 74 6d 70 31 2c | 20 65 70 73 69 6c 6f 6e |rt(tmp1,| epsilon|
|000011d0| 29 3b 0a 58 09 09 71 66 | 72 65 65 28 74 6d 70 31 |);.X..qf|ree(tmp1|
|000011e0| 29 3b 0a 58 09 09 72 65 | 74 75 72 6e 20 72 3b 0a |);.X..re|turn r;.|
|000011f0| 58 09 7d 0a 58 09 68 79 | 70 6f 74 20 3d 20 71 68 |X.}.X.hy|pot = qh|
|00001200| 79 70 6f 74 28 63 2d 3e | 72 65 61 6c 2c 20 63 2d |ypot(c->|real, c-|
|00001210| 3e 69 6d 61 67 2c 20 65 | 70 73 69 6c 6f 6e 29 3b |>imag, e|psilon);|
|00001220| 0a 58 09 74 6d 70 31 20 | 3d 20 71 61 64 64 28 68 |.X.tmp1 |= qadd(h|
|00001230| 79 70 6f 74 2c 20 63 2d | 3e 72 65 61 6c 29 3b 0a |ypot, c-|>real);.|
|00001240| 58 09 74 6d 70 32 20 3d | 20 71 73 63 61 6c 65 28 |X.tmp2 =| qscale(|
|00001250| 74 6d 70 31 2c 20 2d 31 | 4c 29 3b 0a 58 09 71 66 |tmp1, -1|L);.X.qf|
|00001260| 72 65 65 28 74 6d 70 31 | 29 3b 0a 58 09 72 2d 3e |ree(tmp1|);.X.r->|
|00001270| 72 65 61 6c 20 3d 20 71 | 73 71 72 74 28 74 6d 70 |real = q|sqrt(tmp|
|00001280| 32 2c 20 65 70 73 69 6c | 6f 6e 29 3b 0a 58 09 71 |2, epsil|on);.X.q|
|00001290| 66 72 65 65 28 74 6d 70 | 32 29 3b 0a 58 09 74 6d |free(tmp|2);.X.tm|
|000012a0| 70 31 20 3d 20 71 73 75 | 62 28 68 79 70 6f 74 2c |p1 = qsu|b(hypot,|
|000012b0| 20 63 2d 3e 72 65 61 6c | 29 3b 0a 58 09 71 66 72 | c->real|);.X.qfr|
|000012c0| 65 65 28 68 79 70 6f 74 | 29 3b 0a 58 09 74 6d 70 |ee(hypot|);.X.tmp|
|000012d0| 32 20 3d 20 71 73 63 61 | 6c 65 28 74 6d 70 31 2c |2 = qsca|le(tmp1,|
|000012e0| 20 2d 31 4c 29 3b 0a 58 | 09 71 66 72 65 65 28 74 | -1L);.X|.qfree(t|
|000012f0| 6d 70 31 29 3b 0a 58 09 | 74 6d 70 31 20 3d 20 71 |mp1);.X.|tmp1 = q|
|00001300| 73 71 72 74 28 74 6d 70 | 32 2c 20 65 70 73 69 6c |sqrt(tmp|2, epsil|
|00001310| 6f 6e 29 3b 0a 58 09 71 | 66 72 65 65 28 74 6d 70 |on);.X.q|free(tmp|
|00001320| 32 29 3b 0a 58 09 69 66 | 20 28 71 69 73 6e 65 67 |2);.X.if| (qisneg|
|00001330| 28 63 2d 3e 69 6d 61 67 | 29 29 20 7b 0a 58 09 09 |(c->imag|)) {.X..|
|00001340| 74 6d 70 32 20 3d 20 71 | 6e 65 67 28 74 6d 70 31 |tmp2 = q|neg(tmp1|
|00001350| 29 3b 0a 58 09 09 71 66 | 72 65 65 28 74 6d 70 31 |);.X..qf|ree(tmp1|
|00001360| 29 3b 0a 58 09 09 74 6d | 70 31 20 3d 20 74 6d 70 |);.X..tm|p1 = tmp|
|00001370| 32 3b 0a 58 09 7d 0a 58 | 09 72 2d 3e 69 6d 61 67 |2;.X.}.X|.r->imag|
|00001380| 20 3d 20 74 6d 70 31 3b | 0a 58 09 72 65 74 75 72 | = tmp1;|.X.retur|
|00001390| 6e 20 72 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 0a |n r;.X}.|X.X.X/*.|
|000013a0| 58 20 2a 20 54 61 6b 65 | 20 74 68 65 20 4e 74 68 |X * Take| the Nth|
|000013b0| 20 72 6f 6f 74 20 6f 66 | 20 61 20 63 6f 6d 70 6c | root of| a compl|
|000013c0| 65 78 20 6e 75 6d 62 65 | 72 2c 20 77 68 65 72 65 |ex numbe|r, where|
|000013d0| 20 4e 20 69 73 20 61 20 | 70 6f 73 69 74 69 76 65 | N is a |positive|
|000013e0| 20 69 6e 74 65 67 65 72 | 2e 0a 58 20 2a 20 45 61 | integer|..X * Ea|
|000013f0| 63 68 20 63 6f 6d 70 6f | 6e 65 6e 74 20 6f 66 20 |ch compo|nent of |
|00001400| 74 68 65 20 72 65 73 75 | 6c 74 20 69 73 20 77 69 |the resu|lt is wi|
|00001410| 74 68 69 6e 20 74 68 65 | 20 73 70 65 63 69 66 69 |thin the| specifi|
|00001420| 65 64 20 65 72 72 6f 72 | 2e 0a 58 20 2a 2f 0a 58 |ed error|..X */.X|
|00001430| 43 4f 4d 50 4c 45 58 20 | 2a 0a 58 63 72 6f 6f 74 |COMPLEX |*.Xcroot|
|00001440| 28 63 2c 20 71 2c 20 65 | 70 73 69 6c 6f 6e 29 0a |(c, q, e|psilon).|
|00001450| 58 09 43 4f 4d 50 4c 45 | 58 20 2a 63 3b 0a 58 09 |X.COMPLE|X *c;.X.|
|00001460| 4e 55 4d 42 45 52 20 2a | 71 2c 20 2a 65 70 73 69 |NUMBER *|q, *epsi|
|00001470| 6c 6f 6e 3b 0a 58 7b 0a | 58 09 43 4f 4d 50 4c 45 |lon;.X{.|X.COMPLE|
|00001480| 58 20 2a 72 3b 0a 58 09 | 4e 55 4d 42 45 52 20 2a |X *r;.X.|NUMBER *|
|00001490| 61 32 70 62 32 2c 20 2a | 72 6f 6f 74 2c 20 2a 74 |a2pb2, *|root, *t|
|000014a0| 6d 70 31 2c 20 2a 74 6d | 70 32 2c 20 2a 65 70 73 |mp1, *tm|p2, *eps|
|000014b0| 69 6c 6f 6e 32 3b 0a 58 | 0a 58 09 69 66 20 28 71 |ilon2;.X|.X.if (q|
|000014c0| 69 73 6e 65 67 28 71 29 | 20 7c 7c 20 71 69 73 7a |isneg(q)| || qisz|
|000014d0| 65 72 6f 28 71 29 20 7c | 7c 20 71 69 73 66 72 61 |ero(q) ||| qisfra|
|000014e0| 63 28 71 29 29 0a 58 09 | 09 65 72 72 6f 72 28 22 |c(q)).X.|.error("|
|000014f0| 54 61 6b 69 6e 67 20 62 | 61 64 20 72 6f 6f 74 20 |Taking b|ad root |
|00001500| 6f 66 20 63 6f 6d 70 6c | 65 78 20 6e 75 6d 62 65 |of compl|ex numbe|
|00001510| 72 22 29 3b 0a 58 09 69 | 66 20 28 63 69 73 6f 6e |r");.X.i|f (cison|
|00001520| 65 28 63 29 20 7c 7c 20 | 71 69 73 6f 6e 65 28 71 |e(c) || |qisone(q|
|00001530| 29 29 0a 58 09 09 72 65 | 74 75 72 6e 20 63 6c 69 |)).X..re|turn cli|
|00001540| 6e 6b 28 63 29 3b 0a 58 | 09 69 66 20 28 71 69 73 |nk(c);.X|.if (qis|
|00001550| 74 77 6f 28 71 29 29 0a | 58 09 09 72 65 74 75 72 |two(q)).|X..retur|
|00001560| 6e 20 63 73 71 72 74 28 | 63 2c 20 65 70 73 69 6c |n csqrt(|c, epsil|
|00001570| 6f 6e 29 3b 0a 58 09 72 | 20 3d 20 63 6f 6d 61 6c |on);.X.r| = comal|
|00001580| 6c 6f 63 28 29 3b 0a 58 | 09 69 66 20 28 63 69 73 |loc();.X|.if (cis|
|00001590| 72 65 61 6c 28 63 29 20 | 26 26 20 21 71 69 73 6e |real(c) |&& !qisn|
|000015a0| 65 67 28 63 2d 3e 72 65 | 61 6c 29 29 20 7b 0a 58 |eg(c->re|al)) {.X|
|000015b0| 09 09 72 2d 3e 72 65 61 | 6c 20 3d 20 71 72 6f 6f |..r->rea|l = qroo|
|000015c0| 74 28 63 2d 3e 72 65 61 | 6c 2c 20 71 2c 20 65 70 |t(c->rea|l, q, ep|
|000015d0| 73 69 6c 6f 6e 29 3b 0a | 58 09 09 72 65 74 75 72 |silon);.|X..retur|
|000015e0| 6e 20 72 3b 0a 58 09 7d | 0a 58 09 2f 2a 0a 58 09 |n r;.X.}|.X./*.X.|
|000015f0| 20 2a 20 43 61 6c 63 75 | 6c 61 74 65 20 74 68 65 | * Calcu|late the|
|00001600| 20 72 6f 6f 74 20 75 73 | 69 6e 67 20 74 68 65 20 | root us|ing the |
|00001610| 66 6f 72 6d 75 6c 61 3a | 0a 58 09 20 2a 09 63 72 |formula:|.X. *.cr|
|00001620| 6f 6f 74 28 61 20 2b 20 | 62 69 2c 20 6e 29 20 3d |oot(a + |bi, n) =|
|00001630| 0a 58 09 20 2a 09 09 63 | 70 6f 6c 61 72 28 71 72 |.X. *..c|polar(qr|
|00001640| 6f 6f 74 28 61 5e 32 20 | 2b 20 62 5e 32 2c 20 32 |oot(a^2 |+ b^2, 2|
|00001650| 20 2a 20 6e 29 2c 20 71 | 61 74 61 6e 32 28 62 2c | * n), q|atan2(b,|
|00001660| 20 61 29 20 2f 20 6e 29 | 2e 0a 58 09 20 2a 2f 0a | a) / n)|..X. */.|
|00001670| 58 09 65 70 73 69 6c 6f | 6e 32 20 3d 20 71 73 63 |X.epsilo|n2 = qsc|
|00001680| 61 6c 65 28 65 70 73 69 | 6c 6f 6e 2c 20 2d 38 4c |ale(epsi|lon, -8L|
|00001690| 29 3b 0a 58 09 74 6d 70 | 31 20 3d 20 71 73 71 75 |);.X.tmp|1 = qsqu|
|000016a0| 61 72 65 28 63 2d 3e 72 | 65 61 6c 29 3b 0a 58 09 |are(c->r|eal);.X.|
|000016b0| 74 6d 70 32 20 3d 20 71 | 73 71 75 61 72 65 28 63 |tmp2 = q|square(c|
|000016c0| 2d 3e 69 6d 61 67 29 3b | 0a 58 09 61 32 70 62 32 |->imag);|.X.a2pb2|
|000016d0| 20 3d 20 71 61 64 64 28 | 74 6d 70 31 2c 20 74 6d | = qadd(|tmp1, tm|
|000016e0| 70 32 29 3b 0a 58 09 71 | 66 72 65 65 28 74 6d 70 |p2);.X.q|free(tmp|
|000016f0| 31 29 3b 0a 58 09 71 66 | 72 65 65 28 74 6d 70 32 |1);.X.qf|ree(tmp2|
|00001700| 29 3b 0a 58 09 74 6d 70 | 31 20 3d 20 71 73 63 61 |);.X.tmp|1 = qsca|
|00001710| 6c 65 28 71 2c 20 31 4c | 29 3b 0a 58 09 72 6f 6f |le(q, 1L|);.X.roo|
|00001720| 74 20 3d 20 71 72 6f 6f | 74 28 61 32 70 62 32 2c |t = qroo|t(a2pb2,|
|00001730| 20 74 6d 70 31 2c 20 65 | 70 73 69 6c 6f 6e 32 29 | tmp1, e|psilon2)|
|00001740| 3b 0a 58 09 71 66 72 65 | 65 28 61 32 70 62 32 29 |;.X.qfre|e(a2pb2)|
|00001750| 3b 0a 58 09 71 66 72 65 | 65 28 74 6d 70 31 29 3b |;.X.qfre|e(tmp1);|
|00001760| 0a 58 09 74 6d 70 31 20 | 3d 20 71 61 74 61 6e 32 |.X.tmp1 |= qatan2|
|00001770| 28 63 2d 3e 69 6d 61 67 | 2c 20 63 2d 3e 72 65 61 |(c->imag|, c->rea|
|00001780| 6c 2c 20 65 70 73 69 6c | 6f 6e 32 29 3b 0a 58 09 |l, epsil|on2);.X.|
|00001790| 71 66 72 65 65 28 65 70 | 73 69 6c 6f 6e 32 29 3b |qfree(ep|silon2);|
|000017a0| 0a 58 09 74 6d 70 32 20 | 3d 20 71 64 69 76 28 74 |.X.tmp2 |= qdiv(t|
|000017b0| 6d 70 31 2c 20 71 29 3b | 0a 58 09 71 66 72 65 65 |mp1, q);|.X.qfree|
|000017c0| 28 74 6d 70 31 29 3b 0a | 58 09 72 20 3d 20 63 70 |(tmp1);.|X.r = cp|
|000017d0| 6f 6c 61 72 28 72 6f 6f | 74 2c 20 74 6d 70 32 2c |olar(roo|t, tmp2,|
|000017e0| 20 65 70 73 69 6c 6f 6e | 29 3b 0a 58 09 71 66 72 | epsilon|);.X.qfr|
|000017f0| 65 65 28 72 6f 6f 74 29 | 3b 0a 58 09 71 66 72 65 |ee(root)|;.X.qfre|
|00001800| 65 28 74 6d 70 32 29 3b | 0a 58 09 72 65 74 75 72 |e(tmp2);|.X.retur|
|00001810| 6e 20 72 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 0a |n r;.X}.|X.X.X/*.|
|00001820| 58 20 2a 20 43 61 6c 63 | 75 6c 61 74 65 20 74 68 |X * Calc|ulate th|
|00001830| 65 20 63 6f 6d 70 6c 65 | 78 20 65 78 70 6f 6e 65 |e comple|x expone|
|00001840| 6e 74 69 61 6c 20 66 75 | 6e 63 74 69 6f 6e 20 74 |ntial fu|nction t|
|00001850| 6f 20 74 68 65 20 64 65 | 73 69 72 65 64 20 61 63 |o the de|sired ac|
|00001860| 63 75 72 61 63 79 2e 0a | 58 20 2a 20 57 65 20 75 |curacy..|X * We u|
|00001870| 73 65 20 74 68 65 20 66 | 6f 72 6d 75 6c 61 3a 0a |se the f|ormula:.|
|00001880| 58 20 2a 09 65 78 70 28 | 61 20 2b 20 62 69 29 20 |X *.exp(|a + bi) |
|00001890| 3d 20 65 78 70 28 61 29 | 20 2a 20 28 63 6f 73 28 |= exp(a)| * (cos(|
|000018a0| 62 29 20 2b 20 69 20 2a | 20 73 69 6e 28 62 29 29 |b) + i *| sin(b))|
|000018b0| 2e 0a 58 20 2a 2f 0a 58 | 43 4f 4d 50 4c 45 58 20 |..X */.X|COMPLEX |
|000018c0| 2a 0a 58 63 65 78 70 28 | 63 2c 20 65 70 73 69 6c |*.Xcexp(|c, epsil|
|000018d0| 6f 6e 29 0a 58 09 43 4f | 4d 50 4c 45 58 20 2a 63 |on).X.CO|MPLEX *c|
|000018e0| 3b 0a 58 09 4e 55 4d 42 | 45 52 20 2a 65 70 73 69 |;.X.NUMB|ER *epsi|
|000018f0| 6c 6f 6e 3b 0a 58 7b 0a | 58 09 43 4f 4d 50 4c 45 |lon;.X{.|X.COMPLE|
|00001900| 58 20 2a 72 3b 0a 58 09 | 4e 55 4d 42 45 52 20 2a |X *r;.X.|NUMBER *|
|00001910| 74 6d 70 31 2c 20 2a 74 | 6d 70 32 2c 20 2a 65 70 |tmp1, *t|mp2, *ep|
|00001920| 73 69 6c 6f 6e 32 3b 0a | 58 0a 58 09 69 66 20 28 |silon2;.|X.X.if (|
|00001930| 63 69 73 7a 65 72 6f 28 | 63 29 29 0a 58 09 09 72 |ciszero(|c)).X..r|
|00001940| 65 74 75 72 6e 20 63 6c | 69 6e 6b 28 26 5f 63 6f |eturn cl|ink(&_co|
|00001950| 6e 65 5f 29 3b 0a 58 09 | 72 20 3d 20 63 6f 6d 61 |ne_);.X.|r = coma|
|00001960| 6c 6c 6f 63 28 29 3b 0a | 58 09 69 66 20 28 63 69 |lloc();.|X.if (ci|
|00001970| 73 72 65 61 6c 28 63 29 | 29 20 7b 0a 58 09 09 72 |sreal(c)|) {.X..r|
|00001980| 2d 3e 72 65 61 6c 20 3d | 20 71 65 78 70 28 63 2d |->real =| qexp(c-|
|00001990| 3e 72 65 61 6c 2c 20 65 | 70 73 69 6c 6f 6e 29 3b |>real, e|psilon);|
|000019a0| 0a 58 09 09 72 65 74 75 | 72 6e 20 72 3b 0a 58 09 |.X..retu|rn r;.X.|
|000019b0| 7d 0a 58 09 65 70 73 69 | 6c 6f 6e 32 20 3d 20 71 |}.X.epsi|lon2 = q|
|000019c0| 73 63 61 6c 65 28 65 70 | 73 69 6c 6f 6e 2c 20 2d |scale(ep|silon, -|
|000019d0| 32 4c 29 3b 0a 58 09 72 | 2d 3e 72 65 61 6c 20 3d |2L);.X.r|->real =|
|000019e0| 20 71 63 6f 73 28 63 2d | 3e 69 6d 61 67 2c 20 65 | qcos(c-|>imag, e|
|000019f0| 70 73 69 6c 6f 6e 32 29 | 3b 0a 58 09 72 2d 3e 69 |psilon2)|;.X.r->i|
|00001a00| 6d 61 67 20 3d 20 71 6c | 65 67 74 6f 6c 65 67 28 |mag = ql|egtoleg(|
|00001a10| 72 2d 3e 72 65 61 6c 2c | 20 65 70 73 69 6c 6f 6e |r->real,| epsilon|
|00001a20| 32 2c 20 5f 73 69 6e 69 | 73 6e 65 67 5f 29 3b 0a |2, _sini|sneg_);.|
|00001a30| 58 09 69 66 20 28 71 69 | 73 7a 65 72 6f 28 63 2d |X.if (qi|szero(c-|
|00001a40| 3e 72 65 61 6c 29 29 20 | 7b 0a 58 09 09 71 66 72 |>real)) |{.X..qfr|
|00001a50| 65 65 28 65 70 73 69 6c | 6f 6e 32 29 3b 0a 58 09 |ee(epsil|on2);.X.|
|00001a60| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 09 7d 0a 58 |.return |r;.X.}.X|
|00001a70| 09 74 6d 70 31 20 3d 20 | 71 65 78 70 28 63 2d 3e |.tmp1 = |qexp(c->|
|00001a80| 72 65 61 6c 2c 20 65 70 | 73 69 6c 6f 6e 32 29 3b |real, ep|silon2);|
|00001a90| 0a 58 09 71 66 72 65 65 | 28 65 70 73 69 6c 6f 6e |.X.qfree|(epsilon|
|00001aa0| 32 29 3b 0a 58 09 74 6d | 70 32 20 3d 20 71 6d 75 |2);.X.tm|p2 = qmu|
|00001ab0| 6c 28 72 2d 3e 72 65 61 | 6c 2c 20 74 6d 70 31 29 |l(r->rea|l, tmp1)|
|00001ac0| 3b 0a 58 09 71 66 72 65 | 65 28 72 2d 3e 72 65 61 |;.X.qfre|e(r->rea|
|00001ad0| 6c 29 3b 0a 58 09 72 2d | 3e 72 65 61 6c 20 3d 20 |l);.X.r-|>real = |
|00001ae0| 74 6d 70 32 3b 0a 58 09 | 74 6d 70 32 20 3d 20 71 |tmp2;.X.|tmp2 = q|
|00001af0| 6d 75 6c 28 72 2d 3e 69 | 6d 61 67 2c 20 74 6d 70 |mul(r->i|mag, tmp|
|00001b00| 31 29 3b 0a 58 09 71 66 | 72 65 65 28 72 2d 3e 69 |1);.X.qf|ree(r->i|
|00001b10| 6d 61 67 29 3b 0a 58 09 | 71 66 72 65 65 28 74 6d |mag);.X.|qfree(tm|
|00001b20| 70 31 29 3b 0a 58 09 72 | 2d 3e 69 6d 61 67 20 3d |p1);.X.r|->imag =|
|00001b30| 20 74 6d 70 32 3b 0a 58 | 09 72 65 74 75 72 6e 20 | tmp2;.X|.return |
|00001b40| 72 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |r;.X}.X.|X.X/*.X |
|00001b50| 2a 20 43 61 6c 63 75 6c | 61 74 65 20 74 68 65 20 |* Calcul|ate the |
|00001b60| 6e 61 74 75 72 61 6c 20 | 6c 6f 67 61 72 69 74 68 |natural |logarith|
|00001b70| 6d 20 6f 66 20 61 20 63 | 6f 6d 70 6c 65 78 20 6e |m of a c|omplex n|
|00001b80| 75 6d 62 65 72 20 77 69 | 74 68 69 6e 20 74 68 65 |umber wi|thin the|
|00001b90| 20 73 70 65 63 69 66 69 | 65 64 0a 58 20 2a 20 65 | specifi|ed.X * e|
|00001ba0| 72 72 6f 72 2e 20 20 57 | 65 20 75 73 65 20 74 68 |rror.  W|e use th|
|00001bb0| 65 20 66 6f 72 6d 75 6c | 61 3a 0a 58 20 2a 09 6c |e formul|a:.X *.l|
|00001bc0| 6e 28 61 20 2b 20 62 69 | 29 20 3d 20 6c 6e 28 61 |n(a + bi|) = ln(a|
|00001bd0| 5e 32 20 2b 20 62 5e 32 | 29 20 2f 20 32 20 2b 20 |^2 + b^2|) / 2 + |
|00001be0| 69 20 2a 20 61 74 61 6e | 32 28 62 2c 20 61 29 2e |i * atan|2(b, a).|
|00001bf0| 0a 58 20 2a 2f 0a 58 43 | 4f 4d 50 4c 45 58 20 2a |.X */.XC|OMPLEX *|
|00001c00| 0a 58 63 6c 6e 28 63 2c | 20 65 70 73 69 6c 6f 6e |.Xcln(c,| epsilon|
|00001c10| 29 0a 58 09 43 4f 4d 50 | 4c 45 58 20 2a 63 3b 0a |).X.COMP|LEX *c;.|
|00001c20| 58 09 4e 55 4d 42 45 52 | 20 2a 65 70 73 69 6c 6f |X.NUMBER| *epsilo|
|00001c30| 6e 3b 0a 58 7b 0a 58 09 | 43 4f 4d 50 4c 45 58 20 |n;.X{.X.|COMPLEX |
|00001c40| 2a 72 3b 0a 58 09 4e 55 | 4d 42 45 52 20 2a 61 32 |*r;.X.NU|MBER *a2|
|00001c50| 62 32 2c 20 2a 74 6d 70 | 31 2c 20 2a 74 6d 70 32 |b2, *tmp|1, *tmp2|
|00001c60| 3b 0a 58 0a 58 09 69 66 | 20 28 63 69 73 7a 65 72 |;.X.X.if| (ciszer|
|00001c70| 6f 28 63 29 29 0a 58 09 | 09 65 72 72 6f 72 28 22 |o(c)).X.|.error("|
|00001c80| 4c 6f 67 61 72 69 74 68 | 6d 20 6f 66 20 7a 65 72 |Logarith|m of zer|
|00001c90| 6f 22 29 3b 0a 58 09 69 | 66 20 28 63 69 73 6f 6e |o");.X.i|f (cison|
|00001ca0| 65 28 63 29 29 0a 58 09 | 09 72 65 74 75 72 6e 20 |e(c)).X.|.return |
|00001cb0| 63 6c 69 6e 6b 28 26 5f | 63 7a 65 72 6f 5f 29 3b |clink(&_|czero_);|
|00001cc0| 0a 58 09 72 20 3d 20 63 | 6f 6d 61 6c 6c 6f 63 28 |.X.r = c|omalloc(|
|00001cd0| 29 3b 0a 58 09 69 66 20 | 28 63 69 73 72 65 61 6c |);.X.if |(cisreal|
|00001ce0| 28 63 29 20 26 26 20 21 | 71 69 73 6e 65 67 28 63 |(c) && !|qisneg(c|
|00001cf0| 2d 3e 72 65 61 6c 29 29 | 20 7b 0a 58 09 09 72 2d |->real))| {.X..r-|
|00001d00| 3e 72 65 61 6c 20 3d 20 | 71 6c 6e 28 63 2d 3e 72 |>real = |qln(c->r|
|00001d10| 65 61 6c 2c 20 65 70 73 | 69 6c 6f 6e 29 3b 0a 58 |eal, eps|ilon);.X|
|00001d20| 09 09 72 65 74 75 72 6e | 20 72 3b 0a 58 09 7d 0a |..return| r;.X.}.|
|00001d30| 58 09 74 6d 70 31 20 3d | 20 71 73 71 75 61 72 65 |X.tmp1 =| qsquare|
|00001d40| 28 63 2d 3e 72 65 61 6c | 29 3b 0a 58 09 74 6d 70 |(c->real|);.X.tmp|
|00001d50| 32 20 3d 20 71 73 71 75 | 61 72 65 28 63 2d 3e 69 |2 = qsqu|are(c->i|
|00001d60| 6d 61 67 29 3b 0a 58 09 | 61 32 62 32 20 3d 20 71 |mag);.X.|a2b2 = q|
|00001d70| 61 64 64 28 74 6d 70 31 | 2c 20 74 6d 70 32 29 3b |add(tmp1|, tmp2);|
|00001d80| 0a 58 09 71 66 72 65 65 | 28 74 6d 70 31 29 3b 0a |.X.qfree|(tmp1);.|
|00001d90| 58 09 71 66 72 65 65 28 | 74 6d 70 32 29 3b 0a 58 |X.qfree(|tmp2);.X|
|00001da0| 09 74 6d 70 31 20 3d 20 | 71 6c 6e 28 61 32 62 32 |.tmp1 = |qln(a2b2|
|00001db0| 2c 20 65 70 73 69 6c 6f | 6e 29 3b 0a 58 09 71 66 |, epsilo|n);.X.qf|
|00001dc0| 72 65 65 28 61 32 62 32 | 29 3b 0a 58 09 72 2d 3e |ree(a2b2|);.X.r->|
|00001dd0| 72 65 61 6c 20 3d 20 71 | 73 63 61 6c 65 28 74 6d |real = q|scale(tm|
|00001de0| 70 31 2c 20 2d 31 4c 29 | 3b 0a 58 09 71 66 72 65 |p1, -1L)|;.X.qfre|
|00001df0| 65 28 74 6d 70 31 29 3b | 0a 58 09 72 2d 3e 69 6d |e(tmp1);|.X.r->im|
|00001e00| 61 67 20 3d 20 71 61 74 | 61 6e 32 28 63 2d 3e 69 |ag = qat|an2(c->i|
|00001e10| 6d 61 67 2c 20 63 2d 3e | 72 65 61 6c 2c 20 65 70 |mag, c->|real, ep|
|00001e20| 73 69 6c 6f 6e 29 3b 0a | 58 09 72 65 74 75 72 6e |silon);.|X.return|
|00001e30| 20 72 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 0a 58 | r;.X}.X|.X.X/*.X|
|00001e40| 20 2a 20 43 61 6c 63 75 | 6c 61 74 65 20 74 68 65 | * Calcu|late the|
|00001e50| 20 63 6f 6d 70 6c 65 78 | 20 63 6f 73 69 6e 65 20 | complex| cosine |
|00001e60| 77 69 74 68 69 6e 20 74 | 68 65 20 73 70 65 63 69 |within t|he speci|
|00001e70| 66 69 65 64 20 61 63 63 | 75 72 61 63 79 2e 0a 58 |fied acc|uracy..X|
|00001e80| 20 2a 20 54 68 69 73 20 | 75 73 65 73 20 74 68 65 | * This |uses the|
|00001e90| 20 66 6f 72 6d 75 6c 61 | 3a 0a 58 20 2a 09 63 6f | formula|:.X *.co|
|00001ea0| 73 28 61 20 2b 20 62 69 | 29 20 3d 20 63 6f 73 28 |s(a + bi|) = cos(|
|00001eb0| 61 29 20 2a 20 63 6f 73 | 68 28 62 29 20 2d 20 73 |a) * cos|h(b) - s|
|00001ec0| 69 6e 28 61 29 20 2a 20 | 73 69 6e 68 28 62 29 20 |in(a) * |sinh(b) |
|00001ed0| 2a 20 69 2e 0a 58 20 2a | 2f 0a 58 43 4f 4d 50 4c |* i..X *|/.XCOMPL|
|00001ee0| 45 58 20 2a 0a 58 63 63 | 6f 73 28 63 2c 20 65 70 |EX *.Xcc|os(c, ep|
|00001ef0| 73 69 6c 6f 6e 29 0a 58 | 09 43 4f 4d 50 4c 45 58 |silon).X|.COMPLEX|
|00001f00| 20 2a 63 3b 0a 58 09 4e | 55 4d 42 45 52 20 2a 65 | *c;.X.N|UMBER *e|
|00001f10| 70 73 69 6c 6f 6e 3b 0a | 58 7b 0a 58 09 43 4f 4d |psilon;.|X{.X.COM|
|00001f20| 50 4c 45 58 20 2a 72 3b | 0a 58 09 4e 55 4d 42 45 |PLEX *r;|.X.NUMBE|
|00001f30| 52 20 2a 63 6f 73 76 61 | 6c 2c 20 2a 63 6f 73 68 |R *cosva|l, *cosh|
|00001f40| 76 61 6c 2c 20 2a 74 6d | 70 31 2c 20 2a 74 6d 70 |val, *tm|p1, *tmp|
|00001f50| 32 2c 20 2a 74 6d 70 33 | 2c 20 2a 65 70 73 69 6c |2, *tmp3|, *epsil|
|00001f60| 6f 6e 32 3b 0a 58 09 69 | 6e 74 20 6e 65 67 69 6d |on2;.X.i|nt negim|
|00001f70| 61 67 3b 0a 58 0a 58 09 | 69 66 20 28 63 69 73 7a |ag;.X.X.|if (cisz|
|00001f80| 65 72 6f 28 63 29 29 0a | 58 09 09 72 65 74 75 72 |ero(c)).|X..retur|
|00001f90| 6e 20 63 6c 69 6e 6b 28 | 26 5f 63 6f 6e 65 5f 29 |n clink(|&_cone_)|
|00001fa0| 3b 0a 58 09 72 20 3d 20 | 63 6f 6d 61 6c 6c 6f 63 |;.X.r = |comalloc|
|00001fb0| 28 29 3b 0a 58 09 69 66 | 20 28 63 69 73 72 65 61 |();.X.if| (cisrea|
|00001fc0| 6c 28 63 29 29 20 7b 0a | 58 09 09 72 2d 3e 72 65 |l(c)) {.|X..r->re|
|00001fd0| 61 6c 20 3d 20 71 63 6f | 73 28 63 2d 3e 72 65 61 |al = qco|s(c->rea|
|00001fe0| 6c 2c 20 65 70 73 69 6c | 6f 6e 29 3b 0a 58 09 09 |l, epsil|on);.X..|
|00001ff0| 72 65 74 75 72 6e 20 72 | 3b 0a 58 09 7d 0a 58 09 |return r|;.X.}.X.|
|00002000| 69 66 20 28 71 69 73 7a | 65 72 6f 28 63 2d 3e 72 |if (qisz|ero(c->r|
|00002010| 65 61 6c 29 29 20 7b 0a | 58 09 09 72 2d 3e 72 65 |eal)) {.|X..r->re|
|00002020| 61 6c 20 3d 20 71 63 6f | 73 68 28 63 2d 3e 69 6d |al = qco|sh(c->im|
|00002030| 61 67 2c 20 65 70 73 69 | 6c 6f 6e 29 3b 0a 58 09 |ag, epsi|lon);.X.|
|00002040| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 09 7d 0a 58 |.return |r;.X.}.X|
|00002050| 09 65 70 73 69 6c 6f 6e | 32 20 3d 20 71 73 63 61 |.epsilon|2 = qsca|
|00002060| 6c 65 28 65 70 73 69 6c | 6f 6e 2c 20 2d 32 4c 29 |le(epsil|on, -2L)|
|00002070| 3b 0a 58 09 63 6f 73 68 | 76 61 6c 20 3d 20 71 63 |;.X.cosh|val = qc|
|00002080| 6f 73 68 28 63 2d 3e 69 | 6d 61 67 2c 20 65 70 73 |osh(c->i|mag, eps|
|00002090| 69 6c 6f 6e 32 29 3b 0a | 58 09 63 6f 73 76 61 6c |ilon2);.|X.cosval|
|000020a0| 20 3d 20 71 63 6f 73 28 | 63 2d 3e 72 65 61 6c 2c | = qcos(|c->real,|
|000020b0| 20 65 70 73 69 6c 6f 6e | 32 29 3b 0a 58 09 6e 65 | epsilon|2);.X.ne|
|000020c0| 67 69 6d 61 67 20 3d 20 | 21 5f 73 69 6e 69 73 6e |gimag = |!_sinisn|
|000020d0| 65 67 5f 3b 0a 58 09 69 | 66 20 28 71 69 73 6e 65 |eg_;.X.i|f (qisne|
|000020e0| 67 28 63 2d 3e 69 6d 61 | 67 29 29 0a 58 09 09 6e |g(c->ima|g)).X..n|
|000020f0| 65 67 69 6d 61 67 20 3d | 20 21 6e 65 67 69 6d 61 |egimag =| !negima|
|00002100| 67 3b 0a 58 09 72 2d 3e | 72 65 61 6c 20 3d 20 71 |g;.X.r->|real = q|
|00002110| 6d 75 6c 28 63 6f 73 76 | 61 6c 2c 20 63 6f 73 68 |mul(cosv|al, cosh|
|00002120| 76 61 6c 29 3b 0a 58 09 | 2f 2a 0a 58 09 20 2a 20 |val);.X.|/*.X. * |
|00002130| 43 61 6c 63 75 6c 61 74 | 65 20 74 68 65 20 69 6d |Calculat|e the im|
|00002140| 61 67 69 6e 61 72 79 20 | 70 61 72 74 20 75 73 69 |aginary |part usi|
|00002150| 6e 67 20 74 68 65 20 66 | 6f 72 6d 75 6c 61 3a 0a |ng the f|ormula:.|
|00002160| 58 09 20 2a 09 73 69 6e | 28 61 29 20 2a 20 73 69 |X. *.sin|(a) * si|
|00002170| 6e 68 28 62 29 20 3d 20 | 73 71 72 74 28 28 31 20 |nh(b) = |sqrt((1 |
|00002180| 2d 20 61 5e 32 29 20 2a | 20 28 62 5e 32 20 2d 20 |- a^2) *| (b^2 - |
|00002190| 31 29 29 2e 0a 58 09 20 | 2a 2f 0a 58 09 74 6d 70 |1))..X. |*/.X.tmp|
|000021a0| 31 20 3d 20 71 73 71 75 | 61 72 65 28 63 6f 73 76 |1 = qsqu|are(cosv|
|000021b0| 61 6c 29 3b 0a 58 09 71 | 66 72 65 65 28 63 6f 73 |al);.X.q|free(cos|
|000021c0| 76 61 6c 29 3b 0a 58 09 | 74 6d 70 32 20 3d 20 71 |val);.X.|tmp2 = q|
|000021d0| 64 65 63 28 74 6d 70 31 | 29 3b 0a 58 09 71 66 72 |dec(tmp1|);.X.qfr|
|000021e0| 65 65 28 74 6d 70 31 29 | 3b 0a 58 09 74 6d 70 31 |ee(tmp1)|;.X.tmp1|
|000021f0| 20 3d 20 71 6e 65 67 28 | 74 6d 70 32 29 3b 0a 58 | = qneg(|tmp2);.X|
|00002200| 09 71 66 72 65 65 28 74 | 6d 70 32 29 3b 0a 58 09 |.qfree(t|mp2);.X.|
|00002210| 74 6d 70 32 20 3d 20 71 | 73 71 75 61 72 65 28 63 |tmp2 = q|square(c|
|00002220| 6f 73 68 76 61 6c 29 3b | 0a 58 09 71 66 72 65 65 |oshval);|.X.qfree|
|00002230| 28 63 6f 73 68 76 61 6c | 29 3b 0a 58 09 74 6d 70 |(coshval|);.X.tmp|
|00002240| 33 20 3d 20 71 64 65 63 | 28 74 6d 70 32 29 3b 0a |3 = qdec|(tmp2);.|
|00002250| 58 09 71 66 72 65 65 28 | 74 6d 70 32 29 3b 0a 58 |X.qfree(|tmp2);.X|
|00002260| 09 74 6d 70 32 20 3d 20 | 71 6d 75 6c 28 74 6d 70 |.tmp2 = |qmul(tmp|
|00002270| 31 2c 20 74 6d 70 33 29 | 3b 0a 58 09 71 66 72 65 |1, tmp3)|;.X.qfre|
|00002280| 65 28 74 6d 70 31 29 3b | 0a 58 09 71 66 72 65 65 |e(tmp1);|.X.qfree|
|00002290| 28 74 6d 70 33 29 3b 0a | 58 09 72 2d 3e 69 6d 61 |(tmp3);.|X.r->ima|
|000022a0| 67 20 3d 20 71 73 71 72 | 74 28 74 6d 70 32 2c 20 |g = qsqr|t(tmp2, |
|000022b0| 65 70 73 69 6c 6f 6e 32 | 29 3b 0a 58 09 71 66 72 |epsilon2|);.X.qfr|
|000022c0| 65 65 28 74 6d 70 32 29 | 3b 0a 58 09 71 66 72 65 |ee(tmp2)|;.X.qfre|
|000022d0| 65 28 65 70 73 69 6c 6f | 6e 32 29 3b 0a 58 09 69 |e(epsilo|n2);.X.i|
|000022e0| 66 20 28 6e 65 67 69 6d | 61 67 29 20 7b 0a 58 09 |f (negim|ag) {.X.|
|000022f0| 09 74 6d 70 31 20 3d 20 | 71 6e 65 67 28 72 2d 3e |.tmp1 = |qneg(r->|
|00002300| 69 6d 61 67 29 3b 0a 58 | 09 09 71 66 72 65 65 28 |imag);.X|..qfree(|
|00002310| 72 2d 3e 69 6d 61 67 29 | 3b 0a 58 09 09 72 2d 3e |r->imag)|;.X..r->|
|00002320| 69 6d 61 67 20 3d 20 74 | 6d 70 31 3b 0a 58 09 7d |imag = t|mp1;.X.}|
|00002330| 0a 58 09 72 65 74 75 72 | 6e 20 72 3b 0a 58 7d 0a |.X.retur|n r;.X}.|
|00002340| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 43 61 6c 63 |X.X.X/*.|X * Calc|
|00002350| 75 6c 61 74 65 20 74 68 | 65 20 63 6f 6d 70 6c 65 |ulate th|e comple|
|00002360| 78 20 73 69 6e 65 20 77 | 69 74 68 69 6e 20 74 68 |x sine w|ithin th|
|00002370| 65 20 73 70 65 63 69 66 | 69 65 64 20 61 63 63 75 |e specif|ied accu|
|00002380| 72 61 63 79 2e 0a 58 20 | 2a 20 54 68 69 73 20 75 |racy..X |* This u|
|00002390| 73 65 73 20 74 68 65 20 | 66 6f 72 6d 75 6c 61 3a |ses the |formula:|
|000023a0| 0a 58 20 2a 09 73 69 6e | 28 61 20 2b 20 62 69 29 |.X *.sin|(a + bi)|
|000023b0| 20 3d 20 73 69 6e 28 61 | 29 20 2a 20 63 6f 73 68 | = sin(a|) * cosh|
|000023c0| 28 62 29 20 2b 20 63 6f | 73 28 61 29 20 2a 20 73 |(b) + co|s(a) * s|
|000023d0| 69 6e 68 28 62 29 20 2a | 20 69 2e 0a 58 20 2a 2f |inh(b) *| i..X */|
|000023e0| 0a 58 43 4f 4d 50 4c 45 | 58 20 2a 0a 58 63 73 69 |.XCOMPLE|X *.Xcsi|
|000023f0| 6e 28 63 2c 20 65 70 73 | 69 6c 6f 6e 29 0a 58 09 |n(c, eps|ilon).X.|
|00002400| 43 4f 4d 50 4c 45 58 20 | 2a 63 3b 0a 58 09 4e 55 |COMPLEX |*c;.X.NU|
|00002410| 4d 42 45 52 20 2a 65 70 | 73 69 6c 6f 6e 3b 0a 58 |MBER *ep|silon;.X|
|00002420| 7b 0a 58 09 43 4f 4d 50 | 4c 45 58 20 2a 72 3b 0a |{.X.COMP|LEX *r;.|
|00002430| 58 0a 58 09 4e 55 4d 42 | 45 52 20 2a 63 6f 73 76 |X.X.NUMB|ER *cosv|
|00002440| 61 6c 2c 20 2a 63 6f 73 | 68 76 61 6c 2c 20 2a 74 |al, *cos|hval, *t|
|00002450| 6d 70 31 2c 20 2a 74 6d | 70 32 2c 20 2a 65 70 73 |mp1, *tm|p2, *eps|
|00002460| 69 6c 6f 6e 32 3b 0a 58 | 0a 58 09 69 66 20 28 63 |ilon2;.X|.X.if (c|
|00002470| 69 73 7a 65 72 6f 28 63 | 29 29 0a 58 09 09 72 65 |iszero(c|)).X..re|
|00002480| 74 75 72 6e 20 63 6c 69 | 6e 6b 28 26 5f 63 7a 65 |turn cli|nk(&_cze|
|00002490| 72 6f 5f 29 3b 0a 58 09 | 72 20 3d 20 63 6f 6d 61 |ro_);.X.|r = coma|
|000024a0| 6c 6c 6f 63 28 29 3b 0a | 58 09 69 66 20 28 63 69 |lloc();.|X.if (ci|
|000024b0| 73 72 65 61 6c 28 63 29 | 29 20 7b 0a 58 09 09 72 |sreal(c)|) {.X..r|
|000024c0| 2d 3e 72 65 61 6c 20 3d | 20 71 73 69 6e 28 63 2d |->real =| qsin(c-|
|000024d0| 3e 72 65 61 6c 2c 20 65 | 70 73 69 6c 6f 6e 29 3b |>real, e|psilon);|
|000024e0| 0a 58 09 09 72 65 74 75 | 72 6e 20 72 3b 0a 58 09 |.X..retu|rn r;.X.|
|000024f0| 7d 0a 58 09 69 66 20 28 | 71 69 73 7a 65 72 6f 28 |}.X.if (|qiszero(|
|00002500| 63 2d 3e 72 65 61 6c 29 | 29 20 7b 0a 58 09 09 72 |c->real)|) {.X..r|
|00002510| 2d 3e 69 6d 61 67 20 3d | 20 71 73 69 6e 68 28 63 |->imag =| qsinh(c|
|00002520| 2d 3e 69 6d 61 67 2c 20 | 65 70 73 69 6c 6f 6e 29 |->imag, |epsilon)|
|00002530| 3b 0a 58 09 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |;.X..ret|urn r;.X|
|00002540| 09 7d 0a 58 09 65 70 73 | 69 6c 6f 6e 32 20 3d 20 |.}.X.eps|ilon2 = |
|00002550| 71 73 63 61 6c 65 28 65 | 70 73 69 6c 6f 6e 2c 20 |qscale(e|psilon, |
|00002560| 2d 32 4c 29 3b 0a 58 09 | 63 6f 73 68 76 61 6c 20 |-2L);.X.|coshval |
|00002570| 3d 20 71 63 6f 73 68 28 | 63 2d 3e 69 6d 61 67 2c |= qcosh(|c->imag,|
|00002580| 20 65 70 73 69 6c 6f 6e | 32 29 3b 0a 58 09 63 6f | epsilon|2);.X.co|
|00002590| 73 76 61 6c 20 3d 20 71 | 63 6f 73 28 63 2d 3e 72 |sval = q|cos(c->r|
|000025a0| 65 61 6c 2c 20 65 70 73 | 69 6c 6f 6e 32 29 3b 0a |eal, eps|ilon2);.|
|000025b0| 58 09 74 6d 70 31 20 3d | 20 71 6c 65 67 74 6f 6c |X.tmp1 =| qlegtol|
|000025c0| 65 67 28 63 6f 73 76 61 | 6c 2c 20 65 70 73 69 6c |eg(cosva|l, epsil|
|000025d0| 6f 6e 32 2c 20 5f 73 69 | 6e 69 73 6e 65 67 5f 29 |on2, _si|nisneg_)|
|000025e0| 3b 0a 58 09 72 2d 3e 72 | 65 61 6c 20 3d 20 71 6d |;.X.r->r|eal = qm|
|000025f0| 75 6c 28 74 6d 70 31 2c | 20 63 6f 73 68 76 61 6c |ul(tmp1,| coshval|
|00002600| 29 3b 0a 58 09 71 66 72 | 65 65 28 74 6d 70 31 29 |);.X.qfr|ee(tmp1)|
|00002610| 3b 0a 58 09 74 6d 70 31 | 20 3d 20 71 73 71 75 61 |;.X.tmp1| = qsqua|
|00002620| 72 65 28 63 6f 73 68 76 | 61 6c 29 3b 0a 58 09 71 |re(coshv|al);.X.q|
|00002630| 66 72 65 65 28 63 6f 73 | 68 76 61 6c 29 3b 0a 58 |free(cos|hval);.X|
|00002640| 09 74 6d 70 32 20 3d 20 | 71 64 65 63 28 74 6d 70 |.tmp2 = |qdec(tmp|
|00002650| 31 29 3b 0a 58 09 71 66 | 72 65 65 28 74 6d 70 31 |1);.X.qf|ree(tmp1|
|00002660| 29 3b 0a 58 09 74 6d 70 | 31 20 3d 20 71 73 71 72 |);.X.tmp|1 = qsqr|
|00002670| 74 28 74 6d 70 32 2c 20 | 65 70 73 69 6c 6f 6e 32 |t(tmp2, |epsilon2|
|00002680| 29 3b 0a 58 09 71 66 72 | 65 65 28 74 6d 70 32 29 |);.X.qfr|ee(tmp2)|
|00002690| 3b 0a 58 09 72 2d 3e 69 | 6d 61 67 20 3d 20 71 6d |;.X.r->i|mag = qm|
|000026a0| 75 6c 28 74 6d 70 31 2c | 20 63 6f 73 76 61 6c 29 |ul(tmp1,| cosval)|
|000026b0| 3b 0a 58 09 71 66 72 65 | 65 28 74 6d 70 31 29 3b |;.X.qfre|e(tmp1);|
|000026c0| 0a 58 09 71 66 72 65 65 | 28 63 6f 73 76 61 6c 29 |.X.qfree|(cosval)|
|000026d0| 3b 0a 58 09 69 66 20 28 | 71 69 73 6e 65 67 28 63 |;.X.if (|qisneg(c|
|000026e0| 2d 3e 69 6d 61 67 29 29 | 20 7b 0a 58 09 09 74 6d |->imag))| {.X..tm|
|000026f0| 70 31 20 3d 20 71 6e 65 | 67 28 72 2d 3e 69 6d 61 |p1 = qne|g(r->ima|
|00002700| 67 29 3b 0a 58 09 09 71 | 66 72 65 65 28 72 2d 3e |g);.X..q|free(r->|
|00002710| 69 6d 61 67 29 3b 0a 58 | 09 09 72 2d 3e 69 6d 61 |imag);.X|..r->ima|
|00002720| 67 20 3d 20 74 6d 70 31 | 3b 0a 58 09 7d 0a 58 09 |g = tmp1|;.X.}.X.|
|00002730| 72 65 74 75 72 6e 20 72 | 3b 0a 58 7d 0a 58 0a 58 |return r|;.X}.X.X|
|00002740| 0a 58 2f 2a 0a 58 20 2a | 20 43 6f 6e 76 65 72 74 |.X/*.X *| Convert|
|00002750| 20 61 20 6e 75 6d 62 65 | 72 20 66 72 6f 6d 20 70 | a numbe|r from p|
|00002760| 6f 6c 61 72 20 63 6f 6f | 72 64 69 6e 61 74 65 73 |olar coo|rdinates|
|00002770| 20 74 6f 20 6e 6f 72 6d | 61 6c 20 63 6f 6d 70 6c | to norm|al compl|
|00002780| 65 78 20 6e 75 6d 62 65 | 72 20 66 6f 72 6d 0a 58 |ex numbe|r form.X|
|00002790| 20 2a 20 77 69 74 68 69 | 6e 20 74 68 65 20 73 70 | * withi|n the sp|
|000027a0| 65 63 69 66 69 65 64 20 | 61 63 63 75 72 61 63 79 |ecified |accuracy|
|000027b0| 2e 20 20 54 68 69 73 20 | 70 72 6f 64 75 63 65 73 |.  This |produces|
|000027c0| 20 74 68 65 20 76 61 6c | 75 65 3a 0a 58 20 2a 09 | the val|ue:.X *.|
|000027d0| 71 31 20 2a 20 63 6f 73 | 28 71 32 29 20 2b 20 71 |q1 * cos|(q2) + q|
|000027e0| 31 20 2a 20 73 69 6e 28 | 71 32 29 20 2a 20 69 2e |1 * sin(|q2) * i.|
|000027f0| 0a 58 20 2a 2f 0a 58 43 | 4f 4d 50 4c 45 58 20 2a |.X */.XC|OMPLEX *|
|00002800| 0a 58 63 70 6f 6c 61 72 | 28 71 31 2c 20 71 32 2c |.Xcpolar|(q1, q2,|
|00002810| 20 65 70 73 69 6c 6f 6e | 29 0a 58 09 4e 55 4d 42 | epsilon|).X.NUMB|
|00002820| 45 52 20 2a 71 31 2c 20 | 2a 71 32 2c 20 2a 65 70 |ER *q1, |*q2, *ep|
|00002830| 73 69 6c 6f 6e 3b 0a 58 | 7b 0a 58 09 43 4f 4d 50 |silon;.X|{.X.COMP|
|00002840| 4c 45 58 20 2a 72 3b 0a | 58 09 4e 55 4d 42 45 52 |LEX *r;.|X.NUMBER|
|00002850| 20 2a 74 6d 70 2c 20 2a | 65 70 73 69 6c 6f 6e 32 | *tmp, *|epsilon2|
|00002860| 3b 0a 58 09 6c 6f 6e 67 | 20 73 63 61 6c 65 3b 0a |;.X.long| scale;.|
|00002870| 58 0a 58 09 72 20 3d 20 | 63 6f 6d 61 6c 6c 6f 63 |X.X.r = |comalloc|
|00002880| 28 29 3b 0a 58 09 69 66 | 20 28 71 69 73 7a 65 72 |();.X.if| (qiszer|
|00002890| 6f 28 71 31 29 20 7c 7c | 20 71 69 73 7a 65 72 6f |o(q1) ||| qiszero|
|000028a0| 28 71 32 29 29 20 7b 0a | 58 09 09 72 2d 3e 72 65 |(q2)) {.|X..r->re|
|000028b0| 61 6c 20 3d 20 71 6c 69 | 6e 6b 28 71 31 29 3b 0a |al = qli|nk(q1);.|
|000028c0| 58 09 09 72 65 74 75 72 | 6e 20 72 3b 0a 58 09 7d |X..retur|n r;.X.}|
|000028d0| 0a 58 09 65 70 73 69 6c | 6f 6e 32 20 3d 20 65 70 |.X.epsil|on2 = ep|
|000028e0| 73 69 6c 6f 6e 3b 0a 58 | 09 69 66 20 28 21 71 69 |silon;.X|.if (!qi|
|000028f0| 73 75 6e 69 74 28 71 31 | 29 29 20 7b 0a 58 09 09 |sunit(q1|)) {.X..|
|00002900| 73 63 61 6c 65 20 3d 20 | 7a 68 69 67 68 62 69 74 |scale = |zhighbit|
|00002910| 28 71 31 2d 3e 6e 75 6d | 29 20 2d 20 7a 68 69 67 |(q1->num|) - zhig|
|00002920| 68 62 69 74 28 71 31 2d | 3e 64 65 6e 29 20 2b 20 |hbit(q1-|>den) + |
|00002930| 31 3b 0a 58 09 09 69 66 | 20 28 73 63 61 6c 65 20 |1;.X..if| (scale |
|00002940| 3e 20 30 29 0a 58 09 09 | 09 65 70 73 69 6c 6f 6e |> 0).X..|.epsilon|
|00002950| 32 20 3d 20 71 73 63 61 | 6c 65 28 65 70 73 69 6c |2 = qsca|le(epsil|
|00002960| 6f 6e 2c 20 2d 73 63 61 | 6c 65 29 3b 0a 58 09 7d |on, -sca|le);.X.}|
|00002970| 0a 58 09 72 2d 3e 72 65 | 61 6c 20 3d 20 71 63 6f |.X.r->re|al = qco|
|00002980| 73 28 71 32 2c 20 65 70 | 73 69 6c 6f 6e 32 29 3b |s(q2, ep|silon2);|
|00002990| 0a 58 09 72 2d 3e 69 6d | 61 67 20 3d 20 71 6c 65 |.X.r->im|ag = qle|
|000029a0| 67 74 6f 6c 65 67 28 72 | 2d 3e 72 65 61 6c 2c 20 |gtoleg(r|->real, |
|000029b0| 65 70 73 69 6c 6f 6e 32 | 2c 20 5f 73 69 6e 69 73 |epsilon2|, _sinis|
|000029c0| 6e 65 67 5f 29 3b 0a 58 | 09 69 66 20 28 65 70 73 |neg_);.X|.if (eps|
|000029d0| 69 6c 6f 6e 32 20 21 3d | 20 65 70 73 69 6c 6f 6e |ilon2 !=| epsilon|
|000029e0| 29 0a 58 09 09 71 66 72 | 65 65 28 65 70 73 69 6c |).X..qfr|ee(epsil|
|000029f0| 6f 6e 32 29 3b 0a 58 09 | 69 66 20 28 71 69 73 6f |on2);.X.|if (qiso|
|00002a00| 6e 65 28 71 31 29 29 0a | 58 09 09 72 65 74 75 72 |ne(q1)).|X..retur|
|00002a10| 6e 20 72 3b 0a 58 09 74 | 6d 70 20 3d 20 71 6d 75 |n r;.X.t|mp = qmu|
|00002a20| 6c 28 72 2d 3e 72 65 61 | 6c 2c 20 71 31 29 3b 0a |l(r->rea|l, q1);.|
|00002a30| 58 09 71 66 72 65 65 28 | 72 2d 3e 72 65 61 6c 29 |X.qfree(|r->real)|
|00002a40| 3b 0a 58 09 72 2d 3e 72 | 65 61 6c 20 3d 20 74 6d |;.X.r->r|eal = tm|
|00002a50| 70 3b 0a 58 09 74 6d 70 | 20 3d 20 71 6d 75 6c 28 |p;.X.tmp| = qmul(|
|00002a60| 72 2d 3e 69 6d 61 67 2c | 20 71 31 29 3b 0a 58 09 |r->imag,| q1);.X.|
|00002a70| 71 66 72 65 65 28 72 2d | 3e 69 6d 61 67 29 3b 0a |qfree(r-|>imag);.|
|00002a80| 58 09 72 2d 3e 69 6d 61 | 67 20 3d 20 74 6d 70 3b |X.r->ima|g = tmp;|
|00002a90| 0a 58 09 72 65 74 75 72 | 6e 20 72 3b 0a 58 7d 0a |.X.retur|n r;.X}.|
|00002aa0| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 52 61 69 73 |X.X.X/*.|X * Rais|
|00002ab0| 65 20 6f 6e 65 20 63 6f | 6d 70 6c 65 78 20 6e 75 |e one co|mplex nu|
|00002ac0| 6d 62 65 72 20 74 6f 20 | 74 68 65 20 70 6f 77 65 |mber to |the powe|
|00002ad0| 72 20 6f 66 20 61 6e 6f | 74 68 65 72 20 6f 6e 65 |r of ano|ther one|
|00002ae0| 20 74 6f 20 77 69 74 68 | 69 6e 20 74 68 65 0a 58 | to with|in the.X|
|00002af0| 20 2a 20 73 70 65 63 69 | 66 69 65 64 20 65 72 72 | * speci|fied err|
|00002b00| 6f 72 2e 0a 58 20 2a 2f | 0a 58 43 4f 4d 50 4c 45 |or..X */|.XCOMPLE|
|00002b10| 58 20 2a 0a 58 63 70 6f | 77 65 72 28 63 31 2c 20 |X *.Xcpo|wer(c1, |
|00002b20| 63 32 2c 20 65 70 73 69 | 6c 6f 6e 29 0a 58 09 43 |c2, epsi|lon).X.C|
|00002b30| 4f 4d 50 4c 45 58 20 2a | 63 31 2c 20 2a 63 32 3b |OMPLEX *|c1, *c2;|
|00002b40| 0a 58 09 4e 55 4d 42 45 | 52 20 2a 65 70 73 69 6c |.X.NUMBE|R *epsil|
|00002b50| 6f 6e 3b 0a 58 7b 0a 58 | 09 43 4f 4d 50 4c 45 58 |on;.X{.X|.COMPLEX|
|00002b60| 20 2a 74 6d 70 31 2c 20 | 2a 74 6d 70 32 3b 0a 58 | *tmp1, |*tmp2;.X|
|00002b70| 09 4e 55 4d 42 45 52 20 | 2a 65 70 73 69 6c 6f 6e |.NUMBER |*epsilon|
|00002b80| 32 3b 0a 58 0a 58 09 69 | 66 20 28 63 69 73 72 65 |2;.X.X.i|f (cisre|
|00002b90| 61 6c 28 63 32 29 20 26 | 26 20 71 69 73 69 6e 74 |al(c2) &|& qisint|
|00002ba0| 28 63 32 2d 3e 72 65 61 | 6c 29 29 0a 58 09 09 72 |(c2->rea|l)).X..r|
|00002bb0| 65 74 75 72 6e 20 63 70 | 6f 77 69 28 63 31 2c 20 |eturn cp|owi(c1, |
|00002bc0| 63 32 2d 3e 72 65 61 6c | 29 3b 0a 58 09 69 66 20 |c2->real|);.X.if |
|00002bd0| 28 63 69 73 6f 6e 65 28 | 63 31 29 20 7c 7c 20 63 |(cisone(|c1) || c|
|00002be0| 69 73 7a 65 72 6f 28 63 | 31 29 29 0a 58 09 09 72 |iszero(c|1)).X..r|
|00002bf0| 65 74 75 72 6e 20 63 6c | 69 6e 6b 28 63 31 29 3b |eturn cl|ink(c1);|
|00002c00| 0a 58 09 65 70 73 69 6c | 6f 6e 32 20 3d 20 71 73 |.X.epsil|on2 = qs|
|00002c10| 63 61 6c 65 28 65 70 73 | 69 6c 6f 6e 2c 20 2d 34 |cale(eps|ilon, -4|
|00002c20| 4c 29 3b 0a 58 09 74 6d | 70 31 20 3d 20 63 6c 6e |L);.X.tm|p1 = cln|
|00002c30| 28 63 31 2c 20 65 70 73 | 69 6c 6f 6e 32 29 3b 0a |(c1, eps|ilon2);.|
|00002c40| 58 09 74 6d 70 32 20 3d | 20 63 6d 75 6c 28 74 6d |X.tmp2 =| cmul(tm|
|00002c50| 70 31 2c 20 63 32 29 3b | 0a 58 09 63 6f 6d 66 72 |p1, c2);|.X.comfr|
|00002c60| 65 65 28 74 6d 70 31 29 | 3b 0a 58 09 74 6d 70 31 |ee(tmp1)|;.X.tmp1|
|00002c70| 20 3d 20 63 65 78 70 28 | 74 6d 70 32 2c 20 65 70 | = cexp(|tmp2, ep|
|00002c80| 73 69 6c 6f 6e 29 3b 0a | 58 09 63 6f 6d 66 72 65 |silon);.|X.comfre|
|00002c90| 65 28 74 6d 70 32 29 3b | 0a 58 09 71 66 72 65 65 |e(tmp2);|.X.qfree|
|00002ca0| 28 65 70 73 69 6c 6f 6e | 32 29 3b 0a 58 09 72 65 |(epsilon|2);.X.re|
|00002cb0| 74 75 72 6e 20 74 6d 70 | 31 3b 0a 58 7d 0a 58 0a |turn tmp|1;.X}.X.|
|00002cc0| 58 0a 58 2f 2a 0a 58 20 | 2a 20 50 72 69 6e 74 20 |X.X/*.X |* Print |
|00002cd0| 61 20 63 6f 6d 70 6c 65 | 78 20 6e 75 6d 62 65 72 |a comple|x number|
|00002ce0| 2e 0a 58 20 2a 2f 0a 58 | 76 6f 69 64 0a 58 63 6f |..X */.X|void.Xco|
|00002cf0| 6d 70 72 69 6e 74 28 63 | 29 0a 58 09 43 4f 4d 50 |mprint(c|).X.COMP|
|00002d00| 4c 45 58 20 2a 63 3b 0a | 58 7b 0a 58 09 4e 55 4d |LEX *c;.|X{.X.NUM|
|00002d10| 42 45 52 20 71 74 6d 70 | 3b 0a 58 0a 58 09 69 66 |BER qtmp|;.X.X.if|
|00002d20| 20 28 5f 6f 75 74 6d 6f | 64 65 5f 20 3d 3d 20 4d | (_outmo|de_ == M|
|00002d30| 4f 44 45 5f 46 52 41 43 | 29 20 7b 0a 58 09 09 63 |ODE_FRAC|) {.X..c|
|00002d40| 70 72 69 6e 74 66 72 28 | 63 29 3b 0a 58 09 09 72 |printfr(|c);.X..r|
|00002d50| 65 74 75 72 6e 3b 0a 58 | 09 7d 0a 58 09 69 66 20 |eturn;.X|.}.X.if |
|00002d60| 28 21 71 69 73 7a 65 72 | 6f 28 63 2d 3e 72 65 61 |(!qiszer|o(c->rea|
|00002d70| 6c 29 20 7c 7c 20 71 69 | 73 7a 65 72 6f 28 63 2d |l) || qi|szero(c-|
|00002d80| 3e 69 6d 61 67 29 29 0a | 58 09 09 71 70 72 69 6e |>imag)).|X..qprin|
|00002d90| 74 6e 75 6d 28 63 2d 3e | 72 65 61 6c 2c 20 4d 4f |tnum(c->|real, MO|
|00002da0| 44 45 5f 44 45 46 41 55 | 4c 54 29 3b 0a 58 09 71 |DE_DEFAU|LT);.X.q|
|00002db0| 74 6d 70 20 3d 20 63 2d | 3e 69 6d 61 67 5b 30 5d |tmp = c-|>imag[0]|
|00002dc0| 3b 0a 58 09 69 66 20 28 | 71 69 73 7a 65 72 6f 28 |;.X.if (|qiszero(|
|00002dd0| 26 71 74 6d 70 29 29 0a | 58 09 09 72 65 74 75 72 |&qtmp)).|X..retur|
|00002de0| 6e 3b 0a 58 09 69 66 20 | 28 21 71 69 73 7a 65 72 |n;.X.if |(!qiszer|
|00002df0| 6f 28 63 2d 3e 72 65 61 | 6c 29 20 26 26 20 21 71 |o(c->rea|l) && !q|
|00002e00| 69 73 6e 65 67 28 26 71 | 74 6d 70 29 29 0a 58 09 |isneg(&q|tmp)).X.|
|00002e10| 09 6d 61 74 68 5f 63 68 | 72 28 27 2b 27 29 3b 0a |.math_ch|r('+');.|
|00002e20| 58 09 69 66 20 28 71 69 | 73 6e 65 67 28 26 71 74 |X.if (qi|sneg(&qt|
|00002e30| 6d 70 29 29 20 7b 0a 58 | 09 09 6d 61 74 68 5f 63 |mp)) {.X|..math_c|
|00002e40| 68 72 28 27 2d 27 29 3b | 0a 58 09 09 71 74 6d 70 |hr('-');|.X..qtmp|
|00002e50| 2e 6e 75 6d 2e 73 69 67 | 6e 20 3d 20 30 3b 0a 58 |.num.sig|n = 0;.X|
|00002e60| 09 7d 0a 58 09 71 70 72 | 69 6e 74 6e 75 6d 28 26 |.}.X.qpr|intnum(&|
|00002e70| 71 74 6d 70 2c 20 4d 4f | 44 45 5f 44 45 46 41 55 |qtmp, MO|DE_DEFAU|
|00002e80| 4c 54 29 3b 0a 58 09 6d | 61 74 68 5f 63 68 72 28 |LT);.X.m|ath_chr(|
|00002e90| 27 69 27 29 3b 0a 58 7d | 0a 58 0a 58 2f 2a 20 45 |'i');.X}|.X.X/* E|
|00002ea0| 4e 44 20 43 4f 44 45 20 | 2a 2f 0a 45 4e 44 5f 4f |ND CODE |*/.END_O|
|00002eb0| 46 5f 46 49 4c 45 0a 69 | 66 20 74 65 73 74 20 31 |F_FILE.i|f test 1|
|00002ec0| 30 33 33 32 20 2d 6e 65 | 20 60 77 63 20 2d 63 20 |0332 -ne| `wc -c |
|00002ed0| 3c 27 63 6f 6d 66 75 6e | 63 2e 63 27 60 3b 20 74 |<'comfun|c.c'`; t|
|00002ee0| 68 65 6e 0a 20 20 20 20 | 65 63 68 6f 20 73 68 61 |hen.    |echo sha|
|00002ef0| 72 3a 20 5c 22 27 63 6f | 6d 66 75 6e 63 2e 63 27 |r: \"'co|mfunc.c'|
|00002f00| 5c 22 20 75 6e 70 61 63 | 6b 65 64 20 77 69 74 68 |\" unpac|ked with|
|00002f10| 20 77 72 6f 6e 67 20 73 | 69 7a 65 21 0a 66 69 0a | wrong s|ize!.fi.|
|00002f20| 23 20 65 6e 64 20 6f 66 | 20 27 63 6f 6d 66 75 6e |# end of| 'comfun|
|00002f30| 63 2e 63 27 0a 66 69 0a | 69 66 20 74 65 73 74 20 |c.c'.fi.|if test |
|00002f40| 2d 66 20 27 63 6f 6d 6d | 61 74 68 2e 63 27 20 2d |-f 'comm|ath.c' -|
|00002f50| 61 20 22 24 7b 31 7d 22 | 20 21 3d 20 22 2d 63 22 |a "${1}"| != "-c"|
|00002f60| 20 3b 20 74 68 65 6e 20 | 0a 20 20 65 63 68 6f 20 | ; then |.  echo |
|00002f70| 73 68 61 72 3a 20 57 69 | 6c 6c 20 6e 6f 74 20 63 |shar: Wi|ll not c|
|00002f80| 6c 6f 62 62 65 72 20 65 | 78 69 73 74 69 6e 67 20 |lobber e|xisting |
|00002f90| 66 69 6c 65 20 5c 22 27 | 63 6f 6d 6d 61 74 68 2e |file \"'|commath.|
|00002fa0| 63 27 5c 22 0a 65 6c 73 | 65 0a 65 63 68 6f 20 73 |c'\".els|e.echo s|
|00002fb0| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|00002fc0| 5c 22 27 63 6f 6d 6d 61 | 74 68 2e 63 27 5c 22 20 |\"'comma|th.c'\" |
|00002fd0| 5c 28 39 36 35 34 20 63 | 68 61 72 61 63 74 65 72 |\(9654 c|haracter|
|00002fe0| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|00002ff0| 20 3e 27 63 6f 6d 6d 61 | 74 68 2e 63 27 20 3c 3c | >'comma|th.c' <<|
|00003000| 27 45 4e 44 5f 4f 46 5f | 46 49 4c 45 27 0a 58 2f |'END_OF_|FILE'.X/|
|00003010| 2a 0a 58 20 2a 20 43 6f | 70 79 72 69 67 68 74 20 |*.X * Co|pyright |
|00003020| 28 63 29 20 31 39 39 32 | 20 44 61 76 69 64 20 49 |(c) 1992| David I|
|00003030| 2e 20 42 65 6c 6c 0a 58 | 20 2a 20 50 65 72 6d 69 |. Bell.X| * Permi|
|00003040| 73 73 69 6f 6e 20 69 73 | 20 67 72 61 6e 74 65 64 |ssion is| granted|
|00003050| 20 74 6f 20 75 73 65 2c | 20 64 69 73 74 72 69 62 | to use,| distrib|
|00003060| 75 74 65 2c 20 6f 72 20 | 6d 6f 64 69 66 79 20 74 |ute, or |modify t|
|00003070| 68 69 73 20 73 6f 75 72 | 63 65 2c 0a 58 20 2a 20 |his sour|ce,.X * |
|00003080| 70 72 6f 76 69 64 65 64 | 20 74 68 61 74 20 74 68 |provided| that th|
|00003090| 69 73 20 63 6f 70 79 72 | 69 67 68 74 20 6e 6f 74 |is copyr|ight not|
|000030a0| 69 63 65 20 72 65 6d 61 | 69 6e 73 20 69 6e 74 61 |ice rema|ins inta|
|000030b0| 63 74 2e 0a 58 20 2a 0a | 58 20 2a 20 45 78 74 65 |ct..X *.|X * Exte|
|000030c0| 6e 64 65 64 20 70 72 65 | 63 69 73 69 6f 6e 20 63 |nded pre|cision c|
|000030d0| 6f 6d 70 6c 65 78 20 61 | 72 69 74 68 6d 65 74 69 |omplex a|rithmeti|
|000030e0| 63 20 70 72 69 6d 69 74 | 69 76 65 20 72 6f 75 74 |c primit|ive rout|
|000030f0| 69 6e 65 73 0a 58 20 2a | 2f 0a 58 0a 58 23 69 6e |ines.X *|/.X.X#in|
|00003100| 63 6c 75 64 65 20 3c 73 | 74 64 69 6f 2e 68 3e 0a |clude <s|tdio.h>.|
|00003110| 58 23 69 6e 63 6c 75 64 | 65 20 22 6d 61 74 68 2e |X#includ|e "math.|
|00003120| 68 22 0a 58 0a 58 0a 58 | 43 4f 4d 50 4c 45 58 20 |h".X.X.X|COMPLEX |
|00003130| 5f 63 7a 65 72 6f 5f 20 | 3d 09 09 7b 20 26 5f 71 |_czero_ |=..{ &_q|
|00003140| 7a 65 72 6f 5f 2c 20 26 | 5f 71 7a 65 72 6f 5f 2c |zero_, &|_qzero_,|
|00003150| 20 31 20 7d 3b 0a 58 43 | 4f 4d 50 4c 45 58 20 5f | 1 };.XC|OMPLEX _|
|00003160| 63 6f 6e 65 5f 20 3d 09 | 09 7b 20 26 5f 71 6f 6e |cone_ =.|.{ &_qon|
|00003170| 65 5f 2c 20 26 5f 71 7a | 65 72 6f 5f 2c 20 31 20 |e_, &_qz|ero_, 1 |
|00003180| 7d 3b 0a 58 73 74 61 74 | 69 63 20 43 4f 4d 50 4c |};.Xstat|ic COMPL|
|00003190| 45 58 20 5f 63 6e 65 67 | 6f 6e 65 5f 20 3d 09 7b |EX _cneg|one_ =.{|
|000031a0| 20 26 5f 71 6e 65 67 6f | 6e 65 5f 2c 20 26 5f 71 | &_qnego|ne_, &_q|
|000031b0| 7a 65 72 6f 5f 2c 20 31 | 20 7d 3b 0a 58 0a 58 23 |zero_, 1| };.X.X#|
|000031c0| 69 66 20 30 0a 58 43 4f | 4d 50 4c 45 58 20 5f 63 |if 0.XCO|MPLEX _c|
|000031d0| 6f 6e 65 69 5f 20 3d 09 | 7b 20 26 5f 71 7a 65 72 |onei_ =.|{ &_qzer|
|000031e0| 6f 5f 2c 20 26 5f 71 6f | 6e 65 5f 2c 20 31 20 7d |o_, &_qo|ne_, 1 }|
|000031f0| 3b 0a 58 23 65 6e 64 69 | 66 0a 58 0a 58 0a 58 2f |;.X#endi|f.X.X.X/|
|00003200| 2a 0a 58 20 2a 20 46 72 | 65 65 20 6c 69 73 74 20 |*.X * Fr|ee list |
|00003210| 66 6f 72 20 63 6f 6d 70 | 6c 65 78 20 6e 75 6d 62 |for comp|lex numb|
|00003220| 65 72 73 2e 0a 58 20 2a | 2f 0a 58 73 74 61 74 69 |ers..X *|/.Xstati|
|00003230| 63 20 46 52 45 45 4c 49 | 53 54 20 66 72 65 65 6c |c FREELI|ST freel|
|00003240| 69 73 74 20 3d 20 7b 0a | 58 09 73 69 7a 65 6f 66 |ist = {.|X.sizeof|
|00003250| 28 43 4f 4d 50 4c 45 58 | 29 2c 09 2f 2a 20 73 69 |(COMPLEX|),./* si|
|00003260| 7a 65 20 6f 66 20 61 6e | 20 69 74 65 6d 20 2a 2f |ze of an| item */|
|00003270| 0a 58 09 31 30 30 09 09 | 09 2f 2a 20 6e 75 6d 62 |.X.100..|./* numb|
|00003280| 65 72 20 6f 66 20 66 72 | 65 65 20 69 74 65 6d 73 |er of fr|ee items|
|00003290| 20 74 6f 20 6b 65 65 70 | 20 2a 2f 0a 58 7d 3b 0a | to keep| */.X};.|
|000032a0| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 41 64 64 20 |X.X.X/*.|X * Add |
|000032b0| 74 77 6f 20 63 6f 6d 70 | 6c 65 78 20 6e 75 6d 62 |two comp|lex numb|
|000032c0| 65 72 73 2e 0a 58 20 2a | 2f 0a 58 43 4f 4d 50 4c |ers..X *|/.XCOMPL|
|000032d0| 45 58 20 2a 0a 58 63 61 | 64 64 28 63 31 2c 20 63 |EX *.Xca|dd(c1, c|
|000032e0| 32 29 0a 58 09 43 4f 4d | 50 4c 45 58 20 2a 63 31 |2).X.COM|PLEX *c1|
|000032f0| 2c 20 2a 63 32 3b 0a 58 | 7b 0a 58 09 43 4f 4d 50 |, *c2;.X|{.X.COMP|
|00003300| 4c 45 58 20 2a 72 3b 0a | 58 0a 58 09 69 66 20 28 |LEX *r;.|X.X.if (|
|00003310| 63 69 73 7a 65 72 6f 28 | 63 31 29 29 0a 58 09 09 |ciszero(|c1)).X..|
|00003320| 72 65 74 75 72 6e 20 63 | 6c 69 6e 6b 28 63 32 29 |return c|link(c2)|
|00003330| 3b 0a 58 09 69 66 20 28 | 63 69 73 7a 65 72 6f 28 |;.X.if (|ciszero(|
|00003340| 63 32 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 63 |c2)).X..|return c|
|00003350| 6c 69 6e 6b 28 63 31 29 | 3b 0a 58 09 72 20 3d 20 |link(c1)|;.X.r = |
|00003360| 63 6f 6d 61 6c 6c 6f 63 | 28 29 3b 0a 58 09 69 66 |comalloc|();.X.if|
|00003370| 20 28 21 71 69 73 7a 65 | 72 6f 28 63 31 2d 3e 72 | (!qisze|ro(c1->r|
|00003380| 65 61 6c 29 20 7c 7c 20 | 21 71 69 73 7a 65 72 6f |eal) || |!qiszero|
|00003390| 28 63 32 2d 3e 72 65 61 | 6c 29 29 0a 58 09 09 72 |(c2->rea|l)).X..r|
|000033a0| 2d 3e 72 65 61 6c 20 3d | 20 71 61 64 64 28 63 31 |->real =| qadd(c1|
|000033b0| 2d 3e 72 65 61 6c 2c 20 | 63 32 2d 3e 72 65 61 6c |->real, |c2->real|
|000033c0| 29 3b 0a 58 09 69 66 20 | 28 21 71 69 73 7a 65 72 |);.X.if |(!qiszer|
|000033d0| 6f 28 63 31 2d 3e 69 6d | 61 67 29 20 7c 7c 20 21 |o(c1->im|ag) || !|
|000033e0| 71 69 73 7a 65 72 6f 28 | 63 32 2d 3e 69 6d 61 67 |qiszero(|c2->imag|
|000033f0| 29 29 0a 58 09 09 72 2d | 3e 69 6d 61 67 20 3d 20 |)).X..r-|>imag = |
|00003400| 71 61 64 64 28 63 31 2d | 3e 69 6d 61 67 2c 20 63 |qadd(c1-|>imag, c|
|00003410| 32 2d 3e 69 6d 61 67 29 | 3b 0a 58 09 72 65 74 75 |2->imag)|;.X.retu|
|00003420| 72 6e 20 72 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a |rn r;.X}|.X.X.X/*|
|00003430| 0a 58 20 2a 20 53 75 62 | 74 72 61 63 74 20 74 77 |.X * Sub|tract tw|
|00003440| 6f 20 63 6f 6d 70 6c 65 | 78 20 6e 75 6d 62 65 72 |o comple|x number|
|00003450| 73 2e 0a 58 20 2a 2f 0a | 58 43 4f 4d 50 4c 45 58 |s..X */.|XCOMPLEX|
|00003460| 20 2a 0a 58 63 73 75 62 | 28 63 31 2c 20 63 32 29 | *.Xcsub|(c1, c2)|
|00003470| 0a 58 09 43 4f 4d 50 4c | 45 58 20 2a 63 31 2c 20 |.X.COMPL|EX *c1, |
|00003480| 2a 63 32 3b 0a 58 7b 0a | 58 09 43 4f 4d 50 4c 45 |*c2;.X{.|X.COMPLE|
|00003490| 58 20 2a 72 3b 0a 58 0a | 58 09 69 66 20 28 28 63 |X *r;.X.|X.if ((c|
|000034a0| 31 2d 3e 72 65 61 6c 20 | 3d 3d 20 63 32 2d 3e 72 |1->real |== c2->r|
|000034b0| 65 61 6c 29 20 26 26 20 | 28 63 31 2d 3e 69 6d 61 |eal) && |(c1->ima|
|000034c0| 67 20 3d 3d 20 63 32 2d | 3e 69 6d 61 67 29 29 0a |g == c2-|>imag)).|
|000034d0| 58 09 09 72 65 74 75 72 | 6e 20 63 6c 69 6e 6b 28 |X..retur|n clink(|
|000034e0| 26 5f 63 7a 65 72 6f 5f | 29 3b 0a 58 09 69 66 20 |&_czero_|);.X.if |
|000034f0| 28 63 69 73 7a 65 72 6f | 28 63 32 29 29 0a 58 09 |(ciszero|(c2)).X.|
|00003500| 09 72 65 74 75 72 6e 20 | 63 6c 69 6e 6b 28 63 31 |.return |clink(c1|
|00003510| 29 3b 0a 58 09 72 20 3d | 20 63 6f 6d 61 6c 6c 6f |);.X.r =| comallo|
|00003520| 63 28 29 3b 0a 58 09 69 | 66 20 28 21 71 69 73 7a |c();.X.i|f (!qisz|
|00003530| 65 72 6f 28 63 31 2d 3e | 72 65 61 6c 29 20 7c 7c |ero(c1->|real) |||
|00003540| 20 21 71 69 73 7a 65 72 | 6f 28 63 32 2d 3e 72 65 | !qiszer|o(c2->re|
|00003550| 61 6c 29 29 0a 58 09 09 | 72 2d 3e 72 65 61 6c 20 |al)).X..|r->real |
|00003560| 3d 20 71 73 75 62 28 63 | 31 2d 3e 72 65 61 6c 2c |= qsub(c|1->real,|
|00003570| 20 63 32 2d 3e 72 65 61 | 6c 29 3b 0a 58 09 69 66 | c2->rea|l);.X.if|
|00003580| 20 28 21 71 69 73 7a 65 | 72 6f 28 63 31 2d 3e 69 | (!qisze|ro(c1->i|
|00003590| 6d 61 67 29 20 7c 7c 20 | 21 71 69 73 7a 65 72 6f |mag) || |!qiszero|
|000035a0| 28 63 32 2d 3e 69 6d 61 | 67 29 29 0a 58 09 09 72 |(c2->ima|g)).X..r|
|000035b0| 2d 3e 69 6d 61 67 20 3d | 20 71 73 75 62 28 63 31 |->imag =| qsub(c1|
|000035c0| 2d 3e 69 6d 61 67 2c 20 | 63 32 2d 3e 69 6d 61 67 |->imag, |c2->imag|
|000035d0| 29 3b 0a 58 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |);.X.ret|urn r;.X|
|000035e0| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 4d 75 |}.X.X.X/|*.X * Mu|
|000035f0| 6c 74 69 70 6c 79 20 74 | 77 6f 20 63 6f 6d 70 6c |ltiply t|wo compl|
|00003600| 65 78 20 6e 75 6d 62 65 | 72 73 2e 0a 58 20 2a 20 |ex numbe|rs..X * |
|00003610| 54 68 69 73 20 73 61 76 | 65 73 20 6f 6e 65 20 6d |This sav|es one m|
|00003620| 75 6c 74 69 70 6c 69 63 | 61 74 69 6f 6e 20 6f 76 |ultiplic|ation ov|
|00003630| 65 72 20 74 68 65 20 6f | 62 76 69 6f 75 73 20 61 |er the o|bvious a|
|00003640| 6c 67 6f 72 69 74 68 6d | 20 62 79 0a 58 20 2a 20 |lgorithm| by.X * |
|00003650| 74 72 61 64 69 6e 67 20 | 69 74 20 66 6f 72 20 73 |trading |it for s|
|00003660| 65 76 65 72 61 6c 20 65 | 78 74 72 61 20 61 64 64 |everal e|xtra add|
|00003670| 69 74 69 6f 6e 73 2c 20 | 61 73 20 66 6f 6c 6c 6f |itions, |as follo|
|00003680| 77 73 2e 20 20 4c 65 74 | 0a 58 20 2a 09 71 31 20 |ws.  Let|.X *.q1 |
|00003690| 3d 20 28 61 20 2b 20 62 | 29 20 2a 20 28 63 20 2b |= (a + b|) * (c +|
|000036a0| 20 64 29 0a 58 20 2a 09 | 71 32 20 3d 20 61 20 2a | d).X *.|q2 = a *|
|000036b0| 20 63 0a 58 20 2a 09 71 | 33 20 3d 20 62 20 2a 20 | c.X *.q|3 = b * |
|000036c0| 64 0a 58 20 2a 20 54 68 | 65 6e 20 28 61 2b 62 69 |d.X * Th|en (a+bi|
|000036d0| 29 20 2a 20 28 63 2b 64 | 69 29 20 3d 20 28 71 32 |) * (c+d|i) = (q2|
|000036e0| 20 2d 20 71 33 29 20 2b | 20 28 71 31 20 2d 20 71 | - q3) +| (q1 - q|
|000036f0| 32 20 2d 20 71 33 29 69 | 2e 0a 58 20 2a 2f 0a 58 |2 - q3)i|..X */.X|
|00003700| 43 4f 4d 50 4c 45 58 20 | 2a 0a 58 63 6d 75 6c 28 |COMPLEX |*.Xcmul(|
|00003710| 63 31 2c 20 63 32 29 0a | 58 09 43 4f 4d 50 4c 45 |c1, c2).|X.COMPLE|
|00003720| 58 20 2a 63 31 2c 20 2a | 63 32 3b 0a 58 7b 0a 58 |X *c1, *|c2;.X{.X|
|00003730| 09 43 4f 4d 50 4c 45 58 | 20 2a 72 3b 0a 58 09 4e |.COMPLEX| *r;.X.N|
|00003740| 55 4d 42 45 52 20 2a 71 | 31 2c 20 2a 71 32 2c 20 |UMBER *q|1, *q2, |
|00003750| 2a 71 33 2c 20 2a 71 34 | 3b 0a 58 0a 58 09 69 66 |*q3, *q4|;.X.X.if|
|00003760| 20 28 63 69 73 7a 65 72 | 6f 28 63 31 29 20 7c 7c | (ciszer|o(c1) |||
|00003770| 20 63 69 73 7a 65 72 6f | 28 63 32 29 29 0a 58 09 | ciszero|(c2)).X.|
|00003780| 09 72 65 74 75 72 6e 20 | 63 6c 69 6e 6b 28 26 5f |.return |clink(&_|
|00003790| 63 7a 65 72 6f 5f 29 3b | 0a 58 09 69 66 20 28 63 |czero_);|.X.if (c|
|000037a0| 69 73 6f 6e 65 28 63 31 | 29 29 0a 58 09 09 72 65 |isone(c1|)).X..re|
|000037b0| 74 75 72 6e 20 63 6c 69 | 6e 6b 28 63 32 29 3b 0a |turn cli|nk(c2);.|
|000037c0| 58 09 69 66 20 28 63 69 | 73 6f 6e 65 28 63 32 29 |X.if (ci|sone(c2)|
|000037d0| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 63 6c 69 6e |).X..ret|urn clin|
|000037e0| 6b 28 63 31 29 3b 0a 58 | 09 69 66 20 28 63 69 73 |k(c1);.X|.if (cis|
|000037f0| 72 65 61 6c 28 63 32 29 | 29 0a 58 09 09 72 65 74 |real(c2)|).X..ret|
|00003800| 75 72 6e 20 63 6d 75 6c | 71 28 63 31 2c 20 63 32 |urn cmul|q(c1, c2|
|00003810| 2d 3e 72 65 61 6c 29 3b | 0a 58 09 69 66 20 28 63 |->real);|.X.if (c|
|00003820| 69 73 72 65 61 6c 28 63 | 31 29 29 0a 58 09 09 72 |isreal(c|1)).X..r|
|00003830| 65 74 75 72 6e 20 63 6d | 75 6c 71 28 63 32 2c 20 |eturn cm|ulq(c2, |
|00003840| 63 31 2d 3e 72 65 61 6c | 29 3b 0a 58 09 2f 2a 0a |c1->real|);.X./*.|
|00003850| 58 09 20 2a 20 4e 65 65 | 64 20 74 6f 20 64 6f 20 |X. * Nee|d to do |
|00003860| 74 68 65 20 66 75 6c 6c | 20 63 61 6c 63 75 6c 61 |the full| calcula|
|00003870| 74 69 6f 6e 2e 0a 58 09 | 20 2a 2f 0a 58 09 72 20 |tion..X.| */.X.r |
|00003880| 3d 20 63 6f 6d 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |= comall|oc();.X.|
|00003890| 71 32 20 3d 20 71 61 64 | 64 28 63 31 2d 3e 72 65 |q2 = qad|d(c1->re|
|000038a0| 61 6c 2c 20 63 31 2d 3e | 69 6d 61 67 29 3b 0a 58 |al, c1->|imag);.X|
|000038b0| 09 71 33 20 3d 20 71 61 | 64 64 28 63 32 2d 3e 72 |.q3 = qa|dd(c2->r|
|000038c0| 65 61 6c 2c 20 63 32 2d | 3e 69 6d 61 67 29 3b 0a |eal, c2-|>imag);.|
|000038d0| 58 09 71 31 20 3d 20 71 | 6d 75 6c 28 71 32 2c 20 |X.q1 = q|mul(q2, |
|000038e0| 71 33 29 3b 0a 58 09 71 | 66 72 65 65 28 71 32 29 |q3);.X.q|free(q2)|
|000038f0| 3b 0a 58 09 71 66 72 65 | 65 28 71 33 29 3b 0a 58 |;.X.qfre|e(q3);.X|
|00003900| 09 71 32 20 3d 20 71 6d | 75 6c 28 63 31 2d 3e 72 |.q2 = qm|ul(c1->r|
|00003910| 65 61 6c 2c 20 63 32 2d | 3e 72 65 61 6c 29 3b 0a |eal, c2-|>real);.|
|00003920| 58 09 71 33 20 3d 20 71 | 6d 75 6c 28 63 31 2d 3e |X.q3 = q|mul(c1->|
|00003930| 69 6d 61 67 2c 20 63 32 | 2d 3e 69 6d 61 67 29 3b |imag, c2|->imag);|
|00003940| 0a 58 09 71 34 20 3d 20 | 71 61 64 64 28 71 32 2c |.X.q4 = |qadd(q2,|
|00003950| 20 71 33 29 3b 0a 58 09 | 72 2d 3e 72 65 61 6c 20 | q3);.X.|r->real |
|00003960| 3d 20 71 73 75 62 28 71 | 32 2c 20 71 33 29 3b 0a |= qsub(q|2, q3);.|
|00003970| 58 09 72 2d 3e 69 6d 61 | 67 20 3d 20 71 73 75 62 |X.r->ima|g = qsub|
|00003980| 28 71 31 2c 20 71 34 29 | 3b 0a 58 09 71 66 72 65 |(q1, q4)|;.X.qfre|
|00003990| 65 28 71 31 29 3b 0a 58 | 09 71 66 72 65 65 28 71 |e(q1);.X|.qfree(q|
|000039a0| 32 29 3b 0a 58 09 71 66 | 72 65 65 28 71 33 29 3b |2);.X.qf|ree(q3);|
|000039b0| 0a 58 09 71 66 72 65 65 | 28 71 34 29 3b 0a 58 09 |.X.qfree|(q4);.X.|
|000039c0| 72 65 74 75 72 6e 20 72 | 3b 0a 58 7d 0a 58 0a 58 |return r|;.X}.X.X|
|000039d0| 0a 58 2f 2a 0a 58 20 2a | 20 53 71 75 61 72 65 20 |.X/*.X *| Square |
|000039e0| 61 20 63 6f 6d 70 6c 65 | 78 20 6e 75 6d 62 65 72 |a comple|x number|
|000039f0| 2e 0a 58 20 2a 2f 0a 58 | 43 4f 4d 50 4c 45 58 20 |..X */.X|COMPLEX |
|00003a00| 2a 0a 58 63 73 71 75 61 | 72 65 28 63 29 0a 58 09 |*.Xcsqua|re(c).X.|
|00003a10| 43 4f 4d 50 4c 45 58 20 | 2a 63 3b 0a 58 7b 0a 58 |COMPLEX |*c;.X{.X|
|00003a20| 09 43 4f 4d 50 4c 45 58 | 20 2a 72 3b 0a 58 09 4e |.COMPLEX| *r;.X.N|
|00003a30| 55 4d 42 45 52 20 2a 71 | 31 2c 20 2a 71 32 3b 0a |UMBER *q|1, *q2;.|
|00003a40| 58 0a 58 09 69 66 20 28 | 63 69 73 7a 65 72 6f 28 |X.X.if (|ciszero(|
|00003a50| 63 29 29 0a 58 09 09 72 | 65 74 75 72 6e 20 63 6c |c)).X..r|eturn cl|
|00003a60| 69 6e 6b 28 26 5f 63 7a | 65 72 6f 5f 29 3b 0a 58 |ink(&_cz|ero_);.X|
|00003a70| 09 69 66 20 28 63 69 73 | 72 75 6e 69 74 28 63 29 |.if (cis|runit(c)|
|00003a80| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 63 6c 69 6e |).X..ret|urn clin|
|00003a90| 6b 28 26 5f 63 6f 6e 65 | 5f 29 3b 0a 58 09 69 66 |k(&_cone|_);.X.if|
|00003aa0| 20 28 63 69 73 69 75 6e | 69 74 28 63 29 29 0a 58 | (cisiun|it(c)).X|
|00003ab0| 09 09 72 65 74 75 72 6e | 20 63 6c 69 6e 6b 28 26 |..return| clink(&|
|00003ac0| 5f 63 6e 65 67 6f 6e 65 | 5f 29 3b 0a 58 09 72 20 |_cnegone|_);.X.r |
|00003ad0| 3d 20 63 6f 6d 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |= comall|oc();.X.|
|00003ae0| 69 66 20 28 63 69 73 72 | 65 61 6c 28 63 29 29 20 |if (cisr|eal(c)) |
|00003af0| 7b 0a 58 09 09 72 2d 3e | 72 65 61 6c 20 3d 20 71 |{.X..r->|real = q|
|00003b00| 73 71 75 61 72 65 28 63 | 2d 3e 72 65 61 6c 29 3b |square(c|->real);|
|00003b10| 0a 58 09 09 72 65 74 75 | 72 6e 20 72 3b 0a 58 09 |.X..retu|rn r;.X.|
|00003b20| 7d 0a 58 09 69 66 20 28 | 63 69 73 69 6d 61 67 28 |}.X.if (|cisimag(|
|00003b30| 63 29 29 20 7b 0a 58 09 | 09 71 31 20 3d 20 71 73 |c)) {.X.|.q1 = qs|
|00003b40| 71 75 61 72 65 28 63 2d | 3e 69 6d 61 67 29 3b 0a |quare(c-|>imag);.|
|00003b50| 58 09 09 72 2d 3e 72 65 | 61 6c 20 3d 20 71 6e 65 |X..r->re|al = qne|
|00003b60| 67 28 71 31 29 3b 0a 58 | 09 09 71 66 72 65 65 28 |g(q1);.X|..qfree(|
|00003b70| 71 31 29 3b 0a 58 09 09 | 72 65 74 75 72 6e 20 72 |q1);.X..|return r|
|00003b80| 3b 0a 58 09 7d 0a 58 09 | 71 31 20 3d 20 71 73 71 |;.X.}.X.|q1 = qsq|
|00003b90| 75 61 72 65 28 63 2d 3e | 72 65 61 6c 29 3b 0a 58 |uare(c->|real);.X|
|00003ba0| 09 71 32 20 3d 20 71 73 | 71 75 61 72 65 28 63 2d |.q2 = qs|quare(c-|
|00003bb0| 3e 69 6d 61 67 29 3b 0a | 58 09 72 2d 3e 72 65 61 |>imag);.|X.r->rea|
|00003bc0| 6c 20 3d 20 71 73 75 62 | 28 71 31 2c 20 71 32 29 |l = qsub|(q1, q2)|
|00003bd0| 3b 0a 58 09 71 66 72 65 | 65 28 71 31 29 3b 0a 58 |;.X.qfre|e(q1);.X|
|00003be0| 09 71 66 72 65 65 28 71 | 32 29 3b 0a 58 09 71 31 |.qfree(q|2);.X.q1|
|00003bf0| 20 3d 20 71 6d 75 6c 28 | 63 2d 3e 72 65 61 6c 2c | = qmul(|c->real,|
|00003c00| 20 63 2d 3e 69 6d 61 67 | 29 3b 0a 58 09 72 2d 3e | c->imag|);.X.r->|
|00003c10| 69 6d 61 67 20 3d 20 71 | 73 63 61 6c 65 28 71 31 |imag = q|scale(q1|
|00003c20| 2c 20 31 4c 29 3b 0a 58 | 09 71 66 72 65 65 28 71 |, 1L);.X|.qfree(q|
|00003c30| 31 29 3b 0a 58 09 72 65 | 74 75 72 6e 20 72 3b 0a |1);.X.re|turn r;.|
|00003c40| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 44 |X}.X.X.X|/*.X * D|
|00003c50| 69 76 69 64 65 20 74 77 | 6f 20 63 6f 6d 70 6c 65 |ivide tw|o comple|
|00003c60| 78 20 6e 75 6d 62 65 72 | 73 2e 0a 58 20 2a 2f 0a |x number|s..X */.|
|00003c70| 58 43 4f 4d 50 4c 45 58 | 20 2a 0a 58 63 64 69 76 |XCOMPLEX| *.Xcdiv|
|00003c80| 28 63 31 2c 20 63 32 29 | 0a 58 09 43 4f 4d 50 4c |(c1, c2)|.X.COMPL|
|00003c90| 45 58 20 2a 63 31 2c 20 | 2a 63 32 3b 0a 58 7b 0a |EX *c1, |*c2;.X{.|
|00003ca0| 58 09 43 4f 4d 50 4c 45 | 58 20 2a 72 3b 0a 58 09 |X.COMPLE|X *r;.X.|
|00003cb0| 4e 55 4d 42 45 52 20 2a | 71 31 2c 20 2a 71 32 2c |NUMBER *|q1, *q2,|
|00003cc0| 20 2a 71 33 2c 20 2a 64 | 65 6e 3b 0a 58 0a 58 09 | *q3, *d|en;.X.X.|
|00003cd0| 69 66 20 28 63 69 73 7a | 65 72 6f 28 63 32 29 29 |if (cisz|ero(c2))|
|00003ce0| 0a 58 09 09 65 72 72 6f | 72 28 22 44 69 76 69 73 |.X..erro|r("Divis|
|00003cf0| 69 6f 6e 20 62 79 20 7a | 65 72 6f 22 29 3b 0a 58 |ion by z|ero");.X|
|00003d00| 09 69 66 20 28 28 63 31 | 2d 3e 72 65 61 6c 20 3d |.if ((c1|->real =|
|00003d10| 3d 20 63 32 2d 3e 72 65 | 61 6c 29 20 26 26 20 28 |= c2->re|al) && (|
|00003d20| 63 31 2d 3e 69 6d 61 67 | 20 3d 3d 20 63 32 2d 3e |c1->imag| == c2->|
|00003d30| 69 6d 61 67 29 29 0a 58 | 09 09 72 65 74 75 72 6e |imag)).X|..return|
|00003d40| 20 63 6c 69 6e 6b 28 26 | 5f 63 6f 6e 65 5f 29 3b | clink(&|_cone_);|
|00003d50| 0a 58 09 72 20 3d 20 63 | 6f 6d 61 6c 6c 6f 63 28 |.X.r = c|omalloc(|
|00003d60| 29 3b 0a 58 09 69 66 20 | 28 63 69 73 72 65 61 6c |);.X.if |(cisreal|
|00003d70| 28 63 31 29 20 26 26 20 | 63 69 73 72 65 61 6c 28 |(c1) && |cisreal(|
|00003d80| 63 32 29 29 20 7b 0a 58 | 09 09 72 2d 3e 72 65 61 |c2)) {.X|..r->rea|
|00003d90| 6c 20 3d 20 71 64 69 76 | 28 63 31 2d 3e 72 65 61 |l = qdiv|(c1->rea|
|00003da0| 6c 2c 20 63 32 2d 3e 72 | 65 61 6c 29 3b 0a 58 09 |l, c2->r|eal);.X.|
|00003db0| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 09 7d 0a 58 |.return |r;.X.}.X|
|00003dc0| 09 69 66 20 28 63 69 73 | 69 6d 61 67 28 63 31 29 |.if (cis|imag(c1)|
|00003dd0| 20 26 26 20 63 69 73 69 | 6d 61 67 28 63 32 29 29 | && cisi|mag(c2))|
|00003de0| 20 7b 0a 58 09 09 72 2d | 3e 72 65 61 6c 20 3d 20 | {.X..r-|>real = |
|00003df0| 71 64 69 76 28 63 31 2d | 3e 69 6d 61 67 2c 20 63 |qdiv(c1-|>imag, c|
|00003e00| 32 2d 3e 69 6d 61 67 29 | 3b 0a 58 09 09 72 65 74 |2->imag)|;.X..ret|
|00003e10| 75 72 6e 20 72 3b 0a 58 | 09 7d 0a 58 09 69 66 20 |urn r;.X|.}.X.if |
|00003e20| 28 63 69 73 69 6d 61 67 | 28 63 31 29 20 26 26 20 |(cisimag|(c1) && |
|00003e30| 63 69 73 72 65 61 6c 28 | 63 32 29 29 20 7b 0a 58 |cisreal(|c2)) {.X|
|00003e40| 09 09 72 2d 3e 69 6d 61 | 67 20 3d 20 71 64 69 76 |..r->ima|g = qdiv|
|00003e50| 28 63 31 2d 3e 69 6d 61 | 67 2c 20 63 32 2d 3e 72 |(c1->ima|g, c2->r|
|00003e60| 65 61 6c 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 20 |eal);.X.|.return |
|00003e70| 72 3b 0a 58 09 7d 0a 58 | 09 69 66 20 28 63 69 73 |r;.X.}.X|.if (cis|
|00003e80| 72 65 61 6c 28 63 31 29 | 20 26 26 20 63 69 73 69 |real(c1)| && cisi|
|00003e90| 6d 61 67 28 63 32 29 29 | 20 7b 0a 58 09 09 71 31 |mag(c2))| {.X..q1|
|00003ea0| 20 3d 20 71 64 69 76 28 | 63 31 2d 3e 72 65 61 6c | = qdiv(|c1->real|
|00003eb0| 2c 20 63 32 2d 3e 69 6d | 61 67 29 3b 0a 58 09 09 |, c2->im|ag);.X..|
|00003ec0| 72 2d 3e 69 6d 61 67 20 | 3d 20 71 6e 65 67 28 71 |r->imag |= qneg(q|
|00003ed0| 31 29 3b 0a 58 09 09 71 | 66 72 65 65 28 71 31 29 |1);.X..q|free(q1)|
|00003ee0| 3b 0a 58 09 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |;.X..ret|urn r;.X|
|00003ef0| 09 7d 0a 58 09 69 66 20 | 28 63 69 73 72 65 61 6c |.}.X.if |(cisreal|
|00003f00| 28 63 32 29 29 20 7b 0a | 58 09 09 72 2d 3e 72 65 |(c2)) {.|X..r->re|
|00003f10| 61 6c 20 3d 20 71 64 69 | 76 28 63 31 2d 3e 72 65 |al = qdi|v(c1->re|
|00003f20| 61 6c 2c 20 63 32 2d 3e | 72 65 61 6c 29 3b 0a 58 |al, c2->|real);.X|
|00003f30| 09 09 72 2d 3e 69 6d 61 | 67 20 3d 20 71 64 69 76 |..r->ima|g = qdiv|
|00003f40| 28 63 31 2d 3e 69 6d 61 | 67 2c 20 63 32 2d 3e 72 |(c1->ima|g, c2->r|
|00003f50| 65 61 6c 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 20 |eal);.X.|.return |
|00003f60| 72 3b 0a 58 09 7d 0a 58 | 09 71 31 20 3d 20 71 73 |r;.X.}.X|.q1 = qs|
|00003f70| 71 75 61 72 65 28 63 32 | 2d 3e 72 65 61 6c 29 3b |quare(c2|->real);|
|00003f80| 0a 58 09 71 32 20 3d 20 | 71 73 71 75 61 72 65 28 |.X.q2 = |qsquare(|
|00003f90| 63 32 2d 3e 69 6d 61 67 | 29 3b 0a 58 09 64 65 6e |c2->imag|);.X.den|
|00003fa0| 20 3d 20 71 61 64 64 28 | 71 31 2c 20 71 32 29 3b | = qadd(|q1, q2);|
|00003fb0| 0a 58 09 71 66 72 65 65 | 28 71 31 29 3b 0a 58 09 |.X.qfree|(q1);.X.|
|00003fc0| 71 66 72 65 65 28 71 32 | 29 3b 0a 58 09 71 31 20 |qfree(q2|);.X.q1 |
|00003fd0| 3d 20 71 6d 75 6c 28 63 | 31 2d 3e 72 65 61 6c 2c |= qmul(c|1->real,|
|00003fe0| 20 63 32 2d 3e 72 65 61 | 6c 29 3b 0a 58 09 71 32 | c2->rea|l);.X.q2|
|00003ff0| 20 3d 20 71 6d 75 6c 28 | 63 31 2d 3e 69 6d 61 67 | = qmul(|c1->imag|
|00004000| 2c 20 63 32 2d 3e 69 6d | 61 67 29 3b 0a 58 09 71 |, c2->im|ag);.X.q|
|00004010| 33 20 3d 20 71 61 64 64 | 28 71 31 2c 20 71 32 29 |3 = qadd|(q1, q2)|
|00004020| 3b 0a 58 09 71 66 72 65 | 65 28 71 31 29 3b 0a 58 |;.X.qfre|e(q1);.X|
|00004030| 09 71 66 72 65 65 28 71 | 32 29 3b 0a 58 09 72 2d |.qfree(q|2);.X.r-|
|00004040| 3e 72 65 61 6c 20 3d 20 | 71 64 69 76 28 71 33 2c |>real = |qdiv(q3,|
|00004050| 20 64 65 6e 29 3b 0a 58 | 09 71 66 72 65 65 28 71 | den);.X|.qfree(q|
|00004060| 33 29 3b 0a 58 09 71 31 | 20 3d 20 71 6d 75 6c 28 |3);.X.q1| = qmul(|
|00004070| 63 31 2d 3e 72 65 61 6c | 2c 20 63 32 2d 3e 69 6d |c1->real|, c2->im|
|00004080| 61 67 29 3b 0a 58 09 71 | 32 20 3d 20 71 6d 75 6c |ag);.X.q|2 = qmul|
|00004090| 28 63 31 2d 3e 69 6d 61 | 67 2c 20 63 32 2d 3e 72 |(c1->ima|g, c2->r|
|000040a0| 65 61 6c 29 3b 0a 58 09 | 71 33 20 3d 20 71 73 75 |eal);.X.|q3 = qsu|
|000040b0| 62 28 71 32 2c 20 71 31 | 29 3b 0a 58 09 71 66 72 |b(q2, q1|);.X.qfr|
|000040c0| 65 65 28 71 31 29 3b 0a | 58 09 71 66 72 65 65 28 |ee(q1);.|X.qfree(|
|000040d0| 71 32 29 3b 0a 58 09 72 | 2d 3e 69 6d 61 67 20 3d |q2);.X.r|->imag =|
|000040e0| 20 71 64 69 76 28 71 33 | 2c 20 64 65 6e 29 3b 0a | qdiv(q3|, den);.|
|000040f0| 58 09 71 66 72 65 65 28 | 71 33 29 3b 0a 58 09 71 |X.qfree(|q3);.X.q|
|00004100| 66 72 65 65 28 64 65 6e | 29 3b 0a 58 09 72 65 74 |free(den|);.X.ret|
|00004110| 75 72 6e 20 72 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |urn r;.X|}.X.X.X/|
|00004120| 2a 0a 58 20 2a 20 49 6e | 76 65 72 74 20 61 20 63 |*.X * In|vert a c|
|00004130| 6f 6d 70 6c 65 78 20 6e | 75 6d 62 65 72 2e 0a 58 |omplex n|umber..X|
|00004140| 20 2a 2f 0a 58 43 4f 4d | 50 4c 45 58 20 2a 0a 58 | */.XCOM|PLEX *.X|
|00004150| 63 69 6e 76 28 63 29 0a | 58 09 43 4f 4d 50 4c 45 |cinv(c).|X.COMPLE|
|00004160| 58 20 2a 63 3b 0a 58 7b | 0a 58 09 43 4f 4d 50 4c |X *c;.X{|.X.COMPL|
|00004170| 45 58 20 2a 72 3b 0a 58 | 09 4e 55 4d 42 45 52 20 |EX *r;.X|.NUMBER |
|00004180| 2a 71 31 2c 20 2a 71 32 | 2c 20 2a 64 65 6e 3b 0a |*q1, *q2|, *den;.|
|00004190| 58 0a 58 09 69 66 20 28 | 63 69 73 7a 65 72 6f 28 |X.X.if (|ciszero(|
|000041a0| 63 29 29 0a 58 09 09 65 | 72 72 6f 72 28 22 49 6e |c)).X..e|rror("In|
|000041b0| 76 65 72 74 69 6e 67 20 | 7a 65 72 6f 22 29 3b 0a |verting |zero");.|
|000041c0| 58 09 72 20 3d 20 63 6f | 6d 61 6c 6c 6f 63 28 29 |X.r = co|malloc()|
|000041d0| 3b 0a 58 09 69 66 20 28 | 63 69 73 72 65 61 6c 28 |;.X.if (|cisreal(|
|000041e0| 63 29 29 20 7b 0a 58 09 | 09 72 2d 3e 72 65 61 6c |c)) {.X.|.r->real|
|000041f0| 20 3d 20 71 69 6e 76 28 | 63 2d 3e 72 65 61 6c 29 | = qinv(|c->real)|
|00004200| 3b 0a 58 09 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |;.X..ret|urn r;.X|
|00004210| 09 7d 0a 58 09 69 66 20 | 28 63 69 73 69 6d 61 67 |.}.X.if |(cisimag|
|00004220| 28 63 29 29 20 7b 0a 58 | 09 09 71 31 20 3d 20 71 |(c)) {.X|..q1 = q|
|00004230| 69 6e 76 28 63 2d 3e 69 | 6d 61 67 29 3b 0a 58 09 |inv(c->i|mag);.X.|
|00004240| 09 72 2d 3e 69 6d 61 67 | 20 3d 20 71 6e 65 67 28 |.r->imag| = qneg(|
|00004250| 71 31 29 3b 0a 58 09 09 | 71 66 72 65 65 28 71 31 |q1);.X..|qfree(q1|
|00004260| 29 3b 0a 58 09 09 72 65 | 74 75 72 6e 20 72 3b 0a |);.X..re|turn r;.|
|00004270| 58 09 7d 0a 58 09 71 31 | 20 3d 20 71 73 71 75 61 |X.}.X.q1| = qsqua|
|00004280| 72 65 28 63 2d 3e 72 65 | 61 6c 29 3b 0a 58 09 71 |re(c->re|al);.X.q|
|00004290| 32 20 3d 20 71 73 71 75 | 61 72 65 28 63 2d 3e 69 |2 = qsqu|are(c->i|
|000042a0| 6d 61 67 29 3b 0a 58 09 | 64 65 6e 20 3d 20 71 61 |mag);.X.|den = qa|
|000042b0| 64 64 28 71 31 2c 20 71 | 32 29 3b 0a 58 09 71 66 |dd(q1, q|2);.X.qf|
|000042c0| 72 65 65 28 71 31 29 3b | 0a 58 09 71 66 72 65 65 |ree(q1);|.X.qfree|
|000042d0| 28 71 32 29 3b 0a 58 09 | 72 2d 3e 72 65 61 6c 20 |(q2);.X.|r->real |
|000042e0| 3d 20 71 64 69 76 28 63 | 2d 3e 72 65 61 6c 2c 20 |= qdiv(c|->real, |
|000042f0| 64 65 6e 29 3b 0a 58 09 | 71 31 20 3d 20 71 64 69 |den);.X.|q1 = qdi|
|00004300| 76 28 63 2d 3e 69 6d 61 | 67 2c 20 64 65 6e 29 3b |v(c->ima|g, den);|
|00004310| 0a 58 09 72 2d 3e 69 6d | 61 67 20 3d 20 71 6e 65 |.X.r->im|ag = qne|
|00004320| 67 28 71 31 29 3b 0a 58 | 09 71 66 72 65 65 28 71 |g(q1);.X|.qfree(q|
|00004330| 31 29 3b 0a 58 09 71 66 | 72 65 65 28 64 65 6e 29 |1);.X.qf|ree(den)|
|00004340| 3b 0a 58 09 72 65 74 75 | 72 6e 20 72 3b 0a 58 7d |;.X.retu|rn r;.X}|
|00004350| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 4e 65 67 |.X.X.X/*|.X * Neg|
|00004360| 61 74 65 20 61 20 63 6f | 6d 70 6c 65 78 20 6e 75 |ate a co|mplex nu|
|00004370| 6d 62 65 72 2e 0a 58 20 | 2a 2f 0a 58 43 4f 4d 50 |mber..X |*/.XCOMP|
|00004380| 4c 45 58 20 2a 0a 58 63 | 6e 65 67 28 63 29 0a 58 |LEX *.Xc|neg(c).X|
|00004390| 09 43 4f 4d 50 4c 45 58 | 20 2a 63 3b 0a 58 7b 0a |.COMPLEX| *c;.X{.|
|000043a0| 58 09 43 4f 4d 50 4c 45 | 58 20 2a 72 3b 0a 58 0a |X.COMPLE|X *r;.X.|
|000043b0| 58 09 69 66 20 28 63 69 | 73 7a 65 72 6f 28 63 29 |X.if (ci|szero(c)|
|000043c0| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 63 6c 69 6e |).X..ret|urn clin|
|000043d0| 6b 28 26 5f 63 7a 65 72 | 6f 5f 29 3b 0a 58 09 72 |k(&_czer|o_);.X.r|
|000043e0| 20 3d 20 63 6f 6d 61 6c | 6c 6f 63 28 29 3b 0a 58 | = comal|loc();.X|
|000043f0| 09 69 66 20 28 21 71 69 | 73 7a 65 72 6f 28 63 2d |.if (!qi|szero(c-|
|00004400| 3e 72 65 61 6c 29 29 0a | 58 09 09 72 2d 3e 72 65 |>real)).|X..r->re|
|00004410| 61 6c 20 3d 20 71 6e 65 | 67 28 63 2d 3e 72 65 61 |al = qne|g(c->rea|
|00004420| 6c 29 3b 0a 58 09 69 66 | 20 28 21 71 69 73 7a 65 |l);.X.if| (!qisze|
|00004430| 72 6f 28 63 2d 3e 69 6d | 61 67 29 29 0a 58 09 09 |ro(c->im|ag)).X..|
|00004440| 72 2d 3e 69 6d 61 67 20 | 3d 20 71 6e 65 67 28 63 |r->imag |= qneg(c|
|00004450| 2d 3e 69 6d 61 67 29 3b | 0a 58 09 72 65 74 75 72 |->imag);|.X.retur|
|00004460| 6e 20 72 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 0a |n r;.X}.|X.X.X/*.|
|00004470| 58 20 2a 20 54 61 6b 65 | 20 74 68 65 20 69 6e 74 |X * Take| the int|
|00004480| 65 67 65 72 20 70 61 72 | 74 20 6f 66 20 61 20 63 |eger par|t of a c|
|00004490| 6f 6d 70 6c 65 78 20 6e | 75 6d 62 65 72 2e 0a 58 |omplex n|umber..X|
|000044a0| 20 2a 20 54 68 69 73 20 | 6d 65 61 6e 73 20 74 61 | * This |means ta|
|000044b0| 6b 65 20 74 68 65 20 69 | 6e 74 65 67 65 72 20 70 |ke the i|nteger p|
|000044c0| 61 72 74 20 6f 66 20 62 | 6f 74 68 20 63 6f 6d 70 |art of b|oth comp|
|000044d0| 6f 6e 65 6e 74 73 2e 0a | 58 20 2a 2f 0a 58 43 4f |onents..|X */.XCO|
|000044e0| 4d 50 4c 45 58 20 2a 0a | 58 63 69 6e 74 28 63 29 |MPLEX *.|Xcint(c)|
|000044f0| 0a 58 09 43 4f 4d 50 4c | 45 58 20 2a 63 3b 0a 58 |.X.COMPL|EX *c;.X|
|00004500| 7b 0a 58 09 43 4f 4d 50 | 4c 45 58 20 2a 72 3b 0a |{.X.COMP|LEX *r;.|
|00004510| 58 0a 58 09 69 66 20 28 | 63 69 73 69 6e 74 28 63 |X.X.if (|cisint(c|
|00004520| 29 29 0a 58 09 09 72 65 | 74 75 72 6e 20 63 6c 69 |)).X..re|turn cli|
|00004530| 6e 6b 28 63 29 3b 0a 58 | 09 72 20 3d 20 63 6f 6d |nk(c);.X|.r = com|
|00004540| 61 6c 6c 6f 63 28 29 3b | 0a 58 09 72 2d 3e 72 65 |alloc();|.X.r->re|
|00004550| 61 6c 20 3d 20 71 69 6e | 74 28 63 2d 3e 72 65 61 |al = qin|t(c->rea|
|00004560| 6c 29 3b 0a 58 09 72 2d | 3e 69 6d 61 67 20 3d 20 |l);.X.r-|>imag = |
|00004570| 71 69 6e 74 28 63 2d 3e | 69 6d 61 67 29 3b 0a 58 |qint(c->|imag);.X|
|00004580| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 7d 0a 58 0a |.return |r;.X}.X.|
|00004590| 58 0a 58 2f 2a 0a 58 20 | 2a 20 54 61 6b 65 20 74 |X.X/*.X |* Take t|
|000045a0| 68 65 20 66 72 61 63 74 | 69 6f 6e 61 6c 20 70 61 |he fract|ional pa|
|000045b0| 72 74 20 6f 66 20 61 20 | 63 6f 6d 70 6c 65 78 20 |rt of a |complex |
|000045c0| 6e 75 6d 62 65 72 2e 0a | 58 20 2a 20 54 68 69 73 |number..|X * This|
|000045d0| 20 6d 65 61 6e 73 20 74 | 61 6b 65 20 74 68 65 20 | means t|ake the |
|000045e0| 66 72 61 63 74 69 6f 6e | 61 6c 20 70 61 72 74 20 |fraction|al part |
|000045f0| 6f 66 20 62 6f 74 68 20 | 63 6f 6d 70 6f 6e 65 6e |of both |componen|
|00004600| 74 73 2e 0a 58 20 2a 2f | 0a 58 43 4f 4d 50 4c 45 |ts..X */|.XCOMPLE|
|00004610| 58 20 2a 0a 58 63 66 72 | 61 63 28 63 29 0a 58 09 |X *.Xcfr|ac(c).X.|
|00004620| 43 4f 4d 50 4c 45 58 20 | 2a 63 3b 0a 58 7b 0a 58 |COMPLEX |*c;.X{.X|
|00004630| 09 43 4f 4d 50 4c 45 58 | 20 2a 72 3b 0a 58 0a 58 |.COMPLEX| *r;.X.X|
|00004640| 09 69 66 20 28 63 69 73 | 69 6e 74 28 63 29 29 0a |.if (cis|int(c)).|
|00004650| 58 09 09 72 65 74 75 72 | 6e 20 63 6c 69 6e 6b 28 |X..retur|n clink(|
|00004660| 26 5f 63 7a 65 72 6f 5f | 29 3b 0a 58 09 72 20 3d |&_czero_|);.X.r =|
|00004670| 20 63 6f 6d 61 6c 6c 6f | 63 28 29 3b 0a 58 09 72 | comallo|c();.X.r|
|00004680| 2d 3e 72 65 61 6c 20 3d | 20 71 66 72 61 63 28 63 |->real =| qfrac(c|
|00004690| 2d 3e 72 65 61 6c 29 3b | 0a 58 09 72 2d 3e 69 6d |->real);|.X.r->im|
|000046a0| 61 67 20 3d 20 71 66 72 | 61 63 28 63 2d 3e 69 6d |ag = qfr|ac(c->im|
|000046b0| 61 67 29 3b 0a 58 09 72 | 65 74 75 72 6e 20 72 3b |ag);.X.r|eturn r;|
|000046c0| 0a 58 7d 0a 58 0a 58 0a | 58 23 69 66 20 30 0a 58 |.X}.X.X.|X#if 0.X|
|000046d0| 2f 2a 0a 58 20 2a 20 54 | 61 6b 65 20 74 68 65 20 |/*.X * T|ake the |
|000046e0| 63 6f 6e 6a 75 67 61 74 | 65 20 6f 66 20 61 20 63 |conjugat|e of a c|
|000046f0| 6f 6d 70 6c 65 78 20 6e | 75 6d 62 65 72 2e 0a 58 |omplex n|umber..X|
|00004700| 20 2a 20 54 68 69 73 20 | 6e 65 67 61 74 65 73 20 | * This |negates |
|00004710| 74 68 65 20 63 6f 6d 70 | 6c 65 78 20 70 61 72 74 |the comp|lex part|
|00004720| 2e 0a 58 20 2a 2f 0a 58 | 43 4f 4d 50 4c 45 58 20 |..X */.X|COMPLEX |
|00004730| 2a 0a 58 63 63 6f 6e 6a | 28 63 29 0a 58 09 43 4f |*.Xcconj|(c).X.CO|
|00004740| 4d 50 4c 45 58 20 2a 63 | 3b 0a 58 7b 0a 58 09 43 |MPLEX *c|;.X{.X.C|
|00004750| 4f 4d 50 4c 45 58 20 2a | 72 3b 0a 58 0a 58 09 69 |OMPLEX *|r;.X.X.i|
|00004760| 66 20 28 63 69 73 72 65 | 61 6c 28 63 29 29 0a 58 |f (cisre|al(c)).X|
|00004770| 09 09 72 65 74 75 72 6e | 20 63 6c 69 6e 6b 28 63 |..return| clink(c|
|00004780| 29 3b 0a 58 09 72 20 3d | 20 63 6f 6d 61 6c 6c 6f |);.X.r =| comallo|
|00004790| 63 28 29 3b 0a 58 09 69 | 66 20 28 21 71 69 73 7a |c();.X.i|f (!qisz|
|000047a0| 65 72 6f 28 63 2d 3e 72 | 65 61 6c 29 29 0a 58 09 |ero(c->r|eal)).X.|
|000047b0| 09 72 2d 3e 72 65 61 6c | 20 3d 20 71 6c 69 6e 6b |.r->real| = qlink|
|000047c0| 28 63 2d 3e 72 65 61 6c | 29 3b 0a 58 09 72 2d 3e |(c->real|);.X.r->|
|000047d0| 69 6d 61 67 20 3d 20 71 | 6e 65 67 28 63 2d 3e 69 |imag = q|neg(c->i|
|000047e0| 6d 61 67 29 3b 0a 58 09 | 72 65 74 75 72 6e 20 72 |mag);.X.|return r|
|000047f0| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00004800| 20 52 65 74 75 72 6e 20 | 74 68 65 20 72 65 61 6c | Return |the real|
|00004810| 20 70 61 72 74 20 6f 66 | 20 61 20 63 6f 6d 70 6c | part of| a compl|
|00004820| 65 78 20 6e 75 6d 62 65 | 72 2e 0a 58 20 2a 2f 0a |ex numbe|r..X */.|
|00004830| 58 43 4f 4d 50 4c 45 58 | 20 2a 0a 58 63 72 65 61 |XCOMPLEX| *.Xcrea|
|00004840| 6c 28 63 29 0a 58 09 43 | 4f 4d 50 4c 45 58 20 2a |l(c).X.C|OMPLEX *|
|00004850| 63 3b 0a 58 7b 0a 58 09 | 43 4f 4d 50 4c 45 58 20 |c;.X{.X.|COMPLEX |
|00004860| 2a 72 3b 0a 58 0a 58 09 | 69 66 20 28 63 69 73 72 |*r;.X.X.|if (cisr|
|00004870| 65 61 6c 28 63 29 29 0a | 58 09 09 72 65 74 75 72 |eal(c)).|X..retur|
|00004880| 6e 20 63 6c 69 6e 6b 28 | 63 29 3b 0a 58 09 72 20 |n clink(|c);.X.r |
|00004890| 3d 20 63 6f 6d 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |= comall|oc();.X.|
|000048a0| 69 66 20 28 21 71 69 73 | 7a 65 72 6f 28 63 2d 3e |if (!qis|zero(c->|
|000048b0| 72 65 61 6c 29 29 0a 58 | 09 09 72 2d 3e 72 65 61 |real)).X|..r->rea|
|000048c0| 6c 20 3d 20 71 6c 69 6e | 6b 28 63 2d 3e 72 65 61 |l = qlin|k(c->rea|
|000048d0| 6c 29 3b 0a 58 09 72 65 | 74 75 72 6e 20 72 3b 0a |l);.X.re|turn r;.|
|000048e0| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 52 |X}.X.X.X|/*.X * R|
|000048f0| 65 74 75 72 6e 20 74 68 | 65 20 69 6d 61 67 69 6e |eturn th|e imagin|
|00004900| 61 72 79 20 70 61 72 74 | 20 6f 66 20 61 20 63 6f |ary part| of a co|
|00004910| 6d 70 6c 65 78 20 6e 75 | 6d 62 65 72 20 61 73 20 |mplex nu|mber as |
|00004920| 61 20 72 65 61 6c 2e 0a | 58 20 2a 2f 0a 58 43 4f |a real..|X */.XCO|
|00004930| 4d 50 4c 45 58 20 2a 0a | 58 63 69 6d 61 67 28 63 |MPLEX *.|Xcimag(c|
|00004940| 29 0a 58 09 43 4f 4d 50 | 4c 45 58 20 2a 63 3b 0a |).X.COMP|LEX *c;.|
|00004950| 58 7b 0a 58 09 43 4f 4d | 50 4c 45 58 20 2a 72 3b |X{.X.COM|PLEX *r;|
|00004960| 0a 58 0a 58 09 69 66 20 | 28 63 69 73 72 65 61 6c |.X.X.if |(cisreal|
|00004970| 28 63 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 63 |(c)).X..|return c|
|00004980| 6c 69 6e 6b 28 26 5f 63 | 7a 65 72 6f 5f 29 3b 0a |link(&_c|zero_);.|
|00004990| 58 09 72 20 3d 20 63 6f | 6d 61 6c 6c 6f 63 28 29 |X.r = co|malloc()|
|000049a0| 3b 0a 58 09 72 2d 3e 72 | 65 61 6c 20 3d 20 71 6c |;.X.r->r|eal = ql|
|000049b0| 69 6e 6b 28 63 2d 3e 69 | 6d 61 67 29 3b 0a 58 09 |ink(c->i|mag);.X.|
|000049c0| 72 65 74 75 72 6e 20 72 | 3b 0a 58 7d 0a 58 23 65 |return r|;.X}.X#e|
|000049d0| 6e 64 69 66 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |ndif.X.X|.X/*.X *|
|000049e0| 20 41 64 64 20 61 20 72 | 65 61 6c 20 6e 75 6d 62 | Add a r|eal numb|
|000049f0| 65 72 20 74 6f 20 61 20 | 63 6f 6d 70 6c 65 78 20 |er to a |complex |
|00004a00| 6e 75 6d 62 65 72 2e 0a | 58 20 2a 2f 0a 58 43 4f |number..|X */.XCO|
|00004a10| 4d 50 4c 45 58 20 2a 0a | 58 63 61 64 64 71 28 63 |MPLEX *.|Xcaddq(c|
|00004a20| 2c 20 71 29 0a 58 09 43 | 4f 4d 50 4c 45 58 20 2a |, q).X.C|OMPLEX *|
|00004a30| 63 3b 0a 58 09 4e 55 4d | 42 45 52 20 2a 71 3b 0a |c;.X.NUM|BER *q;.|
|00004a40| 58 7b 0a 58 09 43 4f 4d | 50 4c 45 58 20 2a 72 3b |X{.X.COM|PLEX *r;|
|00004a50| 0a 58 0a 58 09 69 66 20 | 28 71 69 73 7a 65 72 6f |.X.X.if |(qiszero|
|00004a60| 28 71 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 63 |(q)).X..|return c|
|00004a70| 6c 69 6e 6b 28 63 29 3b | 0a 58 09 72 20 3d 20 63 |link(c);|.X.r = c|
|00004a80| 6f 6d 61 6c 6c 6f 63 28 | 29 3b 0a 58 09 72 2d 3e |omalloc(|);.X.r->|
|00004a90| 72 65 61 6c 20 3d 20 71 | 61 64 64 28 63 2d 3e 72 |real = q|add(c->r|
|00004aa0| 65 61 6c 2c 20 71 29 3b | 0a 58 09 72 2d 3e 69 6d |eal, q);|.X.r->im|
|00004ab0| 61 67 20 3d 20 71 6c 69 | 6e 6b 28 63 2d 3e 69 6d |ag = qli|nk(c->im|
|00004ac0| 61 67 29 3b 0a 58 09 72 | 65 74 75 72 6e 20 72 3b |ag);.X.r|eturn r;|
|00004ad0| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |.X}.X.X.|X/*.X * |
|00004ae0| 53 75 62 74 72 61 63 74 | 20 61 20 72 65 61 6c 20 |Subtract| a real |
|00004af0| 6e 75 6d 62 65 72 20 66 | 72 6f 6d 20 61 20 63 6f |number f|rom a co|
|00004b00| 6d 70 6c 65 78 20 6e 75 | 6d 62 65 72 2e 0a 58 20 |mplex nu|mber..X |
|00004b10| 2a 2f 0a 58 43 4f 4d 50 | 4c 45 58 20 2a 0a 58 63 |*/.XCOMP|LEX *.Xc|
|00004b20| 73 75 62 71 28 63 2c 20 | 71 29 0a 58 09 43 4f 4d |subq(c, |q).X.COM|
|00004b30| 50 4c 45 58 20 2a 63 3b | 0a 58 09 4e 55 4d 42 45 |PLEX *c;|.X.NUMBE|
|00004b40| 52 20 2a 71 3b 0a 58 7b | 0a 58 09 43 4f 4d 50 4c |R *q;.X{|.X.COMPL|
|00004b50| 45 58 20 2a 72 3b 0a 58 | 0a 58 09 69 66 20 28 71 |EX *r;.X|.X.if (q|
|00004b60| 69 73 7a 65 72 6f 28 71 | 29 29 0a 58 09 09 72 65 |iszero(q|)).X..re|
|00004b70| 74 75 72 6e 20 63 6c 69 | 6e 6b 28 63 29 3b 0a 58 |turn cli|nk(c);.X|
|00004b80| 09 72 20 3d 20 63 6f 6d | 61 6c 6c 6f 63 28 29 3b |.r = com|alloc();|
|00004b90| 0a 58 09 72 2d 3e 72 65 | 61 6c 20 3d 20 71 73 75 |.X.r->re|al = qsu|
|00004ba0| 62 28 63 2d 3e 72 65 61 | 6c 2c 20 71 29 3b 0a 58 |b(c->rea|l, q);.X|
|00004bb0| 09 72 2d 3e 69 6d 61 67 | 20 3d 20 71 6c 69 6e 6b |.r->imag| = qlink|
|00004bc0| 28 63 2d 3e 69 6d 61 67 | 29 3b 0a 58 09 72 65 74 |(c->imag|);.X.ret|
|00004bd0| 75 72 6e 20 72 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |urn r;.X|}.X.X.X/|
|00004be0| 2a 0a 58 20 2a 20 53 68 | 69 66 74 20 74 68 65 20 |*.X * Sh|ift the |
|00004bf0| 63 6f 6d 70 6f 6e 65 6e | 74 73 20 6f 66 20 61 20 |componen|ts of a |
|00004c00| 63 6f 6d 70 6c 65 78 20 | 6e 75 6d 62 65 72 20 6c |complex |number l|
|00004c10| 65 66 74 20 62 79 20 74 | 68 65 20 73 70 65 63 69 |eft by t|he speci|
|00004c20| 66 69 65 64 0a 58 20 2a | 20 6e 75 6d 62 65 72 20 |fied.X *| number |
|00004c30| 6f 66 20 62 69 74 73 2e | 20 20 4e 65 67 61 74 69 |of bits.|  Negati|
|00004c40| 76 65 20 76 61 6c 75 65 | 73 20 73 68 69 66 74 20 |ve value|s shift |
|00004c50| 74 6f 20 74 68 65 20 72 | 69 67 68 74 2e 0a 58 20 |to the r|ight..X |
|00004c60| 2a 2f 0a 58 43 4f 4d 50 | 4c 45 58 20 2a 0a 58 63 |*/.XCOMP|LEX *.Xc|
|00004c70| 73 68 69 66 74 28 63 2c | 20 6e 29 0a 58 09 43 4f |shift(c,| n).X.CO|
|00004c80| 4d 50 4c 45 58 20 2a 63 | 3b 0a 58 09 6c 6f 6e 67 |MPLEX *c|;.X.long|
|00004c90| 20 6e 3b 0a 58 7b 0a 58 | 09 43 4f 4d 50 4c 45 58 | n;.X{.X|.COMPLEX|
|00004ca0| 20 2a 72 3b 0a 58 0a 58 | 09 69 66 20 28 63 69 73 | *r;.X.X|.if (cis|
|00004cb0| 7a 65 72 6f 28 63 29 20 | 7c 7c 20 28 6e 20 3d 3d |zero(c) ||| (n ==|
|00004cc0| 20 30 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 63 | 0)).X..|return c|
|00004cd0| 6c 69 6e 6b 28 63 29 3b | 0a 58 09 72 20 3d 20 63 |link(c);|.X.r = c|
|00004ce0| 6f 6d 61 6c 6c 6f 63 28 | 29 3b 0a 58 09 72 2d 3e |omalloc(|);.X.r->|
|00004cf0| 72 65 61 6c 20 3d 20 71 | 73 68 69 66 74 28 63 2d |real = q|shift(c-|
|00004d00| 3e 72 65 61 6c 2c 20 6e | 29 3b 0a 58 09 72 2d 3e |>real, n|);.X.r->|
|00004d10| 69 6d 61 67 20 3d 20 71 | 73 68 69 66 74 28 63 2d |imag = q|shift(c-|
|00004d20| 3e 69 6d 61 67 2c 20 6e | 29 3b 0a 58 09 72 65 74 |>imag, n|);.X.ret|
|00004d30| 75 72 6e 20 72 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |urn r;.X|}.X.X.X/|
|00004d40| 2a 0a 58 20 2a 20 53 63 | 61 6c 65 20 61 20 63 6f |*.X * Sc|ale a co|
|00004d50| 6d 70 6c 65 78 20 6e 75 | 6d 62 65 72 20 62 79 20 |mplex nu|mber by |
|00004d60| 61 20 70 6f 77 65 72 20 | 6f 66 20 74 77 6f 2e 0a |a power |of two..|
|00004d70| 58 20 2a 2f 0a 58 43 4f | 4d 50 4c 45 58 20 2a 0a |X */.XCO|MPLEX *.|
|00004d80| 58 63 73 63 61 6c 65 28 | 63 2c 20 6e 29 0a 58 09 |Xcscale(|c, n).X.|
|00004d90| 43 4f 4d 50 4c 45 58 20 | 2a 63 3b 0a 58 09 6c 6f |COMPLEX |*c;.X.lo|
|00004da0| 6e 67 20 6e 3b 0a 58 7b | 0a 58 09 43 4f 4d 50 4c |ng n;.X{|.X.COMPL|
|00004db0| 45 58 20 2a 72 3b 0a 58 | 0a 58 09 69 66 20 28 63 |EX *r;.X|.X.if (c|
|00004dc0| 69 73 7a 65 72 6f 28 63 | 29 20 7c 7c 20 28 6e 20 |iszero(c|) || (n |
|00004dd0| 3d 3d 20 30 29 29 0a 58 | 09 09 72 65 74 75 72 6e |== 0)).X|..return|
|00004de0| 20 63 6c 69 6e 6b 28 63 | 29 3b 0a 58 09 72 20 3d | clink(c|);.X.r =|
|00004df0| 20 63 6f 6d 61 6c 6c 6f | 63 28 29 3b 0a 58 09 72 | comallo|c();.X.r|
|00004e00| 2d 3e 72 65 61 6c 20 3d | 20 71 73 63 61 6c 65 28 |->real =| qscale(|
|00004e10| 63 2d 3e 72 65 61 6c 2c | 20 6e 29 3b 0a 58 09 72 |c->real,| n);.X.r|
|00004e20| 2d 3e 69 6d 61 67 20 3d | 20 71 73 63 61 6c 65 28 |->imag =| qscale(|
|00004e30| 63 2d 3e 69 6d 61 67 2c | 20 6e 29 3b 0a 58 09 72 |c->imag,| n);.X.r|
|00004e40| 65 74 75 72 6e 20 72 3b | 0a 58 7d 0a 58 0a 58 0a |eturn r;|.X}.X.X.|
|00004e50| 58 2f 2a 0a 58 20 2a 20 | 4d 75 6c 74 69 70 6c 79 |X/*.X * |Multiply|
|00004e60| 20 61 20 63 6f 6d 70 6c | 65 78 20 6e 75 6d 62 65 | a compl|ex numbe|
|00004e70| 72 20 62 79 20 61 20 72 | 65 61 6c 20 6e 75 6d 62 |r by a r|eal numb|
|00004e80| 65 72 2e 0a 58 20 2a 2f | 0a 58 43 4f 4d 50 4c 45 |er..X */|.XCOMPLE|
|00004e90| 58 20 2a 0a 58 63 6d 75 | 6c 71 28 63 2c 20 71 29 |X *.Xcmu|lq(c, q)|
|00004ea0| 0a 58 09 43 4f 4d 50 4c | 45 58 20 2a 63 3b 0a 58 |.X.COMPL|EX *c;.X|
|00004eb0| 09 4e 55 4d 42 45 52 20 | 2a 71 3b 0a 58 7b 0a 58 |.NUMBER |*q;.X{.X|
|00004ec0| 09 43 4f 4d 50 4c 45 58 | 20 2a 72 3b 0a 58 0a 58 |.COMPLEX| *r;.X.X|
|00004ed0| 09 69 66 20 28 71 69 73 | 7a 65 72 6f 28 71 29 29 |.if (qis|zero(q))|
|00004ee0| 0a 58 09 09 72 65 74 75 | 72 6e 20 63 6c 69 6e 6b |.X..retu|rn clink|
|00004ef0| 28 26 5f 63 7a 65 72 6f | 5f 29 3b 0a 58 09 69 66 |(&_czero|_);.X.if|
|00004f00| 20 28 71 69 73 6f 6e 65 | 28 71 29 29 0a 58 09 09 | (qisone|(q)).X..|
|00004f10| 72 65 74 75 72 6e 20 63 | 6c 69 6e 6b 28 63 29 3b |return c|link(c);|
|00004f20| 0a 58 09 69 66 20 28 71 | 69 73 6e 65 67 6f 6e 65 |.X.if (q|isnegone|
|00004f30| 28 71 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 63 |(q)).X..|return c|
|00004f40| 6e 65 67 28 63 29 3b 0a | 58 09 72 20 3d 20 63 6f |neg(c);.|X.r = co|
|00004f50| 6d 61 6c 6c 6f 63 28 29 | 3b 0a 58 09 72 2d 3e 72 |malloc()|;.X.r->r|
|00004f60| 65 61 6c 20 3d 20 71 6d | 75 6c 28 63 2d 3e 72 65 |eal = qm|ul(c->re|
|00004f70| 61 6c 2c 20 71 29 3b 0a | 58 09 72 2d 3e 69 6d 61 |al, q);.|X.r->ima|
|00004f80| 67 20 3d 20 71 6d 75 6c | 28 63 2d 3e 69 6d 61 67 |g = qmul|(c->imag|
|00004f90| 2c 20 71 29 3b 0a 58 09 | 72 65 74 75 72 6e 20 72 |, q);.X.|return r|
|00004fa0| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00004fb0| 20 44 69 76 69 64 65 20 | 61 20 63 6f 6d 70 6c 65 | Divide |a comple|
|00004fc0| 78 20 6e 75 6d 62 65 72 | 20 62 79 20 61 20 72 65 |x number| by a re|
|00004fd0| 61 6c 20 6e 75 6d 62 65 | 72 2e 0a 58 20 2a 2f 0a |al numbe|r..X */.|
|00004fe0| 58 43 4f 4d 50 4c 45 58 | 20 2a 0a 58 63 64 69 76 |XCOMPLEX| *.Xcdiv|
|00004ff0| 71 28 63 2c 20 71 29 0a | 58 09 43 4f 4d 50 4c 45 |q(c, q).|X.COMPLE|
|00005000| 58 20 2a 63 3b 0a 58 09 | 4e 55 4d 42 45 52 20 2a |X *c;.X.|NUMBER *|
|00005010| 71 3b 0a 58 7b 0a 58 09 | 43 4f 4d 50 4c 45 58 20 |q;.X{.X.|COMPLEX |
|00005020| 2a 72 3b 0a 58 0a 58 09 | 69 66 20 28 71 69 73 7a |*r;.X.X.|if (qisz|
|00005030| 65 72 6f 28 71 29 29 0a | 58 09 09 65 72 72 6f 72 |ero(q)).|X..error|
|00005040| 28 22 44 69 76 69 73 69 | 6f 6e 20 62 79 20 7a 65 |("Divisi|on by ze|
|00005050| 72 6f 22 29 3b 0a 58 09 | 69 66 20 28 71 69 73 6f |ro");.X.|if (qiso|
|00005060| 6e 65 28 71 29 29 0a 58 | 09 09 72 65 74 75 72 6e |ne(q)).X|..return|
|00005070| 20 63 6c 69 6e 6b 28 63 | 29 3b 0a 58 09 69 66 20 | clink(c|);.X.if |
|00005080| 28 71 69 73 6e 65 67 6f | 6e 65 28 71 29 29 0a 58 |(qisnego|ne(q)).X|
|00005090| 09 09 72 65 74 75 72 6e | 20 63 6e 65 67 28 63 29 |..return| cneg(c)|
|000050a0| 3b 0a 58 09 72 20 3d 20 | 63 6f 6d 61 6c 6c 6f 63 |;.X.r = |comalloc|
|000050b0| 28 29 3b 0a 58 09 72 2d | 3e 72 65 61 6c 20 3d 20 |();.X.r-|>real = |
|000050c0| 71 64 69 76 28 63 2d 3e | 72 65 61 6c 2c 20 71 29 |qdiv(c->|real, q)|
|000050d0| 3b 0a 58 09 72 2d 3e 69 | 6d 61 67 20 3d 20 71 64 |;.X.r->i|mag = qd|
|000050e0| 69 76 28 63 2d 3e 69 6d | 61 67 2c 20 71 29 3b 0a |iv(c->im|ag, q);.|
|000050f0| 58 09 72 65 74 75 72 6e | 20 72 3b 0a 58 7d 0a 58 |X.return| r;.X}.X|
|00005100| 0a 58 0a 58 2f 2a 0a 58 | 20 2a 20 54 61 6b 65 20 |.X.X/*.X| * Take |
|00005110| 74 68 65 20 69 6e 74 65 | 67 65 72 20 71 75 6f 74 |the inte|ger quot|
|00005120| 69 65 6e 74 20 6f 66 20 | 61 20 63 6f 6d 70 6c 65 |ient of |a comple|
|00005130| 78 20 6e 75 6d 62 65 72 | 20 62 79 20 61 20 72 65 |x number| by a re|
|00005140| 61 6c 20 6e 75 6d 62 65 | 72 2e 0a 58 20 2a 20 54 |al numbe|r..X * T|
|00005150| 68 69 73 20 69 73 20 64 | 65 66 69 6e 65 64 20 74 |his is d|efined t|
|00005160| 6f 20 62 65 20 74 68 65 | 20 72 65 73 75 6c 74 20 |o be the| result |
|00005170| 6f 66 20 64 6f 69 6e 67 | 20 74 68 65 20 71 75 6f |of doing| the quo|
|00005180| 74 69 65 6e 74 20 66 6f | 72 20 65 61 63 68 20 63 |tient fo|r each c|
|00005190| 6f 6d 70 6f 6e 65 6e 74 | 2e 0a 58 20 2a 2f 0a 58 |omponent|..X */.X|
|000051a0| 43 4f 4d 50 4c 45 58 20 | 2a 0a 58 63 71 75 6f 71 |COMPLEX |*.Xcquoq|
|000051b0| 28 63 2c 20 71 29 0a 58 | 09 43 4f 4d 50 4c 45 58 |(c, q).X|.COMPLEX|
|000051c0| 20 2a 63 3b 0a 58 09 4e | 55 4d 42 45 52 20 2a 71 | *c;.X.N|UMBER *q|
|000051d0| 3b 0a 58 7b 0a 58 09 43 | 4f 4d 50 4c 45 58 20 2a |;.X{.X.C|OMPLEX *|
|000051e0| 72 3b 0a 58 0a 58 09 69 | 66 20 28 71 69 73 7a 65 |r;.X.X.i|f (qisze|
|000051f0| 72 6f 28 71 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |ro(q)).X|..error(|
|00005200| 22 44 69 76 69 73 69 6f | 6e 20 62 79 20 7a 65 72 |"Divisio|n by zer|
|00005210| 6f 22 29 3b 0a 58 09 72 | 20 3d 20 63 6f 6d 61 6c |o");.X.r| = comal|
|00005220| 6c 6f 63 28 29 3b 0a 58 | 09 72 2d 3e 72 65 61 6c |loc();.X|.r->real|
|00005230| 20 3d 20 71 71 75 6f 28 | 63 2d 3e 72 65 61 6c 2c | = qquo(|c->real,|
|00005240| 20 71 29 3b 0a 58 09 72 | 2d 3e 69 6d 61 67 20 3d | q);.X.r|->imag =|
|00005250| 20 71 71 75 6f 28 63 2d | 3e 69 6d 61 67 2c 20 71 | qquo(c-|>imag, q|
|00005260| 29 3b 0a 58 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |);.X.ret|urn r;.X|
|00005270| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 54 61 |}.X.X.X/|*.X * Ta|
|00005280| 6b 65 20 74 68 65 20 6d | 6f 64 75 6c 75 73 20 6f |ke the m|odulus o|
|00005290| 66 20 61 20 63 6f 6d 70 | 6c 65 78 20 6e 75 6d 62 |f a comp|lex numb|
|000052a0| 65 72 20 62 79 20 61 20 | 72 65 61 6c 20 6e 75 6d |er by a |real num|
|000052b0| 62 65 72 2e 0a 58 20 2a | 20 54 68 69 73 20 69 73 |ber..X *| This is|
|000052c0| 20 64 65 66 69 6e 65 64 | 20 74 6f 20 62 65 20 74 | defined| to be t|
|000052d0| 68 65 20 72 65 73 75 6c | 74 20 6f 66 20 64 6f 69 |he resul|t of doi|
|000052e0| 6e 67 20 74 68 65 20 6d | 6f 64 75 6c 6f 20 66 6f |ng the m|odulo fo|
|000052f0| 72 20 65 61 63 68 20 63 | 6f 6d 70 6f 6e 65 6e 74 |r each c|omponent|
|00005300| 2e 0a 58 20 2a 2f 0a 58 | 43 4f 4d 50 4c 45 58 20 |..X */.X|COMPLEX |
|00005310| 2a 0a 58 63 6d 6f 64 71 | 28 63 2c 20 71 29 0a 58 |*.Xcmodq|(c, q).X|
|00005320| 09 43 4f 4d 50 4c 45 58 | 20 2a 63 3b 0a 58 09 4e |.COMPLEX| *c;.X.N|
|00005330| 55 4d 42 45 52 20 2a 71 | 3b 0a 58 7b 0a 58 09 43 |UMBER *q|;.X{.X.C|
|00005340| 4f 4d 50 4c 45 58 20 2a | 72 3b 0a 58 0a 58 09 69 |OMPLEX *|r;.X.X.i|
|00005350| 66 20 28 71 69 73 7a 65 | 72 6f 28 71 29 29 0a 58 |f (qisze|ro(q)).X|
|00005360| 09 09 65 72 72 6f 72 28 | 22 44 69 76 69 73 69 6f |..error(|"Divisio|
|00005370| 6e 20 62 79 20 7a 65 72 | 6f 22 29 3b 0a 58 09 72 |n by zer|o");.X.r|
|00005380| 20 3d 20 63 6f 6d 61 6c | 6c 6f 63 28 29 3b 0a 58 | = comal|loc();.X|
|00005390| 09 72 2d 3e 72 65 61 6c | 20 3d 20 71 6d 6f 64 28 |.r->real| = qmod(|
|000053a0| 63 2d 3e 72 65 61 6c 2c | 20 71 29 3b 0a 58 09 72 |c->real,| q);.X.r|
|000053b0| 2d 3e 69 6d 61 67 20 3d | 20 71 6d 6f 64 28 63 2d |->imag =| qmod(c-|
|000053c0| 3e 69 6d 61 67 2c 20 71 | 29 3b 0a 58 09 72 65 74 |>imag, q|);.X.ret|
|000053d0| 75 72 6e 20 72 3b 0a 58 | 7d 0a 58 0a 58 0a 58 23 |urn r;.X|}.X.X.X#|
|000053e0| 69 66 20 30 0a 58 2f 2a | 0a 58 20 2a 20 43 6f 6e |if 0.X/*|.X * Con|
|000053f0| 73 74 72 75 63 74 20 61 | 20 63 6f 6d 70 6c 65 78 |struct a| complex|
|00005400| 20 6e 75 6d 62 65 72 20 | 67 69 76 65 6e 20 74 68 | number |given th|
|00005410| 65 20 72 65 61 6c 20 61 | 6e 64 20 69 6d 61 67 69 |e real a|nd imagi|
|00005420| 6e 61 72 79 20 63 6f 6d | 70 6f 6e 65 6e 74 73 2e |nary com|ponents.|
|00005430| 0a 58 20 2a 2f 0a 58 43 | 4f 4d 50 4c 45 58 20 2a |.X */.XC|OMPLEX *|
|00005440| 0a 58 71 71 74 6f 63 28 | 71 31 2c 20 71 32 29 0a |.Xqqtoc(|q1, q2).|
|00005450| 58 09 4e 55 4d 42 45 52 | 20 2a 71 31 2c 20 2a 71 |X.NUMBER| *q1, *q|
|00005460| 32 3b 0a 58 7b 0a 58 09 | 43 4f 4d 50 4c 45 58 20 |2;.X{.X.|COMPLEX |
|00005470| 2a 72 3b 0a 58 0a 58 09 | 69 66 20 28 71 69 73 7a |*r;.X.X.|if (qisz|
|00005480| 65 72 6f 28 71 31 29 20 | 26 26 20 71 69 73 7a 65 |ero(q1) |&& qisze|
|00005490| 72 6f 28 71 32 29 29 0a | 58 09 09 72 65 74 75 72 |ro(q2)).|X..retur|
|000054a0| 6e 20 63 6c 69 6e 6b 28 | 26 5f 63 7a 65 72 6f 5f |n clink(|&_czero_|
|000054b0| 29 3b 0a 58 09 72 20 3d | 20 63 6f 6d 61 6c 6c 6f |);.X.r =| comallo|
|000054c0| 63 28 29 3b 0a 58 09 69 | 66 20 28 21 71 69 73 7a |c();.X.i|f (!qisz|
|000054d0| 65 72 6f 28 71 31 29 29 | 0a 58 09 09 72 2d 3e 72 |ero(q1))|.X..r->r|
|000054e0| 65 61 6c 20 3d 20 71 6c | 69 6e 6b 28 71 31 29 3b |eal = ql|ink(q1);|
|000054f0| 0a 58 09 69 66 20 28 21 | 71 69 73 7a 65 72 6f 28 |.X.if (!|qiszero(|
|00005500| 71 32 29 29 0a 58 09 09 | 72 2d 3e 69 6d 61 67 20 |q2)).X..|r->imag |
|00005510| 3d 20 71 6c 69 6e 6b 28 | 71 32 29 3b 0a 58 09 72 |= qlink(|q2);.X.r|
|00005520| 65 74 75 72 6e 20 72 3b | 0a 58 7d 0a 58 23 65 6e |eturn r;|.X}.X#en|
|00005530| 64 69 66 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |dif.X.X.|X/*.X * |
|00005540| 43 6f 6d 70 61 72 65 20 | 74 77 6f 20 63 6f 6d 70 |Compare |two comp|
|00005550| 6c 65 78 20 6e 75 6d 62 | 65 72 73 20 66 6f 72 20 |lex numb|ers for |
|00005560| 65 71 75 61 6c 69 74 79 | 2c 20 72 65 74 75 72 6e |equality|, return|
|00005570| 69 6e 67 20 46 41 4c 53 | 45 20 69 66 20 74 68 65 |ing FALS|E if the|
|00005580| 79 20 61 72 65 20 65 71 | 75 61 6c 2c 0a 58 20 2a |y are eq|ual,.X *|
|00005590| 20 61 6e 64 20 54 52 55 | 45 20 69 66 20 74 68 65 | and TRU|E if the|
|000055a0| 79 20 64 69 66 66 65 72 | 2e 0a 58 20 2a 2f 0a 58 |y differ|..X */.X|
|000055b0| 42 4f 4f 4c 0a 58 63 63 | 6d 70 28 63 31 2c 20 63 |BOOL.Xcc|mp(c1, c|
|000055c0| 32 29 0a 58 09 43 4f 4d | 50 4c 45 58 20 2a 63 31 |2).X.COM|PLEX *c1|
|000055d0| 2c 20 2a 63 32 3b 0a 58 | 7b 0a 58 09 42 4f 4f 4c |, *c2;.X|{.X.BOOL|
|000055e0| 20 69 3b 0a 58 0a 58 09 | 69 20 3d 20 71 63 6d 70 | i;.X.X.|i = qcmp|
|000055f0| 28 63 31 2d 3e 72 65 61 | 6c 2c 20 63 32 2d 3e 72 |(c1->rea|l, c2->r|
|00005600| 65 61 6c 29 3b 0a 58 09 | 69 66 20 28 21 69 29 0a |eal);.X.|if (!i).|
|00005610| 58 09 09 69 20 3d 20 71 | 63 6d 70 28 63 31 2d 3e |X..i = q|cmp(c1->|
|00005620| 69 6d 61 67 2c 20 63 32 | 2d 3e 69 6d 61 67 29 3b |imag, c2|->imag);|
|00005630| 0a 58 09 72 65 74 75 72 | 6e 20 69 3b 0a 58 7d 0a |.X.retur|n i;.X}.|
|00005640| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 41 6c 6c 6f |X.X.X/*.|X * Allo|
|00005650| 63 61 74 65 20 61 20 6e | 65 77 20 63 6f 6d 70 6c |cate a n|ew compl|
|00005660| 65 78 20 6e 75 6d 62 65 | 72 2e 0a 58 20 2a 2f 0a |ex numbe|r..X */.|
|00005670| 58 43 4f 4d 50 4c 45 58 | 20 2a 0a 58 63 6f 6d 61 |XCOMPLEX| *.Xcoma|
|00005680| 6c 6c 6f 63 28 29 0a 58 | 7b 0a 58 09 43 4f 4d 50 |lloc().X|{.X.COMP|
|00005690| 4c 45 58 20 2a 72 3b 0a | 58 0a 58 09 72 20 3d 20 |LEX *r;.|X.X.r = |
|000056a0| 28 43 4f 4d 50 4c 45 58 | 20 2a 29 20 61 6c 6c 6f |(COMPLEX| *) allo|
|000056b0| 63 69 74 65 6d 28 26 66 | 72 65 65 6c 69 73 74 29 |citem(&f|reelist)|
|000056c0| 3b 0a 58 09 69 66 20 28 | 72 20 3d 3d 20 4e 55 4c |;.X.if (|r == NUL|
|000056d0| 4c 29 0a 58 09 09 65 72 | 72 6f 72 28 22 43 61 6e |L).X..er|ror("Can|
|000056e0| 6e 6f 74 20 61 6c 6c 6f | 63 61 74 65 20 63 6f 6d |not allo|cate com|
|000056f0| 70 6c 65 78 20 6e 75 6d | 62 65 72 22 29 3b 0a 58 |plex num|ber");.X|
|00005700| 09 72 2d 3e 6c 69 6e 6b | 73 20 3d 20 31 3b 0a 58 |.r->link|s = 1;.X|
|00005710| 09 72 2d 3e 72 65 61 6c | 20 3d 20 71 6c 69 6e 6b |.r->real| = qlink|
|00005720| 28 26 5f 71 7a 65 72 6f | 5f 29 3b 0a 58 09 72 2d |(&_qzero|_);.X.r-|
|00005730| 3e 69 6d 61 67 20 3d 20 | 71 6c 69 6e 6b 28 26 5f |>imag = |qlink(&_|
|00005740| 71 7a 65 72 6f 5f 29 3b | 0a 58 09 72 65 74 75 72 |qzero_);|.X.retur|
|00005750| 6e 20 72 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 0a |n r;.X}.|X.X.X/*.|
|00005760| 58 20 2a 20 46 72 65 65 | 20 61 20 63 6f 6d 70 6c |X * Free| a compl|
|00005770| 65 78 20 6e 75 6d 62 65 | 72 2e 0a 58 20 2a 2f 0a |ex numbe|r..X */.|
|00005780| 58 76 6f 69 64 0a 58 63 | 6f 6d 66 72 65 65 28 63 |Xvoid.Xc|omfree(c|
|00005790| 29 0a 58 09 43 4f 4d 50 | 4c 45 58 20 2a 63 3b 0a |).X.COMP|LEX *c;.|
|000057a0| 58 7b 0a 58 09 69 66 20 | 28 2d 2d 28 63 2d 3e 6c |X{.X.if |(--(c->l|
|000057b0| 69 6e 6b 73 29 20 3e 20 | 30 29 0a 58 09 09 72 65 |inks) > |0).X..re|
|000057c0| 74 75 72 6e 3b 0a 58 09 | 71 66 72 65 65 28 63 2d |turn;.X.|qfree(c-|
|000057d0| 3e 72 65 61 6c 29 3b 0a | 58 09 71 66 72 65 65 28 |>real);.|X.qfree(|
|000057e0| 63 2d 3e 69 6d 61 67 29 | 3b 0a 58 09 66 72 65 65 |c->imag)|;.X.free|
|000057f0| 69 74 65 6d 28 26 66 72 | 65 65 6c 69 73 74 2c 20 |item(&fr|eelist, |
|00005800| 28 46 52 45 45 49 54 45 | 4d 20 2a 29 20 63 29 3b |(FREEITE|M *) c);|
|00005810| 0a 58 7d 0a 58 0a 58 2f | 2a 20 45 4e 44 20 43 4f |.X}.X.X/|* END CO|
|00005820| 44 45 20 2a 2f 0a 45 4e | 44 5f 4f 46 5f 46 49 4c |DE */.EN|D_OF_FIL|
|00005830| 45 0a 69 66 20 74 65 73 | 74 20 39 36 35 34 20 2d |E.if tes|t 9654 -|
|00005840| 6e 65 20 60 77 63 20 2d | 63 20 3c 27 63 6f 6d 6d |ne `wc -|c <'comm|
|00005850| 61 74 68 2e 63 27 60 3b | 20 74 68 65 6e 0a 20 20 |ath.c'`;| then.  |
|00005860| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 27 |  echo s|har: \"'|
|00005870| 63 6f 6d 6d 61 74 68 2e | 63 27 5c 22 20 75 6e 70 |commath.|c'\" unp|
|00005880| 61 63 6b 65 64 20 77 69 | 74 68 20 77 72 6f 6e 67 |acked wi|th wrong|
|00005890| 20 73 69 7a 65 21 0a 66 | 69 0a 23 20 65 6e 64 20 | size!.f|i.# end |
|000058a0| 6f 66 20 27 63 6f 6d 6d | 61 74 68 2e 63 27 0a 66 |of 'comm|ath.c'.f|
|000058b0| 69 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 66 69 |i.if tes|t -f 'fi|
|000058c0| 6c 65 2e 63 27 20 2d 61 | 20 22 24 7b 31 7d 22 20 |le.c' -a| "${1}" |
|000058d0| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|000058e0| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c |  echo s|har: Wil|
|000058f0| 6c 20 6e 6f 74 20 63 6c | 6f 62 62 65 72 20 65 78 |l not cl|obber ex|
|00005900| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 27 66 |isting f|ile \"'f|
|00005910| 69 6c 65 2e 63 27 5c 22 | 0a 65 6c 73 65 0a 65 63 |ile.c'\"|.else.ec|
|00005920| 68 6f 20 73 68 61 72 3a | 20 45 78 74 72 61 63 74 |ho shar:| Extract|
|00005930| 69 6e 67 20 5c 22 27 66 | 69 6c 65 2e 63 27 5c 22 |ing \"'f|ile.c'\"|
|00005940| 20 5c 28 31 30 34 30 37 | 20 63 68 61 72 61 63 74 | \(10407| charact|
|00005950| 65 72 73 5c 29 0a 73 65 | 64 20 22 73 2f 5e 58 2f |ers\).se|d "s/^X/|
|00005960| 2f 22 20 3e 27 66 69 6c | 65 2e 63 27 20 3c 3c 27 |/" >'fil|e.c' <<'|
|00005970| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 27 0a 58 2f 2a |END_OF_F|ILE'.X/*|
|00005980| 0a 58 20 2a 20 43 6f 70 | 79 72 69 67 68 74 20 28 |.X * Cop|yright (|
|00005990| 63 29 20 31 39 39 32 20 | 44 61 76 69 64 20 49 2e |c) 1992 |David I.|
|000059a0| 20 42 65 6c 6c 0a 58 20 | 2a 20 50 65 72 6d 69 73 | Bell.X |* Permis|
|000059b0| 73 69 6f 6e 20 69 73 20 | 67 72 61 6e 74 65 64 20 |sion is |granted |
|000059c0| 74 6f 20 75 73 65 2c 20 | 64 69 73 74 72 69 62 75 |to use, |distribu|
|000059d0| 74 65 2c 20 6f 72 20 6d | 6f 64 69 66 79 20 74 68 |te, or m|odify th|
|000059e0| 69 73 20 73 6f 75 72 63 | 65 2c 0a 58 20 2a 20 70 |is sourc|e,.X * p|
|000059f0| 72 6f 76 69 64 65 64 20 | 74 68 61 74 20 74 68 69 |rovided |that thi|
|00005a00| 73 20 63 6f 70 79 72 69 | 67 68 74 20 6e 6f 74 69 |s copyri|ght noti|
|00005a10| 63 65 20 72 65 6d 61 69 | 6e 73 20 69 6e 74 61 63 |ce remai|ns intac|
|00005a20| 74 2e 0a 58 20 2a 0a 58 | 20 2a 20 46 69 6c 65 20 |t..X *.X| * File |
|00005a30| 49 2f 4f 20 72 6f 75 74 | 69 6e 65 73 20 63 61 6c |I/O rout|ines cal|
|00005a40| 6c 61 62 6c 65 20 62 79 | 20 75 73 65 72 73 2e 0a |lable by| users..|
|00005a50| 58 20 2a 2f 0a 58 0a 58 | 23 69 6e 63 6c 75 64 65 |X */.X.X|#include|
|00005a60| 20 22 73 74 64 61 72 67 | 2e 68 22 0a 58 23 69 6e | "stdarg|.h".X#in|
|00005a70| 63 6c 75 64 65 20 22 63 | 61 6c 63 2e 68 22 0a 58 |clude "c|alc.h".X|
|00005a80| 0a 58 0a 58 23 64 65 66 | 69 6e 65 09 52 45 41 44 |.X.X#def|ine.READ|
|00005a90| 53 49 5a 45 09 31 30 32 | 34 09 2f 2a 20 62 75 66 |SIZE.102|4./* buf|
|00005aa0| 66 65 72 20 73 69 7a 65 | 20 66 6f 72 20 72 65 61 |fer size| for rea|
|00005ab0| 64 69 6e 67 20 2a 2f 0a | 58 0a 58 2f 2a 0a 58 20 |ding */.|X.X/*.X |
|00005ac0| 2a 20 44 65 66 69 6e 69 | 74 69 6f 6e 20 6f 66 20 |* Defini|tion of |
|00005ad0| 6f 70 65 6e 65 64 20 66 | 69 6c 65 73 2e 0a 58 20 |opened f|iles..X |
|00005ae0| 2a 2f 0a 58 74 79 70 65 | 64 65 66 20 73 74 72 75 |*/.Xtype|def stru|
|00005af0| 63 74 20 7b 0a 58 09 46 | 49 4c 45 49 44 20 69 64 |ct {.X.F|ILEID id|
|00005b00| 3b 09 09 2f 2a 20 69 64 | 20 74 6f 20 69 64 65 6e |;../* id| to iden|
|00005b10| 74 69 66 79 20 74 68 69 | 73 20 66 69 6c 65 20 2a |tify thi|s file *|
|00005b20| 2f 0a 58 09 46 49 4c 45 | 20 2a 66 70 3b 09 09 2f |/.X.FILE| *fp;../|
|00005b30| 2a 20 72 65 61 6c 20 66 | 69 6c 65 20 73 74 72 75 |* real f|ile stru|
|00005b40| 63 74 75 72 65 20 66 6f | 72 20 49 2f 4f 20 2a 2f |cture fo|r I/O */|
|00005b50| 0a 58 09 63 68 61 72 20 | 2a 6e 61 6d 65 3b 09 09 |.X.char |*name;..|
|00005b60| 2f 2a 20 66 69 6c 65 20 | 6e 61 6d 65 20 2a 2f 0a |/* file |name */.|
|00005b70| 58 09 42 4f 4f 4c 20 72 | 65 61 64 69 6e 67 3b 09 |X.BOOL r|eading;.|
|00005b80| 09 2f 2a 20 54 52 55 45 | 20 69 66 20 6f 70 65 6e |./* TRUE| if open|
|00005b90| 65 64 20 66 6f 72 20 72 | 65 61 64 69 6e 67 20 2a |ed for r|eading *|
|00005ba0| 2f 0a 58 09 42 4f 4f 4c | 20 77 72 69 74 69 6e 67 |/.X.BOOL| writing|
|00005bb0| 3b 09 09 2f 2a 20 54 52 | 55 45 20 69 66 20 6f 70 |;../* TR|UE if op|
|00005bc0| 65 6e 65 64 20 66 6f 72 | 20 77 72 69 74 69 6e 67 |ened for| writing|
|00005bd0| 20 2a 2f 0a 58 09 63 68 | 61 72 20 2a 6d 6f 64 65 | */.X.ch|ar *mode|
|00005be0| 3b 09 09 2f 2a 20 6f 70 | 65 6e 20 6d 6f 64 65 20 |;../* op|en mode |
|00005bf0| 2a 2f 0a 58 7d 20 46 49 | 4c 45 49 4f 3b 0a 58 0a |*/.X} FI|LEIO;.X.|
|00005c00| 58 0a 58 2f 2a 0a 58 20 | 2a 20 54 61 62 6c 65 20 |X.X/*.X |* Table |
|00005c10| 6f 66 20 6f 70 65 6e 65 | 64 20 66 69 6c 65 73 2e |of opene|d files.|
|00005c20| 0a 58 20 2a 20 54 68 65 | 20 66 69 72 73 74 20 74 |.X * The| first t|
|00005c30| 68 72 65 65 20 65 6e 74 | 72 69 65 73 20 61 6c 77 |hree ent|ries alw|
|00005c40| 61 79 73 20 63 6f 72 72 | 65 73 70 6f 6e 64 20 74 |ays corr|espond t|
|00005c50| 6f 20 73 74 64 69 6e 2c | 20 73 74 64 6f 75 74 2c |o stdin,| stdout,|
|00005c60| 20 61 6e 64 20 73 74 64 | 65 72 72 2c 0a 58 20 2a | and std|err,.X *|
|00005c70| 20 61 6e 64 20 63 61 6e | 6e 6f 74 20 62 65 20 63 | and can|not be c|
|00005c80| 6c 6f 73 65 64 2e 20 20 | 54 68 65 69 72 20 66 69 |losed.  |Their fi|
|00005c90| 6c 65 20 69 64 73 20 61 | 72 65 20 61 6c 77 61 79 |le ids a|re alway|
|00005ca0| 73 20 30 2c 20 31 2c 20 | 61 6e 64 20 32 2e 0a 58 |s 0, 1, |and 2..X|
|00005cb0| 20 2a 2f 0a 58 73 74 61 | 74 69 63 20 46 49 4c 45 | */.Xsta|tic FILE|
|00005cc0| 49 4f 20 66 69 6c 65 73 | 5b 4d 41 58 46 49 4c 45 |IO files|[MAXFILE|
|00005cd0| 53 5d 20 3d 20 7b 0a 58 | 09 46 49 4c 45 49 44 5f |S] = {.X|.FILEID_|
|00005ce0| 53 54 44 49 4e 2c 20 20 | 73 74 64 69 6e 2c 20 20 |STDIN,  |stdin,  |
|00005cf0| 22 28 73 74 64 69 6e 29 | 22 2c 20 20 54 52 55 45 |"(stdin)|",  TRUE|
|00005d00| 2c 20 46 41 4c 53 45 2c | 20 22 72 65 61 64 69 6e |, FALSE,| "readin|
|00005d10| 67 22 2c 0a 58 09 46 49 | 4c 45 49 44 5f 53 54 44 |g",.X.FI|LEID_STD|
|00005d20| 4f 55 54 2c 20 73 74 64 | 6f 75 74 2c 20 22 28 73 |OUT, std|out, "(s|
|00005d30| 74 64 6f 75 74 29 22 2c | 20 46 41 4c 53 45 2c 20 |tdout)",| FALSE, |
|00005d40| 54 52 55 45 2c 20 22 77 | 72 69 74 69 6e 67 22 2c |TRUE, "w|riting",|
|00005d50| 0a 58 09 46 49 4c 45 49 | 44 5f 53 54 44 45 52 52 |.X.FILEI|D_STDERR|
|00005d60| 2c 20 73 74 64 65 72 72 | 2c 20 22 28 73 74 64 65 |, stderr|, "(stde|
|00005d70| 72 72 29 22 2c 20 46 41 | 4c 53 45 2c 20 54 52 55 |rr)", FA|LSE, TRU|
|00005d80| 45 2c 20 22 77 72 69 74 | 69 6e 67 22 0a 58 7d 3b |E, "writ|ing".X};|
|00005d90| 0a 58 0a 58 73 74 61 74 | 69 63 20 46 49 4c 45 49 |.X.Xstat|ic FILEI|
|00005da0| 44 20 6c 61 73 74 69 64 | 20 3d 20 46 49 4c 45 49 |D lastid| = FILEI|
|00005db0| 44 5f 53 54 44 45 52 52 | 3b 09 09 2f 2a 20 6c 61 |D_STDERR|;../* la|
|00005dc0| 73 74 20 61 6c 6c 6f 63 | 61 74 65 64 20 66 69 6c |st alloc|ated fil|
|00005dd0| 65 20 69 64 20 2a 2f 0a | 58 0a 58 0a 58 0a 58 2f |e id */.|X.X.X.X/|
|00005de0| 2a 0a 58 20 2a 20 4f 70 | 65 6e 20 74 68 65 20 73 |*.X * Op|en the s|
|00005df0| 70 65 63 69 66 69 65 64 | 20 66 69 6c 65 20 6e 61 |pecified| file na|
|00005e00| 6d 65 20 66 6f 72 20 72 | 65 61 64 69 6e 67 20 6f |me for r|eading o|
|00005e10| 72 20 77 72 69 74 69 6e | 67 20 61 73 20 64 65 74 |r writin|g as det|
|00005e20| 65 72 6d 69 6e 65 64 20 | 62 79 0a 58 20 2a 20 74 |ermined |by.X * t|
|00005e30| 68 65 20 73 70 65 63 69 | 66 69 65 64 20 6d 6f 64 |he speci|fied mod|
|00005e40| 65 20 28 22 72 22 2c 20 | 22 77 22 2c 20 6f 72 20 |e ("r", |"w", or |
|00005e50| 22 61 22 29 2e 20 20 52 | 65 74 75 72 6e 73 20 61 |"a").  R|eturns a|
|00005e60| 20 66 69 6c 65 20 69 64 | 20 77 68 69 63 68 20 63 | file id| which c|
|00005e70| 61 6e 20 62 65 0a 58 20 | 2a 20 75 73 65 64 20 74 |an be.X |* used t|
|00005e80| 6f 20 64 6f 20 49 2f 4f | 20 74 6f 20 74 68 65 20 |o do I/O| to the |
|00005e90| 66 69 6c 65 2c 20 6f 72 | 20 65 6c 73 65 20 46 49 |file, or| else FI|
|00005ea0| 4c 45 49 44 5f 4e 4f 4e | 45 20 69 66 20 74 68 65 |LEID_NON|E if the|
|00005eb0| 20 6f 70 65 6e 20 66 61 | 69 6c 65 64 2e 0a 58 20 | open fa|iled..X |
|00005ec0| 2a 20 41 62 6f 72 74 73 | 20 77 69 74 68 20 61 6e |* Aborts| with an|
|00005ed0| 20 65 72 72 6f 72 20 69 | 66 20 74 6f 6f 20 6d 61 | error i|f too ma|
|00005ee0| 6e 79 20 66 69 6c 65 73 | 20 61 72 65 20 6f 70 65 |ny files| are ope|
|00005ef0| 6e 65 64 20 6f 72 20 74 | 68 65 20 6d 6f 64 65 20 |ned or t|he mode |
|00005f00| 69 73 20 69 6c 6c 65 67 | 61 6c 2e 0a 58 20 2a 2f |is illeg|al..X */|
|00005f10| 0a 58 46 49 4c 45 49 44 | 0a 58 6f 70 65 6e 69 64 |.XFILEID|.Xopenid|
|00005f20| 28 6e 61 6d 65 2c 20 6d | 6f 64 65 29 0a 58 09 63 |(name, m|ode).X.c|
|00005f30| 68 61 72 20 2a 6e 61 6d | 65 3b 09 09 2f 2a 20 66 |har *nam|e;../* f|
|00005f40| 69 6c 65 20 6e 61 6d 65 | 20 2a 2f 0a 58 09 63 68 |ile name| */.X.ch|
|00005f50| 61 72 20 2a 6d 6f 64 65 | 3b 09 09 2f 2a 20 6f 70 |ar *mode|;../* op|
|00005f60| 65 6e 20 6d 6f 64 65 20 | 2a 2f 0a 58 7b 0a 58 09 |en mode |*/.X{.X.|
|00005f70| 46 49 4c 45 49 4f 20 2a | 66 69 6f 70 3b 09 09 2f |FILEIO *|fiop;../|
|00005f80| 2a 20 66 69 6c 65 20 73 | 74 72 75 63 74 75 72 65 |* file s|tructure|
|00005f90| 20 2a 2f 0a 58 09 46 49 | 4c 45 49 44 20 69 64 3b | */.X.FI|LEID id;|
|00005fa0| 09 09 2f 2a 20 6e 65 77 | 20 66 69 6c 65 20 69 64 |../* new| file id|
|00005fb0| 20 2a 2f 0a 58 09 69 6e | 74 20 63 6f 75 6e 74 3b | */.X.in|t count;|
|00005fc0| 0a 58 0a 58 09 69 66 20 | 28 28 28 2a 6d 6f 64 65 |.X.X.if |(((*mode|
|00005fd0| 20 21 3d 20 27 72 27 29 | 20 26 26 20 28 2a 6d 6f | != 'r')| && (*mo|
|00005fe0| 64 65 20 21 3d 20 27 77 | 27 29 20 26 26 20 28 2a |de != 'w|') && (*|
|00005ff0| 6d 6f 64 65 20 21 3d 20 | 27 61 27 29 29 20 7c 7c |mode != |'a')) |||
|00006000| 20 6d 6f 64 65 5b 31 5d | 29 0a 58 09 09 65 72 72 | mode[1]|).X..err|
|00006010| 6f 72 28 22 49 6c 6c 65 | 67 61 6c 20 6d 6f 64 65 |or("Ille|gal mode|
|00006020| 20 66 6f 72 20 66 6f 70 | 65 6e 22 29 3b 0a 58 0a | for fop|en");.X.|
|00006030| 58 09 63 6f 75 6e 74 20 | 3d 20 4d 41 58 46 49 4c |X.count |= MAXFIL|
|00006040| 45 53 3b 0a 58 09 64 6f | 20 7b 0a 58 09 09 69 66 |ES;.X.do| {.X..if|
|00006050| 20 28 2d 2d 63 6f 75 6e | 74 20 3c 20 30 29 0a 58 | (--coun|t < 0).X|
|00006060| 09 09 09 65 72 72 6f 72 | 28 22 54 6f 6f 20 6d 61 |...error|("Too ma|
|00006070| 6e 79 20 6f 70 65 6e 20 | 66 69 6c 65 73 22 29 3b |ny open |files");|
|00006080| 0a 58 09 09 69 64 20 3d | 20 2b 2b 6c 61 73 74 69 |.X..id =| ++lasti|
|00006090| 64 3b 0a 58 09 09 66 69 | 6f 70 20 3d 20 26 66 69 |d;.X..fi|op = &fi|
|000060a0| 6c 65 73 5b 69 64 20 25 | 20 4d 41 58 46 49 4c 45 |les[id %| MAXFILE|
|000060b0| 53 5d 3b 0a 58 0a 58 09 | 7d 20 77 68 69 6c 65 20 |S];.X.X.|} while |
|000060c0| 28 66 69 6f 70 2d 3e 72 | 65 61 64 69 6e 67 20 7c |(fiop->r|eading ||
|000060d0| 7c 20 66 69 6f 70 2d 3e | 77 72 69 74 69 6e 67 29 || fiop->|writing)|
|000060e0| 3b 0a 58 0a 58 09 66 69 | 6f 70 2d 3e 6e 61 6d 65 |;.X.X.fi|op->name|
|000060f0| 20 3d 20 28 63 68 61 72 | 20 2a 29 6d 61 6c 6c 6f | = (char| *)mallo|
|00006100| 63 28 73 74 72 6c 65 6e | 28 6e 61 6d 65 29 20 2b |c(strlen|(name) +|
|00006110| 20 31 29 3b 0a 58 09 69 | 66 20 28 66 69 6f 70 2d | 1);.X.i|f (fiop-|
|00006120| 3e 6e 61 6d 65 20 3d 3d | 20 4e 55 4c 4c 29 20 7b |>name ==| NULL) {|
|00006130| 0a 58 09 09 6c 61 73 74 | 69 64 2d 2d 3b 0a 58 09 |.X..last|id--;.X.|
|00006140| 09 65 72 72 6f 72 28 22 | 4e 6f 20 6d 65 6d 6f 72 |.error("|No memor|
|00006150| 79 20 66 6f 72 20 66 69 | 6c 65 6e 61 6d 65 22 29 |y for fi|lename")|
|00006160| 3b 0a 58 09 7d 0a 58 09 | 73 74 72 63 70 79 28 66 |;.X.}.X.|strcpy(f|
|00006170| 69 6f 70 2d 3e 6e 61 6d | 65 2c 20 6e 61 6d 65 29 |iop->nam|e, name)|
|00006180| 3b 0a 58 0a 58 09 66 69 | 6f 70 2d 3e 66 70 20 3d |;.X.X.fi|op->fp =|
|00006190| 20 66 5f 6f 70 65 6e 28 | 6e 61 6d 65 2c 20 6d 6f | f_open(|name, mo|
|000061a0| 64 65 29 3b 0a 58 09 69 | 66 20 28 66 69 6f 70 2d |de);.X.i|f (fiop-|
|000061b0| 3e 66 70 20 3d 3d 20 4e | 55 4c 4c 29 20 7b 0a 58 |>fp == N|ULL) {.X|
|000061c0| 09 09 66 72 65 65 28 66 | 69 6f 70 2d 3e 6e 61 6d |..free(f|iop->nam|
|000061d0| 65 29 3b 0a 58 09 09 66 | 69 6f 70 2d 3e 6e 61 6d |e);.X..f|iop->nam|
|000061e0| 65 20 3d 20 4e 55 4c 4c | 3b 0a 58 09 09 6c 61 73 |e = NULL|;.X..las|
|000061f0| 74 69 64 2d 2d 3b 0a 58 | 09 09 72 65 74 75 72 6e |tid--;.X|..return|
|00006200| 20 46 49 4c 45 49 44 5f | 4e 4f 4e 45 3b 0a 58 09 | FILEID_|NONE;.X.|
|00006210| 7d 0a 58 0a 58 09 73 77 | 69 74 63 68 20 28 2a 6d |}.X.X.sw|itch (*m|
|00006220| 6f 64 65 29 20 7b 0a 58 | 09 09 63 61 73 65 20 27 |ode) {.X|..case '|
|00006230| 72 27 3a 0a 58 09 09 09 | 66 69 6f 70 2d 3e 6d 6f |r':.X...|fiop->mo|
|00006240| 64 65 20 3d 20 22 72 65 | 61 64 69 6e 67 22 3b 0a |de = "re|ading";.|
|00006250| 58 09 09 09 66 69 6f 70 | 2d 3e 72 65 61 64 69 6e |X...fiop|->readin|
|00006260| 67 20 3d 20 54 52 55 45 | 3b 0a 58 09 09 09 62 72 |g = TRUE|;.X...br|
|00006270| 65 61 6b 3b 0a 58 09 09 | 63 61 73 65 20 27 77 27 |eak;.X..|case 'w'|
|00006280| 3a 0a 58 09 09 09 66 69 | 6f 70 2d 3e 6d 6f 64 65 |:.X...fi|op->mode|
|00006290| 20 3d 20 22 77 72 69 74 | 69 6e 67 22 3b 0a 58 09 | = "writ|ing";.X.|
|000062a0| 09 09 66 69 6f 70 2d 3e | 77 72 69 74 69 6e 67 20 |..fiop->|writing |
|000062b0| 3d 20 54 52 55 45 3b 0a | 58 09 09 09 62 72 65 61 |= TRUE;.|X...brea|
|000062c0| 6b 3b 0a 58 09 09 63 61 | 73 65 20 27 61 27 3a 0a |k;.X..ca|se 'a':.|
|000062d0| 58 09 09 09 66 69 6f 70 | 2d 3e 6d 6f 64 65 20 3d |X...fiop|->mode =|
|000062e0| 20 22 61 70 70 65 6e 64 | 69 6e 67 22 3b 0a 58 09 | "append|ing";.X.|
|000062f0| 09 09 66 69 6f 70 2d 3e | 77 72 69 74 69 6e 67 20 |..fiop->|writing |
|00006300| 3d 20 54 52 55 45 3b 0a | 58 09 09 09 62 72 65 61 |= TRUE;.|X...brea|
|00006310| 6b 3b 0a 58 09 7d 0a 58 | 0a 58 09 66 69 6f 70 2d |k;.X.}.X|.X.fiop-|
|00006320| 3e 69 64 20 3d 20 69 64 | 3b 0a 58 0a 58 09 72 65 |>id = id|;.X.X.re|
|00006330| 74 75 72 6e 20 69 64 3b | 0a 58 7d 0a 58 0a 58 0a |turn id;|.X}.X.X.|
|00006340| 58 2f 2a 0a 58 20 2a 20 | 46 69 6e 64 20 74 68 65 |X/*.X * |Find the|
|00006350| 20 66 69 6c 65 20 49 2f | 4f 20 73 74 72 75 63 74 | file I/|O struct|
|00006360| 75 72 65 20 66 6f 72 20 | 74 68 65 20 73 70 65 63 |ure for |the spec|
|00006370| 69 66 69 65 64 20 66 69 | 6c 65 20 69 64 2c 20 61 |ified fi|le id, a|
|00006380| 6e 64 20 76 65 72 69 66 | 79 20 74 68 61 74 0a 58 |nd verif|y that.X|
|00006390| 20 2a 20 69 74 20 69 73 | 20 6f 70 65 6e 65 64 20 | * it is| opened |
|000063a0| 69 6e 20 74 68 65 20 72 | 65 71 75 69 72 65 64 20 |in the r|equired |
|000063b0| 6d 61 6e 6e 65 72 20 28 | 27 72 27 20 66 6f 72 20 |manner (|'r' for |
|000063c0| 72 65 61 64 69 6e 67 20 | 6f 72 20 27 77 27 20 66 |reading |or 'w' f|
|000063d0| 6f 72 20 77 72 69 74 69 | 6e 67 29 2e 0a 58 20 2a |or writi|ng)..X *|
|000063e0| 20 49 66 20 6d 6f 64 65 | 20 69 73 20 30 2c 20 74 | If mode| is 0, t|
|000063f0| 68 65 6e 20 6e 6f 20 6f | 70 65 6e 20 63 68 65 63 |hen no o|pen chec|
+--------+-------------------------+-------------------------+--------+--------+