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

< 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 35 |Subject:| v26i035|
|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 | 39 2f 32 31 0a 53 65 6e |r, Part0|9/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 35 |ume 26, |Issue 35|
|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 39 0a 0a 23 21 20 2f |alc/part|09..#! /|
|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 39 20 28 6f 66 20 32 | 31 29 2e 22 0a 23 20 43 | 9 (of 2|1).".# C|
|000002e0| 6f 6e 74 65 6e 74 73 3a | 20 20 71 6d 61 74 68 2e |ontents:|  qmath.|
|000002f0| 63 20 71 74 72 61 6e 73 | 2e 63 0a 23 20 57 72 61 |c qtrans|.c.# Wra|
|00000300| 70 70 65 64 20 62 79 20 | 64 62 65 6c 6c 40 65 6c |pped by |dbell@el|
|00000310| 6d 20 6f 6e 20 54 75 65 | 20 46 65 62 20 32 35 20 |m on Tue| Feb 25 |
|00000320| 31 35 3a 32 31 3a 30 34 | 20 31 39 39 32 0a 50 41 |15:21:04| 1992.PA|
|00000330| 54 48 3d 2f 62 69 6e 3a | 2f 75 73 72 2f 62 69 6e |TH=/bin:|/usr/bin|
|00000340| 3a 2f 75 73 72 2f 75 63 | 62 20 3b 20 65 78 70 6f |:/usr/uc|b ; expo|
|00000350| 72 74 20 50 41 54 48 0a | 69 66 20 74 65 73 74 20 |rt PATH.|if test |
|00000360| 2d 66 20 27 71 6d 61 74 | 68 2e 63 27 20 2d 61 20 |-f 'qmat|h.c' -a |
|00000370| 22 24 7b 31 7d 22 20 21 | 3d 20 22 2d 63 22 20 3b |"${1}" !|= "-c" ;|
|00000380| 20 74 68 65 6e 20 0a 20 | 20 65 63 68 6f 20 73 68 | then . | echo sh|
|00000390| 61 72 3a 20 57 69 6c 6c | 20 6e 6f 74 20 63 6c 6f |ar: Will| not clo|
|000003a0| 62 62 65 72 20 65 78 69 | 73 74 69 6e 67 20 66 69 |bber exi|sting fi|
|000003b0| 6c 65 20 5c 22 27 71 6d | 61 74 68 2e 63 27 5c 22 |le \"'qm|ath.c'\"|
|000003c0| 0a 65 6c 73 65 0a 65 63 | 68 6f 20 73 68 61 72 3a |.else.ec|ho shar:|
|000003d0| 20 45 78 74 72 61 63 74 | 69 6e 67 20 5c 22 27 71 | Extract|ing \"'q|
|000003e0| 6d 61 74 68 2e 63 27 5c | 22 20 5c 28 32 30 31 34 |math.c'\|" \(2014|
|000003f0| 34 20 63 68 61 72 61 63 | 74 65 72 73 5c 29 0a 73 |4 charac|ters\).s|
|00000400| 65 64 20 22 73 2f 5e 58 | 2f 2f 22 20 3e 27 71 6d |ed "s/^X|//" >'qm|
|00000410| 61 74 68 2e 63 27 20 3c | 3c 27 45 4e 44 5f 4f 46 |ath.c' <|<'END_OF|
|00000420| 5f 46 49 4c 45 27 0a 58 | 2f 2a 0a 58 20 2a 20 43 |_FILE'.X|/*.X * C|
|00000430| 6f 70 79 72 69 67 68 74 | 20 28 63 29 20 31 39 39 |opyright| (c) 199|
|00000440| 32 20 44 61 76 69 64 20 | 49 2e 20 42 65 6c 6c 0a |2 David |I. Bell.|
|00000450| 58 20 2a 20 50 65 72 6d | 69 73 73 69 6f 6e 20 69 |X * Perm|ission i|
|00000460| 73 20 67 72 61 6e 74 65 | 64 20 74 6f 20 75 73 65 |s grante|d to use|
|00000470| 2c 20 64 69 73 74 72 69 | 62 75 74 65 2c 20 6f 72 |, distri|bute, or|
|00000480| 20 6d 6f 64 69 66 79 20 | 74 68 69 73 20 73 6f 75 | modify |this sou|
|00000490| 72 63 65 2c 0a 58 20 2a | 20 70 72 6f 76 69 64 65 |rce,.X *| provide|
|000004a0| 64 20 74 68 61 74 20 74 | 68 69 73 20 63 6f 70 79 |d that t|his copy|
|000004b0| 72 69 67 68 74 20 6e 6f | 74 69 63 65 20 72 65 6d |right no|tice rem|
|000004c0| 61 69 6e 73 20 69 6e 74 | 61 63 74 2e 0a 58 20 2a |ains int|act..X *|
|000004d0| 0a 58 20 2a 20 45 78 74 | 65 6e 64 65 64 20 70 72 |.X * Ext|ended pr|
|000004e0| 65 63 69 73 69 6f 6e 20 | 72 61 74 69 6f 6e 61 6c |ecision |rational|
|000004f0| 20 61 72 69 74 68 6d 65 | 74 69 63 20 70 72 69 6d | arithme|tic prim|
|00000500| 69 74 69 76 65 20 72 6f | 75 74 69 6e 65 73 0a 58 |itive ro|utines.X|
|00000510| 20 2a 2f 0a 58 0a 58 23 | 69 6e 63 6c 75 64 65 20 | */.X.X#|include |
|00000520| 3c 73 74 64 69 6f 2e 68 | 3e 0a 58 23 69 6e 63 6c |<stdio.h|>.X#incl|
|00000530| 75 64 65 20 22 6d 61 74 | 68 2e 68 22 0a 58 0a 58 |ude "mat|h.h".X.X|
|00000540| 0a 58 4e 55 4d 42 45 52 | 20 5f 71 7a 65 72 6f 5f |.XNUMBER| _qzero_|
|00000550| 20 3d 09 7b 20 7b 20 5f | 7a 65 72 6f 76 61 6c 5f | =.{ { _|zeroval_|
|00000560| 2c 20 31 2c 20 30 20 7d | 2c 20 7b 20 5f 6f 6e 65 |, 1, 0 }|, { _one|
|00000570| 76 61 6c 5f 2c 20 31 2c | 20 30 20 7d 2c 20 31 20 |val_, 1,| 0 }, 1 |
|00000580| 7d 3b 0a 58 4e 55 4d 42 | 45 52 20 5f 71 6f 6e 65 |};.XNUMB|ER _qone|
|00000590| 5f 20 3d 09 09 7b 20 7b | 20 5f 6f 6e 65 76 61 6c |_ =..{ {| _oneval|
|000005a0| 5f 2c 20 31 2c 20 30 20 | 7d 2c 20 7b 20 5f 6f 6e |_, 1, 0 |}, { _on|
|000005b0| 65 76 61 6c 5f 2c 20 31 | 2c 20 30 20 7d 2c 20 31 |eval_, 1|, 0 }, 1|
|000005c0| 20 7d 3b 0a 58 73 74 61 | 74 69 63 20 4e 55 4d 42 | };.Xsta|tic NUMB|
|000005d0| 45 52 20 5f 71 74 77 6f | 5f 20 3d 09 7b 20 7b 20 |ER _qtwo|_ =.{ { |
|000005e0| 5f 74 77 6f 76 61 6c 5f | 2c 20 31 2c 20 30 20 7d |_twoval_|, 1, 0 }|
|000005f0| 2c 20 7b 20 5f 6f 6e 65 | 76 61 6c 5f 2c 20 31 2c |, { _one|val_, 1,|
|00000600| 20 30 20 7d 2c 20 31 20 | 7d 3b 0a 58 73 74 61 74 | 0 }, 1 |};.Xstat|
|00000610| 69 63 20 4e 55 4d 42 45 | 52 20 5f 71 74 65 6e 5f |ic NUMBE|R _qten_|
|00000620| 20 3d 09 7b 20 7b 20 5f | 74 65 6e 76 61 6c 5f 2c | =.{ { _|tenval_,|
|00000630| 20 31 2c 20 30 20 7d 2c | 20 7b 20 5f 6f 6e 65 76 | 1, 0 },| { _onev|
|00000640| 61 6c 5f 2c 20 31 2c 20 | 30 20 7d 2c 20 31 20 7d |al_, 1, |0 }, 1 }|
|00000650| 3b 0a 58 4e 55 4d 42 45 | 52 20 5f 71 6e 65 67 6f |;.XNUMBE|R _qnego|
|00000660| 6e 65 5f 20 3d 09 7b 20 | 7b 20 5f 6f 6e 65 76 61 |ne_ =.{ |{ _oneva|
|00000670| 6c 5f 2c 20 31 2c 20 31 | 20 7d 2c 20 7b 20 5f 6f |l_, 1, 1| }, { _o|
|00000680| 6e 65 76 61 6c 5f 2c 20 | 31 2c 20 30 20 7d 2c 20 |neval_, |1, 0 }, |
|00000690| 31 20 7d 3b 0a 58 4e 55 | 4d 42 45 52 20 5f 71 6f |1 };.XNU|MBER _qo|
|000006a0| 6e 65 68 61 6c 66 5f 20 | 3d 09 7b 20 7b 20 5f 6f |nehalf_ |=.{ { _o|
|000006b0| 6e 65 76 61 6c 5f 2c 20 | 31 2c 20 30 20 7d 2c 20 |neval_, |1, 0 }, |
|000006c0| 7b 20 5f 74 77 6f 76 61 | 6c 5f 2c 20 31 2c 20 30 |{ _twova|l_, 1, 0|
|000006d0| 20 7d 2c 20 31 20 7d 3b | 0a 58 0a 58 0a 58 23 69 | }, 1 };|.X.X.X#i|
|000006e0| 66 20 30 0a 58 73 74 61 | 74 69 63 20 63 68 61 72 |f 0.Xsta|tic char|
|000006f0| 20 2a 61 62 6f 72 74 6d | 73 67 20 3d 20 22 43 61 | *abortm|sg = "Ca|
|00000700| 6c 63 75 6c 61 74 69 6f | 6e 20 61 62 6f 72 74 65 |lculatio|n aborte|
|00000710| 64 22 3b 0a 58 23 65 6e | 64 69 66 0a 58 73 74 61 |d";.X#en|dif.Xsta|
|00000720| 74 69 63 20 63 68 61 72 | 20 2a 6d 65 6d 6d 73 67 |tic char| *memmsg|
|00000730| 20 3d 20 22 4e 6f 74 20 | 65 6e 6f 75 67 68 20 6d | = "Not |enough m|
|00000740| 65 6d 6f 72 79 22 3b 0a | 58 0a 58 0a 58 2f 2a 0a |emory";.|X.X.X/*.|
|00000750| 58 20 2a 20 43 72 65 61 | 74 65 20 61 6e 6f 74 68 |X * Crea|te anoth|
|00000760| 65 72 20 63 6f 70 79 20 | 6f 66 20 61 20 6e 75 6d |er copy |of a num|
|00000770| 62 65 72 2e 0a 58 20 2a | 09 71 32 20 3d 20 71 63 |ber..X *|.q2 = qc|
|00000780| 6f 70 79 28 71 31 29 3b | 0a 58 20 2a 2f 0a 58 4e |opy(q1);|.X */.XN|
|00000790| 55 4d 42 45 52 20 2a 0a | 58 71 63 6f 70 79 28 71 |UMBER *.|Xqcopy(q|
|000007a0| 29 0a 58 09 72 65 67 69 | 73 74 65 72 20 4e 55 4d |).X.regi|ster NUM|
|000007b0| 42 45 52 20 2a 71 3b 0a | 58 7b 0a 58 09 72 65 67 |BER *q;.|X{.X.reg|
|000007c0| 69 73 74 65 72 20 4e 55 | 4d 42 45 52 20 2a 72 3b |ister NU|MBER *r;|
|000007d0| 0a 58 0a 58 09 72 20 3d | 20 71 61 6c 6c 6f 63 28 |.X.X.r =| qalloc(|
|000007e0| 29 3b 0a 58 09 72 2d 3e | 6e 75 6d 2e 73 69 67 6e |);.X.r->|num.sign|
|000007f0| 20 3d 20 71 2d 3e 6e 75 | 6d 2e 73 69 67 6e 3b 0a | = q->nu|m.sign;.|
|00000800| 58 09 69 66 20 28 21 69 | 73 75 6e 69 74 28 71 2d |X.if (!i|sunit(q-|
|00000810| 3e 6e 75 6d 29 29 20 7b | 0a 58 09 09 72 2d 3e 6e |>num)) {|.X..r->n|
|00000820| 75 6d 2e 6c 65 6e 20 3d | 20 71 2d 3e 6e 75 6d 2e |um.len =| q->num.|
|00000830| 6c 65 6e 3b 0a 58 09 09 | 72 2d 3e 6e 75 6d 2e 76 |len;.X..|r->num.v|
|00000840| 20 3d 20 61 6c 6c 6f 63 | 28 72 2d 3e 6e 75 6d 2e | = alloc|(r->num.|
|00000850| 6c 65 6e 29 3b 0a 58 09 | 09 63 6f 70 79 76 61 6c |len);.X.|.copyval|
|00000860| 28 71 2d 3e 6e 75 6d 2c | 20 72 2d 3e 6e 75 6d 29 |(q->num,| r->num)|
|00000870| 3b 0a 58 09 7d 0a 58 09 | 69 66 20 28 21 69 73 75 |;.X.}.X.|if (!isu|
|00000880| 6e 69 74 28 71 2d 3e 64 | 65 6e 29 29 20 7b 0a 58 |nit(q->d|en)) {.X|
|00000890| 09 09 72 2d 3e 64 65 6e | 2e 6c 65 6e 20 3d 20 71 |..r->den|.len = q|
|000008a0| 2d 3e 64 65 6e 2e 6c 65 | 6e 3b 0a 58 09 09 72 2d |->den.le|n;.X..r-|
|000008b0| 3e 64 65 6e 2e 76 20 3d | 20 61 6c 6c 6f 63 28 72 |>den.v =| alloc(r|
|000008c0| 2d 3e 64 65 6e 2e 6c 65 | 6e 29 3b 0a 58 09 09 63 |->den.le|n);.X..c|
|000008d0| 6f 70 79 76 61 6c 28 71 | 2d 3e 64 65 6e 2c 20 72 |opyval(q|->den, r|
|000008e0| 2d 3e 64 65 6e 29 3b 0a | 58 09 7d 0a 58 09 72 65 |->den);.|X.}.X.re|
|000008f0| 74 75 72 6e 20 72 3b 0a | 58 7d 0a 58 0a 58 0a 58 |turn r;.|X}.X.X.X|
|00000900| 2f 2a 0a 58 20 2a 20 43 | 6f 6e 76 65 72 74 20 61 |/*.X * C|onvert a|
|00000910| 20 6e 75 6d 62 65 72 20 | 74 6f 20 61 20 6e 6f 72 | number |to a nor|
|00000920| 6d 61 6c 20 69 6e 74 65 | 67 65 72 2e 0a 58 20 2a |mal inte|ger..X *|
|00000930| 09 69 20 3d 20 71 74 6f | 69 28 71 29 3b 0a 58 20 |.i = qto|i(q);.X |
|00000940| 2a 2f 0a 58 6c 6f 6e 67 | 0a 58 71 74 6f 69 28 71 |*/.Xlong|.Xqtoi(q|
|00000950| 29 0a 58 09 72 65 67 69 | 73 74 65 72 20 4e 55 4d |).X.regi|ster NUM|
|00000960| 42 45 52 20 2a 71 3b 0a | 58 7b 0a 58 09 6c 6f 6e |BER *q;.|X{.X.lon|
|00000970| 67 20 69 3b 0a 58 09 5a | 56 41 4c 55 45 20 72 65 |g i;.X.Z|VALUE re|
|00000980| 73 3b 0a 58 0a 58 09 69 | 66 20 28 71 69 73 69 6e |s;.X.X.i|f (qisin|
|00000990| 74 28 71 29 29 0a 58 09 | 09 72 65 74 75 72 6e 20 |t(q)).X.|.return |
|000009a0| 7a 74 6f 69 28 71 2d 3e | 6e 75 6d 29 3b 0a 58 09 |ztoi(q->|num);.X.|
|000009b0| 7a 71 75 6f 28 71 2d 3e | 6e 75 6d 2c 20 71 2d 3e |zquo(q->|num, q->|
|000009c0| 64 65 6e 2c 20 26 72 65 | 73 29 3b 0a 58 09 69 20 |den, &re|s);.X.i |
|000009d0| 3d 20 7a 74 6f 69 28 72 | 65 73 29 3b 0a 58 09 66 |= ztoi(r|es);.X.f|
|000009e0| 72 65 65 68 28 72 65 73 | 2e 76 29 3b 0a 58 09 72 |reeh(res|.v);.X.r|
|000009f0| 65 74 75 72 6e 20 69 3b | 0a 58 7d 0a 58 0a 58 0a |eturn i;|.X}.X.X.|
|00000a00| 58 2f 2a 0a 58 20 2a 20 | 43 6f 6e 76 65 72 74 20 |X/*.X * |Convert |
|00000a10| 61 20 6e 6f 72 6d 61 6c | 20 69 6e 74 65 67 65 72 |a normal| integer|
|00000a20| 20 69 6e 74 6f 20 61 20 | 6e 75 6d 62 65 72 2e 0a | into a |number..|
|00000a30| 58 20 2a 09 71 20 3d 20 | 69 74 6f 71 28 69 29 3b |X *.q = |itoq(i);|
|00000a40| 0a 58 20 2a 2f 0a 58 4e | 55 4d 42 45 52 20 2a 0a |.X */.XN|UMBER *.|
|00000a50| 58 69 74 6f 71 28 69 29 | 0a 58 09 6c 6f 6e 67 20 |Xitoq(i)|.X.long |
|00000a60| 69 3b 0a 58 7b 0a 58 09 | 72 65 67 69 73 74 65 72 |i;.X{.X.|register|
|00000a70| 20 4e 55 4d 42 45 52 20 | 2a 71 3b 0a 58 0a 58 09 | NUMBER |*q;.X.X.|
|00000a80| 69 66 20 28 28 69 20 3e | 3d 20 2d 31 29 20 26 26 |if ((i >|= -1) &&|
|00000a90| 20 28 69 20 3c 3d 20 31 | 30 29 29 20 7b 0a 58 09 | (i <= 1|0)) {.X.|
|00000aa0| 09 73 77 69 74 63 68 20 | 28 28 69 6e 74 29 20 69 |.switch |((int) i|
|00000ab0| 29 20 7b 0a 58 09 09 09 | 63 61 73 65 20 30 3a 20 |) {.X...|case 0: |
|00000ac0| 71 20 3d 20 26 5f 71 7a | 65 72 6f 5f 3b 20 62 72 |q = &_qz|ero_; br|
|00000ad0| 65 61 6b 3b 0a 58 09 09 | 09 63 61 73 65 20 31 3a |eak;.X..|.case 1:|
|00000ae0| 20 71 20 3d 20 26 5f 71 | 6f 6e 65 5f 3b 20 62 72 | q = &_q|one_; br|
|00000af0| 65 61 6b 3b 0a 58 09 09 | 09 63 61 73 65 20 32 3a |eak;.X..|.case 2:|
|00000b00| 20 71 20 3d 20 26 5f 71 | 74 77 6f 5f 3b 20 62 72 | q = &_q|two_; br|
|00000b10| 65 61 6b 3b 0a 58 09 09 | 09 63 61 73 65 20 31 30 |eak;.X..|.case 10|
|00000b20| 3a 20 71 20 3d 20 26 5f | 71 74 65 6e 5f 3b 20 62 |: q = &_|qten_; b|
|00000b30| 72 65 61 6b 3b 0a 58 09 | 09 09 63 61 73 65 20 2d |reak;.X.|..case -|
|00000b40| 31 3a 20 71 20 3d 20 26 | 5f 71 6e 65 67 6f 6e 65 |1: q = &|_qnegone|
|00000b50| 5f 3b 20 62 72 65 61 6b | 3b 0a 58 09 09 09 64 65 |_; break|;.X...de|
|00000b60| 66 61 75 6c 74 3a 20 71 | 20 3d 20 4e 55 4c 4c 3b |fault: q| = NULL;|
|00000b70| 0a 58 09 09 7d 0a 58 09 | 09 69 66 20 28 71 29 0a |.X..}.X.|.if (q).|
|00000b80| 58 09 09 09 72 65 74 75 | 72 6e 20 71 6c 69 6e 6b |X...retu|rn qlink|
|00000b90| 28 71 29 3b 0a 58 09 7d | 0a 58 09 71 20 3d 20 71 |(q);.X.}|.X.q = q|
|00000ba0| 61 6c 6c 6f 63 28 29 3b | 0a 58 09 69 74 6f 7a 28 |alloc();|.X.itoz(|
|00000bb0| 69 2c 20 26 71 2d 3e 6e | 75 6d 29 3b 0a 58 09 72 |i, &q->n|um);.X.r|
|00000bc0| 65 74 75 72 6e 20 71 3b | 0a 58 7d 0a 58 0a 58 0a |eturn q;|.X}.X.X.|
|00000bd0| 58 2f 2a 0a 58 20 2a 20 | 43 72 65 61 74 65 20 61 |X/*.X * |Create a|
|00000be0| 20 6e 75 6d 62 65 72 20 | 66 72 6f 6d 20 74 68 65 | number |from the|
|00000bf0| 20 67 69 76 65 6e 20 69 | 6e 74 65 67 72 61 6c 20 | given i|ntegral |
|00000c00| 6e 75 6d 65 72 61 74 6f | 72 20 61 6e 64 20 64 65 |numerato|r and de|
|00000c10| 6e 6f 6d 69 6e 61 74 6f | 72 2e 0a 58 20 2a 09 71 |nominato|r..X *.q|
|00000c20| 20 3d 20 69 69 74 6f 71 | 28 69 6e 75 6d 2c 20 69 | = iitoq|(inum, i|
|00000c30| 64 65 6e 29 3b 0a 58 20 | 2a 2f 0a 58 4e 55 4d 42 |den);.X |*/.XNUMB|
|00000c40| 45 52 20 2a 0a 58 69 69 | 74 6f 71 28 69 6e 75 6d |ER *.Xii|toq(inum|
|00000c50| 2c 20 69 64 65 6e 29 0a | 58 09 6c 6f 6e 67 20 69 |, iden).|X.long i|
|00000c60| 6e 75 6d 2c 20 69 64 65 | 6e 3b 0a 58 7b 0a 58 09 |num, ide|n;.X{.X.|
|00000c70| 72 65 67 69 73 74 65 72 | 20 4e 55 4d 42 45 52 20 |register| NUMBER |
|00000c80| 2a 71 3b 0a 58 09 6c 6f | 6e 67 20 64 3b 0a 58 09 |*q;.X.lo|ng d;.X.|
|00000c90| 42 4f 4f 4c 20 73 69 67 | 6e 3b 0a 58 0a 58 09 69 |BOOL sig|n;.X.X.i|
|00000ca0| 66 20 28 69 64 65 6e 20 | 3d 3d 20 30 29 0a 58 09 |f (iden |== 0).X.|
|00000cb0| 09 65 72 72 6f 72 28 22 | 44 69 76 69 73 69 6f 6e |.error("|Division|
|00000cc0| 20 62 79 20 7a 65 72 6f | 22 29 3b 0a 58 09 69 66 | by zero|");.X.if|
|00000cd0| 20 28 69 6e 75 6d 20 3d | 3d 20 30 29 0a 58 09 09 | (inum =|= 0).X..|
|00000ce0| 72 65 74 75 72 6e 20 71 | 6c 69 6e 6b 28 26 5f 71 |return q|link(&_q|
|00000cf0| 7a 65 72 6f 5f 29 3b 0a | 58 09 73 69 67 6e 20 3d |zero_);.|X.sign =|
|00000d00| 20 30 3b 0a 58 09 69 66 | 20 28 69 6e 75 6d 20 3c | 0;.X.if| (inum <|
|00000d10| 20 30 29 20 7b 0a 58 09 | 09 73 69 67 6e 20 3d 20 | 0) {.X.|.sign = |
|00000d20| 31 3b 0a 58 09 09 69 6e | 75 6d 20 3d 20 2d 69 6e |1;.X..in|um = -in|
|00000d30| 75 6d 3b 0a 58 09 7d 0a | 58 09 69 66 20 28 69 64 |um;.X.}.|X.if (id|
|00000d40| 65 6e 20 3c 20 30 29 20 | 7b 0a 58 09 09 73 69 67 |en < 0) |{.X..sig|
|00000d50| 6e 20 3d 20 31 20 2d 20 | 73 69 67 6e 3b 0a 58 09 |n = 1 - |sign;.X.|
|00000d60| 09 69 64 65 6e 20 3d 20 | 2d 69 64 65 6e 3b 0a 58 |.iden = |-iden;.X|
|00000d70| 09 7d 0a 58 09 64 20 3d | 20 69 69 67 63 64 28 69 |.}.X.d =| iigcd(i|
|00000d80| 6e 75 6d 2c 20 69 64 65 | 6e 29 3b 0a 58 09 69 6e |num, ide|n);.X.in|
|00000d90| 75 6d 20 2f 3d 20 64 3b | 0a 58 09 69 64 65 6e 20 |um /= d;|.X.iden |
|00000da0| 2f 3d 20 64 3b 0a 58 09 | 69 66 20 28 69 64 65 6e |/= d;.X.|if (iden|
|00000db0| 20 3d 3d 20 31 29 0a 58 | 09 09 72 65 74 75 72 6e | == 1).X|..return|
|00000dc0| 20 69 74 6f 71 28 73 69 | 67 6e 20 3f 20 2d 69 6e | itoq(si|gn ? -in|
|00000dd0| 75 6d 20 3a 20 69 6e 75 | 6d 29 3b 0a 58 09 71 20 |um : inu|m);.X.q |
|00000de0| 3d 20 71 61 6c 6c 6f 63 | 28 29 3b 0a 58 09 69 66 |= qalloc|();.X.if|
|00000df0| 20 28 69 6e 75 6d 20 21 | 3d 20 31 29 0a 58 09 09 | (inum !|= 1).X..|
|00000e00| 69 74 6f 7a 28 69 6e 75 | 6d 2c 20 26 71 2d 3e 6e |itoz(inu|m, &q->n|
|00000e10| 75 6d 29 3b 0a 58 09 69 | 74 6f 7a 28 69 64 65 6e |um);.X.i|toz(iden|
|00000e20| 2c 20 26 71 2d 3e 64 65 | 6e 29 3b 0a 58 09 71 2d |, &q->de|n);.X.q-|
|00000e30| 3e 6e 75 6d 2e 73 69 67 | 6e 20 3d 20 73 69 67 6e |>num.sig|n = sign|
|00000e40| 3b 0a 58 09 72 65 74 75 | 72 6e 20 71 3b 0a 58 7d |;.X.retu|rn q;.X}|
|00000e50| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 41 64 64 |.X.X.X/*|.X * Add|
|00000e60| 20 74 77 6f 20 6e 75 6d | 62 65 72 73 20 74 6f 20 | two num|bers to |
|00000e70| 65 61 63 68 20 6f 74 68 | 65 72 2e 0a 58 20 2a 09 |each oth|er..X *.|
|00000e80| 71 33 20 3d 20 71 61 64 | 64 28 71 31 2c 20 71 32 |q3 = qad|d(q1, q2|
|00000e90| 29 3b 0a 58 20 2a 2f 0a | 58 4e 55 4d 42 45 52 20 |);.X */.|XNUMBER |
|00000ea0| 2a 0a 58 71 61 64 64 28 | 71 31 2c 20 71 32 29 0a |*.Xqadd(|q1, q2).|
|00000eb0| 58 09 72 65 67 69 73 74 | 65 72 20 4e 55 4d 42 45 |X.regist|er NUMBE|
|00000ec0| 52 20 2a 71 31 2c 20 2a | 71 32 3b 0a 58 7b 0a 58 |R *q1, *|q2;.X{.X|
|00000ed0| 09 4e 55 4d 42 45 52 20 | 2a 72 3b 0a 58 09 5a 56 |.NUMBER |*r;.X.ZV|
|00000ee0| 41 4c 55 45 20 74 31 2c | 20 74 32 2c 20 74 65 6d |ALUE t1,| t2, tem|
|00000ef0| 70 2c 20 64 31 2c 20 64 | 32 2c 20 76 70 64 31 2c |p, d1, d|2, vpd1,|
|00000f00| 20 75 70 64 31 3b 0a 58 | 0a 58 09 72 20 3d 20 71 | upd1;.X|.X.r = q|
|00000f10| 61 6c 6c 6f 63 28 29 3b | 0a 58 09 2f 2a 0a 58 09 |alloc();|.X./*.X.|
|00000f20| 20 2a 20 49 66 20 65 69 | 74 68 65 72 20 6e 75 6d | * If ei|ther num|
|00000f30| 62 65 72 20 69 73 20 61 | 6e 20 69 6e 74 65 67 65 |ber is a|n intege|
|00000f40| 72 2c 20 74 68 65 6e 20 | 74 68 65 20 72 65 73 75 |r, then |the resu|
|00000f50| 6c 74 20 69 73 20 65 61 | 73 79 2e 0a 58 09 20 2a |lt is ea|sy..X. *|
|00000f60| 2f 0a 58 09 69 66 20 28 | 71 69 73 69 6e 74 28 71 |/.X.if (|qisint(q|
|00000f70| 31 29 20 26 26 20 71 69 | 73 69 6e 74 28 71 32 29 |1) && qi|sint(q2)|
|00000f80| 29 20 7b 0a 58 09 09 7a | 61 64 64 28 71 31 2d 3e |) {.X..z|add(q1->|
|00000f90| 6e 75 6d 2c 20 71 32 2d | 3e 6e 75 6d 2c 20 26 72 |num, q2-|>num, &r|
|00000fa0| 2d 3e 6e 75 6d 29 3b 0a | 58 09 09 72 65 74 75 72 |->num);.|X..retur|
|00000fb0| 6e 20 72 3b 0a 58 09 7d | 0a 58 09 69 66 20 28 71 |n r;.X.}|.X.if (q|
|00000fc0| 69 73 69 6e 74 28 71 32 | 29 29 20 7b 0a 58 09 09 |isint(q2|)) {.X..|
|00000fd0| 7a 6d 75 6c 28 71 31 2d | 3e 64 65 6e 2c 20 71 32 |zmul(q1-|>den, q2|
|00000fe0| 2d 3e 6e 75 6d 2c 20 26 | 74 65 6d 70 29 3b 0a 58 |->num, &|temp);.X|
|00000ff0| 09 09 7a 61 64 64 28 71 | 31 2d 3e 6e 75 6d 2c 20 |..zadd(q|1->num, |
|00001000| 74 65 6d 70 2c 20 26 72 | 2d 3e 6e 75 6d 29 3b 0a |temp, &r|->num);.|
|00001010| 58 09 09 66 72 65 65 68 | 28 74 65 6d 70 2e 76 29 |X..freeh|(temp.v)|
|00001020| 3b 0a 58 09 09 7a 63 6f | 70 79 28 71 31 2d 3e 64 |;.X..zco|py(q1->d|
|00001030| 65 6e 2c 20 26 72 2d 3e | 64 65 6e 29 3b 0a 58 09 |en, &r->|den);.X.|
|00001040| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 09 7d 0a 58 |.return |r;.X.}.X|
|00001050| 09 69 66 20 28 71 69 73 | 69 6e 74 28 71 31 29 29 |.if (qis|int(q1))|
|00001060| 20 7b 0a 58 09 09 7a 6d | 75 6c 28 71 32 2d 3e 64 | {.X..zm|ul(q2->d|
|00001070| 65 6e 2c 20 71 31 2d 3e | 6e 75 6d 2c 20 26 74 65 |en, q1->|num, &te|
|00001080| 6d 70 29 3b 0a 58 09 09 | 7a 61 64 64 28 71 32 2d |mp);.X..|zadd(q2-|
|00001090| 3e 6e 75 6d 2c 20 74 65 | 6d 70 2c 20 26 72 2d 3e |>num, te|mp, &r->|
|000010a0| 6e 75 6d 29 3b 0a 58 09 | 09 66 72 65 65 68 28 74 |num);.X.|.freeh(t|
|000010b0| 65 6d 70 2e 76 29 3b 0a | 58 09 09 7a 63 6f 70 79 |emp.v);.|X..zcopy|
|000010c0| 28 71 32 2d 3e 64 65 6e | 2c 20 26 72 2d 3e 64 65 |(q2->den|, &r->de|
|000010d0| 6e 29 3b 0a 58 09 09 72 | 65 74 75 72 6e 20 72 3b |n);.X..r|eturn r;|
|000010e0| 0a 58 09 7d 0a 58 09 2f | 2a 0a 58 09 20 2a 20 42 |.X.}.X./|*.X. * B|
|000010f0| 6f 74 68 20 61 72 67 75 | 6d 65 6e 74 73 20 61 72 |oth argu|ments ar|
|00001100| 65 20 74 72 75 65 20 66 | 72 61 63 74 69 6f 6e 73 |e true f|ractions|
|00001110| 2c 20 73 6f 20 77 65 20 | 6e 65 65 64 20 6d 6f 72 |, so we |need mor|
|00001120| 65 20 77 6f 72 6b 2e 0a | 58 09 20 2a 20 49 66 20 |e work..|X. * If |
|00001130| 74 68 65 20 64 65 6e 6f | 6d 69 6e 61 74 6f 72 73 |the deno|minators|
|00001140| 20 61 72 65 20 72 65 6c | 61 74 69 76 65 6c 79 20 | are rel|atively |
|00001150| 70 72 69 6d 65 2c 20 74 | 68 65 6e 20 74 68 65 20 |prime, t|hen the |
|00001160| 61 6e 73 77 65 72 20 69 | 73 20 74 68 65 0a 58 09 |answer i|s the.X.|
|00001170| 20 2a 20 73 74 72 61 69 | 67 68 74 66 6f 72 77 61 | * strai|ghtforwa|
|00001180| 72 64 20 63 72 6f 73 73 | 20 70 72 6f 64 75 63 74 |rd cross| product|
|00001190| 20 72 65 73 75 6c 74 20 | 77 69 74 68 20 6e 6f 20 | result |with no |
|000011a0| 6e 65 65 64 20 66 6f 72 | 20 72 65 64 75 63 74 69 |need for| reducti|
|000011b0| 6f 6e 2e 0a 58 09 20 2a | 2f 0a 58 09 7a 67 63 64 |on..X. *|/.X.zgcd|
|000011c0| 28 71 31 2d 3e 64 65 6e | 2c 20 71 32 2d 3e 64 65 |(q1->den|, q2->de|
|000011d0| 6e 2c 20 26 64 31 29 3b | 0a 58 09 69 66 20 28 69 |n, &d1);|.X.if (i|
|000011e0| 73 75 6e 69 74 28 64 31 | 29 29 20 7b 0a 58 09 09 |sunit(d1|)) {.X..|
|000011f0| 66 72 65 65 68 28 64 31 | 2e 76 29 3b 0a 58 09 09 |freeh(d1|.v);.X..|
|00001200| 7a 6d 75 6c 28 71 31 2d | 3e 6e 75 6d 2c 20 71 32 |zmul(q1-|>num, q2|
|00001210| 2d 3e 64 65 6e 2c 20 26 | 74 31 29 3b 0a 58 09 09 |->den, &|t1);.X..|
|00001220| 7a 6d 75 6c 28 71 31 2d | 3e 64 65 6e 2c 20 71 32 |zmul(q1-|>den, q2|
|00001230| 2d 3e 6e 75 6d 2c 20 26 | 74 32 29 3b 0a 58 09 09 |->num, &|t2);.X..|
|00001240| 7a 61 64 64 28 74 31 2c | 20 74 32 2c 20 26 72 2d |zadd(t1,| t2, &r-|
|00001250| 3e 6e 75 6d 29 3b 0a 58 | 09 09 66 72 65 65 68 28 |>num);.X|..freeh(|
|00001260| 74 31 2e 76 29 3b 0a 58 | 09 09 66 72 65 65 68 28 |t1.v);.X|..freeh(|
|00001270| 74 32 2e 76 29 3b 0a 58 | 09 09 7a 6d 75 6c 28 71 |t2.v);.X|..zmul(q|
|00001280| 31 2d 3e 64 65 6e 2c 20 | 71 32 2d 3e 64 65 6e 2c |1->den, |q2->den,|
|00001290| 20 26 72 2d 3e 64 65 6e | 29 3b 0a 58 09 09 72 65 | &r->den|);.X..re|
|000012a0| 74 75 72 6e 20 72 3b 0a | 58 09 7d 0a 58 09 2f 2a |turn r;.|X.}.X./*|
|000012b0| 0a 58 09 20 2a 20 54 68 | 65 20 63 61 6c 63 75 6c |.X. * Th|e calcul|
|000012c0| 61 74 69 6f 6e 20 69 73 | 20 6e 6f 77 20 6d 6f 72 |ation is| now mor|
|000012d0| 65 20 63 6f 6d 70 6c 69 | 63 61 74 65 64 2e 0a 58 |e compli|cated..X|
|000012e0| 09 20 2a 20 53 65 65 20 | 4b 6e 75 74 68 20 56 6f |. * See |Knuth Vo|
|000012f0| 6c 20 32 20 66 6f 72 20 | 64 65 74 61 69 6c 73 2e |l 2 for |details.|
|00001300| 0a 58 09 20 2a 2f 0a 58 | 09 7a 71 75 6f 28 71 32 |.X. */.X|.zquo(q2|
|00001310| 2d 3e 64 65 6e 2c 20 64 | 31 2c 20 26 76 70 64 31 |->den, d|1, &vpd1|
|00001320| 29 3b 0a 58 09 7a 71 75 | 6f 28 71 31 2d 3e 64 65 |);.X.zqu|o(q1->de|
|00001330| 6e 2c 20 64 31 2c 20 26 | 75 70 64 31 29 3b 0a 58 |n, d1, &|upd1);.X|
|00001340| 09 7a 6d 75 6c 28 71 31 | 2d 3e 6e 75 6d 2c 20 76 |.zmul(q1|->num, v|
|00001350| 70 64 31 2c 20 26 74 31 | 29 3b 0a 58 09 7a 6d 75 |pd1, &t1|);.X.zmu|
|00001360| 6c 28 71 32 2d 3e 6e 75 | 6d 2c 20 75 70 64 31 2c |l(q2->nu|m, upd1,|
|00001370| 20 26 74 32 29 3b 0a 58 | 09 7a 61 64 64 28 74 31 | &t2);.X|.zadd(t1|
|00001380| 2c 20 74 32 2c 20 26 74 | 65 6d 70 29 3b 0a 58 09 |, t2, &t|emp);.X.|
|00001390| 66 72 65 65 68 28 74 31 | 2e 76 29 3b 0a 58 09 66 |freeh(t1|.v);.X.f|
|000013a0| 72 65 65 68 28 74 32 2e | 76 29 3b 0a 58 09 66 72 |reeh(t2.|v);.X.fr|
|000013b0| 65 65 68 28 76 70 64 31 | 2e 76 29 3b 0a 58 09 7a |eeh(vpd1|.v);.X.z|
|000013c0| 67 63 64 28 74 65 6d 70 | 2c 20 64 31 2c 20 26 64 |gcd(temp|, d1, &d|
|000013d0| 32 29 3b 0a 58 09 66 72 | 65 65 68 28 64 31 2e 76 |2);.X.fr|eeh(d1.v|
|000013e0| 29 3b 0a 58 09 69 66 20 | 28 69 73 75 6e 69 74 28 |);.X.if |(isunit(|
|000013f0| 64 32 29 29 20 7b 0a 58 | 09 09 66 72 65 65 68 28 |d2)) {.X|..freeh(|
|00001400| 64 32 2e 76 29 3b 0a 58 | 09 09 72 2d 3e 6e 75 6d |d2.v);.X|..r->num|
|00001410| 20 3d 20 74 65 6d 70 3b | 0a 58 09 09 7a 6d 75 6c | = temp;|.X..zmul|
|00001420| 28 75 70 64 31 2c 20 71 | 32 2d 3e 64 65 6e 2c 20 |(upd1, q|2->den, |
|00001430| 26 72 2d 3e 64 65 6e 29 | 3b 0a 58 09 09 66 72 65 |&r->den)|;.X..fre|
|00001440| 65 68 28 75 70 64 31 2e | 76 29 3b 0a 58 09 09 72 |eh(upd1.|v);.X..r|
|00001450| 65 74 75 72 6e 20 72 3b | 0a 58 09 7d 0a 58 09 7a |eturn r;|.X.}.X.z|
|00001460| 71 75 6f 28 74 65 6d 70 | 2c 20 64 32 2c 20 26 72 |quo(temp|, d2, &r|
|00001470| 2d 3e 6e 75 6d 29 3b 0a | 58 09 66 72 65 65 68 28 |->num);.|X.freeh(|
|00001480| 74 65 6d 70 2e 76 29 3b | 0a 58 09 7a 71 75 6f 28 |temp.v);|.X.zquo(|
|00001490| 71 32 2d 3e 64 65 6e 2c | 20 64 32 2c 20 26 74 65 |q2->den,| d2, &te|
|000014a0| 6d 70 29 3b 0a 58 09 66 | 72 65 65 68 28 64 32 2e |mp);.X.f|reeh(d2.|
|000014b0| 76 29 3b 0a 58 09 7a 6d | 75 6c 28 74 65 6d 70 2c |v);.X.zm|ul(temp,|
|000014c0| 20 75 70 64 31 2c 20 26 | 72 2d 3e 64 65 6e 29 3b | upd1, &|r->den);|
|000014d0| 0a 58 09 66 72 65 65 68 | 28 74 65 6d 70 2e 76 29 |.X.freeh|(temp.v)|
|000014e0| 3b 0a 58 09 66 72 65 65 | 68 28 75 70 64 31 2e 76 |;.X.free|h(upd1.v|
|000014f0| 29 3b 0a 58 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |);.X.ret|urn r;.X|
|00001500| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 53 75 |}.X.X.X/|*.X * Su|
|00001510| 62 74 72 61 63 74 20 6f | 6e 65 20 6e 75 6d 62 65 |btract o|ne numbe|
|00001520| 72 20 66 72 6f 6d 20 61 | 6e 6f 74 68 65 72 2e 0a |r from a|nother..|
|00001530| 58 20 2a 09 71 33 20 3d | 20 71 73 75 62 28 71 31 |X *.q3 =| qsub(q1|
|00001540| 2c 20 71 32 29 3b 0a 58 | 20 2a 2f 0a 58 4e 55 4d |, q2);.X| */.XNUM|
|00001550| 42 45 52 20 2a 0a 58 71 | 73 75 62 28 71 31 2c 20 |BER *.Xq|sub(q1, |
|00001560| 71 32 29 0a 58 09 72 65 | 67 69 73 74 65 72 20 4e |q2).X.re|gister N|
|00001570| 55 4d 42 45 52 20 2a 71 | 31 2c 20 2a 71 32 3b 0a |UMBER *q|1, *q2;.|
|00001580| 58 7b 0a 58 09 4e 55 4d | 42 45 52 20 74 2c 20 2a |X{.X.NUM|BER t, *|
|00001590| 72 3b 0a 58 0a 58 09 69 | 66 20 28 71 31 20 3d 3d |r;.X.X.i|f (q1 ==|
|000015a0| 20 71 32 29 0a 58 09 09 | 72 65 74 75 72 6e 20 71 | q2).X..|return q|
|000015b0| 6c 69 6e 6b 28 26 5f 71 | 7a 65 72 6f 5f 29 3b 0a |link(&_q|zero_);.|
|000015c0| 58 09 69 66 20 28 71 69 | 73 69 6e 74 28 71 31 29 |X.if (qi|sint(q1)|
|000015d0| 20 26 26 20 71 69 73 69 | 6e 74 28 71 32 29 29 20 | && qisi|nt(q2)) |
|000015e0| 7b 0a 58 09 09 72 20 3d | 20 71 61 6c 6c 6f 63 28 |{.X..r =| qalloc(|
|000015f0| 29 3b 0a 58 09 09 7a 73 | 75 62 28 71 31 2d 3e 6e |);.X..zs|ub(q1->n|
|00001600| 75 6d 2c 20 71 32 2d 3e | 6e 75 6d 2c 20 26 72 2d |um, q2->|num, &r-|
|00001610| 3e 6e 75 6d 29 3b 0a 58 | 09 09 72 65 74 75 72 6e |>num);.X|..return|
|00001620| 20 72 3b 0a 58 09 7d 0a | 58 09 74 20 3d 20 2a 71 | r;.X.}.|X.t = *q|
|00001630| 32 3b 0a 58 09 74 2e 6e | 75 6d 2e 73 69 67 6e 20 |2;.X.t.n|um.sign |
|00001640| 3d 20 21 74 2e 6e 75 6d | 2e 73 69 67 6e 3b 0a 58 |= !t.num|.sign;.X|
|00001650| 09 72 65 74 75 72 6e 20 | 71 61 64 64 28 71 31 2c |.return |qadd(q1,|
|00001660| 20 26 74 29 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a | &t);.X}|.X.X.X/*|
|00001670| 0a 58 20 2a 20 49 6e 63 | 72 65 6d 65 6e 74 20 61 |.X * Inc|rement a|
|00001680| 20 6e 75 6d 62 65 72 20 | 62 79 20 6f 6e 65 2e 0a | number |by one..|
|00001690| 58 20 2a 2f 0a 58 4e 55 | 4d 42 45 52 20 2a 0a 58 |X */.XNU|MBER *.X|
|000016a0| 71 69 6e 63 28 71 29 0a | 58 09 4e 55 4d 42 45 52 |qinc(q).|X.NUMBER|
|000016b0| 20 2a 71 3b 0a 58 7b 0a | 58 09 4e 55 4d 42 45 52 | *q;.X{.|X.NUMBER|
|000016c0| 20 2a 72 3b 0a 58 0a 58 | 09 72 20 3d 20 71 61 6c | *r;.X.X|.r = qal|
|000016d0| 6c 6f 63 28 29 3b 0a 58 | 09 69 66 20 28 71 69 73 |loc();.X|.if (qis|
|000016e0| 69 6e 74 28 71 29 29 20 | 7b 0a 58 09 09 7a 61 64 |int(q)) |{.X..zad|
|000016f0| 64 28 71 2d 3e 6e 75 6d | 2c 20 5f 6f 6e 65 5f 2c |d(q->num|, _one_,|
|00001700| 20 26 72 2d 3e 6e 75 6d | 29 3b 0a 58 09 09 72 65 | &r->num|);.X..re|
|00001710| 74 75 72 6e 20 72 3b 0a | 58 09 7d 0a 58 09 7a 61 |turn r;.|X.}.X.za|
|00001720| 64 64 28 71 2d 3e 6e 75 | 6d 2c 20 71 2d 3e 64 65 |dd(q->nu|m, q->de|
|00001730| 6e 2c 20 26 72 2d 3e 6e | 75 6d 29 3b 0a 58 09 7a |n, &r->n|um);.X.z|
|00001740| 63 6f 70 79 28 71 2d 3e | 64 65 6e 2c 20 26 72 2d |copy(q->|den, &r-|
|00001750| 3e 64 65 6e 29 3b 0a 58 | 09 72 65 74 75 72 6e 20 |>den);.X|.return |
|00001760| 72 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |r;.X}.X.|X.X/*.X |
|00001770| 2a 20 44 65 63 72 65 6d | 65 6e 74 20 61 20 6e 75 |* Decrem|ent a nu|
|00001780| 6d 62 65 72 20 62 79 20 | 6f 6e 65 2e 0a 58 20 2a |mber by |one..X *|
|00001790| 2f 0a 58 4e 55 4d 42 45 | 52 20 2a 0a 58 71 64 65 |/.XNUMBE|R *.Xqde|
|000017a0| 63 28 71 29 0a 58 09 4e | 55 4d 42 45 52 20 2a 71 |c(q).X.N|UMBER *q|
|000017b0| 3b 0a 58 7b 0a 58 09 4e | 55 4d 42 45 52 20 2a 72 |;.X{.X.N|UMBER *r|
|000017c0| 3b 0a 58 0a 58 09 72 20 | 3d 20 71 61 6c 6c 6f 63 |;.X.X.r |= qalloc|
|000017d0| 28 29 3b 0a 58 09 69 66 | 20 28 71 69 73 69 6e 74 |();.X.if| (qisint|
|000017e0| 28 71 29 29 20 7b 0a 58 | 09 09 7a 73 75 62 28 71 |(q)) {.X|..zsub(q|
|000017f0| 2d 3e 6e 75 6d 2c 20 5f | 6f 6e 65 5f 2c 20 26 72 |->num, _|one_, &r|
|00001800| 2d 3e 6e 75 6d 29 3b 0a | 58 09 09 72 65 74 75 72 |->num);.|X..retur|
|00001810| 6e 20 72 3b 0a 58 09 7d | 0a 58 09 7a 73 75 62 28 |n r;.X.}|.X.zsub(|
|00001820| 71 2d 3e 6e 75 6d 2c 20 | 71 2d 3e 64 65 6e 2c 20 |q->num, |q->den, |
|00001830| 26 72 2d 3e 6e 75 6d 29 | 3b 0a 58 09 7a 63 6f 70 |&r->num)|;.X.zcop|
|00001840| 79 28 71 2d 3e 64 65 6e | 2c 20 26 72 2d 3e 64 65 |y(q->den|, &r->de|
|00001850| 6e 29 3b 0a 58 09 72 65 | 74 75 72 6e 20 72 3b 0a |n);.X.re|turn r;.|
|00001860| 58 7d 0a 58 0a 58 0a 58 | 23 69 66 64 65 66 20 43 |X}.X.X.X|#ifdef C|
|00001870| 4f 44 45 0a 58 2f 2a 0a | 58 20 2a 20 41 64 64 20 |ODE.X/*.|X * Add |
|00001880| 61 20 6e 6f 72 6d 61 6c | 20 73 6d 61 6c 6c 20 69 |a normal| small i|
|00001890| 6e 74 65 67 65 72 20 76 | 61 6c 75 65 20 74 6f 20 |nteger v|alue to |
|000018a0| 61 6e 20 61 72 62 69 74 | 72 61 72 79 20 6e 75 6d |an arbit|rary num|
|000018b0| 62 65 72 2e 0a 58 20 2a | 2f 0a 58 4e 55 4d 42 45 |ber..X *|/.XNUMBE|
|000018c0| 52 20 2a 0a 58 71 61 64 | 64 69 28 71 31 2c 20 6e |R *.Xqad|di(q1, n|
|000018d0| 29 0a 58 09 4e 55 4d 42 | 45 52 20 2a 71 31 3b 0a |).X.NUMB|ER *q1;.|
|000018e0| 58 09 6c 6f 6e 67 20 6e | 3b 0a 58 7b 0a 58 09 4e |X.long n|;.X{.X.N|
|000018f0| 55 4d 42 45 52 20 61 64 | 64 6e 75 6d 3b 09 09 2f |UMBER ad|dnum;../|
|00001900| 2a 20 74 65 6d 70 6f 72 | 61 72 79 20 6e 75 6d 62 |* tempor|ary numb|
|00001910| 65 72 20 2a 2f 0a 58 09 | 48 41 4c 46 20 61 64 64 |er */.X.|HALF add|
|00001920| 76 61 6c 5b 32 5d 3b 09 | 09 2f 2a 20 76 61 6c 75 |val[2];.|./* valu|
|00001930| 65 20 6f 66 20 73 6d 61 | 6c 6c 20 6e 75 6d 62 65 |e of sma|ll numbe|
|00001940| 72 20 2a 2f 0a 58 09 42 | 4f 4f 4c 20 6e 65 67 3b |r */.X.B|OOL neg;|
|00001950| 09 09 2f 2a 20 54 52 55 | 45 20 69 66 20 6e 75 6d |../* TRU|E if num|
|00001960| 62 65 72 20 69 73 20 6e | 65 67 20 2a 2f 0a 58 0a |ber is n|eg */.X.|
|00001970| 58 09 69 66 20 28 6e 20 | 3d 3d 20 30 29 0a 58 09 |X.if (n |== 0).X.|
|00001980| 09 72 65 74 75 72 6e 20 | 71 6c 69 6e 6b 28 71 31 |.return |qlink(q1|
|00001990| 29 3b 0a 58 09 69 66 20 | 28 6e 20 3d 3d 20 31 29 |);.X.if |(n == 1)|
|000019a0| 0a 58 09 09 72 65 74 75 | 72 6e 20 71 69 6e 63 28 |.X..retu|rn qinc(|
|000019b0| 71 31 29 3b 0a 58 09 69 | 66 20 28 6e 20 3d 3d 20 |q1);.X.i|f (n == |
|000019c0| 2d 31 29 0a 58 09 09 72 | 65 74 75 72 6e 20 71 64 |-1).X..r|eturn qd|
|000019d0| 65 63 28 71 31 29 3b 0a | 58 09 69 66 20 28 71 69 |ec(q1);.|X.if (qi|
|000019e0| 73 7a 65 72 6f 28 71 31 | 29 29 0a 58 09 09 72 65 |szero(q1|)).X..re|
|000019f0| 74 75 72 6e 20 69 74 6f | 71 28 6e 29 3b 0a 58 09 |turn ito|q(n);.X.|
|00001a00| 61 64 64 6e 75 6d 2e 6e | 75 6d 2e 73 69 67 6e 20 |addnum.n|um.sign |
|00001a10| 3d 20 30 3b 0a 58 09 61 | 64 64 6e 75 6d 2e 6e 75 |= 0;.X.a|ddnum.nu|
|00001a20| 6d 2e 6c 65 6e 20 3d 20 | 31 3b 0a 58 09 61 64 64 |m.len = |1;.X.add|
|00001a30| 6e 75 6d 2e 6e 75 6d 2e | 76 20 3d 20 61 64 64 76 |num.num.|v = addv|
|00001a40| 61 6c 3b 0a 58 09 61 64 | 64 6e 75 6d 2e 64 65 6e |al;.X.ad|dnum.den|
|00001a50| 20 3d 20 5f 6f 6e 65 5f | 3b 0a 58 09 6e 65 67 20 | = _one_|;.X.neg |
|00001a60| 3d 20 28 6e 20 3c 20 30 | 29 3b 0a 58 09 69 66 20 |= (n < 0|);.X.if |
|00001a70| 28 6e 65 67 29 0a 58 09 | 09 6e 20 3d 20 2d 6e 3b |(neg).X.|.n = -n;|
|00001a80| 0a 58 09 61 64 64 76 61 | 6c 5b 30 5d 20 3d 20 28 |.X.addva|l[0] = (|
|00001a90| 48 41 4c 46 29 20 6e 3b | 0a 58 09 6e 20 3d 20 28 |HALF) n;|.X.n = (|
|00001aa0| 28 28 46 55 4c 4c 29 20 | 6e 29 20 3e 3e 20 42 41 |((FULL) |n) >> BA|
|00001ab0| 53 45 42 29 3b 0a 58 09 | 69 66 20 28 6e 29 20 7b |SEB);.X.|if (n) {|
|00001ac0| 0a 58 09 09 61 64 64 76 | 61 6c 5b 31 5d 20 3d 20 |.X..addv|al[1] = |
|00001ad0| 28 48 41 4c 46 29 20 6e | 3b 0a 58 09 09 61 64 64 |(HALF) n|;.X..add|
|00001ae0| 6e 75 6d 2e 6e 75 6d 2e | 6c 65 6e 20 3d 20 32 3b |num.num.|len = 2;|
|00001af0| 0a 58 09 7d 0a 58 09 69 | 66 20 28 6e 65 67 29 0a |.X.}.X.i|f (neg).|
|00001b00| 58 09 09 72 65 74 75 72 | 6e 20 71 73 75 62 28 71 |X..retur|n qsub(q|
|00001b10| 31 2c 20 26 61 64 64 6e | 75 6d 29 3b 0a 58 09 65 |1, &addn|um);.X.e|
|00001b20| 6c 73 65 0a 58 09 09 72 | 65 74 75 72 6e 20 71 61 |lse.X..r|eturn qa|
|00001b30| 64 64 28 71 31 2c 20 26 | 61 64 64 6e 75 6d 29 3b |dd(q1, &|addnum);|
|00001b40| 0a 58 7d 0a 58 23 65 6e | 64 69 66 0a 58 0a 58 0a |.X}.X#en|dif.X.X.|
|00001b50| 58 2f 2a 0a 58 20 2a 20 | 4d 75 6c 74 69 70 6c 79 |X/*.X * |Multiply|
|00001b60| 20 74 77 6f 20 6e 75 6d | 62 65 72 73 2e 0a 58 20 | two num|bers..X |
|00001b70| 2a 09 71 33 20 3d 20 71 | 6d 75 6c 28 71 31 2c 20 |*.q3 = q|mul(q1, |
|00001b80| 71 32 29 3b 0a 58 20 2a | 2f 0a 58 4e 55 4d 42 45 |q2);.X *|/.XNUMBE|
|00001b90| 52 20 2a 0a 58 71 6d 75 | 6c 28 71 31 2c 20 71 32 |R *.Xqmu|l(q1, q2|
|00001ba0| 29 0a 58 09 72 65 67 69 | 73 74 65 72 20 4e 55 4d |).X.regi|ster NUM|
|00001bb0| 42 45 52 20 2a 71 31 2c | 20 2a 71 32 3b 0a 58 7b |BER *q1,| *q2;.X{|
|00001bc0| 0a 58 09 4e 55 4d 42 45 | 52 20 2a 72 3b 09 09 09 |.X.NUMBE|R *r;...|
|00001bd0| 2f 2a 20 72 65 74 75 72 | 6e 65 64 20 76 61 6c 75 |/* retur|ned valu|
|00001be0| 65 20 2a 2f 0a 58 09 5a | 56 41 4c 55 45 20 6e 31 |e */.X.Z|VALUE n1|
|00001bf0| 2c 20 6e 32 2c 20 64 31 | 2c 20 64 32 3b 09 09 2f |, n2, d1|, d2;../|
|00001c00| 2a 20 6e 75 6d 65 72 61 | 74 6f 72 73 20 61 6e 64 |* numera|tors and|
|00001c10| 20 64 65 6e 6f 6d 69 6e | 61 74 6f 72 73 20 2a 2f | denomin|ators */|
|00001c20| 0a 58 09 5a 56 41 4c 55 | 45 20 74 6d 70 3b 0a 58 |.X.ZVALU|E tmp;.X|
|00001c30| 0a 58 09 69 66 20 28 71 | 69 73 7a 65 72 6f 28 71 |.X.if (q|iszero(q|
|00001c40| 31 29 20 7c 7c 20 71 69 | 73 7a 65 72 6f 28 71 32 |1) || qi|szero(q2|
|00001c50| 29 29 0a 58 09 09 72 65 | 74 75 72 6e 20 71 6c 69 |)).X..re|turn qli|
|00001c60| 6e 6b 28 26 5f 71 7a 65 | 72 6f 5f 29 3b 0a 58 09 |nk(&_qze|ro_);.X.|
|00001c70| 69 66 20 28 71 69 73 6f | 6e 65 28 71 31 29 29 0a |if (qiso|ne(q1)).|
|00001c80| 58 09 09 72 65 74 75 72 | 6e 20 71 6c 69 6e 6b 28 |X..retur|n qlink(|
|00001c90| 71 32 29 3b 0a 58 09 69 | 66 20 28 71 69 73 6f 6e |q2);.X.i|f (qison|
|00001ca0| 65 28 71 32 29 29 0a 58 | 09 09 72 65 74 75 72 6e |e(q2)).X|..return|
|00001cb0| 20 71 6c 69 6e 6b 28 71 | 31 29 3b 0a 58 09 69 66 | qlink(q|1);.X.if|
|00001cc0| 20 28 71 69 73 69 6e 74 | 28 71 31 29 20 26 26 20 | (qisint|(q1) && |
|00001cd0| 71 69 73 69 6e 74 28 71 | 32 29 29 20 7b 09 2f 2a |qisint(q|2)) {./*|
|00001ce0| 20 65 61 73 79 20 72 65 | 73 75 6c 74 73 20 69 66 | easy re|sults if|
|00001cf0| 20 69 6e 74 65 67 65 72 | 73 20 2a 2f 0a 58 09 09 | integer|s */.X..|
|00001d00| 72 20 3d 20 71 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |r = qall|oc();.X.|
|00001d10| 09 7a 6d 75 6c 28 71 31 | 2d 3e 6e 75 6d 2c 20 71 |.zmul(q1|->num, q|
|00001d20| 32 2d 3e 6e 75 6d 2c 20 | 26 72 2d 3e 6e 75 6d 29 |2->num, |&r->num)|
|00001d30| 3b 0a 58 09 09 72 65 74 | 75 72 6e 20 72 3b 0a 58 |;.X..ret|urn r;.X|
|00001d40| 09 7d 0a 58 09 6e 31 20 | 3d 20 71 31 2d 3e 6e 75 |.}.X.n1 |= q1->nu|
|00001d50| 6d 3b 0a 58 09 6e 32 20 | 3d 20 71 32 2d 3e 6e 75 |m;.X.n2 |= q2->nu|
|00001d60| 6d 3b 0a 58 09 64 31 20 | 3d 20 71 31 2d 3e 64 65 |m;.X.d1 |= q1->de|
|00001d70| 6e 3b 0a 58 09 64 32 20 | 3d 20 71 32 2d 3e 64 65 |n;.X.d2 |= q2->de|
|00001d80| 6e 3b 0a 58 09 69 66 20 | 28 69 73 7a 65 72 6f 28 |n;.X.if |(iszero(|
|00001d90| 64 31 29 20 7c 7c 20 69 | 73 7a 65 72 6f 28 64 32 |d1) || i|szero(d2|
|00001da0| 29 29 0a 58 09 09 65 72 | 72 6f 72 28 22 44 69 76 |)).X..er|ror("Div|
|00001db0| 69 73 69 6f 6e 20 62 79 | 20 7a 65 72 6f 22 29 3b |ision by| zero");|
|00001dc0| 0a 58 09 69 66 20 28 69 | 73 7a 65 72 6f 28 6e 31 |.X.if (i|szero(n1|
|00001dd0| 29 20 7c 7c 20 69 73 7a | 65 72 6f 28 6e 32 29 29 |) || isz|ero(n2))|
|00001de0| 0a 58 09 09 72 65 74 75 | 72 6e 20 71 6c 69 6e 6b |.X..retu|rn qlink|
|00001df0| 28 26 5f 71 7a 65 72 6f | 5f 29 3b 0a 58 09 69 66 |(&_qzero|_);.X.if|
|00001e00| 20 28 21 69 73 75 6e 69 | 74 28 6e 31 29 20 26 26 | (!isuni|t(n1) &&|
|00001e10| 20 21 69 73 75 6e 69 74 | 28 64 32 29 29 20 7b 09 | !isunit|(d2)) {.|
|00001e20| 2f 2a 20 70 6f 73 73 69 | 62 6c 79 20 72 65 64 75 |/* possi|bly redu|
|00001e30| 63 65 20 2a 2f 0a 58 09 | 09 7a 67 63 64 28 6e 31 |ce */.X.|.zgcd(n1|
|00001e40| 2c 20 64 32 2c 20 26 74 | 6d 70 29 3b 0a 58 09 09 |, d2, &t|mp);.X..|
|00001e50| 69 66 20 28 21 69 73 75 | 6e 69 74 28 74 6d 70 29 |if (!isu|nit(tmp)|
|00001e60| 29 20 7b 0a 58 09 09 09 | 7a 71 75 6f 28 71 31 2d |) {.X...|zquo(q1-|
|00001e70| 3e 6e 75 6d 2c 20 74 6d | 70 2c 20 26 6e 31 29 3b |>num, tm|p, &n1);|
|00001e80| 0a 58 09 09 09 7a 71 75 | 6f 28 71 32 2d 3e 64 65 |.X...zqu|o(q2->de|
|00001e90| 6e 2c 20 74 6d 70 2c 20 | 26 64 32 29 3b 0a 58 09 |n, tmp, |&d2);.X.|
|00001ea0| 09 7d 0a 58 09 09 66 72 | 65 65 68 28 74 6d 70 2e |.}.X..fr|eeh(tmp.|
|00001eb0| 76 29 3b 0a 58 09 7d 0a | 58 09 69 66 20 28 21 69 |v);.X.}.|X.if (!i|
|00001ec0| 73 75 6e 69 74 28 6e 32 | 29 20 26 26 20 21 69 73 |sunit(n2|) && !is|
|00001ed0| 75 6e 69 74 28 64 31 29 | 29 20 7b 09 2f 2a 20 61 |unit(d1)|) {./* a|
|00001ee0| 67 61 69 6e 20 70 6f 73 | 73 69 62 6c 79 20 72 65 |gain pos|sibly re|
|00001ef0| 64 75 63 65 20 2a 2f 0a | 58 09 09 7a 67 63 64 28 |duce */.|X..zgcd(|
|00001f00| 6e 32 2c 20 64 31 2c 20 | 26 74 6d 70 29 3b 0a 58 |n2, d1, |&tmp);.X|
|00001f10| 09 09 69 66 20 28 21 69 | 73 75 6e 69 74 28 74 6d |..if (!i|sunit(tm|
|00001f20| 70 29 29 20 7b 0a 58 09 | 09 09 7a 71 75 6f 28 71 |p)) {.X.|..zquo(q|
|00001f30| 32 2d 3e 6e 75 6d 2c 20 | 74 6d 70 2c 20 26 6e 32 |2->num, |tmp, &n2|
|00001f40| 29 3b 0a 58 09 09 09 7a | 71 75 6f 28 71 31 2d 3e |);.X...z|quo(q1->|
|00001f50| 64 65 6e 2c 20 74 6d 70 | 2c 20 26 64 31 29 3b 0a |den, tmp|, &d1);.|
|00001f60| 58 09 09 7d 0a 58 09 09 | 66 72 65 65 68 28 74 6d |X..}.X..|freeh(tm|
|00001f70| 70 2e 76 29 3b 0a 58 09 | 7d 0a 58 09 72 20 3d 20 |p.v);.X.|}.X.r = |
|00001f80| 71 61 6c 6c 6f 63 28 29 | 3b 0a 58 09 7a 6d 75 6c |qalloc()|;.X.zmul|
|00001f90| 28 6e 31 2c 20 6e 32 2c | 20 26 72 2d 3e 6e 75 6d |(n1, n2,| &r->num|
|00001fa0| 29 3b 0a 58 09 7a 6d 75 | 6c 28 64 31 2c 20 64 32 |);.X.zmu|l(d1, d2|
|00001fb0| 2c 20 26 72 2d 3e 64 65 | 6e 29 3b 0a 58 09 69 66 |, &r->de|n);.X.if|
|00001fc0| 20 28 71 31 2d 3e 6e 75 | 6d 2e 76 20 21 3d 20 6e | (q1->nu|m.v != n|
|00001fd0| 31 2e 76 29 0a 58 09 09 | 66 72 65 65 68 28 6e 31 |1.v).X..|freeh(n1|
|00001fe0| 2e 76 29 3b 0a 58 09 69 | 66 20 28 71 31 2d 3e 64 |.v);.X.i|f (q1->d|
|00001ff0| 65 6e 2e 76 20 21 3d 20 | 64 31 2e 76 29 0a 58 09 |en.v != |d1.v).X.|
|00002000| 09 66 72 65 65 68 28 64 | 31 2e 76 29 3b 0a 58 09 |.freeh(d|1.v);.X.|
|00002010| 69 66 20 28 71 32 2d 3e | 6e 75 6d 2e 76 20 21 3d |if (q2->|num.v !=|
|00002020| 20 6e 32 2e 76 29 0a 58 | 09 09 66 72 65 65 68 28 | n2.v).X|..freeh(|
|00002030| 6e 32 2e 76 29 3b 0a 58 | 09 69 66 20 28 71 32 2d |n2.v);.X|.if (q2-|
|00002040| 3e 64 65 6e 2e 76 20 21 | 3d 20 64 32 2e 76 29 0a |>den.v !|= d2.v).|
|00002050| 58 09 09 66 72 65 65 68 | 28 64 32 2e 76 29 3b 0a |X..freeh|(d2.v);.|
|00002060| 58 09 72 65 74 75 72 6e | 20 72 3b 0a 58 7d 0a 58 |X.return| r;.X}.X|
|00002070| 0a 58 0a 58 23 69 66 20 | 30 0a 58 2f 2a 0a 58 20 |.X.X#if |0.X/*.X |
|00002080| 2a 20 4d 75 6c 74 69 70 | 6c 79 20 61 20 6e 75 6d |* Multip|ly a num|
|00002090| 62 65 72 20 62 79 20 61 | 20 73 6d 61 6c 6c 20 69 |ber by a| small i|
|000020a0| 6e 74 65 67 65 72 2e 0a | 58 20 2a 09 71 32 20 3d |nteger..|X *.q2 =|
|000020b0| 20 71 6d 75 6c 69 28 71 | 31 2c 20 6e 29 3b 0a 58 | qmuli(q|1, n);.X|
|000020c0| 20 2a 2f 0a 58 4e 55 4d | 42 45 52 20 2a 0a 58 71 | */.XNUM|BER *.Xq|
|000020d0| 6d 75 6c 69 28 71 2c 20 | 6e 29 0a 58 09 4e 55 4d |muli(q, |n).X.NUM|
|000020e0| 42 45 52 20 2a 71 3b 0a | 58 09 6c 6f 6e 67 20 6e |BER *q;.|X.long n|
|000020f0| 3b 0a 58 7b 0a 58 09 4e | 55 4d 42 45 52 20 2a 72 |;.X{.X.N|UMBER *r|
|00002100| 3b 0a 58 09 6c 6f 6e 67 | 20 64 3b 09 09 09 2f 2a |;.X.long| d;.../*|
|00002110| 20 67 63 64 20 6f 66 20 | 6d 75 6c 74 69 70 6c 69 | gcd of |multipli|
|00002120| 65 72 20 61 6e 64 20 64 | 65 6e 6f 6d 69 6e 61 74 |er and d|enominat|
|00002130| 6f 72 20 2a 2f 0a 58 09 | 69 6e 74 20 73 69 67 6e |or */.X.|int sign|
|00002140| 3b 0a 58 0a 58 09 69 66 | 20 28 28 6e 20 3d 3d 20 |;.X.X.if| ((n == |
|00002150| 30 29 20 7c 7c 20 71 69 | 73 7a 65 72 6f 28 71 29 |0) || qi|szero(q)|
|00002160| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 71 6c 69 6e |).X..ret|urn qlin|
|00002170| 6b 28 26 5f 71 7a 65 72 | 6f 5f 29 3b 0a 58 09 69 |k(&_qzer|o_);.X.i|
|00002180| 66 20 28 6e 20 3d 3d 20 | 31 29 0a 58 09 09 72 65 |f (n == |1).X..re|
|00002190| 74 75 72 6e 20 71 6c 69 | 6e 6b 28 71 29 3b 0a 58 |turn qli|nk(q);.X|
|000021a0| 09 72 20 3d 20 71 61 6c | 6c 6f 63 28 29 3b 0a 58 |.r = qal|loc();.X|
|000021b0| 09 69 66 20 28 71 69 73 | 69 6e 74 28 71 29 29 20 |.if (qis|int(q)) |
|000021c0| 7b 0a 58 09 09 7a 6d 75 | 6c 69 28 71 2d 3e 6e 75 |{.X..zmu|li(q->nu|
|000021d0| 6d 2c 20 6e 2c 20 26 72 | 2d 3e 6e 75 6d 29 3b 0a |m, n, &r|->num);.|
|000021e0| 58 09 09 72 65 74 75 72 | 6e 20 72 3b 0a 58 09 7d |X..retur|n r;.X.}|
|000021f0| 0a 58 09 73 69 67 6e 20 | 3d 20 31 3b 0a 58 09 69 |.X.sign |= 1;.X.i|
|00002200| 66 20 28 6e 20 3c 20 30 | 29 20 7b 0a 58 09 09 6e |f (n < 0|) {.X..n|
|00002210| 20 3d 20 2d 6e 3b 0a 58 | 09 09 73 69 67 6e 20 3d | = -n;.X|..sign =|
|00002220| 20 2d 31 3b 0a 58 09 7d | 0a 58 09 64 20 3d 20 7a | -1;.X.}|.X.d = z|
|00002230| 6d 6f 64 69 28 71 2d 3e | 64 65 6e 2c 20 6e 29 3b |modi(q->|den, n);|
|00002240| 0a 58 09 64 20 3d 20 69 | 69 67 63 64 28 64 2c 20 |.X.d = i|igcd(d, |
|00002250| 6e 29 3b 0a 58 09 7a 6d | 75 6c 69 28 71 2d 3e 6e |n);.X.zm|uli(q->n|
|00002260| 75 6d 2c 20 28 6e 20 2a | 20 73 69 67 6e 29 20 2f |um, (n *| sign) /|
|00002270| 20 64 2c 20 26 72 2d 3e | 6e 75 6d 29 3b 0a 58 09 | d, &r->|num);.X.|
|00002280| 28 76 6f 69 64 29 20 7a | 64 69 76 69 28 71 2d 3e |(void) z|divi(q->|
|00002290| 64 65 6e 2c 20 64 2c 20 | 26 72 2d 3e 64 65 6e 29 |den, d, |&r->den)|
|000022a0| 3b 0a 58 09 72 65 74 75 | 72 6e 20 72 3b 0a 58 7d |;.X.retu|rn r;.X}|
|000022b0| 0a 58 23 65 6e 64 69 66 | 0a 58 0a 58 0a 58 2f 2a |.X#endif|.X.X.X/*|
|000022c0| 0a 58 20 2a 20 44 69 76 | 69 64 65 20 74 77 6f 20 |.X * Div|ide two |
|000022d0| 6e 75 6d 62 65 72 73 20 | 28 61 73 20 66 72 61 63 |numbers |(as frac|
|000022e0| 74 69 6f 6e 73 29 2e 0a | 58 20 2a 09 71 33 20 3d |tions)..|X *.q3 =|
|000022f0| 20 71 64 69 76 28 71 31 | 2c 20 71 32 29 3b 0a 58 | qdiv(q1|, q2);.X|
|00002300| 20 2a 2f 0a 58 4e 55 4d | 42 45 52 20 2a 0a 58 71 | */.XNUM|BER *.Xq|
|00002310| 64 69 76 28 71 31 2c 20 | 71 32 29 0a 58 09 72 65 |div(q1, |q2).X.re|
|00002320| 67 69 73 74 65 72 20 4e | 55 4d 42 45 52 20 2a 71 |gister N|UMBER *q|
|00002330| 31 2c 20 2a 71 32 3b 0a | 58 7b 0a 58 09 4e 55 4d |1, *q2;.|X{.X.NUM|
|00002340| 42 45 52 20 74 65 6d 70 | 3b 0a 58 0a 58 09 69 66 |BER temp|;.X.X.if|
|00002350| 20 28 71 31 20 3d 3d 20 | 71 32 29 20 7b 0a 58 09 | (q1 == |q2) {.X.|
|00002360| 09 69 66 20 28 71 69 73 | 7a 65 72 6f 28 71 32 29 |.if (qis|zero(q2)|
|00002370| 29 0a 58 09 09 09 65 72 | 72 6f 72 28 22 44 69 76 |).X...er|ror("Div|
|00002380| 69 73 69 6f 6e 20 62 79 | 20 7a 65 72 6f 22 29 3b |ision by| zero");|
|00002390| 0a 58 09 09 72 65 74 75 | 72 6e 20 71 6c 69 6e 6b |.X..retu|rn qlink|
|000023a0| 28 26 5f 71 6f 6e 65 5f | 29 3b 0a 58 09 7d 0a 58 |(&_qone_|);.X.}.X|
|000023b0| 09 69 66 20 28 71 69 73 | 6f 6e 65 28 71 31 29 29 |.if (qis|one(q1))|
|000023c0| 0a 58 09 09 72 65 74 75 | 72 6e 20 71 69 6e 76 28 |.X..retu|rn qinv(|
|000023d0| 71 32 29 3b 0a 58 09 74 | 65 6d 70 2e 6e 75 6d 20 |q2);.X.t|emp.num |
|000023e0| 3d 20 71 32 2d 3e 64 65 | 6e 3b 0a 58 09 74 65 6d |= q2->de|n;.X.tem|
|000023f0| 70 2e 64 65 6e 20 3d 20 | 71 32 2d 3e 6e 75 6d 3b |p.den = |q2->num;|
|00002400| 0a 58 09 74 65 6d 70 2e | 6e 75 6d 2e 73 69 67 6e |.X.temp.|num.sign|
|00002410| 20 3d 20 74 65 6d 70 2e | 64 65 6e 2e 73 69 67 6e | = temp.|den.sign|
|00002420| 3b 0a 58 09 74 65 6d 70 | 2e 64 65 6e 2e 73 69 67 |;.X.temp|.den.sig|
|00002430| 6e 20 3d 20 30 3b 0a 58 | 09 74 65 6d 70 2e 6c 69 |n = 0;.X|.temp.li|
|00002440| 6e 6b 73 20 3d 20 31 3b | 0a 58 09 72 65 74 75 72 |nks = 1;|.X.retur|
|00002450| 6e 20 71 6d 75 6c 28 71 | 31 2c 20 26 74 65 6d 70 |n qmul(q|1, &temp|
|00002460| 29 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |);.X}.X.|X.X/*.X |
|00002470| 2a 20 44 69 76 69 64 65 | 20 61 20 6e 75 6d 62 65 |* Divide| a numbe|
|00002480| 72 20 62 79 20 61 20 73 | 6d 61 6c 6c 20 69 6e 74 |r by a s|mall int|
|00002490| 65 67 65 72 2e 0a 58 20 | 2a 09 71 32 20 3d 20 71 |eger..X |*.q2 = q|
|000024a0| 64 69 76 69 28 71 31 2c | 20 6e 29 3b 0a 58 20 2a |divi(q1,| n);.X *|
|000024b0| 2f 0a 58 4e 55 4d 42 45 | 52 20 2a 0a 58 71 64 69 |/.XNUMBE|R *.Xqdi|
|000024c0| 76 69 28 71 2c 20 6e 29 | 0a 58 09 4e 55 4d 42 45 |vi(q, n)|.X.NUMBE|
|000024d0| 52 20 2a 71 3b 0a 58 09 | 6c 6f 6e 67 20 6e 3b 0a |R *q;.X.|long n;.|
|000024e0| 58 7b 0a 58 09 4e 55 4d | 42 45 52 20 2a 72 3b 0a |X{.X.NUM|BER *r;.|
|000024f0| 58 09 6c 6f 6e 67 20 64 | 3b 09 09 09 2f 2a 20 67 |X.long d|;.../* g|
|00002500| 63 64 20 6f 66 20 64 69 | 76 69 73 6f 72 20 61 6e |cd of di|visor an|
|00002510| 64 20 6e 75 6d 65 72 61 | 74 6f 72 20 2a 2f 0a 58 |d numera|tor */.X|
|00002520| 09 69 6e 74 20 73 69 67 | 6e 3b 0a 58 0a 58 09 69 |.int sig|n;.X.X.i|
|00002530| 66 20 28 6e 20 3d 3d 20 | 30 29 0a 58 09 09 65 72 |f (n == |0).X..er|
|00002540| 72 6f 72 28 22 44 69 76 | 69 73 69 6f 6e 20 62 79 |ror("Div|ision by|
|00002550| 20 7a 65 72 6f 22 29 3b | 0a 58 09 69 66 20 28 28 | zero");|.X.if ((|
|00002560| 6e 20 3d 3d 20 31 29 20 | 7c 7c 20 71 69 73 7a 65 |n == 1) ||| qisze|
|00002570| 72 6f 28 71 29 29 0a 58 | 09 09 72 65 74 75 72 6e |ro(q)).X|..return|
|00002580| 20 71 6c 69 6e 6b 28 71 | 29 3b 0a 58 09 73 69 67 | qlink(q|);.X.sig|
|00002590| 6e 20 3d 20 31 3b 0a 58 | 09 69 66 20 28 6e 20 3c |n = 1;.X|.if (n <|
|000025a0| 20 30 29 20 7b 0a 58 09 | 09 6e 20 3d 20 2d 6e 3b | 0) {.X.|.n = -n;|
|000025b0| 0a 58 09 09 73 69 67 6e | 20 3d 20 2d 31 3b 0a 58 |.X..sign| = -1;.X|
|000025c0| 09 7d 0a 58 09 72 20 3d | 20 71 61 6c 6c 6f 63 28 |.}.X.r =| qalloc(|
|000025d0| 29 3b 0a 58 09 64 20 3d | 20 7a 6d 6f 64 69 28 71 |);.X.d =| zmodi(q|
|000025e0| 2d 3e 6e 75 6d 2c 20 6e | 29 3b 0a 58 09 64 20 3d |->num, n|);.X.d =|
|000025f0| 20 69 69 67 63 64 28 64 | 2c 20 6e 29 3b 0a 58 09 | iigcd(d|, n);.X.|
|00002600| 28 76 6f 69 64 29 20 7a | 64 69 76 69 28 71 2d 3e |(void) z|divi(q->|
|00002610| 6e 75 6d 2c 20 64 20 2a | 20 73 69 67 6e 2c 20 26 |num, d *| sign, &|
|00002620| 72 2d 3e 6e 75 6d 29 3b | 0a 58 09 7a 6d 75 6c 69 |r->num);|.X.zmuli|
|00002630| 28 71 2d 3e 64 65 6e 2c | 20 6e 20 2f 20 64 2c 20 |(q->den,| n / d, |
|00002640| 26 72 2d 3e 64 65 6e 29 | 3b 0a 58 09 72 65 74 75 |&r->den)|;.X.retu|
|00002650| 72 6e 20 72 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a |rn r;.X}|.X.X.X/*|
|00002660| 0a 58 20 2a 20 52 65 74 | 75 72 6e 20 74 68 65 20 |.X * Ret|urn the |
|00002670| 71 75 6f 74 69 65 6e 74 | 20 77 68 65 6e 20 6f 6e |quotient| when on|
|00002680| 65 20 6e 75 6d 62 65 72 | 20 69 73 20 64 69 76 69 |e number| is divi|
|00002690| 64 65 64 20 62 79 20 61 | 6e 6f 74 68 65 72 2e 0a |ded by a|nother..|
|000026a0| 58 20 2a 20 54 68 69 73 | 20 77 6f 72 6b 73 20 66 |X * This| works f|
|000026b0| 6f 72 20 66 72 61 63 74 | 69 6f 6e 61 6c 20 76 61 |or fract|ional va|
|000026c0| 6c 75 65 73 20 61 6c 73 | 6f 2c 20 61 6e 64 20 69 |lues als|o, and i|
|000026d0| 6e 20 61 6c 6c 20 63 61 | 73 65 73 3a 0a 58 20 2a |n all ca|ses:.X *|
|000026e0| 09 71 71 75 6f 28 71 31 | 2c 20 71 32 29 20 3d 20 |.qquo(q1|, q2) = |
|000026f0| 69 6e 74 28 71 31 20 2f | 20 71 32 29 2e 0a 58 20 |int(q1 /| q2)..X |
|00002700| 2a 2f 0a 58 4e 55 4d 42 | 45 52 20 2a 0a 58 71 71 |*/.XNUMB|ER *.Xqq|
|00002710| 75 6f 28 71 31 2c 20 71 | 32 29 0a 58 09 72 65 67 |uo(q1, q|2).X.reg|
|00002720| 69 73 74 65 72 20 4e 55 | 4d 42 45 52 20 2a 71 31 |ister NU|MBER *q1|
|00002730| 2c 20 2a 71 32 3b 0a 58 | 7b 0a 58 09 5a 56 41 4c |, *q2;.X|{.X.ZVAL|
|00002740| 55 45 20 6e 75 6d 2c 20 | 64 65 6e 2c 20 72 65 73 |UE num, |den, res|
|00002750| 3b 0a 58 09 4e 55 4d 42 | 45 52 20 2a 71 3b 0a 58 |;.X.NUMB|ER *q;.X|
|00002760| 0a 58 09 69 66 20 28 69 | 73 75 6e 69 74 28 71 31 |.X.if (i|sunit(q1|
|00002770| 2d 3e 6e 75 6d 29 29 0a | 58 09 09 6e 75 6d 20 3d |->num)).|X..num =|
|00002780| 20 71 32 2d 3e 64 65 6e | 3b 0a 58 09 65 6c 73 65 | q2->den|;.X.else|
|00002790| 20 69 66 20 28 69 73 75 | 6e 69 74 28 71 32 2d 3e | if (isu|nit(q2->|
|000027a0| 64 65 6e 29 29 0a 58 09 | 09 6e 75 6d 20 3d 20 71 |den)).X.|.num = q|
|000027b0| 31 2d 3e 6e 75 6d 3b 0a | 58 09 65 6c 73 65 0a 58 |1->num;.|X.else.X|
|000027c0| 09 09 7a 6d 75 6c 28 71 | 31 2d 3e 6e 75 6d 2c 20 |..zmul(q|1->num, |
|000027d0| 71 32 2d 3e 64 65 6e 2c | 20 26 6e 75 6d 29 3b 0a |q2->den,| &num);.|
|000027e0| 58 09 69 66 20 28 69 73 | 75 6e 69 74 28 71 31 2d |X.if (is|unit(q1-|
|000027f0| 3e 64 65 6e 29 29 0a 58 | 09 09 64 65 6e 20 3d 20 |>den)).X|..den = |
|00002800| 71 32 2d 3e 6e 75 6d 3b | 0a 58 09 65 6c 73 65 20 |q2->num;|.X.else |
|00002810| 69 66 20 28 69 73 75 6e | 69 74 28 71 32 2d 3e 6e |if (isun|it(q2->n|
|00002820| 75 6d 29 29 0a 58 09 09 | 64 65 6e 20 3d 20 71 31 |um)).X..|den = q1|
|00002830| 2d 3e 64 65 6e 3b 0a 58 | 09 65 6c 73 65 0a 58 09 |->den;.X|.else.X.|
|00002840| 09 7a 6d 75 6c 28 71 31 | 2d 3e 64 65 6e 2c 20 71 |.zmul(q1|->den, q|
|00002850| 32 2d 3e 6e 75 6d 2c 20 | 26 64 65 6e 29 3b 0a 58 |2->num, |&den);.X|
|00002860| 09 7a 71 75 6f 28 6e 75 | 6d 2c 20 64 65 6e 2c 20 |.zquo(nu|m, den, |
|00002870| 26 72 65 73 29 3b 0a 58 | 09 69 66 20 28 28 6e 75 |&res);.X|.if ((nu|
|00002880| 6d 2e 76 20 21 3d 20 71 | 32 2d 3e 64 65 6e 2e 76 |m.v != q|2->den.v|
|00002890| 29 20 26 26 20 28 6e 75 | 6d 2e 76 20 21 3d 20 71 |) && (nu|m.v != q|
|000028a0| 31 2d 3e 6e 75 6d 2e 76 | 29 29 0a 58 09 09 66 72 |1->num.v|)).X..fr|
|000028b0| 65 65 68 28 6e 75 6d 2e | 76 29 3b 0a 58 09 69 66 |eeh(num.|v);.X.if|
|000028c0| 20 28 28 64 65 6e 2e 76 | 20 21 3d 20 71 32 2d 3e | ((den.v| != q2->|
|000028d0| 6e 75 6d 2e 76 29 20 26 | 26 20 28 64 65 6e 2e 76 |num.v) &|& (den.v|
|000028e0| 20 21 3d 20 71 31 2d 3e | 64 65 6e 2e 76 29 29 0a | != q1->|den.v)).|
|000028f0| 58 09 09 66 72 65 65 68 | 28 64 65 6e 2e 76 29 3b |X..freeh|(den.v);|
|00002900| 0a 58 09 69 66 20 28 69 | 73 7a 65 72 6f 28 72 65 |.X.if (i|szero(re|
|00002910| 73 29 29 20 7b 0a 58 09 | 09 66 72 65 65 68 28 72 |s)) {.X.|.freeh(r|
|00002920| 65 73 2e 76 29 3b 0a 58 | 09 09 72 65 74 75 72 6e |es.v);.X|..return|
|00002930| 20 71 6c 69 6e 6b 28 26 | 5f 71 7a 65 72 6f 5f 29 | qlink(&|_qzero_)|
|00002940| 3b 0a 58 09 7d 0a 58 09 | 69 66 20 28 69 73 75 6e |;.X.}.X.|if (isun|
|00002950| 69 74 28 72 65 73 29 29 | 20 7b 0a 58 09 09 71 20 |it(res))| {.X..q |
|00002960| 3d 20 28 72 65 73 2e 73 | 69 67 6e 20 3f 20 26 5f |= (res.s|ign ? &_|
|00002970| 71 6e 65 67 6f 6e 65 5f | 20 3a 20 26 5f 71 6f 6e |qnegone_| : &_qon|
|00002980| 65 5f 29 3b 0a 58 09 09 | 66 72 65 65 68 28 72 65 |e_);.X..|freeh(re|
|00002990| 73 2e 76 29 3b 0a 58 09 | 09 72 65 74 75 72 6e 20 |s.v);.X.|.return |
|000029a0| 71 6c 69 6e 6b 28 71 29 | 3b 0a 58 09 7d 0a 58 09 |qlink(q)|;.X.}.X.|
|000029b0| 71 20 3d 20 71 61 6c 6c | 6f 63 28 29 3b 0a 58 09 |q = qall|oc();.X.|
|000029c0| 71 2d 3e 6e 75 6d 20 3d | 20 72 65 73 3b 0a 58 09 |q->num =| res;.X.|
|000029d0| 72 65 74 75 72 6e 20 71 | 3b 0a 58 7d 0a 58 0a 58 |return q|;.X}.X.X|
|000029e0| 0a 58 2f 2a 0a 58 20 2a | 20 52 65 74 75 72 6e 20 |.X/*.X *| Return |
|000029f0| 74 68 65 20 61 62 73 6f | 6c 75 74 65 20 76 61 6c |the abso|lute val|
|00002a00| 75 65 20 6f 66 20 61 20 | 6e 75 6d 62 65 72 2e 0a |ue of a |number..|
|00002a10| 58 20 2a 09 71 32 20 3d | 20 71 61 62 73 28 71 31 |X *.q2 =| qabs(q1|
|00002a20| 29 3b 0a 58 20 2a 2f 0a | 58 4e 55 4d 42 45 52 20 |);.X */.|XNUMBER |
|00002a30| 2a 0a 58 71 61 62 73 28 | 71 29 0a 58 09 72 65 67 |*.Xqabs(|q).X.reg|
|00002a40| 69 73 74 65 72 20 4e 55 | 4d 42 45 52 20 2a 71 3b |ister NU|MBER *q;|
|00002a50| 0a 58 7b 0a 58 09 72 65 | 67 69 73 74 65 72 20 4e |.X{.X.re|gister N|
|00002a60| 55 4d 42 45 52 20 2a 72 | 3b 0a 58 0a 58 09 69 66 |UMBER *r|;.X.X.if|
|00002a70| 20 28 71 2d 3e 6e 75 6d | 2e 73 69 67 6e 20 3d 3d | (q->num|.sign ==|
|00002a80| 20 30 29 0a 58 09 09 72 | 65 74 75 72 6e 20 71 6c | 0).X..r|eturn ql|
|00002a90| 69 6e 6b 28 71 29 3b 0a | 58 09 72 20 3d 20 71 61 |ink(q);.|X.r = qa|
|00002aa0| 6c 6c 6f 63 28 29 3b 0a | 58 09 69 66 20 28 21 69 |lloc();.|X.if (!i|
|00002ab0| 73 75 6e 69 74 28 71 2d | 3e 6e 75 6d 29 29 0a 58 |sunit(q-|>num)).X|
|00002ac0| 09 09 7a 63 6f 70 79 28 | 71 2d 3e 6e 75 6d 2c 20 |..zcopy(|q->num, |
|00002ad0| 26 72 2d 3e 6e 75 6d 29 | 3b 0a 58 09 69 66 20 28 |&r->num)|;.X.if (|
|00002ae0| 21 69 73 75 6e 69 74 28 | 71 2d 3e 64 65 6e 29 29 |!isunit(|q->den))|
|00002af0| 0a 58 09 09 7a 63 6f 70 | 79 28 71 2d 3e 64 65 6e |.X..zcop|y(q->den|
|00002b00| 2c 20 26 72 2d 3e 64 65 | 6e 29 3b 0a 58 09 72 2d |, &r->de|n);.X.r-|
|00002b10| 3e 6e 75 6d 2e 73 69 67 | 6e 20 3d 20 30 3b 0a 58 |>num.sig|n = 0;.X|
|00002b20| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 7d 0a 58 0a |.return |r;.X}.X.|
|00002b30| 58 0a 58 2f 2a 0a 58 20 | 2a 20 4e 65 67 61 74 65 |X.X/*.X |* Negate|
|00002b40| 20 61 20 6e 75 6d 62 65 | 72 2e 0a 58 20 2a 09 71 | a numbe|r..X *.q|
|00002b50| 32 20 3d 20 71 6e 65 67 | 28 71 31 29 3b 0a 58 20 |2 = qneg|(q1);.X |
|00002b60| 2a 2f 0a 58 4e 55 4d 42 | 45 52 20 2a 0a 58 71 6e |*/.XNUMB|ER *.Xqn|
|00002b70| 65 67 28 71 29 0a 58 09 | 72 65 67 69 73 74 65 72 |eg(q).X.|register|
|00002b80| 20 4e 55 4d 42 45 52 20 | 2a 71 3b 0a 58 7b 0a 58 | NUMBER |*q;.X{.X|
|00002b90| 09 72 65 67 69 73 74 65 | 72 20 4e 55 4d 42 45 52 |.registe|r NUMBER|
|00002ba0| 20 2a 72 3b 0a 58 0a 58 | 09 69 66 20 28 71 69 73 | *r;.X.X|.if (qis|
|00002bb0| 7a 65 72 6f 28 71 29 29 | 0a 58 09 09 72 65 74 75 |zero(q))|.X..retu|
|00002bc0| 72 6e 20 71 6c 69 6e 6b | 28 26 5f 71 7a 65 72 6f |rn qlink|(&_qzero|
|00002bd0| 5f 29 3b 0a 58 09 72 20 | 3d 20 71 61 6c 6c 6f 63 |_);.X.r |= qalloc|
|00002be0| 28 29 3b 0a 58 09 69 66 | 20 28 21 69 73 75 6e 69 |();.X.if| (!isuni|
|00002bf0| 74 28 71 2d 3e 6e 75 6d | 29 29 0a 58 09 09 7a 63 |t(q->num|)).X..zc|
|00002c00| 6f 70 79 28 71 2d 3e 6e | 75 6d 2c 20 26 72 2d 3e |opy(q->n|um, &r->|
|00002c10| 6e 75 6d 29 3b 0a 58 09 | 69 66 20 28 21 69 73 75 |num);.X.|if (!isu|
|00002c20| 6e 69 74 28 71 2d 3e 64 | 65 6e 29 29 0a 58 09 09 |nit(q->d|en)).X..|
|00002c30| 7a 63 6f 70 79 28 71 2d | 3e 64 65 6e 2c 20 26 72 |zcopy(q-|>den, &r|
|00002c40| 2d 3e 64 65 6e 29 3b 0a | 58 09 72 2d 3e 6e 75 6d |->den);.|X.r->num|
|00002c50| 2e 73 69 67 6e 20 3d 20 | 21 71 2d 3e 6e 75 6d 2e |.sign = |!q->num.|
|00002c60| 73 69 67 6e 3b 0a 58 09 | 72 65 74 75 72 6e 20 72 |sign;.X.|return r|
|00002c70| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00002c80| 20 52 65 74 75 72 6e 20 | 74 68 65 20 73 69 67 6e | Return |the sign|
|00002c90| 20 6f 66 20 61 20 6e 75 | 6d 62 65 72 20 28 2d 31 | of a nu|mber (-1|
|00002ca0| 2c 20 30 20 6f 72 20 31 | 29 0a 58 20 2a 2f 0a 58 |, 0 or 1|).X */.X|
|00002cb0| 4e 55 4d 42 45 52 20 2a | 0a 58 71 73 69 67 6e 28 |NUMBER *|.Xqsign(|
|00002cc0| 71 29 0a 58 09 4e 55 4d | 42 45 52 20 2a 71 3b 0a |q).X.NUM|BER *q;.|
|00002cd0| 58 7b 0a 58 09 69 66 20 | 28 71 69 73 7a 65 72 6f |X{.X.if |(qiszero|
|00002ce0| 28 71 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 71 |(q)).X..|return q|
|00002cf0| 6c 69 6e 6b 28 26 5f 71 | 7a 65 72 6f 5f 29 3b 0a |link(&_q|zero_);.|
|00002d00| 58 09 69 66 20 28 71 69 | 73 6e 65 67 28 71 29 29 |X.if (qi|sneg(q))|
|00002d10| 0a 58 09 09 72 65 74 75 | 72 6e 20 71 6c 69 6e 6b |.X..retu|rn qlink|
|00002d20| 28 26 5f 71 6e 65 67 6f | 6e 65 5f 29 3b 0a 58 09 |(&_qnego|ne_);.X.|
|00002d30| 72 65 74 75 72 6e 20 71 | 6c 69 6e 6b 28 26 5f 71 |return q|link(&_q|
|00002d40| 6f 6e 65 5f 29 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |one_);.X|}.X.X.X/|
|00002d50| 2a 0a 58 20 2a 20 49 6e | 76 65 72 74 20 61 20 6e |*.X * In|vert a n|
|00002d60| 75 6d 62 65 72 2e 0a 58 | 20 2a 09 71 32 20 3d 20 |umber..X| *.q2 = |
|00002d70| 71 69 6e 76 28 71 31 29 | 3b 0a 58 20 2a 2f 0a 58 |qinv(q1)|;.X */.X|
|00002d80| 4e 55 4d 42 45 52 20 2a | 0a 58 71 69 6e 76 28 71 |NUMBER *|.Xqinv(q|
|00002d90| 29 0a 58 09 72 65 67 69 | 73 74 65 72 20 4e 55 4d |).X.regi|ster NUM|
|00002da0| 42 45 52 20 2a 71 3b 0a | 58 7b 0a 58 09 72 65 67 |BER *q;.|X{.X.reg|
|00002db0| 69 73 74 65 72 20 4e 55 | 4d 42 45 52 20 2a 72 3b |ister NU|MBER *r;|
|00002dc0| 0a 58 0a 58 09 69 66 20 | 28 71 69 73 75 6e 69 74 |.X.X.if |(qisunit|
|00002dd0| 28 71 29 29 20 7b 0a 58 | 09 09 72 20 3d 20 28 71 |(q)) {.X|..r = (q|
|00002de0| 69 73 6e 65 67 28 71 29 | 20 3f 20 26 5f 71 6e 65 |isneg(q)| ? &_qne|
|00002df0| 67 6f 6e 65 5f 20 3a 20 | 26 5f 71 6f 6e 65 5f 29 |gone_ : |&_qone_)|
|00002e00| 3b 0a 58 09 09 72 65 74 | 75 72 6e 20 71 6c 69 6e |;.X..ret|urn qlin|
|00002e10| 6b 28 72 29 3b 0a 58 09 | 7d 0a 58 09 69 66 20 28 |k(r);.X.|}.X.if (|
|00002e20| 71 69 73 7a 65 72 6f 28 | 71 29 29 0a 58 09 09 65 |qiszero(|q)).X..e|
|00002e30| 72 72 6f 72 28 22 44 69 | 76 69 73 69 6f 6e 20 62 |rror("Di|vision b|
|00002e40| 79 20 7a 65 72 6f 22 29 | 3b 0a 58 09 72 20 3d 20 |y zero")|;.X.r = |
|00002e50| 71 61 6c 6c 6f 63 28 29 | 3b 0a 58 09 69 66 20 28 |qalloc()|;.X.if (|
|00002e60| 21 69 73 75 6e 69 74 28 | 71 2d 3e 6e 75 6d 29 29 |!isunit(|q->num))|
|00002e70| 0a 58 09 09 7a 63 6f 70 | 79 28 71 2d 3e 6e 75 6d |.X..zcop|y(q->num|
|00002e80| 2c 20 26 72 2d 3e 64 65 | 6e 29 3b 0a 58 09 69 66 |, &r->de|n);.X.if|
|00002e90| 20 28 21 69 73 75 6e 69 | 74 28 71 2d 3e 64 65 6e | (!isuni|t(q->den|
|00002ea0| 29 29 0a 58 09 09 7a 63 | 6f 70 79 28 71 2d 3e 64 |)).X..zc|opy(q->d|
|00002eb0| 65 6e 2c 20 26 72 2d 3e | 6e 75 6d 29 3b 0a 58 09 |en, &r->|num);.X.|
|00002ec0| 72 2d 3e 6e 75 6d 2e 73 | 69 67 6e 20 3d 20 71 2d |r->num.s|ign = q-|
|00002ed0| 3e 6e 75 6d 2e 73 69 67 | 6e 3b 0a 58 09 72 2d 3e |>num.sig|n;.X.r->|
|00002ee0| 64 65 6e 2e 73 69 67 6e | 20 3d 20 30 3b 0a 58 09 |den.sign| = 0;.X.|
|00002ef0| 72 65 74 75 72 6e 20 72 | 3b 0a 58 7d 0a 58 0a 58 |return r|;.X}.X.X|
|00002f00| 0a 58 2f 2a 0a 58 20 2a | 20 52 65 74 75 72 6e 20 |.X/*.X *| Return |
|00002f10| 6a 75 73 74 20 74 68 65 | 20 6e 75 6d 65 72 61 74 |just the| numerat|
|00002f20| 6f 72 20 6f 66 20 61 20 | 6e 75 6d 62 65 72 2e 0a |or of a |number..|
|00002f30| 58 20 2a 09 71 32 20 3d | 20 71 6e 75 6d 28 71 31 |X *.q2 =| qnum(q1|
|00002f40| 29 3b 0a 58 20 2a 2f 0a | 58 4e 55 4d 42 45 52 20 |);.X */.|XNUMBER |
|00002f50| 2a 0a 58 71 6e 75 6d 28 | 71 29 0a 58 09 72 65 67 |*.Xqnum(|q).X.reg|
|00002f60| 69 73 74 65 72 20 4e 55 | 4d 42 45 52 20 2a 71 3b |ister NU|MBER *q;|
|00002f70| 0a 58 7b 0a 58 09 72 65 | 67 69 73 74 65 72 20 4e |.X{.X.re|gister N|
|00002f80| 55 4d 42 45 52 20 2a 72 | 3b 0a 58 0a 58 09 69 66 |UMBER *r|;.X.X.if|
|00002f90| 20 28 71 69 73 69 6e 74 | 28 71 29 29 0a 58 09 09 | (qisint|(q)).X..|
|00002fa0| 72 65 74 75 72 6e 20 71 | 6c 69 6e 6b 28 71 29 3b |return q|link(q);|
|00002fb0| 0a 58 09 69 66 20 28 69 | 73 75 6e 69 74 28 71 2d |.X.if (i|sunit(q-|
|00002fc0| 3e 6e 75 6d 29 29 20 7b | 0a 58 09 09 72 20 3d 20 |>num)) {|.X..r = |
|00002fd0| 28 71 69 73 6e 65 67 28 | 71 29 20 3f 20 26 5f 71 |(qisneg(|q) ? &_q|
|00002fe0| 6e 65 67 6f 6e 65 5f 20 | 3a 20 26 5f 71 6f 6e 65 |negone_ |: &_qone|
|00002ff0| 5f 29 3b 0a 58 09 09 72 | 65 74 75 72 6e 20 71 6c |_);.X..r|eturn ql|
|00003000| 69 6e 6b 28 72 29 3b 0a | 58 09 7d 0a 58 09 72 20 |ink(r);.|X.}.X.r |
|00003010| 3d 20 71 61 6c 6c 6f 63 | 28 29 3b 0a 58 09 7a 63 |= qalloc|();.X.zc|
|00003020| 6f 70 79 28 71 2d 3e 6e | 75 6d 2c 20 26 72 2d 3e |opy(q->n|um, &r->|
|00003030| 6e 75 6d 29 3b 0a 58 09 | 72 65 74 75 72 6e 20 72 |num);.X.|return r|
|00003040| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00003050| 20 52 65 74 75 72 6e 20 | 6a 75 73 74 20 74 68 65 | Return |just the|
|00003060| 20 64 65 6e 6f 6d 69 6e | 61 74 6f 72 20 6f 66 20 | denomin|ator of |
|00003070| 61 20 6e 75 6d 62 65 72 | 2e 0a 58 20 2a 09 71 32 |a number|..X *.q2|
|00003080| 20 3d 20 71 64 65 6e 28 | 71 31 29 3b 0a 58 20 2a | = qden(|q1);.X *|
|00003090| 2f 0a 58 4e 55 4d 42 45 | 52 20 2a 0a 58 71 64 65 |/.XNUMBE|R *.Xqde|
|000030a0| 6e 28 71 29 0a 58 09 72 | 65 67 69 73 74 65 72 20 |n(q).X.r|egister |
|000030b0| 4e 55 4d 42 45 52 20 2a | 71 3b 0a 58 7b 0a 58 09 |NUMBER *|q;.X{.X.|
|000030c0| 72 65 67 69 73 74 65 72 | 20 4e 55 4d 42 45 52 20 |register| NUMBER |
|000030d0| 2a 72 3b 0a 58 0a 58 09 | 69 66 20 28 71 69 73 69 |*r;.X.X.|if (qisi|
|000030e0| 6e 74 28 71 29 29 0a 58 | 09 09 72 65 74 75 72 6e |nt(q)).X|..return|
|000030f0| 20 71 6c 69 6e 6b 28 26 | 5f 71 6f 6e 65 5f 29 3b | qlink(&|_qone_);|
|00003100| 0a 58 09 72 20 3d 20 71 | 61 6c 6c 6f 63 28 29 3b |.X.r = q|alloc();|
|00003110| 0a 58 09 7a 63 6f 70 79 | 28 71 2d 3e 64 65 6e 2c |.X.zcopy|(q->den,|
|00003120| 20 26 72 2d 3e 6e 75 6d | 29 3b 0a 58 09 72 65 74 | &r->num|);.X.ret|
|00003130| 75 72 6e 20 72 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |urn r;.X|}.X.X.X/|
|00003140| 2a 0a 58 20 2a 20 52 65 | 74 75 72 6e 20 74 68 65 |*.X * Re|turn the|
|00003150| 20 66 72 61 63 74 69 6f | 6e 61 6c 20 70 61 72 74 | fractio|nal part|
|00003160| 20 6f 66 20 61 20 6e 75 | 6d 62 65 72 2e 0a 58 20 | of a nu|mber..X |
|00003170| 2a 09 71 32 20 3d 20 71 | 66 72 61 63 28 71 31 29 |*.q2 = q|frac(q1)|
|00003180| 3b 0a 58 20 2a 2f 0a 58 | 4e 55 4d 42 45 52 20 2a |;.X */.X|NUMBER *|
|00003190| 0a 58 71 66 72 61 63 28 | 71 29 0a 58 09 72 65 67 |.Xqfrac(|q).X.reg|
|000031a0| 69 73 74 65 72 20 4e 55 | 4d 42 45 52 20 2a 71 3b |ister NU|MBER *q;|
|000031b0| 0a 58 7b 0a 58 09 72 65 | 67 69 73 74 65 72 20 4e |.X{.X.re|gister N|
|000031c0| 55 4d 42 45 52 20 2a 72 | 3b 0a 58 0a 58 09 69 66 |UMBER *r|;.X.X.if|
|000031d0| 20 28 71 69 73 69 6e 74 | 28 71 29 29 0a 58 09 09 | (qisint|(q)).X..|
|000031e0| 72 65 74 75 72 6e 20 71 | 6c 69 6e 6b 28 26 5f 71 |return q|link(&_q|
|000031f0| 7a 65 72 6f 5f 29 3b 0a | 58 09 69 66 20 28 28 71 |zero_);.|X.if ((q|
|00003200| 2d 3e 6e 75 6d 2e 6c 65 | 6e 20 3c 20 71 2d 3e 64 |->num.le|n < q->d|
|00003210| 65 6e 2e 6c 65 6e 29 20 | 7c 7c 20 28 28 71 2d 3e |en.len) ||| ((q->|
|00003220| 6e 75 6d 2e 6c 65 6e 20 | 3d 3d 20 71 2d 3e 64 65 |num.len |== q->de|
|00003230| 6e 2e 6c 65 6e 29 20 26 | 26 0a 58 09 09 28 71 2d |n.len) &|&.X..(q-|
|00003240| 3e 6e 75 6d 2e 76 5b 71 | 2d 3e 6e 75 6d 2e 6c 65 |>num.v[q|->num.le|
|00003250| 6e 20 2d 20 31 5d 20 3c | 20 71 2d 3e 64 65 6e 2e |n - 1] <| q->den.|
|00003260| 76 5b 71 2d 3e 6e 75 6d | 2e 6c 65 6e 20 2d 20 31 |v[q->num|.len - 1|
|00003270| 5d 29 29 29 0a 58 09 09 | 09 72 65 74 75 72 6e 20 |]))).X..|.return |
|00003280| 71 6c 69 6e 6b 28 71 29 | 3b 0a 58 09 72 20 3d 20 |qlink(q)|;.X.r = |
|00003290| 71 61 6c 6c 6f 63 28 29 | 3b 0a 58 09 7a 6d 6f 64 |qalloc()|;.X.zmod|
|000032a0| 28 71 2d 3e 6e 75 6d 2c | 20 71 2d 3e 64 65 6e 2c |(q->num,| q->den,|
|000032b0| 20 26 72 2d 3e 6e 75 6d | 29 3b 0a 58 09 7a 63 6f | &r->num|);.X.zco|
|000032c0| 70 79 28 71 2d 3e 64 65 | 6e 2c 20 26 72 2d 3e 64 |py(q->de|n, &r->d|
|000032d0| 65 6e 29 3b 0a 58 09 72 | 2d 3e 6e 75 6d 2e 73 69 |en);.X.r|->num.si|
|000032e0| 67 6e 20 3d 20 71 2d 3e | 6e 75 6d 2e 73 69 67 6e |gn = q->|num.sign|
|000032f0| 3b 0a 58 09 72 65 74 75 | 72 6e 20 72 3b 0a 58 7d |;.X.retu|rn r;.X}|
|00003300| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 52 65 74 |.X.X.X/*|.X * Ret|
|00003310| 75 72 6e 20 74 68 65 20 | 69 6e 74 65 67 72 61 6c |urn the |integral|
|00003320| 20 70 61 72 74 20 6f 66 | 20 61 20 6e 75 6d 62 65 | part of| a numbe|
|00003330| 72 2e 0a 58 20 2a 09 71 | 32 20 3d 20 71 69 6e 74 |r..X *.q|2 = qint|
|00003340| 28 71 31 29 3b 0a 58 20 | 2a 2f 0a 58 4e 55 4d 42 |(q1);.X |*/.XNUMB|
|00003350| 45 52 20 2a 0a 58 71 69 | 6e 74 28 71 29 0a 58 09 |ER *.Xqi|nt(q).X.|
|00003360| 72 65 67 69 73 74 65 72 | 20 4e 55 4d 42 45 52 20 |register| NUMBER |
|00003370| 2a 71 3b 0a 58 7b 0a 58 | 09 72 65 67 69 73 74 65 |*q;.X{.X|.registe|
|00003380| 72 20 4e 55 4d 42 45 52 | 20 2a 72 3b 0a 58 0a 58 |r NUMBER| *r;.X.X|
|00003390| 09 69 66 20 28 71 69 73 | 69 6e 74 28 71 29 29 0a |.if (qis|int(q)).|
|000033a0| 58 09 09 72 65 74 75 72 | 6e 20 71 6c 69 6e 6b 28 |X..retur|n qlink(|
|000033b0| 71 29 3b 0a 58 09 69 66 | 20 28 28 71 2d 3e 6e 75 |q);.X.if| ((q->nu|
|000033c0| 6d 2e 6c 65 6e 20 3c 20 | 71 2d 3e 64 65 6e 2e 6c |m.len < |q->den.l|
|000033d0| 65 6e 29 20 7c 7c 20 28 | 28 71 2d 3e 6e 75 6d 2e |en) || (|(q->num.|
|000033e0| 6c 65 6e 20 3d 3d 20 71 | 2d 3e 64 65 6e 2e 6c 65 |len == q|->den.le|
|000033f0| 6e 29 20 26 26 0a 58 09 | 09 28 71 2d 3e 6e 75 6d |n) &&.X.|.(q->num|
|00003400| 2e 76 5b 71 2d 3e 6e 75 | 6d 2e 6c 65 6e 20 2d 20 |.v[q->nu|m.len - |
|00003410| 31 5d 20 3c 20 71 2d 3e | 64 65 6e 2e 76 5b 71 2d |1] < q->|den.v[q-|
|00003420| 3e 6e 75 6d 2e 6c 65 6e | 20 2d 20 31 5d 29 29 29 |>num.len| - 1])))|
|00003430| 0a 58 09 09 09 72 65 74 | 75 72 6e 20 71 6c 69 6e |.X...ret|urn qlin|
|00003440| 6b 28 26 5f 71 7a 65 72 | 6f 5f 29 3b 0a 58 09 72 |k(&_qzer|o_);.X.r|
|00003450| 20 3d 20 71 61 6c 6c 6f | 63 28 29 3b 0a 58 09 7a | = qallo|c();.X.z|
|00003460| 71 75 6f 28 71 2d 3e 6e | 75 6d 2c 20 71 2d 3e 64 |quo(q->n|um, q->d|
|00003470| 65 6e 2c 20 26 72 2d 3e | 6e 75 6d 29 3b 0a 58 09 |en, &r->|num);.X.|
|00003480| 72 65 74 75 72 6e 20 72 | 3b 0a 58 7d 0a 58 0a 58 |return r|;.X}.X.X|
|00003490| 0a 58 2f 2a 0a 58 20 2a | 20 43 6f 6d 70 75 74 65 |.X/*.X *| Compute|
|000034a0| 20 74 68 65 20 73 71 75 | 61 72 65 20 6f 66 20 61 | the squ|are of a|
|000034b0| 20 6e 75 6d 62 65 72 2e | 0a 58 20 2a 2f 0a 58 4e | number.|.X */.XN|
|000034c0| 55 4d 42 45 52 20 2a 0a | 58 71 73 71 75 61 72 65 |UMBER *.|Xqsquare|
|000034d0| 28 71 29 0a 58 09 72 65 | 67 69 73 74 65 72 20 4e |(q).X.re|gister N|
|000034e0| 55 4d 42 45 52 20 2a 71 | 3b 0a 58 7b 0a 58 09 5a |UMBER *q|;.X{.X.Z|
|000034f0| 56 41 4c 55 45 20 6e 75 | 6d 2c 20 64 65 6e 3b 0a |VALUE nu|m, den;.|
|00003500| 58 0a 58 09 69 66 20 28 | 71 69 73 7a 65 72 6f 28 |X.X.if (|qiszero(|
|00003510| 71 29 29 0a 58 09 09 72 | 65 74 75 72 6e 20 71 6c |q)).X..r|eturn ql|
|00003520| 69 6e 6b 28 26 5f 71 7a | 65 72 6f 5f 29 3b 0a 58 |ink(&_qz|ero_);.X|
|00003530| 09 69 66 20 28 71 69 73 | 75 6e 69 74 28 71 29 29 |.if (qis|unit(q))|
|00003540| 0a 58 09 09 72 65 74 75 | 72 6e 20 71 6c 69 6e 6b |.X..retu|rn qlink|
|00003550| 28 26 5f 71 6f 6e 65 5f | 29 3b 0a 58 09 6e 75 6d |(&_qone_|);.X.num|
|00003560| 20 3d 20 71 2d 3e 6e 75 | 6d 3b 0a 58 09 64 65 6e | = q->nu|m;.X.den|
|00003570| 20 3d 20 71 2d 3e 64 65 | 6e 3b 0a 58 09 71 20 3d | = q->de|n;.X.q =|
|00003580| 20 71 61 6c 6c 6f 63 28 | 29 3b 0a 58 09 69 66 20 | qalloc(|);.X.if |
|00003590| 28 21 69 73 75 6e 69 74 | 28 6e 75 6d 29 29 0a 58 |(!isunit|(num)).X|
|000035a0| 09 09 7a 73 71 75 61 72 | 65 28 6e 75 6d 2c 20 26 |..zsquar|e(num, &|
|000035b0| 71 2d 3e 6e 75 6d 29 3b | 0a 58 09 69 66 20 28 21 |q->num);|.X.if (!|
|000035c0| 69 73 75 6e 69 74 28 64 | 65 6e 29 29 0a 58 09 09 |isunit(d|en)).X..|
|000035d0| 7a 73 71 75 61 72 65 28 | 64 65 6e 2c 20 26 71 2d |zsquare(|den, &q-|
|000035e0| 3e 64 65 6e 29 3b 0a 58 | 09 72 65 74 75 72 6e 20 |>den);.X|.return |
|000035f0| 71 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |q;.X}.X.|X.X/*.X |
|00003600| 2a 20 53 68 69 66 74 20 | 61 6e 20 69 6e 74 65 67 |* Shift |an integ|
|00003610| 65 72 20 62 79 20 61 20 | 67 69 76 65 6e 20 6e 75 |er by a |given nu|
|00003620| 6d 62 65 72 20 6f 66 20 | 62 69 74 73 2e 20 54 68 |mber of |bits. Th|
|00003630| 69 73 20 6d 75 6c 74 69 | 70 6c 69 65 73 20 74 68 |is multi|plies th|
|00003640| 65 20 6e 75 6d 62 65 72 | 0a 58 20 2a 20 62 79 20 |e number|.X * by |
|00003650| 74 68 65 20 61 70 70 72 | 6f 70 72 69 61 74 65 20 |the appr|opriate |
|00003660| 70 6f 77 65 72 20 6f 66 | 20 74 77 6f 2e 20 20 50 |power of| two.  P|
|00003670| 6f 73 69 74 69 76 65 20 | 6e 75 6d 62 65 72 73 20 |ositive |numbers |
|00003680| 73 68 69 66 74 20 6c 65 | 66 74 2c 20 6e 65 67 61 |shift le|ft, nega|
|00003690| 74 69 76 65 0a 58 20 2a | 20 6f 6e 65 73 20 73 68 |tive.X *| ones sh|
|000036a0| 69 66 74 20 72 69 67 68 | 74 2e 20 20 4c 6f 77 20 |ift righ|t.  Low |
|000036b0| 62 69 74 73 20 61 72 65 | 20 74 72 75 6e 63 61 74 |bits are| truncat|
|000036c0| 65 64 20 77 68 65 6e 20 | 73 68 69 66 74 69 6e 67 |ed when |shifting|
|000036d0| 20 72 69 67 68 74 2e 0a | 58 20 2a 2f 0a 58 4e 55 | right..|X */.XNU|
|000036e0| 4d 42 45 52 20 2a 0a 58 | 71 73 68 69 66 74 28 71 |MBER *.X|qshift(q|
|000036f0| 2c 20 6e 29 0a 58 09 4e | 55 4d 42 45 52 20 2a 71 |, n).X.N|UMBER *q|
|00003700| 3b 0a 58 09 6c 6f 6e 67 | 20 6e 3b 0a 58 7b 0a 58 |;.X.long| n;.X{.X|
|00003710| 09 72 65 67 69 73 74 65 | 72 20 4e 55 4d 42 45 52 |.registe|r NUMBER|
|00003720| 20 2a 72 3b 0a 58 0a 58 | 09 69 66 20 28 71 69 73 | *r;.X.X|.if (qis|
|00003730| 66 72 61 63 28 71 29 29 | 0a 58 09 09 65 72 72 6f |frac(q))|.X..erro|
|00003740| 72 28 22 53 68 69 66 74 | 20 6f 66 20 6e 6f 6e 2d |r("Shift| of non-|
|00003750| 69 6e 74 65 67 65 72 22 | 29 3b 0a 58 09 69 66 20 |integer"|);.X.if |
|00003760| 28 71 69 73 7a 65 72 6f | 28 71 29 20 7c 7c 20 28 |(qiszero|(q) || (|
|00003770| 6e 20 3d 3d 20 30 29 29 | 0a 58 09 09 72 65 74 75 |n == 0))|.X..retu|
|00003780| 72 6e 20 71 6c 69 6e 6b | 28 71 29 3b 0a 58 09 69 |rn qlink|(q);.X.i|
|00003790| 66 20 28 6e 20 3c 3d 20 | 2d 28 71 2d 3e 6e 75 6d |f (n <= |-(q->num|
|000037a0| 2e 6c 65 6e 20 2a 20 42 | 41 53 45 42 29 29 0a 58 |.len * B|ASEB)).X|
|000037b0| 09 09 72 65 74 75 72 6e | 20 71 6c 69 6e 6b 28 26 |..return| qlink(&|
|000037c0| 5f 71 7a 65 72 6f 5f 29 | 3b 0a 58 09 72 20 3d 20 |_qzero_)|;.X.r = |
|000037d0| 71 61 6c 6c 6f 63 28 29 | 3b 0a 58 09 7a 73 68 69 |qalloc()|;.X.zshi|
|000037e0| 66 74 28 71 2d 3e 6e 75 | 6d 2c 20 6e 2c 20 26 72 |ft(q->nu|m, n, &r|
|000037f0| 2d 3e 6e 75 6d 29 3b 0a | 58 09 72 65 74 75 72 6e |->num);.|X.return|
|00003800| 20 72 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 0a 58 | r;.X}.X|.X.X/*.X|
|00003810| 20 2a 20 53 63 61 6c 65 | 20 61 20 6e 75 6d 62 65 | * Scale| a numbe|
|00003820| 72 20 62 79 20 61 20 70 | 6f 77 65 72 20 6f 66 20 |r by a p|ower of |
|00003830| 74 77 6f 2c 20 61 73 20 | 69 6e 3a 0a 58 20 2a 09 |two, as |in:.X *.|
|00003840| 61 6e 73 20 3d 20 71 20 | 2a 20 32 5e 6e 2e 0a 58 |ans = q |* 2^n..X|
|00003850| 20 2a 20 54 68 69 73 20 | 69 73 20 73 69 6d 69 6c | * This |is simil|
|00003860| 61 72 20 74 6f 20 73 68 | 69 66 74 69 6e 67 2c 20 |ar to sh|ifting, |
|00003870| 65 78 63 65 70 74 20 74 | 68 61 74 20 66 72 61 63 |except t|hat frac|
|00003880| 74 69 6f 6e 73 20 77 6f | 72 6b 2e 0a 58 20 2a 2f |tions wo|rk..X */|
|00003890| 0a 58 4e 55 4d 42 45 52 | 20 2a 0a 58 71 73 63 61 |.XNUMBER| *.Xqsca|
|000038a0| 6c 65 28 71 2c 20 70 6f | 77 29 0a 58 09 4e 55 4d |le(q, po|w).X.NUM|
|000038b0| 42 45 52 20 2a 71 3b 0a | 58 09 6c 6f 6e 67 20 70 |BER *q;.|X.long p|
|000038c0| 6f 77 3b 0a 58 7b 0a 58 | 09 6c 6f 6e 67 20 6e 75 |ow;.X{.X|.long nu|
|000038d0| 6d 73 68 69 66 74 2c 20 | 64 65 6e 73 68 69 66 74 |mshift, |denshift|
|000038e0| 2c 20 74 6d 70 3b 0a 58 | 09 4e 55 4d 42 45 52 20 |, tmp;.X|.NUMBER |
|000038f0| 2a 72 3b 0a 58 0a 58 09 | 69 66 20 28 71 69 73 7a |*r;.X.X.|if (qisz|
|00003900| 65 72 6f 28 71 29 20 7c | 7c 20 28 70 6f 77 20 3d |ero(q) ||| (pow =|
|00003910| 3d 20 30 29 29 0a 58 09 | 09 72 65 74 75 72 6e 20 |= 0)).X.|.return |
|00003920| 71 6c 69 6e 6b 28 71 29 | 3b 0a 58 09 69 66 20 28 |qlink(q)|;.X.if (|
|00003930| 28 70 6f 77 20 3e 20 31 | 30 30 30 30 30 30 4c 29 |(pow > 1|000000L)|
|00003940| 20 7c 7c 20 28 70 6f 77 | 20 3c 20 2d 31 30 30 30 | || (pow| < -1000|
|00003950| 30 30 30 4c 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |000L)).X|..error(|
|00003960| 22 56 65 72 79 20 6c 61 | 72 67 65 20 73 63 61 6c |"Very la|rge scal|
|00003970| 65 20 76 61 6c 75 65 22 | 29 3b 0a 58 09 6e 75 6d |e value"|);.X.num|
|00003980| 73 68 69 66 74 20 3d 20 | 69 73 6f 64 64 28 71 2d |shift = |isodd(q-|
|00003990| 3e 6e 75 6d 29 20 3f 20 | 30 20 3a 20 7a 6c 6f 77 |>num) ? |0 : zlow|
|000039a0| 62 69 74 28 71 2d 3e 6e | 75 6d 29 3b 0a 58 09 64 |bit(q->n|um);.X.d|
|000039b0| 65 6e 73 68 69 66 74 20 | 3d 20 69 73 6f 64 64 28 |enshift |= isodd(|
|000039c0| 71 2d 3e 64 65 6e 29 20 | 3f 20 30 20 3a 20 7a 6c |q->den) |? 0 : zl|
|000039d0| 6f 77 62 69 74 28 71 2d | 3e 64 65 6e 29 3b 0a 58 |owbit(q-|>den);.X|
|000039e0| 09 69 66 20 28 70 6f 77 | 20 3e 20 30 29 20 7b 0a |.if (pow| > 0) {.|
|000039f0| 58 09 09 74 6d 70 20 3d | 20 70 6f 77 3b 0a 58 09 |X..tmp =| pow;.X.|
|00003a00| 09 69 66 20 28 74 6d 70 | 20 3e 20 64 65 6e 73 68 |.if (tmp| > densh|
|00003a10| 69 66 74 29 0a 58 09 09 | 74 6d 70 20 3d 20 64 65 |ift).X..|tmp = de|
|00003a20| 6e 73 68 69 66 74 3b 0a | 58 09 09 64 65 6e 73 68 |nshift;.|X..densh|
|00003a30| 69 66 74 20 3d 20 2d 74 | 6d 70 3b 0a 58 09 09 6e |ift = -t|mp;.X..n|
|00003a40| 75 6d 73 68 69 66 74 20 | 3d 20 28 70 6f 77 20 2d |umshift |= (pow -|
|00003a50| 20 74 6d 70 29 3b 0a 58 | 09 7d 20 65 6c 73 65 20 | tmp);.X|.} else |
|00003a60| 7b 0a 58 09 09 70 6f 77 | 20 3d 20 2d 70 6f 77 3b |{.X..pow| = -pow;|
|00003a70| 0a 58 09 09 74 6d 70 20 | 3d 20 70 6f 77 3b 0a 58 |.X..tmp |= pow;.X|
|00003a80| 09 09 69 66 20 28 74 6d | 70 20 3e 20 6e 75 6d 73 |..if (tm|p > nums|
|00003a90| 68 69 66 74 29 0a 58 09 | 09 09 74 6d 70 20 3d 20 |hift).X.|..tmp = |
|00003aa0| 6e 75 6d 73 68 69 66 74 | 3b 0a 58 09 09 6e 75 6d |numshift|;.X..num|
|00003ab0| 73 68 69 66 74 20 3d 20 | 2d 74 6d 70 3b 0a 58 09 |shift = |-tmp;.X.|
|00003ac0| 09 64 65 6e 73 68 69 66 | 74 20 3d 20 28 70 6f 77 |.denshif|t = (pow|
|00003ad0| 20 2d 20 74 6d 70 29 3b | 0a 58 09 7d 0a 58 09 72 | - tmp);|.X.}.X.r|
|00003ae0| 20 3d 20 71 61 6c 6c 6f | 63 28 29 3b 0a 58 09 69 | = qallo|c();.X.i|
|00003af0| 66 20 28 6e 75 6d 73 68 | 69 66 74 29 0a 58 09 09 |f (numsh|ift).X..|
|00003b00| 7a 73 68 69 66 74 28 71 | 2d 3e 6e 75 6d 2c 20 6e |zshift(q|->num, n|
|00003b10| 75 6d 73 68 69 66 74 2c | 20 26 72 2d 3e 6e 75 6d |umshift,| &r->num|
|00003b20| 29 3b 0a 58 09 65 6c 73 | 65 0a 58 09 09 7a 63 6f |);.X.els|e.X..zco|
|00003b30| 70 79 28 71 2d 3e 6e 75 | 6d 2c 20 26 72 2d 3e 6e |py(q->nu|m, &r->n|
|00003b40| 75 6d 29 3b 0a 58 09 69 | 66 20 28 64 65 6e 73 68 |um);.X.i|f (densh|
|00003b50| 69 66 74 29 0a 58 09 09 | 7a 73 68 69 66 74 28 71 |ift).X..|zshift(q|
|00003b60| 2d 3e 64 65 6e 2c 20 64 | 65 6e 73 68 69 66 74 2c |->den, d|enshift,|
|00003b70| 20 26 72 2d 3e 64 65 6e | 29 3b 0a 58 09 65 6c 73 | &r->den|);.X.els|
|00003b80| 65 0a 58 09 09 7a 63 6f | 70 79 28 71 2d 3e 64 65 |e.X..zco|py(q->de|
|00003b90| 6e 2c 20 26 72 2d 3e 64 | 65 6e 29 3b 0a 58 09 72 |n, &r->d|en);.X.r|
|00003ba0| 65 74 75 72 6e 20 72 3b | 0a 58 7d 0a 58 0a 58 0a |eturn r;|.X}.X.X.|
|00003bb0| 58 2f 2a 0a 58 20 2a 20 | 52 65 74 75 72 6e 20 74 |X/*.X * |Return t|
|00003bc0| 68 65 20 6d 69 6e 69 6d | 75 6d 20 6f 66 20 74 77 |he minim|um of tw|
|00003bd0| 6f 20 6e 75 6d 62 65 72 | 73 2e 0a 58 20 2a 2f 0a |o number|s..X */.|
|00003be0| 58 4e 55 4d 42 45 52 20 | 2a 0a 58 71 6d 69 6e 28 |XNUMBER |*.Xqmin(|
|00003bf0| 71 31 2c 20 71 32 29 0a | 58 09 4e 55 4d 42 45 52 |q1, q2).|X.NUMBER|
|00003c00| 20 2a 71 31 2c 20 2a 71 | 32 3b 0a 58 7b 0a 58 09 | *q1, *q|2;.X{.X.|
|00003c10| 69 66 20 28 71 72 65 6c | 28 71 31 2c 20 71 32 29 |if (qrel|(q1, q2)|
|00003c20| 20 3e 20 30 29 0a 58 09 | 09 71 31 20 3d 20 71 32 | > 0).X.|.q1 = q2|
|00003c30| 3b 0a 58 09 72 65 74 75 | 72 6e 20 71 6c 69 6e 6b |;.X.retu|rn qlink|
|00003c40| 28 71 31 29 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a |(q1);.X}|.X.X.X/*|
|00003c50| 0a 58 20 2a 20 52 65 74 | 75 72 6e 20 74 68 65 20 |.X * Ret|urn the |
|00003c60| 6d 61 78 69 6d 75 6d 20 | 6f 66 20 74 77 6f 20 6e |maximum |of two n|
|00003c70| 75 6d 62 65 72 73 2e 0a | 58 20 2a 2f 0a 58 4e 55 |umbers..|X */.XNU|
|00003c80| 4d 42 45 52 20 2a 0a 58 | 71 6d 61 78 28 71 31 2c |MBER *.X|qmax(q1,|
|00003c90| 20 71 32 29 0a 58 09 4e | 55 4d 42 45 52 20 2a 71 | q2).X.N|UMBER *q|
|00003ca0| 31 2c 20 2a 71 32 3b 0a | 58 7b 0a 58 09 69 66 20 |1, *q2;.|X{.X.if |
|00003cb0| 28 71 72 65 6c 28 71 31 | 2c 20 71 32 29 20 3c 20 |(qrel(q1|, q2) < |
|00003cc0| 30 29 0a 58 09 09 71 31 | 20 3d 20 71 32 3b 0a 58 |0).X..q1| = q2;.X|
|00003cd0| 09 72 65 74 75 72 6e 20 | 71 6c 69 6e 6b 28 71 31 |.return |qlink(q1|
|00003ce0| 29 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |);.X}.X.|X.X/*.X |
|00003cf0| 2a 20 50 65 72 66 6f 72 | 6d 20 74 68 65 20 6c 6f |* Perfor|m the lo|
|00003d00| 67 69 63 61 6c 20 4f 52 | 20 6f 66 20 74 77 6f 20 |gical OR| of two |
|00003d10| 69 6e 74 65 67 65 72 73 | 2e 0a 58 20 2a 2f 0a 58 |integers|..X */.X|
|00003d20| 4e 55 4d 42 45 52 20 2a | 0a 58 71 6f 72 28 71 31 |NUMBER *|.Xqor(q1|
|00003d30| 2c 20 71 32 29 0a 58 09 | 4e 55 4d 42 45 52 20 2a |, q2).X.|NUMBER *|
|00003d40| 71 31 2c 20 2a 71 32 3b | 0a 58 7b 0a 58 09 72 65 |q1, *q2;|.X{.X.re|
|00003d50| 67 69 73 74 65 72 20 4e | 55 4d 42 45 52 20 2a 72 |gister N|UMBER *r|
|00003d60| 3b 0a 58 0a 58 09 69 66 | 20 28 71 69 73 66 72 61 |;.X.X.if| (qisfra|
|00003d70| 63 28 71 31 29 20 7c 7c | 20 71 69 73 66 72 61 63 |c(q1) ||| qisfrac|
|00003d80| 28 71 32 29 29 0a 58 09 | 09 65 72 72 6f 72 28 22 |(q2)).X.|.error("|
|00003d90| 4e 6f 6e 2d 69 6e 74 65 | 67 65 72 73 20 66 6f 72 |Non-inte|gers for|
|00003da0| 20 6c 6f 67 69 63 61 6c | 20 6f 72 22 29 3b 0a 58 | logical| or");.X|
|00003db0| 09 69 66 20 28 28 71 31 | 20 3d 3d 20 71 32 29 20 |.if ((q1| == q2) |
|00003dc0| 7c 7c 20 71 69 73 7a 65 | 72 6f 28 71 32 29 29 0a ||| qisze|ro(q2)).|
|00003dd0| 58 09 09 72 65 74 75 72 | 6e 20 71 6c 69 6e 6b 28 |X..retur|n qlink(|
|00003de0| 71 31 29 3b 0a 58 09 69 | 66 20 28 71 69 73 7a 65 |q1);.X.i|f (qisze|
|00003df0| 72 6f 28 71 31 29 29 0a | 58 09 09 72 65 74 75 72 |ro(q1)).|X..retur|
|00003e00| 6e 20 71 6c 69 6e 6b 28 | 71 32 29 3b 0a 58 09 72 |n qlink(|q2);.X.r|
|00003e10| 20 3d 20 71 61 6c 6c 6f | 63 28 29 3b 0a 58 09 7a | = qallo|c();.X.z|
|00003e20| 6f 72 28 71 31 2d 3e 6e | 75 6d 2c 20 71 32 2d 3e |or(q1->n|um, q2->|
|00003e30| 6e 75 6d 2c 20 26 72 2d | 3e 6e 75 6d 29 3b 0a 58 |num, &r-|>num);.X|
|00003e40| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 7d 0a 58 0a |.return |r;.X}.X.|
|00003e50| 58 0a 58 2f 2a 0a 58 20 | 2a 20 50 65 72 66 6f 72 |X.X/*.X |* Perfor|
|00003e60| 6d 20 74 68 65 20 6c 6f | 67 69 63 61 6c 20 41 4e |m the lo|gical AN|
|00003e70| 44 20 6f 66 20 74 77 6f | 20 69 6e 74 65 67 65 72 |D of two| integer|
|00003e80| 73 2e 0a 58 20 2a 2f 0a | 58 4e 55 4d 42 45 52 20 |s..X */.|XNUMBER |
|00003e90| 2a 0a 58 71 61 6e 64 28 | 71 31 2c 20 71 32 29 0a |*.Xqand(|q1, q2).|
|00003ea0| 58 09 4e 55 4d 42 45 52 | 20 2a 71 31 2c 20 2a 71 |X.NUMBER| *q1, *q|
|00003eb0| 32 3b 0a 58 7b 0a 58 09 | 72 65 67 69 73 74 65 72 |2;.X{.X.|register|
|00003ec0| 20 4e 55 4d 42 45 52 20 | 2a 72 3b 0a 58 09 5a 56 | NUMBER |*r;.X.ZV|
|00003ed0| 41 4c 55 45 20 72 65 73 | 3b 0a 58 0a 58 09 69 66 |ALUE res|;.X.X.if|
|00003ee0| 20 28 71 69 73 66 72 61 | 63 28 71 31 29 20 7c 7c | (qisfra|c(q1) |||
|00003ef0| 20 71 69 73 66 72 61 63 | 28 71 32 29 29 0a 58 09 | qisfrac|(q2)).X.|
|00003f00| 09 65 72 72 6f 72 28 22 | 4e 6f 6e 2d 69 6e 74 65 |.error("|Non-inte|
|00003f10| 67 65 72 73 20 66 6f 72 | 20 6c 6f 67 69 63 61 6c |gers for| logical|
|00003f20| 20 61 6e 64 22 29 3b 0a | 58 09 69 66 20 28 71 31 | and");.|X.if (q1|
|00003f30| 20 3d 3d 20 71 32 29 0a | 58 09 09 72 65 74 75 72 | == q2).|X..retur|
|00003f40| 6e 20 71 6c 69 6e 6b 28 | 71 31 29 3b 0a 58 09 69 |n qlink(|q1);.X.i|
|00003f50| 66 20 28 71 69 73 7a 65 | 72 6f 28 71 31 29 20 7c |f (qisze|ro(q1) ||
|00003f60| 7c 20 71 69 73 7a 65 72 | 6f 28 71 32 29 29 0a 58 || qiszer|o(q2)).X|
|00003f70| 09 09 72 65 74 75 72 6e | 20 71 6c 69 6e 6b 28 26 |..return| qlink(&|
|00003f80| 5f 71 7a 65 72 6f 5f 29 | 3b 0a 58 09 7a 61 6e 64 |_qzero_)|;.X.zand|
|00003f90| 28 71 31 2d 3e 6e 75 6d | 2c 20 71 32 2d 3e 6e 75 |(q1->num|, q2->nu|
|00003fa0| 6d 2c 20 26 72 65 73 29 | 3b 0a 58 09 69 66 20 28 |m, &res)|;.X.if (|
|00003fb0| 69 73 7a 65 72 6f 28 72 | 65 73 29 29 20 7b 0a 58 |iszero(r|es)) {.X|
|00003fc0| 09 09 66 72 65 65 68 28 | 72 65 73 2e 76 29 3b 0a |..freeh(|res.v);.|
|00003fd0| 58 09 09 72 65 74 75 72 | 6e 20 71 6c 69 6e 6b 28 |X..retur|n qlink(|
|00003fe0| 26 5f 71 7a 65 72 6f 5f | 29 3b 0a 58 09 7d 0a 58 |&_qzero_|);.X.}.X|
|00003ff0| 09 72 20 3d 20 71 61 6c | 6c 6f 63 28 29 3b 0a 58 |.r = qal|loc();.X|
|00004000| 09 72 2d 3e 6e 75 6d 20 | 3d 20 72 65 73 3b 0a 58 |.r->num |= res;.X|
|00004010| 09 72 65 74 75 72 6e 20 | 72 3b 0a 58 7d 0a 58 0a |.return |r;.X}.X.|
|00004020| 58 0a 58 2f 2a 0a 58 20 | 2a 20 50 65 72 66 6f 72 |X.X/*.X |* Perfor|
|00004030| 6d 20 74 68 65 20 6c 6f | 67 69 63 61 6c 20 58 4f |m the lo|gical XO|
|00004040| 52 20 6f 66 20 74 77 6f | 20 69 6e 74 65 67 65 72 |R of two| integer|
|00004050| 73 2e 0a 58 20 2a 2f 0a | 58 4e 55 4d 42 45 52 20 |s..X */.|XNUMBER |
|00004060| 2a 0a 58 71 78 6f 72 28 | 71 31 2c 20 71 32 29 0a |*.Xqxor(|q1, q2).|
|00004070| 58 09 4e 55 4d 42 45 52 | 20 2a 71 31 2c 20 2a 71 |X.NUMBER| *q1, *q|
|00004080| 32 3b 0a 58 7b 0a 58 09 | 72 65 67 69 73 74 65 72 |2;.X{.X.|register|
|00004090| 20 4e 55 4d 42 45 52 20 | 2a 72 3b 0a 58 09 5a 56 | NUMBER |*r;.X.ZV|
|000040a0| 41 4c 55 45 20 72 65 73 | 3b 0a 58 0a 58 09 69 66 |ALUE res|;.X.X.if|
|000040b0| 20 28 71 69 73 66 72 61 | 63 28 71 31 29 20 7c 7c | (qisfra|c(q1) |||
|000040c0| 20 71 69 73 66 72 61 63 | 28 71 32 29 29 0a 58 09 | qisfrac|(q2)).X.|
|000040d0| 09 65 72 72 6f 72 28 22 | 4e 6f 6e 2d 69 6e 74 65 |.error("|Non-inte|
|000040e0| 67 65 72 73 20 66 6f 72 | 20 6c 6f 67 69 63 61 6c |gers for| logical|
|000040f0| 20 78 6f 72 22 29 3b 0a | 58 09 69 66 20 28 71 31 | xor");.|X.if (q1|
|00004100| 20 3d 3d 20 71 32 29 0a | 58 09 09 72 65 74 75 72 | == q2).|X..retur|
|00004110| 6e 20 71 6c 69 6e 6b 28 | 26 5f 71 7a 65 72 6f 5f |n qlink(|&_qzero_|
|00004120| 29 3b 0a 58 09 69 66 20 | 28 71 69 73 7a 65 72 6f |);.X.if |(qiszero|
|00004130| 28 71 31 29 29 0a 58 09 | 09 72 65 74 75 72 6e 20 |(q1)).X.|.return |
|00004140| 71 6c 69 6e 6b 28 71 32 | 29 3b 0a 58 09 69 66 20 |qlink(q2|);.X.if |
|00004150| 28 71 69 73 7a 65 72 6f | 28 71 32 29 29 0a 58 09 |(qiszero|(q2)).X.|
|00004160| 09 72 65 74 75 72 6e 20 | 71 6c 69 6e 6b 28 71 31 |.return |qlink(q1|
|00004170| 29 3b 0a 58 09 7a 78 6f | 72 28 71 31 2d 3e 6e 75 |);.X.zxo|r(q1->nu|
|00004180| 6d 2c 20 71 32 2d 3e 6e | 75 6d 2c 20 26 72 65 73 |m, q2->n|um, &res|
|00004190| 29 3b 0a 58 09 69 66 20 | 28 69 73 7a 65 72 6f 28 |);.X.if |(iszero(|
|000041a0| 72 65 73 29 29 20 7b 0a | 58 09 09 66 72 65 65 68 |res)) {.|X..freeh|
|000041b0| 28 72 65 73 2e 76 29 3b | 0a 58 09 09 72 65 74 75 |(res.v);|.X..retu|
|000041c0| 72 6e 20 71 6c 69 6e 6b | 28 26 5f 71 7a 65 72 6f |rn qlink|(&_qzero|
|000041d0| 5f 29 3b 0a 58 09 7d 0a | 58 09 72 20 3d 20 71 61 |_);.X.}.|X.r = qa|
|000041e0| 6c 6c 6f 63 28 29 3b 0a | 58 09 72 2d 3e 6e 75 6d |lloc();.|X.r->num|
|000041f0| 20 3d 20 72 65 73 3b 0a | 58 09 72 65 74 75 72 6e | = res;.|X.return|
|00004200| 20 72 3b 0a 58 7d 0a 58 | 0a 58 0a 58 23 69 66 20 | r;.X}.X|.X.X#if |
|00004210| 30 0a 58 2f 2a 0a 58 20 | 2a 20 52 65 74 75 72 6e |0.X/*.X |* Return|
|00004220| 20 74 68 65 20 6e 75 6d | 62 65 72 20 77 68 6f 73 | the num|ber whos|
|00004230| 65 20 62 69 6e 61 72 79 | 20 72 65 70 72 65 73 65 |e binary| represe|
|00004240| 6e 74 61 74 69 6f 6e 20 | 6f 6e 6c 79 20 68 61 73 |ntation |only has|
|00004250| 20 74 68 65 20 73 70 65 | 63 69 66 69 65 64 0a 58 | the spe|cified.X|
|00004260| 20 2a 20 62 69 74 20 73 | 65 74 20 28 63 6f 75 6e | * bit s|et (coun|
|00004270| 74 69 6e 67 20 66 72 6f | 6d 20 7a 65 72 6f 29 2e |ting fro|m zero).|
|00004280| 20 20 54 68 69 73 20 74 | 68 75 73 20 70 72 6f 64 |  This t|hus prod|
|00004290| 75 63 65 73 20 61 20 67 | 69 76 65 6e 20 70 6f 77 |uces a g|iven pow|
|000042a0| 65 72 20 6f 66 20 74 77 | 6f 2e 0a 58 20 2a 2f 0a |er of tw|o..X */.|
|000042b0| 58 4e 55 4d 42 45 52 20 | 2a 0a 58 71 62 69 74 76 |XNUMBER |*.Xqbitv|
|000042c0| 61 6c 75 65 28 6e 29 0a | 58 09 6c 6f 6e 67 20 6e |alue(n).|X.long n|
|000042d0| 3b 0a 58 7b 0a 58 09 72 | 65 67 69 73 74 65 72 20 |;.X{.X.r|egister |
|000042e0| 4e 55 4d 42 45 52 20 2a | 72 3b 0a 58 0a 58 09 69 |NUMBER *|r;.X.X.i|
|000042f0| 66 20 28 6e 20 3c 3d 20 | 30 29 0a 58 09 09 72 65 |f (n <= |0).X..re|
|00004300| 74 75 72 6e 20 71 6c 69 | 6e 6b 28 26 5f 71 6f 6e |turn qli|nk(&_qon|
|00004310| 65 5f 29 3b 0a 58 09 72 | 20 3d 20 71 61 6c 6c 6f |e_);.X.r| = qallo|
|00004320| 63 28 29 3b 0a 58 09 7a | 62 69 74 76 61 6c 75 65 |c();.X.z|bitvalue|
|00004330| 28 6e 2c 20 26 72 2d 3e | 6e 75 6d 29 3b 0a 58 09 |(n, &r->|num);.X.|
|00004340| 72 65 74 75 72 6e 20 72 | 3b 0a 58 7d 0a 58 0a 58 |return r|;.X}.X.X|
|00004350| 0a 58 2f 2a 0a 58 20 2a | 20 54 65 73 74 20 74 6f |.X/*.X *| Test to|
|00004360| 20 73 65 65 20 69 66 20 | 74 68 65 20 73 70 65 63 | see if |the spec|
|00004370| 69 66 69 65 64 20 62 69 | 74 20 6f 66 20 61 20 6e |ified bi|t of a n|
|00004380| 75 6d 62 65 72 20 69 73 | 20 6f 6e 20 28 63 6f 75 |umber is| on (cou|
|00004390| 6e 74 65 64 20 66 72 6f | 6d 20 7a 65 72 6f 29 2e |nted fro|m zero).|
|000043a0| 0a 58 20 2a 20 52 65 74 | 75 72 6e 73 20 54 52 55 |.X * Ret|urns TRU|
|000043b0| 45 20 69 66 20 74 68 65 | 20 62 69 74 20 69 73 20 |E if the| bit is |
|000043c0| 73 65 74 2c 20 6f 72 20 | 46 41 4c 53 45 20 69 66 |set, or |FALSE if|
|000043d0| 20 69 74 20 69 73 20 6e | 6f 74 2e 0a 58 20 2a 09 | it is n|ot..X *.|
|000043e0| 69 20 3d 20 71 62 69 74 | 74 65 73 74 28 71 2c 20 |i = qbit|test(q, |
|000043f0| 6e 29 3b 0a 58 20 2a 2f | 0a 58 42 4f 4f 4c 0a 58 |n);.X */|.XBOOL.X|
|00004400| 71 62 69 74 74 65 73 74 | 28 71 2c 20 6e 29 0a 58 |qbittest|(q, n).X|
|00004410| 09 72 65 67 69 73 74 65 | 72 20 4e 55 4d 42 45 52 |.registe|r NUMBER|
|00004420| 20 2a 71 3b 0a 58 09 6c | 6f 6e 67 20 6e 3b 0a 58 | *q;.X.l|ong n;.X|
|00004430| 7b 0a 58 09 69 6e 74 20 | 78 2c 20 79 3b 0a 58 0a |{.X.int |x, y;.X.|
|00004440| 58 09 69 66 20 28 28 6e | 20 3c 20 30 29 20 7c 7c |X.if ((n| < 0) |||
|00004450| 20 28 6e 20 3e 3d 20 28 | 71 2d 3e 6e 75 6d 2e 6c | (n >= (|q->num.l|
|00004460| 65 6e 20 2a 20 42 41 53 | 45 42 29 29 29 0a 58 09 |en * BAS|EB))).X.|
|00004470| 09 72 65 74 75 72 6e 20 | 46 41 4c 53 45 3b 0a 58 |.return |FALSE;.X|
|00004480| 09 78 20 3d 20 71 2d 3e | 6e 75 6d 2e 76 5b 6e 20 |.x = q->|num.v[n |
|00004490| 2f 20 42 41 53 45 42 5d | 3b 0a 58 09 79 20 3d 20 |/ BASEB]|;.X.y = |
|000044a0| 28 31 20 3c 3c 20 28 6e | 20 25 20 42 41 53 45 42 |(1 << (n| % BASEB|
|000044b0| 29 29 3b 0a 58 09 72 65 | 74 75 72 6e 20 28 28 78 |));.X.re|turn ((x|
|000044c0| 20 26 20 79 29 20 21 3d | 20 30 29 3b 0a 58 7d 0a | & y) !=| 0);.X}.|
|000044d0| 58 23 65 6e 64 69 66 0a | 58 0a 58 0a 58 2f 2a 0a |X#endif.|X.X.X/*.|
|000044e0| 58 20 2a 20 52 65 74 75 | 72 6e 20 74 68 65 20 70 |X * Retu|rn the p|
|000044f0| 72 65 63 69 73 69 6f 6e | 20 6f 66 20 61 20 6e 75 |recision| of a nu|
|00004500| 6d 62 65 72 20 28 75 73 | 75 61 6c 6c 79 20 66 6f |mber (us|ually fo|
|00004510| 72 20 65 78 61 6d 69 6e | 69 6e 67 20 61 6e 20 65 |r examin|ing an e|
|00004520| 70 73 69 6c 6f 6e 20 76 | 61 6c 75 65 29 2e 0a 58 |psilon v|alue)..X|
|00004530| 20 2a 20 54 68 69 73 20 | 69 73 20 74 68 65 20 6c | * This |is the l|
|00004540| 61 72 67 65 73 74 20 70 | 6f 77 65 72 20 6f 66 20 |argest p|ower of |
|00004550| 74 77 6f 20 77 68 6f 73 | 65 20 72 65 63 69 70 72 |two whos|e recipr|
|00004560| 6f 63 61 6c 20 69 73 20 | 6e 6f 74 20 73 6d 61 6c |ocal is |not smal|
|00004570| 6c 65 72 20 69 6e 20 61 | 62 73 6f 6c 75 74 65 0a |ler in a|bsolute.|
|00004580| 58 20 2a 20 76 61 6c 75 | 65 20 74 68 61 6e 20 74 |X * valu|e than t|
|00004590| 68 65 20 73 70 65 63 69 | 66 69 65 64 20 6e 75 6d |he speci|fied num|
|000045a0| 62 65 72 2e 20 20 46 6f | 72 20 65 78 61 6d 70 6c |ber.  Fo|r exampl|
|000045b0| 65 2c 20 71 62 69 74 70 | 72 65 63 28 31 2f 31 30 |e, qbitp|rec(1/10|
|000045c0| 30 29 20 3d 20 36 2e 0a | 58 20 2a 20 4e 75 6d 62 |0) = 6..|X * Numb|
|000045d0| 65 72 73 20 6c 61 72 67 | 65 72 20 74 68 61 6e 20 |ers larg|er than |
|000045e0| 6f 6e 65 20 68 61 76 65 | 20 61 20 70 72 65 63 69 |one have| a preci|
|000045f0| 73 69 6f 6e 20 6f 66 20 | 7a 65 72 6f 2e 0a 58 20 |sion of |zero..X |
|00004600| 2a 2f 0a 58 6c 6f 6e 67 | 0a 58 71 70 72 65 63 69 |*/.Xlong|.Xqpreci|
|00004610| 73 69 6f 6e 28 71 29 0a | 58 09 4e 55 4d 42 45 52 |sion(q).|X.NUMBER|
|00004620| 20 2a 71 3b 0a 58 7b 0a | 58 09 6c 6f 6e 67 20 72 | *q;.X{.|X.long r|
|00004630| 3b 0a 58 0a 58 09 69 66 | 20 28 71 69 73 69 6e 74 |;.X.X.if| (qisint|
|00004640| 28 71 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 30 |(q)).X..|return 0|
|00004650| 3b 0a 58 09 69 66 20 28 | 69 73 75 6e 69 74 28 71 |;.X.if (|isunit(q|
|00004660| 2d 3e 6e 75 6d 29 29 0a | 58 09 09 72 65 74 75 72 |->num)).|X..retur|
|00004670| 6e 20 7a 68 69 67 68 62 | 69 74 28 71 2d 3e 64 65 |n zhighb|it(q->de|
|00004680| 6e 29 3b 0a 58 09 72 20 | 3d 20 7a 68 69 67 68 62 |n);.X.r |= zhighb|
|00004690| 69 74 28 71 2d 3e 64 65 | 6e 29 20 2d 20 7a 68 69 |it(q->de|n) - zhi|
|000046a0| 67 68 62 69 74 28 71 2d | 3e 6e 75 6d 29 20 2d 20 |ghbit(q-|>num) - |
|000046b0| 31 3b 0a 58 09 69 66 20 | 28 72 20 3c 20 30 29 0a |1;.X.if |(r < 0).|
|000046c0| 58 09 09 72 20 3d 20 30 | 3b 0a 58 09 72 65 74 75 |X..r = 0|;.X.retu|
|000046d0| 72 6e 20 72 3b 0a 58 7d | 0a 58 0a 58 0a 58 23 69 |rn r;.X}|.X.X.X#i|
|000046e0| 66 20 30 0a 58 2f 2a 0a | 58 20 2a 20 52 65 74 75 |f 0.X/*.|X * Retu|
|000046f0| 72 6e 20 61 6e 20 69 6e | 74 65 67 65 72 20 69 6e |rn an in|teger in|
|00004700| 64 69 63 61 74 69 6e 67 | 20 74 68 65 20 73 69 67 |dicating| the sig|
|00004710| 6e 20 6f 66 20 61 20 6e | 75 6d 62 65 72 20 28 2d |n of a n|umber (-|
|00004720| 31 2c 20 30 2c 20 6f 72 | 20 31 29 2e 0a 58 20 2a |1, 0, or| 1)..X *|
|00004730| 09 69 20 3d 20 71 74 73 | 74 28 71 29 3b 0a 58 20 |.i = qts|t(q);.X |
|00004740| 2a 2f 0a 58 46 4c 41 47 | 0a 58 71 74 65 73 74 28 |*/.XFLAG|.Xqtest(|
|00004750| 71 29 0a 58 09 72 65 67 | 69 73 74 65 72 20 4e 55 |q).X.reg|ister NU|
|00004760| 4d 42 45 52 20 2a 71 3b | 0a 58 7b 0a 58 09 69 66 |MBER *q;|.X{.X.if|
|00004770| 20 28 21 7a 74 65 73 74 | 28 71 2d 3e 6e 75 6d 29 | (!ztest|(q->num)|
|00004780| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 30 3b 0a 58 |).X..ret|urn 0;.X|
|00004790| 09 69 66 20 28 71 2d 3e | 6e 75 6d 2e 73 69 67 6e |.if (q->|num.sign|
|000047a0| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 2d 31 3b 0a |).X..ret|urn -1;.|
|000047b0| 58 09 72 65 74 75 72 6e | 20 31 3b 0a 58 7d 0a 58 |X.return| 1;.X}.X|
|000047c0| 23 65 6e 64 69 66 0a 58 | 0a 58 0a 58 2f 2a 0a 58 |#endif.X|.X.X/*.X|
|000047d0| 20 2a 20 43 6f 6d 70 61 | 72 65 20 74 77 6f 20 6e | * Compa|re two n|
|000047e0| 75 6d 62 65 72 73 20 61 | 6e 64 20 72 65 74 75 72 |umbers a|nd retur|
|000047f0| 6e 20 61 6e 20 69 6e 74 | 65 67 65 72 20 69 6e 64 |n an int|eger ind|
|00004800| 69 63 61 74 69 6e 67 20 | 74 68 65 69 72 20 72 65 |icating |their re|
|00004810| 6c 61 74 69 76 65 20 73 | 69 7a 65 2e 0a 58 20 2a |lative s|ize..X *|
|00004820| 09 69 20 3d 20 71 72 65 | 6c 28 71 31 2c 20 71 32 |.i = qre|l(q1, q2|
|00004830| 29 3b 0a 58 20 2a 2f 0a | 58 46 4c 41 47 0a 58 71 |);.X */.|XFLAG.Xq|
|00004840| 72 65 6c 28 71 31 2c 20 | 71 32 29 0a 58 09 72 65 |rel(q1, |q2).X.re|
|00004850| 67 69 73 74 65 72 20 4e | 55 4d 42 45 52 20 2a 71 |gister N|UMBER *q|
|00004860| 31 2c 20 2a 71 32 3b 0a | 58 7b 0a 58 09 5a 56 41 |1, *q2;.|X{.X.ZVA|
|00004870| 4c 55 45 20 7a 31 2c 20 | 7a 32 3b 0a 58 09 6c 6f |LUE z1, |z2;.X.lo|
|00004880| 6e 67 20 77 63 31 2c 20 | 77 63 32 3b 0a 58 09 69 |ng wc1, |wc2;.X.i|
|00004890| 6e 74 20 73 69 67 6e 3b | 0a 58 09 69 6e 74 20 7a |nt sign;|.X.int z|
|000048a0| 31 66 20 3d 20 30 2c 20 | 7a 32 66 20 3d 20 30 3b |1f = 0, |z2f = 0;|
|000048b0| 0a 58 0a 58 09 69 66 20 | 28 71 31 20 3d 3d 20 71 |.X.X.if |(q1 == q|
|000048c0| 32 29 0a 58 09 09 72 65 | 74 75 72 6e 20 30 3b 0a |2).X..re|turn 0;.|
|000048d0| 58 09 73 69 67 6e 20 3d | 20 71 32 2d 3e 6e 75 6d |X.sign =| q2->num|
|000048e0| 2e 73 69 67 6e 20 2d 20 | 71 31 2d 3e 6e 75 6d 2e |.sign - |q1->num.|
|000048f0| 73 69 67 6e 3b 0a 58 09 | 69 66 20 28 73 69 67 6e |sign;.X.|if (sign|
|00004900| 29 0a 58 09 09 72 65 74 | 75 72 6e 20 73 69 67 6e |).X..ret|urn sign|
|00004910| 3b 0a 58 09 69 66 20 28 | 71 69 73 7a 65 72 6f 28 |;.X.if (|qiszero(|
|00004920| 71 32 29 29 0a 58 09 09 | 72 65 74 75 72 6e 20 21 |q2)).X..|return !|
|00004930| 71 69 73 7a 65 72 6f 28 | 71 31 29 3b 0a 58 09 69 |qiszero(|q1);.X.i|
|00004940| 66 20 28 71 69 73 7a 65 | 72 6f 28 71 31 29 29 0a |f (qisze|ro(q1)).|
|00004950| 58 09 09 72 65 74 75 72 | 6e 20 2d 31 3b 0a 58 09 |X..retur|n -1;.X.|
|00004960| 2f 2a 0a 58 09 20 2a 20 | 4d 61 6b 65 20 61 20 71 |/*.X. * |Make a q|
|00004970| 75 69 63 6b 20 63 6f 6d | 70 61 72 69 73 6f 6e 20 |uick com|parison |
|00004980| 62 79 20 63 61 6c 63 75 | 6c 61 74 69 6e 67 20 74 |by calcu|lating t|
|00004990| 68 65 20 6e 75 6d 62 65 | 72 20 6f 66 20 77 6f 72 |he numbe|r of wor|
|000049a0| 64 73 20 72 65 73 75 6c | 74 69 6e 67 20 61 73 0a |ds resul|ting as.|
|000049b0| 58 09 20 2a 20 69 66 20 | 77 65 20 6d 75 6c 74 69 |X. * if |we multi|
|000049c0| 70 6c 69 65 64 20 74 68 | 72 6f 75 67 68 20 62 79 |plied th|rough by|
|000049d0| 20 74 68 65 20 64 65 6e | 6f 6d 69 6e 61 74 6f 72 | the den|ominator|
|000049e0| 73 2c 20 61 6e 64 20 74 | 68 65 6e 20 63 6f 6d 70 |s, and t|hen comp|
|000049f0| 61 72 69 6e 67 20 74 68 | 65 0a 58 09 20 2a 20 77 |aring th|e.X. * w|
|00004a00| 6f 72 64 20 63 6f 75 6e | 74 73 2e 0a 58 09 20 2a |ord coun|ts..X. *|
|00004a10| 2f 0a 58 09 73 69 67 6e | 20 3d 20 31 3b 0a 58 09 |/.X.sign| = 1;.X.|
|00004a20| 69 66 20 28 71 69 73 6e | 65 67 28 71 31 29 29 0a |if (qisn|eg(q1)).|
|00004a30| 58 09 09 73 69 67 6e 20 | 3d 20 2d 31 3b 0a 58 09 |X..sign |= -1;.X.|
|00004a40| 77 63 31 20 3d 20 71 31 | 2d 3e 6e 75 6d 2e 6c 65 |wc1 = q1|->num.le|
|00004a50| 6e 20 2b 20 71 32 2d 3e | 64 65 6e 2e 6c 65 6e 3b |n + q2->|den.len;|
|00004a60| 0a 58 09 77 63 32 20 3d | 20 71 32 2d 3e 6e 75 6d |.X.wc2 =| q2->num|
|00004a70| 2e 6c 65 6e 20 2b 20 71 | 31 2d 3e 64 65 6e 2e 6c |.len + q|1->den.l|
|00004a80| 65 6e 3b 0a 58 09 69 66 | 20 28 77 63 31 20 3c 20 |en;.X.if| (wc1 < |
|00004a90| 77 63 32 20 2d 20 31 29 | 0a 58 09 09 72 65 74 75 |wc2 - 1)|.X..retu|
|00004aa0| 72 6e 20 2d 73 69 67 6e | 3b 0a 58 09 69 66 20 28 |rn -sign|;.X.if (|
|00004ab0| 77 63 32 20 3c 20 77 63 | 31 20 2d 20 31 29 0a 58 |wc2 < wc|1 - 1).X|
|00004ac0| 09 09 72 65 74 75 72 6e | 20 73 69 67 6e 3b 0a 58 |..return| sign;.X|
|00004ad0| 09 2f 2a 0a 58 09 20 2a | 20 51 75 69 63 6b 20 63 |./*.X. *| Quick c|
|00004ae0| 68 65 63 6b 20 66 61 69 | 6c 65 64 2c 20 6d 75 73 |heck fai|led, mus|
|00004af0| 74 20 61 63 74 75 61 6c | 6c 79 20 64 6f 20 74 68 |t actual|ly do th|
|00004b00| 65 20 66 75 6c 6c 20 63 | 6f 6d 70 61 72 69 73 6f |e full c|ompariso|
|00004b10| 6e 2e 0a 58 09 20 2a 2f | 0a 58 09 69 66 20 28 69 |n..X. */|.X.if (i|
|00004b20| 73 75 6e 69 74 28 71 32 | 2d 3e 64 65 6e 29 29 0a |sunit(q2|->den)).|
|00004b30| 58 09 09 7a 31 20 3d 20 | 71 31 2d 3e 6e 75 6d 3b |X..z1 = |q1->num;|
|00004b40| 0a 58 09 65 6c 73 65 20 | 69 66 20 28 69 73 6f 6e |.X.else |if (ison|
|00004b50| 65 28 71 31 2d 3e 6e 75 | 6d 29 29 0a 58 09 09 7a |e(q1->nu|m)).X..z|
|00004b60| 31 20 3d 20 71 32 2d 3e | 64 65 6e 3b 0a 58 09 65 |1 = q2->|den;.X.e|
|00004b70| 6c 73 65 20 7b 0a 58 09 | 09 7a 31 66 20 3d 20 31 |lse {.X.|.z1f = 1|
|00004b80| 3b 0a 58 09 09 7a 6d 75 | 6c 28 71 31 2d 3e 6e 75 |;.X..zmu|l(q1->nu|
|00004b90| 6d 2c 20 71 32 2d 3e 64 | 65 6e 2c 20 26 7a 31 29 |m, q2->d|en, &z1)|
|00004ba0| 3b 0a 58 09 7d 0a 58 09 | 69 66 20 28 69 73 75 6e |;.X.}.X.|if (isun|
|00004bb0| 69 74 28 71 31 2d 3e 64 | 65 6e 29 29 0a 58 09 09 |it(q1->d|en)).X..|
|00004bc0| 7a 32 20 3d 20 71 32 2d | 3e 6e 75 6d 3b 0a 58 09 |z2 = q2-|>num;.X.|
|00004bd0| 65 6c 73 65 20 69 66 20 | 28 69 73 6f 6e 65 28 71 |else if |(isone(q|
|00004be0| 32 2d 3e 6e 75 6d 29 29 | 0a 58 09 09 7a 32 20 3d |2->num))|.X..z2 =|
|00004bf0| 20 71 31 2d 3e 64 65 6e | 3b 0a 58 09 65 6c 73 65 | q1->den|;.X.else|
|00004c00| 20 7b 0a 58 09 09 7a 32 | 66 20 3d 20 31 3b 0a 58 | {.X..z2|f = 1;.X|
|00004c10| 09 09 7a 6d 75 6c 28 71 | 32 2d 3e 6e 75 6d 2c 20 |..zmul(q|2->num, |
|00004c20| 71 31 2d 3e 64 65 6e 2c | 20 26 7a 32 29 3b 0a 58 |q1->den,| &z2);.X|
|00004c30| 09 7d 0a 58 09 73 69 67 | 6e 20 3d 20 7a 72 65 6c |.}.X.sig|n = zrel|
|00004c40| 28 7a 31 2c 20 7a 32 29 | 3b 0a 58 09 69 66 20 28 |(z1, z2)|;.X.if (|
|00004c50| 7a 31 66 29 0a 58 09 09 | 66 72 65 65 68 28 7a 31 |z1f).X..|freeh(z1|
|00004c60| 2e 76 29 3b 0a 58 09 69 | 66 20 28 7a 32 66 29 0a |.v);.X.i|f (z2f).|
|00004c70| 58 09 09 66 72 65 65 68 | 28 7a 32 2e 76 29 3b 0a |X..freeh|(z2.v);.|
|00004c80| 58 09 72 65 74 75 72 6e | 20 73 69 67 6e 3b 0a 58 |X.return| sign;.X|
|00004c90| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 43 6f |}.X.X.X/|*.X * Co|
|00004ca0| 6d 70 61 72 65 20 74 77 | 6f 20 6e 75 6d 62 65 72 |mpare tw|o number|
|00004cb0| 73 20 74 6f 20 73 65 65 | 20 69 66 20 74 68 65 79 |s to see| if they|
|00004cc0| 20 61 72 65 20 65 71 75 | 61 6c 2e 0a 58 20 2a 20 | are equ|al..X * |
|00004cd0| 54 68 69 73 20 64 69 66 | 66 65 72 73 20 66 72 6f |This dif|fers fro|
|00004ce0| 6d 20 71 72 65 6c 20 69 | 6e 20 74 68 61 74 20 74 |m qrel i|n that t|
|00004cf0| 68 65 20 6e 75 6d 62 65 | 72 73 20 61 72 65 20 6e |he numbe|rs are n|
|00004d00| 6f 74 20 6f 72 64 65 72 | 65 64 2e 0a 58 20 2a 20 |ot order|ed..X * |
|00004d10| 52 65 74 75 72 6e 73 20 | 54 52 55 45 20 69 66 20 |Returns |TRUE if |
|00004d20| 74 68 65 79 20 64 69 66 | 66 65 72 2e 0a 58 20 2a |they dif|fer..X *|
|00004d30| 2f 0a 58 42 4f 4f 4c 0a | 58 71 63 6d 70 28 71 31 |/.XBOOL.|Xqcmp(q1|
|00004d40| 2c 20 71 32 29 0a 58 09 | 72 65 67 69 73 74 65 72 |, q2).X.|register|
|00004d50| 20 4e 55 4d 42 45 52 20 | 2a 71 31 2c 20 2a 71 32 | NUMBER |*q1, *q2|
|00004d60| 3b 0a 58 7b 0a 58 09 69 | 66 20 28 71 31 20 3d 3d |;.X{.X.i|f (q1 ==|
|00004d70| 20 71 32 29 0a 58 09 09 | 72 65 74 75 72 6e 20 46 | q2).X..|return F|
|00004d80| 41 4c 53 45 3b 0a 58 09 | 69 66 20 28 28 71 31 2d |ALSE;.X.|if ((q1-|
|00004d90| 3e 6e 75 6d 2e 73 69 67 | 6e 20 21 3d 20 71 32 2d |>num.sig|n != q2-|
|00004da0| 3e 6e 75 6d 2e 73 69 67 | 6e 29 20 7c 7c 20 28 71 |>num.sig|n) || (q|
|00004db0| 31 2d 3e 6e 75 6d 2e 6c | 65 6e 20 21 3d 20 71 32 |1->num.l|en != q2|
|00004dc0| 2d 3e 6e 75 6d 2e 6c 65 | 6e 29 20 7c 7c 0a 58 09 |->num.le|n) ||.X.|
|00004dd0| 09 28 71 32 2d 3e 64 65 | 6e 2e 6c 65 6e 20 21 3d |.(q2->de|n.len !=|
|00004de0| 20 71 32 2d 3e 64 65 6e | 2e 6c 65 6e 29 20 7c 7c | q2->den|.len) |||
|00004df0| 20 28 2a 71 31 2d 3e 6e | 75 6d 2e 76 20 21 3d 20 | (*q1->n|um.v != |
|00004e00| 2a 71 32 2d 3e 6e 75 6d | 2e 76 29 20 7c 7c 0a 58 |*q2->num|.v) ||.X|
|00004e10| 09 09 28 2a 71 31 2d 3e | 64 65 6e 2e 76 20 21 3d |..(*q1->|den.v !=|
|00004e20| 20 2a 71 32 2d 3e 64 65 | 6e 2e 76 29 29 0a 58 09 | *q2->de|n.v)).X.|
|00004e30| 09 09 72 65 74 75 72 6e | 20 54 52 55 45 3b 0a 58 |..return| TRUE;.X|
|00004e40| 09 69 66 20 28 7a 63 6d | 70 28 71 31 2d 3e 6e 75 |.if (zcm|p(q1->nu|
|00004e50| 6d 2c 20 71 32 2d 3e 6e | 75 6d 29 29 0a 58 09 09 |m, q2->n|um)).X..|
|00004e60| 72 65 74 75 72 6e 20 54 | 52 55 45 3b 0a 58 09 69 |return T|RUE;.X.i|
|00004e70| 66 20 28 71 69 73 69 6e | 74 28 71 31 29 29 0a 58 |f (qisin|t(q1)).X|
|00004e80| 09 09 72 65 74 75 72 6e | 20 46 41 4c 53 45 3b 0a |..return| FALSE;.|
|00004e90| 58 09 72 65 74 75 72 6e | 20 7a 63 6d 70 28 71 31 |X.return| zcmp(q1|
|00004ea0| 2d 3e 64 65 6e 2c 20 71 | 32 2d 3e 64 65 6e 29 3b |->den, q|2->den);|
|00004eb0| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |.X}.X.X.|X/*.X * |
|00004ec0| 43 6f 6d 70 61 72 65 20 | 61 20 6e 75 6d 62 65 72 |Compare |a number|
|00004ed0| 20 61 67 61 69 6e 73 74 | 20 61 20 6e 6f 72 6d 61 | against| a norma|
|00004ee0| 6c 20 73 6d 61 6c 6c 20 | 69 6e 74 65 67 65 72 2e |l small |integer.|
|00004ef0| 0a 58 20 2a 20 52 65 74 | 75 72 6e 73 20 31 2c 20 |.X * Ret|urns 1, |
|00004f00| 30 2c 20 6f 72 20 2d 31 | 2c 20 61 63 63 6f 72 64 |0, or -1|, accord|
|00004f10| 69 6e 67 20 74 6f 20 77 | 68 65 74 68 65 72 20 74 |ing to w|hether t|
|00004f20| 68 65 20 66 69 72 73 74 | 20 6e 75 6d 62 65 72 20 |he first| number |
|00004f30| 69 73 20 67 72 65 61 74 | 65 72 2c 0a 58 20 2a 20 |is great|er,.X * |
|00004f40| 65 71 75 61 6c 2c 20 6f | 72 20 6c 65 73 73 20 74 |equal, o|r less t|
|00004f50| 68 61 6e 20 74 68 65 20 | 73 65 63 6f 6e 64 20 6e |han the |second n|
|00004f60| 75 6d 62 65 72 2e 0a 58 | 20 2a 09 6e 20 3d 20 71 |umber..X| *.n = q|
|00004f70| 72 65 6c 69 28 71 2c 20 | 6e 29 3b 0a 58 20 2a 2f |reli(q, |n);.X */|
|00004f80| 0a 58 46 4c 41 47 0a 58 | 71 72 65 6c 69 28 71 2c |.XFLAG.X|qreli(q,|
|00004f90| 20 6e 29 0a 58 09 4e 55 | 4d 42 45 52 20 2a 71 3b | n).X.NU|MBER *q;|
|00004fa0| 0a 58 09 6c 6f 6e 67 20 | 6e 3b 0a 58 7b 0a 58 09 |.X.long |n;.X{.X.|
|00004fb0| 69 6e 74 20 73 69 67 6e | 3b 0a 58 09 5a 56 41 4c |int sign|;.X.ZVAL|
|00004fc0| 55 45 20 6e 75 6d 3b 0a | 58 09 48 41 4c 46 20 68 |UE num;.|X.HALF h|
|00004fd0| 32 5b 32 5d 3b 0a 58 09 | 4e 55 4d 42 45 52 20 71 |2[2];.X.|NUMBER q|
|00004fe0| 32 3b 0a 58 0a 58 09 73 | 69 67 6e 20 3d 20 7a 74 |2;.X.X.s|ign = zt|
|00004ff0| 65 73 74 28 71 2d 3e 6e | 75 6d 29 3b 09 09 2f 2a |est(q->n|um);../*|
|00005000| 20 64 6f 20 74 72 69 76 | 69 61 6c 20 73 69 67 6e | do triv|ial sign|
|00005010| 20 63 68 65 63 6b 73 20 | 2a 2f 0a 58 09 69 66 20 | checks |*/.X.if |
|00005020| 28 73 69 67 6e 20 3d 3d | 20 30 29 20 7b 0a 58 09 |(sign ==| 0) {.X.|
|00005030| 09 69 66 20 28 6e 20 3e | 20 30 29 0a 58 09 09 09 |.if (n >| 0).X...|
|00005040| 72 65 74 75 72 6e 20 2d | 31 3b 0a 58 09 09 72 65 |return -|1;.X..re|
|00005050| 74 75 72 6e 20 28 6e 20 | 3c 20 30 29 3b 0a 58 09 |turn (n |< 0);.X.|
|00005060| 7d 0a 58 09 69 66 20 28 | 28 73 69 67 6e 20 3c 20 |}.X.if (|(sign < |
|00005070| 30 29 20 26 26 20 28 6e | 20 3e 3d 20 30 29 29 0a |0) && (n| >= 0)).|
|00005080| 58 09 09 72 65 74 75 72 | 6e 20 2d 31 3b 0a 58 09 |X..retur|n -1;.X.|
|00005090| 69 66 20 28 28 73 69 67 | 6e 20 3e 20 30 29 20 26 |if ((sig|n > 0) &|
|000050a0| 26 20 28 6e 20 3c 3d 20 | 30 29 29 0a 58 09 09 72 |& (n <= |0)).X..r|
|000050b0| 65 74 75 72 6e 20 31 3b | 0a 58 09 6e 20 2a 3d 20 |eturn 1;|.X.n *= |
|000050c0| 73 69 67 6e 3b 0a 58 09 | 69 66 20 28 6e 20 3d 3d |sign;.X.|if (n ==|
|000050d0| 20 31 29 20 7b 09 09 09 | 2f 2a 20 71 75 69 63 6b | 1) {...|/* quick|
|000050e0| 20 63 68 65 63 6b 20 61 | 67 61 69 6e 73 74 20 31 | check a|gainst 1|
|000050f0| 20 6f 72 20 2d 31 20 2a | 2f 0a 58 09 09 6e 75 6d | or -1 *|/.X..num|
|00005100| 20 3d 20 71 2d 3e 6e 75 | 6d 3b 0a 58 09 09 6e 75 | = q->nu|m;.X..nu|
|00005110| 6d 2e 73 69 67 6e 20 3d | 20 30 3b 0a 58 09 09 72 |m.sign =| 0;.X..r|
|00005120| 65 74 75 72 6e 20 28 73 | 69 67 6e 20 2a 20 7a 72 |eturn (s|ign * zr|
|00005130| 65 6c 28 6e 75 6d 2c 20 | 71 2d 3e 64 65 6e 29 29 |el(num, |q->den))|
|00005140| 3b 0a 58 09 7d 0a 58 09 | 6e 75 6d 2e 73 69 67 6e |;.X.}.X.|num.sign|
|00005150| 20 3d 20 28 73 69 67 6e | 20 3c 20 30 29 3b 0a 58 | = (sign| < 0);.X|
|00005160| 09 6e 75 6d 2e 6c 65 6e | 20 3d 20 31 20 2b 20 28 |.num.len| = 1 + (|
|00005170| 6e 20 3e 3d 20 42 41 53 | 45 29 3b 0a 58 09 6e 75 |n >= BAS|E);.X.nu|
|00005180| 6d 2e 76 20 3d 20 68 32 | 3b 0a 58 09 68 32 5b 30 |m.v = h2|;.X.h2[0|
|00005190| 5d 20 3d 20 28 6e 20 26 | 20 42 41 53 45 31 29 3b |] = (n &| BASE1);|
|000051a0| 0a 58 09 68 32 5b 31 5d | 20 3d 20 28 6e 20 3e 3e |.X.h2[1]| = (n >>|
|000051b0| 20 42 41 53 45 42 29 3b | 0a 58 09 69 66 20 28 69 | BASEB);|.X.if (i|
|000051c0| 73 75 6e 69 74 28 71 2d | 3e 64 65 6e 29 29 09 2f |sunit(q-|>den))./|
|000051d0| 2a 20 69 6e 74 65 67 65 | 72 20 63 6f 6d 70 61 72 |* intege|r compar|
|000051e0| 65 20 69 66 20 6e 6f 20 | 64 65 6e 6f 6d 69 6e 61 |e if no |denomina|
|000051f0| 74 6f 72 20 2a 2f 0a 58 | 09 09 72 65 74 75 72 6e |tor */.X|..return|
|00005200| 20 7a 72 65 6c 28 71 2d | 3e 6e 75 6d 2c 20 6e 75 | zrel(q-|>num, nu|
|00005210| 6d 29 3b 0a 58 09 71 32 | 2e 6e 75 6d 20 3d 20 6e |m);.X.q2|.num = n|
|00005220| 75 6d 3b 0a 58 09 71 32 | 2e 64 65 6e 20 3d 20 5f |um;.X.q2|.den = _|
|00005230| 6f 6e 65 5f 3b 0a 58 09 | 71 32 2e 6c 69 6e 6b 73 |one_;.X.|q2.links|
|00005240| 20 3d 20 31 3b 0a 58 09 | 72 65 74 75 72 6e 20 71 | = 1;.X.|return q|
|00005250| 72 65 6c 28 71 2c 20 26 | 71 32 29 3b 09 2f 2a 20 |rel(q, &|q2);./* |
|00005260| 66 75 6c 6c 20 66 72 61 | 63 74 69 6f 6e 61 6c 20 |full fra|ctional |
|00005270| 63 6f 6d 70 61 72 65 20 | 2a 2f 0a 58 7d 0a 58 0a |compare |*/.X}.X.|
|00005280| 58 0a 58 2f 2a 0a 58 20 | 2a 20 43 6f 6d 70 61 72 |X.X/*.X |* Compar|
|00005290| 65 20 61 20 6e 75 6d 62 | 65 72 20 61 67 61 69 6e |e a numb|er again|
|000052a0| 73 74 20 61 20 73 6d 61 | 6c 6c 20 69 6e 74 65 67 |st a sma|ll integ|
|000052b0| 65 72 20 74 6f 20 73 65 | 65 20 69 66 20 74 68 65 |er to se|e if the|
|000052c0| 79 20 61 72 65 20 65 71 | 75 61 6c 2e 0a 58 20 2a |y are eq|ual..X *|
|000052d0| 20 52 65 74 75 72 6e 73 | 20 54 52 55 45 20 69 66 | Returns| TRUE if|
|000052e0| 20 74 68 65 79 20 64 69 | 66 66 65 72 2e 0a 58 20 | they di|ffer..X |
|000052f0| 2a 2f 0a 58 42 4f 4f 4c | 0a 58 71 63 6d 70 69 28 |*/.XBOOL|.Xqcmpi(|
|00005300| 71 2c 20 6e 29 0a 58 09 | 4e 55 4d 42 45 52 20 2a |q, n).X.|NUMBER *|
|00005310| 71 3b 0a 58 09 6c 6f 6e | 67 20 6e 3b 0a 58 7b 0a |q;.X.lon|g n;.X{.|
|00005320| 58 09 6c 6f 6e 67 20 6c | 65 6e 3b 0a 58 0a 58 09 |X.long l|en;.X.X.|
|00005330| 6c 65 6e 20 3d 20 71 2d | 3e 6e 75 6d 2e 6c 65 6e |len = q-|>num.len|
|00005340| 3b 0a 58 09 69 66 20 28 | 28 6c 65 6e 20 3e 20 32 |;.X.if (|(len > 2|
|00005350| 29 20 7c 7c 20 71 69 73 | 66 72 61 63 28 71 29 20 |) || qis|frac(q) |
|00005360| 7c 7c 20 28 71 2d 3e 6e | 75 6d 2e 73 69 67 6e 20 ||| (q->n|um.sign |
|00005370| 21 3d 20 28 6e 20 3c 20 | 30 29 29 29 0a 58 09 09 |!= (n < |0))).X..|
|00005380| 72 65 74 75 72 6e 20 54 | 52 55 45 3b 0a 58 09 69 |return T|RUE;.X.i|
|00005390| 66 20 28 6e 20 3c 20 30 | 29 0a 58 09 09 6e 20 3d |f (n < 0|).X..n =|
|000053a0| 20 2d 6e 3b 0a 58 09 69 | 66 20 28 28 28 48 41 4c | -n;.X.i|f (((HAL|
|000053b0| 46 29 28 6e 29 29 20 21 | 3d 20 71 2d 3e 6e 75 6d |F)(n)) !|= q->num|
|000053c0| 2e 76 5b 30 5d 29 0a 58 | 09 09 72 65 74 75 72 6e |.v[0]).X|..return|
|000053d0| 20 54 52 55 45 3b 0a 58 | 09 6e 20 3d 20 28 28 46 | TRUE;.X|.n = ((F|
|000053e0| 55 4c 4c 29 20 6e 29 20 | 3e 3e 20 42 41 53 45 42 |ULL) n) |>> BASEB|
|000053f0| 3b 0a 58 09 72 65 74 75 | 72 6e 20 28 28 28 6e 20 |;.X.retu|rn (((n |
|00005400| 21 3d 20 30 29 20 21 3d | 20 28 6c 65 6e 20 3d 3d |!= 0) !=| (len ==|
|00005410| 20 32 29 29 20 7c 7c 20 | 28 6e 20 21 3d 20 71 2d | 2)) || |(n != q-|
|00005420| 3e 6e 75 6d 2e 76 5b 31 | 5d 29 29 3b 0a 58 7d 0a |>num.v[1|]));.X}.|
|00005430| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 4e 75 6d 62 |X.X.X/*.|X * Numb|
|00005440| 65 72 20 6e 6f 64 65 20 | 61 6c 6c 6f 63 61 74 69 |er node |allocati|
|00005450| 6f 6e 20 72 6f 75 74 69 | 6e 65 73 0a 58 20 2a 2f |on routi|nes.X */|
|00005460| 0a 58 0a 58 23 64 65 66 | 69 6e 65 09 4e 4e 41 4c |.X.X#def|ine.NNAL|
|00005470| 4c 4f 43 09 31 30 30 30 | 0a 58 0a 58 75 6e 69 6f |LOC.1000|.X.Xunio|
|00005480| 6e 20 61 6c 6c 6f 63 4e | 6f 64 65 20 7b 0a 58 09 |n allocN|ode {.X.|
|00005490| 4e 55 4d 42 45 52 09 6e | 75 6d 3b 0a 58 09 75 6e |NUMBER.n|um;.X.un|
|000054a0| 69 6f 6e 20 61 6c 6c 6f | 63 4e 6f 64 65 09 2a 6c |ion allo|cNode.*l|
|000054b0| 69 6e 6b 3b 0a 58 7d 3b | 0a 58 0a 58 73 74 61 74 |ink;.X};|.X.Xstat|
|000054c0| 69 63 20 75 6e 69 6f 6e | 20 61 6c 6c 6f 63 4e 6f |ic union| allocNo|
|000054d0| 64 65 09 2a 66 72 65 65 | 4e 75 6d 3b 0a 58 0a 58 |de.*free|Num;.X.X|
|000054e0| 0a 58 4e 55 4d 42 45 52 | 20 2a 0a 58 71 61 6c 6c |.XNUMBER| *.Xqall|
|000054f0| 6f 63 28 29 0a 58 7b 0a | 58 09 72 65 67 69 73 74 |oc().X{.|X.regist|
|00005500| 65 72 20 75 6e 69 6f 6e | 20 61 6c 6c 6f 63 4e 6f |er union| allocNo|
|00005510| 64 65 20 2a 74 65 6d 70 | 3b 0a 58 0a 58 09 69 66 |de *temp|;.X.X.if|
|00005520| 20 28 21 66 72 65 65 4e | 75 6d 29 20 7b 0a 58 09 | (!freeN|um) {.X.|
|00005530| 09 66 72 65 65 4e 75 6d | 20 3d 20 28 75 6e 69 6f |.freeNum| = (unio|
|00005540| 6e 20 61 6c 6c 6f 63 4e | 6f 64 65 20 2a 29 0a 58 |n allocN|ode *).X|
|00005550| 09 09 6d 61 6c 6c 6f 63 | 28 73 69 7a 65 6f 66 20 |..malloc|(sizeof |
|00005560| 28 4e 55 4d 42 45 52 29 | 20 2a 20 4e 4e 41 4c 4c |(NUMBER)| * NNALL|
|00005570| 4f 43 29 3b 0a 58 09 09 | 69 66 20 28 66 72 65 65 |OC);.X..|if (free|
|00005580| 4e 75 6d 20 3d 3d 20 4e | 55 4c 4c 29 0a 58 09 09 |Num == N|ULL).X..|
|00005590| 09 65 72 72 6f 72 28 6d | 65 6d 6d 73 67 29 3b 0a |.error(m|emmsg);.|
|000055a0| 58 09 09 74 65 6d 70 20 | 3d 20 66 72 65 65 4e 75 |X..temp |= freeNu|
|000055b0| 6d 3b 0a 58 09 09 77 68 | 69 6c 65 20 28 74 65 6d |m;.X..wh|ile (tem|
|000055c0| 70 20 21 3d 20 66 72 65 | 65 4e 75 6d 20 2b 20 4e |p != fre|eNum + N|
|000055d0| 4e 41 4c 4c 4f 43 20 2d | 20 32 29 20 7b 0a 58 09 |NALLOC -| 2) {.X.|
|000055e0| 09 09 74 65 6d 70 2d 3e | 6c 69 6e 6b 20 3d 20 74 |..temp->|link = t|
|000055f0| 65 6d 70 2b 31 3b 0a 58 | 09 09 09 2b 2b 74 65 6d |emp+1;.X|...++tem|
|00005600| 70 3b 0a 58 09 09 7d 0a | 58 09 7d 0a 58 09 74 65 |p;.X..}.|X.}.X.te|
|00005610| 6d 70 20 3d 20 66 72 65 | 65 4e 75 6d 3b 0a 58 09 |mp = fre|eNum;.X.|
|00005620| 66 72 65 65 4e 75 6d 20 | 3d 20 74 65 6d 70 2d 3e |freeNum |= temp->|
|00005630| 6c 69 6e 6b 3b 0a 58 09 | 74 65 6d 70 2d 3e 6e 75 |link;.X.|temp->nu|
|00005640| 6d 2e 6c 69 6e 6b 73 20 | 3d 20 31 3b 0a 58 09 74 |m.links |= 1;.X.t|
|00005650| 65 6d 70 2d 3e 6e 75 6d | 2e 6e 75 6d 20 3d 20 5f |emp->num|.num = _|
|00005660| 6f 6e 65 5f 3b 0a 58 09 | 74 65 6d 70 2d 3e 6e 75 |one_;.X.|temp->nu|
|00005670| 6d 2e 64 65 6e 20 3d 20 | 5f 6f 6e 65 5f 3b 0a 58 |m.den = |_one_;.X|
|00005680| 09 72 65 74 75 72 6e 20 | 26 74 65 6d 70 2d 3e 6e |.return |&temp->n|
|00005690| 75 6d 3b 0a 58 7d 0a 58 | 0a 58 0a 58 76 6f 69 64 |um;.X}.X|.X.Xvoid|
|000056a0| 0a 58 71 66 72 65 65 6e | 75 6d 28 71 29 0a 58 09 |.Xqfreen|um(q).X.|
|000056b0| 72 65 67 69 73 74 65 72 | 20 4e 55 4d 42 45 52 20 |register| NUMBER |
|000056c0| 2a 71 3b 0a 58 7b 0a 58 | 09 75 6e 69 6f 6e 20 61 |*q;.X{.X|.union a|
|000056d0| 6c 6c 6f 63 4e 6f 64 65 | 20 2a 61 3b 0a 58 0a 58 |llocNode| *a;.X.X|
|000056e0| 09 69 66 20 28 71 20 3d | 3d 20 4e 55 4c 4c 29 0a |.if (q =|= NULL).|
|000056f0| 58 09 09 72 65 74 75 72 | 6e 3b 0a 58 09 66 72 65 |X..retur|n;.X.fre|
|00005700| 65 68 28 71 2d 3e 6e 75 | 6d 2e 76 29 3b 0a 58 09 |eh(q->nu|m.v);.X.|
|00005710| 66 72 65 65 68 28 71 2d | 3e 64 65 6e 2e 76 29 3b |freeh(q-|>den.v);|
|00005720| 0a 58 09 61 20 3d 20 28 | 75 6e 69 6f 6e 20 61 6c |.X.a = (|union al|
|00005730| 6c 6f 63 4e 6f 64 65 20 | 2a 29 20 71 3b 0a 58 09 |locNode |*) q;.X.|
|00005740| 61 2d 3e 6c 69 6e 6b 20 | 3d 20 66 72 65 65 4e 75 |a->link |= freeNu|
|00005750| 6d 3b 0a 58 09 66 72 65 | 65 4e 75 6d 20 3d 20 61 |m;.X.fre|eNum = a|
|00005760| 3b 0a 58 7d 0a 58 0a 58 | 2f 2a 20 45 4e 44 20 43 |;.X}.X.X|/* END C|
|00005770| 4f 44 45 20 2a 2f 0a 45 | 4e 44 5f 4f 46 5f 46 49 |ODE */.E|ND_OF_FI|
|00005780| 4c 45 0a 69 66 20 74 65 | 73 74 20 32 30 31 34 34 |LE.if te|st 20144|
|00005790| 20 2d 6e 65 20 60 77 63 | 20 2d 63 20 3c 27 71 6d | -ne `wc| -c <'qm|
|000057a0| 61 74 68 2e 63 27 60 3b | 20 74 68 65 6e 0a 20 20 |ath.c'`;| then.  |
|000057b0| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 27 |  echo s|har: \"'|
|000057c0| 71 6d 61 74 68 2e 63 27 | 5c 22 20 75 6e 70 61 63 |qmath.c'|\" unpac|
|000057d0| 6b 65 64 20 77 69 74 68 | 20 77 72 6f 6e 67 20 73 |ked with| wrong s|
|000057e0| 69 7a 65 21 0a 66 69 0a | 23 20 65 6e 64 20 6f 66 |ize!.fi.|# end of|
|000057f0| 20 27 71 6d 61 74 68 2e | 63 27 0a 66 69 0a 69 66 | 'qmath.|c'.fi.if|
|00005800| 20 74 65 73 74 20 2d 66 | 20 27 71 74 72 61 6e 73 | test -f| 'qtrans|
|00005810| 2e 63 27 20 2d 61 20 22 | 24 7b 31 7d 22 20 21 3d |.c' -a "|${1}" !=|
|00005820| 20 22 2d 63 22 20 3b 20 | 74 68 65 6e 20 0a 20 20 | "-c" ; |then .  |
|00005830| 65 63 68 6f 20 73 68 61 | 72 3a 20 57 69 6c 6c 20 |echo sha|r: Will |
|00005840| 6e 6f 74 20 63 6c 6f 62 | 62 65 72 20 65 78 69 73 |not clob|ber exis|
|00005850| 74 69 6e 67 20 66 69 6c | 65 20 5c 22 27 71 74 72 |ting fil|e \"'qtr|
|00005860| 61 6e 73 2e 63 27 5c 22 | 0a 65 6c 73 65 0a 65 63 |ans.c'\"|.else.ec|
|00005870| 68 6f 20 73 68 61 72 3a | 20 45 78 74 72 61 63 74 |ho shar:| Extract|
|00005880| 69 6e 67 20 5c 22 27 71 | 74 72 61 6e 73 2e 63 27 |ing \"'q|trans.c'|
|00005890| 5c 22 20 5c 28 32 30 35 | 36 35 20 63 68 61 72 61 |\" \(205|65 chara|
|000058a0| 63 74 65 72 73 5c 29 0a | 73 65 64 20 22 73 2f 5e |cters\).|sed "s/^|
|000058b0| 58 2f 2f 22 20 3e 27 71 | 74 72 61 6e 73 2e 63 27 |X//" >'q|trans.c'|
|000058c0| 20 3c 3c 27 45 4e 44 5f | 4f 46 5f 46 49 4c 45 27 | <<'END_|OF_FILE'|
|000058d0| 0a 58 2f 2a 0a 58 20 2a | 20 43 6f 70 79 72 69 67 |.X/*.X *| Copyrig|
|000058e0| 68 74 20 28 63 29 20 31 | 39 39 32 20 44 61 76 69 |ht (c) 1|992 Davi|
|000058f0| 64 20 49 2e 20 42 65 6c | 6c 0a 58 20 2a 20 50 65 |d I. Bel|l.X * Pe|
|00005900| 72 6d 69 73 73 69 6f 6e | 20 69 73 20 67 72 61 6e |rmission| is gran|
|00005910| 74 65 64 20 74 6f 20 75 | 73 65 2c 20 64 69 73 74 |ted to u|se, dist|
|00005920| 72 69 62 75 74 65 2c 20 | 6f 72 20 6d 6f 64 69 66 |ribute, |or modif|
|00005930| 79 20 74 68 69 73 20 73 | 6f 75 72 63 65 2c 0a 58 |y this s|ource,.X|
|00005940| 20 2a 20 70 72 6f 76 69 | 64 65 64 20 74 68 61 74 | * provi|ded that|
|00005950| 20 74 68 69 73 20 63 6f | 70 79 72 69 67 68 74 20 | this co|pyright |
|00005960| 6e 6f 74 69 63 65 20 72 | 65 6d 61 69 6e 73 20 69 |notice r|emains i|
|00005970| 6e 74 61 63 74 2e 0a 58 | 20 2a 0a 58 20 2a 20 54 |ntact..X| *.X * T|
|00005980| 72 61 6e 73 63 65 6e 64 | 65 6e 74 61 6c 20 66 75 |ranscend|ental fu|
|00005990| 6e 63 74 69 6f 6e 73 20 | 66 6f 72 20 72 65 61 6c |nctions |for real|
|000059a0| 20 6e 75 6d 62 65 72 73 | 2e 0a 58 20 2a 20 54 68 | numbers|..X * Th|
|000059b0| 65 73 65 20 61 72 65 20 | 73 69 6e 2c 20 63 6f 73 |ese are |sin, cos|
|000059c0| 2c 20 65 78 70 2c 20 6c | 6e 2c 20 70 6f 77 65 72 |, exp, l|n, power|
|000059d0| 2c 20 63 6f 73 68 2c 20 | 73 69 6e 68 2e 0a 58 20 |, cosh, |sinh..X |
|000059e0| 2a 2f 0a 58 0a 58 23 69 | 6e 63 6c 75 64 65 20 22 |*/.X.X#i|nclude "|
|000059f0| 6d 61 74 68 2e 68 22 0a | 58 0a 58 42 4f 4f 4c 20 |math.h".|X.XBOOL |
|00005a00| 5f 73 69 6e 69 73 6e 65 | 67 5f 3b 09 2f 2a 20 77 |_sinisne|g_;./* w|
|00005a10| 68 65 74 68 65 72 20 73 | 69 6e 28 78 29 20 3c 20 |hether s|in(x) < |
|00005a20| 30 20 28 73 65 74 20 62 | 79 20 63 6f 73 28 78 29 |0 (set b|y cos(x)|
|00005a30| 29 20 2a 2f 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |) */.X.X|.X/*.X *|
|00005a40| 20 43 61 6c 63 75 6c 61 | 74 65 20 74 68 65 20 63 | Calcula|te the c|
|00005a50| 6f 73 69 6e 65 20 6f 66 | 20 61 20 6e 75 6d 62 65 |osine of| a numbe|
|00005a60| 72 20 77 69 74 68 20 61 | 6e 20 61 63 63 75 72 61 |r with a|n accura|
|00005a70| 63 79 20 77 69 74 68 69 | 6e 20 65 70 73 69 6c 6f |cy withi|n epsilo|
|00005a80| 6e 2e 0a 58 20 2a 20 54 | 68 69 73 20 61 6c 73 6f |n..X * T|his also|
|00005a90| 20 73 61 76 65 73 20 74 | 68 65 20 73 69 67 6e 20 | saves t|he sign |
|00005aa0| 6f 66 20 74 68 65 20 63 | 6f 72 72 65 73 70 6f 6e |of the c|orrespon|
|00005ab0| 64 69 6e 67 20 73 69 6e | 20 66 75 6e 63 74 69 6f |ding sin| functio|
|00005ac0| 6e 2e 0a 58 20 2a 2f 0a | 58 4e 55 4d 42 45 52 20 |n..X */.|XNUMBER |
|00005ad0| 2a 0a 58 71 63 6f 73 28 | 71 2c 20 65 70 73 69 6c |*.Xqcos(|q, epsil|
|00005ae0| 6f 6e 29 0a 58 09 4e 55 | 4d 42 45 52 20 2a 71 2c |on).X.NU|MBER *q,|
|00005af0| 20 2a 65 70 73 69 6c 6f | 6e 3b 0a 58 7b 0a 58 09 | *epsilo|n;.X{.X.|
|00005b00| 4e 55 4d 42 45 52 20 2a | 74 65 72 6d 2c 20 2a 73 |NUMBER *|term, *s|
|00005b10| 75 6d 2c 20 2a 71 73 71 | 2c 20 2a 65 70 73 69 6c |um, *qsq|, *epsil|
|00005b20| 6f 6e 32 2c 20 2a 74 6d | 70 3b 0a 58 09 46 55 4c |on2, *tm|p;.X.FUL|
|00005b30| 4c 20 6e 2c 20 69 3b 0a | 58 09 6c 6f 6e 67 20 73 |L n, i;.|X.long s|
|00005b40| 63 61 6c 65 2c 20 62 69 | 74 73 2c 20 62 69 74 73 |cale, bi|ts, bits|
|00005b50| 32 3b 0a 58 0a 58 09 5f | 73 69 6e 69 73 6e 65 67 |2;.X.X._|sinisneg|
|00005b60| 5f 20 3d 20 71 69 73 6e | 65 67 28 71 29 3b 0a 58 |_ = qisn|eg(q);.X|
|00005b70| 09 69 66 20 28 71 69 73 | 6e 65 67 28 65 70 73 69 |.if (qis|neg(epsi|
|00005b80| 6c 6f 6e 29 20 7c 7c 20 | 71 69 73 7a 65 72 6f 28 |lon) || |qiszero(|
|00005b90| 65 70 73 69 6c 6f 6e 29 | 29 0a 58 09 09 65 72 72 |epsilon)|).X..err|
|00005ba0| 6f 72 28 22 49 6c 6c 65 | 67 61 6c 20 65 70 73 69 |or("Ille|gal epsi|
|00005bb0| 6c 6f 6e 20 76 61 6c 75 | 65 20 66 6f 72 20 63 6f |lon valu|e for co|
|00005bc0| 73 69 6e 65 22 29 3b 0a | 58 09 69 66 20 28 71 69 |sine");.|X.if (qi|
|00005bd0| 73 7a 65 72 6f 28 71 29 | 29 0a 58 09 09 72 65 74 |szero(q)|).X..ret|
|00005be0| 75 72 6e 20 71 6c 69 6e | 6b 28 26 5f 71 6f 6e 65 |urn qlin|k(&_qone|
|00005bf0| 5f 29 3b 0a 58 09 62 69 | 74 73 20 3d 20 71 70 72 |_);.X.bi|ts = qpr|
|00005c00| 65 63 69 73 69 6f 6e 28 | 65 70 73 69 6c 6f 6e 29 |ecision(|epsilon)|
|00005c10| 20 2b 20 31 3b 0a 58 09 | 65 70 73 69 6c 6f 6e 20 | + 1;.X.|epsilon |
|00005c20| 3d 20 71 73 63 61 6c 65 | 28 65 70 73 69 6c 6f 6e |= qscale|(epsilon|
|00005c30| 2c 20 2d 34 4c 29 3b 0a | 58 09 2f 2a 0a 58 09 20 |, -4L);.|X./*.X. |
|00005c40| 2a 20 49 66 20 74 68 65 | 20 61 72 67 75 6d 65 6e |* If the| argumen|
|00005c50| 74 20 69 73 20 6c 61 72 | 67 65 72 20 74 68 61 6e |t is lar|ger than|
|00005c60| 20 6f 6e 65 2c 20 74 68 | 65 6e 20 64 69 76 69 64 | one, th|en divid|
|00005c70| 65 20 69 74 20 62 79 20 | 61 20 70 6f 77 65 72 20 |e it by |a power |
|00005c80| 6f 66 20 74 77 6f 0a 58 | 09 20 2a 20 73 6f 20 74 |of two.X|. * so t|
|00005c90| 68 61 74 20 69 74 20 69 | 73 20 6f 6e 65 20 6f 72 |hat it i|s one or|
|00005ca0| 20 6c 65 73 73 2e 20 20 | 54 68 69 73 20 77 69 6c | less.  |This wil|
|00005cb0| 6c 20 6d 61 6b 65 20 74 | 68 65 20 73 65 72 69 65 |l make t|he serie|
|00005cc0| 73 20 63 6f 6e 76 65 72 | 67 65 20 71 75 69 63 6b |s conver|ge quick|
|00005cd0| 6c 79 2e 0a 58 09 20 2a | 20 57 65 20 77 69 6c 6c |ly..X. *| We will|
|00005ce0| 20 65 78 74 72 61 70 6f | 6c 61 74 65 20 74 68 65 | extrapo|late the|
|00005cf0| 20 72 65 73 75 6c 74 20 | 66 6f 72 20 74 68 65 20 | result |for the |
|00005d00| 6f 72 69 67 69 6e 61 6c | 20 61 72 67 75 6d 65 6e |original| argumen|
|00005d10| 74 20 61 66 74 65 72 77 | 61 72 64 73 2e 0a 58 09 |t afterw|ards..X.|
|00005d20| 20 2a 2f 0a 58 09 73 63 | 61 6c 65 20 3d 20 7a 68 | */.X.sc|ale = zh|
|00005d30| 69 67 68 62 69 74 28 71 | 2d 3e 6e 75 6d 29 20 2d |ighbit(q|->num) -|
|00005d40| 20 7a 68 69 67 68 62 69 | 74 28 71 2d 3e 64 65 6e | zhighbi|t(q->den|
|00005d50| 29 20 2b 20 31 3b 0a 58 | 09 69 66 20 28 73 63 61 |) + 1;.X|.if (sca|
|00005d60| 6c 65 20 3c 20 30 29 0a | 58 09 09 73 63 61 6c 65 |le < 0).|X..scale|
|00005d70| 20 3d 20 30 3b 0a 58 09 | 69 66 20 28 73 63 61 6c | = 0;.X.|if (scal|
|00005d80| 65 20 3e 20 30 29 20 7b | 0a 58 09 09 71 20 3d 20 |e > 0) {|.X..q = |
|00005d90| 71 73 63 61 6c 65 28 71 | 2c 20 2d 73 63 61 6c 65 |qscale(q|, -scale|
|00005da0| 29 3b 0a 58 09 09 74 6d | 70 20 3d 20 71 73 63 61 |);.X..tm|p = qsca|
|00005db0| 6c 65 28 65 70 73 69 6c | 6f 6e 2c 20 2d 73 63 61 |le(epsil|on, -sca|
|00005dc0| 6c 65 29 3b 0a 58 09 09 | 71 66 72 65 65 28 65 70 |le);.X..|qfree(ep|
|00005dd0| 73 69 6c 6f 6e 29 3b 0a | 58 09 09 65 70 73 69 6c |silon);.|X..epsil|
|00005de0| 6f 6e 20 3d 20 74 6d 70 | 3b 0a 58 09 7d 0a 58 09 |on = tmp|;.X.}.X.|
|00005df0| 65 70 73 69 6c 6f 6e 32 | 20 3d 20 71 73 63 61 6c |epsilon2| = qscal|
|00005e00| 65 28 65 70 73 69 6c 6f | 6e 2c 20 2d 34 4c 29 3b |e(epsilo|n, -4L);|
|00005e10| 0a 58 09 71 66 72 65 65 | 28 65 70 73 69 6c 6f 6e |.X.qfree|(epsilon|
|00005e20| 29 3b 0a 58 09 62 69 74 | 73 32 20 3d 20 71 70 72 |);.X.bit|s2 = qpr|
|00005e30| 65 63 69 73 69 6f 6e 28 | 65 70 73 69 6c 6f 6e 32 |ecision(|epsilon2|
|00005e40| 29 20 2b 20 31 30 3b 0a | 58 09 2f 2a 0a 58 09 20 |) + 10;.|X./*.X. |
|00005e50| 2a 20 4e 6f 77 20 75 73 | 65 20 74 68 65 20 54 61 |* Now us|e the Ta|
|00005e60| 79 6c 6f 72 20 73 65 72 | 69 65 73 20 65 78 70 61 |ylor ser|ies expa|
|00005e70| 6e 73 69 6f 6e 20 74 6f | 20 63 61 6c 63 75 6c 61 |nsion to| calcula|
|00005e80| 74 65 20 74 68 65 20 63 | 6f 73 69 6e 65 2e 0a 58 |te the c|osine..X|
|00005e90| 09 20 2a 20 4b 65 65 70 | 20 75 73 69 6e 67 20 61 |. * Keep| using a|
|00005ea0| 70 70 72 6f 78 69 6d 61 | 74 69 6f 6e 73 20 73 6f |pproxima|tions so|
|00005eb0| 20 74 68 61 74 20 74 68 | 65 20 66 72 61 63 74 69 | that th|e fracti|
|00005ec0| 6f 6e 73 20 64 6f 6e 27 | 74 20 67 65 74 20 74 6f |ons don'|t get to|
|00005ed0| 6f 20 6c 61 72 67 65 2e | 0a 58 09 20 2a 2f 0a 58 |o large.|.X. */.X|
|00005ee0| 09 71 73 71 20 3d 20 71 | 73 71 75 61 72 65 28 71 |.qsq = q|square(q|
|00005ef0| 29 3b 0a 58 09 69 66 20 | 28 73 63 61 6c 65 20 3e |);.X.if |(scale >|
|00005f00| 20 30 29 0a 58 09 09 71 | 66 72 65 65 28 71 29 3b | 0).X..q|free(q);|
|00005f10| 0a 58 09 74 65 72 6d 20 | 3d 20 71 6c 69 6e 6b 28 |.X.term |= qlink(|
|00005f20| 26 5f 71 6f 6e 65 5f 29 | 3b 0a 58 09 73 75 6d 20 |&_qone_)|;.X.sum |
|00005f30| 3d 20 71 6c 69 6e 6b 28 | 26 5f 71 6f 6e 65 5f 29 |= qlink(|&_qone_)|
|00005f40| 3b 0a 58 09 6e 20 3d 20 | 30 3b 0a 58 09 77 68 69 |;.X.n = |0;.X.whi|
|00005f50| 6c 65 20 28 71 72 65 6c | 28 74 65 72 6d 2c 20 65 |le (qrel|(term, e|
|00005f60| 70 73 69 6c 6f 6e 32 29 | 20 3e 20 30 29 20 7b 0a |psilon2)| > 0) {.|
|00005f70| 58 09 09 69 20 3d 20 2b | 2b 6e 3b 0a 58 09 09 69 |X..i = +|+n;.X..i|
|00005f80| 20 2a 3d 20 2b 2b 6e 3b | 0a 58 09 09 74 6d 70 20 | *= ++n;|.X..tmp |
|00005f90| 3d 20 71 6d 75 6c 28 74 | 65 72 6d 2c 20 71 73 71 |= qmul(t|erm, qsq|
|00005fa0| 29 3b 0a 58 09 09 71 66 | 72 65 65 28 74 65 72 6d |);.X..qf|ree(term|
|00005fb0| 29 3b 0a 58 09 09 74 65 | 72 6d 20 3d 20 71 64 69 |);.X..te|rm = qdi|
|00005fc0| 76 69 28 74 6d 70 2c 20 | 28 6c 6f 6e 67 29 20 69 |vi(tmp, |(long) i|
|00005fd0| 29 3b 0a 58 09 09 71 66 | 72 65 65 28 74 6d 70 29 |);.X..qf|ree(tmp)|
|00005fe0| 3b 0a 58 09 09 74 6d 70 | 20 3d 20 71 62 72 6f 75 |;.X..tmp| = qbrou|
|00005ff0| 6e 64 28 74 65 72 6d 2c | 20 62 69 74 73 32 29 3b |nd(term,| bits2);|
|00006000| 0a 58 09 09 71 66 72 65 | 65 28 74 65 72 6d 29 3b |.X..qfre|e(term);|
|00006010| 0a 58 09 09 74 65 72 6d | 20 3d 20 74 6d 70 3b 0a |.X..term| = tmp;.|
|00006020| 58 09 09 69 66 20 28 6e | 20 26 20 32 29 0a 58 09 |X..if (n| & 2).X.|
|00006030| 09 09 74 6d 70 20 3d 20 | 71 73 75 62 28 73 75 6d |..tmp = |qsub(sum|
|00006040| 2c 20 74 65 72 6d 29 3b | 0a 58 09 09 65 6c 73 65 |, term);|.X..else|
|00006050| 0a 58 09 09 09 74 6d 70 | 20 3d 20 71 61 64 64 28 |.X...tmp| = qadd(|
|00006060| 73 75 6d 2c 20 74 65 72 | 6d 29 3b 0a 58 09 09 71 |sum, ter|m);.X..q|
|00006070| 66 72 65 65 28 73 75 6d | 29 3b 0a 58 09 09 73 75 |free(sum|);.X..su|
|00006080| 6d 20 3d 20 71 62 72 6f | 75 6e 64 28 74 6d 70 2c |m = qbro|und(tmp,|
|00006090| 20 62 69 74 73 32 29 3b | 0a 58 09 09 71 66 72 65 | bits2);|.X..qfre|
|000060a0| 65 28 74 6d 70 29 3b 0a | 58 09 7d 0a 58 09 71 66 |e(tmp);.|X.}.X.qf|
|000060b0| 72 65 65 28 74 65 72 6d | 29 3b 0a 58 09 71 66 72 |ree(term|);.X.qfr|
|000060c0| 65 65 28 71 73 71 29 3b | 0a 58 09 71 66 72 65 65 |ee(qsq);|.X.qfree|
|000060d0| 28 65 70 73 69 6c 6f 6e | 32 29 3b 0a 58 09 2f 2a |(epsilon|2);.X./*|
|000060e0| 0a 58 09 20 2a 20 4e 6f | 77 20 73 63 61 6c 65 20 |.X. * No|w scale |
|000060f0| 62 61 63 6b 20 75 70 20 | 74 6f 20 74 68 65 20 6f |back up |to the o|
|00006100| 72 69 67 69 6e 61 6c 20 | 76 61 6c 75 65 20 6f 66 |riginal |value of|
|00006110| 20 78 20 62 79 20 75 73 | 69 6e 67 20 74 68 65 20 | x by us|ing the |
|00006120| 66 6f 72 6d 75 6c 61 3a | 0a 58 09 20 2a 09 63 6f |formula:|.X. *.co|
|00006130| 73 28 32 20 2a 20 78 29 | 20 3d 20 32 20 2a 20 28 |s(2 * x)| = 2 * (|
|00006140| 63 6f 73 28 78 29 20 5e | 20 32 29 20 2d 20 31 2e |cos(x) ^| 2) - 1.|
|00006150| 0a 58 09 20 2a 2f 0a 58 | 09 77 68 69 6c 65 20 28 |.X. */.X|.while (|
|00006160| 2d 2d 73 63 61 6c 65 20 | 3e 3d 20 30 29 20 7b 0a |--scale |>= 0) {.|
|00006170| 58 09 09 69 66 20 28 71 | 69 73 6e 65 67 28 73 75 |X..if (q|isneg(su|
|00006180| 6d 29 29 0a 58 09 09 09 | 5f 73 69 6e 69 73 6e 65 |m)).X...|_sinisne|
|00006190| 67 5f 20 3d 20 21 5f 73 | 69 6e 69 73 6e 65 67 5f |g_ = !_s|inisneg_|
|000061a0| 3b 0a 58 09 09 74 6d 70 | 20 3d 20 71 73 71 75 61 |;.X..tmp| = qsqua|
|000061b0| 72 65 28 73 75 6d 29 3b | 0a 58 09 09 71 66 72 65 |re(sum);|.X..qfre|
|000061c0| 65 28 73 75 6d 29 3b 0a | 58 09 09 73 75 6d 20 3d |e(sum);.|X..sum =|
|000061d0| 20 71 73 63 61 6c 65 28 | 74 6d 70 2c 20 31 4c 29 | qscale(|tmp, 1L)|
|000061e0| 3b 0a 58 09 09 71 66 72 | 65 65 28 74 6d 70 29 3b |;.X..qfr|ee(tmp);|
|000061f0| 0a 58 09 09 74 6d 70 20 | 3d 20 71 64 65 63 28 73 |.X..tmp |= qdec(s|
|00006200| 75 6d 29 3b 0a 58 09 09 | 71 66 72 65 65 28 73 75 |um);.X..|qfree(su|
|00006210| 6d 29 3b 0a 58 09 09 73 | 75 6d 20 3d 20 71 62 72 |m);.X..s|um = qbr|
|00006220| 6f 75 6e 64 28 74 6d 70 | 2c 20 62 69 74 73 32 29 |ound(tmp|, bits2)|
|00006230| 3b 0a 58 09 09 71 66 72 | 65 65 28 74 6d 70 29 3b |;.X..qfr|ee(tmp);|
|00006240| 0a 58 09 7d 0a 58 09 74 | 6d 70 20 3d 20 71 62 72 |.X.}.X.t|mp = qbr|
|00006250| 6f 75 6e 64 28 73 75 6d | 2c 20 62 69 74 73 29 3b |ound(sum|, bits);|
|00006260| 0a 58 09 71 66 72 65 65 | 28 73 75 6d 29 3b 0a 58 |.X.qfree|(sum);.X|
|00006270| 09 72 65 74 75 72 6e 20 | 74 6d 70 3b 0a 58 7d 0a |.return |tmp;.X}.|
|00006280| 58 0a 58 0a 58 2f 2a 0a | 58 20 2a 20 43 61 6c 63 |X.X.X/*.|X * Calc|
|00006290| 75 6c 61 74 65 20 74 68 | 65 20 73 69 6e 65 20 6f |ulate th|e sine o|
|000062a0| 66 20 61 20 6e 75 6d 62 | 65 72 20 77 69 74 68 20 |f a numb|er with |
|000062b0| 61 6e 20 61 63 63 75 72 | 61 63 79 20 77 69 74 68 |an accur|acy with|
|000062c0| 69 6e 20 65 70 73 69 6c | 6f 6e 2e 0a 58 20 2a 20 |in epsil|on..X * |
|000062d0| 54 68 69 73 20 69 73 20 | 63 61 6c 63 75 6c 61 74 |This is |calculat|
|000062e0| 65 64 20 75 73 69 6e 67 | 20 74 68 65 20 66 6f 72 |ed using| the for|
|000062f0| 6d 75 6c 61 3a 0a 58 20 | 2a 09 73 69 6e 28 78 29 |mula:.X |*.sin(x)|
|00006300| 5e 32 20 2b 20 63 6f 73 | 28 78 29 5e 32 20 3d 20 |^2 + cos|(x)^2 = |
|00006310| 31 2e 0a 58 20 2a 20 54 | 68 65 20 6f 6e 6c 79 20 |1..X * T|he only |
|00006320| 74 72 69 63 6b 79 20 62 | 69 74 20 69 73 20 72 65 |tricky b|it is re|
|00006330| 73 6f 6c 76 69 6e 67 20 | 74 68 65 20 73 69 67 6e |solving |the sign|
|00006340| 20 6f 66 20 74 68 65 20 | 72 65 73 75 6c 74 2e 0a | of the |result..|
|00006350| 58 20 2a 20 46 75 74 75 | 72 65 3a 20 55 73 65 20 |X * Futu|re: Use |
|00006360| 73 69 6e 28 33 2a 78 29 | 20 3d 20 33 2a 73 69 6e |sin(3*x)| = 3*sin|
|00006370| 28 78 29 20 2d 20 34 2a | 73 69 6e 28 78 29 5e 33 |(x) - 4*|sin(x)^3|
|00006380| 2e 0a 58 20 2a 2f 0a 58 | 4e 55 4d 42 45 52 20 2a |..X */.X|NUMBER *|
|00006390| 0a 58 71 73 69 6e 28 71 | 2c 20 65 70 73 69 6c 6f |.Xqsin(q|, epsilo|
|000063a0| 6e 29 0a 58 09 4e 55 4d | 42 45 52 20 2a 71 2c 20 |n).X.NUM|BER *q, |
|000063b0| 2a 65 70 73 69 6c 6f 6e | 3b 0a 58 7b 0a 58 09 4e |*epsilon|;.X{.X.N|
|000063c0| 55 4d 42 45 52 20 2a 74 | 6d 70 31 2c 20 2a 74 6d |UMBER *t|mp1, *tm|
|000063d0| 70 32 2c 20 2a 65 70 73 | 69 6c 6f 6e 32 3b 0a 58 |p2, *eps|ilon2;.X|
|000063e0| 0a 58 09 69 66 20 28 71 | 69 73 6e 65 67 28 65 70 |.X.if (q|isneg(ep|
|000063f0| 73 69 6c 6f 6e 29 20 7c | 7c 20 71 69 73 7a 65 72 |silon) ||| qiszer|
+--------+-------------------------+-------------------------+--------+--------+