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

< 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 34 30 |Subject:| v26i040|
|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 31 | 34 2f 32 31 0a 53 65 6e |r, Part1|4/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 34 30 |ume 26, |Issue 40|
|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 | 31 34 0a 0a 23 21 20 2f |alc/part|14..#! /|
|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 31 34 20 28 6f 66 20 | 32 31 29 2e 22 0a 23 20 | 14 (of |21).".# |
|000002e0| 43 6f 6e 74 65 6e 74 73 | 3a 20 20 6d 61 74 66 75 |Contents|:  matfu|
|000002f0| 6e 63 2e 63 0a 23 20 57 | 72 61 70 70 65 64 20 62 |nc.c.# W|rapped b|
|00000300| 79 20 64 62 65 6c 6c 40 | 65 6c 6d 20 6f 6e 20 54 |y dbell@|elm on T|
|00000310| 75 65 20 46 65 62 20 32 | 35 20 31 35 3a 32 31 3a |ue Feb 2|5 15:21:|
|00000320| 31 32 20 31 39 39 32 0a | 50 41 54 48 3d 2f 62 69 |12 1992.|PATH=/bi|
|00000330| 6e 3a 2f 75 73 72 2f 62 | 69 6e 3a 2f 75 73 72 2f |n:/usr/b|in:/usr/|
|00000340| 75 63 62 20 3b 20 65 78 | 70 6f 72 74 20 50 41 54 |ucb ; ex|port PAT|
|00000350| 48 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 6d 61 |H.if tes|t -f 'ma|
|00000360| 74 66 75 6e 63 2e 63 27 | 20 2d 61 20 22 24 7b 31 |tfunc.c'| -a "${1|
|00000370| 7d 22 20 21 3d 20 22 2d | 63 22 20 3b 20 74 68 65 |}" != "-|c" ; the|
|00000380| 6e 20 0a 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |n .  ech|o shar: |
|00000390| 57 69 6c 6c 20 6e 6f 74 | 20 63 6c 6f 62 62 65 72 |Will not| clobber|
|000003a0| 20 65 78 69 73 74 69 6e | 67 20 66 69 6c 65 20 5c | existin|g file \|
|000003b0| 22 27 6d 61 74 66 75 6e | 63 2e 63 27 5c 22 0a 65 |"'matfun|c.c'\".e|
|000003c0| 6c 73 65 0a 65 63 68 6f | 20 73 68 61 72 3a 20 45 |lse.echo| shar: E|
|000003d0| 78 74 72 61 63 74 69 6e | 67 20 5c 22 27 6d 61 74 |xtractin|g \"'mat|
|000003e0| 66 75 6e 63 2e 63 27 5c | 22 20 5c 28 32 38 30 37 |func.c'\|" \(2807|
|000003f0| 35 20 63 68 61 72 61 63 | 74 65 72 73 5c 29 0a 73 |5 charac|ters\).s|
|00000400| 65 64 20 22 73 2f 5e 58 | 2f 2f 22 20 3e 27 6d 61 |ed "s/^X|//" >'ma|
|00000410| 74 66 75 6e 63 2e 63 27 | 20 3c 3c 27 45 4e 44 5f |tfunc.c'| <<'END_|
|00000420| 4f 46 5f 46 49 4c 45 27 | 0a 58 2f 2a 0a 58 20 2a |OF_FILE'|.X/*.X *|
|00000430| 20 43 6f 70 79 72 69 67 | 68 74 20 28 63 29 20 31 | Copyrig|ht (c) 1|
|00000440| 39 39 32 20 44 61 76 69 | 64 20 49 2e 20 42 65 6c |992 Davi|d I. Bel|
|00000450| 6c 0a 58 20 2a 20 50 65 | 72 6d 69 73 73 69 6f 6e |l.X * Pe|rmission|
|00000460| 20 69 73 20 67 72 61 6e | 74 65 64 20 74 6f 20 75 | is gran|ted to u|
|00000470| 73 65 2c 20 64 69 73 74 | 72 69 62 75 74 65 2c 20 |se, dist|ribute, |
|00000480| 6f 72 20 6d 6f 64 69 66 | 79 20 74 68 69 73 20 73 |or modif|y this s|
|00000490| 6f 75 72 63 65 2c 0a 58 | 20 2a 20 70 72 6f 76 69 |ource,.X| * provi|
|000004a0| 64 65 64 20 74 68 61 74 | 20 74 68 69 73 20 63 6f |ded that| this co|
|000004b0| 70 79 72 69 67 68 74 20 | 6e 6f 74 69 63 65 20 72 |pyright |notice r|
|000004c0| 65 6d 61 69 6e 73 20 69 | 6e 74 61 63 74 2e 0a 58 |emains i|ntact..X|
|000004d0| 20 2a 0a 58 20 2a 20 45 | 78 74 65 6e 64 65 64 20 | *.X * E|xtended |
|000004e0| 70 72 65 63 69 73 69 6f | 6e 20 72 61 74 69 6f 6e |precisio|n ration|
|000004f0| 61 6c 20 61 72 69 74 68 | 6d 65 74 69 63 20 6d 61 |al arith|metic ma|
|00000500| 74 72 69 78 20 66 75 6e | 63 74 69 6f 6e 73 2e 0a |trix fun|ctions..|
|00000510| 58 20 2a 20 4d 61 74 72 | 69 63 65 73 20 63 61 6e |X * Matr|ices can|
|00000520| 20 63 6f 6e 74 61 69 6e | 20 61 72 62 69 74 72 61 | contain| arbitra|
|00000530| 72 79 20 74 79 70 65 73 | 20 6f 66 20 65 6c 65 6d |ry types| of elem|
|00000540| 65 6e 74 73 2e 0a 58 20 | 2a 2f 0a 58 0a 58 23 69 |ents..X |*/.X.X#i|
|00000550| 6e 63 6c 75 64 65 20 22 | 63 61 6c 63 2e 68 22 0a |nclude "|calc.h".|
|00000560| 58 0a 58 0a 58 23 69 66 | 20 30 0a 58 73 74 61 74 |X.X.X#if| 0.Xstat|
|00000570| 69 63 20 63 68 61 72 20 | 2a 61 62 6f 72 74 6d 73 |ic char |*abortms|
|00000580| 67 20 3d 20 22 43 61 6c | 63 75 6c 61 74 69 6f 6e |g = "Cal|culation|
|00000590| 20 61 62 6f 72 74 65 64 | 22 3b 0a 58 73 74 61 74 | aborted|";.Xstat|
|000005a0| 69 63 20 63 68 61 72 20 | 2a 6d 65 6d 6d 73 67 20 |ic char |*memmsg |
|000005b0| 3d 20 22 4e 6f 74 20 65 | 6e 6f 75 67 68 20 6d 65 |= "Not e|nough me|
|000005c0| 6d 6f 72 79 22 3b 0a 58 | 23 65 6e 64 69 66 0a 58 |mory";.X|#endif.X|
|000005d0| 0a 58 0a 58 73 74 61 74 | 69 63 20 76 6f 69 64 20 |.X.Xstat|ic void |
|000005e0| 6d 61 74 73 77 61 70 72 | 6f 77 28 29 2c 20 6d 61 |matswapr|ow(), ma|
|000005f0| 74 73 75 62 72 6f 77 28 | 29 2c 20 6d 61 74 6d 75 |tsubrow(|), matmu|
|00000600| 6c 72 6f 77 28 29 3b 0a | 58 23 69 66 20 30 0a 58 |lrow();.|X#if 0.X|
|00000610| 73 74 61 74 69 63 20 76 | 6f 69 64 20 6d 61 74 61 |static v|oid mata|
|00000620| 64 64 72 6f 77 28 29 3b | 0a 58 23 65 6e 64 69 66 |ddrow();|.X#endif|
|00000630| 0a 58 0a 58 73 74 61 74 | 69 63 20 4d 41 54 52 49 |.X.Xstat|ic MATRI|
|00000640| 58 20 2a 6d 61 74 69 64 | 65 6e 74 28 29 3b 0a 58 |X *matid|ent();.X|
|00000650| 0a 58 0a 58 0a 58 2f 2a | 0a 58 20 2a 20 41 64 64 |.X.X.X/*|.X * Add|
|00000660| 20 74 77 6f 20 63 6f 6d | 70 61 74 69 62 6c 65 20 | two com|patible |
|00000670| 6d 61 74 72 69 63 65 73 | 2e 0a 58 20 2a 2f 0a 58 |matrices|..X */.X|
|00000680| 4d 41 54 52 49 58 20 2a | 0a 58 6d 61 74 61 64 64 |MATRIX *|.Xmatadd|
|00000690| 28 6d 31 2c 20 6d 32 29 | 0a 58 09 4d 41 54 52 49 |(m1, m2)|.X.MATRI|
|000006a0| 58 20 2a 6d 31 2c 20 2a | 6d 32 3b 0a 58 7b 0a 58 |X *m1, *|m2;.X{.X|
|000006b0| 09 69 6e 74 20 64 69 6d | 3b 0a 58 0a 58 09 6c 6f |.int dim|;.X.X.lo|
|000006c0| 6e 67 20 6d 69 6e 31 2c | 20 6d 69 6e 32 2c 20 6d |ng min1,| min2, m|
|000006d0| 61 78 31 2c 20 6d 61 78 | 32 2c 20 69 6e 64 65 78 |ax1, max|2, index|
|000006e0| 3b 0a 58 09 56 41 4c 55 | 45 20 2a 76 31 2c 20 2a |;.X.VALU|E *v1, *|
|000006f0| 76 32 2c 20 2a 76 72 65 | 73 3b 0a 58 09 4d 41 54 |v2, *vre|s;.X.MAT|
|00000700| 52 49 58 20 2a 72 65 73 | 3b 0a 58 09 4d 41 54 52 |RIX *res|;.X.MATR|
|00000710| 49 58 20 74 6d 70 3b 0a | 58 0a 58 09 69 66 20 28 |IX tmp;.|X.X.if (|
|00000720| 6d 31 2d 3e 6d 5f 64 69 | 6d 20 21 3d 20 6d 32 2d |m1->m_di|m != m2-|
|00000730| 3e 6d 5f 64 69 6d 29 0a | 58 09 09 65 72 72 6f 72 |>m_dim).|X..error|
|00000740| 28 22 49 6e 63 6f 6d 70 | 61 74 69 62 6c 65 20 6d |("Incomp|atible m|
|00000750| 61 74 72 69 78 20 64 69 | 6d 65 6e 73 69 6f 6e 73 |atrix di|mensions|
|00000760| 20 66 6f 72 20 61 64 64 | 22 29 3b 0a 58 09 74 6d | for add|");.X.tm|
|00000770| 70 2e 6d 5f 64 69 6d 20 | 3d 20 6d 31 2d 3e 6d 5f |p.m_dim |= m1->m_|
|00000780| 64 69 6d 3b 0a 58 09 74 | 6d 70 2e 6d 5f 73 69 7a |dim;.X.t|mp.m_siz|
|00000790| 65 20 3d 20 6d 31 2d 3e | 6d 5f 73 69 7a 65 3b 0a |e = m1->|m_size;.|
|000007a0| 58 09 66 6f 72 20 28 64 | 69 6d 20 3d 20 30 3b 20 |X.for (d|im = 0; |
|000007b0| 64 69 6d 20 3c 20 6d 31 | 2d 3e 6d 5f 64 69 6d 3b |dim < m1|->m_dim;|
|000007c0| 20 64 69 6d 2b 2b 29 20 | 7b 0a 58 09 09 6d 69 6e | dim++) |{.X..min|
|000007d0| 31 20 3d 20 6d 31 2d 3e | 6d 5f 6d 69 6e 5b 64 69 |1 = m1->|m_min[di|
|000007e0| 6d 5d 3b 0a 58 09 09 6d | 61 78 31 20 3d 20 6d 31 |m];.X..m|ax1 = m1|
|000007f0| 2d 3e 6d 5f 6d 61 78 5b | 64 69 6d 5d 3b 0a 58 09 |->m_max[|dim];.X.|
|00000800| 09 6d 69 6e 32 20 3d 20 | 6d 32 2d 3e 6d 5f 6d 69 |.min2 = |m2->m_mi|
|00000810| 6e 5b 64 69 6d 5d 3b 0a | 58 09 09 6d 61 78 32 20 |n[dim];.|X..max2 |
|00000820| 3d 20 6d 32 2d 3e 6d 5f | 6d 61 78 5b 64 69 6d 5d |= m2->m_|max[dim]|
|00000830| 3b 0a 58 09 09 69 66 20 | 28 28 6d 69 6e 31 20 26 |;.X..if |((min1 &|
|00000840| 26 20 6d 69 6e 32 20 26 | 26 20 28 6d 69 6e 31 20 |& min2 &|& (min1 |
|00000850| 21 3d 20 6d 69 6e 32 29 | 29 20 7c 7c 20 28 28 6d |!= min2)|) || ((m|
|00000860| 61 78 31 2d 6d 69 6e 31 | 29 20 21 3d 20 28 6d 61 |ax1-min1|) != (ma|
|00000870| 78 32 2d 6d 69 6e 32 29 | 29 29 0a 58 09 09 09 65 |x2-min2)|)).X...e|
|00000880| 72 72 6f 72 28 22 49 6e | 63 6f 6d 70 61 74 69 62 |rror("In|compatib|
|00000890| 6c 65 20 6d 61 74 72 69 | 78 20 62 6f 75 6e 64 73 |le matri|x bounds|
|000008a0| 20 66 6f 72 20 61 64 64 | 22 29 3b 0a 58 09 09 74 | for add|");.X..t|
|000008b0| 6d 70 2e 6d 5f 6d 69 6e | 5b 64 69 6d 5d 20 3d 20 |mp.m_min|[dim] = |
|000008c0| 28 6d 69 6e 31 20 3f 20 | 6d 69 6e 31 20 3a 20 6d |(min1 ? |min1 : m|
|000008d0| 69 6e 32 29 3b 0a 58 09 | 09 74 6d 70 2e 6d 5f 6d |in2);.X.|.tmp.m_m|
|000008e0| 61 78 5b 64 69 6d 5d 20 | 3d 20 74 6d 70 2e 6d 5f |ax[dim] |= tmp.m_|
|000008f0| 6d 69 6e 5b 64 69 6d 5d | 20 2b 20 28 6d 61 78 31 |min[dim]| + (max1|
|00000900| 20 2d 20 6d 69 6e 31 29 | 3b 0a 58 09 7d 0a 58 09 | - min1)|;.X.}.X.|
|00000910| 72 65 73 20 3d 20 6d 61 | 74 61 6c 6c 6f 63 28 6d |res = ma|talloc(m|
|00000920| 31 2d 3e 6d 5f 73 69 7a | 65 29 3b 0a 58 09 2a 72 |1->m_siz|e);.X.*r|
|00000930| 65 73 20 3d 20 74 6d 70 | 3b 0a 58 09 76 31 20 3d |es = tmp|;.X.v1 =|
|00000940| 20 6d 31 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 | m1->m_t|able;.X.|
|00000950| 76 32 20 3d 20 6d 32 2d | 3e 6d 5f 74 61 62 6c 65 |v2 = m2-|>m_table|
|00000960| 3b 0a 58 09 76 72 65 73 | 20 3d 20 72 65 73 2d 3e |;.X.vres| = res->|
|00000970| 6d 5f 74 61 62 6c 65 3b | 0a 58 09 66 6f 72 20 28 |m_table;|.X.for (|
|00000980| 69 6e 64 65 78 20 3d 20 | 6d 31 2d 3e 6d 5f 73 69 |index = |m1->m_si|
|00000990| 7a 65 3b 20 69 6e 64 65 | 78 20 3e 20 30 3b 20 69 |ze; inde|x > 0; i|
|000009a0| 6e 64 65 78 2d 2d 29 0a | 58 09 09 61 64 64 76 61 |ndex--).|X..addva|
|000009b0| 6c 75 65 28 76 31 2b 2b | 2c 20 76 32 2b 2b 2c 20 |lue(v1++|, v2++, |
|000009c0| 76 72 65 73 2b 2b 29 3b | 0a 58 09 72 65 74 75 72 |vres++);|.X.retur|
|000009d0| 6e 20 72 65 73 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |n res;.X|}.X.X.X/|
|000009e0| 2a 0a 58 20 2a 20 53 75 | 62 74 72 61 63 74 20 74 |*.X * Su|btract t|
|000009f0| 77 6f 20 63 6f 6d 70 61 | 74 69 62 6c 65 20 6d 61 |wo compa|tible ma|
|00000a00| 74 72 69 63 65 73 2e 0a | 58 20 2a 2f 0a 58 4d 41 |trices..|X */.XMA|
|00000a10| 54 52 49 58 20 2a 0a 58 | 6d 61 74 73 75 62 28 6d |TRIX *.X|matsub(m|
|00000a20| 31 2c 20 6d 32 29 0a 58 | 09 4d 41 54 52 49 58 20 |1, m2).X|.MATRIX |
|00000a30| 2a 6d 31 2c 20 2a 6d 32 | 3b 0a 58 7b 0a 58 09 69 |*m1, *m2|;.X{.X.i|
|00000a40| 6e 74 20 64 69 6d 3b 0a | 58 09 6c 6f 6e 67 20 6d |nt dim;.|X.long m|
|00000a50| 69 6e 31 2c 20 6d 69 6e | 32 2c 20 6d 61 78 31 2c |in1, min|2, max1,|
|00000a60| 20 6d 61 78 32 2c 20 69 | 6e 64 65 78 3b 0a 58 09 | max2, i|ndex;.X.|
|00000a70| 56 41 4c 55 45 20 2a 76 | 31 2c 20 2a 76 32 2c 20 |VALUE *v|1, *v2, |
|00000a80| 2a 76 72 65 73 3b 0a 58 | 09 4d 41 54 52 49 58 20 |*vres;.X|.MATRIX |
|00000a90| 2a 72 65 73 3b 0a 58 09 | 4d 41 54 52 49 58 20 74 |*res;.X.|MATRIX t|
|00000aa0| 6d 70 3b 0a 58 0a 58 09 | 69 66 20 28 6d 31 2d 3e |mp;.X.X.|if (m1->|
|00000ab0| 6d 5f 64 69 6d 20 21 3d | 20 6d 32 2d 3e 6d 5f 64 |m_dim !=| m2->m_d|
|00000ac0| 69 6d 29 0a 58 09 09 65 | 72 72 6f 72 28 22 49 6e |im).X..e|rror("In|
|00000ad0| 63 6f 6d 70 61 74 69 62 | 6c 65 20 6d 61 74 72 69 |compatib|le matri|
|00000ae0| 78 20 64 69 6d 65 6e 73 | 69 6f 6e 73 20 66 6f 72 |x dimens|ions for|
|00000af0| 20 73 75 62 22 29 3b 0a | 58 09 74 6d 70 2e 6d 5f | sub");.|X.tmp.m_|
|00000b00| 64 69 6d 20 3d 20 6d 31 | 2d 3e 6d 5f 64 69 6d 3b |dim = m1|->m_dim;|
|00000b10| 0a 58 09 74 6d 70 2e 6d | 5f 73 69 7a 65 20 3d 20 |.X.tmp.m|_size = |
|00000b20| 6d 31 2d 3e 6d 5f 73 69 | 7a 65 3b 0a 58 09 66 6f |m1->m_si|ze;.X.fo|
|00000b30| 72 20 28 64 69 6d 20 3d | 20 30 3b 20 64 69 6d 20 |r (dim =| 0; dim |
|00000b40| 3c 20 6d 31 2d 3e 6d 5f | 64 69 6d 3b 20 64 69 6d |< m1->m_|dim; dim|
|00000b50| 2b 2b 29 20 7b 0a 58 09 | 09 6d 69 6e 31 20 3d 20 |++) {.X.|.min1 = |
|00000b60| 6d 31 2d 3e 6d 5f 6d 69 | 6e 5b 64 69 6d 5d 3b 0a |m1->m_mi|n[dim];.|
|00000b70| 58 09 09 6d 61 78 31 20 | 3d 20 6d 31 2d 3e 6d 5f |X..max1 |= m1->m_|
|00000b80| 6d 61 78 5b 64 69 6d 5d | 3b 0a 58 09 09 6d 69 6e |max[dim]|;.X..min|
|00000b90| 32 20 3d 20 6d 32 2d 3e | 6d 5f 6d 69 6e 5b 64 69 |2 = m2->|m_min[di|
|00000ba0| 6d 5d 3b 0a 58 09 09 6d | 61 78 32 20 3d 20 6d 32 |m];.X..m|ax2 = m2|
|00000bb0| 2d 3e 6d 5f 6d 61 78 5b | 64 69 6d 5d 3b 0a 58 09 |->m_max[|dim];.X.|
|00000bc0| 09 69 66 20 28 28 6d 69 | 6e 31 20 26 26 20 6d 69 |.if ((mi|n1 && mi|
|00000bd0| 6e 32 20 26 26 20 28 6d | 69 6e 31 20 21 3d 20 6d |n2 && (m|in1 != m|
|00000be0| 69 6e 32 29 29 20 7c 7c | 20 28 28 6d 61 78 31 2d |in2)) ||| ((max1-|
|00000bf0| 6d 69 6e 31 29 20 21 3d | 20 28 6d 61 78 32 2d 6d |min1) !=| (max2-m|
|00000c00| 69 6e 32 29 29 29 0a 58 | 09 09 09 65 72 72 6f 72 |in2))).X|...error|
|00000c10| 28 22 49 6e 63 6f 6d 70 | 61 74 69 62 6c 65 20 6d |("Incomp|atible m|
|00000c20| 61 74 72 69 78 20 62 6f | 75 6e 64 73 20 66 6f 72 |atrix bo|unds for|
|00000c30| 20 73 75 62 22 29 3b 0a | 58 09 09 74 6d 70 2e 6d | sub");.|X..tmp.m|
|00000c40| 5f 6d 69 6e 5b 64 69 6d | 5d 20 3d 20 28 6d 69 6e |_min[dim|] = (min|
|00000c50| 31 20 3f 20 6d 69 6e 31 | 20 3a 20 6d 69 6e 32 29 |1 ? min1| : min2)|
|00000c60| 3b 0a 58 09 09 74 6d 70 | 2e 6d 5f 6d 61 78 5b 64 |;.X..tmp|.m_max[d|
|00000c70| 69 6d 5d 20 3d 20 74 6d | 70 2e 6d 5f 6d 69 6e 5b |im] = tm|p.m_min[|
|00000c80| 64 69 6d 5d 20 2b 20 28 | 6d 61 78 31 20 2d 20 6d |dim] + (|max1 - m|
|00000c90| 69 6e 31 29 3b 0a 58 09 | 7d 0a 58 09 72 65 73 20 |in1);.X.|}.X.res |
|00000ca0| 3d 20 6d 61 74 61 6c 6c | 6f 63 28 6d 31 2d 3e 6d |= matall|oc(m1->m|
|00000cb0| 5f 73 69 7a 65 29 3b 0a | 58 09 2a 72 65 73 20 3d |_size);.|X.*res =|
|00000cc0| 20 74 6d 70 3b 0a 58 09 | 76 31 20 3d 20 6d 31 2d | tmp;.X.|v1 = m1-|
|00000cd0| 3e 6d 5f 74 61 62 6c 65 | 3b 0a 58 09 76 32 20 3d |>m_table|;.X.v2 =|
|00000ce0| 20 6d 32 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 | m2->m_t|able;.X.|
|00000cf0| 76 72 65 73 20 3d 20 72 | 65 73 2d 3e 6d 5f 74 61 |vres = r|es->m_ta|
|00000d00| 62 6c 65 3b 0a 58 09 66 | 6f 72 20 28 69 6e 64 65 |ble;.X.f|or (inde|
|00000d10| 78 20 3d 20 6d 31 2d 3e | 6d 5f 73 69 7a 65 3b 20 |x = m1->|m_size; |
|00000d20| 69 6e 64 65 78 20 3e 20 | 30 3b 20 69 6e 64 65 78 |index > |0; index|
|00000d30| 2d 2d 29 0a 58 09 09 73 | 75 62 76 61 6c 75 65 28 |--).X..s|ubvalue(|
|00000d40| 76 31 2b 2b 2c 20 76 32 | 2b 2b 2c 20 76 72 65 73 |v1++, v2|++, vres|
|00000d50| 2b 2b 29 3b 0a 58 09 72 | 65 74 75 72 6e 20 72 65 |++);.X.r|eturn re|
|00000d60| 73 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |s;.X}.X.|X.X/*.X |
|00000d70| 2a 20 50 72 6f 64 75 63 | 65 20 74 68 65 20 6e 65 |* Produc|e the ne|
|00000d80| 67 61 74 69 76 65 20 6f | 66 20 61 20 6d 61 74 72 |gative o|f a matr|
|00000d90| 69 78 2e 0a 58 20 2a 2f | 0a 58 4d 41 54 52 49 58 |ix..X */|.XMATRIX|
|00000da0| 20 2a 0a 58 6d 61 74 6e | 65 67 28 6d 29 0a 58 09 | *.Xmatn|eg(m).X.|
|00000db0| 4d 41 54 52 49 58 20 2a | 6d 3b 0a 58 7b 0a 58 09 |MATRIX *|m;.X{.X.|
|00000dc0| 72 65 67 69 73 74 65 72 | 20 56 41 4c 55 45 20 2a |register| VALUE *|
|00000dd0| 76 61 6c 2c 20 2a 76 72 | 65 73 3b 0a 58 09 6c 6f |val, *vr|es;.X.lo|
|00000de0| 6e 67 20 69 6e 64 65 78 | 3b 0a 58 09 4d 41 54 52 |ng index|;.X.MATR|
|00000df0| 49 58 20 2a 72 65 73 3b | 0a 58 0a 58 09 72 65 73 |IX *res;|.X.X.res|
|00000e00| 20 3d 20 6d 61 74 61 6c | 6c 6f 63 28 6d 2d 3e 6d | = matal|loc(m->m|
|00000e10| 5f 73 69 7a 65 29 3b 0a | 58 09 2a 72 65 73 20 3d |_size);.|X.*res =|
|00000e20| 20 2a 6d 3b 0a 58 09 76 | 61 6c 20 3d 20 6d 2d 3e | *m;.X.v|al = m->|
|00000e30| 6d 5f 74 61 62 6c 65 3b | 0a 58 09 76 72 65 73 20 |m_table;|.X.vres |
|00000e40| 3d 20 72 65 73 2d 3e 6d | 5f 74 61 62 6c 65 3b 0a |= res->m|_table;.|
|00000e50| 58 09 66 6f 72 20 28 69 | 6e 64 65 78 20 3d 20 6d |X.for (i|ndex = m|
|00000e60| 2d 3e 6d 5f 73 69 7a 65 | 3b 20 69 6e 64 65 78 20 |->m_size|; index |
|00000e70| 3e 20 30 3b 20 69 6e 64 | 65 78 2d 2d 29 0a 58 09 |> 0; ind|ex--).X.|
|00000e80| 09 6e 65 67 76 61 6c 75 | 65 28 76 61 6c 2b 2b 2c |.negvalu|e(val++,|
|00000e90| 20 76 72 65 73 2b 2b 29 | 3b 0a 58 09 72 65 74 75 | vres++)|;.X.retu|
|00000ea0| 72 6e 20 72 65 73 3b 0a | 58 7d 0a 58 0a 58 0a 58 |rn res;.|X}.X.X.X|
|00000eb0| 2f 2a 0a 58 20 2a 20 4d | 75 6c 74 69 70 6c 79 20 |/*.X * M|ultiply |
|00000ec0| 74 77 6f 20 63 6f 6d 70 | 61 74 69 62 6c 65 20 6d |two comp|atible m|
|00000ed0| 61 74 72 69 63 65 73 2e | 0a 58 20 2a 2f 0a 58 4d |atrices.|.X */.XM|
|00000ee0| 41 54 52 49 58 20 2a 0a | 58 6d 61 74 6d 75 6c 28 |ATRIX *.|Xmatmul(|
|00000ef0| 6d 31 2c 20 6d 32 29 0a | 58 09 4d 41 54 52 49 58 |m1, m2).|X.MATRIX|
|00000f00| 20 2a 6d 31 2c 20 2a 6d | 32 3b 0a 58 7b 0a 58 09 | *m1, *m|2;.X{.X.|
|00000f10| 72 65 67 69 73 74 65 72 | 20 4d 41 54 52 49 58 20 |register| MATRIX |
|00000f20| 2a 72 65 73 3b 0a 58 09 | 6c 6f 6e 67 20 69 31 2c |*res;.X.|long i1,|
|00000f30| 20 69 32 2c 20 6d 61 78 | 31 2c 20 6d 61 78 32 2c | i2, max|1, max2,|
|00000f40| 20 69 6e 64 65 78 2c 20 | 6d 61 78 69 6e 64 65 78 | index, |maxindex|
|00000f50| 3b 0a 58 09 56 41 4c 55 | 45 20 2a 76 31 2c 20 2a |;.X.VALU|E *v1, *|
|00000f60| 76 32 3b 0a 58 09 56 41 | 4c 55 45 20 73 75 6d 2c |v2;.X.VA|LUE sum,|
|00000f70| 20 74 6d 70 31 2c 20 74 | 6d 70 32 3b 0a 58 0a 58 | tmp1, t|mp2;.X.X|
|00000f80| 09 69 66 20 28 28 6d 31 | 2d 3e 6d 5f 64 69 6d 20 |.if ((m1|->m_dim |
|00000f90| 21 3d 20 32 29 20 7c 7c | 20 28 6d 32 2d 3e 6d 5f |!= 2) ||| (m2->m_|
|00000fa0| 64 69 6d 20 21 3d 20 32 | 29 29 0a 58 09 09 65 72 |dim != 2|)).X..er|
|00000fb0| 72 6f 72 28 22 4d 61 74 | 72 69 78 20 64 69 6d 65 |ror("Mat|rix dime|
|00000fc0| 6e 73 69 6f 6e 20 6d 75 | 73 74 20 62 65 20 74 77 |nsion mu|st be tw|
|00000fd0| 6f 20 66 6f 72 20 6d 75 | 6c 22 29 3b 0a 58 09 69 |o for mu|l");.X.i|
|00000fe0| 66 20 28 28 6d 31 2d 3e | 6d 5f 6d 61 78 5b 31 5d |f ((m1->|m_max[1]|
|00000ff0| 20 2d 20 6d 31 2d 3e 6d | 5f 6d 69 6e 5b 31 5d 29 | - m1->m|_min[1])|
|00001000| 20 21 3d 20 28 6d 32 2d | 3e 6d 5f 6d 61 78 5b 30 | != (m2-|>m_max[0|
|00001010| 5d 20 2d 20 6d 32 2d 3e | 6d 5f 6d 69 6e 5b 30 5d |] - m2->|m_min[0]|
|00001020| 29 29 0a 58 09 09 65 72 | 72 6f 72 28 22 49 6e 63 |)).X..er|ror("Inc|
|00001030| 6f 6d 70 61 74 69 62 6c | 65 20 62 6f 75 6e 64 73 |ompatibl|e bounds|
|00001040| 20 66 6f 72 20 6d 61 74 | 72 69 78 20 6d 75 6c 22 | for mat|rix mul"|
|00001050| 29 3b 0a 58 09 6d 61 78 | 31 20 3d 20 28 6d 31 2d |);.X.max|1 = (m1-|
|00001060| 3e 6d 5f 6d 61 78 5b 30 | 5d 20 2d 20 6d 31 2d 3e |>m_max[0|] - m1->|
|00001070| 6d 5f 6d 69 6e 5b 30 5d | 20 2b 20 31 29 3b 0a 58 |m_min[0]| + 1);.X|
|00001080| 09 6d 61 78 32 20 3d 20 | 28 6d 32 2d 3e 6d 5f 6d |.max2 = |(m2->m_m|
|00001090| 61 78 5b 31 5d 20 2d 20 | 6d 32 2d 3e 6d 5f 6d 69 |ax[1] - |m2->m_mi|
|000010a0| 6e 5b 31 5d 20 2b 20 31 | 29 3b 0a 58 09 6d 61 78 |n[1] + 1|);.X.max|
|000010b0| 69 6e 64 65 78 20 3d 20 | 28 6d 31 2d 3e 6d 5f 6d |index = |(m1->m_m|
|000010c0| 61 78 5b 31 5d 20 2d 20 | 6d 31 2d 3e 6d 5f 6d 69 |ax[1] - |m1->m_mi|
|000010d0| 6e 5b 31 5d 20 2b 20 31 | 29 3b 0a 58 09 72 65 73 |n[1] + 1|);.X.res|
|000010e0| 20 3d 20 6d 61 74 61 6c | 6c 6f 63 28 6d 61 78 31 | = matal|loc(max1|
|000010f0| 20 2a 20 6d 61 78 32 29 | 3b 0a 58 09 72 65 73 2d | * max2)|;.X.res-|
|00001100| 3e 6d 5f 64 69 6d 20 3d | 20 32 3b 0a 58 09 72 65 |>m_dim =| 2;.X.re|
|00001110| 73 2d 3e 6d 5f 6d 69 6e | 5b 30 5d 20 3d 20 6d 31 |s->m_min|[0] = m1|
|00001120| 2d 3e 6d 5f 6d 69 6e 5b | 30 5d 3b 0a 58 09 72 65 |->m_min[|0];.X.re|
|00001130| 73 2d 3e 6d 5f 6d 61 78 | 5b 30 5d 20 3d 20 6d 31 |s->m_max|[0] = m1|
|00001140| 2d 3e 6d 5f 6d 61 78 5b | 30 5d 3b 0a 58 09 72 65 |->m_max[|0];.X.re|
|00001150| 73 2d 3e 6d 5f 6d 69 6e | 5b 31 5d 20 3d 20 6d 32 |s->m_min|[1] = m2|
|00001160| 2d 3e 6d 5f 6d 69 6e 5b | 31 5d 3b 0a 58 09 72 65 |->m_min[|1];.X.re|
|00001170| 73 2d 3e 6d 5f 6d 61 78 | 5b 31 5d 20 3d 20 6d 32 |s->m_max|[1] = m2|
|00001180| 2d 3e 6d 5f 6d 61 78 5b | 31 5d 3b 0a 58 09 66 6f |->m_max[|1];.X.fo|
|00001190| 72 20 28 69 31 20 3d 20 | 30 3b 20 69 31 20 3c 20 |r (i1 = |0; i1 < |
|000011a0| 6d 61 78 31 3b 20 69 31 | 2b 2b 29 20 7b 0a 58 09 |max1; i1|++) {.X.|
|000011b0| 09 66 6f 72 20 28 69 32 | 20 3d 20 30 3b 20 69 32 |.for (i2| = 0; i2|
|000011c0| 20 3c 20 6d 61 78 32 3b | 20 69 32 2b 2b 29 20 7b | < max2;| i2++) {|
|000011d0| 0a 58 09 09 09 73 75 6d | 2e 76 5f 6e 75 6d 20 3d |.X...sum|.v_num =|
|000011e0| 20 71 6c 69 6e 6b 28 26 | 5f 71 7a 65 72 6f 5f 29 | qlink(&|_qzero_)|
|000011f0| 3b 0a 58 09 09 09 73 75 | 6d 2e 76 5f 74 79 70 65 |;.X...su|m.v_type|
|00001200| 20 3d 20 56 5f 4e 55 4d | 3b 0a 58 09 09 09 76 31 | = V_NUM|;.X...v1|
|00001210| 20 3d 20 26 6d 31 2d 3e | 6d 5f 74 61 62 6c 65 5b | = &m1->|m_table[|
|00001220| 69 31 20 2a 20 6d 61 78 | 69 6e 64 65 78 5d 3b 0a |i1 * max|index];.|
|00001230| 58 09 09 09 76 32 20 3d | 20 26 6d 32 2d 3e 6d 5f |X...v2 =| &m2->m_|
|00001240| 74 61 62 6c 65 5b 69 32 | 5d 3b 0a 58 09 09 09 66 |table[i2|];.X...f|
|00001250| 6f 72 20 28 69 6e 64 65 | 78 20 3d 20 30 3b 20 69 |or (inde|x = 0; i|
|00001260| 6e 64 65 78 20 3c 20 6d | 61 78 69 6e 64 65 78 3b |ndex < m|axindex;|
|00001270| 20 69 6e 64 65 78 2b 2b | 29 20 7b 0a 58 09 09 09 | index++|) {.X...|
|00001280| 09 6d 75 6c 76 61 6c 75 | 65 28 76 31 2c 20 76 32 |.mulvalu|e(v1, v2|
|00001290| 2c 20 26 74 6d 70 31 29 | 3b 0a 58 09 09 09 09 61 |, &tmp1)|;.X....a|
|000012a0| 64 64 76 61 6c 75 65 28 | 26 73 75 6d 2c 20 26 74 |ddvalue(|&sum, &t|
|000012b0| 6d 70 31 2c 20 26 74 6d | 70 32 29 3b 0a 58 09 09 |mp1, &tm|p2);.X..|
|000012c0| 09 09 66 72 65 65 76 61 | 6c 75 65 28 26 74 6d 70 |..freeva|lue(&tmp|
|000012d0| 31 29 3b 0a 58 09 09 09 | 09 66 72 65 65 76 61 6c |1);.X...|.freeval|
|000012e0| 75 65 28 26 73 75 6d 29 | 3b 0a 58 09 09 09 09 73 |ue(&sum)|;.X....s|
|000012f0| 75 6d 20 3d 20 74 6d 70 | 32 3b 0a 58 09 09 09 09 |um = tmp|2;.X....|
|00001300| 76 31 2b 2b 3b 0a 58 09 | 09 09 09 76 32 20 2b 3d |v1++;.X.|...v2 +=|
|00001310| 20 6d 61 78 32 3b 0a 58 | 09 09 09 7d 0a 58 09 09 | max2;.X|...}.X..|
|00001320| 09 69 6e 64 65 78 20 3d | 20 28 69 31 20 2a 20 6d |.index =| (i1 * m|
|00001330| 61 78 32 29 20 2b 20 69 | 32 3b 0a 58 09 09 09 72 |ax2) + i|2;.X...r|
|00001340| 65 73 2d 3e 6d 5f 74 61 | 62 6c 65 5b 69 6e 64 65 |es->m_ta|ble[inde|
|00001350| 78 5d 20 3d 20 73 75 6d | 3b 0a 58 09 09 7d 0a 58 |x] = sum|;.X..}.X|
|00001360| 09 7d 0a 58 09 72 65 74 | 75 72 6e 20 72 65 73 3b |.}.X.ret|urn res;|
|00001370| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |.X}.X.X.|X/*.X * |
|00001380| 53 71 75 61 72 65 20 61 | 20 6d 61 74 72 69 78 2e |Square a| matrix.|
|00001390| 0a 58 20 2a 2f 0a 58 4d | 41 54 52 49 58 20 2a 0a |.X */.XM|ATRIX *.|
|000013a0| 58 6d 61 74 73 71 75 61 | 72 65 28 6d 29 0a 58 09 |Xmatsqua|re(m).X.|
|000013b0| 4d 41 54 52 49 58 20 2a | 6d 3b 0a 58 7b 0a 58 09 |MATRIX *|m;.X{.X.|
|000013c0| 72 65 67 69 73 74 65 72 | 20 4d 41 54 52 49 58 20 |register| MATRIX |
|000013d0| 2a 72 65 73 3b 0a 58 09 | 6c 6f 6e 67 20 69 31 2c |*res;.X.|long i1,|
|000013e0| 20 69 32 2c 20 6d 61 78 | 2c 20 69 6e 64 65 78 3b | i2, max|, index;|
|000013f0| 0a 58 09 56 41 4c 55 45 | 20 2a 76 31 2c 20 2a 76 |.X.VALUE| *v1, *v|
|00001400| 32 3b 0a 58 09 56 41 4c | 55 45 20 73 75 6d 2c 20 |2;.X.VAL|UE sum, |
|00001410| 74 6d 70 31 2c 20 74 6d | 70 32 3b 0a 58 0a 58 09 |tmp1, tm|p2;.X.X.|
|00001420| 69 66 20 28 6d 2d 3e 6d | 5f 64 69 6d 20 21 3d 20 |if (m->m|_dim != |
|00001430| 32 29 0a 58 09 09 65 72 | 72 6f 72 28 22 4d 61 74 |2).X..er|ror("Mat|
|00001440| 72 69 78 20 64 69 6d 65 | 6e 73 69 6f 6e 20 6d 75 |rix dime|nsion mu|
|00001450| 73 74 20 62 65 20 74 77 | 6f 20 66 6f 72 20 73 71 |st be tw|o for sq|
|00001460| 75 61 72 65 22 29 3b 0a | 58 09 69 66 20 28 28 6d |uare");.|X.if ((m|
|00001470| 2d 3e 6d 5f 6d 61 78 5b | 30 5d 20 2d 20 6d 2d 3e |->m_max[|0] - m->|
|00001480| 6d 5f 6d 69 6e 5b 30 5d | 29 20 21 3d 20 28 6d 2d |m_min[0]|) != (m-|
|00001490| 3e 6d 5f 6d 61 78 5b 31 | 5d 20 2d 20 6d 2d 3e 6d |>m_max[1|] - m->m|
|000014a0| 5f 6d 69 6e 5b 31 5d 29 | 29 0a 58 09 09 65 72 72 |_min[1])|).X..err|
|000014b0| 6f 72 28 22 53 71 75 61 | 72 69 6e 67 20 6e 6f 6e |or("Squa|ring non|
|000014c0| 2d 73 71 75 61 72 65 20 | 6d 61 74 72 69 78 22 29 |-square |matrix")|
|000014d0| 3b 0a 58 09 6d 61 78 20 | 3d 20 28 6d 2d 3e 6d 5f |;.X.max |= (m->m_|
|000014e0| 6d 61 78 5b 30 5d 20 2d | 20 6d 2d 3e 6d 5f 6d 69 |max[0] -| m->m_mi|
|000014f0| 6e 5b 30 5d 20 2b 20 31 | 29 3b 0a 58 09 72 65 73 |n[0] + 1|);.X.res|
|00001500| 20 3d 20 6d 61 74 61 6c | 6c 6f 63 28 6d 61 78 20 | = matal|loc(max |
|00001510| 2a 20 6d 61 78 29 3b 0a | 58 09 72 65 73 2d 3e 6d |* max);.|X.res->m|
|00001520| 5f 64 69 6d 20 3d 20 32 | 3b 0a 58 09 72 65 73 2d |_dim = 2|;.X.res-|
|00001530| 3e 6d 5f 6d 69 6e 5b 30 | 5d 20 3d 20 6d 2d 3e 6d |>m_min[0|] = m->m|
|00001540| 5f 6d 69 6e 5b 30 5d 3b | 0a 58 09 72 65 73 2d 3e |_min[0];|.X.res->|
|00001550| 6d 5f 6d 61 78 5b 30 5d | 20 3d 20 6d 2d 3e 6d 5f |m_max[0]| = m->m_|
|00001560| 6d 61 78 5b 30 5d 3b 0a | 58 09 72 65 73 2d 3e 6d |max[0];.|X.res->m|
|00001570| 5f 6d 69 6e 5b 31 5d 20 | 3d 20 6d 2d 3e 6d 5f 6d |_min[1] |= m->m_m|
|00001580| 69 6e 5b 31 5d 3b 0a 58 | 09 72 65 73 2d 3e 6d 5f |in[1];.X|.res->m_|
|00001590| 6d 61 78 5b 31 5d 20 3d | 20 6d 2d 3e 6d 5f 6d 61 |max[1] =| m->m_ma|
|000015a0| 78 5b 31 5d 3b 0a 58 09 | 66 6f 72 20 28 69 31 20 |x[1];.X.|for (i1 |
|000015b0| 3d 20 30 3b 20 69 31 20 | 3c 20 6d 61 78 3b 20 69 |= 0; i1 |< max; i|
|000015c0| 31 2b 2b 29 20 7b 0a 58 | 09 09 66 6f 72 20 28 69 |1++) {.X|..for (i|
|000015d0| 32 20 3d 20 30 3b 20 69 | 32 20 3c 20 6d 61 78 3b |2 = 0; i|2 < max;|
|000015e0| 20 69 32 2b 2b 29 20 7b | 0a 58 09 09 09 73 75 6d | i2++) {|.X...sum|
|000015f0| 2e 76 5f 6e 75 6d 20 3d | 20 71 6c 69 6e 6b 28 26 |.v_num =| qlink(&|
|00001600| 5f 71 7a 65 72 6f 5f 29 | 3b 0a 58 09 09 09 73 75 |_qzero_)|;.X...su|
|00001610| 6d 2e 76 5f 74 79 70 65 | 20 3d 20 56 5f 4e 55 4d |m.v_type| = V_NUM|
|00001620| 3b 0a 58 09 09 09 76 31 | 20 3d 20 26 6d 2d 3e 6d |;.X...v1| = &m->m|
|00001630| 5f 74 61 62 6c 65 5b 69 | 31 20 2a 20 6d 61 78 5d |_table[i|1 * max]|
|00001640| 3b 0a 58 09 09 09 76 32 | 20 3d 20 26 6d 2d 3e 6d |;.X...v2| = &m->m|
|00001650| 5f 74 61 62 6c 65 5b 69 | 32 5d 3b 0a 58 09 09 09 |_table[i|2];.X...|
|00001660| 66 6f 72 20 28 69 6e 64 | 65 78 20 3d 20 30 3b 20 |for (ind|ex = 0; |
|00001670| 69 6e 64 65 78 20 3c 20 | 6d 61 78 3b 20 69 6e 64 |index < |max; ind|
|00001680| 65 78 2b 2b 29 20 7b 0a | 58 09 09 09 09 6d 75 6c |ex++) {.|X....mul|
|00001690| 76 61 6c 75 65 28 76 31 | 2c 20 76 32 2c 20 26 74 |value(v1|, v2, &t|
|000016a0| 6d 70 31 29 3b 0a 58 09 | 09 09 09 61 64 64 76 61 |mp1);.X.|...addva|
|000016b0| 6c 75 65 28 26 73 75 6d | 2c 20 26 74 6d 70 31 2c |lue(&sum|, &tmp1,|
|000016c0| 20 26 74 6d 70 32 29 3b | 0a 58 09 09 09 09 66 72 | &tmp2);|.X....fr|
|000016d0| 65 65 76 61 6c 75 65 28 | 26 74 6d 70 31 29 3b 0a |eevalue(|&tmp1);.|
|000016e0| 58 09 09 09 09 66 72 65 | 65 76 61 6c 75 65 28 26 |X....fre|evalue(&|
|000016f0| 73 75 6d 29 3b 0a 58 09 | 09 09 09 73 75 6d 20 3d |sum);.X.|...sum =|
|00001700| 20 74 6d 70 32 3b 0a 58 | 09 09 09 09 76 31 2b 2b | tmp2;.X|....v1++|
|00001710| 3b 0a 58 09 09 09 09 76 | 32 20 2b 3d 20 6d 61 78 |;.X....v|2 += max|
|00001720| 3b 0a 58 09 09 09 7d 0a | 58 09 09 09 69 6e 64 65 |;.X...}.|X...inde|
|00001730| 78 20 3d 20 28 69 31 20 | 2a 20 6d 61 78 29 20 2b |x = (i1 |* max) +|
|00001740| 20 69 32 3b 0a 58 09 09 | 09 72 65 73 2d 3e 6d 5f | i2;.X..|.res->m_|
|00001750| 74 61 62 6c 65 5b 69 6e | 64 65 78 5d 20 3d 20 73 |table[in|dex] = s|
|00001760| 75 6d 3b 0a 58 09 09 7d | 0a 58 09 7d 0a 58 09 72 |um;.X..}|.X.}.X.r|
|00001770| 65 74 75 72 6e 20 72 65 | 73 3b 0a 58 7d 0a 58 0a |eturn re|s;.X}.X.|
|00001780| 58 0a 58 2f 2a 0a 58 20 | 2a 20 43 6f 6d 70 75 74 |X.X/*.X |* Comput|
|00001790| 65 20 74 68 65 20 72 65 | 73 75 6c 74 20 6f 66 20 |e the re|sult of |
|000017a0| 72 61 69 73 69 6e 67 20 | 61 20 73 71 75 61 72 65 |raising |a square|
|000017b0| 20 6d 61 74 72 69 78 20 | 74 6f 20 61 6e 20 69 6e | matrix |to an in|
|000017c0| 74 65 67 65 72 20 70 6f | 77 65 72 2e 0a 58 20 2a |teger po|wer..X *|
|000017d0| 20 4e 65 67 61 74 69 76 | 65 20 70 6f 77 65 72 73 | Negativ|e powers|
|000017e0| 20 6d 65 61 6e 20 74 68 | 65 20 70 6f 73 69 74 69 | mean th|e positi|
|000017f0| 76 65 20 70 6f 77 65 72 | 20 6f 66 20 74 68 65 20 |ve power| of the |
|00001800| 69 6e 76 65 72 73 65 2e | 0a 58 20 2a 20 4e 6f 74 |inverse.|.X * Not|
|00001810| 65 3a 20 54 68 69 73 20 | 63 61 6c 63 75 6c 61 74 |e: This |calculat|
|00001820| 69 6f 6e 20 63 6f 75 6c | 64 20 73 6f 6d 65 64 61 |ion coul|d someda|
|00001830| 79 20 62 65 20 69 6d 70 | 72 6f 76 65 64 20 66 6f |y be imp|roved fo|
|00001840| 72 20 6c 61 72 67 65 20 | 70 6f 77 65 72 73 0a 58 |r large |powers.X|
|00001850| 20 2a 20 62 79 20 75 73 | 69 6e 67 20 74 68 65 20 | * by us|ing the |
|00001860| 63 68 61 72 61 63 74 65 | 72 69 73 74 69 63 20 70 |characte|ristic p|
|00001870| 6f 6c 79 6e 6f 6d 69 61 | 6c 20 6f 66 20 74 68 65 |olynomia|l of the|
|00001880| 20 6d 61 74 72 69 78 2e | 0a 58 20 2a 2f 0a 58 4d | matrix.|.X */.XM|
|00001890| 41 54 52 49 58 20 2a 0a | 58 6d 61 74 70 6f 77 69 |ATRIX *.|Xmatpowi|
|000018a0| 28 6d 2c 20 71 29 0a 58 | 09 4d 41 54 52 49 58 20 |(m, q).X|.MATRIX |
|000018b0| 2a 6d 3b 09 09 2f 2a 20 | 6d 61 74 72 69 78 20 74 |*m;../* |matrix t|
|000018c0| 6f 20 62 65 20 72 61 69 | 73 65 64 20 2a 2f 0a 58 |o be rai|sed */.X|
|000018d0| 09 4e 55 4d 42 45 52 20 | 2a 71 3b 09 09 2f 2a 20 |.NUMBER |*q;../* |
|000018e0| 70 6f 77 65 72 20 74 6f | 20 72 61 69 73 65 20 69 |power to| raise i|
|000018f0| 74 20 74 6f 20 2a 2f 0a | 58 7b 0a 58 09 4d 41 54 |t to */.|X{.X.MAT|
|00001900| 52 49 58 20 2a 72 65 73 | 2c 20 2a 74 6d 70 3b 0a |RIX *res|, *tmp;.|
|00001910| 58 09 6c 6f 6e 67 20 70 | 6f 77 65 72 3b 09 09 2f |X.long p|ower;../|
|00001920| 2a 20 70 6f 77 65 72 20 | 74 6f 20 72 61 69 73 65 |* power |to raise|
|00001930| 20 74 6f 20 2a 2f 0a 58 | 09 75 6e 73 69 67 6e 65 | to */.X|.unsigne|
|00001940| 64 20 6c 6f 6e 67 20 62 | 69 74 3b 09 2f 2a 20 63 |d long b|it;./* c|
|00001950| 75 72 72 65 6e 74 20 62 | 69 74 20 76 61 6c 75 65 |urrent b|it value|
|00001960| 20 2a 2f 0a 58 0a 58 09 | 69 66 20 28 6d 2d 3e 6d | */.X.X.|if (m->m|
|00001970| 5f 64 69 6d 20 21 3d 20 | 32 29 0a 58 09 09 65 72 |_dim != |2).X..er|
|00001980| 72 6f 72 28 22 4d 61 74 | 72 69 78 20 64 69 6d 65 |ror("Mat|rix dime|
|00001990| 6e 73 69 6f 6e 20 6d 75 | 73 74 20 62 65 20 74 77 |nsion mu|st be tw|
|000019a0| 6f 20 66 6f 72 20 70 6f | 77 65 72 22 29 3b 0a 58 |o for po|wer");.X|
|000019b0| 09 69 66 20 28 28 6d 2d | 3e 6d 5f 6d 61 78 5b 30 |.if ((m-|>m_max[0|
|000019c0| 5d 20 2d 20 6d 2d 3e 6d | 5f 6d 69 6e 5b 30 5d 29 |] - m->m|_min[0])|
|000019d0| 20 21 3d 20 28 6d 2d 3e | 6d 5f 6d 61 78 5b 31 5d | != (m->|m_max[1]|
|000019e0| 20 2d 20 6d 2d 3e 6d 5f | 6d 69 6e 5b 31 5d 29 29 | - m->m_|min[1]))|
|000019f0| 0a 58 09 09 65 72 72 6f | 72 28 22 52 61 69 73 69 |.X..erro|r("Raisi|
|00001a00| 6e 67 20 6e 6f 6e 2d 73 | 71 75 61 72 65 20 6d 61 |ng non-s|quare ma|
|00001a10| 74 72 69 78 20 74 6f 20 | 61 20 70 6f 77 65 72 22 |trix to |a power"|
|00001a20| 29 3b 0a 58 09 69 66 20 | 28 71 69 73 66 72 61 63 |);.X.if |(qisfrac|
|00001a30| 28 71 29 29 0a 58 09 09 | 65 72 72 6f 72 28 22 52 |(q)).X..|error("R|
|00001a40| 61 69 73 69 6e 67 20 6d | 61 74 72 69 78 20 74 6f |aising m|atrix to|
|00001a50| 20 6e 6f 6e 2d 69 6e 74 | 65 67 72 61 6c 20 70 6f | non-int|egral po|
|00001a60| 77 65 72 22 29 3b 0a 58 | 09 69 66 20 28 69 73 62 |wer");.X|.if (isb|
|00001a70| 69 67 28 71 2d 3e 6e 75 | 6d 29 29 0a 58 09 09 65 |ig(q->nu|m)).X..e|
|00001a80| 72 72 6f 72 28 22 52 61 | 69 73 69 6e 67 20 6d 61 |rror("Ra|ising ma|
|00001a90| 74 72 69 78 20 74 6f 20 | 76 65 72 79 20 6c 61 72 |trix to |very lar|
|00001aa0| 67 65 20 70 6f 77 65 72 | 22 29 3b 0a 58 09 70 6f |ge power|");.X.po|
|00001ab0| 77 65 72 20 3d 20 28 69 | 73 74 69 6e 79 28 71 2d |wer = (i|stiny(q-|
|00001ac0| 3e 6e 75 6d 29 20 3f 20 | 7a 31 74 6f 6c 28 71 2d |>num) ? |z1tol(q-|
|00001ad0| 3e 6e 75 6d 29 20 3a 20 | 7a 32 74 6f 6c 28 71 2d |>num) : |z2tol(q-|
|00001ae0| 3e 6e 75 6d 29 29 3b 0a | 58 09 69 66 20 28 71 69 |>num));.|X.if (qi|
|00001af0| 73 6e 65 67 28 71 29 29 | 0a 58 09 09 70 6f 77 65 |sneg(q))|.X..powe|
|00001b00| 72 20 3d 20 2d 70 6f 77 | 65 72 3b 0a 58 09 2f 2a |r = -pow|er;.X./*|
|00001b10| 0a 58 09 20 2a 20 48 61 | 6e 64 6c 65 20 73 6f 6d |.X. * Ha|ndle som|
|00001b20| 65 20 6c 6f 77 20 70 6f | 77 65 72 73 20 73 70 65 |e low po|wers spe|
|00001b30| 63 69 61 6c 6c 79 0a 58 | 09 20 2a 2f 0a 58 09 69 |cially.X|. */.X.i|
|00001b40| 66 20 28 28 70 6f 77 65 | 72 20 3c 3d 20 34 29 20 |f ((powe|r <= 4) |
|00001b50| 26 26 20 28 70 6f 77 65 | 72 20 3e 3d 20 2d 32 29 |&& (powe|r >= -2)|
|00001b60| 29 20 7b 0a 58 09 09 73 | 77 69 74 63 68 20 28 28 |) {.X..s|witch ((|
|00001b70| 69 6e 74 29 20 70 6f 77 | 65 72 29 20 7b 0a 58 09 |int) pow|er) {.X.|
|00001b80| 09 09 63 61 73 65 20 30 | 3a 0a 58 09 09 09 09 72 |..case 0|:.X....r|
|00001b90| 65 74 75 72 6e 20 6d 61 | 74 69 64 65 6e 74 28 6d |eturn ma|tident(m|
|00001ba0| 29 3b 0a 58 09 09 09 63 | 61 73 65 20 31 3a 0a 58 |);.X...c|ase 1:.X|
|00001bb0| 09 09 09 09 72 65 74 75 | 72 6e 20 6d 61 74 63 6f |....retu|rn matco|
|00001bc0| 70 79 28 6d 29 3b 0a 58 | 09 09 09 63 61 73 65 20 |py(m);.X|...case |
|00001bd0| 2d 31 3a 0a 58 09 09 09 | 09 72 65 74 75 72 6e 20 |-1:.X...|.return |
|00001be0| 6d 61 74 69 6e 76 28 6d | 29 3b 0a 58 09 09 09 63 |matinv(m|);.X...c|
|00001bf0| 61 73 65 20 32 3a 0a 58 | 09 09 09 09 72 65 74 75 |ase 2:.X|....retu|
|00001c00| 72 6e 20 6d 61 74 73 71 | 75 61 72 65 28 6d 29 3b |rn matsq|uare(m);|
|00001c10| 0a 58 09 09 09 63 61 73 | 65 20 2d 32 3a 0a 58 09 |.X...cas|e -2:.X.|
|00001c20| 09 09 09 74 6d 70 20 3d | 20 6d 61 74 69 6e 76 28 |...tmp =| matinv(|
|00001c30| 6d 29 3b 0a 58 09 09 09 | 09 72 65 73 20 3d 20 6d |m);.X...|.res = m|
|00001c40| 61 74 73 71 75 61 72 65 | 28 74 6d 70 29 3b 0a 58 |atsquare|(tmp);.X|
|00001c50| 09 09 09 09 6d 61 74 66 | 72 65 65 28 74 6d 70 29 |....matf|ree(tmp)|
|00001c60| 3b 0a 58 09 09 09 09 72 | 65 74 75 72 6e 20 72 65 |;.X....r|eturn re|
|00001c70| 73 3b 0a 58 09 09 09 63 | 61 73 65 20 33 3a 0a 58 |s;.X...c|ase 3:.X|
|00001c80| 09 09 09 09 74 6d 70 20 | 3d 20 6d 61 74 73 71 75 |....tmp |= matsqu|
|00001c90| 61 72 65 28 6d 29 3b 0a | 58 09 09 09 09 72 65 73 |are(m);.|X....res|
|00001ca0| 20 3d 20 6d 61 74 6d 75 | 6c 28 6d 2c 20 74 6d 70 | = matmu|l(m, tmp|
|00001cb0| 29 3b 0a 58 09 09 09 09 | 6d 61 74 66 72 65 65 28 |);.X....|matfree(|
|00001cc0| 74 6d 70 29 3b 0a 58 09 | 09 09 09 72 65 74 75 72 |tmp);.X.|...retur|
|00001cd0| 6e 20 72 65 73 3b 0a 58 | 09 09 09 63 61 73 65 20 |n res;.X|...case |
|00001ce0| 34 3a 0a 58 09 09 09 09 | 74 6d 70 20 3d 20 6d 61 |4:.X....|tmp = ma|
|00001cf0| 74 73 71 75 61 72 65 28 | 6d 29 3b 0a 58 09 09 09 |tsquare(|m);.X...|
|00001d00| 09 72 65 73 20 3d 20 6d | 61 74 73 71 75 61 72 65 |.res = m|atsquare|
|00001d10| 28 74 6d 70 29 3b 0a 58 | 09 09 09 09 6d 61 74 66 |(tmp);.X|....matf|
|00001d20| 72 65 65 28 74 6d 70 29 | 3b 0a 58 09 09 09 09 72 |ree(tmp)|;.X....r|
|00001d30| 65 74 75 72 6e 20 72 65 | 73 3b 0a 58 09 09 7d 0a |eturn re|s;.X..}.|
|00001d40| 58 09 7d 0a 58 09 69 66 | 20 28 70 6f 77 65 72 20 |X.}.X.if| (power |
|00001d50| 3c 20 30 29 20 7b 0a 58 | 09 09 6d 20 3d 20 6d 61 |< 0) {.X|..m = ma|
|00001d60| 74 69 6e 76 28 6d 29 3b | 0a 58 09 09 70 6f 77 65 |tinv(m);|.X..powe|
|00001d70| 72 20 3d 20 2d 70 6f 77 | 65 72 3b 0a 58 09 7d 0a |r = -pow|er;.X.}.|
|00001d80| 58 09 2f 2a 0a 58 09 20 | 2a 20 43 6f 6d 70 75 74 |X./*.X. |* Comput|
|00001d90| 65 20 74 68 65 20 70 6f | 77 65 72 20 62 79 20 73 |e the po|wer by s|
|00001da0| 71 75 61 72 69 6e 67 20 | 61 6e 64 20 6d 75 6c 74 |quaring |and mult|
|00001db0| 69 70 6c 79 69 6e 67 2e | 0a 58 09 20 2a 20 54 68 |iplying.|.X. * Th|
|00001dc0| 69 73 20 75 73 65 73 20 | 74 68 65 20 6c 65 66 74 |is uses |the left|
|00001dd0| 20 74 6f 20 72 69 67 68 | 74 20 6d 65 74 68 6f 64 | to righ|t method|
|00001de0| 20 6f 66 20 70 6f 77 65 | 72 20 72 61 69 73 69 6e | of powe|r raisin|
|00001df0| 67 2e 0a 58 09 20 2a 2f | 0a 58 09 62 69 74 20 3d |g..X. */|.X.bit =|
|00001e00| 20 54 4f 50 46 55 4c 4c | 3b 0a 58 09 77 68 69 6c | TOPFULL|;.X.whil|
|00001e10| 65 20 28 28 62 69 74 20 | 26 20 70 6f 77 65 72 29 |e ((bit |& power)|
|00001e20| 20 3d 3d 20 30 29 0a 58 | 09 09 62 69 74 20 3e 3e | == 0).X|..bit >>|
|00001e30| 3d 20 31 4c 3b 0a 58 09 | 62 69 74 20 3e 3e 3d 20 |= 1L;.X.|bit >>= |
|00001e40| 31 4c 3b 0a 58 09 72 65 | 73 20 3d 20 6d 61 74 73 |1L;.X.re|s = mats|
|00001e50| 71 75 61 72 65 28 6d 29 | 3b 0a 58 09 69 66 20 28 |quare(m)|;.X.if (|
|00001e60| 62 69 74 20 26 20 70 6f | 77 65 72 29 20 7b 0a 58 |bit & po|wer) {.X|
|00001e70| 09 09 74 6d 70 20 3d 20 | 6d 61 74 6d 75 6c 28 72 |..tmp = |matmul(r|
|00001e80| 65 73 2c 20 6d 29 3b 0a | 58 09 09 6d 61 74 66 72 |es, m);.|X..matfr|
|00001e90| 65 65 28 72 65 73 29 3b | 0a 58 09 09 72 65 73 20 |ee(res);|.X..res |
|00001ea0| 3d 20 74 6d 70 3b 0a 58 | 09 7d 0a 58 09 62 69 74 |= tmp;.X|.}.X.bit|
|00001eb0| 20 3e 3e 3d 20 31 4c 3b | 0a 58 09 77 68 69 6c 65 | >>= 1L;|.X.while|
|00001ec0| 20 28 62 69 74 29 20 7b | 0a 58 09 09 74 6d 70 20 | (bit) {|.X..tmp |
|00001ed0| 3d 20 6d 61 74 73 71 75 | 61 72 65 28 72 65 73 29 |= matsqu|are(res)|
|00001ee0| 3b 0a 58 09 09 6d 61 74 | 66 72 65 65 28 72 65 73 |;.X..mat|free(res|
|00001ef0| 29 3b 0a 58 09 09 72 65 | 73 20 3d 20 74 6d 70 3b |);.X..re|s = tmp;|
|00001f00| 0a 58 09 09 69 66 20 28 | 62 69 74 20 26 20 70 6f |.X..if (|bit & po|
|00001f10| 77 65 72 29 20 7b 0a 58 | 09 09 09 74 6d 70 20 3d |wer) {.X|...tmp =|
|00001f20| 20 6d 61 74 6d 75 6c 28 | 72 65 73 2c 20 6d 29 3b | matmul(|res, m);|
|00001f30| 0a 58 09 09 09 6d 61 74 | 66 72 65 65 28 72 65 73 |.X...mat|free(res|
|00001f40| 29 3b 0a 58 09 09 09 72 | 65 73 20 3d 20 74 6d 70 |);.X...r|es = tmp|
|00001f50| 3b 0a 58 09 09 7d 0a 58 | 09 09 62 69 74 20 3e 3e |;.X..}.X|..bit >>|
|00001f60| 3d 20 31 4c 3b 0a 58 09 | 7d 0a 58 09 69 66 20 28 |= 1L;.X.|}.X.if (|
|00001f70| 71 69 73 6e 65 67 28 71 | 29 29 0a 58 09 09 6d 61 |qisneg(q|)).X..ma|
|00001f80| 74 66 72 65 65 28 6d 29 | 3b 0a 58 09 72 65 74 75 |tfree(m)|;.X.retu|
|00001f90| 72 6e 20 72 65 73 3b 0a | 58 7d 0a 58 0a 58 0a 58 |rn res;.|X}.X.X.X|
|00001fa0| 2f 2a 0a 58 20 2a 20 43 | 61 6c 63 75 6c 61 74 65 |/*.X * C|alculate|
|00001fb0| 20 74 68 65 20 63 72 6f | 73 73 20 70 72 6f 64 75 | the cro|ss produ|
|00001fc0| 63 74 20 6f 66 20 74 77 | 6f 20 6f 6e 65 20 64 69 |ct of tw|o one di|
|00001fd0| 6d 65 6e 73 69 6f 6e 61 | 6c 20 6d 61 74 72 69 63 |mensiona|l matric|
|00001fe0| 65 73 20 65 61 63 68 0a | 58 20 2a 20 77 69 74 68 |es each.|X * with|
|00001ff0| 20 74 68 72 65 65 20 63 | 6f 6d 70 6f 6e 65 6e 74 | three c|omponent|
|00002000| 73 2e 0a 58 20 2a 09 6d | 33 20 3d 20 6d 61 74 63 |s..X *.m|3 = matc|
|00002010| 72 6f 73 73 28 6d 31 2c | 20 6d 32 29 3b 0a 58 20 |ross(m1,| m2);.X |
|00002020| 2a 2f 0a 58 4d 41 54 52 | 49 58 20 2a 0a 58 6d 61 |*/.XMATR|IX *.Xma|
|00002030| 74 63 72 6f 73 73 28 6d | 31 2c 20 6d 32 29 0a 58 |tcross(m|1, m2).X|
|00002040| 09 4d 41 54 52 49 58 20 | 2a 6d 31 2c 20 2a 6d 32 |.MATRIX |*m1, *m2|
|00002050| 3b 0a 58 7b 0a 58 09 4d | 41 54 52 49 58 20 2a 72 |;.X{.X.M|ATRIX *r|
|00002060| 65 73 3b 0a 58 09 56 41 | 4c 55 45 20 2a 76 31 2c |es;.X.VA|LUE *v1,|
|00002070| 20 2a 76 32 2c 20 2a 76 | 72 3b 0a 58 09 56 41 4c | *v2, *v|r;.X.VAL|
|00002080| 55 45 20 74 6d 70 31 2c | 20 74 6d 70 32 3b 0a 58 |UE tmp1,| tmp2;.X|
|00002090| 0a 58 09 69 66 20 28 28 | 6d 31 2d 3e 6d 5f 64 69 |.X.if ((|m1->m_di|
|000020a0| 6d 20 21 3d 20 31 29 20 | 7c 7c 20 28 6d 32 2d 3e |m != 1) ||| (m2->|
|000020b0| 6d 5f 64 69 6d 20 21 3d | 20 31 29 29 0a 58 09 09 |m_dim !=| 1)).X..|
|000020c0| 65 72 72 6f 72 28 22 4d | 61 74 72 69 78 20 6e 6f |error("M|atrix no|
|000020d0| 74 20 31 64 20 66 6f 72 | 20 63 72 6f 73 73 20 70 |t 1d for| cross p|
|000020e0| 72 6f 64 75 63 74 22 29 | 3b 0a 58 09 69 66 20 28 |roduct")|;.X.if (|
|000020f0| 28 6d 31 2d 3e 6d 5f 73 | 69 7a 65 20 21 3d 20 33 |(m1->m_s|ize != 3|
|00002100| 29 20 7c 7c 20 28 6d 32 | 2d 3e 6d 5f 73 69 7a 65 |) || (m2|->m_size|
|00002110| 20 21 3d 20 33 29 29 0a | 58 09 09 65 72 72 6f 72 | != 3)).|X..error|
|00002120| 28 22 4d 61 74 72 69 78 | 20 6e 6f 74 20 73 69 7a |("Matrix| not siz|
|00002130| 65 20 33 20 66 6f 72 20 | 63 72 6f 73 73 20 70 72 |e 3 for |cross pr|
|00002140| 6f 64 75 63 74 22 29 3b | 0a 58 09 72 65 73 20 3d |oduct");|.X.res =|
|00002150| 20 6d 61 74 61 6c 6c 6f | 63 28 33 4c 29 3b 0a 58 | matallo|c(3L);.X|
|00002160| 09 72 65 73 2d 3e 6d 5f | 64 69 6d 20 3d 20 31 3b |.res->m_|dim = 1;|
|00002170| 0a 58 09 72 65 73 2d 3e | 6d 5f 6d 69 6e 5b 30 5d |.X.res->|m_min[0]|
|00002180| 20 3d 20 30 3b 0a 58 09 | 72 65 73 2d 3e 6d 5f 6d | = 0;.X.|res->m_m|
|00002190| 61 78 5b 30 5d 20 3d 20 | 32 3b 0a 58 09 76 31 20 |ax[0] = |2;.X.v1 |
|000021a0| 3d 20 6d 31 2d 3e 6d 5f | 74 61 62 6c 65 3b 0a 58 |= m1->m_|table;.X|
|000021b0| 09 76 32 20 3d 20 6d 32 | 2d 3e 6d 5f 74 61 62 6c |.v2 = m2|->m_tabl|
|000021c0| 65 3b 0a 58 09 76 72 20 | 3d 20 72 65 73 2d 3e 6d |e;.X.vr |= res->m|
|000021d0| 5f 74 61 62 6c 65 3b 0a | 58 09 6d 75 6c 76 61 6c |_table;.|X.mulval|
|000021e0| 75 65 28 76 31 20 2b 20 | 31 2c 20 76 32 20 2b 20 |ue(v1 + |1, v2 + |
|000021f0| 32 2c 20 26 74 6d 70 31 | 29 3b 0a 58 09 6d 75 6c |2, &tmp1|);.X.mul|
|00002200| 76 61 6c 75 65 28 76 31 | 20 2b 20 32 2c 20 76 32 |value(v1| + 2, v2|
|00002210| 20 2b 20 31 2c 20 26 74 | 6d 70 32 29 3b 0a 58 09 | + 1, &t|mp2);.X.|
|00002220| 73 75 62 76 61 6c 75 65 | 28 26 74 6d 70 31 2c 20 |subvalue|(&tmp1, |
|00002230| 26 74 6d 70 32 2c 20 76 | 72 20 2b 20 30 29 3b 0a |&tmp2, v|r + 0);.|
|00002240| 58 09 66 72 65 65 76 61 | 6c 75 65 28 26 74 6d 70 |X.freeva|lue(&tmp|
|00002250| 31 29 3b 0a 58 09 66 72 | 65 65 76 61 6c 75 65 28 |1);.X.fr|eevalue(|
|00002260| 26 74 6d 70 32 29 3b 0a | 58 09 6d 75 6c 76 61 6c |&tmp2);.|X.mulval|
|00002270| 75 65 28 76 31 20 2b 20 | 32 2c 20 76 32 20 2b 20 |ue(v1 + |2, v2 + |
|00002280| 30 2c 20 26 74 6d 70 31 | 29 3b 0a 58 09 6d 75 6c |0, &tmp1|);.X.mul|
|00002290| 76 61 6c 75 65 28 76 31 | 20 2b 20 30 2c 20 76 32 |value(v1| + 0, v2|
|000022a0| 20 2b 20 32 2c 20 26 74 | 6d 70 32 29 3b 0a 58 09 | + 2, &t|mp2);.X.|
|000022b0| 73 75 62 76 61 6c 75 65 | 28 26 74 6d 70 31 2c 20 |subvalue|(&tmp1, |
|000022c0| 26 74 6d 70 32 2c 20 76 | 72 20 2b 20 31 29 3b 0a |&tmp2, v|r + 1);.|
|000022d0| 58 09 66 72 65 65 76 61 | 6c 75 65 28 26 74 6d 70 |X.freeva|lue(&tmp|
|000022e0| 31 29 3b 0a 58 09 66 72 | 65 65 76 61 6c 75 65 28 |1);.X.fr|eevalue(|
|000022f0| 26 74 6d 70 32 29 3b 0a | 58 09 6d 75 6c 76 61 6c |&tmp2);.|X.mulval|
|00002300| 75 65 28 76 31 20 2b 20 | 30 2c 20 76 32 20 2b 20 |ue(v1 + |0, v2 + |
|00002310| 31 2c 20 26 74 6d 70 31 | 29 3b 0a 58 09 6d 75 6c |1, &tmp1|);.X.mul|
|00002320| 76 61 6c 75 65 28 76 31 | 20 2b 20 31 2c 20 76 32 |value(v1| + 1, v2|
|00002330| 20 2b 20 30 2c 20 26 74 | 6d 70 32 29 3b 0a 58 09 | + 0, &t|mp2);.X.|
|00002340| 73 75 62 76 61 6c 75 65 | 28 26 74 6d 70 31 2c 20 |subvalue|(&tmp1, |
|00002350| 26 74 6d 70 32 2c 20 76 | 72 20 2b 20 32 29 3b 0a |&tmp2, v|r + 2);.|
|00002360| 58 09 66 72 65 65 76 61 | 6c 75 65 28 26 74 6d 70 |X.freeva|lue(&tmp|
|00002370| 31 29 3b 0a 58 09 66 72 | 65 65 76 61 6c 75 65 28 |1);.X.fr|eevalue(|
|00002380| 26 74 6d 70 32 29 3b 0a | 58 09 72 65 74 75 72 6e |&tmp2);.|X.return|
|00002390| 20 72 65 73 3b 0a 58 7d | 0a 58 0a 58 0a 58 2f 2a | res;.X}|.X.X.X/*|
|000023a0| 0a 58 20 2a 20 52 65 74 | 75 72 6e 20 74 68 65 20 |.X * Ret|urn the |
|000023b0| 64 6f 74 20 70 72 6f 64 | 75 63 74 20 6f 66 20 74 |dot prod|uct of t|
|000023c0| 77 6f 20 6d 61 74 72 69 | 63 65 73 2e 0a 58 20 2a |wo matri|ces..X *|
|000023d0| 09 72 65 73 75 6c 74 20 | 3d 20 6d 61 74 64 6f 74 |.result |= matdot|
|000023e0| 28 6d 31 2c 20 6d 32 29 | 3b 0a 58 20 2a 2f 0a 58 |(m1, m2)|;.X */.X|
|000023f0| 56 41 4c 55 45 0a 58 6d | 61 74 64 6f 74 28 6d 31 |VALUE.Xm|atdot(m1|
|00002400| 2c 20 6d 32 29 0a 58 09 | 4d 41 54 52 49 58 20 2a |, m2).X.|MATRIX *|
|00002410| 6d 31 2c 20 2a 6d 32 3b | 0a 58 7b 0a 58 09 56 41 |m1, *m2;|.X{.X.VA|
|00002420| 4c 55 45 20 2a 76 31 2c | 20 2a 76 32 3b 0a 58 09 |LUE *v1,| *v2;.X.|
|00002430| 56 41 4c 55 45 20 72 65 | 73 75 6c 74 2c 20 74 6d |VALUE re|sult, tm|
|00002440| 70 31 2c 20 74 6d 70 32 | 3b 0a 58 09 6c 6f 6e 67 |p1, tmp2|;.X.long|
|00002450| 20 6c 65 6e 3b 0a 58 0a | 58 09 69 66 20 28 28 6d | len;.X.|X.if ((m|
|00002460| 31 2d 3e 6d 5f 64 69 6d | 20 21 3d 20 31 29 20 7c |1->m_dim| != 1) ||
|00002470| 7c 20 28 6d 32 2d 3e 6d | 5f 64 69 6d 20 21 3d 20 || (m2->m|_dim != |
|00002480| 31 29 29 0a 58 09 09 65 | 72 72 6f 72 28 22 4d 61 |1)).X..e|rror("Ma|
|00002490| 74 72 69 78 20 6e 6f 74 | 20 31 64 20 66 6f 72 20 |trix not| 1d for |
|000024a0| 64 6f 74 20 70 72 6f 64 | 75 63 74 22 29 3b 0a 58 |dot prod|uct");.X|
|000024b0| 09 69 66 20 28 6d 31 2d | 3e 6d 5f 73 69 7a 65 20 |.if (m1-|>m_size |
|000024c0| 21 3d 20 6d 32 2d 3e 6d | 5f 73 69 7a 65 29 0a 58 |!= m2->m|_size).X|
|000024d0| 09 09 65 72 72 6f 72 28 | 22 49 6e 63 6f 6d 70 61 |..error(|"Incompa|
|000024e0| 74 69 62 6c 65 20 6d 61 | 74 72 69 78 20 73 69 7a |tible ma|trix siz|
|000024f0| 65 73 20 66 6f 72 20 64 | 6f 74 20 70 72 6f 64 75 |es for d|ot produ|
|00002500| 63 74 22 29 3b 0a 58 09 | 76 31 20 3d 20 6d 31 2d |ct");.X.|v1 = m1-|
|00002510| 3e 6d 5f 74 61 62 6c 65 | 3b 0a 58 09 76 32 20 3d |>m_table|;.X.v2 =|
|00002520| 20 6d 32 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 | m2->m_t|able;.X.|
|00002530| 6d 75 6c 76 61 6c 75 65 | 28 76 31 2c 20 76 32 2c |mulvalue|(v1, v2,|
|00002540| 20 26 72 65 73 75 6c 74 | 29 3b 0a 58 09 6c 65 6e | &result|);.X.len|
|00002550| 20 3d 20 6d 31 2d 3e 6d | 5f 73 69 7a 65 3b 0a 58 | = m1->m|_size;.X|
|00002560| 09 77 68 69 6c 65 20 28 | 2d 2d 6c 65 6e 20 3e 20 |.while (|--len > |
|00002570| 30 29 20 7b 0a 58 09 09 | 6d 75 6c 76 61 6c 75 65 |0) {.X..|mulvalue|
|00002580| 28 2b 2b 76 31 2c 20 2b | 2b 76 32 2c 20 26 74 6d |(++v1, +|+v2, &tm|
|00002590| 70 31 29 3b 0a 58 09 09 | 61 64 64 76 61 6c 75 65 |p1);.X..|addvalue|
|000025a0| 28 26 72 65 73 75 6c 74 | 2c 20 26 74 6d 70 31 2c |(&result|, &tmp1,|
|000025b0| 20 26 74 6d 70 32 29 3b | 0a 58 09 09 66 72 65 65 | &tmp2);|.X..free|
|000025c0| 76 61 6c 75 65 28 26 74 | 6d 70 31 29 3b 0a 58 09 |value(&t|mp1);.X.|
|000025d0| 09 66 72 65 65 76 61 6c | 75 65 28 26 72 65 73 75 |.freeval|ue(&resu|
|000025e0| 6c 74 29 3b 0a 58 09 09 | 72 65 73 75 6c 74 20 3d |lt);.X..|result =|
|000025f0| 20 74 6d 70 32 3b 0a 58 | 09 7d 0a 58 09 72 65 74 | tmp2;.X|.}.X.ret|
|00002600| 75 72 6e 20 72 65 73 75 | 6c 74 3b 0a 58 7d 0a 58 |urn resu|lt;.X}.X|
|00002610| 0a 58 0a 58 2f 2a 0a 58 | 20 2a 20 53 63 61 6c 65 |.X.X/*.X| * Scale|
|00002620| 20 74 68 65 20 65 6c 65 | 6d 65 6e 74 73 20 6f 66 | the ele|ments of|
|00002630| 20 61 20 6d 61 74 72 69 | 78 20 62 79 20 61 20 73 | a matri|x by a s|
|00002640| 70 65 63 69 66 69 65 64 | 20 70 6f 77 65 72 20 6f |pecified| power o|
|00002650| 66 20 74 77 6f 2e 0a 58 | 20 2a 2f 0a 58 4d 41 54 |f two..X| */.XMAT|
|00002660| 52 49 58 20 2a 0a 58 6d | 61 74 73 63 61 6c 65 28 |RIX *.Xm|atscale(|
|00002670| 6d 2c 20 6e 29 0a 58 09 | 4d 41 54 52 49 58 20 2a |m, n).X.|MATRIX *|
|00002680| 6d 3b 09 09 2f 2a 20 6d | 61 74 72 69 78 20 74 6f |m;../* m|atrix to|
|00002690| 20 62 65 20 73 63 61 6c | 65 64 20 2a 2f 0a 58 09 | be scal|ed */.X.|
|000026a0| 6c 6f 6e 67 20 6e 3b 0a | 58 7b 0a 58 09 72 65 67 |long n;.|X{.X.reg|
|000026b0| 69 73 74 65 72 20 56 41 | 4c 55 45 20 2a 76 61 6c |ister VA|LUE *val|
|000026c0| 2c 20 2a 76 72 65 73 3b | 0a 58 09 56 41 4c 55 45 |, *vres;|.X.VALUE|
|000026d0| 20 6e 75 6d 3b 0a 58 09 | 6c 6f 6e 67 20 69 6e 64 | num;.X.|long ind|
|000026e0| 65 78 3b 0a 58 09 4d 41 | 54 52 49 58 20 2a 72 65 |ex;.X.MA|TRIX *re|
|000026f0| 73 3b 09 09 2f 2a 20 72 | 65 73 75 6c 74 69 6e 67 |s;../* r|esulting|
|00002700| 20 6d 61 74 72 69 78 20 | 2a 2f 0a 58 0a 58 09 69 | matrix |*/.X.X.i|
|00002710| 66 20 28 6e 20 3d 3d 20 | 30 29 0a 58 09 09 72 65 |f (n == |0).X..re|
|00002720| 74 75 72 6e 20 6d 61 74 | 63 6f 70 79 28 6d 29 3b |turn mat|copy(m);|
|00002730| 0a 58 09 6e 75 6d 2e 76 | 5f 74 79 70 65 20 3d 20 |.X.num.v|_type = |
|00002740| 56 5f 4e 55 4d 3b 0a 58 | 09 6e 75 6d 2e 76 5f 6e |V_NUM;.X|.num.v_n|
|00002750| 75 6d 20 3d 20 69 74 6f | 71 28 6e 29 3b 0a 58 09 |um = ito|q(n);.X.|
|00002760| 72 65 73 20 3d 20 6d 61 | 74 61 6c 6c 6f 63 28 6d |res = ma|talloc(m|
|00002770| 2d 3e 6d 5f 73 69 7a 65 | 29 3b 0a 58 09 2a 72 65 |->m_size|);.X.*re|
|00002780| 73 20 3d 20 2a 6d 3b 0a | 58 09 76 61 6c 20 3d 20 |s = *m;.|X.val = |
|00002790| 6d 2d 3e 6d 5f 74 61 62 | 6c 65 3b 0a 58 09 76 72 |m->m_tab|le;.X.vr|
|000027a0| 65 73 20 3d 20 72 65 73 | 2d 3e 6d 5f 74 61 62 6c |es = res|->m_tabl|
|000027b0| 65 3b 0a 58 09 66 6f 72 | 20 28 69 6e 64 65 78 20 |e;.X.for| (index |
|000027c0| 3d 20 6d 2d 3e 6d 5f 73 | 69 7a 65 3b 20 69 6e 64 |= m->m_s|ize; ind|
|000027d0| 65 78 20 3e 20 30 3b 20 | 69 6e 64 65 78 2d 2d 29 |ex > 0; |index--)|
|000027e0| 0a 58 09 09 73 63 61 6c | 65 76 61 6c 75 65 28 76 |.X..scal|evalue(v|
|000027f0| 61 6c 2b 2b 2c 20 26 6e | 75 6d 2c 20 76 72 65 73 |al++, &n|um, vres|
|00002800| 2b 2b 29 3b 0a 58 09 71 | 66 72 65 65 28 6e 75 6d |++);.X.q|free(num|
|00002810| 2e 76 5f 6e 75 6d 29 3b | 0a 58 09 72 65 74 75 72 |.v_num);|.X.retur|
|00002820| 6e 20 72 65 73 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |n res;.X|}.X.X.X/|
|00002830| 2a 0a 58 20 2a 20 53 68 | 69 66 74 20 74 68 65 20 |*.X * Sh|ift the |
|00002840| 65 6c 65 6d 65 6e 74 73 | 20 6f 66 20 61 20 6d 61 |elements| of a ma|
|00002850| 74 72 69 78 20 62 79 20 | 74 68 65 20 73 70 65 63 |trix by |the spec|
|00002860| 69 66 69 65 64 20 6e 75 | 6d 62 65 72 20 6f 66 20 |ified nu|mber of |
|00002870| 62 69 74 73 2e 0a 58 20 | 2a 20 50 6f 73 69 74 69 |bits..X |* Positi|
|00002880| 76 65 20 73 68 69 66 74 | 20 6d 65 61 6e 73 20 6c |ve shift| means l|
|00002890| 65 66 74 77 61 72 64 73 | 2c 20 6e 65 67 61 74 69 |eftwards|, negati|
|000028a0| 76 65 20 73 68 69 66 74 | 20 72 69 67 68 74 77 61 |ve shift| rightwa|
|000028b0| 72 64 73 2e 0a 58 20 2a | 2f 0a 58 4d 41 54 52 49 |rds..X *|/.XMATRI|
|000028c0| 58 20 2a 0a 58 6d 61 74 | 73 68 69 66 74 28 6d 2c |X *.Xmat|shift(m,|
|000028d0| 20 6e 29 0a 58 09 4d 41 | 54 52 49 58 20 2a 6d 3b | n).X.MA|TRIX *m;|
|000028e0| 09 09 2f 2a 20 6d 61 74 | 72 69 78 20 74 6f 20 62 |../* mat|rix to b|
|000028f0| 65 20 73 63 61 6c 65 64 | 20 2a 2f 0a 58 09 6c 6f |e scaled| */.X.lo|
|00002900| 6e 67 20 6e 3b 0a 58 7b | 0a 58 09 72 65 67 69 73 |ng n;.X{|.X.regis|
|00002910| 74 65 72 20 56 41 4c 55 | 45 20 2a 76 61 6c 2c 20 |ter VALU|E *val, |
|00002920| 2a 76 72 65 73 3b 0a 58 | 09 56 41 4c 55 45 20 6e |*vres;.X|.VALUE n|
|00002930| 75 6d 3b 0a 58 09 6c 6f | 6e 67 20 69 6e 64 65 78 |um;.X.lo|ng index|
|00002940| 3b 0a 58 09 4d 41 54 52 | 49 58 20 2a 72 65 73 3b |;.X.MATR|IX *res;|
|00002950| 09 09 2f 2a 20 72 65 73 | 75 6c 74 69 6e 67 20 6d |../* res|ulting m|
|00002960| 61 74 72 69 78 20 2a 2f | 0a 58 0a 58 09 69 66 20 |atrix */|.X.X.if |
|00002970| 28 6e 20 3d 3d 20 30 29 | 0a 58 09 09 72 65 74 75 |(n == 0)|.X..retu|
|00002980| 72 6e 20 6d 61 74 63 6f | 70 79 28 6d 29 3b 0a 58 |rn matco|py(m);.X|
|00002990| 09 6e 75 6d 2e 76 5f 74 | 79 70 65 20 3d 20 56 5f |.num.v_t|ype = V_|
|000029a0| 4e 55 4d 3b 0a 58 09 6e | 75 6d 2e 76 5f 6e 75 6d |NUM;.X.n|um.v_num|
|000029b0| 20 3d 20 69 74 6f 71 28 | 6e 29 3b 0a 58 09 72 65 | = itoq(|n);.X.re|
|000029c0| 73 20 3d 20 6d 61 74 61 | 6c 6c 6f 63 28 6d 2d 3e |s = mata|lloc(m->|
|000029d0| 6d 5f 73 69 7a 65 29 3b | 0a 58 09 2a 72 65 73 20 |m_size);|.X.*res |
|000029e0| 3d 20 2a 6d 3b 0a 58 09 | 76 61 6c 20 3d 20 6d 2d |= *m;.X.|val = m-|
|000029f0| 3e 6d 5f 74 61 62 6c 65 | 3b 0a 58 09 76 72 65 73 |>m_table|;.X.vres|
|00002a00| 20 3d 20 72 65 73 2d 3e | 6d 5f 74 61 62 6c 65 3b | = res->|m_table;|
|00002a10| 0a 58 09 66 6f 72 20 28 | 69 6e 64 65 78 20 3d 20 |.X.for (|index = |
|00002a20| 6d 2d 3e 6d 5f 73 69 7a | 65 3b 20 69 6e 64 65 78 |m->m_siz|e; index|
|00002a30| 20 3e 20 30 3b 20 69 6e | 64 65 78 2d 2d 29 0a 58 | > 0; in|dex--).X|
|00002a40| 09 09 73 68 69 66 74 76 | 61 6c 75 65 28 76 61 6c |..shiftv|alue(val|
|00002a50| 2b 2b 2c 20 26 6e 75 6d | 2c 20 46 41 4c 53 45 2c |++, &num|, FALSE,|
|00002a60| 20 76 72 65 73 2b 2b 29 | 3b 0a 58 09 71 66 72 65 | vres++)|;.X.qfre|
|00002a70| 65 28 6e 75 6d 2e 76 5f | 6e 75 6d 29 3b 0a 58 09 |e(num.v_|num);.X.|
|00002a80| 72 65 74 75 72 6e 20 72 | 65 73 3b 0a 58 7d 0a 58 |return r|es;.X}.X|
|00002a90| 0a 58 0a 58 2f 2a 0a 58 | 20 2a 20 4d 75 6c 74 69 |.X.X/*.X| * Multi|
|00002aa0| 70 6c 79 20 74 68 65 20 | 65 6c 65 6d 65 6e 74 73 |ply the |elements|
|00002ab0| 20 6f 66 20 61 20 6d 61 | 74 72 69 78 20 62 79 20 | of a ma|trix by |
|00002ac0| 61 20 73 70 65 63 69 66 | 69 65 64 20 76 61 6c 75 |a specif|ied valu|
|00002ad0| 65 2e 0a 58 20 2a 2f 0a | 58 4d 41 54 52 49 58 20 |e..X */.|XMATRIX |
|00002ae0| 2a 0a 58 6d 61 74 6d 75 | 6c 76 61 6c 28 6d 2c 20 |*.Xmatmu|lval(m, |
|00002af0| 76 70 29 0a 58 09 4d 41 | 54 52 49 58 20 2a 6d 3b |vp).X.MA|TRIX *m;|
|00002b00| 09 09 2f 2a 20 6d 61 74 | 72 69 78 20 74 6f 20 62 |../* mat|rix to b|
|00002b10| 65 20 6d 75 6c 74 69 70 | 6c 69 65 64 20 2a 2f 0a |e multip|lied */.|
|00002b20| 58 09 56 41 4c 55 45 20 | 2a 76 70 3b 09 09 2f 2a |X.VALUE |*vp;../*|
|00002b30| 20 76 61 6c 75 65 20 74 | 6f 20 6d 75 6c 74 69 70 | value t|o multip|
|00002b40| 6c 79 20 62 79 20 2a 2f | 0a 58 7b 0a 58 09 72 65 |ly by */|.X{.X.re|
|00002b50| 67 69 73 74 65 72 20 56 | 41 4c 55 45 20 2a 76 61 |gister V|ALUE *va|
|00002b60| 6c 2c 20 2a 76 72 65 73 | 3b 0a 58 09 6c 6f 6e 67 |l, *vres|;.X.long|
|00002b70| 20 69 6e 64 65 78 3b 0a | 58 09 4d 41 54 52 49 58 | index;.|X.MATRIX|
|00002b80| 20 2a 72 65 73 3b 0a 58 | 0a 58 09 72 65 73 20 3d | *res;.X|.X.res =|
|00002b90| 20 6d 61 74 61 6c 6c 6f | 63 28 6d 2d 3e 6d 5f 73 | matallo|c(m->m_s|
|00002ba0| 69 7a 65 29 3b 0a 58 09 | 2a 72 65 73 20 3d 20 2a |ize);.X.|*res = *|
|00002bb0| 6d 3b 0a 58 09 76 61 6c | 20 3d 20 6d 2d 3e 6d 5f |m;.X.val| = m->m_|
|00002bc0| 74 61 62 6c 65 3b 0a 58 | 09 76 72 65 73 20 3d 20 |table;.X|.vres = |
|00002bd0| 72 65 73 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 |res->m_t|able;.X.|
|00002be0| 66 6f 72 20 28 69 6e 64 | 65 78 20 3d 20 6d 2d 3e |for (ind|ex = m->|
|00002bf0| 6d 5f 73 69 7a 65 3b 20 | 69 6e 64 65 78 20 3e 20 |m_size; |index > |
|00002c00| 30 3b 20 69 6e 64 65 78 | 2d 2d 29 0a 58 09 09 6d |0; index|--).X..m|
|00002c10| 75 6c 76 61 6c 75 65 28 | 76 61 6c 2b 2b 2c 20 76 |ulvalue(|val++, v|
|00002c20| 70 2c 20 76 72 65 73 2b | 2b 29 3b 0a 58 09 72 65 |p, vres+|+);.X.re|
|00002c30| 74 75 72 6e 20 72 65 73 | 3b 0a 58 7d 0a 58 0a 58 |turn res|;.X}.X.X|
|00002c40| 0a 58 2f 2a 0a 58 20 2a | 20 44 69 76 69 64 65 20 |.X/*.X *| Divide |
|00002c50| 74 68 65 20 65 6c 65 6d | 65 6e 74 73 20 6f 66 20 |the elem|ents of |
|00002c60| 61 20 6d 61 74 72 69 78 | 20 62 79 20 61 20 73 70 |a matrix| by a sp|
|00002c70| 65 63 69 66 69 65 64 20 | 76 61 6c 75 65 2c 20 6b |ecified |value, k|
|00002c80| 65 65 70 69 6e 67 0a 58 | 20 2a 20 6f 6e 6c 79 20 |eeping.X| * only |
|00002c90| 74 68 65 20 69 6e 74 65 | 67 65 72 20 71 75 6f 74 |the inte|ger quot|
|00002ca0| 69 65 6e 74 2e 0a 58 20 | 2a 2f 0a 58 4d 41 54 52 |ient..X |*/.XMATR|
|00002cb0| 49 58 20 2a 0a 58 6d 61 | 74 71 75 6f 76 61 6c 28 |IX *.Xma|tquoval(|
|00002cc0| 6d 2c 20 76 70 29 0a 58 | 09 4d 41 54 52 49 58 20 |m, vp).X|.MATRIX |
|00002cd0| 2a 6d 3b 09 09 2f 2a 20 | 6d 61 74 72 69 78 20 74 |*m;../* |matrix t|
|00002ce0| 6f 20 62 65 20 64 69 76 | 69 64 65 64 20 2a 2f 0a |o be div|ided */.|
|00002cf0| 58 09 56 41 4c 55 45 20 | 2a 76 70 3b 09 09 2f 2a |X.VALUE |*vp;../*|
|00002d00| 20 76 61 6c 75 65 20 74 | 6f 20 64 69 76 69 64 65 | value t|o divide|
|00002d10| 20 62 79 20 2a 2f 0a 58 | 7b 0a 58 09 72 65 67 69 | by */.X|{.X.regi|
|00002d20| 73 74 65 72 20 56 41 4c | 55 45 20 2a 76 61 6c 2c |ster VAL|UE *val,|
|00002d30| 20 2a 76 72 65 73 3b 0a | 58 09 6c 6f 6e 67 20 69 | *vres;.|X.long i|
|00002d40| 6e 64 65 78 3b 0a 58 09 | 4d 41 54 52 49 58 20 2a |ndex;.X.|MATRIX *|
|00002d50| 72 65 73 3b 0a 58 0a 58 | 09 69 66 20 28 28 76 70 |res;.X.X|.if ((vp|
|00002d60| 2d 3e 76 5f 74 79 70 65 | 20 3d 3d 20 56 5f 4e 55 |->v_type| == V_NU|
|00002d70| 4d 29 20 26 26 20 71 69 | 73 7a 65 72 6f 28 76 70 |M) && qi|szero(vp|
|00002d80| 2d 3e 76 5f 6e 75 6d 29 | 29 0a 58 09 09 65 72 72 |->v_num)|).X..err|
|00002d90| 6f 72 28 22 44 69 76 69 | 73 69 6f 6e 20 62 79 20 |or("Divi|sion by |
|00002da0| 7a 65 72 6f 22 29 3b 0a | 58 09 72 65 73 20 3d 20 |zero");.|X.res = |
|00002db0| 6d 61 74 61 6c 6c 6f 63 | 28 6d 2d 3e 6d 5f 73 69 |matalloc|(m->m_si|
|00002dc0| 7a 65 29 3b 0a 58 09 2a | 72 65 73 20 3d 20 2a 6d |ze);.X.*|res = *m|
|00002dd0| 3b 0a 58 09 76 61 6c 20 | 3d 20 6d 2d 3e 6d 5f 74 |;.X.val |= m->m_t|
|00002de0| 61 62 6c 65 3b 0a 58 09 | 76 72 65 73 20 3d 20 72 |able;.X.|vres = r|
|00002df0| 65 73 2d 3e 6d 5f 74 61 | 62 6c 65 3b 0a 58 09 66 |es->m_ta|ble;.X.f|
|00002e00| 6f 72 20 28 69 6e 64 65 | 78 20 3d 20 6d 2d 3e 6d |or (inde|x = m->m|
|00002e10| 5f 73 69 7a 65 3b 20 69 | 6e 64 65 78 20 3e 20 30 |_size; i|ndex > 0|
|00002e20| 3b 20 69 6e 64 65 78 2d | 2d 29 0a 58 09 09 71 75 |; index-|-).X..qu|
|00002e30| 6f 76 61 6c 75 65 28 76 | 61 6c 2b 2b 2c 20 76 70 |ovalue(v|al++, vp|
|00002e40| 2c 20 76 72 65 73 2b 2b | 29 3b 0a 58 09 72 65 74 |, vres++|);.X.ret|
|00002e50| 75 72 6e 20 72 65 73 3b | 0a 58 7d 0a 58 0a 58 0a |urn res;|.X}.X.X.|
|00002e60| 58 2f 2a 0a 58 20 2a 20 | 44 69 76 69 64 65 20 74 |X/*.X * |Divide t|
|00002e70| 68 65 20 65 6c 65 6d 65 | 6e 74 73 20 6f 66 20 61 |he eleme|nts of a|
|00002e80| 20 6d 61 74 72 69 78 20 | 62 79 20 61 20 73 70 65 | matrix |by a spe|
|00002e90| 63 69 66 69 65 64 20 76 | 61 6c 75 65 2c 20 6b 65 |cified v|alue, ke|
|00002ea0| 65 70 69 6e 67 0a 58 20 | 2a 20 6f 6e 6c 79 20 74 |eping.X |* only t|
|00002eb0| 68 65 20 72 65 6d 61 69 | 6e 64 65 72 20 6f 66 20 |he remai|nder of |
|00002ec0| 74 68 65 20 64 69 76 69 | 73 69 6f 6e 2e 0a 58 20 |the divi|sion..X |
|00002ed0| 2a 2f 0a 58 4d 41 54 52 | 49 58 20 2a 0a 58 6d 61 |*/.XMATR|IX *.Xma|
|00002ee0| 74 6d 6f 64 76 61 6c 28 | 6d 2c 20 76 70 29 0a 58 |tmodval(|m, vp).X|
|00002ef0| 09 4d 41 54 52 49 58 20 | 2a 6d 3b 09 09 2f 2a 20 |.MATRIX |*m;../* |
|00002f00| 6d 61 74 72 69 78 20 74 | 6f 20 62 65 20 64 69 76 |matrix t|o be div|
|00002f10| 69 64 65 64 20 2a 2f 0a | 58 09 56 41 4c 55 45 20 |ided */.|X.VALUE |
|00002f20| 2a 76 70 3b 09 09 2f 2a | 20 76 61 6c 75 65 20 74 |*vp;../*| value t|
|00002f30| 6f 20 64 69 76 69 64 65 | 20 62 79 20 2a 2f 0a 58 |o divide| by */.X|
|00002f40| 7b 0a 58 09 72 65 67 69 | 73 74 65 72 20 56 41 4c |{.X.regi|ster VAL|
|00002f50| 55 45 20 2a 76 61 6c 2c | 20 2a 76 72 65 73 3b 0a |UE *val,| *vres;.|
|00002f60| 58 09 6c 6f 6e 67 20 69 | 6e 64 65 78 3b 0a 58 09 |X.long i|ndex;.X.|
|00002f70| 4d 41 54 52 49 58 20 2a | 72 65 73 3b 0a 58 0a 58 |MATRIX *|res;.X.X|
|00002f80| 09 69 66 20 28 28 76 70 | 2d 3e 76 5f 74 79 70 65 |.if ((vp|->v_type|
|00002f90| 20 3d 3d 20 56 5f 4e 55 | 4d 29 20 26 26 20 71 69 | == V_NU|M) && qi|
|00002fa0| 73 7a 65 72 6f 28 76 70 | 2d 3e 76 5f 6e 75 6d 29 |szero(vp|->v_num)|
|00002fb0| 29 0a 58 09 09 65 72 72 | 6f 72 28 22 44 69 76 69 |).X..err|or("Divi|
|00002fc0| 73 69 6f 6e 20 62 79 20 | 7a 65 72 6f 22 29 3b 0a |sion by |zero");.|
|00002fd0| 58 09 72 65 73 20 3d 20 | 6d 61 74 61 6c 6c 6f 63 |X.res = |matalloc|
|00002fe0| 28 6d 2d 3e 6d 5f 73 69 | 7a 65 29 3b 0a 58 09 2a |(m->m_si|ze);.X.*|
|00002ff0| 72 65 73 20 3d 20 2a 6d | 3b 0a 58 09 76 61 6c 20 |res = *m|;.X.val |
|00003000| 3d 20 6d 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 |= m->m_t|able;.X.|
|00003010| 76 72 65 73 20 3d 20 72 | 65 73 2d 3e 6d 5f 74 61 |vres = r|es->m_ta|
|00003020| 62 6c 65 3b 0a 58 09 66 | 6f 72 20 28 69 6e 64 65 |ble;.X.f|or (inde|
|00003030| 78 20 3d 20 6d 2d 3e 6d | 5f 73 69 7a 65 3b 20 69 |x = m->m|_size; i|
|00003040| 6e 64 65 78 20 3e 20 30 | 3b 20 69 6e 64 65 78 2d |ndex > 0|; index-|
|00003050| 2d 29 0a 58 09 09 6d 6f | 64 76 61 6c 75 65 28 76 |-).X..mo|dvalue(v|
|00003060| 61 6c 2b 2b 2c 20 76 70 | 2c 20 76 72 65 73 2b 2b |al++, vp|, vres++|
|00003070| 29 3b 0a 58 09 72 65 74 | 75 72 6e 20 72 65 73 3b |);.X.ret|urn res;|
|00003080| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 0a 58 20 2a 20 |.X}.X.X.|X/*.X * |
|00003090| 54 72 61 6e 73 70 6f 73 | 65 20 74 68 65 20 65 6c |Transpos|e the el|
|000030a0| 65 6d 65 6e 74 73 20 6f | 66 20 61 20 74 77 6f 20 |ements o|f a two |
|000030b0| 64 69 6d 65 6e 73 69 6f | 6e 61 6c 20 6d 61 74 72 |dimensio|nal matr|
|000030c0| 69 78 2e 0a 58 20 2a 2f | 0a 58 4d 41 54 52 49 58 |ix..X */|.XMATRIX|
|000030d0| 20 2a 0a 58 6d 61 74 74 | 72 61 6e 73 28 6d 29 0a | *.Xmatt|rans(m).|
|000030e0| 58 09 4d 41 54 52 49 58 | 20 2a 6d 3b 09 09 09 2f |X.MATRIX| *m;.../|
|000030f0| 2a 20 6d 61 74 72 69 78 | 20 74 6f 20 62 65 20 74 |* matrix| to be t|
|00003100| 72 61 6e 73 70 6f 73 65 | 64 20 2a 2f 0a 58 7b 0a |ranspose|d */.X{.|
|00003110| 58 09 72 65 67 69 73 74 | 65 72 20 56 41 4c 55 45 |X.regist|er VALUE|
|00003120| 20 2a 76 31 2c 20 2a 76 | 32 3b 09 2f 2a 20 63 75 | *v1, *v|2;./* cu|
|00003130| 72 72 65 6e 74 20 76 61 | 6c 75 65 73 20 2a 2f 0a |rrent va|lues */.|
|00003140| 58 09 6c 6f 6e 67 20 72 | 6f 77 73 2c 20 63 6f 6c |X.long r|ows, col|
|00003150| 73 3b 09 09 2f 2a 20 72 | 6f 77 73 20 61 6e 64 20 |s;../* r|ows and |
|00003160| 63 6f 6c 75 6d 6e 73 20 | 69 6e 20 6f 6c 64 20 6d |columns |in old m|
|00003170| 61 74 72 69 78 20 2a 2f | 0a 58 09 6c 6f 6e 67 20 |atrix */|.X.long |
|00003180| 72 6f 77 2c 20 63 6f 6c | 3b 09 09 09 2f 2a 20 63 |row, col|;.../* c|
|00003190| 75 72 72 65 6e 74 20 72 | 6f 77 20 61 6e 64 20 63 |urrent r|ow and c|
|000031a0| 6f 6c 75 6d 6e 20 2a 2f | 0a 58 09 4d 41 54 52 49 |olumn */|.X.MATRI|
|000031b0| 58 20 2a 72 65 73 3b 0a | 58 0a 58 09 69 66 20 28 |X *res;.|X.X.if (|
|000031c0| 6d 2d 3e 6d 5f 64 69 6d | 20 21 3d 20 32 29 0a 58 |m->m_dim| != 2).X|
|000031d0| 09 09 65 72 72 6f 72 28 | 22 4d 61 74 72 69 78 20 |..error(|"Matrix |
|000031e0| 64 69 6d 65 6e 73 69 6f | 6e 20 6d 75 73 74 20 62 |dimensio|n must b|
|000031f0| 65 20 74 77 6f 20 66 6f | 72 20 74 72 61 6e 73 70 |e two fo|r transp|
|00003200| 6f 73 65 22 29 3b 0a 58 | 09 72 65 73 20 3d 20 6d |ose");.X|.res = m|
|00003210| 61 74 61 6c 6c 6f 63 28 | 6d 2d 3e 6d 5f 73 69 7a |atalloc(|m->m_siz|
|00003220| 65 29 3b 0a 58 09 72 65 | 73 2d 3e 6d 5f 64 69 6d |e);.X.re|s->m_dim|
|00003230| 20 3d 20 32 3b 0a 58 09 | 72 65 73 2d 3e 6d 5f 6d | = 2;.X.|res->m_m|
|00003240| 69 6e 5b 30 5d 20 3d 20 | 6d 2d 3e 6d 5f 6d 69 6e |in[0] = |m->m_min|
|00003250| 5b 31 5d 3b 0a 58 09 72 | 65 73 2d 3e 6d 5f 6d 61 |[1];.X.r|es->m_ma|
|00003260| 78 5b 30 5d 20 3d 20 6d | 2d 3e 6d 5f 6d 61 78 5b |x[0] = m|->m_max[|
|00003270| 31 5d 3b 0a 58 09 72 65 | 73 2d 3e 6d 5f 6d 69 6e |1];.X.re|s->m_min|
|00003280| 5b 31 5d 20 3d 20 6d 2d | 3e 6d 5f 6d 69 6e 5b 30 |[1] = m-|>m_min[0|
|00003290| 5d 3b 0a 58 09 72 65 73 | 2d 3e 6d 5f 6d 61 78 5b |];.X.res|->m_max[|
|000032a0| 31 5d 20 3d 20 6d 2d 3e | 6d 5f 6d 61 78 5b 30 5d |1] = m->|m_max[0]|
|000032b0| 3b 0a 58 09 72 6f 77 73 | 20 3d 20 28 6d 2d 3e 6d |;.X.rows| = (m->m|
|000032c0| 5f 6d 61 78 5b 30 5d 20 | 2d 20 6d 2d 3e 6d 5f 6d |_max[0] |- m->m_m|
|000032d0| 69 6e 5b 30 5d 20 2b 20 | 31 29 3b 0a 58 09 63 6f |in[0] + |1);.X.co|
|000032e0| 6c 73 20 3d 20 28 6d 2d | 3e 6d 5f 6d 61 78 5b 31 |ls = (m-|>m_max[1|
|000032f0| 5d 20 2d 20 6d 2d 3e 6d | 5f 6d 69 6e 5b 31 5d 20 |] - m->m|_min[1] |
|00003300| 2b 20 31 29 3b 0a 58 09 | 76 31 20 3d 20 72 65 73 |+ 1);.X.|v1 = res|
|00003310| 2d 3e 6d 5f 74 61 62 6c | 65 3b 0a 58 09 66 6f 72 |->m_tabl|e;.X.for|
|00003320| 20 28 72 6f 77 20 3d 20 | 30 3b 20 72 6f 77 20 3c | (row = |0; row <|
|00003330| 20 72 6f 77 73 3b 20 72 | 6f 77 2b 2b 29 20 7b 0a | rows; r|ow++) {.|
|00003340| 58 09 09 76 32 20 3d 20 | 26 6d 2d 3e 6d 5f 74 61 |X..v2 = |&m->m_ta|
|00003350| 62 6c 65 5b 72 6f 77 5d | 3b 0a 58 09 09 66 6f 72 |ble[row]|;.X..for|
|00003360| 20 28 63 6f 6c 20 3d 20 | 30 3b 20 63 6f 6c 20 3c | (col = |0; col <|
|00003370| 20 63 6f 6c 73 3b 20 63 | 6f 6c 2b 2b 29 20 7b 0a | cols; c|ol++) {.|
|00003380| 58 09 09 09 63 6f 70 79 | 76 61 6c 75 65 28 76 32 |X...copy|value(v2|
|00003390| 2c 20 76 31 29 3b 0a 58 | 09 09 09 76 31 2b 2b 3b |, v1);.X|...v1++;|
|000033a0| 0a 58 09 09 09 76 32 20 | 2b 3d 20 72 6f 77 73 3b |.X...v2 |+= rows;|
|000033b0| 0a 58 09 09 7d 0a 58 09 | 7d 0a 58 09 72 65 74 75 |.X..}.X.|}.X.retu|
|000033c0| 72 6e 20 72 65 73 3b 0a | 58 7d 0a 58 0a 58 0a 58 |rn res;.|X}.X.X.X|
|000033d0| 2f 2a 0a 58 20 2a 20 50 | 72 6f 64 75 63 65 20 61 |/*.X * P|roduce a|
|000033e0| 20 6d 61 74 72 69 78 20 | 77 69 74 68 20 76 61 6c | matrix |with val|
|000033f0| 75 65 73 20 61 6c 6c 20 | 6f 66 20 77 68 69 63 68 |ues all |of which|
|00003400| 20 61 72 65 20 63 6f 6e | 6a 75 67 61 74 65 64 2e | are con|jugated.|
|00003410| 0a 58 20 2a 2f 0a 58 4d | 41 54 52 49 58 20 2a 0a |.X */.XM|ATRIX *.|
|00003420| 58 6d 61 74 63 6f 6e 6a | 28 6d 29 0a 58 09 4d 41 |Xmatconj|(m).X.MA|
|00003430| 54 52 49 58 20 2a 6d 3b | 0a 58 7b 0a 58 09 72 65 |TRIX *m;|.X{.X.re|
|00003440| 67 69 73 74 65 72 20 56 | 41 4c 55 45 20 2a 76 61 |gister V|ALUE *va|
|00003450| 6c 2c 20 2a 76 72 65 73 | 3b 0a 58 09 6c 6f 6e 67 |l, *vres|;.X.long|
|00003460| 20 69 6e 64 65 78 3b 0a | 58 09 4d 41 54 52 49 58 | index;.|X.MATRIX|
|00003470| 20 2a 72 65 73 3b 0a 58 | 0a 58 09 72 65 73 20 3d | *res;.X|.X.res =|
|00003480| 20 6d 61 74 61 6c 6c 6f | 63 28 6d 2d 3e 6d 5f 73 | matallo|c(m->m_s|
|00003490| 69 7a 65 29 3b 0a 58 09 | 2a 72 65 73 20 3d 20 2a |ize);.X.|*res = *|
|000034a0| 6d 3b 0a 58 09 76 61 6c | 20 3d 20 6d 2d 3e 6d 5f |m;.X.val| = m->m_|
|000034b0| 74 61 62 6c 65 3b 0a 58 | 09 76 72 65 73 20 3d 20 |table;.X|.vres = |
|000034c0| 72 65 73 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 |res->m_t|able;.X.|
|000034d0| 66 6f 72 20 28 69 6e 64 | 65 78 20 3d 20 6d 2d 3e |for (ind|ex = m->|
|000034e0| 6d 5f 73 69 7a 65 3b 20 | 69 6e 64 65 78 20 3e 20 |m_size; |index > |
|000034f0| 30 3b 20 69 6e 64 65 78 | 2d 2d 29 0a 58 09 09 63 |0; index|--).X..c|
|00003500| 6f 6e 6a 76 61 6c 75 65 | 28 76 61 6c 2b 2b 2c 20 |onjvalue|(val++, |
|00003510| 76 72 65 73 2b 2b 29 3b | 0a 58 09 72 65 74 75 72 |vres++);|.X.retur|
|00003520| 6e 20 72 65 73 3b 0a 58 | 7d 0a 58 0a 58 0a 58 2f |n res;.X|}.X.X.X/|
|00003530| 2a 0a 58 20 2a 20 50 72 | 6f 64 75 63 65 20 61 20 |*.X * Pr|oduce a |
|00003540| 6d 61 74 72 69 78 20 77 | 69 74 68 20 76 61 6c 75 |matrix w|ith valu|
|00003550| 65 73 20 61 6c 6c 20 6f | 66 20 77 68 69 63 68 20 |es all o|f which |
|00003560| 68 61 76 65 20 62 65 65 | 6e 20 72 6f 75 6e 64 65 |have bee|n rounde|
|00003570| 64 20 74 6f 20 74 68 65 | 0a 58 20 2a 20 73 70 65 |d to the|.X * spe|
|00003580| 63 69 66 69 65 64 20 6e | 75 6d 62 65 72 20 6f 66 |cified n|umber of|
|00003590| 20 64 65 63 69 6d 61 6c | 20 70 6c 61 63 65 73 2e | decimal| places.|
|000035a0| 0a 58 20 2a 2f 0a 58 4d | 41 54 52 49 58 20 2a 0a |.X */.XM|ATRIX *.|
|000035b0| 58 6d 61 74 72 6f 75 6e | 64 28 6d 2c 20 70 6c 61 |Xmatroun|d(m, pla|
|000035c0| 63 65 73 29 0a 58 09 4d | 41 54 52 49 58 20 2a 6d |ces).X.M|ATRIX *m|
|000035d0| 3b 0a 58 09 6c 6f 6e 67 | 20 70 6c 61 63 65 73 3b |;.X.long| places;|
|000035e0| 0a 58 7b 0a 58 09 72 65 | 67 69 73 74 65 72 20 56 |.X{.X.re|gister V|
|000035f0| 41 4c 55 45 20 2a 76 61 | 6c 2c 20 2a 76 72 65 73 |ALUE *va|l, *vres|
|00003600| 3b 0a 58 09 56 41 4c 55 | 45 20 74 6d 70 3b 0a 58 |;.X.VALU|E tmp;.X|
|00003610| 09 6c 6f 6e 67 20 69 6e | 64 65 78 3b 0a 58 09 4d |.long in|dex;.X.M|
|00003620| 41 54 52 49 58 20 2a 72 | 65 73 3b 0a 58 0a 58 09 |ATRIX *r|es;.X.X.|
|00003630| 69 66 20 28 70 6c 61 63 | 65 73 20 3c 20 30 29 0a |if (plac|es < 0).|
|00003640| 58 09 09 65 72 72 6f 72 | 28 22 4e 65 67 61 74 69 |X..error|("Negati|
|00003650| 76 65 20 6e 75 6d 62 65 | 72 20 6f 66 20 70 6c 61 |ve numbe|r of pla|
|00003660| 63 65 73 20 66 6f 72 20 | 6d 61 74 72 6f 75 6e 64 |ces for |matround|
|00003670| 22 29 3b 0a 58 09 72 65 | 73 20 3d 20 6d 61 74 61 |");.X.re|s = mata|
|00003680| 6c 6c 6f 63 28 6d 2d 3e | 6d 5f 73 69 7a 65 29 3b |lloc(m->|m_size);|
|00003690| 0a 58 09 2a 72 65 73 20 | 3d 20 2a 6d 3b 0a 58 09 |.X.*res |= *m;.X.|
|000036a0| 74 6d 70 2e 76 5f 74 79 | 70 65 20 3d 20 56 5f 49 |tmp.v_ty|pe = V_I|
|000036b0| 4e 54 3b 0a 58 09 74 6d | 70 2e 76 5f 69 6e 74 20 |NT;.X.tm|p.v_int |
|000036c0| 3d 20 70 6c 61 63 65 73 | 3b 0a 58 09 76 61 6c 20 |= places|;.X.val |
|000036d0| 3d 20 6d 2d 3e 6d 5f 74 | 61 62 6c 65 3b 0a 58 09 |= m->m_t|able;.X.|
|000036e0| 76 72 65 73 20 3d 20 72 | 65 73 2d 3e 6d 5f 74 61 |vres = r|es->m_ta|
|000036f0| 62 6c 65 3b 0a 58 09 66 | 6f 72 20 28 69 6e 64 65 |ble;.X.f|or (inde|
|00003700| 78 20 3d 20 6d 2d 3e 6d | 5f 73 69 7a 65 3b 20 69 |x = m->m|_size; i|
|00003710| 6e 64 65 78 20 3e 20 30 | 3b 20 69 6e 64 65 78 2d |ndex > 0|; index-|
|00003720| 2d 29 0a 58 09 09 72 6f | 75 6e 64 76 61 6c 75 65 |-).X..ro|undvalue|
|00003730| 28 76 61 6c 2b 2b 2c 20 | 26 74 6d 70 2c 20 76 72 |(val++, |&tmp, vr|
|00003740| 65 73 2b 2b 29 3b 0a 58 | 09 72 65 74 75 72 6e 20 |es++);.X|.return |
|00003750| 72 65 73 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 0a |res;.X}.|X.X.X/*.|
|00003760| 58 20 2a 20 50 72 6f 64 | 75 63 65 20 61 20 6d 61 |X * Prod|uce a ma|
|00003770| 74 72 69 78 20 77 69 74 | 68 20 76 61 6c 75 65 73 |trix wit|h values|
|00003780| 20 61 6c 6c 20 6f 66 20 | 77 68 69 63 68 20 68 61 | all of |which ha|
|00003790| 76 65 20 62 65 65 6e 20 | 72 6f 75 6e 64 65 64 20 |ve been |rounded |
|000037a0| 74 6f 20 74 68 65 0a 58 | 20 2a 20 73 70 65 63 69 |to the.X| * speci|
|000037b0| 66 69 65 64 20 6e 75 6d | 62 65 72 20 6f 66 20 62 |fied num|ber of b|
|000037c0| 69 6e 61 72 79 20 70 6c | 61 63 65 73 2e 0a 58 20 |inary pl|aces..X |
|000037d0| 2a 2f 0a 58 4d 41 54 52 | 49 58 20 2a 0a 58 6d 61 |*/.XMATR|IX *.Xma|
|000037e0| 74 62 72 6f 75 6e 64 28 | 6d 2c 20 70 6c 61 63 65 |tbround(|m, place|
|000037f0| 73 29 0a 58 09 4d 41 54 | 52 49 58 20 2a 6d 3b 0a |s).X.MAT|RIX *m;.|
|00003800| 58 09 6c 6f 6e 67 20 70 | 6c 61 63 65 73 3b 0a 58 |X.long p|laces;.X|
|00003810| 7b 0a 58 09 72 65 67 69 | 73 74 65 72 20 56 41 4c |{.X.regi|ster VAL|
|00003820| 55 45 20 2a 76 61 6c 2c | 20 2a 76 72 65 73 3b 0a |UE *val,| *vres;.|
|00003830| 58 09 56 41 4c 55 45 20 | 74 6d 70 3b 0a 58 09 6c |X.VALUE |tmp;.X.l|
|00003840| 6f 6e 67 20 69 6e 64 65 | 78 3b 0a 58 09 4d 41 54 |ong inde|x;.X.MAT|
|00003850| 52 49 58 20 2a 72 65 73 | 3b 0a 58 0a 58 09 69 66 |RIX *res|;.X.X.if|
|00003860| 20 28 70 6c 61 63 65 73 | 20 3c 20 30 29 0a 58 09 | (places| < 0).X.|
|00003870| 09 65 72 72 6f 72 28 22 | 4e 65 67 61 74 69 76 65 |.error("|Negative|
|00003880| 20 6e 75 6d 62 65 72 20 | 6f 66 20 70 6c 61 63 65 | number |of place|
|00003890| 73 20 66 6f 72 20 6d 61 | 74 62 72 6f 75 6e 64 22 |s for ma|tbround"|
|000038a0| 29 3b 0a 58 09 72 65 73 | 20 3d 20 6d 61 74 61 6c |);.X.res| = matal|
|000038b0| 6c 6f 63 28 6d 2d 3e 6d | 5f 73 69 7a 65 29 3b 0a |loc(m->m|_size);.|
|000038c0| 58 09 2a 72 65 73 20 3d | 20 2a 6d 3b 0a 58 09 74 |X.*res =| *m;.X.t|
|000038d0| 6d 70 2e 76 5f 74 79 70 | 65 20 3d 20 56 5f 49 4e |mp.v_typ|e = V_IN|
|000038e0| 54 3b 0a 58 09 74 6d 70 | 2e 76 5f 69 6e 74 20 3d |T;.X.tmp|.v_int =|
|000038f0| 20 70 6c 61 63 65 73 3b | 0a 58 09 76 61 6c 20 3d | places;|.X.val =|
|00003900| 20 6d 2d 3e 6d 5f 74 61 | 62 6c 65 3b 0a 58 09 76 | m->m_ta|ble;.X.v|
|00003910| 72 65 73 20 3d 20 72 65 | 73 2d 3e 6d 5f 74 61 62 |res = re|s->m_tab|
|00003920| 6c 65 3b 0a 58 09 66 6f | 72 20 28 69 6e 64 65 78 |le;.X.fo|r (index|
|00003930| 20 3d 20 6d 2d 3e 6d 5f | 73 69 7a 65 3b 20 69 6e | = m->m_|size; in|
|00003940| 64 65 78 20 3e 20 30 3b | 20 69 6e 64 65 78 2d 2d |dex > 0;| index--|
|00003950| 29 0a 58 09 09 62 72 6f | 75 6e 64 76 61 6c 75 65 |).X..bro|undvalue|
|00003960| 28 76 61 6c 2b 2b 2c 20 | 26 74 6d 70 2c 20 76 72 |(val++, |&tmp, vr|
|00003970| 65 73 2b 2b 29 3b 0a 58 | 09 72 65 74 75 72 6e 20 |es++);.X|.return |
|00003980| 72 65 73 3b 0a 58 7d 0a | 58 0a 58 0a 58 2f 2a 0a |res;.X}.|X.X.X/*.|
|00003990| 58 20 2a 20 50 72 6f 64 | 75 63 65 20 61 20 6d 61 |X * Prod|uce a ma|
|000039a0| 74 72 69 78 20 77 69 74 | 68 20 76 61 6c 75 65 73 |trix wit|h values|
|000039b0| 20 61 6c 6c 20 6f 66 20 | 77 68 69 63 68 20 68 61 | all of |which ha|
|000039c0| 76 65 20 62 65 65 6e 20 | 74 72 75 6e 63 61 74 65 |ve been |truncate|
|000039d0| 64 20 74 6f 20 69 6e 74 | 65 67 65 72 73 2e 0a 58 |d to int|egers..X|
|000039e0| 20 2a 2f 0a 58 4d 41 54 | 52 49 58 20 2a 0a 58 6d | */.XMAT|RIX *.Xm|
|000039f0| 61 74 69 6e 74 28 6d 29 | 0a 58 09 4d 41 54 52 49 |atint(m)|.X.MATRI|
|00003a00| 58 20 2a 6d 3b 0a 58 7b | 0a 58 09 72 65 67 69 73 |X *m;.X{|.X.regis|
|00003a10| 74 65 72 20 56 41 4c 55 | 45 20 2a 76 61 6c 2c 20 |ter VALU|E *val, |
|00003a20| 2a 76 72 65 73 3b 0a 58 | 09 6c 6f 6e 67 20 69 6e |*vres;.X|.long in|
|00003a30| 64 65 78 3b 0a 58 09 4d | 41 54 52 49 58 20 2a 72 |dex;.X.M|ATRIX *r|
|00003a40| 65 73 3b 0a 58 0a 58 09 | 72 65 73 20 3d 20 6d 61 |es;.X.X.|res = ma|
|00003a50| 74 61 6c 6c 6f 63 28 6d | 2d 3e 6d 5f 73 69 7a 65 |talloc(m|->m_size|
|00003a60| 29 3b 0a 58 09 2a 72 65 | 73 20 3d 20 2a 6d 3b 0a |);.X.*re|s = *m;.|
|00003a70| 58 09 76 61 6c 20 3d 20 | 6d 2d 3e 6d 5f 74 61 62 |X.val = |m->m_tab|
|00003a80| 6c 65 3b 0a 58 09 76 72 | 65 73 20 3d 20 72 65 73 |le;.X.vr|es = res|
|00003a90| 2d 3e 6d 5f 74 61 62 6c | 65 3b 0a 58 09 66 6f 72 |->m_tabl|e;.X.for|
|00003aa0| 20 28 69 6e 64 65 78 20 | 3d 20 6d 2d 3e 6d 5f 73 | (index |= m->m_s|
|00003ab0| 69 7a 65 3b 20 69 6e 64 | 65 78 20 3e 20 30 3b 20 |ize; ind|ex > 0; |
|00003ac0| 69 6e 64 65 78 2d 2d 29 | 0a 58 09 09 69 6e 74 76 |index--)|.X..intv|
|00003ad0| 61 6c 75 65 28 76 61 6c | 2b 2b 2c 20 76 72 65 73 |alue(val|++, vres|
|00003ae0| 2b 2b 29 3b 0a 58 09 72 | 65 74 75 72 6e 20 72 65 |++);.X.r|eturn re|
|00003af0| 73 3b 0a 58 7d 0a 58 0a | 58 0a 58 2f 2a 0a 58 20 |s;.X}.X.|X.X/*.X |
|00003b00| 2a 20 50 72 6f 64 75 63 | 65 20 61 20 6d 61 74 72 |* Produc|e a matr|
|00003b10| 69 78 20 77 69 74 68 20 | 76 61 6c 75 65 73 20 61 |ix with |values a|
|00003b20| 6c 6c 20 6f 66 20 77 68 | 69 63 68 20 68 61 76 65 |ll of wh|ich have|
|00003b30| 20 6f 6e 6c 79 20 74 68 | 65 20 66 72 61 63 74 69 | only th|e fracti|
|00003b40| 6f 6e 20 70 61 72 74 20 | 6c 65 66 74 2e 0a 58 20 |on part |left..X |
|00003b50| 2a 2f 0a 58 4d 41 54 52 | 49 58 20 2a 0a 58 6d 61 |*/.XMATR|IX *.Xma|
|00003b60| 74 66 72 61 63 28 6d 29 | 0a 58 09 4d 41 54 52 49 |tfrac(m)|.X.MATRI|
|00003b70| 58 20 2a 6d 3b 0a 58 7b | 0a 58 09 72 65 67 69 73 |X *m;.X{|.X.regis|
|00003b80| 74 65 72 20 56 41 4c 55 | 45 20 2a 76 61 6c 2c 20 |ter VALU|E *val, |
|00003b90| 2a 76 72 65 73 3b 0a 58 | 09 6c 6f 6e 67 20 69 6e |*vres;.X|.long in|
|00003ba0| 64 65 78 3b 0a 58 09 4d | 41 54 52 49 58 20 2a 72 |dex;.X.M|ATRIX *r|
|00003bb0| 65 73 3b 0a 58 0a 58 09 | 72 65 73 20 3d 20 6d 61 |es;.X.X.|res = ma|
|00003bc0| 74 61 6c 6c 6f 63 28 6d | 2d 3e 6d 5f 73 69 7a 65 |talloc(m|->m_size|
|00003bd0| 29 3b 0a 58 09 2a 72 65 | 73 20 3d 20 2a 6d 3b 0a |);.X.*re|s = *m;.|
|00003be0| 58 09 76 61 6c 20 3d 20 | 6d 2d 3e 6d 5f 74 61 62 |X.val = |m->m_tab|
|00003bf0| 6c 65 3b 0a 58 09 76 72 | 65 73 20 3d 20 72 65 73 |le;.X.vr|es = res|
|00003c00| 2d 3e 6d 5f 74 61 62 6c | 65 3b 0a 58 09 66 6f 72 |->m_tabl|e;.X.for|
|00003c10| 20 28 69 6e 64 65 78 20 | 3d 20 6d 2d 3e 6d 5f 73 | (index |= m->m_s|
|00003c20| 69 7a 65 3b 20 69 6e 64 | 65 78 20 3e 20 30 3b 20 |ize; ind|ex > 0; |
|00003c30| 69 6e 64 65 78 2d 2d 29 | 0a 58 09 09 66 72 61 63 |index--)|.X..frac|
|00003c40| 76 61 6c 75 65 28 76 61 | 6c 2b 2b 2c 20 76 72 65 |value(va|l++, vre|
|00003c50| 73 2b 2b 29 3b 0a 58 09 | 72 65 74 75 72 6e 20 72 |s++);.X.|return r|
|00003c60| 65 73 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 0a 58 |es;.X}.X|.X.X/*.X|
|00003c70| 20 2a 20 53 65 61 72 63 | 68 20 61 20 6d 61 74 72 | * Searc|h a matr|
|00003c80| 69 78 20 66 6f 72 20 74 | 68 65 20 73 70 65 63 69 |ix for t|he speci|
|00003c90| 66 69 65 64 20 76 61 6c | 75 65 2c 20 73 74 61 72 |fied val|ue, star|
|00003ca0| 74 69 6e 67 20 77 69 74 | 68 20 74 68 65 20 73 70 |ting wit|h the sp|
|00003cb0| 65 63 69 66 69 65 64 20 | 69 6e 64 65 78 2e 0a 58 |ecified |index..X|
|00003cc0| 20 2a 20 52 65 74 75 72 | 6e 73 20 74 68 65 20 69 | * Retur|ns the i|
|00003cd0| 6e 64 65 78 20 6f 66 20 | 74 68 65 20 66 6f 75 6e |ndex of |the foun|
|00003ce0| 64 20 76 61 6c 75 65 2c | 20 6f 72 20 2d 31 20 69 |d value,| or -1 i|
|00003cf0| 66 20 74 68 65 20 76 61 | 6c 75 65 20 77 61 73 20 |f the va|lue was |
|00003d00| 6e 6f 74 20 66 6f 75 6e | 64 2e 0a 58 20 2a 2f 0a |not foun|d..X */.|
|00003d10| 58 6c 6f 6e 67 0a 58 6d | 61 74 73 65 61 72 63 68 |Xlong.Xm|atsearch|
|00003d20| 28 6d 2c 20 76 70 2c 20 | 69 6e 64 65 78 29 0a 58 |(m, vp, |index).X|
|00003d30| 09 4d 41 54 52 49 58 20 | 2a 6d 3b 0a 58 09 56 41 |.MATRIX |*m;.X.VA|
|00003d40| 4c 55 45 20 2a 76 70 3b | 0a 58 09 6c 6f 6e 67 20 |LUE *vp;|.X.long |
|00003d50| 69 6e 64 65 78 3b 0a 58 | 7b 0a 58 09 72 65 67 69 |index;.X|{.X.regi|
|00003d60| 73 74 65 72 20 56 41 4c | 55 45 20 2a 76 61 6c 3b |ster VAL|UE *val;|
|00003d70| 0a 58 0a 58 09 69 66 20 | 28 69 6e 64 65 78 20 3c |.X.X.if |(index <|
|00003d80| 20 30 29 0a 58 09 09 69 | 6e 64 65 78 20 3d 20 30 | 0).X..i|ndex = 0|
|00003d90| 3b 0a 58 09 76 61 6c 20 | 3d 20 26 6d 2d 3e 6d 5f |;.X.val |= &m->m_|
|00003da0| 74 61 62 6c 65 5b 69 6e | 64 65 78 5d 3b 0a 58 09 |table[in|dex];.X.|
|00003db0| 77 68 69 6c 65 20 28 69 | 6e 64 65 78 20 3c 20 6d |while (i|ndex < m|
|00003dc0| 2d 3e 6d 5f 73 69 7a 65 | 29 20 7b 0a 58 09 09 69 |->m_size|) {.X..i|
|00003dd0| 66 20 28 21 63 6f 6d 70 | 61 72 65 76 61 6c 75 65 |f (!comp|arevalue|
|00003de0| 28 76 70 2c 20 76 61 6c | 29 29 0a 58 09 09 09 72 |(vp, val|)).X...r|
|00003df0| 65 74 75 72 6e 20 69 6e | 64 65 78 3b 0a 58 09 09 |eturn in|dex;.X..|
|00003e00| 69 6e 64 65 78 2b 2b 3b | 0a 58 09 09 76 61 6c 2b |index++;|.X..val+|
|00003e10| 2b 3b 0a 58 09 7d 0a 58 | 09 72 65 74 75 72 6e 20 |+;.X.}.X|.return |
|00003e20| 2d 31 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 0a 58 |-1;.X}.X|.X.X/*.X|
|00003e30| 20 2a 20 53 65 61 72 63 | 68 20 61 20 6d 61 74 72 | * Searc|h a matr|
|00003e40| 69 78 20 62 61 63 6b 77 | 61 72 64 73 20 66 6f 72 |ix backw|ards for|
|00003e50| 20 74 68 65 20 73 70 65 | 63 69 66 69 65 64 20 76 | the spe|cified v|
|00003e60| 61 6c 75 65 2c 20 73 74 | 61 72 74 69 6e 67 20 77 |alue, st|arting w|
|00003e70| 69 74 68 20 74 68 65 0a | 58 20 2a 20 73 70 65 63 |ith the.|X * spec|
|00003e80| 69 66 69 65 64 20 69 6e | 64 65 78 2e 20 20 52 65 |ified in|dex.  Re|
|00003e90| 74 75 72 6e 73 20 74 68 | 65 20 69 6e 64 65 78 20 |turns th|e index |
|00003ea0| 6f 66 20 74 68 65 20 66 | 6f 75 6e 64 20 76 61 6c |of the f|ound val|
|00003eb0| 75 65 2c 20 6f 72 20 2d | 31 20 69 66 20 74 68 65 |ue, or -|1 if the|
|00003ec0| 0a 58 20 2a 20 76 61 6c | 75 65 20 77 61 73 20 6e |.X * val|ue was n|
|00003ed0| 6f 74 20 66 6f 75 6e 64 | 2e 0a 58 20 2a 2f 0a 58 |ot found|..X */.X|
|00003ee0| 6c 6f 6e 67 0a 58 6d 61 | 74 72 73 65 61 72 63 68 |long.Xma|trsearch|
|00003ef0| 28 6d 2c 20 76 70 2c 20 | 69 6e 64 65 78 29 0a 58 |(m, vp, |index).X|
|00003f00| 09 4d 41 54 52 49 58 20 | 2a 6d 3b 0a 58 09 56 41 |.MATRIX |*m;.X.VA|
|00003f10| 4c 55 45 20 2a 76 70 3b | 0a 58 09 6c 6f 6e 67 20 |LUE *vp;|.X.long |
|00003f20| 69 6e 64 65 78 3b 0a 58 | 7b 0a 58 09 72 65 67 69 |index;.X|{.X.regi|
|00003f30| 73 74 65 72 20 56 41 4c | 55 45 20 2a 76 61 6c 3b |ster VAL|UE *val;|
|00003f40| 0a 58 0a 58 09 69 66 20 | 28 69 6e 64 65 78 20 3e |.X.X.if |(index >|
|00003f50| 3d 20 6d 2d 3e 6d 5f 73 | 69 7a 65 29 0a 58 09 09 |= m->m_s|ize).X..|
|00003f60| 69 6e 64 65 78 20 3d 20 | 6d 2d 3e 6d 5f 73 69 7a |index = |m->m_siz|
|00003f70| 65 20 2d 20 31 3b 0a 58 | 09 76 61 6c 20 3d 20 26 |e - 1;.X|.val = &|
|00003f80| 6d 2d 3e 6d 5f 74 61 62 | 6c 65 5b 69 6e 64 65 78 |m->m_tab|le[index|
|00003f90| 5d 3b 0a 58 09 77 68 69 | 6c 65 20 28 69 6e 64 65 |];.X.whi|le (inde|
|00003fa0| 78 20 3e 3d 20 30 29 20 | 7b 0a 58 09 09 69 66 20 |x >= 0) |{.X..if |
|00003fb0| 28 21 63 6f 6d 70 61 72 | 65 76 61 6c 75 65 28 76 |(!compar|evalue(v|
|00003fc0| 70 2c 20 76 61 6c 29 29 | 0a 58 09 09 09 72 65 74 |p, val))|.X...ret|
|00003fd0| 75 72 6e 20 69 6e 64 65 | 78 3b 0a 58 09 09 69 6e |urn inde|x;.X..in|
|00003fe0| 64 65 78 2d 2d 3b 0a 58 | 09 09 76 61 6c 2d 2d 3b |dex--;.X|..val--;|
|00003ff0| 0a 58 09 7d 0a 58 09 72 | 65 74 75 72 6e 20 2d 31 |.X.}.X.r|eturn -1|
|00004000| 3b 0a 58 7d 0a 58 0a 58 | 0a 58 2f 2a 0a 58 20 2a |;.X}.X.X|.X/*.X *|
|00004010| 20 46 69 6c 6c 20 61 6c | 6c 20 6f 66 20 74 68 65 | Fill al|l of the|
|00004020| 20 65 6c 65 6d 65 6e 74 | 73 20 6f 66 20 61 20 6d | element|s of a m|
|00004030| 61 74 72 69 78 20 77 69 | 74 68 20 6f 6e 65 20 6f |atrix wi|th one o|
|00004040| 66 20 74 77 6f 20 73 70 | 65 63 69 66 69 65 64 20 |f two sp|ecified |
|00004050| 76 61 6c 75 65 73 2e 0a | 58 20 2a 20 41 6c 6c 20 |values..|X * All |
|00004060| 65 6e 74 72 69 65 73 20 | 61 72 65 20 66 69 6c 6c |entries |are fill|
|00004070| 65 64 20 77 69 74 68 20 | 74 68 65 20 66 69 72 73 |ed with |the firs|
|00004080| 74 20 73 70 65 63 69 66 | 69 65 64 20 76 61 6c 75 |t specif|ied valu|
|00004090| 65 2c 20 65 78 63 65 70 | 74 20 74 68 61 74 20 69 |e, excep|t that i|
|000040a0| 66 0a 58 20 2a 20 74 68 | 65 20 6d 61 74 72 69 78 |f.X * th|e matrix|
|000040b0| 20 69 73 20 73 71 75 61 | 72 65 20 61 6e 64 20 74 | is squa|re and t|
|000040c0| 68 65 20 73 65 63 6f 6e | 64 20 76 61 6c 75 65 20 |he secon|d value |
|000040d0| 70 6f 69 6e 74 65 72 20 | 69 73 20 6e 6f 6e 2d 4e |pointer |is non-N|
|000040e0| 55 4c 4c 2c 20 74 68 65 | 6e 0a 58 20 2a 20 61 6c |ULL, the|n.X * al|
|000040f0| 6c 20 64 69 61 67 6f 6e | 61 6c 20 65 6e 74 72 69 |l diagon|al entri|
|00004100| 65 73 20 61 72 65 20 66 | 69 6c 6c 65 64 20 77 69 |es are f|illed wi|
|00004110| 74 68 20 74 68 65 20 73 | 65 63 6f 6e 64 20 76 61 |th the s|econd va|
|00004120| 6c 75 65 2e 20 20 54 68 | 69 73 20 72 6f 75 74 69 |lue.  Th|is routi|
|00004130| 6e 65 0a 58 20 2a 20 61 | 66 66 65 63 74 73 20 74 |ne.X * a|ffects t|
|00004140| 68 65 20 73 75 70 70 6c | 69 65 64 20 6d 61 74 72 |he suppl|ied matr|
|00004150| 69 78 20 64 69 72 65 63 | 74 6c 79 2c 20 61 6e 64 |ix direc|tly, and|
|00004160| 20 64 6f 65 73 6e 27 74 | 20 72 65 74 75 72 6e 20 | doesn't| return |
|00004170| 61 20 63 6f 70 79 2e 0a | 58 20 2a 2f 0a 58 76 6f |a copy..|X */.Xvo|
|00004180| 69 64 0a 58 6d 61 74 66 | 69 6c 6c 28 6d 2c 20 76 |id.Xmatf|ill(m, v|
|00004190| 31 2c 20 76 32 29 0a 58 | 09 4d 41 54 52 49 58 20 |1, v2).X|.MATRIX |
|000041a0| 2a 6d 3b 09 09 2f 2a 20 | 6d 61 74 72 69 78 20 74 |*m;../* |matrix t|
|000041b0| 6f 20 62 65 20 66 69 6c | 6c 65 64 20 2a 2f 0a 58 |o be fil|led */.X|
|000041c0| 09 56 41 4c 55 45 20 2a | 76 31 3b 09 09 2f 2a 20 |.VALUE *|v1;../* |
|000041d0| 76 61 6c 75 65 20 74 6f | 20 66 69 6c 6c 20 6d 6f |value to| fill mo|
|000041e0| 73 74 20 6f 66 20 6d 61 | 74 72 69 78 20 77 69 74 |st of ma|trix wit|
|000041f0| 68 20 2a 2f 0a 58 09 56 | 41 4c 55 45 20 2a 76 32 |h */.X.V|ALUE *v2|
|00004200| 3b 09 09 2f 2a 20 76 61 | 6c 75 65 20 66 6f 72 20 |;../* va|lue for |
|00004210| 64 69 61 67 6f 6e 61 6c | 20 65 6e 74 72 69 65 73 |diagonal| entries|
|00004220| 20 28 6f 72 20 4e 55 4c | 4c 29 20 2a 2f 0a 58 7b | (or NUL|L) */.X{|
|00004230| 0a 58 09 72 65 67 69 73 | 74 65 72 20 56 41 4c 55 |.X.regis|ter VALU|
|00004240| 45 20 2a 76 61 6c 3b 0a | 58 09 6c 6f 6e 67 20 72 |E *val;.|X.long r|
|00004250| 6f 77 2c 20 63 6f 6c 3b | 0a 58 09 6c 6f 6e 67 20 |ow, col;|.X.long |
|00004260| 72 6f 77 73 3b 0a 58 09 | 6c 6f 6e 67 20 69 6e 64 |rows;.X.|long ind|
|00004270| 65 78 3b 0a 58 0a 58 09 | 69 66 20 28 76 32 20 26 |ex;.X.X.|if (v2 &|
|00004280| 26 20 28 28 6d 2d 3e 6d | 5f 64 69 6d 20 21 3d 20 |& ((m->m|_dim != |
|00004290| 32 29 20 7c 7c 0a 58 09 | 09 28 28 6d 2d 3e 6d 5f |2) ||.X.|.((m->m_|
|000042a0| 6d 61 78 5b 30 5d 20 2d | 20 6d 2d 3e 6d 5f 6d 69 |max[0] -| m->m_mi|
|000042b0| 6e 5b 30 5d 29 20 21 3d | 20 28 6d 2d 3e 6d 5f 6d |n[0]) !=| (m->m_m|
|000042c0| 61 78 5b 31 5d 20 2d 20 | 6d 2d 3e 6d 5f 6d 69 6e |ax[1] - |m->m_min|
|000042d0| 5b 31 5d 29 29 29 29 0a | 58 09 09 09 65 72 72 6f |[1])))).|X...erro|
|000042e0| 72 28 22 46 69 6c 6c 69 | 6e 67 20 64 69 61 67 6f |r("Filli|ng diago|
|000042f0| 6e 61 6c 73 20 6f 66 20 | 6e 6f 6e 2d 73 71 75 61 |nals of |non-squa|
|00004300| 72 65 20 6d 61 74 72 69 | 78 22 29 3b 0a 58 09 76 |re matri|x");.X.v|
|00004310| 61 6c 20 3d 20 6d 2d 3e | 6d 5f 74 61 62 6c 65 3b |al = m->|m_table;|
|00004320| 0a 58 09 66 6f 72 20 28 | 69 6e 64 65 78 20 3d 20 |.X.for (|index = |
|00004330| 6d 2d 3e 6d 5f 73 69 7a | 65 3b 20 69 6e 64 65 78 |m->m_siz|e; index|
|00004340| 20 3e 20 30 3b 20 69 6e | 64 65 78 2d 2d 29 0a 58 | > 0; in|dex--).X|
|00004350| 09 09 66 72 65 65 76 61 | 6c 75 65 28 76 61 6c 2b |..freeva|lue(val+|
|00004360| 2b 29 3b 0a 58 09 76 61 | 6c 20 3d 20 6d 2d 3e 6d |+);.X.va|l = m->m|
|00004370| 5f 74 61 62 6c 65 3b 0a | 58 09 69 66 20 28 76 32 |_table;.|X.if (v2|
|00004380| 20 3d 3d 20 4e 55 4c 4c | 29 20 7b 0a 58 09 09 66 | == NULL|) {.X..f|
|00004390| 6f 72 20 28 69 6e 64 65 | 78 20 3d 20 6d 2d 3e 6d |or (inde|x = m->m|
|000043a0| 5f 73 69 7a 65 3b 20 69 | 6e 64 65 78 20 3e 20 30 |_size; i|ndex > 0|
|000043b0| 3b 20 69 6e 64 65 78 2d | 2d 29 0a 58 09 09 09 63 |; index-|-).X...c|
|000043c0| 6f 70 79 76 61 6c 75 65 | 28 76 31 2c 20 76 61 6c |opyvalue|(v1, val|
|000043d0| 2b 2b 29 3b 0a 58 09 09 | 72 65 74 75 72 6e 3b 0a |++);.X..|return;.|
|000043e0| 58 09 7d 0a 58 09 72 6f | 77 73 20 3d 20 6d 2d 3e |X.}.X.ro|ws = m->|
|000043f0| 6d 5f 6d 61 78 5b 30 5d | 20 2d 20 6d 2d 3e 6d 5f |m_max[0]| - m->m_|
|00004400| 6d 69 6e 5b 30 5d 20 2b | 20 31 3b 0a 58 09 66 6f |min[0] +| 1;.X.fo|
|00004410| 72 20 28 72 6f 77 20 3d | 20 30 3b 20 72 6f 77 20 |r (row =| 0; row |
|00004420| 3c 20 72 6f 77 73 3b 20 | 72 6f 77 2b 2b 29 20 7b |< rows; |row++) {|
|00004430| 0a 58 09 09 66 6f 72 20 | 28 63 6f 6c 20 3d 20 30 |.X..for |(col = 0|
|00004440| 3b 20 63 6f 6c 20 3c 20 | 72 6f 77 73 3b 20 63 6f |; col < |rows; co|
|00004450| 6c 2b 2b 29 20 7b 0a 58 | 09 09 09 63 6f 70 79 76 |l++) {.X|...copyv|
|00004460| 61 6c 75 65 28 28 28 72 | 6f 77 20 21 3d 20 63 6f |alue(((r|ow != co|
|00004470| 6c 29 20 3f 20 76 31 20 | 3a 20 76 32 29 2c 20 76 |l) ? v1 |: v2), v|
|00004480| 61 6c 2b 2b 29 3b 0a 58 | 09 09 7d 0a 58 09 7d 0a |al++);.X|..}.X.}.|
|00004490| 58 7d 0a 58 0a 58 0a 58 | 2f 2a 0a 58 20 2a 20 53 |X}.X.X.X|/*.X * S|
|000044a0| 65 74 20 61 20 63 6f 70 | 79 20 6f 66 20 61 20 73 |et a cop|y of a s|
|000044b0| 71 75 61 72 65 20 6d 61 | 74 72 69 78 20 74 6f 20 |quare ma|trix to |
|000044c0| 74 68 65 20 69 64 65 6e | 74 69 74 79 20 6d 61 74 |the iden|tity mat|
|000044d0| 72 69 78 2e 0a 58 20 2a | 2f 0a 58 73 74 61 74 69 |rix..X *|/.Xstati|
|000044e0| 63 20 4d 41 54 52 49 58 | 20 2a 0a 58 6d 61 74 69 |c MATRIX| *.Xmati|
|000044f0| 64 65 6e 74 28 6d 29 0a | 58 09 4d 41 54 52 49 58 |dent(m).|X.MATRIX|
|00004500| 20 2a 6d 3b 0a 58 7b 0a | 58 09 72 65 67 69 73 74 | *m;.X{.|X.regist|
|00004510| 65 72 20 56 41 4c 55 45 | 20 2a 76 61 6c 3b 09 2f |er VALUE| *val;./|
|00004520| 2a 20 63 75 72 72 65 6e | 74 20 76 61 6c 75 65 20 |* curren|t value |
|00004530| 2a 2f 0a 58 09 6c 6f 6e | 67 20 72 6f 77 2c 20 63 |*/.X.lon|g row, c|
|00004540| 6f 6c 3b 09 09 2f 2a 20 | 63 75 72 72 65 6e 74 20 |ol;../* |current |
|00004550| 72 6f 77 20 61 6e 64 20 | 63 6f 6c 75 6d 6e 20 2a |row and |column *|
|00004560| 2f 0a 58 09 6c 6f 6e 67 | 20 72 6f 77 73 3b 09 09 |/.X.long| rows;..|
|00004570| 2f 2a 20 6e 75 6d 62 65 | 72 20 6f 66 20 72 6f 77 |/* numbe|r of row|
|00004580| 73 20 2a 2f 0a 58 09 4d | 41 54 52 49 58 20 2a 72 |s */.X.M|ATRIX *r|
|00004590| 65 73 3b 09 09 2f 2a 20 | 72 65 73 75 6c 74 69 6e |es;../* |resultin|
|000045a0| 67 20 6d 61 74 72 69 78 | 20 2a 2f 0a 58 0a 58 09 |g matrix| */.X.X.|
|000045b0| 69 66 20 28 6d 2d 3e 6d | 5f 64 69 6d 20 21 3d 20 |if (m->m|_dim != |
|000045c0| 32 29 0a 58 09 09 65 72 | 72 6f 72 28 22 4d 61 74 |2).X..er|ror("Mat|
|000045d0| 72 69 78 20 64 69 6d 65 | 6e 73 69 6f 6e 20 6d 75 |rix dime|nsion mu|
|000045e0| 73 74 20 62 65 20 74 77 | 6f 20 66 6f 72 20 73 65 |st be tw|o for se|
|000045f0| 74 74 69 6e 67 20 74 6f | 20 69 64 65 6e 74 69 74 |tting to| identit|
|00004600| 79 22 29 3b 0a 58 09 69 | 66 20 28 28 6d 2d 3e 6d |y");.X.i|f ((m->m|
|00004610| 5f 6d 61 78 5b 30 5d 20 | 2d 20 6d 2d 3e 6d 5f 6d |_max[0] |- m->m_m|
|00004620| 69 6e 5b 30 5d 29 20 21 | 3d 20 28 6d 2d 3e 6d 5f |in[0]) !|= (m->m_|
|00004630| 6d 61 78 5b 31 5d 20 2d | 20 6d 2d 3e 6d 5f 6d 69 |max[1] -| m->m_mi|
|00004640| 6e 5b 31 5d 29 29 0a 58 | 09 09 65 72 72 6f 72 28 |n[1])).X|..error(|
|00004650| 22 4d 61 74 72 69 78 20 | 6d 75 73 74 20 62 65 20 |"Matrix |must be |
|00004660| 73 71 75 61 72 65 20 66 | 6f 72 20 73 65 74 74 69 |square f|or setti|
|00004670| 6e 67 20 74 6f 20 69 64 | 65 6e 74 69 74 79 22 29 |ng to id|entity")|
|00004680| 3b 0a 58 09 72 65 73 20 | 3d 20 6d 61 74 61 6c 6c |;.X.res |= matall|
|00004690| 6f 63 28 6d 2d 3e 6d 5f | 73 69 7a 65 29 3b 0a 58 |oc(m->m_|size);.X|
|000046a0| 09 2a 72 65 73 20 3d 20 | 2a 6d 3b 0a 58 09 76 61 |.*res = |*m;.X.va|
|000046b0| 6c 20 3d 20 72 65 73 2d | 3e 6d 5f 74 61 62 6c 65 |l = res-|>m_table|
|000046c0| 3b 0a 58 09 72 6f 77 73 | 20 3d 20 28 72 65 73 2d |;.X.rows| = (res-|
|000046d0| 3e 6d 5f 6d 61 78 5b 30 | 5d 20 2d 20 72 65 73 2d |>m_max[0|] - res-|
|000046e0| 3e 6d 5f 6d 69 6e 5b 30 | 5d 20 2b 20 31 29 3b 0a |>m_min[0|] + 1);.|
|000046f0| 58 09 66 6f 72 20 28 72 | 6f 77 20 3d 20 30 3b 20 |X.for (r|ow = 0; |
|00004700| 72 6f 77 20 3c 20 72 6f | 77 73 3b 20 72 6f 77 2b |row < ro|ws; row+|
|00004710| 2b 29 20 7b 0a 58 09 09 | 66 6f 72 20 28 63 6f 6c |+) {.X..|for (col|
|00004720| 20 3d 20 30 3b 20 63 6f | 6c 20 3c 20 72 6f 77 73 | = 0; co|l < rows|
|00004730| 3b 20 63 6f 6c 2b 2b 29 | 20 7b 0a 58 09 09 09 76 |; col++)| {.X...v|
|00004740| 61 6c 2d 3e 76 5f 74 79 | 70 65 20 3d 20 56 5f 4e |al->v_ty|pe = V_N|
|00004750| 55 4d 3b 0a 58 09 09 09 | 76 61 6c 2d 3e 76 5f 6e |UM;.X...|val->v_n|
|00004760| 75 6d 20 3d 20 28 28 72 | 6f 77 20 3d 3d 20 63 6f |um = ((r|ow == co|
|00004770| 6c 29 20 3f 20 71 6c 69 | 6e 6b 28 26 5f 71 6f 6e |l) ? qli|nk(&_qon|
|00004780| 65 5f 29 20 3a 20 71 6c | 69 6e 6b 28 26 5f 71 7a |e_) : ql|ink(&_qz|
|00004790| 65 72 6f 5f 29 29 3b 0a | 58 09 09 09 76 61 6c 2b |ero_));.|X...val+|
|000047a0| 2b 3b 0a 58 09 09 7d 0a | 58 09 7d 0a 58 09 72 65 |+;.X..}.|X.}.X.re|
|000047b0| 74 75 72 6e 20 72 65 73 | 3b 0a 58 7d 0a 58 0a 58 |turn res|;.X}.X.X|
|000047c0| 0a 58 2f 2a 0a 58 20 2a | 20 43 61 6c 63 75 6c 61 |.X/*.X *| Calcula|
|000047d0| 74 65 20 74 68 65 20 69 | 6e 76 65 72 73 65 20 6f |te the i|nverse o|
|000047e0| 66 20 61 20 6d 61 74 72 | 69 78 20 69 66 20 69 74 |f a matr|ix if it|
|000047f0| 20 65 78 69 73 74 73 2e | 0a 58 20 2a 20 54 68 69 | exists.|.X * Thi|
|00004800| 73 20 69 73 20 64 6f 6e | 65 20 62 79 20 75 73 69 |s is don|e by usi|
|00004810| 6e 67 20 74 72 61 6e 73 | 66 6f 72 6d 61 74 69 6f |ng trans|formatio|
|00004820| 6e 73 20 6f 6e 20 74 68 | 65 20 73 75 70 70 6c 69 |ns on th|e suppli|
|00004830| 65 64 20 6d 61 74 72 69 | 78 20 74 6f 20 63 6f 6e |ed matri|x to con|
|00004840| 76 65 72 74 0a 58 20 2a | 20 69 74 20 74 6f 20 74 |vert.X *| it to t|
|00004850| 68 65 20 69 64 65 6e 74 | 69 74 79 20 6d 61 74 72 |he ident|ity matr|
|00004860| 69 78 2c 20 61 6e 64 20 | 73 69 6d 75 6c 74 61 6e |ix, and |simultan|
|00004870| 65 6f 75 73 6c 79 20 61 | 70 70 6c 79 69 6e 67 20 |eously a|pplying |
|00004880| 74 68 65 20 73 61 6d 65 | 20 73 65 74 20 6f 66 0a |the same| set of.|
|00004890| 58 20 2a 20 74 72 61 6e | 73 66 6f 72 6d 61 74 69 |X * tran|sformati|
|000048a0| 6f 6e 73 20 74 6f 20 74 | 68 65 20 69 64 65 6e 74 |ons to t|he ident|
|000048b0| 69 74 79 20 6d 61 74 72 | 69 78 2e 0a 58 20 2a 2f |ity matr|ix..X */|
|000048c0| 0a 58 4d 41 54 52 49 58 | 20 2a 0a 58 6d 61 74 69 |.XMATRIX| *.Xmati|
|000048d0| 6e 76 28 6d 29 0a 58 09 | 4d 41 54 52 49 58 20 2a |nv(m).X.|MATRIX *|
|000048e0| 6d 3b 0a 58 7b 0a 58 09 | 4d 41 54 52 49 58 20 2a |m;.X{.X.|MATRIX *|
|000048f0| 72 65 73 3b 09 09 2f 2a | 20 6d 61 74 72 69 78 20 |res;../*| matrix |
|00004900| 74 6f 20 62 65 63 6f 6d | 65 20 74 68 65 20 69 6e |to becom|e the in|
|00004910| 76 65 72 73 65 20 2a 2f | 0a 58 09 6c 6f 6e 67 20 |verse */|.X.long |
|00004920| 72 6f 77 73 3b 09 09 2f | 2a 20 6e 75 6d 62 65 72 |rows;../|* number|
|00004930| 20 6f 66 20 72 6f 77 73 | 20 2a 2f 0a 58 09 6c 6f | of rows| */.X.lo|
|00004940| 6e 67 20 63 75 72 3b 09 | 09 2f 2a 20 63 75 72 72 |ng cur;.|./* curr|
|00004950| 65 6e 74 20 72 6f 77 20 | 62 65 69 6e 67 20 77 6f |ent row |being wo|
|00004960| 72 6b 65 64 20 6f 6e 20 | 2a 2f 0a 58 09 6c 6f 6e |rked on |*/.X.lon|
|00004970| 67 20 72 6f 77 2c 20 63 | 6f 6c 3b 09 09 2f 2a 20 |g row, c|ol;../* |
|00004980| 74 65 6d 70 20 72 6f 77 | 20 61 6e 64 20 63 6f 6c |temp row| and col|
|00004990| 75 6d 6e 20 76 61 6c 75 | 65 73 20 2a 2f 0a 58 09 |umn valu|es */.X.|
|000049a0| 56 41 4c 55 45 20 2a 76 | 61 6c 3b 09 09 2f 2a 20 |VALUE *v|al;../* |
|000049b0| 63 75 72 72 65 6e 74 20 | 76 61 6c 75 65 20 69 6e |current |value in|
|000049c0| 20 6d 61 74 72 69 78 2a | 2f 0a 58 09 56 41 4c 55 | matrix*|/.X.VALU|
|000049d0| 45 20 6d 75 6c 76 61 6c | 3b 09 09 2f 2a 20 76 61 |E mulval|;../* va|
|000049e0| 6c 75 65 20 74 6f 20 6d | 75 6c 74 69 70 6c 79 20 |lue to m|ultiply |
|000049f0| 72 6f 77 73 20 62 79 20 | 2a 2f 0a 58 09 56 41 4c |rows by |*/.X.VAL|
|00004a00| 55 45 20 74 6d 70 76 61 | 6c 3b 09 09 2f 2a 20 74 |UE tmpva|l;../* t|
|00004a10| 65 6d 70 6f 72 61 72 79 | 20 76 61 6c 75 65 20 2a |emporary| value *|
|00004a20| 2f 0a 58 0a 58 09 69 66 | 20 28 6d 2d 3e 6d 5f 64 |/.X.X.if| (m->m_d|
|00004a30| 69 6d 20 21 3d 20 32 29 | 0a 58 09 09 65 72 72 6f |im != 2)|.X..erro|
|00004a40| 72 28 22 4d 61 74 72 69 | 78 20 64 69 6d 65 6e 73 |r("Matri|x dimens|
|00004a50| 69 6f 6e 20 6d 75 73 74 | 20 62 65 20 74 77 6f 20 |ion must| be two |
|00004a60| 66 6f 72 20 69 6e 76 65 | 72 73 65 22 29 3b 0a 58 |for inve|rse");.X|
|00004a70| 09 69 66 20 28 28 6d 2d | 3e 6d 5f 6d 61 78 5b 30 |.if ((m-|>m_max[0|
|00004a80| 5d 20 2d 20 6d 2d 3e 6d | 5f 6d 69 6e 5b 30 5d 29 |] - m->m|_min[0])|
|00004a90| 20 21 3d 20 28 6d 2d 3e | 6d 5f 6d 61 78 5b 31 5d | != (m->|m_max[1]|
|00004aa0| 20 2d 20 6d 2d 3e 6d 5f | 6d 69 6e 5b 31 5d 29 29 | - m->m_|min[1]))|
|00004ab0| 0a 58 09 09 65 72 72 6f | 72 28 22 49 6e 76 65 72 |.X..erro|r("Inver|
|00004ac0| 74 69 6e 67 20 6e 6f 6e | 2d 73 71 75 61 72 65 20 |ting non|-square |
|00004ad0| 6d 61 74 72 69 78 22 29 | 3b 0a 58 09 2f 2a 0a 58 |matrix")|;.X./*.X|
|00004ae0| 09 20 2a 20 42 65 67 69 | 6e 20 62 79 20 63 72 65 |. * Begi|n by cre|
|00004af0| 61 74 69 6e 67 20 74 68 | 65 20 69 64 65 6e 74 69 |ating th|e identi|
|00004b00| 74 79 20 6d 61 74 72 69 | 78 20 77 69 74 68 20 74 |ty matri|x with t|
|00004b10| 68 65 20 73 61 6d 65 20 | 61 74 74 72 69 62 75 74 |he same |attribut|
|00004b20| 65 73 2e 0a 58 09 20 2a | 2f 0a 58 09 72 65 73 20 |es..X. *|/.X.res |
|00004b30| 3d 20 6d 61 74 61 6c 6c | 6f 63 28 6d 2d 3e 6d 5f |= matall|oc(m->m_|
|00004b40| 73 69 7a 65 29 3b 0a 58 | 09 2a 72 65 73 20 3d 20 |size);.X|.*res = |
|00004b50| 2a 6d 3b 0a 58 09 72 6f | 77 73 20 3d 20 28 6d 2d |*m;.X.ro|ws = (m-|
|00004b60| 3e 6d 5f 6d 61 78 5b 30 | 5d 20 2d 20 6d 2d 3e 6d |>m_max[0|] - m->m|
|00004b70| 5f 6d 69 6e 5b 30 5d 20 | 2b 20 31 29 3b 0a 58 09 |_min[0] |+ 1);.X.|
|00004b80| 76 61 6c 20 3d 20 72 65 | 73 2d 3e 6d 5f 74 61 62 |val = re|s->m_tab|
|00004b90| 6c 65 3b 0a 58 09 66 6f | 72 20 28 72 6f 77 20 3d |le;.X.fo|r (row =|
|00004ba0| 20 30 3b 20 72 6f 77 20 | 3c 20 72 6f 77 73 3b 20 | 0; row |< rows; |
|00004bb0| 72 6f 77 2b 2b 29 20 7b | 0a 58 09 09 66 6f 72 20 |row++) {|.X..for |
|00004bc0| 28 63 6f 6c 20 3d 20 30 | 3b 20 63 6f 6c 20 3c 20 |(col = 0|; col < |
|00004bd0| 72 6f 77 73 3b 20 63 6f | 6c 2b 2b 29 20 7b 0a 58 |rows; co|l++) {.X|
|00004be0| 09 09 09 69 66 20 28 72 | 6f 77 20 3d 3d 20 63 6f |...if (r|ow == co|
|00004bf0| 6c 29 0a 58 09 09 09 09 | 76 61 6c 2d 3e 76 5f 6e |l).X....|val->v_n|
|00004c00| 75 6d 20 3d 20 71 6c 69 | 6e 6b 28 26 5f 71 6f 6e |um = qli|nk(&_qon|
|00004c10| 65 5f 29 3b 0a 58 09 09 | 09 65 6c 73 65 0a 58 09 |e_);.X..|.else.X.|
|00004c20| 09 09 09 76 61 6c 2d 3e | 76 5f 6e 75 6d 20 3d 20 |...val->|v_num = |
|00004c30| 71 6c 69 6e 6b 28 26 5f | 71 7a 65 72 6f 5f 29 3b |qlink(&_|qzero_);|
|00004c40| 0a 58 09 09 09 76 61 6c | 2d 3e 76 5f 74 79 70 65 |.X...val|->v_type|
|00004c50| 20 3d 20 56 5f 4e 55 4d | 3b 0a 58 09 09 09 76 61 | = V_NUM|;.X...va|
|00004c60| 6c 2b 2b 3b 0a 58 09 09 | 7d 0a 58 09 7d 0a 58 09 |l++;.X..|}.X.}.X.|
|00004c70| 2f 2a 0a 58 09 20 2a 20 | 4e 6f 77 20 6c 6f 6f 70 |/*.X. * |Now loop|
|00004c80| 20 6f 76 65 72 20 65 61 | 63 68 20 72 6f 77 2c 20 | over ea|ch row, |
|00004c90| 61 6e 64 20 65 6c 69 6d | 69 6e 61 74 65 20 61 6c |and elim|inate al|
|00004ca0| 6c 20 65 6e 74 72 69 65 | 73 20 69 6e 20 74 68 65 |l entrie|s in the|
|00004cb0| 0a 58 09 20 2a 20 63 6f | 72 72 65 73 70 6f 6e 64 |.X. * co|rrespond|
|00004cc0| 69 6e 67 20 63 6f 6c 75 | 6d 6e 20 62 79 20 75 73 |ing colu|mn by us|
|00004cd0| 69 6e 67 20 72 6f 77 20 | 6f 70 65 72 61 74 69 6f |ing row |operatio|
|00004ce0| 6e 73 2e 20 20 44 6f 20 | 74 68 65 20 73 61 6d 65 |ns.  Do |the same|
|00004cf0| 0a 58 09 20 2a 20 6f 70 | 65 72 61 74 69 6f 6e 73 |.X. * op|erations|
|00004d00| 20 6f 6e 20 74 68 65 20 | 72 65 73 75 6c 74 69 6e | on the |resultin|
|00004d10| 67 20 6d 61 74 72 69 78 | 2e 20 20 43 6f 70 79 20 |g matrix|.  Copy |
|00004d20| 74 68 65 20 6f 72 69 67 | 69 6e 61 6c 20 6d 61 74 |the orig|inal mat|
|00004d30| 72 69 78 0a 58 09 20 2a | 20 73 6f 20 74 68 61 74 |rix.X. *| so that|
|00004d40| 20 77 65 20 64 6f 6e 27 | 74 20 64 65 73 74 72 6f | we don'|t destro|
|00004d50| 79 20 69 74 2e 0a 58 09 | 20 2a 2f 0a 58 09 6d 20 |y it..X.| */.X.m |
|00004d60| 3d 20 6d 61 74 63 6f 70 | 79 28 6d 29 3b 0a 58 09 |= matcop|y(m);.X.|
|00004d70| 66 6f 72 20 28 63 75 72 | 20 3d 20 30 3b 20 63 75 |for (cur| = 0; cu|
|00004d80| 72 20 3c 20 72 6f 77 73 | 3b 20 63 75 72 2b 2b 29 |r < rows|; cur++)|
|00004d90| 20 7b 0a 58 09 09 2f 2a | 0a 58 09 09 20 2a 20 46 | {.X../*|.X.. * F|
|00004da0| 69 6e 64 20 74 68 65 20 | 66 69 72 73 74 20 6e 6f |ind the |first no|
|00004db0| 6e 7a 65 72 6f 20 76 61 | 6c 75 65 20 69 6e 20 74 |nzero va|lue in t|
|00004dc0| 68 65 20 72 65 73 74 20 | 6f 66 20 74 68 65 20 63 |he rest |of the c|
|00004dd0| 6f 6c 75 6d 6e 0a 58 09 | 09 20 2a 20 64 6f 77 6e |olumn.X.|. * down|
|00004de0| 77 61 72 64 73 20 66 72 | 6f 6d 20 5b 63 75 72 2c |wards fr|om [cur,|
|00004df0| 63 75 72 5d 2e 20 20 49 | 66 20 74 68 65 72 65 20 |cur].  I|f there |
|00004e00| 69 73 20 6e 6f 20 73 75 | 63 68 20 76 61 6c 75 65 |is no su|ch value|
|00004e10| 2c 20 74 68 65 6e 0a 58 | 09 09 20 2a 20 74 68 65 |, then.X|.. * the|
|00004e20| 20 6d 61 74 72 69 78 20 | 69 73 20 6e 6f 74 20 69 | matrix |is not i|
|00004e30| 6e 76 65 72 74 69 62 6c | 65 2e 20 20 49 66 20 74 |nvertibl|e.  If t|
|00004e40| 68 65 20 66 69 72 73 74 | 20 6e 6f 6e 7a 65 72 6f |he first| nonzero|
|00004e50| 20 65 6e 74 72 79 0a 58 | 09 09 20 2a 20 69 73 20 | entry.X|.. * is |
|00004e60| 6e 6f 74 20 74 68 65 20 | 63 75 72 72 65 6e 74 20 |not the |current |
|00004e70| 72 6f 77 2c 20 74 68 65 | 6e 20 73 77 61 70 20 74 |row, the|n swap t|
|00004e80| 68 65 20 74 77 6f 20 72 | 6f 77 73 20 74 6f 20 6d |he two r|ows to m|
|00004e90| 61 6b 65 20 74 68 65 0a | 58 09 09 20 2a 20 63 75 |ake the.|X.. * cu|
|00004ea0| 72 72 65 6e 74 20 6f 6e | 65 20 6e 6f 6e 7a 65 72 |rrent on|e nonzer|
|00004eb0| 6f 2e 0a 58 09 09 20 2a | 2f 0a 58 09 09 72 6f 77 |o..X.. *|/.X..row|
|00004ec0| 20 3d 20 63 75 72 3b 0a | 58 09 09 76 61 6c 20 3d | = cur;.|X..val =|
|00004ed0| 20 26 6d 2d 3e 6d 5f 74 | 61 62 6c 65 5b 28 72 6f | &m->m_t|able[(ro|
|00004ee0| 77 20 2a 20 72 6f 77 73 | 29 20 2b 20 72 6f 77 5d |w * rows|) + row]|
|00004ef0| 3b 0a 58 09 09 77 68 69 | 6c 65 20 28 74 65 73 74 |;.X..whi|le (test|
|00004f00| 76 61 6c 75 65 28 76 61 | 6c 29 20 3d 3d 20 30 29 |value(va|l) == 0)|
|00004f10| 20 7b 0a 58 09 09 09 69 | 66 20 28 2b 2b 72 6f 77 | {.X...i|f (++row|
|00004f20| 20 3e 3d 20 72 6f 77 73 | 29 20 7b 0a 58 09 09 09 | >= rows|) {.X...|
|00004f30| 09 6d 61 74 66 72 65 65 | 28 6d 29 3b 0a 58 09 09 |.matfree|(m);.X..|
|00004f40| 09 09 6d 61 74 66 72 65 | 65 28 72 65 73 29 3b 0a |..matfre|e(res);.|
|00004f50| 58 09 09 09 09 65 72 72 | 6f 72 28 22 4d 61 74 72 |X....err|or("Matr|
|00004f60| 69 78 20 69 73 20 6e 6f | 74 20 69 6e 76 65 72 74 |ix is no|t invert|
|00004f70| 69 62 6c 65 22 29 3b 0a | 58 09 09 09 7d 0a 58 09 |ible");.|X...}.X.|
|00004f80| 09 09 76 61 6c 20 2b 3d | 20 72 6f 77 73 3b 0a 58 |..val +=| rows;.X|
|00004f90| 09 09 7d 0a 58 09 09 69 | 6e 76 65 72 74 76 61 6c |..}.X..i|nvertval|
|00004fa0| 75 65 28 76 61 6c 2c 20 | 26 6d 75 6c 76 61 6c 29 |ue(val, |&mulval)|
|00004fb0| 3b 0a 58 09 09 69 66 20 | 28 72 6f 77 20 21 3d 20 |;.X..if |(row != |
|00004fc0| 63 75 72 29 20 7b 0a 58 | 09 09 09 6d 61 74 73 77 |cur) {.X|...matsw|
|00004fd0| 61 70 72 6f 77 28 6d 2c | 20 72 6f 77 2c 20 63 75 |aprow(m,| row, cu|
|00004fe0| 72 29 3b 0a 58 09 09 09 | 6d 61 74 73 77 61 70 72 |r);.X...|matswapr|
|00004ff0| 6f 77 28 72 65 73 2c 20 | 72 6f 77 2c 20 63 75 72 |ow(res, |row, cur|
|00005000| 29 3b 0a 58 09 09 7d 0a | 58 09 09 2f 2a 0a 58 09 |);.X..}.|X../*.X.|
|00005010| 09 20 2a 20 4e 6f 77 20 | 66 6f 72 20 65 76 65 72 |. * Now |for ever|
|00005020| 79 20 6f 74 68 65 72 20 | 6e 6f 6e 7a 65 72 6f 20 |y other |nonzero |
|00005030| 65 6e 74 72 79 20 69 6e | 20 74 68 65 20 63 75 72 |entry in| the cur|
|00005040| 72 65 6e 74 20 63 6f 6c | 75 6d 6e 2c 20 73 75 62 |rent col|umn, sub|
|00005050| 74 72 61 63 74 0a 58 09 | 09 20 2a 20 74 68 65 20 |tract.X.|. * the |
|00005060| 61 70 70 72 6f 70 72 69 | 61 74 65 20 6d 75 6c 74 |appropri|ate mult|
|00005070| 69 70 6c 65 20 6f 66 20 | 74 68 65 20 63 75 72 72 |iple of |the curr|
|00005080| 65 6e 74 20 72 6f 77 20 | 74 6f 20 66 6f 72 63 65 |ent row |to force|
|00005090| 20 74 68 61 74 20 65 6e | 74 72 79 0a 58 09 09 20 | that en|try.X.. |
|000050a0| 2a 20 74 6f 20 62 65 63 | 6f 6d 65 20 7a 65 72 6f |* to bec|ome zero|
|000050b0| 2e 0a 58 09 09 20 2a 2f | 0a 58 09 09 76 61 6c 20 |..X.. */|.X..val |
|000050c0| 3d 20 26 6d 2d 3e 6d 5f | 74 61 62 6c 65 5b 63 75 |= &m->m_|table[cu|
|000050d0| 72 5d 3b 0a 58 09 09 66 | 6f 72 20 28 72 6f 77 20 |r];.X..f|or (row |
|000050e0| 3d 20 30 3b 20 72 6f 77 | 20 3c 20 72 6f 77 73 3b |= 0; row| < rows;|
|000050f0| 20 72 6f 77 2b 2b 2c 20 | 76 61 6c 20 2b 3d 20 72 | row++, |val += r|
|00005100| 6f 77 73 29 20 7b 0a 58 | 09 09 09 69 66 20 28 28 |ows) {.X|...if ((|
|00005110| 72 6f 77 20 3d 3d 20 63 | 75 72 29 20 7c 7c 20 28 |row == c|ur) || (|
|00005120| 74 65 73 74 76 61 6c 75 | 65 28 76 61 6c 29 20 3d |testvalu|e(val) =|
|00005130| 3d 20 30 29 29 0a 58 09 | 09 09 09 63 6f 6e 74 69 |= 0)).X.|...conti|
|00005140| 6e 75 65 3b 0a 58 09 09 | 09 6d 75 6c 76 61 6c 75 |nue;.X..|.mulvalu|
|00005150| 65 28 76 61 6c 2c 20 26 | 6d 75 6c 76 61 6c 2c 20 |e(val, &|mulval, |
|00005160| 26 74 6d 70 76 61 6c 29 | 3b 0a 58 09 09 09 6d 61 |&tmpval)|;.X...ma|
|00005170| 74 73 75 62 72 6f 77 28 | 6d 2c 20 72 6f 77 2c 20 |tsubrow(|m, row, |
|00005180| 63 75 72 2c 20 26 74 6d | 70 76 61 6c 29 3b 0a 58 |cur, &tm|pval);.X|
|00005190| 09 09 09 6d 61 74 73 75 | 62 72 6f 77 28 72 65 73 |...matsu|brow(res|
|000051a0| 2c 20 72 6f 77 2c 20 63 | 75 72 2c 20 26 74 6d 70 |, row, c|ur, &tmp|
|000051b0| 76 61 6c 29 3b 0a 58 09 | 09 09 66 72 65 65 76 61 |val);.X.|..freeva|
|000051c0| 6c 75 65 28 26 74 6d 70 | 76 61 6c 29 3b 0a 58 09 |lue(&tmp|val);.X.|
|000051d0| 09 7d 0a 58 09 09 66 72 | 65 65 76 61 6c 75 65 28 |.}.X..fr|eevalue(|
|000051e0| 26 6d 75 6c 76 61 6c 29 | 3b 0a 58 09 7d 0a 58 09 |&mulval)|;.X.}.X.|
|000051f0| 2f 2a 0a 58 09 20 2a 20 | 4e 6f 77 20 74 68 65 20 |/*.X. * |Now the |
|00005200| 6f 72 69 67 69 6e 61 6c | 20 6d 61 74 72 69 78 20 |original| matrix |
|00005210| 68 61 73 20 6e 6f 6e 7a | 65 72 6f 20 65 6e 74 72 |has nonz|ero entr|
|00005220| 69 65 73 20 6f 6e 6c 79 | 20 6f 6e 20 69 74 73 20 |ies only| on its |
|00005230| 6d 61 69 6e 20 64 69 61 | 67 6f 6e 61 6c 2e 0a 58 |main dia|gonal..X|
|00005240| 09 20 2a 20 53 63 61 6c | 65 20 74 68 65 20 72 6f |. * Scal|e the ro|
|00005250| 77 73 20 6f 66 20 74 68 | 65 20 72 65 73 75 6c 74 |ws of th|e result|
|00005260| 20 6d 61 74 72 69 78 20 | 62 79 20 74 68 65 20 69 | matrix |by the i|
|00005270| 6e 76 65 72 73 65 20 6f | 66 20 74 68 6f 73 65 20 |nverse o|f those |
|00005280| 65 6e 74 72 69 65 73 2e | 0a 58 09 20 2a 2f 0a 58 |entries.|.X. */.X|
|00005290| 09 76 61 6c 20 3d 20 6d | 2d 3e 6d 5f 74 61 62 6c |.val = m|->m_tabl|
|000052a0| 65 3b 0a 58 09 66 6f 72 | 20 28 72 6f 77 20 3d 20 |e;.X.for| (row = |
|000052b0| 30 3b 20 72 6f 77 20 3c | 20 72 6f 77 73 3b 20 72 |0; row <| rows; r|
|000052c0| 6f 77 2b 2b 29 20 7b 0a | 58 09 09 69 66 20 28 28 |ow++) {.|X..if ((|
|000052d0| 76 61 6c 2d 3e 76 5f 74 | 79 70 65 20 21 3d 20 56 |val->v_t|ype != V|
|000052e0| 5f 4e 55 4d 29 20 7c 7c | 20 21 71 69 73 6f 6e 65 |_NUM) ||| !qisone|
|000052f0| 28 76 61 6c 2d 3e 76 5f | 6e 75 6d 29 29 20 7b 0a |(val->v_|num)) {.|
|00005300| 58 09 09 09 69 6e 76 65 | 72 74 76 61 6c 75 65 28 |X...inve|rtvalue(|
|00005310| 76 61 6c 2c 20 26 6d 75 | 6c 76 61 6c 29 3b 0a 58 |val, &mu|lval);.X|
|00005320| 09 09 09 6d 61 74 6d 75 | 6c 72 6f 77 28 72 65 73 |...matmu|lrow(res|
|00005330| 2c 20 72 6f 77 2c 20 26 | 6d 75 6c 76 61 6c 29 3b |, row, &|mulval);|
|00005340| 0a 58 09 09 09 66 72 65 | 65 76 61 6c 75 65 28 26 |.X...fre|evalue(&|
|00005350| 6d 75 6c 76 61 6c 29 3b | 0a 58 09 09 7d 0a 58 09 |mulval);|.X..}.X.|
|00005360| 09 76 61 6c 20 2b 3d 20 | 28 72 6f 77 73 20 2b 20 |.val += |(rows + |
|00005370| 31 29 3b 0a 58 09 7d 0a | 58 09 6d 61 74 66 72 65 |1);.X.}.|X.matfre|
|00005380| 65 28 6d 29 3b 0a 58 09 | 72 65 74 75 72 6e 20 72 |e(m);.X.|return r|
|00005390| 65 73 3b 0a 58 7d 0a 58 | 0a 58 0a 58 2f 2a 0a 58 |es;.X}.X|.X.X/*.X|
|000053a0| 20 2a 20 43 61 6c 63 75 | 6c 61 74 65 20 74 68 65 | * Calcu|late the|
|000053b0| 20 64 65 74 65 72 6d 69 | 6e 61 6e 74 20 6f 66 20 | determi|nant of |
|000053c0| 61 20 73 71 75 61 72 65 | 20 6d 61 74 72 69 78 2e |a square| matrix.|
|000053d0| 0a 58 20 2a 20 54 68 69 | 73 20 69 73 20 64 6f 6e |.X * Thi|s is don|
|000053e0| 65 20 75 73 69 6e 67 20 | 72 6f 77 20 6f 70 65 72 |e using |row oper|
|000053f0| 61 74 69 6f 6e 73 20 74 | 6f 20 63 72 65 61 74 65 |ations t|o create|
|00005400| 20 61 6e 20 75 70 70 65 | 72 2d 64 69 61 67 6f 6e | an uppe|r-diagon|
|00005410| 61 6c 20 6d 61 74 72 69 | 78 2e 0a 58 20 2a 2f 0a |al matri|x..X */.|
|00005420| 58 56 41 4c 55 45 0a 58 | 6d 61 74 64 65 74 28 6d |XVALUE.X|matdet(m|
|00005430| 29 0a 58 09 4d 41 54 52 | 49 58 20 2a 6d 3b 0a 58 |).X.MATR|IX *m;.X|
|00005440| 7b 0a 58 09 6c 6f 6e 67 | 20 72 6f 77 73 3b 09 09 |{.X.long| rows;..|
|00005450| 2f 2a 20 6e 75 6d 62 65 | 72 20 6f 66 20 72 6f 77 |/* numbe|r of row|
|00005460| 73 20 2a 2f 0a 58 09 6c | 6f 6e 67 20 63 75 72 3b |s */.X.l|ong cur;|
|00005470| 09 09 2f 2a 20 63 75 72 | 72 65 6e 74 20 72 6f 77 |../* cur|rent row|
|00005480| 20 62 65 69 6e 67 20 77 | 6f 72 6b 65 64 20 6f 6e | being w|orked on|
|00005490| 20 2a 2f 0a 58 09 6c 6f | 6e 67 20 72 6f 77 3b 09 | */.X.lo|ng row;.|
|000054a0| 09 2f 2a 20 74 65 6d 70 | 20 72 6f 77 20 76 61 6c |./* temp| row val|
|000054b0| 75 65 73 20 2a 2f 0a 58 | 09 69 6e 74 20 6e 65 67 |ues */.X|.int neg|
|000054c0| 3b 09 09 2f 2a 20 77 68 | 65 74 68 65 72 20 74 6f |;../* wh|ether to|
|000054d0| 20 6e 65 67 61 74 65 20 | 64 65 74 65 72 6d 69 6e | negate |determin|
|000054e0| 61 6e 74 20 2a 2f 0a 58 | 09 56 41 4c 55 45 20 2a |ant */.X|.VALUE *|
|000054f0| 76 61 6c 3b 09 09 2f 2a | 20 63 75 72 72 65 6e 74 |val;../*| current|
|00005500| 20 76 61 6c 75 65 20 2a | 2f 0a 58 09 56 41 4c 55 | value *|/.X.VALU|
|00005510| 45 20 6d 75 6c 76 61 6c | 2c 20 74 6d 70 76 61 6c |E mulval|, tmpval|
|00005520| 3b 09 2f 2a 20 6f 74 68 | 65 72 20 76 61 6c 75 65 |;./* oth|er value|
|00005530| 73 20 2a 2f 0a 58 0a 58 | 09 69 66 20 28 6d 2d 3e |s */.X.X|.if (m->|
|00005540| 6d 5f 64 69 6d 20 21 3d | 20 32 29 0a 58 09 09 65 |m_dim !=| 2).X..e|
|00005550| 72 72 6f 72 28 22 4d 61 | 74 72 69 78 20 64 69 6d |rror("Ma|trix dim|
|00005560| 65 6e 73 69 6f 6e 20 6d | 75 73 74 20 62 65 20 74 |ension m|ust be t|
|00005570| 77 6f 20 66 6f 72 20 64 | 65 74 65 72 6d 69 6e 61 |wo for d|etermina|
|00005580| 6e 74 22 29 3b 0a 58 09 | 69 66 20 28 28 6d 2d 3e |nt");.X.|if ((m->|
|00005590| 6d 5f 6d 61 78 5b 30 5d | 20 2d 20 6d 2d 3e 6d 5f |m_max[0]| - m->m_|
|000055a0| 6d 69 6e 5b 30 5d 29 20 | 21 3d 20 28 6d 2d 3e 6d |min[0]) |!= (m->m|
|000055b0| 5f 6d 61 78 5b 31 5d 20 | 2d 20 6d 2d 3e 6d 5f 6d |_max[1] |- m->m_m|
|000055c0| 69 6e 5b 31 5d 29 29 0a | 58 09 09 65 72 72 6f 72 |in[1])).|X..error|
|000055d0| 28 22 4e 6f 6e 2d 73 71 | 75 61 72 65 20 6d 61 74 |("Non-sq|uare mat|
|000055e0| 72 69 78 20 66 6f 72 20 | 64 65 74 65 72 6d 69 6e |rix for |determin|
|000055f0| 61 6e 74 22 29 3b 0a 58 | 09 2f 2a 0a 58 09 20 2a |ant");.X|./*.X. *|
|00005600| 20 4c 6f 6f 70 20 6f 76 | 65 72 20 65 61 63 68 20 | Loop ov|er each |
|00005610| 72 6f 77 2c 20 61 6e 64 | 20 65 6c 69 6d 69 6e 61 |row, and| elimina|
|00005620| 74 65 20 61 6c 6c 20 6c | 6f 77 65 72 20 65 6e 74 |te all l|ower ent|
|00005630| 72 69 65 73 20 69 6e 20 | 74 68 65 0a 58 09 20 2a |ries in |the.X. *|
|00005640| 20 63 6f 72 72 65 73 70 | 6f 6e 64 69 6e 67 20 63 | corresp|onding c|
|00005650| 6f 6c 75 6d 6e 20 62 79 | 20 75 73 69 6e 67 20 72 |olumn by| using r|
|00005660| 6f 77 20 6f 70 65 72 61 | 74 69 6f 6e 73 2e 20 20 |ow opera|tions.  |
|00005670| 43 6f 70 79 20 74 68 65 | 20 6f 72 69 67 69 6e 61 |Copy the| origina|
|00005680| 6c 0a 58 09 20 2a 20 6d | 61 74 72 69 78 20 73 6f |l.X. * m|atrix so|
|00005690| 20 74 68 61 74 20 77 65 | 20 64 6f 6e 27 74 20 64 | that we| don't d|
|000056a0| 65 73 74 72 6f 79 20 69 | 74 2e 0a 58 09 20 2a 2f |estroy i|t..X. */|
|000056b0| 0a 58 09 6e 65 67 20 3d | 20 30 3b 0a 58 09 6d 20 |.X.neg =| 0;.X.m |
|000056c0| 3d 20 6d 61 74 63 6f 70 | 79 28 6d 29 3b 0a 58 09 |= matcop|y(m);.X.|
|000056d0| 72 6f 77 73 20 3d 20 28 | 6d 2d 3e 6d 5f 6d 61 78 |rows = (|m->m_max|
|000056e0| 5b 30 5d 20 2d 20 6d 2d | 3e 6d 5f 6d 69 6e 5b 30 |[0] - m-|>m_min[0|
|000056f0| 5d 20 2b 20 31 29 3b 0a | 58 09 66 6f 72 20 28 63 |] + 1);.|X.for (c|
|00005700| 75 72 20 3d 20 30 3b 20 | 63 75 72 20 3c 20 72 6f |ur = 0; |cur < ro|
|00005710| 77 73 3b 20 63 75 72 2b | 2b 29 20 7b 0a 58 09 09 |ws; cur+|+) {.X..|
|00005720| 2f 2a 0a 58 09 09 20 2a | 20 46 69 6e 64 20 74 68 |/*.X.. *| Find th|
|00005730| 65 20 66 69 72 73 74 20 | 6e 6f 6e 7a 65 72 6f 20 |e first |nonzero |
|00005740| 76 61 6c 75 65 20 69 6e | 20 74 68 65 20 72 65 73 |value in| the res|
|00005750| 74 20 6f 66 20 74 68 65 | 20 63 6f 6c 75 6d 6e 0a |t of the| column.|
|00005760| 58 09 09 20 2a 20 64 6f | 77 6e 77 61 72 64 73 20 |X.. * do|wnwards |
|00005770| 66 72 6f 6d 20 5b 63 75 | 72 2c 63 75 72 5d 2e 20 |from [cu|r,cur]. |
|00005780| 20 49 66 20 74 68 65 72 | 65 20 69 73 20 6e 6f 20 | If ther|e is no |
|00005790| 73 75 63 68 20 76 61 6c | 75 65 2c 20 74 68 65 6e |such val|ue, then|
|000057a0| 0a 58 09 09 20 2a 20 74 | 68 65 20 64 65 74 65 72 |.X.. * t|he deter|
|000057b0| 6d 69 6e 61 6e 74 20 69 | 73 20 7a 65 72 6f 2e 20 |minant i|s zero. |
|000057c0| 20 49 66 20 74 68 65 20 | 66 69 72 73 74 20 6e 6f | If the |first no|
|000057d0| 6e 7a 65 72 6f 20 65 6e | 74 72 79 20 69 73 20 6e |nzero en|try is n|
|000057e0| 6f 74 0a 58 09 09 20 2a | 20 74 68 65 20 63 75 72 |ot.X.. *| the cur|
|000057f0| 72 65 6e 74 20 72 6f 77 | 2c 20 74 68 65 6e 20 73 |rent row|, then s|
|00005800| 77 61 70 20 74 68 65 20 | 74 77 6f 20 72 6f 77 73 |wap the |two rows|
|00005810| 20 74 6f 20 6d 61 6b 65 | 20 74 68 65 20 63 75 72 | to make| the cur|
|00005820| 72 65 6e 74 0a 58 09 09 | 20 2a 20 6f 6e 65 20 6e |rent.X..| * one n|
|00005830| 6f 6e 7a 65 72 6f 2c 20 | 61 6e 64 20 72 65 6d 65 |onzero, |and reme|
|00005840| 6d 62 65 72 20 74 68 61 | 74 20 74 68 65 20 64 65 |mber tha|t the de|
|00005850| 74 65 72 6d 69 6e 61 6e | 74 20 63 68 61 6e 67 65 |terminan|t change|
|00005860| 73 20 73 69 67 6e 2e 0a | 58 09 09 20 2a 2f 0a 58 |s sign..|X.. */.X|
|00005870| 09 09 72 6f 77 20 3d 20 | 63 75 72 3b 0a 58 09 09 |..row = |cur;.X..|
|00005880| 76 61 6c 20 3d 20 26 6d | 2d 3e 6d 5f 74 61 62 6c |val = &m|->m_tabl|
|00005890| 65 5b 28 72 6f 77 20 2a | 20 72 6f 77 73 29 20 2b |e[(row *| rows) +|
|000058a0| 20 72 6f 77 5d 3b 0a 58 | 09 09 77 68 69 6c 65 20 | row];.X|..while |
|000058b0| 28 74 65 73 74 76 61 6c | 75 65 28 76 61 6c 29 20 |(testval|ue(val) |
|000058c0| 3d 3d 20 30 29 20 7b 0a | 58 09 09 09 69 66 20 28 |== 0) {.|X...if (|
|000058d0| 2b 2b 72 6f 77 20 3e 3d | 20 72 6f 77 73 29 20 7b |++row >=| rows) {|
|000058e0| 0a 58 09 09 09 09 6d 61 | 74 66 72 65 65 28 6d 29 |.X....ma|tfree(m)|
|000058f0| 3b 0a 58 09 09 09 09 6d | 75 6c 76 61 6c 2e 76 5f |;.X....m|ulval.v_|
|00005900| 74 79 70 65 20 3d 20 56 | 5f 4e 55 4d 3b 0a 58 09 |type = V|_NUM;.X.|
|00005910| 09 09 09 6d 75 6c 76 61 | 6c 2e 76 5f 6e 75 6d 20 |...mulva|l.v_num |
|00005920| 3d 20 71 6c 69 6e 6b 28 | 26 5f 71 7a 65 72 6f 5f |= qlink(|&_qzero_|
|00005930| 29 3b 0a 58 09 09 09 09 | 72 65 74 75 72 6e 20 6d |);.X....|return m|
|00005940| 75 6c 76 61 6c 3b 0a 58 | 09 09 09 7d 0a 58 09 09 |ulval;.X|...}.X..|
|00005950| 09 76 61 6c 20 2b 3d 20 | 72 6f 77 73 3b 0a 58 09 |.val += |rows;.X.|
|00005960| 09 7d 0a 58 09 09 69 6e | 76 65 72 74 76 61 6c 75 |.}.X..in|vertvalu|
|00005970| 65 28 76 61 6c 2c 20 26 | 6d 75 6c 76 61 6c 29 3b |e(val, &|mulval);|
|00005980| 0a 58 09 09 69 66 20 28 | 72 6f 77 20 21 3d 20 63 |.X..if (|row != c|
|00005990| 75 72 29 20 7b 0a 58 09 | 09 09 6d 61 74 73 77 61 |ur) {.X.|..matswa|
|000059a0| 70 72 6f 77 28 6d 2c 20 | 72 6f 77 2c 20 63 75 72 |prow(m, |row, cur|
|000059b0| 29 3b 0a 58 09 09 09 6e | 65 67 20 3d 20 21 6e 65 |);.X...n|eg = !ne|
|000059c0| 67 3b 0a 58 09 09 7d 0a | 58 09 09 2f 2a 0a 58 09 |g;.X..}.|X../*.X.|
|000059d0| 09 20 2a 20 4e 6f 77 20 | 66 6f 72 20 65 76 65 72 |. * Now |for ever|
|000059e0| 79 20 6f 74 68 65 72 20 | 6e 6f 6e 7a 65 72 6f 20 |y other |nonzero |
|000059f0| 65 6e 74 72 79 20 6c 6f | 77 65 72 20 64 6f 77 6e |entry lo|wer down|
|00005a00| 20 69 6e 20 74 68 65 20 | 63 75 72 72 65 6e 74 20 | in the |current |
|00005a10| 63 6f 6c 75 6d 6e 2c 0a | 58 09 09 20 2a 20 73 75 |column,.|X.. * su|
|00005a20| 62 74 72 61 63 74 20 74 | 68 65 20 61 70 70 72 6f |btract t|he appro|
|00005a30| 70 72 69 61 74 65 20 6d | 75 6c 74 69 70 6c 65 20 |priate m|ultiple |
|00005a40| 6f 66 20 74 68 65 20 63 | 75 72 72 65 6e 74 20 72 |of the c|urrent r|
|00005a50| 6f 77 20 74 6f 20 66 6f | 72 63 65 20 74 68 61 74 |ow to fo|rce that|
|00005a60| 0a 58 09 09 20 2a 20 65 | 6e 74 72 79 20 74 6f 20 |.X.. * e|ntry to |
|00005a70| 62 65 63 6f 6d 65 20 7a | 65 72 6f 2e 0a 58 09 09 |become z|ero..X..|
|00005a80| 20 2a 2f 0a 58 09 09 72 | 6f 77 20 3d 20 63 75 72 | */.X..r|ow = cur|
|00005a90| 20 2b 20 31 3b 0a 58 09 | 09 76 61 6c 20 3d 20 26 | + 1;.X.|.val = &|
|00005aa0| 6d 2d 3e 6d 5f 74 61 62 | 6c 65 5b 28 72 6f 77 20 |m->m_tab|le[(row |
|00005ab0| 2a 20 72 6f 77 73 29 20 | 2b 20 63 75 72 5d 3b 0a |* rows) |+ cur];.|
|00005ac0| 58 09 09 66 6f 72 20 28 | 3b 20 72 6f 77 20 3c 20 |X..for (|; row < |
|00005ad0| 72 6f 77 73 3b 20 72 6f | 77 2b 2b 2c 20 76 61 6c |rows; ro|w++, val|
|00005ae0| 20 2b 3d 20 72 6f 77 73 | 29 20 7b 0a 58 09 09 09 | += rows|) {.X...|
|00005af0| 69 66 20 28 74 65 73 74 | 76 61 6c 75 65 28 76 61 |if (test|value(va|
|00005b00| 6c 29 20 3d 3d 20 30 29 | 0a 58 09 09 09 09 63 6f |l) == 0)|.X....co|
|00005b10| 6e 74 69 6e 75 65 3b 0a | 58 09 09 09 6d 75 6c 76 |ntinue;.|X...mulv|
|00005b20| 61 6c 75 65 28 76 61 6c | 2c 20 26 6d 75 6c 76 61 |alue(val|, &mulva|
|00005b30| 6c 2c 20 26 74 6d 70 76 | 61 6c 29 3b 0a 58 09 09 |l, &tmpv|al);.X..|
|00005b40| 09 6d 61 74 73 75 62 72 | 6f 77 28 6d 2c 20 72 6f |.matsubr|ow(m, ro|
|00005b50| 77 2c 20 63 75 72 2c 20 | 26 74 6d 70 76 61 6c 29 |w, cur, |&tmpval)|
|00005b60| 3b 0a 58 09 09 09 66 72 | 65 65 76 61 6c 75 65 28 |;.X...fr|eevalue(|
|00005b70| 26 74 6d 70 76 61 6c 29 | 3b 0a 58 09 09 7d 0a 58 |&tmpval)|;.X..}.X|
|00005b80| 09 09 66 72 65 65 76 61 | 6c 75 65 28 26 6d 75 6c |..freeva|lue(&mul|
|00005b90| 76 61 6c 29 3b 0a 58 09 | 7d 0a 58 09 2f 2a 0a 58 |val);.X.|}.X./*.X|
|00005ba0| 09 20 2a 20 4e 6f 77 20 | 74 68 65 20 6d 61 74 72 |. * Now |the matr|
|00005bb0| 69 78 20 69 73 20 75 70 | 70 65 72 2d 64 69 61 67 |ix is up|per-diag|
|00005bc0| 6f 6e 61 6c 2c 20 61 6e | 64 20 74 68 65 20 64 65 |onal, an|d the de|
|00005bd0| 74 65 72 6d 69 6e 61 6e | 74 20 69 73 20 74 68 65 |terminan|t is the|
|00005be0| 0a 58 09 20 2a 20 70 72 | 6f 64 75 63 74 20 6f 66 |.X. * pr|oduct of|
|00005bf0| 20 74 68 65 20 6d 61 69 | 6e 20 64 69 61 67 6f 6e | the mai|n diagon|
|00005c00| 61 6c 20 65 6e 74 72 69 | 65 73 2c 20 61 6e 64 20 |al entri|es, and |
|00005c10| 69 73 20 70 6f 73 73 69 | 62 6c 79 20 6e 65 67 61 |is possi|bly nega|
|00005c20| 74 65 64 2e 0a 58 09 20 | 2a 2f 0a 58 09 76 61 6c |ted..X. |*/.X.val|
|00005c30| 20 3d 20 6d 2d 3e 6d 5f | 74 61 62 6c 65 3b 0a 58 | = m->m_|table;.X|
|00005c40| 09 6d 75 6c 76 61 6c 2e | 76 5f 74 79 70 65 20 3d |.mulval.|v_type =|
|00005c50| 20 56 5f 4e 55 4d 3b 0a | 58 09 6d 75 6c 76 61 6c | V_NUM;.|X.mulval|
|00005c60| 2e 76 5f 6e 75 6d 20 3d | 20 71 6c 69 6e 6b 28 26 |.v_num =| qlink(&|
|00005c70| 5f 71 6f 6e 65 5f 29 3b | 0a 58 09 66 6f 72 20 28 |_qone_);|.X.for (|
|00005c80| 72 6f 77 20 3d 20 30 3b | 20 72 6f 77 20 3c 20 72 |row = 0;| row < r|
|00005c90| 6f 77 73 3b 20 72 6f 77 | 2b 2b 29 20 7b 0a 58 09 |ows; row|++) {.X.|
|00005ca0| 09 6d 75 6c 76 61 6c 75 | 65 28 26 6d 75 6c 76 61 |.mulvalu|e(&mulva|
|00005cb0| 6c 2c 20 76 61 6c 2c 20 | 26 74 6d 70 76 61 6c 29 |l, val, |&tmpval)|
|00005cc0| 3b 0a 58 09 09 66 72 65 | 65 76 61 6c 75 65 28 26 |;.X..fre|evalue(&|
|00005cd0| 6d 75 6c 76 61 6c 29 3b | 0a 58 09 09 6d 75 6c 76 |mulval);|.X..mulv|
|00005ce0| 61 6c 20 3d 20 74 6d 70 | 76 61 6c 3b 0a 58 09 09 |al = tmp|val;.X..|
|00005cf0| 76 61 6c 20 2b 3d 20 28 | 72 6f 77 73 20 2b 20 31 |val += (|rows + 1|
|00005d00| 29 3b 0a 58 09 7d 0a 58 | 09 6d 61 74 66 72 65 65 |);.X.}.X|.matfree|
|00005d10| 28 6d 29 3b 0a 58 09 69 | 66 20 28 6e 65 67 29 20 |(m);.X.i|f (neg) |
|00005d20| 7b 0a 58 09 09 6e 65 67 | 76 61 6c 75 65 28 26 6d |{.X..neg|value(&m|
|00005d30| 75 6c 76 61 6c 2c 20 26 | 74 6d 70 76 61 6c 29 3b |ulval, &|tmpval);|
|00005d40| 0a 58 09 09 66 72 65 65 | 76 61 6c 75 65 28 26 6d |.X..free|value(&m|
|00005d50| 75 6c 76 61 6c 29 3b 0a | 58 09 09 72 65 74 75 72 |ulval);.|X..retur|
|00005d60| 6e 20 74 6d 70 76 61 6c | 3b 0a 58 09 7d 0a 58 09 |n tmpval|;.X.}.X.|
|00005d70| 72 65 74 75 72 6e 20 6d | 75 6c 76 61 6c 3b 0a 58 |return m|ulval;.X|
|00005d80| 7d 0a 58 0a 58 0a 58 2f | 2a 0a 58 20 2a 20 4c 6f |}.X.X.X/|*.X * Lo|
|00005d90| 63 61 6c 20 75 74 69 6c | 69 74 79 20 72 6f 75 74 |cal util|ity rout|
|00005da0| 69 6e 65 20 74 6f 20 73 | 77 61 70 20 74 77 6f 20 |ine to s|wap two |
|00005db0| 72 6f 77 73 20 6f 66 20 | 61 20 73 71 75 61 72 65 |rows of |a square|
|00005dc0| 20 6d 61 74 72 69 78 2e | 0a 58 20 2a 20 4e 6f 20 | matrix.|.X * No |
|00005dd0| 63 68 65 63 6b 73 20 61 | 72 65 20 6d 61 64 65 20 |checks a|re made |
|00005de0| 74 6f 20 76 65 72 69 66 | 79 20 74 68 65 20 6c 65 |to verif|y the le|
|00005df0| 67 61 6c 69 74 79 20 6f | 66 20 74 68 65 20 61 72 |gality o|f the ar|
|00005e00| 67 75 6d 65 6e 74 73 2e | 0a 58 20 2a 2f 0a 58 73 |guments.|.X */.Xs|
|00005e10| 74 61 74 69 63 20 76 6f | 69 64 0a 58 6d 61 74 73 |tatic vo|id.Xmats|
|00005e20| 77 61 70 72 6f 77 28 6d | 2c 20 72 31 2c 20 72 32 |waprow(m|, r1, r2|
|00005e30| 29 0a 58 09 4d 41 54 52 | 49 58 20 2a 6d 3b 0a 58 |).X.MATR|IX *m;.X|
|00005e40| 09 6c 6f 6e 67 20 72 31 | 2c 20 72 32 3b 0a 58 7b |.long r1|, r2;.X{|
|00005e50| 0a 58 09 72 65 67 69 73 | 74 65 72 20 56 41 4c 55 |.X.regis|ter VALU|
|00005e60| 45 20 2a 76 31 2c 20 2a | 76 32 3b 0a 58 09 72 65 |E *v1, *|v2;.X.re|
|00005e70| 67 69 73 74 65 72 20 6c | 6f 6e 67 20 72 6f 77 73 |gister l|ong rows|
|00005e80| 3b 0a 58 09 56 41 4c 55 | 45 20 74 6d 70 3b 0a 58 |;.X.VALU|E tmp;.X|
|00005e90| 0a 58 09 69 66 20 28 72 | 31 20 3d 3d 20 72 32 29 |.X.if (r|1 == r2)|
|00005ea0| 0a 58 09 09 72 65 74 75 | 72 6e 3b 0a 58 09 72 6f |.X..retu|rn;.X.ro|
|00005eb0| 77 73 20 3d 20 28 6d 2d | 3e 6d 5f 6d 61 78 5b 30 |ws = (m-|>m_max[0|
|00005ec0| 5d 20 2d 20 6d 2d 3e 6d | 5f 6d 69 6e 5b 30 5d 20 |] - m->m|_min[0] |
|00005ed0| 2b 20 31 29 3b 0a 58 09 | 76 31 20 3d 20 26 6d 2d |+ 1);.X.|v1 = &m-|
|00005ee0| 3e 6d 5f 74 61 62 6c 65 | 5b 72 31 20 2a 20 72 6f |>m_table|[r1 * ro|
|00005ef0| 77 73 5d 3b 0a 58 09 76 | 32 20 3d 20 26 6d 2d 3e |ws];.X.v|2 = &m->|
|00005f00| 6d 5f 74 61 62 6c 65 5b | 72 32 20 2a 20 72 6f 77 |m_table[|r2 * row|
|00005f10| 73 5d 3b 0a 58 09 77 68 | 69 6c 65 20 28 72 6f 77 |s];.X.wh|ile (row|
|00005f20| 73 2d 2d 20 3e 20 30 29 | 20 7b 0a 58 09 09 74 6d |s-- > 0)| {.X..tm|
|00005f30| 70 20 3d 20 2a 76 31 3b | 0a 58 09 09 2a 76 31 20 |p = *v1;|.X..*v1 |
|00005f40| 3d 20 2a 76 32 3b 0a 58 | 09 09 2a 76 32 20 3d 20 |= *v2;.X|..*v2 = |
|00005f50| 74 6d 70 3b 0a 58 09 09 | 76 31 2b 2b 3b 0a 58 09 |tmp;.X..|v1++;.X.|
|00005f60| 09 76 32 2b 2b 3b 0a 58 | 09 7d 0a 58 7d 0a 58 0a |.v2++;.X|.}.X}.X.|
|00005f70| 58 0a 58 2f 2a 0a 58 20 | 2a 20 4c 6f 63 61 6c 20 |X.X/*.X |* Local |
|00005f80| 75 74 69 6c 69 74 79 20 | 72 6f 75 74 69 6e 65 20 |utility |routine |
|00005f90| 74 6f 20 73 75 62 74 72 | 61 63 74 20 61 20 6d 75 |to subtr|act a mu|
|00005fa0| 6c 74 69 70 6c 65 20 6f | 66 20 6f 6e 65 20 72 6f |ltiple o|f one ro|
|00005fb0| 77 20 74 6f 20 61 6e 6f | 74 68 65 72 20 6f 6e 65 |w to ano|ther one|
|00005fc0| 2e 0a 58 20 2a 20 54 68 | 65 20 72 6f 77 20 74 6f |..X * Th|e row to|
|00005fd0| 20 62 65 20 63 68 61 6e | 67 65 64 20 69 73 20 6f | be chan|ged is o|
|00005fe0| 70 72 6f 77 2c 20 74 68 | 65 20 72 6f 77 20 74 6f |prow, th|e row to|
|00005ff0| 20 62 65 20 73 75 62 74 | 72 61 63 74 65 64 20 69 | be subt|racted i|
|00006000| 73 20 62 61 73 65 72 6f | 77 2e 0a 58 20 2a 20 4e |s basero|w..X * N|
|00006010| 6f 20 63 68 65 63 6b 73 | 20 61 72 65 20 6d 61 64 |o checks| are mad|
|00006020| 65 20 74 6f 20 76 65 72 | 69 66 79 20 74 68 65 20 |e to ver|ify the |
|00006030| 6c 65 67 61 6c 69 74 79 | 20 6f 66 20 74 68 65 20 |legality| of the |
|00006040| 61 72 67 75 6d 65 6e 74 | 73 2e 0a 58 20 2a 2f 0a |argument|s..X */.|
|00006050| 58 73 74 61 74 69 63 20 | 76 6f 69 64 0a 58 6d 61 |Xstatic |void.Xma|
|00006060| 74 73 75 62 72 6f 77 28 | 6d 2c 20 6f 70 72 6f 77 |tsubrow(|m, oprow|
|00006070| 2c 20 62 61 73 65 72 6f | 77 2c 20 6d 75 6c 76 61 |, basero|w, mulva|
|00006080| 6c 29 0a 58 09 4d 41 54 | 52 49 58 20 2a 6d 3b 0a |l).X.MAT|RIX *m;.|
|00006090| 58 09 6c 6f 6e 67 20 6f | 70 72 6f 77 2c 20 62 61 |X.long o|prow, ba|
|000060a0| 73 65 72 6f 77 3b 0a 58 | 09 56 41 4c 55 45 20 2a |serow;.X|.VALUE *|
|000060b0| 6d 75 6c 76 61 6c 3b 0a | 58 7b 0a 58 09 72 65 67 |mulval;.|X{.X.reg|
|000060c0| 69 73 74 65 72 20 56 41 | 4c 55 45 20 2a 76 6f 70 |ister VA|LUE *vop|
|000060d0| 2c 20 2a 76 62 61 73 65 | 3b 0a 58 09 72 65 67 69 |, *vbase|;.X.regi|
|000060e0| 73 74 65 72 20 6c 6f 6e | 67 20 65 6e 74 72 69 65 |ster lon|g entrie|
|000060f0| 73 3b 0a 58 09 56 41 4c | 55 45 20 74 6d 70 31 2c |s;.X.VAL|UE tmp1,|
|00006100| 20 74 6d 70 32 3b 0a 58 | 0a 58 09 65 6e 74 72 69 | tmp2;.X|.X.entri|
|00006110| 65 73 20 3d 20 28 6d 2d | 3e 6d 5f 6d 61 78 5b 30 |es = (m-|>m_max[0|
|00006120| 5d 20 2d 20 6d 2d 3e 6d | 5f 6d 69 6e 5b 30 5d 20 |] - m->m|_min[0] |
|00006130| 2b 20 31 29 3b 0a 58 09 | 76 6f 70 20 3d 20 26 6d |+ 1);.X.|vop = &m|
|00006140| 2d 3e 6d 5f 74 61 62 6c | 65 5b 6f 70 72 6f 77 20 |->m_tabl|e[oprow |
|00006150| 2a 20 65 6e 74 72 69 65 | 73 5d 3b 0a 58 09 76 62 |* entrie|s];.X.vb|
|00006160| 61 73 65 20 3d 20 26 6d | 2d 3e 6d 5f 74 61 62 6c |ase = &m|->m_tabl|
|00006170| 65 5b 62 61 73 65 72 6f | 77 20 2a 20 65 6e 74 72 |e[basero|w * entr|
|00006180| 69 65 73 5d 3b 0a 58 09 | 77 68 69 6c 65 20 28 65 |ies];.X.|while (e|
|00006190| 6e 74 72 69 65 73 2d 2d | 20 3e 20 30 29 20 7b 0a |ntries--| > 0) {.|
|000061a0| 58 09 09 6d 75 6c 76 61 | 6c 75 65 28 76 62 61 73 |X..mulva|lue(vbas|
|000061b0| 65 2c 20 6d 75 6c 76 61 | 6c 2c 20 26 74 6d 70 31 |e, mulva|l, &tmp1|
|000061c0| 29 3b 0a 58 09 09 73 75 | 62 76 61 6c 75 65 28 76 |);.X..su|bvalue(v|
|000061d0| 6f 70 2c 20 26 74 6d 70 | 31 2c 20 26 74 6d 70 32 |op, &tmp|1, &tmp2|
|000061e0| 29 3b 0a 58 09 09 66 72 | 65 65 76 61 6c 75 65 28 |);.X..fr|eevalue(|
|000061f0| 26 74 6d 70 31 29 3b 0a | 58 09 09 66 72 65 65 76 |&tmp1);.|X..freev|
|00006200| 61 6c 75 65 28 76 6f 70 | 29 3b 0a 58 09 09 2a 76 |alue(vop|);.X..*v|
|00006210| 6f 70 20 3d 20 74 6d 70 | 32 3b 0a 58 09 09 76 6f |op = tmp|2;.X..vo|
|00006220| 70 2b 2b 3b 0a 58 09 09 | 76 62 61 73 65 2b 2b 3b |p++;.X..|vbase++;|
|00006230| 0a 58 09 7d 0a 58 7d 0a | 58 0a 58 0a 58 23 69 66 |.X.}.X}.|X.X.X#if|
|00006240| 20 30 0a 58 2f 2a 0a 58 | 20 2a 20 4c 6f 63 61 6c | 0.X/*.X| * Local|
|00006250| 20 75 74 69 6c 69 74 79 | 20 72 6f 75 74 69 6e 65 | utility| routine|
|00006260| 20 74 6f 20 61 64 64 20 | 6f 6e 65 20 72 6f 77 20 | to add |one row |
|00006270| 74 6f 20 61 6e 6f 74 68 | 65 72 20 6f 6e 65 2e 0a |to anoth|er one..|
|00006280| 58 20 2a 20 4e 6f 20 63 | 68 65 63 6b 73 20 61 72 |X * No c|hecks ar|
|00006290| 65 20 6d 61 64 65 20 74 | 6f 20 76 65 72 69 66 79 |e made t|o verify|
|000062a0| 20 74 68 65 20 6c 65 67 | 61 6c 69 74 79 20 6f 66 | the leg|ality of|
|000062b0| 20 74 68 65 20 61 72 67 | 75 6d 65 6e 74 73 2e 0a | the arg|uments..|
|000062c0| 58 20 2a 2f 0a 58 73 74 | 61 74 69 63 20 76 6f 69 |X */.Xst|atic voi|
|000062d0| 64 0a 58 6d 61 74 61 64 | 64 72 6f 77 28 6d 2c 20 |d.Xmatad|drow(m, |
|000062e0| 72 31 2c 20 72 32 29 0a | 58 09 4d 41 54 52 49 58 |r1, r2).|X.MATRIX|
|000062f0| 20 2a 6d 3b 0a 58 09 6c | 6f 6e 67 20 72 31 3b 09 | *m;.X.l|ong r1;.|
|00006300| 2f 2a 20 72 6f 77 20 74 | 6f 20 62 65 20 61 64 64 |/* row t|o be add|
|00006310| 65 64 20 69 6e 74 6f 20 | 2a 2f 0a 58 09 6c 6f 6e |ed into |*/.X.lon|
|00006320| 67 20 72 32 3b 09 2f 2a | 20 72 6f 77 20 74 6f 20 |g r2;./*| row to |
|00006330| 61 64 64 20 2a 2f 0a 58 | 7b 0a 58 09 72 65 67 69 |add */.X|{.X.regi|
|00006340| 73 74 65 72 20 56 41 4c | 55 45 20 2a 76 31 2c 20 |ster VAL|UE *v1, |
|00006350| 2a 76 32 3b 0a 58 09 72 | 65 67 69 73 74 65 72 20 |*v2;.X.r|egister |
|00006360| 6c 6f 6e 67 20 72 6f 77 | 73 3b 0a 58 09 56 41 4c |long row|s;.X.VAL|
|00006370| 55 45 20 74 6d 70 3b 0a | 58 0a 58 09 72 6f 77 73 |UE tmp;.|X.X.rows|
|00006380| 20 3d 20 28 6d 2d 3e 6d | 5f 6d 61 78 5b 30 5d 20 | = (m->m|_max[0] |
|00006390| 2d 20 6d 2d 3e 6d 5f 6d | 69 6e 5b 30 5d 20 2b 20 |- m->m_m|in[0] + |
|000063a0| 31 29 3b 0a 58 09 76 31 | 20 3d 20 26 6d 2d 3e 6d |1);.X.v1| = &m->m|
|000063b0| 5f 74 61 62 6c 65 5b 72 | 31 20 2a 20 72 6f 77 73 |_table[r|1 * rows|
|000063c0| 5d 3b 0a 58 09 76 32 20 | 3d 20 26 6d 2d 3e 6d 5f |];.X.v2 |= &m->m_|
|000063d0| 74 61 62 6c 65 5b 72 32 | 20 2a 20 72 6f 77 73 5d |table[r2| * rows]|
|000063e0| 3b 0a 58 09 77 68 69 6c | 65 20 28 72 6f 77 73 2d |;.X.whil|e (rows-|
|000063f0| 2d 20 3e 20 30 29 20 7b | 0a 58 09 09 61 64 64 76 |- > 0) {|.X..addv|
+--------+-------------------------+-------------------------+--------+--------+