/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / misc / volume26 / newmat04 / part04

< prev

next >

view JSON data

view as text

open on a Mac

open on a PC

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 6d 69 73 63 0a 46 72 |.sources|.misc.Fr|
|00000020| 6f 6d 3a 20 72 6f 62 65 | 72 74 40 61 6d 2e 64 73 |om: robe|rt@am.ds|
|00000030| 69 72 2e 67 6f 76 74 2e | 6e 7a 20 28 52 6f 62 65 |ir.govt.|nz (Robe|
|00000040| 72 74 20 44 61 76 69 65 | 73 29 0a 53 75 62 6a 65 |rt Davie|s).Subje|
|00000050| 63 74 3a 20 20 76 32 36 | 69 30 39 30 3a 20 20 6e |ct:  v26|i090:  n|
|00000060| 65 77 6d 61 74 30 34 20 | 2d 20 41 20 6d 61 74 72 |ewmat04 |- A matr|
|00000070| 69 78 20 70 61 63 6b 61 | 67 65 20 69 6e 20 43 2b |ix packa|ge in C+|
|00000080| 2b 2c 20 50 61 72 74 30 | 34 2f 30 35 0a 4d 65 73 |+, Part0|4/05.Mes|
|00000090| 73 61 67 65 2d 49 44 3a | 20 3c 31 39 39 31 4e 6f |sage-ID:| <1991No|
|000000a0| 76 33 30 2e 32 30 35 33 | 35 33 2e 34 30 32 31 40 |v30.2053|53.4021@|
|000000b0| 73 70 61 72 6b 79 2e 69 | 6d 64 2e 73 74 65 72 6c |sparky.i|md.sterl|
|000000c0| 69 6e 67 2e 63 6f 6d 3e | 0a 58 2d 4d 64 34 2d 53 |ing.com>|.X-Md4-S|
|000000d0| 69 67 6e 61 74 75 72 65 | 3a 20 39 38 30 35 33 32 |ignature|: 980532|
|000000e0| 36 62 30 38 31 39 32 35 | 31 37 33 63 31 31 64 32 |6b081925|173c11d2|
|000000f0| 32 39 35 36 35 65 64 61 | 62 30 0a 44 61 74 65 3a |29565eda|b0.Date:|
|00000100| 20 53 61 74 2c 20 33 30 | 20 4e 6f 76 20 31 39 39 | Sat, 30| Nov 199|
|00000110| 31 20 32 30 3a 35 33 3a | 35 33 20 47 4d 54 0a 41 |1 20:53:|53 GMT.A|
|00000120| 70 70 72 6f 76 65 64 3a | 20 6b 65 6e 74 40 73 70 |pproved:| kent@sp|
|00000130| 61 72 6b 79 2e 69 6d 64 | 2e 73 74 65 72 6c 69 6e |arky.imd|.sterlin|
|00000140| 67 2e 63 6f 6d 0a 0a 53 | 75 62 6d 69 74 74 65 64 |g.com..S|ubmitted|
|00000150| 2d 62 79 3a 20 72 6f 62 | 65 72 74 40 61 6d 2e 64 |-by: rob|ert@am.d|
|00000160| 73 69 72 2e 67 6f 76 74 | 2e 6e 7a 20 28 52 6f 62 |sir.govt|.nz (Rob|
|00000170| 65 72 74 20 44 61 76 69 | 65 73 29 0a 50 6f 73 74 |ert Davi|es).Post|
|00000180| 69 6e 67 2d 6e 75 6d 62 | 65 72 3a 20 56 6f 6c 75 |ing-numb|er: Volu|
|00000190| 6d 65 20 32 36 2c 20 49 | 73 73 75 65 20 39 30 0a |me 26, I|ssue 90.|
|000001a0| 41 72 63 68 69 76 65 2d | 6e 61 6d 65 3a 20 6e 65 |Archive-|name: ne|
|000001b0| 77 6d 61 74 30 34 2f 70 | 61 72 74 30 34 0a 45 6e |wmat04/p|art04.En|
|000001c0| 76 69 72 6f 6e 6d 65 6e | 74 3a 20 43 2b 2b 0a 53 |vironmen|t: C++.S|
|000001d0| 75 70 65 72 73 65 64 65 | 73 3a 20 6e 65 77 6d 61 |upersede|s: newma|
|000001e0| 74 30 32 3a 20 56 6f 6c | 75 6d 65 20 32 31 2c 20 |t02: Vol|ume 21, |
|000001f0| 49 73 73 75 65 20 34 39 | 2d 35 32 0a 0a 23 21 20 |Issue 49|-52..#! |
|00000200| 2f 62 69 6e 2f 73 68 0a | 23 20 54 68 69 73 20 69 |/bin/sh.|# This i|
|00000210| 73 20 61 20 73 68 65 6c | 6c 20 61 72 63 68 69 76 |s a shel|l archiv|
|00000220| 65 2e 20 20 52 65 6d 6f | 76 65 20 61 6e 79 74 68 |e.  Remo|ve anyth|
|00000230| 69 6e 67 20 62 65 66 6f | 72 65 20 74 68 69 73 20 |ing befo|re this |
|00000240| 6c 69 6e 65 2c 20 74 68 | 65 6e 20 75 6e 70 61 63 |line, th|en unpac|
|00000250| 6b 0a 23 20 69 74 20 62 | 79 20 73 61 76 69 6e 67 |k.# it b|y saving|
|00000260| 20 69 74 20 69 6e 74 6f | 20 61 20 66 69 6c 65 20 | it into| a file |
|00000270| 61 6e 64 20 74 79 70 69 | 6e 67 20 22 73 68 20 66 |and typi|ng "sh f|
|00000280| 69 6c 65 22 2e 20 20 54 | 6f 20 6f 76 65 72 77 72 |ile".  T|o overwr|
|00000290| 69 74 65 20 65 78 69 73 | 74 69 6e 67 0a 23 20 66 |ite exis|ting.# f|
|000002a0| 69 6c 65 73 2c 20 74 79 | 70 65 20 22 73 68 20 66 |iles, ty|pe "sh f|
|000002b0| 69 6c 65 20 2d 63 22 2e | 20 20 59 6f 75 20 63 61 |ile -c".|  You ca|
|000002c0| 6e 20 61 6c 73 6f 20 66 | 65 65 64 20 74 68 69 73 |n also f|eed this|
|000002d0| 20 61 73 20 73 74 61 6e | 64 61 72 64 20 69 6e 70 | as stan|dard inp|
|000002e0| 75 74 20 76 69 61 0a 23 | 20 75 6e 73 68 61 72 2c |ut via.#| unshar,|
|000002f0| 20 6f 72 20 62 79 20 74 | 79 70 69 6e 67 20 22 73 | or by t|yping "s|
|00000300| 68 20 3c 66 69 6c 65 22 | 2c 20 65 2e 67 2e 2e 20 |h <file"|, e.g.. |
|00000310| 20 49 66 20 74 68 69 73 | 20 61 72 63 68 69 76 65 | If this| archive|
|00000320| 20 69 73 20 63 6f 6d 70 | 6c 65 74 65 2c 20 79 6f | is comp|lete, yo|
|00000330| 75 0a 23 20 77 69 6c 6c | 20 73 65 65 20 74 68 65 |u.# will| see the|
|00000340| 20 66 6f 6c 6c 6f 77 69 | 6e 67 20 6d 65 73 73 61 | followi|ng messa|
|00000350| 67 65 20 61 74 20 74 68 | 65 20 65 6e 64 3a 0a 23 |ge at th|e end:.#|
|00000360| 09 09 22 45 6e 64 20 6f | 66 20 73 68 65 6c 6c 20 |.."End o|f shell |
|00000370| 61 72 63 68 69 76 65 2e | 22 0a 23 20 43 6f 6e 74 |archive.|".# Cont|
|00000380| 65 6e 74 73 3a 20 20 6e | 65 77 6d 61 74 35 2e 63 |ents:  n|ewmat5.c|
|00000390| 78 78 20 6e 65 77 6d 61 | 74 36 2e 63 78 78 20 6e |xx newma|t6.cxx n|
|000003a0| 65 77 6d 61 74 37 2e 63 | 78 78 20 6e 65 77 6d 61 |ewmat7.c|xx newma|
|000003b0| 74 38 2e 63 78 78 0a 23 | 20 20 20 6e 65 77 6d 61 |t8.cxx.#|   newma|
|000003c0| 74 39 2e 63 78 78 20 6e | 65 77 6d 61 74 72 6d 2e |t9.cxx n|ewmatrm.|
|000003d0| 63 78 78 0a 23 20 57 72 | 61 70 70 65 64 20 62 79 |cxx.# Wr|apped by|
|000003e0| 20 72 6f 62 65 72 74 40 | 6b 65 61 20 6f 6e 20 53 | robert@|kea on S|
|000003f0| 61 74 20 4e 6f 76 20 33 | 30 20 31 38 3a 32 39 3a |at Nov 3|0 18:29:|
|00000400| 34 33 20 31 39 39 31 0a | 50 41 54 48 3d 2f 62 69 |43 1991.|PATH=/bi|
|00000410| 6e 3a 2f 75 73 72 2f 62 | 69 6e 3a 2f 75 73 72 2f |n:/usr/b|in:/usr/|
|00000420| 75 63 62 20 3b 20 65 78 | 70 6f 72 74 20 50 41 54 |ucb ; ex|port PAT|
|00000430| 48 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 6e 65 |H.if tes|t -f 'ne|
|00000440| 77 6d 61 74 35 2e 63 78 | 78 27 20 2d 61 20 22 24 |wmat5.cx|x' -a "$|
|00000450| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|00000460| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen .  e|cho shar|
|00000470| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 63 6c 6f 62 62 |: Will n|ot clobb|
|00000480| 65 72 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |er exist|ing file|
|00000490| 20 5c 22 27 6e 65 77 6d | 61 74 35 2e 63 78 78 27 | \"'newm|at5.cxx'|
|000004a0| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|000004b0| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|000004c0| 27 6e 65 77 6d 61 74 35 | 2e 63 78 78 27 5c 22 20 |'newmat5|.cxx'\" |
|000004d0| 5c 28 37 37 37 35 20 63 | 68 61 72 61 63 74 65 72 |\(7775 c|haracter|
|000004e0| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|000004f0| 20 3e 27 6e 65 77 6d 61 | 74 35 2e 63 78 78 27 20 | >'newma|t5.cxx' |
|00000500| 3c 3c 27 45 4e 44 5f 4f | 46 5f 46 49 4c 45 27 0a |<<'END_O|F_FILE'.|
|00000510| 58 2f 2f 24 24 20 6e 65 | 77 6d 61 74 35 2e 63 78 |X//$$ ne|wmat5.cx|
|00000520| 78 20 20 20 20 20 20 20 | 20 20 54 72 61 6e 73 70 |x       |  Transp|
|00000530| 6f 73 65 2c 20 65 76 61 | 6c 75 61 74 65 20 65 74 |ose, eva|luate et|
|00000540| 63 0a 58 0a 58 2f 2f 20 | 43 6f 70 79 72 69 67 68 |c.X.X// |Copyrigh|
|00000550| 74 20 28 43 29 20 31 39 | 39 31 3a 20 52 20 42 20 |t (C) 19|91: R B |
|00000560| 44 61 76 69 65 73 20 61 | 6e 64 20 44 53 49 52 0a |Davies a|nd DSIR.|
|00000570| 58 0a 58 23 69 6e 63 6c | 75 64 65 20 22 69 6e 63 |X.X#incl|ude "inc|
|00000580| 6c 75 64 65 2e 68 78 78 | 22 0a 58 0a 58 23 69 6e |lude.hxx|".X.X#in|
|00000590| 63 6c 75 64 65 20 22 6e | 65 77 6d 61 74 2e 68 78 |clude "n|ewmat.hx|
|000005a0| 78 22 0a 58 23 69 6e 63 | 6c 75 64 65 20 22 6e 65 |x".X#inc|lude "ne|
|000005b0| 77 6d 61 74 72 63 2e 68 | 78 78 22 0a 58 0a 58 2f |wmatrc.h|xx".X.X/|
|000005c0| 2f 23 64 65 66 69 6e 65 | 20 52 45 50 4f 52 54 20 |/#define| REPORT |
|000005d0| 7b 20 73 74 61 74 69 63 | 20 45 78 65 43 6f 75 6e |{ static| ExeCoun|
|000005e0| 74 65 72 20 45 78 65 43 | 6f 75 6e 74 28 5f 5f 4c |ter ExeC|ount(__L|
|000005f0| 49 4e 45 5f 5f 2c 35 29 | 3b 20 45 78 65 43 6f 75 |INE__,5)|; ExeCou|
|00000600| 6e 74 2b 2b 3b 20 7d 0a | 58 0a 58 23 64 65 66 69 |nt++; }.|X.X#defi|
|00000610| 6e 65 20 52 45 50 4f 52 | 54 20 7b 7d 0a 58 0a 58 |ne REPOR|T {}.X.X|
|00000620| 0a 58 2f 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X/*****|********|
|00000630| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 20 63 61 72 72 |********|*** carr|
|00000640| 79 20 6f 75 74 20 6f 70 | 65 72 61 74 69 6f 6e 73 |y out op|erations|
|00000650| 20 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a | *******|********|
|00000660| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2f |********|*******/|
|00000670| 0a 58 0a 58 0a 58 47 65 | 6e 65 72 61 6c 4d 61 74 |.X.X.XGe|neralMat|
|00000680| 72 69 78 2a 20 47 65 6e | 65 72 61 6c 4d 61 74 72 |rix* Gen|eralMatr|
|00000690| 69 78 3a 3a 54 72 61 6e | 73 70 6f 73 65 28 4d 61 |ix::Tran|spose(Ma|
|000006a0| 74 72 69 78 54 79 70 65 | 20 6d 74 29 0a 58 7b 0a |trixType| mt).X{.|
|000006b0| 58 20 20 20 69 66 20 28 | 28 69 6e 74 29 6d 74 3d |X   if (|(int)mt=|
|000006c0| 3d 30 29 20 6d 74 20 3d | 20 54 79 70 65 28 29 2e |=0) mt =| Type().|
|000006d0| 74 28 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 2f |t();    |       /|
|000006e0| 2f 20 74 79 70 65 20 6f | 66 20 74 72 61 6e 73 70 |/ type o|f transp|
|000006f0| 6f 73 65 64 20 6d 61 74 | 72 69 78 0a 58 20 20 20 |osed mat|rix.X   |
|00000700| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00000710| 6d 31 20 3d 20 6d 74 2e | 4e 65 77 28 6e 63 6f 6c |m1 = mt.|New(ncol|
|00000720| 73 2c 6e 72 6f 77 73 29 | 3b 20 67 6d 31 2d 3e 52 |s,nrows)|; gm1->R|
|00000730| 65 6c 65 61 73 65 41 6e | 64 44 65 6c 65 74 65 28 |eleaseAn|dDelete(|
|00000740| 29 3b 0a 58 0a 58 20 20 | 20 69 66 20 28 6d 74 20 |);.X.X  | if (mt |
|00000750| 3d 3d 20 54 79 70 65 28 | 29 2e 74 28 29 29 0a 58 |== Type(|).t()).X|
|00000760| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 52 45 50 4f |   {.X  |    REPO|
|00000770| 52 54 0a 58 20 20 20 20 | 20 20 66 6f 72 20 28 69 |RT.X    |  for (i|
|00000780| 6e 74 20 69 3d 30 3b 20 | 69 3c 6e 63 6f 6c 73 3b |nt i=0; |i<ncols;|
|00000790| 20 69 2b 2b 29 0a 58 20 | 20 20 20 20 20 7b 0a 58 | i++).X |     {.X|
|000007a0| 09 20 4d 61 74 72 69 78 | 52 6f 77 20 6d 72 28 67 |. Matrix|Row mr(g|
|000007b0| 6d 31 2c 20 53 74 6f 72 | 65 4f 6e 45 78 69 74 2b |m1, Stor|eOnExit+|
|000007c0| 44 69 72 65 63 74 50 61 | 72 74 2c 20 69 29 3b 0a |DirectPa|rt, i);.|
|000007d0| 58 20 20 20 20 20 20 20 | 20 20 4d 61 74 72 69 78 |X       |  Matrix|
|000007e0| 43 6f 6c 20 6d 63 28 74 | 68 69 73 2c 20 6d 72 2e |Col mc(t|his, mr.|
|000007f0| 53 74 6f 72 65 28 29 2c | 20 4c 6f 61 64 4f 6e 45 |Store(),| LoadOnE|
|00000800| 6e 74 72 79 2c 20 69 29 | 3b 0a 58 20 20 20 20 20 |ntry, i)|;.X     |
|00000810| 20 7d 0a 58 20 20 20 7d | 0a 58 20 20 20 65 6c 73 | }.X   }|.X   els|
|00000820| 65 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 52 |e.X   {.|X      R|
|00000830| 45 50 4f 52 54 0a 58 20 | 20 20 20 20 20 4d 61 74 |EPORT.X |     Mat|
|00000840| 72 69 78 52 6f 77 20 6d | 72 28 74 68 69 73 2c 20 |rixRow m|r(this, |
|00000850| 4c 6f 61 64 4f 6e 45 6e | 74 72 79 29 3b 0a 58 20 |LoadOnEn|try);.X |
|00000860| 20 20 20 20 20 4d 61 74 | 72 69 78 43 6f 6c 20 6d |     Mat|rixCol m|
|00000870| 63 28 67 6d 31 2c 20 53 | 74 6f 72 65 4f 6e 45 78 |c(gm1, S|toreOnEx|
|00000880| 69 74 2b 44 69 72 65 63 | 74 50 61 72 74 29 3b 0a |it+Direc|tPart);.|
|00000890| 58 20 20 20 20 20 20 69 | 6e 74 20 69 20 3d 20 6e |X      i|nt i = n|
|000008a0| 72 6f 77 73 3b 0a 58 20 | 20 20 20 20 20 77 68 69 |rows;.X |     whi|
|000008b0| 6c 65 20 28 69 2d 2d 29 | 20 7b 20 6d 63 2e 43 6f |le (i--)| { mc.Co|
|000008c0| 70 79 28 6d 72 29 3b 20 | 6d 72 2e 4e 65 78 74 28 |py(mr); |mr.Next(|
|000008d0| 29 3b 20 6d 63 2e 4e 65 | 78 74 28 29 3b 20 7d 0a |); mc.Ne|xt(); }.|
|000008e0| 58 20 20 20 7d 0a 58 20 | 20 20 74 44 65 6c 65 74 |X   }.X |  tDelet|
|000008f0| 65 28 29 3b 20 72 65 74 | 75 72 6e 20 67 6d 31 3b |e(); ret|urn gm1;|
|00000900| 0a 58 7d 0a 58 0a 58 47 | 65 6e 65 72 61 6c 4d 61 |.X}.X.XG|eneralMa|
|00000910| 74 72 69 78 2a 20 53 79 | 6d 6d 65 74 72 69 63 4d |trix* Sy|mmetricM|
|00000920| 61 74 72 69 78 3a 3a 54 | 72 61 6e 73 70 6f 73 65 |atrix::T|ranspose|
|00000930| 28 4d 61 74 72 69 78 54 | 79 70 65 20 6d 74 29 0a |(MatrixT|ype mt).|
|00000940| 58 7b 20 52 45 50 4f 52 | 54 20 20 72 65 74 75 72 |X{ REPOR|T  retur|
|00000950| 6e 20 45 76 61 6c 75 61 | 74 65 28 6d 74 29 3b 20 |n Evalua|te(mt); |
|00000960| 7d 0a 58 0a 58 0a 58 47 | 65 6e 65 72 61 6c 4d 61 |}.X.X.XG|eneralMa|
|00000970| 74 72 69 78 2a 20 44 69 | 61 67 6f 6e 61 6c 4d 61 |trix* Di|agonalMa|
|00000980| 74 72 69 78 3a 3a 54 72 | 61 6e 73 70 6f 73 65 28 |trix::Tr|anspose(|
|00000990| 4d 61 74 72 69 78 54 79 | 70 65 20 6d 74 29 0a 58 |MatrixTy|pe mt).X|
|000009a0| 7b 20 52 45 50 4f 52 54 | 20 72 65 74 75 72 6e 20 |{ REPORT| return |
|000009b0| 45 76 61 6c 75 61 74 65 | 28 6d 74 29 3b 20 7d 0a |Evaluate|(mt); }.|
|000009c0| 58 0a 58 42 4f 4f 4c 20 | 47 65 6e 65 72 61 6c 4d |X.XBOOL |GeneralM|
|000009d0| 61 74 72 69 78 3a 3a 49 | 73 5a 65 72 6f 28 29 20 |atrix::I|sZero() |
|000009e0| 63 6f 6e 73 74 0a 58 7b | 0a 58 20 20 20 52 45 50 |const.X{|.X   REP|
|000009f0| 4f 52 54 0a 58 20 20 20 | 72 65 61 6c 2a 20 73 3d |ORT.X   |real* s=|
|00000a00| 73 74 6f 72 65 3b 20 69 | 6e 74 20 69 3d 73 74 6f |store; i|nt i=sto|
|00000a10| 72 61 67 65 3b 0a 58 20 | 20 20 77 68 69 6c 65 20 |rage;.X |  while |
|00000a20| 28 69 2d 2d 29 20 7b 20 | 69 66 20 28 2a 73 2b 2b |(i--) { |if (*s++|
|00000a30| 29 20 72 65 74 75 72 6e | 20 46 41 4c 53 45 3b 20 |) return| FALSE; |
|00000a40| 7d 0a 58 20 20 20 72 65 | 74 75 72 6e 20 54 52 55 |}.X   re|turn TRU|
|00000a50| 45 3b 0a 58 7d 0a 58 0a | 58 47 65 6e 65 72 61 6c |E;.X}.X.|XGeneral|
|00000a60| 4d 61 74 72 69 78 2a 20 | 43 6f 6c 75 6d 6e 56 65 |Matrix* |ColumnVe|
|00000a70| 63 74 6f 72 3a 3a 54 72 | 61 6e 73 70 6f 73 65 28 |ctor::Tr|anspose(|
|00000a80| 4d 61 74 72 69 78 54 79 | 70 65 20 6d 74 29 0a 58 |MatrixTy|pe mt).X|
|00000a90| 7b 0a 58 20 20 20 52 45 | 50 4f 52 54 0a 58 20 20 |{.X   RE|PORT.X  |
|00000aa0| 20 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 | General|Matrix* |
|00000ab0| 67 6d 78 20 3d 20 6e 65 | 77 20 52 6f 77 56 65 63 |gmx = ne|w RowVec|
|00000ac0| 74 6f 72 3b 20 4d 61 74 | 72 69 78 45 72 72 6f 72 |tor; Mat|rixError|
|00000ad0| 4e 6f 53 70 61 63 65 28 | 67 6d 78 29 3b 0a 58 20 |NoSpace(|gmx);.X |
|00000ae0| 20 20 67 6d 78 2d 3e 6e | 72 6f 77 73 20 3d 20 31 |  gmx->n|rows = 1|
|00000af0| 3b 20 67 6d 78 2d 3e 6e | 63 6f 6c 73 20 3d 20 67 |; gmx->n|cols = g|
|00000b00| 6d 78 2d 3e 73 74 6f 72 | 61 67 65 20 3d 20 73 74 |mx->stor|age = st|
|00000b10| 6f 72 61 67 65 3b 0a 58 | 20 20 20 72 65 74 75 72 |orage;.X|   retur|
|00000b20| 6e 20 42 6f 72 72 6f 77 | 53 74 6f 72 65 28 67 6d |n Borrow|Store(gm|
|00000b30| 78 2c 6d 74 29 3b 0a 58 | 7d 0a 58 0a 58 47 65 6e |x,mt);.X|}.X.XGen|
|00000b40| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 52 6f 77 56 |eralMatr|ix* RowV|
|00000b50| 65 63 74 6f 72 3a 3a 54 | 72 61 6e 73 70 6f 73 65 |ector::T|ranspose|
|00000b60| 28 4d 61 74 72 69 78 54 | 79 70 65 20 6d 74 29 0a |(MatrixT|ype mt).|
|00000b70| 58 7b 0a 58 20 20 20 52 | 45 50 4f 52 54 0a 58 20 |X{.X   R|EPORT.X |
|00000b80| 20 20 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |  Genera|lMatrix*|
|00000b90| 20 67 6d 78 20 3d 20 6e | 65 77 20 43 6f 6c 75 6d | gmx = n|ew Colum|
|00000ba0| 6e 56 65 63 74 6f 72 3b | 20 4d 61 74 72 69 78 45 |nVector;| MatrixE|
|00000bb0| 72 72 6f 72 4e 6f 53 70 | 61 63 65 28 67 6d 78 29 |rrorNoSp|ace(gmx)|
|00000bc0| 3b 0a 58 20 20 20 67 6d | 78 2d 3e 6e 63 6f 6c 73 |;.X   gm|x->ncols|
|00000bd0| 20 3d 20 31 3b 20 67 6d | 78 2d 3e 6e 72 6f 77 73 | = 1; gm|x->nrows|
|00000be0| 20 3d 20 67 6d 78 2d 3e | 73 74 6f 72 61 67 65 20 | = gmx->|storage |
|00000bf0| 3d 20 73 74 6f 72 61 67 | 65 3b 0a 58 20 20 20 72 |= storag|e;.X   r|
|00000c00| 65 74 75 72 6e 20 42 6f | 72 72 6f 77 53 74 6f 72 |eturn Bo|rrowStor|
|00000c10| 65 28 67 6d 78 2c 6d 74 | 29 3b 0a 58 7d 0a 58 0a |e(gmx,mt|);.X}.X.|
|00000c20| 58 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 |XGeneral|Matrix* |
|00000c30| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 3a 3a 45 |GeneralM|atrix::E|
|00000c40| 76 61 6c 75 61 74 65 28 | 4d 61 74 72 69 78 54 79 |valuate(|MatrixTy|
|00000c50| 70 65 20 6d 74 29 0a 58 | 7b 0a 58 20 20 20 69 66 |pe mt).X|{.X   if|
|00000c60| 20 28 28 69 6e 74 29 6d | 74 3d 3d 30 29 20 7b 20 | ((int)m|t==0) { |
|00000c70| 52 45 50 4f 52 54 20 72 | 65 74 75 72 6e 20 74 68 |REPORT r|eturn th|
|00000c80| 69 73 3b 20 7d 0a 58 20 | 20 20 69 66 20 28 6d 74 |is; }.X |  if (mt|
|00000c90| 3d 3d 74 68 69 73 2d 3e | 54 79 70 65 28 29 29 20 |==this->|Type()) |
|00000ca0| 7b 20 52 45 50 4f 52 54 | 20 72 65 74 75 72 6e 20 |{ REPORT| return |
|00000cb0| 74 68 69 73 3b 20 7d 0a | 58 20 20 20 52 45 50 4f |this; }.|X   REPO|
|00000cc0| 52 54 0a 58 20 20 20 47 | 65 6e 65 72 61 6c 4d 61 |RT.X   G|eneralMa|
|00000cd0| 74 72 69 78 2a 20 67 6d | 78 20 3d 20 6d 74 2e 4e |trix* gm|x = mt.N|
|00000ce0| 65 77 28 6e 72 6f 77 73 | 2c 6e 63 6f 6c 73 29 3b |ew(nrows|,ncols);|
|00000cf0| 0a 58 20 20 20 4d 61 74 | 72 69 78 52 6f 77 20 6d |.X   Mat|rixRow m|
|00000d00| 72 28 74 68 69 73 2c 20 | 4c 6f 61 64 4f 6e 45 6e |r(this, |LoadOnEn|
|00000d10| 74 72 79 29 3b 0a 58 20 | 20 20 4d 61 74 72 69 78 |try);.X |  Matrix|
|00000d20| 52 6f 77 20 6d 72 78 28 | 67 6d 78 2c 20 53 74 6f |Row mrx(|gmx, Sto|
|00000d30| 72 65 4f 6e 45 78 69 74 | 2b 44 69 72 65 63 74 50 |reOnExit|+DirectP|
|00000d40| 61 72 74 29 3b 0a 58 20 | 20 20 69 6e 74 20 69 3d |art);.X |  int i=|
|00000d50| 6e 72 6f 77 73 3b 0a 58 | 20 20 20 77 68 69 6c 65 |nrows;.X|   while|
|00000d60| 20 28 69 2d 2d 29 20 7b | 20 6d 72 78 2e 43 6f 70 | (i--) {| mrx.Cop|
|00000d70| 79 28 6d 72 29 3b 20 6d | 72 78 2e 4e 65 78 74 28 |y(mr); m|rx.Next(|
|00000d80| 29 3b 20 6d 72 2e 4e 65 | 78 74 28 29 3b 20 7d 0a |); mr.Ne|xt(); }.|
|00000d90| 58 20 20 20 74 44 65 6c | 65 74 65 28 29 3b 20 67 |X   tDel|ete(); g|
|00000da0| 6d 78 2d 3e 52 65 6c 65 | 61 73 65 41 6e 64 44 65 |mx->Rele|aseAndDe|
|00000db0| 6c 65 74 65 28 29 3b 20 | 72 65 74 75 72 6e 20 67 |lete(); |return g|
|00000dc0| 6d 78 3b 0a 58 7d 0a 58 | 0a 58 47 65 6e 65 72 61 |mx;.X}.X|.XGenera|
|00000dd0| 6c 4d 61 74 72 69 78 2a | 20 43 6f 6e 73 74 4d 61 |lMatrix*| ConstMa|
|00000de0| 74 72 69 78 3a 3a 45 76 | 61 6c 75 61 74 65 28 4d |trix::Ev|aluate(M|
|00000df0| 61 74 72 69 78 54 79 70 | 65 20 6d 74 29 0a 58 7b |atrixTyp|e mt).X{|
|00000e00| 0a 58 20 20 20 69 66 20 | 28 28 69 6e 74 29 6d 74 |.X   if |((int)mt|
|00000e10| 3d 3d 30 20 7c 7c 20 6d | 74 3d 3d 63 67 6d 2d 3e |==0 || m|t==cgm->|
|00000e20| 54 79 70 65 28 29 29 20 | 7b 20 52 45 50 4f 52 54 |Type()) |{ REPORT|
|00000e30| 20 72 65 74 75 72 6e 20 | 28 47 65 6e 65 72 61 6c | return |(General|
|00000e40| 4d 61 74 72 69 78 2a 29 | 63 67 6d 3b 20 7d 0a 58 |Matrix*)|cgm; }.X|
|00000e50| 0a 58 20 20 20 52 45 50 | 4f 52 54 0a 58 20 20 20 |.X   REP|ORT.X   |
|00000e60| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00000e70| 6d 78 20 3d 20 63 67 6d | 2d 3e 54 79 70 65 28 29 |mx = cgm|->Type()|
|00000e80| 2e 4e 65 77 28 63 67 6d | 2d 3e 4e 72 6f 77 73 28 |.New(cgm|->Nrows(|
|00000e90| 29 2c 63 67 6d 2d 3e 4e | 63 6f 6c 73 28 29 29 3b |),cgm->N|cols());|
|00000ea0| 0a 58 20 20 20 4d 61 74 | 72 69 78 52 6f 77 20 6d |.X   Mat|rixRow m|
|00000eb0| 72 28 28 47 65 6e 65 72 | 61 6c 4d 61 74 72 69 78 |r((Gener|alMatrix|
|00000ec0| 2a 29 63 67 6d 2c 20 4c | 6f 61 64 4f 6e 45 6e 74 |*)cgm, L|oadOnEnt|
|00000ed0| 72 79 29 3b 2f 2f 61 73 | 73 75 6d 65 20 77 6f 6e |ry);//as|sume won|
|00000ee0| 27 74 20 63 68 61 6e 67 | 65 20 74 68 69 73 0a 58 |'t chang|e this.X|
|00000ef0| 20 20 20 4d 61 74 72 69 | 78 52 6f 77 20 6d 72 78 |   Matri|xRow mrx|
|00000f00| 28 67 6d 78 2c 20 53 74 | 6f 72 65 4f 6e 45 78 69 |(gmx, St|oreOnExi|
|00000f10| 74 2b 44 69 72 65 63 74 | 50 61 72 74 29 3b 0a 58 |t+Direct|Part);.X|
|00000f20| 20 20 20 69 6e 74 20 69 | 3d 63 67 6d 2d 3e 4e 72 |   int i|=cgm->Nr|
|00000f30| 6f 77 73 28 29 3b 0a 58 | 20 20 20 77 68 69 6c 65 |ows();.X|   while|
|00000f40| 20 28 69 2d 2d 29 20 7b | 20 6d 72 78 2e 43 6f 70 | (i--) {| mrx.Cop|
|00000f50| 79 28 6d 72 29 3b 20 6d | 72 78 2e 4e 65 78 74 28 |y(mr); m|rx.Next(|
|00000f60| 29 3b 20 6d 72 2e 4e 65 | 78 74 28 29 3b 20 7d 0a |); mr.Ne|xt(); }.|
|00000f70| 58 20 20 20 67 6d 78 2d | 3e 52 65 6c 65 61 73 65 |X   gmx-|>Release|
|00000f80| 41 6e 64 44 65 6c 65 74 | 65 28 29 3b 20 72 65 74 |AndDelet|e(); ret|
|00000f90| 75 72 6e 20 67 6d 78 3b | 20 2f 2f 20 6e 6f 20 74 |urn gmx;| // no t|
|00000fa0| 44 65 6c 65 74 65 0a 58 | 7d 0a 58 0a 58 47 65 6e |Delete.X|}.X.XGen|
|00000fb0| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 53 68 69 66 |eralMatr|ix* Shif|
|00000fc0| 74 65 64 4d 61 74 72 69 | 78 3a 3a 45 76 61 6c 75 |tedMatri|x::Evalu|
|00000fd0| 61 74 65 28 4d 61 74 72 | 69 78 54 79 70 65 20 6d |ate(Matr|ixType m|
|00000fe0| 74 29 0a 58 7b 0a 58 20 | 20 20 47 65 6e 65 72 61 |t).X{.X |  Genera|
|00000ff0| 6c 4d 61 74 72 69 78 2a | 20 67 6d 3d 62 6d 2d 3e |lMatrix*| gm=bm->|
|00001000| 45 76 61 6c 75 61 74 65 | 28 29 3b 0a 58 20 20 20 |Evaluate|();.X   |
|00001010| 69 66 20 28 28 69 6e 74 | 29 6d 74 3d 3d 30 29 0a |if ((int|)mt==0).|
|00001020| 58 20 20 20 20 20 20 7b | 20 4d 61 74 72 69 78 54 |X      {| MatrixT|
|00001030| 79 70 65 20 6d 74 65 71 | 65 6c 28 4d 61 74 72 69 |ype mteq|el(Matri|
|00001040| 78 54 79 70 65 3a 3a 45 | 71 45 6c 29 3b 20 6d 74 |xType::E|qEl); mt|
|00001050| 20 3d 20 67 6d 2d 3e 54 | 79 70 65 28 29 2b 6d 74 | = gm->T|ype()+mt|
|00001060| 65 71 65 6c 3b 20 7d 0a | 58 20 20 20 69 6e 74 20 |eqel; }.|X   int |
|00001070| 6e 72 3d 67 6d 2d 3e 4e | 72 6f 77 73 28 29 3b 20 |nr=gm->N|rows(); |
|00001080| 69 6e 74 20 6e 63 3d 67 | 6d 2d 3e 4e 63 6f 6c 73 |int nc=g|m->Ncols|
|00001090| 28 29 3b 0a 58 20 20 20 | 69 66 20 28 21 28 6d 74 |();.X   |if (!(mt|
|000010a0| 3d 3d 67 6d 2d 3e 54 79 | 70 65 28 29 29 29 0a 58 |==gm->Ty|pe())).X|
|000010b0| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 52 45 50 4f |   {.X  |    REPO|
|000010c0| 52 54 0a 58 20 20 20 20 | 20 20 47 65 6e 65 72 61 |RT.X    |  Genera|
|000010d0| 6c 4d 61 74 72 69 78 2a | 20 67 6d 78 20 3d 20 6d |lMatrix*| gmx = m|
|000010e0| 74 2e 4e 65 77 28 6e 72 | 2c 6e 63 29 3b 0a 58 20 |t.New(nr|,nc);.X |
|000010f0| 20 20 20 20 20 4d 61 74 | 72 69 78 52 6f 77 20 6d |     Mat|rixRow m|
|00001100| 72 28 67 6d 2c 20 4c 6f | 61 64 4f 6e 45 6e 74 72 |r(gm, Lo|adOnEntr|
|00001110| 79 29 3b 20 0a 58 20 20 | 20 20 20 20 4d 61 74 72 |y); .X  |    Matr|
|00001120| 69 78 52 6f 77 20 6d 72 | 78 28 67 6d 78 2c 20 53 |ixRow mr|x(gmx, S|
|00001130| 74 6f 72 65 4f 6e 45 78 | 69 74 2b 44 69 72 65 63 |toreOnEx|it+Direc|
|00001140| 74 50 61 72 74 29 3b 0a | 58 20 20 20 20 20 20 77 |tPart);.|X      w|
|00001150| 68 69 6c 65 20 28 6e 72 | 2d 2d 29 20 7b 20 6d 72 |hile (nr|--) { mr|
|00001160| 78 2e 41 64 64 28 6d 72 | 2c 66 29 3b 20 6d 72 78 |x.Add(mr|,f); mrx|
|00001170| 2e 4e 65 78 74 28 29 3b | 20 6d 72 2e 4e 65 78 74 |.Next();| mr.Next|
|00001180| 28 29 3b 20 7d 0a 58 20 | 20 20 20 20 20 67 6d 78 |(); }.X |     gmx|
|00001190| 2d 3e 52 65 6c 65 61 73 | 65 41 6e 64 44 65 6c 65 |->Releas|eAndDele|
|000011a0| 74 65 28 29 3b 20 67 6d | 2d 3e 74 44 65 6c 65 74 |te(); gm|->tDelet|
|000011b0| 65 28 29 3b 20 72 65 74 | 75 72 6e 20 67 6d 78 3b |e(); ret|urn gmx;|
|000011c0| 0a 58 20 20 20 7d 0a 58 | 20 20 20 65 6c 73 65 20 |.X   }.X|   else |
|000011d0| 69 66 20 28 67 6d 2d 3e | 72 65 75 73 65 28 29 29 |if (gm->|reuse())|
|000011e0| 20 7b 20 52 45 50 4f 52 | 54 20 67 6d 2d 3e 41 64 | { REPOR|T gm->Ad|
|000011f0| 64 28 66 29 3b 20 72 65 | 74 75 72 6e 20 67 6d 3b |d(f); re|turn gm;|
|00001200| 20 7d 0a 58 20 20 20 65 | 6c 73 65 0a 58 20 20 20 | }.X   e|lse.X   |
|00001210| 7b 0a 58 20 20 20 20 20 | 20 52 45 50 4f 52 54 20 |{.X     | REPORT |
|00001220| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00001230| 6d 79 20 3d 20 67 6d 2d | 3e 54 79 70 65 28 29 2e |my = gm-|>Type().|
|00001240| 4e 65 77 28 6e 72 2c 6e | 63 29 3b 0a 58 20 20 20 |New(nr,n|c);.X   |
|00001250| 20 20 20 67 6d 79 2d 3e | 52 65 6c 65 61 73 65 41 |   gmy->|ReleaseA|
|00001260| 6e 64 44 65 6c 65 74 65 | 28 29 3b 20 67 6d 79 2d |ndDelete|(); gmy-|
|00001270| 3e 41 64 64 28 67 6d 2c | 66 29 3b 20 72 65 74 75 |>Add(gm,|f); retu|
|00001280| 72 6e 20 67 6d 79 3b 0a | 58 20 20 20 7d 0a 58 7d |rn gmy;.|X   }.X}|
|00001290| 0a 58 0a 58 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |.X.XGene|ralMatri|
|000012a0| 78 2a 20 53 63 61 6c 65 | 64 4d 61 74 72 69 78 3a |x* Scale|dMatrix:|
|000012b0| 3a 45 76 61 6c 75 61 74 | 65 28 4d 61 74 72 69 78 |:Evaluat|e(Matrix|
|000012c0| 54 79 70 65 20 6d 74 29 | 0a 58 7b 0a 58 20 20 20 |Type mt)|.X{.X   |
|000012d0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|000012e0| 6d 3d 62 6d 2d 3e 45 76 | 61 6c 75 61 74 65 28 29 |m=bm->Ev|aluate()|
|000012f0| 3b 0a 58 20 20 20 69 6e | 74 20 6e 72 3d 67 6d 2d |;.X   in|t nr=gm-|
|00001300| 3e 4e 72 6f 77 73 28 29 | 3b 20 69 6e 74 20 6e 63 |>Nrows()|; int nc|
|00001310| 3d 67 6d 2d 3e 4e 63 6f | 6c 73 28 29 3b 0a 58 20 |=gm->Nco|ls();.X |
|00001320| 20 20 69 66 20 28 28 69 | 6e 74 29 6d 74 3d 3d 30 |  if ((i|nt)mt==0|
|00001330| 29 20 6d 74 20 3d 20 67 | 6d 2d 3e 54 79 70 65 28 |) mt = g|m->Type(|
|00001340| 29 3b 0a 58 20 20 20 69 | 66 20 28 6d 74 3d 3d 67 |);.X   i|f (mt==g|
|00001350| 6d 2d 3e 54 79 70 65 28 | 29 29 0a 58 20 20 20 7b |m->Type(|)).X   {|
|00001360| 0a 58 20 20 20 20 20 20 | 69 66 20 28 67 6d 2d 3e |.X      |if (gm->|
|00001370| 72 65 75 73 65 28 29 29 | 20 7b 20 52 45 50 4f 52 |reuse())| { REPOR|
|00001380| 54 20 67 6d 2d 3e 4d 75 | 6c 74 69 70 6c 79 28 66 |T gm->Mu|ltiply(f|
|00001390| 29 3b 20 72 65 74 75 72 | 6e 20 67 6d 3b 20 7d 0a |); retur|n gm; }.|
|000013a0| 58 20 20 20 20 20 20 65 | 6c 73 65 0a 58 20 20 20 |X      e|lse.X   |
|000013b0| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 20 20 20 52 |   {.X  |       R|
|000013c0| 45 50 4f 52 54 20 47 65 | 6e 65 72 61 6c 4d 61 74 |EPORT Ge|neralMat|
|000013d0| 72 69 78 2a 20 67 6d 78 | 20 3d 20 67 6d 2d 3e 54 |rix* gmx| = gm->T|
|000013e0| 79 70 65 28 29 2e 4e 65 | 77 28 6e 72 2c 6e 63 29 |ype().Ne|w(nr,nc)|
|000013f0| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 67 6d 78 2d |;.X     |    gmx-|
|00001400| 3e 52 65 6c 65 61 73 65 | 41 6e 64 44 65 6c 65 74 |>Release|AndDelet|
|00001410| 65 28 29 3b 20 67 6d 78 | 2d 3e 4d 75 6c 74 69 70 |e(); gmx|->Multip|
|00001420| 6c 79 28 67 6d 2c 66 29 | 3b 20 72 65 74 75 72 6e |ly(gm,f)|; return|
|00001430| 20 67 6d 78 3b 0a 58 20 | 20 20 20 20 20 7d 0a 58 | gmx;.X |     }.X|
|00001440| 20 20 20 7d 0a 58 20 20 | 20 65 6c 73 65 0a 58 20 |   }.X  | else.X |
|00001450| 20 20 7b 0a 58 20 20 20 | 20 20 20 52 45 50 4f 52 |  {.X   |   REPOR|
|00001460| 54 0a 58 20 20 20 20 20 | 20 47 65 6e 65 72 61 6c |T.X     | General|
|00001470| 4d 61 74 72 69 78 2a 20 | 67 6d 78 20 3d 20 6d 74 |Matrix* |gmx = mt|
|00001480| 2e 4e 65 77 28 6e 72 2c | 6e 63 29 3b 0a 58 20 20 |.New(nr,|nc);.X  |
|00001490| 20 20 20 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 |    Matr|ixRow mr|
|000014a0| 28 67 6d 2c 20 4c 6f 61 | 64 4f 6e 45 6e 74 72 79 |(gm, Loa|dOnEntry|
|000014b0| 29 3b 20 0a 58 20 20 20 | 20 20 20 4d 61 74 72 69 |); .X   |   Matri|
|000014c0| 78 52 6f 77 20 6d 72 78 | 28 67 6d 78 2c 20 53 74 |xRow mrx|(gmx, St|
|000014d0| 6f 72 65 4f 6e 45 78 69 | 74 2b 44 69 72 65 63 74 |oreOnExi|t+Direct|
|000014e0| 50 61 72 74 29 3b 0a 58 | 20 20 20 20 20 20 77 68 |Part);.X|      wh|
|000014f0| 69 6c 65 20 28 6e 72 2d | 2d 29 20 7b 20 6d 72 78 |ile (nr-|-) { mrx|
|00001500| 2e 4d 75 6c 74 69 70 6c | 79 28 6d 72 2c 66 29 3b |.Multipl|y(mr,f);|
|00001510| 20 6d 72 78 2e 4e 65 78 | 74 28 29 3b 20 6d 72 2e | mrx.Nex|t(); mr.|
|00001520| 4e 65 78 74 28 29 3b 20 | 7d 0a 58 20 20 20 20 20 |Next(); |}.X     |
|00001530| 20 67 6d 78 2d 3e 52 65 | 6c 65 61 73 65 41 6e 64 | gmx->Re|leaseAnd|
|00001540| 44 65 6c 65 74 65 28 29 | 3b 20 67 6d 2d 3e 74 44 |Delete()|; gm->tD|
|00001550| 65 6c 65 74 65 28 29 3b | 20 72 65 74 75 72 6e 20 |elete();| return |
|00001560| 67 6d 78 3b 0a 58 20 20 | 20 7d 0a 58 7d 0a 58 0a |gmx;.X  | }.X}.X.|
|00001570| 58 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 |XGeneral|Matrix* |
|00001580| 4e 65 67 61 74 65 64 4d | 61 74 72 69 78 3a 3a 45 |NegatedM|atrix::E|
|00001590| 76 61 6c 75 61 74 65 28 | 4d 61 74 72 69 78 54 79 |valuate(|MatrixTy|
|000015a0| 70 65 20 6d 74 29 0a 58 | 7b 0a 58 20 20 20 47 65 |pe mt).X|{.X   Ge|
|000015b0| 6e 65 72 61 6c 4d 61 74 | 72 69 78 2a 20 67 6d 3d |neralMat|rix* gm=|
|000015c0| 62 6d 2d 3e 45 76 61 6c | 75 61 74 65 28 29 3b 0a |bm->Eval|uate();.|
|000015d0| 58 20 20 20 69 66 20 28 | 28 69 6e 74 29 6d 74 3d |X   if (|(int)mt=|
|000015e0| 3d 30 29 20 6d 74 20 3d | 20 67 6d 2d 3e 54 79 70 |=0) mt =| gm->Typ|
|000015f0| 65 28 29 3b 0a 58 20 20 | 20 69 6e 74 20 6e 72 3d |e();.X  | int nr=|
|00001600| 67 6d 2d 3e 4e 72 6f 77 | 73 28 29 3b 20 69 6e 74 |gm->Nrow|s(); int|
|00001610| 20 6e 63 3d 67 6d 2d 3e | 4e 63 6f 6c 73 28 29 3b | nc=gm->|Ncols();|
|00001620| 0a 58 20 20 20 69 66 20 | 28 6d 74 3d 3d 67 6d 2d |.X   if |(mt==gm-|
|00001630| 3e 54 79 70 65 28 29 29 | 0a 58 20 20 20 7b 0a 58 |>Type())|.X   {.X|
|00001640| 20 20 20 20 20 20 69 66 | 20 28 67 6d 2d 3e 72 65 |      if| (gm->re|
|00001650| 75 73 65 28 29 29 20 7b | 20 52 45 50 4f 52 54 20 |use()) {| REPORT |
|00001660| 67 6d 2d 3e 4e 65 67 61 | 74 65 28 29 3b 20 72 65 |gm->Nega|te(); re|
|00001670| 74 75 72 6e 20 67 6d 3b | 20 7d 0a 58 20 20 20 20 |turn gm;| }.X    |
|00001680| 20 20 65 6c 73 65 0a 58 | 20 20 20 20 20 20 7b 0a |  else.X|      {.|
|00001690| 58 20 20 20 20 20 20 20 | 20 20 52 45 50 4f 52 54 |X       |  REPORT|
|000016a0| 0a 58 20 20 20 20 20 20 | 20 20 20 47 65 6e 65 72 |.X      |   Gener|
|000016b0| 61 6c 4d 61 74 72 69 78 | 2a 20 67 6d 78 20 3d 20 |alMatrix|* gmx = |
|000016c0| 67 6d 2d 3e 54 79 70 65 | 28 29 2e 4e 65 77 28 6e |gm->Type|().New(n|
|000016d0| 72 2c 6e 63 29 3b 20 67 | 6d 78 2d 3e 52 65 6c 65 |r,nc); g|mx->Rele|
|000016e0| 61 73 65 41 6e 64 44 65 | 6c 65 74 65 28 29 3b 0a |aseAndDe|lete();.|
|000016f0| 58 20 20 20 20 20 20 20 | 20 20 67 6d 78 2d 3e 4e |X       |  gmx->N|
|00001700| 65 67 61 74 65 28 67 6d | 29 3b 20 72 65 74 75 72 |egate(gm|); retur|
|00001710| 6e 20 67 6d 78 3b 0a 58 | 20 20 20 20 20 20 7d 0a |n gmx;.X|      }.|
|00001720| 58 20 20 20 7d 0a 58 20 | 20 20 65 6c 73 65 0a 58 |X   }.X |  else.X|
|00001730| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 52 45 50 4f |   {.X  |    REPO|
|00001740| 52 54 0a 58 20 20 20 20 | 20 20 47 65 6e 65 72 61 |RT.X    |  Genera|
|00001750| 6c 4d 61 74 72 69 78 2a | 20 67 6d 78 20 3d 20 6d |lMatrix*| gmx = m|
|00001760| 74 2e 4e 65 77 28 6e 72 | 2c 6e 63 29 3b 0a 58 20 |t.New(nr|,nc);.X |
|00001770| 20 20 20 20 20 4d 61 74 | 72 69 78 52 6f 77 20 6d |     Mat|rixRow m|
|00001780| 72 28 67 6d 2c 20 4c 6f | 61 64 4f 6e 45 6e 74 72 |r(gm, Lo|adOnEntr|
|00001790| 79 29 3b 20 0a 58 20 20 | 20 20 20 20 4d 61 74 72 |y); .X  |    Matr|
|000017a0| 69 78 52 6f 77 20 6d 72 | 78 28 67 6d 78 2c 20 53 |ixRow mr|x(gmx, S|
|000017b0| 74 6f 72 65 4f 6e 45 78 | 69 74 2b 44 69 72 65 63 |toreOnEx|it+Direc|
|000017c0| 74 50 61 72 74 29 3b 0a | 58 20 20 20 20 20 20 77 |tPart);.|X      w|
|000017d0| 68 69 6c 65 20 28 6e 72 | 2d 2d 29 20 7b 20 6d 72 |hile (nr|--) { mr|
|000017e0| 78 2e 4e 65 67 61 74 65 | 28 6d 72 29 3b 20 6d 72 |x.Negate|(mr); mr|
|000017f0| 78 2e 4e 65 78 74 28 29 | 3b 20 6d 72 2e 4e 65 78 |x.Next()|; mr.Nex|
|00001800| 74 28 29 3b 20 7d 0a 58 | 20 20 20 20 20 20 67 6d |t(); }.X|      gm|
|00001810| 78 2d 3e 52 65 6c 65 61 | 73 65 41 6e 64 44 65 6c |x->Relea|seAndDel|
|00001820| 65 74 65 28 29 3b 20 67 | 6d 2d 3e 74 44 65 6c 65 |ete(); g|m->tDele|
|00001830| 74 65 28 29 3b 20 72 65 | 74 75 72 6e 20 67 6d 78 |te(); re|turn gmx|
|00001840| 3b 0a 58 20 20 20 7d 0a | 58 7d 20 20 20 0a 58 0a |;.X   }.|X}   .X.|
|00001850| 58 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 |XGeneral|Matrix* |
|00001860| 54 72 61 6e 73 70 6f 73 | 65 64 4d 61 74 72 69 78 |Transpos|edMatrix|
|00001870| 3a 3a 45 76 61 6c 75 61 | 74 65 28 4d 61 74 72 69 |::Evalua|te(Matri|
|00001880| 78 54 79 70 65 20 6d 74 | 29 0a 58 7b 0a 58 20 20 |xType mt|).X{.X  |
|00001890| 20 52 45 50 4f 52 54 0a | 58 20 20 20 47 65 6e 65 | REPORT.|X   Gene|
|000018a0| 72 61 6c 4d 61 74 72 69 | 78 2a 20 67 6d 3d 62 6d |ralMatri|x* gm=bm|
|000018b0| 2d 3e 45 76 61 6c 75 61 | 74 65 28 29 3b 0a 58 20 |->Evalua|te();.X |
|000018c0| 20 20 69 66 20 28 28 69 | 6e 74 29 6d 74 3d 3d 30 |  if ((i|nt)mt==0|
|000018d0| 29 20 6d 74 20 3d 20 67 | 6d 2d 3e 54 79 70 65 28 |) mt = g|m->Type(|
|000018e0| 29 2e 74 28 29 3b 0a 58 | 20 20 20 47 65 6e 65 72 |).t();.X|   Gener|
|000018f0| 61 6c 4d 61 74 72 69 78 | 2a 20 67 6d 78 3d 67 6d |alMatrix|* gmx=gm|
|00001900| 2d 3e 54 72 61 6e 73 70 | 6f 73 65 28 6d 74 29 3b |->Transp|ose(mt);|
|00001910| 0a 58 20 20 20 72 65 74 | 75 72 6e 20 67 6d 78 3b |.X   ret|urn gmx;|
|00001920| 0a 58 7d 0a 58 20 20 20 | 0a 58 47 65 6e 65 72 61 |.X}.X   |.XGenera|
|00001930| 6c 4d 61 74 72 69 78 2a | 20 52 6f 77 65 64 4d 61 |lMatrix*| RowedMa|
|00001940| 74 72 69 78 3a 3a 45 76 | 61 6c 75 61 74 65 28 4d |trix::Ev|aluate(M|
|00001950| 61 74 72 69 78 54 79 70 | 65 20 6d 74 29 0a 58 7b |atrixTyp|e mt).X{|
|00001960| 0a 58 20 20 20 47 65 6e | 65 72 61 6c 4d 61 74 72 |.X   Gen|eralMatr|
|00001970| 69 78 2a 20 67 6d 20 3d | 20 62 6d 2d 3e 45 76 61 |ix* gm =| bm->Eva|
|00001980| 6c 75 61 74 65 28 29 3b | 0a 58 20 20 20 47 65 6e |luate();|.X   Gen|
|00001990| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 67 6d 78 20 |eralMatr|ix* gmx |
|000019a0| 3d 20 6e 65 77 20 52 6f | 77 56 65 63 74 6f 72 3b |= new Ro|wVector;|
|000019b0| 20 4d 61 74 72 69 78 45 | 72 72 6f 72 4e 6f 53 70 | MatrixE|rrorNoSp|
|000019c0| 61 63 65 28 67 6d 78 29 | 3b 0a 58 20 20 20 67 6d |ace(gmx)|;.X   gm|
|000019d0| 78 2d 3e 6e 72 6f 77 73 | 20 3d 20 31 3b 20 67 6d |x->nrows| = 1; gm|
|000019e0| 78 2d 3e 6e 63 6f 6c 73 | 20 3d 20 67 6d 78 2d 3e |x->ncols| = gmx->|
|000019f0| 73 74 6f 72 61 67 65 20 | 3d 20 67 6d 2d 3e 73 74 |storage |= gm->st|
|00001a00| 6f 72 61 67 65 3b 0a 58 | 20 20 20 72 65 74 75 72 |orage;.X|   retur|
|00001a10| 6e 20 67 6d 2d 3e 42 6f | 72 72 6f 77 53 74 6f 72 |n gm->Bo|rrowStor|
|00001a20| 65 28 67 6d 78 2c 6d 74 | 29 3b 0a 58 7d 0a 58 0a |e(gmx,mt|);.X}.X.|
|00001a30| 58 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 |XGeneral|Matrix* |
|00001a40| 43 6f 6c 65 64 4d 61 74 | 72 69 78 3a 3a 45 76 61 |ColedMat|rix::Eva|
|00001a50| 6c 75 61 74 65 28 4d 61 | 74 72 69 78 54 79 70 65 |luate(Ma|trixType|
|00001a60| 20 6d 74 29 0a 58 7b 0a | 58 20 20 20 47 65 6e 65 | mt).X{.|X   Gene|
|00001a70| 72 61 6c 4d 61 74 72 69 | 78 2a 20 67 6d 20 3d 20 |ralMatri|x* gm = |
|00001a80| 62 6d 2d 3e 45 76 61 6c | 75 61 74 65 28 29 3b 0a |bm->Eval|uate();.|
|00001a90| 58 20 20 20 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |X   Gene|ralMatri|
|00001aa0| 78 2a 20 67 6d 78 20 3d | 20 6e 65 77 20 43 6f 6c |x* gmx =| new Col|
|00001ab0| 75 6d 6e 56 65 63 74 6f | 72 3b 20 4d 61 74 72 69 |umnVecto|r; Matri|
|00001ac0| 78 45 72 72 6f 72 4e 6f | 53 70 61 63 65 28 67 6d |xErrorNo|Space(gm|
|00001ad0| 78 29 3b 0a 58 20 20 20 | 67 6d 78 2d 3e 6e 63 6f |x);.X   |gmx->nco|
|00001ae0| 6c 73 20 3d 20 31 3b 20 | 67 6d 78 2d 3e 6e 72 6f |ls = 1; |gmx->nro|
|00001af0| 77 73 20 3d 20 67 6d 78 | 2d 3e 73 74 6f 72 61 67 |ws = gmx|->storag|
|00001b00| 65 20 3d 20 67 6d 2d 3e | 73 74 6f 72 61 67 65 3b |e = gm->|storage;|
|00001b10| 0a 58 20 20 20 72 65 74 | 75 72 6e 20 67 6d 2d 3e |.X   ret|urn gm->|
|00001b20| 42 6f 72 72 6f 77 53 74 | 6f 72 65 28 67 6d 78 2c |BorrowSt|ore(gmx,|
|00001b30| 6d 74 29 3b 0a 58 7d 0a | 58 0a 58 47 65 6e 65 72 |mt);.X}.|X.XGener|
|00001b40| 61 6c 4d 61 74 72 69 78 | 2a 20 44 69 61 67 65 64 |alMatrix|* Diaged|
|00001b50| 4d 61 74 72 69 78 3a 3a | 45 76 61 6c 75 61 74 65 |Matrix::|Evaluate|
|00001b60| 28 4d 61 74 72 69 78 54 | 79 70 65 20 6d 74 29 0a |(MatrixT|ype mt).|
|00001b70| 58 7b 0a 58 20 20 20 47 | 65 6e 65 72 61 6c 4d 61 |X{.X   G|eneralMa|
|00001b80| 74 72 69 78 2a 20 67 6d | 20 3d 20 62 6d 2d 3e 45 |trix* gm| = bm->E|
|00001b90| 76 61 6c 75 61 74 65 28 | 29 3b 0a 58 20 20 20 47 |valuate(|);.X   G|
|00001ba0| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 20 67 6d |eneralMa|trix* gm|
|00001bb0| 78 20 3d 20 6e 65 77 20 | 44 69 61 67 6f 6e 61 6c |x = new |Diagonal|
|00001bc0| 4d 61 74 72 69 78 3b 20 | 4d 61 74 72 69 78 45 72 |Matrix; |MatrixEr|
|00001bd0| 72 6f 72 4e 6f 53 70 61 | 63 65 28 67 6d 78 29 3b |rorNoSpa|ce(gmx);|
|00001be0| 0a 58 20 20 20 67 6d 78 | 2d 3e 6e 72 6f 77 73 20 |.X   gmx|->nrows |
|00001bf0| 3d 20 67 6d 78 2d 3e 6e | 63 6f 6c 73 20 3d 20 67 |= gmx->n|cols = g|
|00001c00| 6d 78 2d 3e 73 74 6f 72 | 61 67 65 20 3d 20 67 6d |mx->stor|age = gm|
|00001c10| 2d 3e 73 74 6f 72 61 67 | 65 3b 0a 58 20 20 20 72 |->storag|e;.X   r|
|00001c20| 65 74 75 72 6e 20 67 6d | 2d 3e 42 6f 72 72 6f 77 |eturn gm|->Borrow|
|00001c30| 53 74 6f 72 65 28 67 6d | 78 2c 6d 74 29 3b 0a 58 |Store(gm|x,mt);.X|
|00001c40| 7d 0a 58 0a 58 47 65 6e | 65 72 61 6c 4d 61 74 72 |}.X.XGen|eralMatr|
|00001c50| 69 78 2a 20 4d 61 74 65 | 64 4d 61 74 72 69 78 3a |ix* Mate|dMatrix:|
|00001c60| 3a 45 76 61 6c 75 61 74 | 65 28 4d 61 74 72 69 78 |:Evaluat|e(Matrix|
|00001c70| 54 79 70 65 20 6d 74 29 | 0a 58 7b 0a 58 20 20 20 |Type mt)|.X{.X   |
|00001c80| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00001c90| 6d 20 3d 20 62 6d 2d 3e | 45 76 61 6c 75 61 74 65 |m = bm->|Evaluate|
|00001ca0| 28 29 3b 0a 58 20 20 20 | 47 65 6e 65 72 61 6c 4d |();.X   |GeneralM|
|00001cb0| 61 74 72 69 78 2a 20 67 | 6d 78 20 3d 20 6e 65 77 |atrix* g|mx = new|
|00001cc0| 20 4d 61 74 72 69 78 3b | 20 4d 61 74 72 69 78 45 | Matrix;| MatrixE|
|00001cd0| 72 72 6f 72 4e 6f 53 70 | 61 63 65 28 67 6d 78 29 |rrorNoSp|ace(gmx)|
|00001ce0| 3b 0a 58 20 20 20 67 6d | 78 2d 3e 6e 72 6f 77 73 |;.X   gm|x->nrows|
|00001cf0| 20 3d 20 6e 72 3b 20 67 | 6d 78 2d 3e 6e 63 6f 6c | = nr; g|mx->ncol|
|00001d00| 73 20 3d 20 6e 63 3b 20 | 67 6d 78 2d 3e 73 74 6f |s = nc; |gmx->sto|
|00001d10| 72 61 67 65 20 3d 20 67 | 6d 2d 3e 73 74 6f 72 61 |rage = g|m->stora|
|00001d20| 67 65 3b 0a 58 20 20 20 | 69 66 20 28 6e 72 2a 6e |ge;.X   |if (nr*n|
|00001d30| 63 20 21 3d 20 67 6d 78 | 2d 3e 73 74 6f 72 61 67 |c != gmx|->storag|
|00001d40| 65 29 0a 58 20 20 20 20 | 20 20 4d 61 74 72 69 78 |e).X    |  Matrix|
|00001d50| 45 72 72 6f 72 28 22 49 | 6e 63 6f 6d 70 61 74 69 |Error("I|ncompati|
|00001d60| 62 6c 65 20 64 69 6d 65 | 6e 73 69 6f 6e 73 20 69 |ble dime|nsions i|
|00001d70| 6e 20 43 6f 70 79 54 6f | 4d 61 74 72 69 78 22 29 |n CopyTo|Matrix")|
|00001d80| 3b 0a 58 20 20 20 72 65 | 74 75 72 6e 20 67 6d 2d |;.X   re|turn gm-|
|00001d90| 3e 42 6f 72 72 6f 77 53 | 74 6f 72 65 28 67 6d 78 |>BorrowS|tore(gmx|
|00001da0| 2c 6d 74 29 3b 0a 58 7d | 0a 58 0a 58 47 65 6e 65 |,mt);.X}|.X.XGene|
|00001db0| 72 61 6c 4d 61 74 72 69 | 78 2a 20 47 65 74 53 75 |ralMatri|x* GetSu|
|00001dc0| 62 4d 61 74 72 69 78 3a | 3a 45 76 61 6c 75 61 74 |bMatrix:|:Evaluat|
|00001dd0| 65 28 4d 61 74 72 69 78 | 54 79 70 65 20 6d 74 29 |e(Matrix|Type mt)|
|00001de0| 0a 58 7b 0a 58 20 20 20 | 52 45 50 4f 52 54 0a 58 |.X{.X   |REPORT.X|
|00001df0| 20 20 20 47 65 6e 65 72 | 61 6c 4d 61 74 72 69 78 |   Gener|alMatrix|
|00001e00| 2a 20 67 6d 20 3d 20 62 | 6d 2d 3e 45 76 61 6c 75 |* gm = b|m->Evalu|
|00001e10| 61 74 65 28 29 3b 0a 58 | 20 20 20 69 66 20 28 72 |ate();.X|   if (r|
|00001e20| 6f 77 5f 73 6b 69 70 2b | 72 6f 77 5f 6e 75 6d 62 |ow_skip+|row_numb|
|00001e30| 65 72 20 3e 20 67 6d 2d | 3e 4e 72 6f 77 73 28 29 |er > gm-|>Nrows()|
|00001e40| 20 7c 7c 20 63 6f 6c 5f | 73 6b 69 70 2b 63 6f 6c | || col_|skip+col|
|00001e50| 5f 6e 75 6d 62 65 72 20 | 3e 20 67 6d 2d 3e 4e 63 |_number |> gm->Nc|
|00001e60| 6f 6c 73 28 29 29 0a 58 | 20 20 20 20 20 20 4d 61 |ols()).X|      Ma|
|00001e70| 74 72 69 78 45 72 72 6f | 72 28 22 53 75 62 4d 61 |trixErro|r("SubMa|
|00001e80| 74 72 69 78 20 64 69 6d | 65 6e 73 69 6f 6e 20 65 |trix dim|ension e|
|00001e90| 72 72 6f 72 22 29 3b 0a | 58 20 20 20 69 66 20 28 |rror");.|X   if (|
|00001ea0| 69 6e 74 28 6d 74 29 3d | 3d 30 29 20 6d 74 20 3d |int(mt)=|=0) mt =|
|00001eb0| 20 4d 61 74 72 69 78 54 | 79 70 65 3a 3a 52 65 63 | MatrixT|ype::Rec|
|00001ec0| 74 3b 0a 58 20 20 20 47 | 65 6e 65 72 61 6c 4d 61 |t;.X   G|eneralMa|
|00001ed0| 74 72 69 78 2a 20 67 6d | 78 20 3d 20 6d 74 2e 4e |trix* gm|x = mt.N|
|00001ee0| 65 77 28 72 6f 77 5f 6e | 75 6d 62 65 72 2c 20 63 |ew(row_n|umber, c|
|00001ef0| 6f 6c 5f 6e 75 6d 62 65 | 72 29 3b 0a 58 20 20 20 |ol_numbe|r);.X   |
|00001f00| 69 6e 74 20 69 20 3d 20 | 72 6f 77 5f 6e 75 6d 62 |int i = |row_numb|
|00001f10| 65 72 3b 0a 58 20 20 20 | 4d 61 74 72 69 78 52 6f |er;.X   |MatrixRo|
|00001f20| 77 20 6d 72 28 67 6d 2c | 20 4c 6f 61 64 4f 6e 45 |w mr(gm,| LoadOnE|
|00001f30| 6e 74 72 79 2c 20 72 6f | 77 5f 73 6b 69 70 29 3b |ntry, ro|w_skip);|
|00001f40| 20 0a 58 20 20 20 4d 61 | 74 72 69 78 52 6f 77 20 | .X   Ma|trixRow |
|00001f50| 6d 72 78 28 67 6d 78 2c | 20 53 74 6f 72 65 4f 6e |mrx(gmx,| StoreOn|
|00001f60| 45 78 69 74 2b 44 69 72 | 65 63 74 50 61 72 74 29 |Exit+Dir|ectPart)|
|00001f70| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 52 6f 77 43 |;.X   Ma|trixRowC|
|00001f80| 6f 6c 20 73 75 62 3b 0a | 58 20 20 20 77 68 69 6c |ol sub;.|X   whil|
|00001f90| 65 20 28 69 2d 2d 29 0a | 58 20 20 20 7b 0a 58 20 |e (i--).|X   {.X |
|00001fa0| 20 20 20 20 20 6d 72 2e | 53 75 62 52 6f 77 43 6f |     mr.|SubRowCo|
|00001fb0| 6c 28 73 75 62 2c 20 63 | 6f 6c 5f 73 6b 69 70 2c |l(sub, c|ol_skip,|
|00001fc0| 20 63 6f 6c 5f 6e 75 6d | 62 65 72 29 3b 20 20 20 | col_num|ber);   |
|00001fd0| 2f 2f 20 70 75 74 20 76 | 61 6c 75 65 73 20 69 6e |// put v|alues in|
|00001fe0| 20 73 75 62 0a 58 20 20 | 20 20 20 20 6d 72 78 2e | sub.X  |    mrx.|
|00001ff0| 43 6f 70 79 28 73 75 62 | 29 3b 20 6d 72 78 2e 4e |Copy(sub|); mrx.N|
|00002000| 65 78 74 28 29 3b 20 6d | 72 2e 4e 65 78 74 28 29 |ext(); m|r.Next()|
|00002010| 3b 0a 58 20 20 20 7d 0a | 58 20 20 20 67 6d 78 2d |;.X   }.|X   gmx-|
|00002020| 3e 52 65 6c 65 61 73 65 | 41 6e 64 44 65 6c 65 74 |>Release|AndDelet|
|00002030| 65 28 29 3b 20 67 6d 2d | 3e 74 44 65 6c 65 74 65 |e(); gm-|>tDelete|
|00002040| 28 29 3b 20 72 65 74 75 | 72 6e 20 67 6d 78 3b 0a |(); retu|rn gmx;.|
|00002050| 58 7d 20 20 20 0a 58 0a | 58 47 65 6e 65 72 61 6c |X}   .X.|XGeneral|
|00002060| 4d 61 74 72 69 78 2a 20 | 52 65 74 75 72 6e 4d 61 |Matrix* |ReturnMa|
|00002070| 74 72 69 78 3a 3a 45 76 | 61 6c 75 61 74 65 28 4d |trix::Ev|aluate(M|
|00002080| 61 74 72 69 78 54 79 70 | 65 20 6d 74 29 0a 58 7b |atrixTyp|e mt).X{|
|00002090| 20 72 65 74 75 72 6e 20 | 67 6d 2d 3e 45 76 61 6c | return |gm->Eval|
|000020a0| 75 61 74 65 28 6d 74 29 | 3b 20 7d 0a 58 0a 58 0a |uate(mt)|; }.X.X.|
|000020b0| 58 76 6f 69 64 20 47 65 | 6e 65 72 61 6c 4d 61 74 |Xvoid Ge|neralMat|
|000020c0| 72 69 78 3a 3a 41 64 64 | 28 47 65 6e 65 72 61 6c |rix::Add|(General|
|000020d0| 4d 61 74 72 69 78 2a 20 | 67 6d 31 2c 20 72 65 61 |Matrix* |gm1, rea|
|000020e0| 6c 20 66 29 0a 58 7b 0a | 58 20 20 20 52 45 50 4f |l f).X{.|X   REPO|
|000020f0| 52 54 0a 58 20 20 20 72 | 65 67 69 73 74 65 72 20 |RT.X   r|egister |
|00002100| 72 65 61 6c 2a 20 73 31 | 3d 67 6d 31 2d 3e 73 74 |real* s1|=gm1->st|
|00002110| 6f 72 65 3b 0a 58 20 20 | 20 72 65 67 69 73 74 65 |ore;.X  | registe|
|00002120| 72 20 72 65 61 6c 2a 20 | 73 3d 73 74 6f 72 65 3b |r real* |s=store;|
|00002130| 20 72 65 67 69 73 74 65 | 72 20 69 6e 74 20 69 3d | registe|r int i=|
|00002140| 73 74 6f 72 61 67 65 3b | 0a 58 20 20 20 77 68 69 |storage;|.X   whi|
|00002150| 6c 65 20 28 69 2d 2d 29 | 20 2a 73 2b 2b 20 3d 20 |le (i--)| *s++ = |
|00002160| 2a 73 31 2b 2b 20 2b 20 | 66 3b 0a 58 7d 0a 58 20 |*s1++ + |f;.X}.X |
|00002170| 20 20 0a 58 76 6f 69 64 | 20 47 65 6e 65 72 61 6c |  .Xvoid| General|
|00002180| 4d 61 74 72 69 78 3a 3a | 41 64 64 28 72 65 61 6c |Matrix::|Add(real|
|00002190| 20 66 29 0a 58 7b 0a 58 | 20 20 20 52 45 50 4f 52 | f).X{.X|   REPOR|
|000021a0| 54 0a 58 20 20 20 72 65 | 67 69 73 74 65 72 20 72 |T.X   re|gister r|
|000021b0| 65 61 6c 2a 20 73 3d 73 | 74 6f 72 65 3b 20 72 65 |eal* s=s|tore; re|
|000021c0| 67 69 73 74 65 72 20 69 | 6e 74 20 69 3d 73 74 6f |gister i|nt i=sto|
|000021d0| 72 61 67 65 3b 20 77 68 | 69 6c 65 20 28 69 2d 2d |rage; wh|ile (i--|
|000021e0| 29 20 2a 73 2b 2b 20 2b | 3d 20 66 3b 0a 58 7d 0a |) *s++ +|= f;.X}.|
|000021f0| 58 20 20 20 0a 58 76 6f | 69 64 20 47 65 6e 65 72 |X   .Xvo|id Gener|
|00002200| 61 6c 4d 61 74 72 69 78 | 3a 3a 4e 65 67 61 74 65 |alMatrix|::Negate|
|00002210| 28 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 |(General|Matrix* |
|00002220| 67 6d 31 29 0a 58 7b 0a | 58 20 20 20 2f 2f 20 63 |gm1).X{.|X   // c|
|00002230| 68 61 6e 67 65 20 73 69 | 67 6e 20 6f 66 20 65 6c |hange si|gn of el|
|00002240| 65 6d 65 6e 74 73 0a 58 | 20 20 20 52 45 50 4f 52 |ements.X|   REPOR|
|00002250| 54 0a 58 20 20 20 72 65 | 61 6c 2a 20 73 31 3d 67 |T.X   re|al* s1=g|
|00002260| 6d 31 2d 3e 73 74 6f 72 | 65 3b 20 72 65 61 6c 2a |m1->stor|e; real*|
|00002270| 20 73 3d 73 74 6f 72 65 | 3b 20 72 65 67 69 73 74 | s=store|; regist|
|00002280| 65 72 20 69 6e 74 20 69 | 3d 73 74 6f 72 61 67 65 |er int i|=storage|
|00002290| 3b 0a 58 20 20 20 77 68 | 69 6c 65 28 69 2d 2d 29 |;.X   wh|ile(i--)|
|000022a0| 20 2a 73 2b 2b 20 3d 20 | 2d 28 2a 73 31 2b 2b 29 | *s++ = |-(*s1++)|
|000022b0| 3b 0a 58 7d 0a 58 20 20 | 20 0a 58 76 6f 69 64 20 |;.X}.X  | .Xvoid |
|000022c0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 3a 3a 4e |GeneralM|atrix::N|
|000022d0| 65 67 61 74 65 28 29 0a | 58 7b 0a 58 20 20 20 52 |egate().|X{.X   R|
|000022e0| 45 50 4f 52 54 0a 58 20 | 20 20 72 65 61 6c 2a 20 |EPORT.X |  real* |
|000022f0| 73 3d 73 74 6f 72 65 3b | 20 72 65 67 69 73 74 65 |s=store;| registe|
|00002300| 72 20 69 6e 74 20 69 3d | 73 74 6f 72 61 67 65 3b |r int i=|storage;|
|00002310| 20 77 68 69 6c 65 28 69 | 2d 2d 29 20 7b 20 2a 73 | while(i|--) { *s|
|00002320| 20 3d 20 2d 28 2a 73 29 | 3b 20 73 2b 2b 3b 20 7d | = -(*s)|; s++; }|
|00002330| 20 0a 58 7d 0a 58 20 20 | 20 0a 58 76 6f 69 64 20 | .X}.X  | .Xvoid |
|00002340| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 3a 3a 4d |GeneralM|atrix::M|
|00002350| 75 6c 74 69 70 6c 79 28 | 47 65 6e 65 72 61 6c 4d |ultiply(|GeneralM|
|00002360| 61 74 72 69 78 2a 20 67 | 6d 31 2c 20 72 65 61 6c |atrix* g|m1, real|
|00002370| 20 66 29 0a 58 7b 0a 58 | 20 20 20 52 45 50 4f 52 | f).X{.X|   REPOR|
|00002380| 54 0a 58 20 20 20 72 65 | 67 69 73 74 65 72 20 72 |T.X   re|gister r|
|00002390| 65 61 6c 2a 20 73 31 3d | 67 6d 31 2d 3e 73 74 6f |eal* s1=|gm1->sto|
|000023a0| 72 65 3b 0a 58 20 20 20 | 72 65 67 69 73 74 65 72 |re;.X   |register|
|000023b0| 20 72 65 61 6c 2a 20 73 | 3d 73 74 6f 72 65 3b 20 | real* s|=store; |
|000023c0| 72 65 67 69 73 74 65 72 | 20 69 6e 74 20 69 3d 73 |register| int i=s|
|000023d0| 74 6f 72 61 67 65 3b 0a | 58 20 20 20 77 68 69 6c |torage;.|X   whil|
|000023e0| 65 20 28 69 2d 2d 29 20 | 2a 73 2b 2b 20 3d 20 2a |e (i--) |*s++ = *|
|000023f0| 73 31 2b 2b 20 2a 20 66 | 3b 0a 58 7d 0a 58 20 20 |s1++ * f|;.X}.X  |
|00002400| 20 0a 58 76 6f 69 64 20 | 47 65 6e 65 72 61 6c 4d | .Xvoid |GeneralM|
|00002410| 61 74 72 69 78 3a 3a 4d | 75 6c 74 69 70 6c 79 28 |atrix::M|ultiply(|
|00002420| 72 65 61 6c 20 66 29 0a | 58 7b 0a 58 20 20 20 52 |real f).|X{.X   R|
|00002430| 45 50 4f 52 54 0a 58 20 | 20 20 72 65 67 69 73 74 |EPORT.X |  regist|
|00002440| 65 72 20 72 65 61 6c 2a | 20 73 3d 73 74 6f 72 65 |er real*| s=store|
|00002450| 3b 20 72 65 67 69 73 74 | 65 72 20 69 6e 74 20 69 |; regist|er int i|
|00002460| 3d 73 74 6f 72 61 67 65 | 3b 20 77 68 69 6c 65 20 |=storage|; while |
|00002470| 28 69 2d 2d 29 20 2a 73 | 2b 2b 20 2a 3d 20 66 3b |(i--) *s|++ *= f;|
|00002480| 0a 58 7d 0a 58 20 20 20 | 0a 45 4e 44 5f 4f 46 5f |.X}.X   |.END_OF_|
|00002490| 46 49 4c 45 0a 69 66 20 | 74 65 73 74 20 37 37 37 |FILE.if |test 777|
|000024a0| 35 20 2d 6e 65 20 60 77 | 63 20 2d 63 20 3c 27 6e |5 -ne `w|c -c <'n|
|000024b0| 65 77 6d 61 74 35 2e 63 | 78 78 27 60 3b 20 74 68 |ewmat5.c|xx'`; th|
|000024c0| 65 6e 0a 20 20 20 20 65 | 63 68 6f 20 73 68 61 72 |en.    e|cho shar|
|000024d0| 3a 20 5c 22 27 6e 65 77 | 6d 61 74 35 2e 63 78 78 |: \"'new|mat5.cxx|
|000024e0| 27 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |'\" unpa|cked wit|
|000024f0| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 66 69 |h wrong |size!.fi|
|00002500| 0a 23 20 65 6e 64 20 6f | 66 20 27 6e 65 77 6d 61 |.# end o|f 'newma|
|00002510| 74 35 2e 63 78 78 27 0a | 66 69 0a 69 66 20 74 65 |t5.cxx'.|fi.if te|
|00002520| 73 74 20 2d 66 20 27 6e | 65 77 6d 61 74 36 2e 63 |st -f 'n|ewmat6.c|
|00002530| 78 78 27 20 2d 61 20 22 | 24 7b 31 7d 22 20 21 3d |xx' -a "|${1}" !=|
|00002540| 20 22 2d 63 22 20 3b 20 | 74 68 65 6e 20 0a 20 20 | "-c" ; |then .  |
|00002550| 65 63 68 6f 20 73 68 61 | 72 3a 20 57 69 6c 6c 20 |echo sha|r: Will |
|00002560| 6e 6f 74 20 63 6c 6f 62 | 62 65 72 20 65 78 69 73 |not clob|ber exis|
|00002570| 74 69 6e 67 20 66 69 6c | 65 20 5c 22 27 6e 65 77 |ting fil|e \"'new|
|00002580| 6d 61 74 36 2e 63 78 78 | 27 5c 22 0a 65 6c 73 65 |mat6.cxx|'\".else|
|00002590| 0a 65 63 68 6f 20 73 68 | 61 72 3a 20 45 78 74 72 |.echo sh|ar: Extr|
|000025a0| 61 63 74 69 6e 67 20 5c | 22 27 6e 65 77 6d 61 74 |acting \|"'newmat|
|000025b0| 36 2e 63 78 78 27 5c 22 | 20 5c 28 37 30 38 34 20 |6.cxx'\"| \(7084 |
|000025c0| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|000025d0| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 6e 65 77 6d | "s/^X//|" >'newm|
|000025e0| 61 74 36 2e 63 78 78 27 | 20 3c 3c 27 45 4e 44 5f |at6.cxx'| <<'END_|
|000025f0| 4f 46 5f 46 49 4c 45 27 | 0a 58 2f 2f 24 24 20 6e |OF_FILE'|.X//$$ n|
|00002600| 65 77 6d 61 74 36 2e 63 | 78 78 20 20 20 20 20 20 |ewmat6.c|xx      |
|00002610| 20 20 20 20 20 20 4f 70 | 65 72 61 74 6f 72 73 2c |      Op|erators,|
|00002620| 20 65 6c 65 6d 65 6e 74 | 20 61 63 63 65 73 73 2c | element| access,|
|00002630| 20 73 75 62 6d 61 74 72 | 69 63 65 73 0a 58 0a 58 | submatr|ices.X.X|
|00002640| 2f 2f 20 43 6f 70 79 72 | 69 67 68 74 20 28 43 29 |// Copyr|ight (C)|
|00002650| 20 31 39 39 31 3a 20 52 | 20 42 20 44 61 76 69 65 | 1991: R| B Davie|
|00002660| 73 20 61 6e 64 20 44 53 | 49 52 0a 58 0a 58 23 69 |s and DS|IR.X.X#i|
|00002670| 6e 63 6c 75 64 65 20 22 | 69 6e 63 6c 75 64 65 2e |nclude "|include.|
|00002680| 68 78 78 22 0a 58 0a 58 | 23 69 6e 63 6c 75 64 65 |hxx".X.X|#include|
|00002690| 20 22 6e 65 77 6d 61 74 | 2e 68 78 78 22 0a 58 23 | "newmat|.hxx".X#|
|000026a0| 69 6e 63 6c 75 64 65 20 | 22 6e 65 77 6d 61 74 72 |include |"newmatr|
|000026b0| 63 2e 68 78 78 22 0a 58 | 0a 58 0a 58 2f 2f 23 64 |c.hxx".X|.X.X//#d|
|000026c0| 65 66 69 6e 65 20 52 45 | 50 4f 52 54 20 7b 20 73 |efine RE|PORT { s|
|000026d0| 74 61 74 69 63 20 45 78 | 65 43 6f 75 6e 74 65 72 |tatic Ex|eCounter|
|000026e0| 20 45 78 65 43 6f 75 6e | 74 28 5f 5f 4c 49 4e 45 | ExeCoun|t(__LINE|
|000026f0| 5f 5f 2c 36 29 3b 20 45 | 78 65 43 6f 75 6e 74 2b |__,6); E|xeCount+|
|00002700| 2b 3b 20 7d 0a 58 0a 58 | 23 64 65 66 69 6e 65 20 |+; }.X.X|#define |
|00002710| 52 45 50 4f 52 54 20 7b | 7d 0a 58 0a 58 2f 2a 2a |REPORT {|}.X.X/**|
|00002720| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002730| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 20 67 65 6e 65 72 61 |********|* genera|
|00002740| 6c 20 75 74 69 6c 69 74 | 69 65 73 20 2a 2a 2a 2a |l utilit|ies ****|
|00002750| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002760| 2a 2a 2a 2a 2a 2f 0a 58 | 0a 58 73 74 61 74 69 63 |*****/.X|.Xstatic|
|00002770| 20 69 6e 74 20 74 72 69 | 73 74 6f 72 65 28 69 6e | int tri|store(in|
|00002780| 74 20 6e 29 20 20 20 20 | 20 20 20 20 20 20 20 20 |t n)    |        |
|00002790| 20 20 20 20 20 20 20 20 | 20 20 2f 2f 20 65 6c 73 |        |  // els|
|000027a0| 20 69 6e 20 74 72 69 61 | 6e 67 75 6c 61 72 20 6d | in tria|ngular m|
|000027b0| 61 74 72 69 78 0a 58 7b | 20 72 65 74 75 72 6e 20 |atrix.X{| return |
|000027c0| 28 6e 2a 28 6e 2b 31 29 | 29 2f 32 3b 20 7d 0a 58 |(n*(n+1)|)/2; }.X|
|000027d0| 0a 58 0a 58 2f 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X.X/***|********|
|000027e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000027f0| 2a 2a 2a 20 6f 70 65 72 | 61 74 6f 72 73 20 2a 2a |*** oper|ators **|
|00002800| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002810| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2f 0a 58 |********|*****/.X|
|00002820| 0a 58 72 65 61 6c 26 20 | 4d 61 74 72 69 78 3a 3a |.Xreal& |Matrix::|
|00002830| 6f 70 65 72 61 74 6f 72 | 28 29 28 69 6e 74 20 6d |operator|()(int m|
|00002840| 2c 20 69 6e 74 20 6e 29 | 0a 58 7b 0a 58 20 20 20 |, int n)|.X{.X   |
|00002850| 69 66 20 28 6d 3c 3d 30 | 20 7c 7c 20 6d 3e 6e 72 |if (m<=0| || m>nr|
|00002860| 6f 77 73 20 7c 7c 20 6e | 3c 3d 30 20 7c 7c 20 6e |ows || n|<=0 || n|
|00002870| 3e 6e 63 6f 6c 73 29 20 | 4d 61 74 72 69 78 45 72 |>ncols) |MatrixEr|
|00002880| 72 6f 72 28 22 49 6e 64 | 65 78 20 6f 75 74 20 6f |ror("Ind|ex out o|
|00002890| 66 20 72 61 6e 67 65 22 | 29 3b 0a 58 20 20 20 72 |f range"|);.X   r|
|000028a0| 65 74 75 72 6e 20 73 74 | 6f 72 65 5b 28 6d 2d 31 |eturn st|ore[(m-1|
|000028b0| 29 2a 6e 63 6f 6c 73 2b | 6e 2d 31 5d 3b 0a 58 7d |)*ncols+|n-1];.X}|
|000028c0| 0a 58 0a 58 72 65 61 6c | 26 20 53 79 6d 6d 65 74 |.X.Xreal|& Symmet|
|000028d0| 72 69 63 4d 61 74 72 69 | 78 3a 3a 6f 70 65 72 61 |ricMatri|x::opera|
|000028e0| 74 6f 72 28 29 28 69 6e | 74 20 6d 2c 20 69 6e 74 |tor()(in|t m, int|
|000028f0| 20 6e 29 0a 58 7b 0a 58 | 20 20 20 69 66 20 28 6d | n).X{.X|   if (m|
|00002900| 3c 3d 30 20 7c 7c 20 6e | 3c 3d 30 20 7c 7c 20 6d |<=0 || n|<=0 || m|
|00002910| 3e 6e 72 6f 77 73 20 7c | 7c 20 6e 3e 6e 63 6f 6c |>nrows ||| n>ncol|
|00002920| 73 29 20 4d 61 74 72 69 | 78 45 72 72 6f 72 28 22 |s) Matri|xError("|
|00002930| 49 6e 64 65 78 20 6f 75 | 74 20 6f 66 20 72 61 6e |Index ou|t of ran|
|00002940| 67 65 22 29 3b 0a 58 20 | 20 20 69 66 20 28 6d 3e |ge");.X |  if (m>|
|00002950| 3d 6e 29 20 72 65 74 75 | 72 6e 20 73 74 6f 72 65 |=n) retu|rn store|
|00002960| 5b 74 72 69 73 74 6f 72 | 65 28 6d 2d 31 29 2b 6e |[tristor|e(m-1)+n|
|00002970| 2d 31 5d 3b 0a 58 20 20 | 20 65 6c 73 65 20 72 65 |-1];.X  | else re|
|00002980| 74 75 72 6e 20 73 74 6f | 72 65 5b 74 72 69 73 74 |turn sto|re[trist|
|00002990| 6f 72 65 28 6e 2d 31 29 | 2b 6d 2d 31 5d 3b 0a 58 |ore(n-1)|+m-1];.X|
|000029a0| 7d 0a 58 0a 58 72 65 61 | 6c 26 20 55 70 70 65 72 |}.X.Xrea|l& Upper|
|000029b0| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|000029c0| 3a 3a 6f 70 65 72 61 74 | 6f 72 28 29 28 69 6e 74 |::operat|or()(int|
|000029d0| 20 6d 2c 20 69 6e 74 20 | 6e 29 0a 58 7b 0a 58 20 | m, int |n).X{.X |
|000029e0| 20 20 69 66 20 28 6d 3c | 3d 30 20 7c 7c 20 6e 3c |  if (m<|=0 || n<|
|000029f0| 6d 20 7c 7c 20 6e 3e 6e | 63 6f 6c 73 29 20 4d 61 |m || n>n|cols) Ma|
|00002a00| 74 72 69 78 45 72 72 6f | 72 28 22 49 6e 64 65 78 |trixErro|r("Index|
|00002a10| 20 6f 75 74 20 6f 66 20 | 72 61 6e 67 65 22 29 3b | out of |range");|
|00002a20| 0a 58 20 20 20 72 65 74 | 75 72 6e 20 73 74 6f 72 |.X   ret|urn stor|
|00002a30| 65 5b 28 6d 2d 31 29 2a | 6e 63 6f 6c 73 2b 6e 2d |e[(m-1)*|ncols+n-|
|00002a40| 31 2d 74 72 69 73 74 6f | 72 65 28 6d 2d 31 29 5d |1-tristo|re(m-1)]|
|00002a50| 3b 0a 58 7d 0a 58 0a 58 | 72 65 61 6c 26 20 4c 6f |;.X}.X.X|real& Lo|
|00002a60| 77 65 72 54 72 69 61 6e | 67 75 6c 61 72 4d 61 74 |werTrian|gularMat|
|00002a70| 72 69 78 3a 3a 6f 70 65 | 72 61 74 6f 72 28 29 28 |rix::ope|rator()(|
|00002a80| 69 6e 74 20 6d 2c 20 69 | 6e 74 20 6e 29 0a 58 7b |int m, i|nt n).X{|
|00002a90| 0a 58 20 20 20 69 66 20 | 28 6e 3c 3d 30 20 7c 7c |.X   if |(n<=0 |||
|00002aa0| 20 6d 3c 6e 20 7c 7c 20 | 6d 3e 6e 72 6f 77 73 29 | m<n || |m>nrows)|
|00002ab0| 20 4d 61 74 72 69 78 45 | 72 72 6f 72 28 22 49 6e | MatrixE|rror("In|
|00002ac0| 64 65 78 20 6f 75 74 20 | 6f 66 20 72 61 6e 67 65 |dex out |of range|
|00002ad0| 22 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 73 |");.X   |return s|
|00002ae0| 74 6f 72 65 5b 74 72 69 | 73 74 6f 72 65 28 6d 2d |tore[tri|store(m-|
|00002af0| 31 29 2b 6e 2d 31 5d 3b | 0a 58 7d 0a 58 0a 58 72 |1)+n-1];|.X}.X.Xr|
|00002b00| 65 61 6c 26 20 44 69 61 | 67 6f 6e 61 6c 4d 61 74 |eal& Dia|gonalMat|
|00002b10| 72 69 78 3a 3a 6f 70 65 | 72 61 74 6f 72 28 29 28 |rix::ope|rator()(|
|00002b20| 69 6e 74 20 6d 2c 20 69 | 6e 74 20 6e 29 0a 58 7b |int m, i|nt n).X{|
|00002b30| 0a 58 20 20 20 69 66 20 | 28 6e 3c 3d 30 20 7c 7c |.X   if |(n<=0 |||
|00002b40| 20 6d 21 3d 6e 20 7c 7c | 20 6d 3e 6e 72 6f 77 73 | m!=n ||| m>nrows|
|00002b50| 20 7c 7c 20 6e 3e 6e 63 | 6f 6c 73 29 20 4d 61 74 | || n>nc|ols) Mat|
|00002b60| 72 69 78 45 72 72 6f 72 | 28 22 49 6e 64 65 78 20 |rixError|("Index |
|00002b70| 6f 75 74 20 6f 66 20 72 | 61 6e 67 65 22 29 3b 0a |out of r|ange");.|
|00002b80| 58 20 20 20 72 65 74 75 | 72 6e 20 73 74 6f 72 65 |X   retu|rn store|
|00002b90| 5b 6e 2d 31 5d 3b 0a 58 | 7d 0a 58 0a 58 72 65 61 |[n-1];.X|}.X.Xrea|
|00002ba0| 6c 26 20 44 69 61 67 6f | 6e 61 6c 4d 61 74 72 69 |l& Diago|nalMatri|
|00002bb0| 78 3a 3a 6f 70 65 72 61 | 74 6f 72 28 29 28 69 6e |x::opera|tor()(in|
|00002bc0| 74 20 6d 29 0a 58 7b 0a | 58 20 20 20 69 66 20 28 |t m).X{.|X   if (|
|00002bd0| 6d 3c 3d 30 20 7c 7c 20 | 6d 3e 6e 72 6f 77 73 29 |m<=0 || |m>nrows)|
|00002be0| 20 4d 61 74 72 69 78 45 | 72 72 6f 72 28 22 49 6e | MatrixE|rror("In|
|00002bf0| 64 65 78 20 6f 75 74 20 | 6f 66 20 72 61 6e 67 65 |dex out |of range|
|00002c00| 22 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 73 |");.X   |return s|
|00002c10| 74 6f 72 65 5b 6d 2d 31 | 5d 3b 0a 58 7d 0a 58 0a |tore[m-1|];.X}.X.|
|00002c20| 58 72 65 61 6c 26 20 43 | 6f 6c 75 6d 6e 56 65 63 |Xreal& C|olumnVec|
|00002c30| 74 6f 72 3a 3a 6f 70 65 | 72 61 74 6f 72 28 29 28 |tor::ope|rator()(|
|00002c40| 69 6e 74 20 6d 29 0a 58 | 7b 0a 58 20 20 20 69 66 |int m).X|{.X   if|
|00002c50| 20 28 6d 3c 3d 30 20 7c | 7c 20 6d 3e 20 6e 72 6f | (m<=0 ||| m> nro|
|00002c60| 77 73 29 20 4d 61 74 72 | 69 78 45 72 72 6f 72 28 |ws) Matr|ixError(|
|00002c70| 22 49 6e 64 65 78 20 6f | 75 74 20 6f 66 20 72 61 |"Index o|ut of ra|
|00002c80| 6e 67 65 22 29 3b 0a 58 | 20 20 20 72 65 74 75 72 |nge");.X|   retur|
|00002c90| 6e 20 73 74 6f 72 65 5b | 6d 2d 31 5d 3b 0a 58 7d |n store[|m-1];.X}|
|00002ca0| 0a 58 0a 58 72 65 61 6c | 26 20 52 6f 77 56 65 63 |.X.Xreal|& RowVec|
|00002cb0| 74 6f 72 3a 3a 6f 70 65 | 72 61 74 6f 72 28 29 28 |tor::ope|rator()(|
|00002cc0| 69 6e 74 20 6e 29 0a 58 | 7b 0a 58 20 20 20 69 66 |int n).X|{.X   if|
|00002cd0| 20 28 6e 3c 3d 30 20 7c | 7c 20 6e 3e 20 6e 63 6f | (n<=0 ||| n> nco|
|00002ce0| 6c 73 29 20 4d 61 74 72 | 69 78 45 72 72 6f 72 28 |ls) Matr|ixError(|
|00002cf0| 22 49 6e 64 65 78 20 6f | 75 74 20 6f 66 20 72 61 |"Index o|ut of ra|
|00002d00| 6e 67 65 22 29 3b 0a 58 | 20 20 20 72 65 74 75 72 |nge");.X|   retur|
|00002d10| 6e 20 73 74 6f 72 65 5b | 6e 2d 31 5d 3b 0a 58 7d |n store[|n-1];.X}|
|00002d20| 0a 58 0a 58 23 69 66 6e | 64 65 66 20 5f 5f 5a 54 |.X.X#ifn|def __ZT|
|00002d30| 43 5f 5f 0a 58 0a 58 72 | 65 61 6c 20 4d 61 74 72 |C__.X.Xr|eal Matr|
|00002d40| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 28 29 28 69 |ix::oper|ator()(i|
|00002d50| 6e 74 20 6d 2c 20 69 6e | 74 20 6e 29 20 63 6f 6e |nt m, in|t n) con|
|00002d60| 73 74 0a 58 7b 0a 58 20 | 20 20 69 66 20 28 6d 3c |st.X{.X |  if (m<|
|00002d70| 3d 30 20 7c 7c 20 6d 3e | 6e 72 6f 77 73 20 7c 7c |=0 || m>|nrows |||
|00002d80| 20 6e 3c 3d 30 20 7c 7c | 20 6e 3e 6e 63 6f 6c 73 | n<=0 ||| n>ncols|
|00002d90| 29 20 4d 61 74 72 69 78 | 45 72 72 6f 72 28 22 49 |) Matrix|Error("I|
|00002da0| 6e 64 65 78 20 6f 75 74 | 20 6f 66 20 72 61 6e 67 |ndex out| of rang|
|00002db0| 65 22 29 3b 0a 58 20 20 | 20 72 65 74 75 72 6e 20 |e");.X  | return |
|00002dc0| 73 74 6f 72 65 5b 28 6d | 2d 31 29 2a 6e 63 6f 6c |store[(m|-1)*ncol|
|00002dd0| 73 2b 6e 2d 31 5d 3b 0a | 58 7d 0a 58 0a 58 72 65 |s+n-1];.|X}.X.Xre|
|00002de0| 61 6c 20 53 79 6d 6d 65 | 74 72 69 63 4d 61 74 72 |al Symme|tricMatr|
|00002df0| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 28 29 28 69 |ix::oper|ator()(i|
|00002e00| 6e 74 20 6d 2c 20 69 6e | 74 20 6e 29 20 63 6f 6e |nt m, in|t n) con|
|00002e10| 73 74 0a 58 7b 0a 58 20 | 20 20 69 66 20 28 6d 3c |st.X{.X |  if (m<|
|00002e20| 3d 30 20 7c 7c 20 6e 3c | 3d 30 20 7c 7c 20 6d 3e |=0 || n<|=0 || m>|
|00002e30| 6e 72 6f 77 73 20 7c 7c | 20 6e 3e 6e 63 6f 6c 73 |nrows ||| n>ncols|
|00002e40| 29 20 4d 61 74 72 69 78 | 45 72 72 6f 72 28 22 49 |) Matrix|Error("I|
|00002e50| 6e 64 65 78 20 6f 75 74 | 20 6f 66 20 72 61 6e 67 |ndex out| of rang|
|00002e60| 65 22 29 3b 0a 58 20 20 | 20 69 66 20 28 6d 3e 3d |e");.X  | if (m>=|
|00002e70| 6e 29 20 72 65 74 75 72 | 6e 20 73 74 6f 72 65 5b |n) retur|n store[|
|00002e80| 74 72 69 73 74 6f 72 65 | 28 6d 2d 31 29 2b 6e 2d |tristore|(m-1)+n-|
|00002e90| 31 5d 3b 0a 58 20 20 20 | 65 6c 73 65 20 72 65 74 |1];.X   |else ret|
|00002ea0| 75 72 6e 20 73 74 6f 72 | 65 5b 74 72 69 73 74 6f |urn stor|e[tristo|
|00002eb0| 72 65 28 6e 2d 31 29 2b | 6d 2d 31 5d 3b 0a 58 7d |re(n-1)+|m-1];.X}|
|00002ec0| 0a 58 0a 58 72 65 61 6c | 20 55 70 70 65 72 54 72 |.X.Xreal| UpperTr|
|00002ed0| 69 61 6e 67 75 6c 61 72 | 4d 61 74 72 69 78 3a 3a |iangular|Matrix::|
|00002ee0| 6f 70 65 72 61 74 6f 72 | 28 29 28 69 6e 74 20 6d |operator|()(int m|
|00002ef0| 2c 20 69 6e 74 20 6e 29 | 20 63 6f 6e 73 74 0a 58 |, int n)| const.X|
|00002f00| 7b 0a 58 20 20 20 69 66 | 20 28 6d 3c 3d 30 20 7c |{.X   if| (m<=0 ||
|00002f10| 7c 20 6e 3c 6d 20 7c 7c | 20 6e 3e 6e 63 6f 6c 73 || n<m ||| n>ncols|
|00002f20| 29 20 4d 61 74 72 69 78 | 45 72 72 6f 72 28 22 49 |) Matrix|Error("I|
|00002f30| 6e 64 65 78 20 6f 75 74 | 20 6f 66 20 72 61 6e 67 |ndex out| of rang|
|00002f40| 65 22 29 3b 0a 58 20 20 | 20 72 65 74 75 72 6e 20 |e");.X  | return |
|00002f50| 73 74 6f 72 65 5b 28 6d | 2d 31 29 2a 6e 63 6f 6c |store[(m|-1)*ncol|
|00002f60| 73 2b 6e 2d 31 2d 74 72 | 69 73 74 6f 72 65 28 6d |s+n-1-tr|istore(m|
|00002f70| 2d 31 29 5d 3b 0a 58 7d | 0a 58 0a 58 72 65 61 6c |-1)];.X}|.X.Xreal|
|00002f80| 20 4c 6f 77 65 72 54 72 | 69 61 6e 67 75 6c 61 72 | LowerTr|iangular|
|00002f90| 4d 61 74 72 69 78 3a 3a | 6f 70 65 72 61 74 6f 72 |Matrix::|operator|
|00002fa0| 28 29 28 69 6e 74 20 6d | 2c 20 69 6e 74 20 6e 29 |()(int m|, int n)|
|00002fb0| 20 63 6f 6e 73 74 0a 58 | 7b 0a 58 20 20 20 69 66 | const.X|{.X   if|
|00002fc0| 20 28 6e 3c 3d 30 20 7c | 7c 20 6d 3c 6e 20 7c 7c | (n<=0 ||| m<n |||
|00002fd0| 20 6d 3e 6e 72 6f 77 73 | 29 20 4d 61 74 72 69 78 | m>nrows|) Matrix|
|00002fe0| 45 72 72 6f 72 28 22 49 | 6e 64 65 78 20 6f 75 74 |Error("I|ndex out|
|00002ff0| 20 6f 66 20 72 61 6e 67 | 65 22 29 3b 0a 58 20 20 | of rang|e");.X  |
|00003000| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 74 72 | return |store[tr|
|00003010| 69 73 74 6f 72 65 28 6d | 2d 31 29 2b 6e 2d 31 5d |istore(m|-1)+n-1]|
|00003020| 3b 0a 58 7d 0a 58 0a 58 | 72 65 61 6c 20 44 69 61 |;.X}.X.X|real Dia|
|00003030| 67 6f 6e 61 6c 4d 61 74 | 72 69 78 3a 3a 6f 70 65 |gonalMat|rix::ope|
|00003040| 72 61 74 6f 72 28 29 28 | 69 6e 74 20 6d 2c 20 69 |rator()(|int m, i|
|00003050| 6e 74 20 6e 29 20 63 6f | 6e 73 74 0a 58 7b 0a 58 |nt n) co|nst.X{.X|
|00003060| 20 20 20 69 66 20 28 6e | 3c 3d 30 20 7c 7c 20 6d |   if (n|<=0 || m|
|00003070| 21 3d 6e 20 7c 7c 20 6d | 3e 6e 72 6f 77 73 20 7c |!=n || m|>nrows ||
|00003080| 7c 20 6e 3e 6e 63 6f 6c | 73 29 20 4d 61 74 72 69 || n>ncol|s) Matri|
|00003090| 78 45 72 72 6f 72 28 22 | 49 6e 64 65 78 20 6f 75 |xError("|Index ou|
|000030a0| 74 20 6f 66 20 72 61 6e | 67 65 22 29 3b 0a 58 20 |t of ran|ge");.X |
|000030b0| 20 20 72 65 74 75 72 6e | 20 73 74 6f 72 65 5b 6e |  return| store[n|
|000030c0| 2d 31 5d 3b 0a 58 7d 0a | 58 0a 58 72 65 61 6c 20 |-1];.X}.|X.Xreal |
|000030d0| 44 69 61 67 6f 6e 61 6c | 4d 61 74 72 69 78 3a 3a |Diagonal|Matrix::|
|000030e0| 6f 70 65 72 61 74 6f 72 | 28 29 28 69 6e 74 20 6d |operator|()(int m|
|000030f0| 29 20 63 6f 6e 73 74 0a | 58 7b 0a 58 20 20 20 69 |) const.|X{.X   i|
|00003100| 66 20 28 6d 3c 3d 30 20 | 7c 7c 20 6d 3e 6e 72 6f |f (m<=0 ||| m>nro|
|00003110| 77 73 29 20 4d 61 74 72 | 69 78 45 72 72 6f 72 28 |ws) Matr|ixError(|
|00003120| 22 49 6e 64 65 78 20 6f | 75 74 20 6f 66 20 72 61 |"Index o|ut of ra|
|00003130| 6e 67 65 22 29 3b 0a 58 | 20 20 20 72 65 74 75 72 |nge");.X|   retur|
|00003140| 6e 20 73 74 6f 72 65 5b | 6d 2d 31 5d 3b 0a 58 7d |n store[|m-1];.X}|
|00003150| 0a 58 0a 58 72 65 61 6c | 20 43 6f 6c 75 6d 6e 56 |.X.Xreal| ColumnV|
|00003160| 65 63 74 6f 72 3a 3a 6f | 70 65 72 61 74 6f 72 28 |ector::o|perator(|
|00003170| 29 28 69 6e 74 20 6d 29 | 20 63 6f 6e 73 74 0a 58 |)(int m)| const.X|
|00003180| 7b 0a 58 20 20 20 69 66 | 20 28 6d 3c 3d 30 20 7c |{.X   if| (m<=0 ||
|00003190| 7c 20 6d 3e 20 6e 72 6f | 77 73 29 20 4d 61 74 72 || m> nro|ws) Matr|
|000031a0| 69 78 45 72 72 6f 72 28 | 22 49 6e 64 65 78 20 6f |ixError(|"Index o|
|000031b0| 75 74 20 6f 66 20 72 61 | 6e 67 65 22 29 3b 0a 58 |ut of ra|nge");.X|
|000031c0| 20 20 20 72 65 74 75 72 | 6e 20 73 74 6f 72 65 5b |   retur|n store[|
|000031d0| 6d 2d 31 5d 3b 0a 58 7d | 0a 58 0a 58 72 65 61 6c |m-1];.X}|.X.Xreal|
|000031e0| 20 52 6f 77 56 65 63 74 | 6f 72 3a 3a 6f 70 65 72 | RowVect|or::oper|
|000031f0| 61 74 6f 72 28 29 28 69 | 6e 74 20 6e 29 20 63 6f |ator()(i|nt n) co|
|00003200| 6e 73 74 0a 58 7b 0a 58 | 20 20 20 69 66 20 28 6e |nst.X{.X|   if (n|
|00003210| 3c 3d 30 20 7c 7c 20 6e | 3e 20 6e 63 6f 6c 73 29 |<=0 || n|> ncols)|
|00003220| 20 4d 61 74 72 69 78 45 | 72 72 6f 72 28 22 49 6e | MatrixE|rror("In|
|00003230| 64 65 78 20 6f 75 74 20 | 6f 66 20 72 61 6e 67 65 |dex out |of range|
|00003240| 22 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 73 |");.X   |return s|
|00003250| 74 6f 72 65 5b 6e 2d 31 | 5d 3b 0a 58 7d 0a 58 0a |tore[n-1|];.X}.X.|
|00003260| 58 23 65 6e 64 69 66 0a | 58 0a 58 42 61 73 65 4d |X#endif.|X.XBaseM|
|00003270| 61 74 72 69 78 3a 3a 6f | 70 65 72 61 74 6f 72 20 |atrix::o|perator |
|00003280| 72 65 61 6c 28 29 0a 58 | 7b 0a 58 20 20 20 52 45 |real().X|{.X   RE|
|00003290| 50 4f 52 54 0a 58 20 20 | 20 47 65 6e 65 72 61 6c |PORT.X  | General|
|000032a0| 4d 61 74 72 69 78 2a 20 | 67 6d 20 3d 20 45 76 61 |Matrix* |gm = Eva|
|000032b0| 6c 75 61 74 65 28 29 3b | 0a 58 20 20 20 69 66 20 |luate();|.X   if |
|000032c0| 28 67 6d 2d 3e 6e 72 6f | 77 73 21 3d 31 20 7c 7c |(gm->nro|ws!=1 |||
|000032d0| 20 67 6d 2d 3e 6e 63 6f | 6c 73 21 3d 31 29 0a 58 | gm->nco|ls!=1).X|
|000032e0| 20 20 20 20 20 20 4d 61 | 74 72 69 78 45 72 72 6f |      Ma|trixErro|
|000032f0| 72 28 22 41 74 74 65 6d | 70 74 20 74 6f 20 63 6f |r("Attem|pt to co|
|00003300| 6e 76 65 72 74 20 6e 6f | 6e 20 31 78 31 20 6d 61 |nvert no|n 1x1 ma|
|00003310| 74 72 69 78 20 74 6f 20 | 73 63 61 6c 61 72 22 29 |trix to |scalar")|
|00003320| 3b 0a 58 20 20 20 72 65 | 61 6c 20 78 20 3d 20 2a |;.X   re|al x = *|
|00003330| 28 67 6d 2d 3e 73 74 6f | 72 65 29 3b 20 67 6d 2d |(gm->sto|re); gm-|
|00003340| 3e 74 44 65 6c 65 74 65 | 28 29 3b 20 72 65 74 75 |>tDelete|(); retu|
|00003350| 72 6e 20 78 3b 0a 58 7d | 0a 58 0a 58 41 64 64 65 |rn x;.X}|.X.XAdde|
|00003360| 64 4d 61 74 72 69 78 20 | 42 61 73 65 4d 61 74 72 |dMatrix |BaseMatr|
|00003370| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 2b 28 42 61 |ix::oper|ator+(Ba|
|00003380| 73 65 4d 61 74 72 69 78 | 26 20 62 6d 29 0a 58 7b |seMatrix|& bm).X{|
|00003390| 20 52 45 50 4f 52 54 20 | 72 65 74 75 72 6e 20 41 | REPORT |return A|
|000033a0| 64 64 65 64 4d 61 74 72 | 69 78 28 74 68 69 73 2c |ddedMatr|ix(this,|
|000033b0| 20 26 62 6d 29 3b 20 7d | 0a 58 0a 58 4d 75 6c 74 | &bm); }|.X.XMult|
|000033c0| 69 70 6c 69 65 64 4d 61 | 74 72 69 78 20 42 61 73 |ipliedMa|trix Bas|
|000033d0| 65 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |eMatrix:|:operato|
|000033e0| 72 2a 28 42 61 73 65 4d | 61 74 72 69 78 26 20 62 |r*(BaseM|atrix& b|
|000033f0| 6d 29 0a 58 7b 20 52 45 | 50 4f 52 54 20 72 65 74 |m).X{ RE|PORT ret|
|00003400| 75 72 6e 20 4d 75 6c 74 | 69 70 6c 69 65 64 4d 61 |urn Mult|ipliedMa|
|00003410| 74 72 69 78 28 74 68 69 | 73 2c 20 26 62 6d 29 3b |trix(thi|s, &bm);|
|00003420| 20 7d 0a 58 0a 58 53 6f | 6c 76 65 64 4d 61 74 72 | }.X.XSo|lvedMatr|
|00003430| 69 78 20 49 6e 76 65 72 | 74 65 64 4d 61 74 72 69 |ix Inver|tedMatri|
|00003440| 78 3a 3a 6f 70 65 72 61 | 74 6f 72 2a 28 42 61 73 |x::opera|tor*(Bas|
|00003450| 65 4d 61 74 72 69 78 26 | 20 62 6d 78 29 0a 58 7b |eMatrix&| bmx).X{|
|00003460| 20 52 45 50 4f 52 54 20 | 72 65 74 75 72 6e 20 53 | REPORT |return S|
|00003470| 6f 6c 76 65 64 4d 61 74 | 72 69 78 28 62 6d 2c 20 |olvedMat|rix(bm, |
|00003480| 26 62 6d 78 29 3b 20 7d | 0a 58 0a 58 53 75 62 74 |&bmx); }|.X.XSubt|
|00003490| 72 61 63 74 65 64 4d 61 | 74 72 69 78 20 42 61 73 |ractedMa|trix Bas|
|000034a0| 65 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |eMatrix:|:operato|
|000034b0| 72 2d 28 42 61 73 65 4d | 61 74 72 69 78 26 20 62 |r-(BaseM|atrix& b|
|000034c0| 6d 29 0a 58 7b 20 52 45 | 50 4f 52 54 20 72 65 74 |m).X{ RE|PORT ret|
|000034d0| 75 72 6e 20 53 75 62 74 | 72 61 63 74 65 64 4d 61 |urn Subt|ractedMa|
|000034e0| 74 72 69 78 28 74 68 69 | 73 2c 20 26 62 6d 29 3b |trix(thi|s, &bm);|
|000034f0| 20 7d 0a 58 0a 58 53 68 | 69 66 74 65 64 4d 61 74 | }.X.XSh|iftedMat|
|00003500| 72 69 78 20 42 61 73 65 | 4d 61 74 72 69 78 3a 3a |rix Base|Matrix::|
|00003510| 6f 70 65 72 61 74 6f 72 | 2b 28 72 65 61 6c 20 66 |operator|+(real f|
|00003520| 29 0a 58 7b 20 52 45 50 | 4f 52 54 20 72 65 74 75 |).X{ REP|ORT retu|
|00003530| 72 6e 20 53 68 69 66 74 | 65 64 4d 61 74 72 69 78 |rn Shift|edMatrix|
|00003540| 28 74 68 69 73 2c 20 66 | 29 3b 20 7d 0a 58 0a 58 |(this, f|); }.X.X|
|00003550| 53 63 61 6c 65 64 4d 61 | 74 72 69 78 20 42 61 73 |ScaledMa|trix Bas|
|00003560| 65 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |eMatrix:|:operato|
|00003570| 72 2a 28 72 65 61 6c 20 | 66 29 0a 58 7b 20 52 45 |r*(real |f).X{ RE|
|00003580| 50 4f 52 54 20 72 65 74 | 75 72 6e 20 53 63 61 6c |PORT ret|urn Scal|
|00003590| 65 64 4d 61 74 72 69 78 | 28 74 68 69 73 2c 20 66 |edMatrix|(this, f|
|000035a0| 29 3b 20 7d 0a 58 0a 58 | 53 63 61 6c 65 64 4d 61 |); }.X.X|ScaledMa|
|000035b0| 74 72 69 78 20 42 61 73 | 65 4d 61 74 72 69 78 3a |trix Bas|eMatrix:|
|000035c0| 3a 6f 70 65 72 61 74 6f | 72 2f 28 72 65 61 6c 20 |:operato|r/(real |
|000035d0| 66 29 0a 58 7b 20 52 45 | 50 4f 52 54 20 72 65 74 |f).X{ RE|PORT ret|
|000035e0| 75 72 6e 20 53 63 61 6c | 65 64 4d 61 74 72 69 78 |urn Scal|edMatrix|
|000035f0| 28 74 68 69 73 2c 20 31 | 2e 30 2f 66 29 3b 20 7d |(this, 1|.0/f); }|
|00003600| 0a 58 0a 58 53 68 69 66 | 74 65 64 4d 61 74 72 69 |.X.XShif|tedMatri|
|00003610| 78 20 42 61 73 65 4d 61 | 74 72 69 78 3a 3a 6f 70 |x BaseMa|trix::op|
|00003620| 65 72 61 74 6f 72 2d 28 | 72 65 61 6c 20 66 29 0a |erator-(|real f).|
|00003630| 58 7b 20 52 45 50 4f 52 | 54 20 72 65 74 75 72 6e |X{ REPOR|T return|
|00003640| 20 53 68 69 66 74 65 64 | 4d 61 74 72 69 78 28 74 | Shifted|Matrix(t|
|00003650| 68 69 73 2c 20 2d 66 29 | 3b 20 7d 0a 58 0a 58 54 |his, -f)|; }.X.XT|
|00003660| 72 61 6e 73 70 6f 73 65 | 64 4d 61 74 72 69 78 20 |ranspose|dMatrix |
|00003670| 42 61 73 65 4d 61 74 72 | 69 78 3a 3a 74 28 29 20 |BaseMatr|ix::t() |
|00003680| 7b 20 52 45 50 4f 52 54 | 20 72 65 74 75 72 6e 20 |{ REPORT| return |
|00003690| 54 72 61 6e 73 70 6f 73 | 65 64 4d 61 74 72 69 78 |Transpos|edMatrix|
|000036a0| 28 74 68 69 73 29 3b 20 | 7d 0a 58 0a 58 4e 65 67 |(this); |}.X.XNeg|
|000036b0| 61 74 65 64 4d 61 74 72 | 69 78 20 42 61 73 65 4d |atedMatr|ix BaseM|
|000036c0| 61 74 72 69 78 3a 3a 6f | 70 65 72 61 74 6f 72 2d |atrix::o|perator-|
|000036d0| 28 29 20 7b 20 52 45 50 | 4f 52 54 20 72 65 74 75 |() { REP|ORT retu|
|000036e0| 72 6e 20 4e 65 67 61 74 | 65 64 4d 61 74 72 69 78 |rn Negat|edMatrix|
|000036f0| 28 74 68 69 73 29 3b 20 | 7d 0a 58 0a 58 49 6e 76 |(this); |}.X.XInv|
|00003700| 65 72 74 65 64 4d 61 74 | 72 69 78 20 42 61 73 65 |ertedMat|rix Base|
|00003710| 4d 61 74 72 69 78 3a 3a | 69 28 29 20 7b 20 52 45 |Matrix::|i() { RE|
|00003720| 50 4f 52 54 20 72 65 74 | 75 72 6e 20 49 6e 76 65 |PORT ret|urn Inve|
|00003730| 72 74 65 64 4d 61 74 72 | 69 78 28 74 68 69 73 29 |rtedMatr|ix(this)|
|00003740| 3b 20 7d 0a 58 0a 58 43 | 6f 6e 73 74 4d 61 74 72 |; }.X.XC|onstMatr|
|00003750| 69 78 20 47 65 6e 65 72 | 61 6c 4d 61 74 72 69 78 |ix Gener|alMatrix|
|00003760| 3a 3a 63 28 29 20 63 6f | 6e 73 74 0a 58 7b 0a 58 |::c() co|nst.X{.X|
|00003770| 20 20 20 69 66 20 28 74 | 61 67 20 21 3d 20 2d 31 |   if (t|ag != -1|
|00003780| 29 20 4d 61 74 72 69 78 | 45 72 72 6f 72 28 22 2e |) Matrix|Error(".|
|00003790| 63 28 29 20 61 70 70 6c | 69 65 64 20 74 6f 20 74 |c() appl|ied to t|
|000037a0| 65 6d 70 6f 72 61 72 79 | 20 6d 61 74 72 69 78 22 |emporary| matrix"|
|000037b0| 29 3b 0a 58 20 20 20 52 | 45 50 4f 52 54 20 72 65 |);.X   R|EPORT re|
|000037c0| 74 75 72 6e 20 43 6f 6e | 73 74 4d 61 74 72 69 78 |turn Con|stMatrix|
|000037d0| 28 74 68 69 73 29 3b 0a | 58 7d 0a 58 0a 58 52 6f |(this);.|X}.X.XRo|
|000037e0| 77 65 64 4d 61 74 72 69 | 78 20 42 61 73 65 4d 61 |wedMatri|x BaseMa|
|000037f0| 74 72 69 78 3a 3a 43 6f | 70 79 54 6f 52 6f 77 28 |trix::Co|pyToRow(|
|00003800| 29 20 7b 20 52 45 50 4f | 52 54 20 72 65 74 75 72 |) { REPO|RT retur|
|00003810| 6e 20 52 6f 77 65 64 4d | 61 74 72 69 78 28 74 68 |n RowedM|atrix(th|
|00003820| 69 73 29 3b 20 7d 0a 58 | 0a 58 43 6f 6c 65 64 4d |is); }.X|.XColedM|
|00003830| 61 74 72 69 78 20 42 61 | 73 65 4d 61 74 72 69 78 |atrix Ba|seMatrix|
|00003840| 3a 3a 43 6f 70 79 54 6f | 43 6f 6c 75 6d 6e 28 29 |::CopyTo|Column()|
|00003850| 20 7b 20 52 45 50 4f 52 | 54 20 72 65 74 75 72 6e | { REPOR|T return|
|00003860| 20 43 6f 6c 65 64 4d 61 | 74 72 69 78 28 74 68 69 | ColedMa|trix(thi|
|00003870| 73 29 3b 20 7d 0a 58 0a | 58 44 69 61 67 65 64 4d |s); }.X.|XDiagedM|
|00003880| 61 74 72 69 78 20 42 61 | 73 65 4d 61 74 72 69 78 |atrix Ba|seMatrix|
|00003890| 3a 3a 43 6f 70 79 54 6f | 44 69 61 67 6f 6e 61 6c |::CopyTo|Diagonal|
|000038a0| 28 29 20 7b 20 52 45 50 | 4f 52 54 20 72 65 74 75 |() { REP|ORT retu|
|000038b0| 72 6e 20 44 69 61 67 65 | 64 4d 61 74 72 69 78 28 |rn Diage|dMatrix(|
|000038c0| 74 68 69 73 29 3b 20 7d | 0a 58 0a 58 4d 61 74 65 |this); }|.X.XMate|
|000038d0| 64 4d 61 74 72 69 78 20 | 42 61 73 65 4d 61 74 72 |dMatrix |BaseMatr|
|000038e0| 69 78 3a 3a 43 6f 70 79 | 54 6f 4d 61 74 72 69 78 |ix::Copy|ToMatrix|
|000038f0| 28 69 6e 74 20 6e 72 78 | 2c 20 69 6e 74 20 6e 63 |(int nrx|, int nc|
|00003900| 78 29 0a 58 7b 20 52 45 | 50 4f 52 54 20 72 65 74 |x).X{ RE|PORT ret|
|00003910| 75 72 6e 20 4d 61 74 65 | 64 4d 61 74 72 69 78 28 |urn Mate|dMatrix(|
|00003920| 74 68 69 73 2c 6e 72 78 | 2c 6e 63 78 29 3b 20 7d |this,nrx|,ncx); }|
|00003930| 0a 58 0a 58 76 6f 69 64 | 20 47 65 6e 65 72 61 6c |.X.Xvoid| General|
|00003940| 4d 61 74 72 69 78 3a 3a | 6f 70 65 72 61 74 6f 72 |Matrix::|operator|
|00003950| 3d 28 72 65 61 6c 20 66 | 29 0a 58 7b 20 52 45 50 |=(real f|).X{ REP|
|00003960| 4f 52 54 20 69 6e 74 20 | 69 3d 73 74 6f 72 61 67 |ORT int |i=storag|
|00003970| 65 3b 20 72 65 61 6c 2a | 20 73 3d 73 74 6f 72 65 |e; real*| s=store|
|00003980| 3b 20 77 68 69 6c 65 20 | 28 69 2d 2d 29 20 7b 20 |; while |(i--) { |
|00003990| 2a 73 2b 2b 20 3d 20 66 | 3b 20 7d 20 7d 0a 58 0a |*s++ = f|; } }.X.|
|000039a0| 58 76 6f 69 64 20 4d 61 | 74 72 69 78 3a 3a 6f 70 |Xvoid Ma|trix::op|
|000039b0| 65 72 61 74 6f 72 3d 28 | 42 61 73 65 4d 61 74 72 |erator=(|BaseMatr|
|000039c0| 69 78 26 20 58 29 0a 58 | 7b 20 52 45 50 4f 52 54 |ix& X).X|{ REPORT|
|000039d0| 20 43 68 65 63 6b 43 6f | 6e 76 65 72 73 69 6f 6e | CheckCo|nversion|
|000039e0| 28 58 29 3b 20 45 71 28 | 58 2c 4d 61 74 72 69 78 |(X); Eq(|X,Matrix|
|000039f0| 54 79 70 65 3a 3a 52 65 | 63 74 29 3b 20 7d 20 0a |Type::Re|ct); } .|
|00003a00| 58 0a 58 76 6f 69 64 20 | 52 6f 77 56 65 63 74 6f |X.Xvoid |RowVecto|
|00003a10| 72 3a 3a 6f 70 65 72 61 | 74 6f 72 3d 28 42 61 73 |r::opera|tor=(Bas|
|00003a20| 65 4d 61 74 72 69 78 26 | 20 58 29 0a 58 7b 0a 58 |eMatrix&| X).X{.X|
|00003a30| 20 20 20 52 45 50 4f 52 | 54 20 43 68 65 63 6b 43 |   REPOR|T CheckC|
|00003a40| 6f 6e 76 65 72 73 69 6f | 6e 28 58 29 3b 20 45 71 |onversio|n(X); Eq|
|00003a50| 28 58 2c 4d 61 74 72 69 | 78 54 79 70 65 3a 3a 52 |(X,Matri|xType::R|
|00003a60| 6f 77 56 29 3b 0a 58 20 | 20 20 69 66 20 28 6e 72 |owV);.X |  if (nr|
|00003a70| 6f 77 73 21 3d 31 29 20 | 4d 61 74 72 69 78 45 72 |ows!=1) |MatrixEr|
|00003a80| 72 6f 72 28 22 49 6c 6c | 65 67 61 6c 20 63 6f 6e |ror("Ill|egal con|
|00003a90| 76 65 72 73 69 6f 6e 20 | 74 6f 20 72 6f 77 20 76 |version |to row v|
|00003aa0| 65 63 74 6f 72 22 29 3b | 0a 58 7d 0a 58 0a 58 76 |ector");|.X}.X.Xv|
|00003ab0| 6f 69 64 20 43 6f 6c 75 | 6d 6e 56 65 63 74 6f 72 |oid Colu|mnVector|
|00003ac0| 3a 3a 6f 70 65 72 61 74 | 6f 72 3d 28 42 61 73 65 |::operat|or=(Base|
|00003ad0| 4d 61 74 72 69 78 26 20 | 58 29 0a 58 7b 0a 58 20 |Matrix& |X).X{.X |
|00003ae0| 20 20 52 45 50 4f 52 54 | 20 43 68 65 63 6b 43 6f |  REPORT| CheckCo|
|00003af0| 6e 76 65 72 73 69 6f 6e | 28 58 29 3b 20 45 71 28 |nversion|(X); Eq(|
|00003b00| 58 2c 4d 61 74 72 69 78 | 54 79 70 65 3a 3a 43 6f |X,Matrix|Type::Co|
|00003b10| 6c 56 29 3b 0a 58 20 20 | 20 69 66 20 28 6e 63 6f |lV);.X  | if (nco|
|00003b20| 6c 73 21 3d 31 29 20 4d | 61 74 72 69 78 45 72 72 |ls!=1) M|atrixErr|
|00003b30| 6f 72 28 22 49 6c 6c 65 | 67 61 6c 20 63 6f 6e 76 |or("Ille|gal conv|
|00003b40| 65 72 73 69 6f 6e 20 74 | 6f 20 63 6f 6c 75 6d 6e |ersion t|o column|
|00003b50| 20 76 65 63 74 6f 72 22 | 29 3b 0a 58 7d 0a 58 0a | vector"|);.X}.X.|
|00003b60| 58 76 6f 69 64 20 53 79 | 6d 6d 65 74 72 69 63 4d |Xvoid Sy|mmetricM|
|00003b70| 61 74 72 69 78 3a 3a 6f | 70 65 72 61 74 6f 72 3d |atrix::o|perator=|
|00003b80| 28 42 61 73 65 4d 61 74 | 72 69 78 26 20 58 29 0a |(BaseMat|rix& X).|
|00003b90| 58 7b 20 52 45 50 4f 52 | 54 20 43 68 65 63 6b 43 |X{ REPOR|T CheckC|
|00003ba0| 6f 6e 76 65 72 73 69 6f | 6e 28 58 29 3b 20 45 71 |onversio|n(X); Eq|
|00003bb0| 28 58 2c 4d 61 74 72 69 | 78 54 79 70 65 3a 3a 53 |(X,Matri|xType::S|
|00003bc0| 79 6d 29 3b 20 7d 0a 58 | 20 0a 58 76 6f 69 64 20 |ym); }.X| .Xvoid |
|00003bd0| 55 70 70 65 72 54 72 69 | 61 6e 67 75 6c 61 72 4d |UpperTri|angularM|
|00003be0| 61 74 72 69 78 3a 3a 6f | 70 65 72 61 74 6f 72 3d |atrix::o|perator=|
|00003bf0| 28 42 61 73 65 4d 61 74 | 72 69 78 26 20 58 29 0a |(BaseMat|rix& X).|
|00003c00| 58 7b 20 52 45 50 4f 52 | 54 20 43 68 65 63 6b 43 |X{ REPOR|T CheckC|
|00003c10| 6f 6e 76 65 72 73 69 6f | 6e 28 58 29 3b 20 45 71 |onversio|n(X); Eq|
|00003c20| 28 58 2c 4d 61 74 72 69 | 78 54 79 70 65 3a 3a 55 |(X,Matri|xType::U|
|00003c30| 54 29 3b 20 7d 0a 58 0a | 58 76 6f 69 64 20 4c 6f |T); }.X.|Xvoid Lo|
|00003c40| 77 65 72 54 72 69 61 6e | 67 75 6c 61 72 4d 61 74 |werTrian|gularMat|
|00003c50| 72 69 78 3a 3a 6f 70 65 | 72 61 74 6f 72 3d 28 42 |rix::ope|rator=(B|
|00003c60| 61 73 65 4d 61 74 72 69 | 78 26 20 58 29 0a 58 7b |aseMatri|x& X).X{|
|00003c70| 20 52 45 50 4f 52 54 20 | 43 68 65 63 6b 43 6f 6e | REPORT |CheckCon|
|00003c80| 76 65 72 73 69 6f 6e 28 | 58 29 3b 20 45 71 28 58 |version(|X); Eq(X|
|00003c90| 2c 4d 61 74 72 69 78 54 | 79 70 65 3a 3a 4c 54 29 |,MatrixT|ype::LT)|
|00003ca0| 3b 20 7d 0a 58 0a 58 76 | 6f 69 64 20 44 69 61 67 |; }.X.Xv|oid Diag|
|00003cb0| 6f 6e 61 6c 4d 61 74 72 | 69 78 3a 3a 6f 70 65 72 |onalMatr|ix::oper|
|00003cc0| 61 74 6f 72 3d 28 42 61 | 73 65 4d 61 74 72 69 78 |ator=(Ba|seMatrix|
|00003cd0| 26 20 58 29 0a 58 7b 20 | 52 45 50 4f 52 54 20 43 |& X).X{ |REPORT C|
|00003ce0| 68 65 63 6b 43 6f 6e 76 | 65 72 73 69 6f 6e 28 58 |heckConv|ersion(X|
|00003cf0| 29 3b 20 45 71 28 58 2c | 4d 61 74 72 69 78 54 79 |); Eq(X,|MatrixTy|
|00003d00| 70 65 3a 3a 44 69 61 67 | 29 3b 20 7d 0a 58 0a 58 |pe::Diag|); }.X.X|
|00003d10| 76 6f 69 64 20 47 65 6e | 65 72 61 6c 4d 61 74 72 |void Gen|eralMatr|
|00003d20| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 3c 3c 28 63 |ix::oper|ator<<(c|
|00003d30| 6f 6e 73 74 20 72 65 61 | 6c 2a 20 72 29 0a 58 7b |onst rea|l* r).X{|
|00003d40| 0a 58 20 20 20 52 45 50 | 4f 52 54 0a 58 20 20 20 |.X   REP|ORT.X   |
|00003d50| 69 6e 74 20 69 20 3d 20 | 73 74 6f 72 61 67 65 3b |int i = |storage;|
|00003d60| 20 72 65 61 6c 2a 20 73 | 3d 73 74 6f 72 65 3b 0a | real* s|=store;.|
|00003d70| 58 20 20 20 77 68 69 6c | 65 28 69 2d 2d 29 20 2a |X   whil|e(i--) *|
|00003d80| 73 2b 2b 20 3d 20 2a 72 | 2b 2b 3b 0a 58 7d 0a 58 |s++ = *r|++;.X}.X|
|00003d90| 0a 58 0a 58 2f 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X.X/***|********|
|00003da0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 20 65 |********|****** e|
|00003db0| 6c 65 6d 65 6e 74 20 61 | 63 63 65 73 73 20 2a 2a |lement a|ccess **|
|00003dc0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003dd0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2f |********|*******/|
|00003de0| 0a 58 0a 58 72 65 61 6c | 26 20 4d 61 74 72 69 78 |.X.Xreal|& Matrix|
|00003df0| 3a 3a 65 6c 65 6d 65 6e | 74 28 69 6e 74 20 6d 2c |::elemen|t(int m,|
|00003e00| 20 69 6e 74 20 6e 29 0a | 58 7b 0a 58 20 20 20 69 | int n).|X{.X   i|
|00003e10| 66 20 28 6d 3c 30 20 7c | 7c 20 6d 3e 3d 20 6e 72 |f (m<0 ||| m>= nr|
|00003e20| 6f 77 73 20 7c 7c 20 6e | 3c 30 20 7c 7c 20 6e 3e |ows || n|<0 || n>|
|00003e30| 3d 20 6e 63 6f 6c 73 29 | 20 4d 61 74 72 69 78 45 |= ncols)| MatrixE|
|00003e40| 72 72 6f 72 28 22 49 6e | 64 65 78 20 6f 75 74 20 |rror("In|dex out |
|00003e50| 6f 66 20 72 61 6e 67 65 | 22 29 3b 0a 58 20 20 20 |of range|");.X   |
|00003e60| 72 65 74 75 72 6e 20 73 | 74 6f 72 65 5b 6d 2a 6e |return s|tore[m*n|
|00003e70| 63 6f 6c 73 2b 6e 5d 3b | 0a 58 7d 0a 58 0a 58 72 |cols+n];|.X}.X.Xr|
|00003e80| 65 61 6c 26 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |eal& Sym|metricMa|
|00003e90| 74 72 69 78 3a 3a 65 6c | 65 6d 65 6e 74 28 69 6e |trix::el|ement(in|
|00003ea0| 74 20 6d 2c 20 69 6e 74 | 20 6e 29 0a 58 7b 0a 58 |t m, int| n).X{.X|
|00003eb0| 20 20 20 69 66 20 28 6d | 3c 30 20 7c 7c 20 6e 3c |   if (m|<0 || n<|
|00003ec0| 30 20 7c 7c 20 6d 20 3e | 3d 20 6e 72 6f 77 73 20 |0 || m >|= nrows |
|00003ed0| 7c 7c 20 6e 3e 3d 6e 63 | 6f 6c 73 29 20 4d 61 74 ||| n>=nc|ols) Mat|
|00003ee0| 72 69 78 45 72 72 6f 72 | 28 22 49 6e 64 65 78 20 |rixError|("Index |
|00003ef0| 6f 75 74 20 6f 66 20 72 | 61 6e 67 65 22 29 3b 0a |out of r|ange");.|
|00003f00| 58 20 20 20 69 66 20 28 | 6d 3e 3d 6e 29 20 72 65 |X   if (|m>=n) re|
|00003f10| 74 75 72 6e 20 73 74 6f | 72 65 5b 74 72 69 73 74 |turn sto|re[trist|
|00003f20| 6f 72 65 28 6d 29 2b 6e | 5d 3b 0a 58 20 20 20 65 |ore(m)+n|];.X   e|
|00003f30| 6c 73 65 20 72 65 74 75 | 72 6e 20 73 74 6f 72 65 |lse retu|rn store|
|00003f40| 5b 74 72 69 73 74 6f 72 | 65 28 6e 29 2b 6d 5d 3b |[tristor|e(n)+m];|
|00003f50| 0a 58 7d 0a 58 0a 58 72 | 65 61 6c 26 20 55 70 70 |.X}.X.Xr|eal& Upp|
|00003f60| 65 72 54 72 69 61 6e 67 | 75 6c 61 72 4d 61 74 72 |erTriang|ularMatr|
|00003f70| 69 78 3a 3a 65 6c 65 6d | 65 6e 74 28 69 6e 74 20 |ix::elem|ent(int |
|00003f80| 6d 2c 20 69 6e 74 20 6e | 29 0a 58 7b 0a 58 20 20 |m, int n|).X{.X  |
|00003f90| 20 69 66 20 28 6d 3c 30 | 20 7c 7c 20 6e 3c 6d 20 | if (m<0| || n<m |
|00003fa0| 7c 7c 20 6e 3e 3d 6e 63 | 6f 6c 73 29 20 4d 61 74 ||| n>=nc|ols) Mat|
|00003fb0| 72 69 78 45 72 72 6f 72 | 28 22 49 6e 64 65 78 20 |rixError|("Index |
|00003fc0| 6f 75 74 20 6f 66 20 72 | 61 6e 67 65 22 29 3b 0a |out of r|ange");.|
|00003fd0| 58 20 20 20 72 65 74 75 | 72 6e 20 73 74 6f 72 65 |X   retu|rn store|
|00003fe0| 5b 6d 2a 6e 63 6f 6c 73 | 2b 6e 2d 74 72 69 73 74 |[m*ncols|+n-trist|
|00003ff0| 6f 72 65 28 6d 29 5d 3b | 0a 58 7d 0a 58 0a 58 72 |ore(m)];|.X}.X.Xr|
|00004000| 65 61 6c 26 20 4c 6f 77 | 65 72 54 72 69 61 6e 67 |eal& Low|erTriang|
|00004010| 75 6c 61 72 4d 61 74 72 | 69 78 3a 3a 65 6c 65 6d |ularMatr|ix::elem|
|00004020| 65 6e 74 28 69 6e 74 20 | 6d 2c 20 69 6e 74 20 6e |ent(int |m, int n|
|00004030| 29 0a 58 7b 0a 58 20 20 | 20 69 66 20 28 6e 3c 30 |).X{.X  | if (n<0|
|00004040| 20 7c 7c 20 6d 3c 6e 20 | 7c 7c 20 6d 3e 3d 6e 72 | || m<n ||| m>=nr|
|00004050| 6f 77 73 29 20 4d 61 74 | 72 69 78 45 72 72 6f 72 |ows) Mat|rixError|
|00004060| 28 22 49 6e 64 65 78 20 | 6f 75 74 20 6f 66 20 72 |("Index |out of r|
|00004070| 61 6e 67 65 22 29 3b 0a | 58 20 20 20 72 65 74 75 |ange");.|X   retu|
|00004080| 72 6e 20 73 74 6f 72 65 | 5b 74 72 69 73 74 6f 72 |rn store|[tristor|
|00004090| 65 28 6d 29 2b 6e 5d 3b | 0a 58 7d 0a 58 0a 58 72 |e(m)+n];|.X}.X.Xr|
|000040a0| 65 61 6c 26 20 44 69 61 | 67 6f 6e 61 6c 4d 61 74 |eal& Dia|gonalMat|
|000040b0| 72 69 78 3a 3a 65 6c 65 | 6d 65 6e 74 28 69 6e 74 |rix::ele|ment(int|
|000040c0| 20 6d 2c 20 69 6e 74 20 | 6e 29 0a 58 7b 0a 58 20 | m, int |n).X{.X |
|000040d0| 20 20 69 66 20 28 6e 3c | 30 20 7c 7c 20 6d 21 3d |  if (n<|0 || m!=|
|000040e0| 6e 20 7c 7c 20 6d 3e 3d | 6e 72 6f 77 73 20 7c 7c |n || m>=|nrows |||
|000040f0| 20 6e 3e 3d 6e 63 6f 6c | 73 29 20 4d 61 74 72 69 | n>=ncol|s) Matri|
|00004100| 78 45 72 72 6f 72 28 22 | 49 6e 64 65 78 20 6f 75 |xError("|Index ou|
|00004110| 74 20 6f 66 20 72 61 6e | 67 65 22 29 3b 0a 58 20 |t of ran|ge");.X |
|00004120| 20 20 72 65 74 75 72 6e | 20 73 74 6f 72 65 5b 6e |  return| store[n|
|00004130| 5d 3b 0a 58 7d 0a 58 0a | 58 72 65 61 6c 26 20 44 |];.X}.X.|Xreal& D|
|00004140| 69 61 67 6f 6e 61 6c 4d | 61 74 72 69 78 3a 3a 65 |iagonalM|atrix::e|
|00004150| 6c 65 6d 65 6e 74 28 69 | 6e 74 20 6d 29 0a 58 7b |lement(i|nt m).X{|
|00004160| 0a 58 20 20 20 69 66 20 | 28 6d 3c 30 20 7c 7c 20 |.X   if |(m<0 || |
|00004170| 6d 3e 3d 6e 72 6f 77 73 | 29 20 4d 61 74 72 69 78 |m>=nrows|) Matrix|
|00004180| 45 72 72 6f 72 28 22 49 | 6e 64 65 78 20 6f 75 74 |Error("I|ndex out|
|00004190| 20 6f 66 20 72 61 6e 67 | 65 22 29 3b 0a 58 20 20 | of rang|e");.X  |
|000041a0| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 6d 5d | return |store[m]|
|000041b0| 3b 0a 58 7d 0a 58 0a 58 | 72 65 61 6c 26 20 43 6f |;.X}.X.X|real& Co|
|000041c0| 6c 75 6d 6e 56 65 63 74 | 6f 72 3a 3a 65 6c 65 6d |lumnVect|or::elem|
|000041d0| 65 6e 74 28 69 6e 74 20 | 6d 29 0a 58 7b 0a 58 20 |ent(int |m).X{.X |
|000041e0| 20 20 69 66 20 28 6d 3c | 30 20 7c 7c 20 6d 3e 3d |  if (m<|0 || m>=|
|000041f0| 20 6e 72 6f 77 73 29 20 | 4d 61 74 72 69 78 45 72 | nrows) |MatrixEr|
|00004200| 72 6f 72 28 22 49 6e 64 | 65 78 20 6f 75 74 20 6f |ror("Ind|ex out o|
|00004210| 66 20 72 61 6e 67 65 22 | 29 3b 0a 58 20 20 20 72 |f range"|);.X   r|
|00004220| 65 74 75 72 6e 20 73 74 | 6f 72 65 5b 6d 5d 3b 0a |eturn st|ore[m];.|
|00004230| 58 7d 0a 58 0a 58 72 65 | 61 6c 26 20 52 6f 77 56 |X}.X.Xre|al& RowV|
|00004240| 65 63 74 6f 72 3a 3a 65 | 6c 65 6d 65 6e 74 28 69 |ector::e|lement(i|
|00004250| 6e 74 20 6e 29 0a 58 7b | 0a 58 20 20 20 69 66 20 |nt n).X{|.X   if |
|00004260| 28 6e 3c 30 20 7c 7c 20 | 6e 3e 3d 20 6e 63 6f 6c |(n<0 || |n>= ncol|
|00004270| 73 29 20 4d 61 74 72 69 | 78 45 72 72 6f 72 28 22 |s) Matri|xError("|
|00004280| 49 6e 64 65 78 20 6f 75 | 74 20 6f 66 20 72 61 6e |Index ou|t of ran|
|00004290| 67 65 22 29 3b 0a 58 20 | 20 20 72 65 74 75 72 6e |ge");.X |  return|
|000042a0| 20 73 74 6f 72 65 5b 6e | 5d 3b 0a 58 7d 0a 58 0a | store[n|];.X}.X.|
|000042b0| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 0a 69 66 20 74 |END_OF_F|ILE.if t|
|000042c0| 65 73 74 20 37 30 38 34 | 20 2d 6e 65 20 60 77 63 |est 7084| -ne `wc|
|000042d0| 20 2d 63 20 3c 27 6e 65 | 77 6d 61 74 36 2e 63 78 | -c <'ne|wmat6.cx|
|000042e0| 78 27 60 3b 20 74 68 65 | 6e 0a 20 20 20 20 65 63 |x'`; the|n.    ec|
|000042f0| 68 6f 20 73 68 61 72 3a | 20 5c 22 27 6e 65 77 6d |ho shar:| \"'newm|
|00004300| 61 74 36 2e 63 78 78 27 | 5c 22 20 75 6e 70 61 63 |at6.cxx'|\" unpac|
|00004310| 6b 65 64 20 77 69 74 68 | 20 77 72 6f 6e 67 20 73 |ked with| wrong s|
|00004320| 69 7a 65 21 0a 66 69 0a | 23 20 65 6e 64 20 6f 66 |ize!.fi.|# end of|
|00004330| 20 27 6e 65 77 6d 61 74 | 36 2e 63 78 78 27 0a 66 | 'newmat|6.cxx'.f|
|00004340| 69 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 6e 65 |i.if tes|t -f 'ne|
|00004350| 77 6d 61 74 37 2e 63 78 | 78 27 20 2d 61 20 22 24 |wmat7.cx|x' -a "$|
|00004360| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|00004370| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen .  e|cho shar|
|00004380| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 63 6c 6f 62 62 |: Will n|ot clobb|
|00004390| 65 72 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |er exist|ing file|
|000043a0| 20 5c 22 27 6e 65 77 6d | 61 74 37 2e 63 78 78 27 | \"'newm|at7.cxx'|
|000043b0| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|000043c0| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|000043d0| 27 6e 65 77 6d 61 74 37 | 2e 63 78 78 27 5c 22 20 |'newmat7|.cxx'\" |
|000043e0| 5c 28 31 31 35 33 35 20 | 63 68 61 72 61 63 74 65 |\(11535 |characte|
|000043f0| 72 73 5c 29 0a 73 65 64 | 20 22 73 2f 5e 58 2f 2f |rs\).sed| "s/^X//|
|00004400| 22 20 3e 27 6e 65 77 6d | 61 74 37 2e 63 78 78 27 |" >'newm|at7.cxx'|
|00004410| 20 3c 3c 27 45 4e 44 5f | 4f 46 5f 46 49 4c 45 27 | <<'END_|OF_FILE'|
|00004420| 0a 58 2f 2f 24 24 20 6e | 65 77 6d 61 74 37 2e 63 |.X//$$ n|ewmat7.c|
|00004430| 78 78 20 20 20 20 20 49 | 6e 76 65 72 74 2c 20 73 |xx     I|nvert, s|
|00004440| 6f 6c 76 65 2c 20 62 69 | 6e 61 72 79 20 6f 70 65 |olve, bi|nary ope|
|00004450| 72 61 74 69 6f 6e 73 0a | 58 0a 58 2f 2f 20 43 6f |rations.|X.X// Co|
|00004460| 70 79 72 69 67 68 74 20 | 28 43 29 20 31 39 39 31 |pyright |(C) 1991|
|00004470| 3a 20 52 20 42 20 44 61 | 76 69 65 73 20 61 6e 64 |: R B Da|vies and|
|00004480| 20 44 53 49 52 0a 58 0a | 58 23 69 6e 63 6c 75 64 | DSIR.X.|X#includ|
|00004490| 65 20 22 69 6e 63 6c 75 | 64 65 2e 68 78 78 22 0a |e "inclu|de.hxx".|
|000044a0| 58 0a 58 23 69 6e 63 6c | 75 64 65 20 22 6e 65 77 |X.X#incl|ude "new|
|000044b0| 6d 61 74 2e 68 78 78 22 | 0a 58 23 69 6e 63 6c 75 |mat.hxx"|.X#inclu|
|000044c0| 64 65 20 22 6e 65 77 6d | 61 74 72 63 2e 68 78 78 |de "newm|atrc.hxx|
|000044d0| 22 0a 58 0a 58 2f 2f 23 | 64 65 66 69 6e 65 20 52 |".X.X//#|define R|
|000044e0| 45 50 4f 52 54 20 7b 20 | 73 74 61 74 69 63 20 45 |EPORT { |static E|
|000044f0| 78 65 43 6f 75 6e 74 65 | 72 20 45 78 65 43 6f 75 |xeCounte|r ExeCou|
|00004500| 6e 74 28 5f 5f 4c 49 4e | 45 5f 5f 2c 37 29 3b 20 |nt(__LIN|E__,7); |
|00004510| 45 78 65 43 6f 75 6e 74 | 2b 2b 3b 20 7d 0a 58 0a |ExeCount|++; }.X.|
|00004520| 58 23 64 65 66 69 6e 65 | 20 52 45 50 4f 52 54 20 |X#define| REPORT |
|00004530| 7b 7d 0a 58 0a 58 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a |{}.X.X.X|/*******|
|00004540| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004550| 2a 2a 2a 2a 2a 2a 20 73 | 6f 6c 76 65 20 72 6f 75 |****** s|olve rou|
|00004560| 74 69 6e 65 73 20 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |tines **|********|
|00004570| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004580| 2a 2a 2a 2a 2f 0a 58 0a | 58 47 65 6e 65 72 61 6c |****/.X.|XGeneral|
|00004590| 4d 61 74 72 69 78 2a 20 | 47 65 6e 65 72 61 6c 4d |Matrix* |GeneralM|
|000045a0| 61 74 72 69 78 3a 3a 4d | 61 6b 65 53 6f 6c 76 65 |atrix::M|akeSolve|
|000045b0| 72 28 29 0a 58 7b 0a 58 | 20 20 20 52 45 50 4f 52 |r().X{.X|   REPOR|
|000045c0| 54 0a 58 20 20 20 47 65 | 6e 65 72 61 6c 4d 61 74 |T.X   Ge|neralMat|
|000045d0| 72 69 78 2a 20 67 6d 20 | 3d 20 6e 65 77 20 43 72 |rix* gm |= new Cr|
|000045e0| 6f 75 74 4d 61 74 72 69 | 78 28 2a 74 68 69 73 29 |outMatri|x(*this)|
|000045f0| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 45 72 72 6f |;.X   Ma|trixErro|
|00004600| 72 4e 6f 53 70 61 63 65 | 28 67 6d 29 3b 20 67 6d |rNoSpace|(gm); gm|
|00004610| 2d 3e 52 65 6c 65 61 73 | 65 41 6e 64 44 65 6c 65 |->Releas|eAndDele|
|00004620| 74 65 28 29 3b 20 72 65 | 74 75 72 6e 20 67 6d 3b |te(); re|turn gm;|
|00004630| 0a 58 7d 0a 58 0a 58 47 | 65 6e 65 72 61 6c 4d 61 |.X}.X.XG|eneralMa|
|00004640| 74 72 69 78 2a 20 4d 61 | 74 72 69 78 3a 3a 4d 61 |trix* Ma|trix::Ma|
|00004650| 6b 65 53 6f 6c 76 65 72 | 28 29 0a 58 7b 0a 58 20 |keSolver|().X{.X |
|00004660| 20 20 52 45 50 4f 52 54 | 0a 58 20 20 20 47 65 6e |  REPORT|.X   Gen|
|00004670| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 67 6d 20 3d |eralMatr|ix* gm =|
|00004680| 20 6e 65 77 20 43 72 6f | 75 74 4d 61 74 72 69 78 | new Cro|utMatrix|
|00004690| 28 2a 74 68 69 73 29 3b | 0a 58 20 20 20 4d 61 74 |(*this);|.X   Mat|
|000046a0| 72 69 78 45 72 72 6f 72 | 4e 6f 53 70 61 63 65 28 |rixError|NoSpace(|
|000046b0| 67 6d 29 3b 20 67 6d 2d | 3e 52 65 6c 65 61 73 65 |gm); gm-|>Release|
|000046c0| 41 6e 64 44 65 6c 65 74 | 65 28 29 3b 20 72 65 74 |AndDelet|e(); ret|
|000046d0| 75 72 6e 20 67 6d 3b 0a | 58 7d 0a 58 0a 58 76 6f |urn gm;.|X}.X.Xvo|
|000046e0| 69 64 20 43 72 6f 75 74 | 4d 61 74 72 69 78 3a 3a |id Crout|Matrix::|
|000046f0| 53 6f 6c 76 65 72 28 4d | 61 74 72 69 78 52 6f 77 |Solver(M|atrixRow|
|00004700| 43 6f 6c 26 20 6d 63 6f | 75 74 2c 20 63 6f 6e 73 |Col& mco|ut, cons|
|00004710| 74 20 4d 61 74 72 69 78 | 52 6f 77 43 6f 6c 26 20 |t Matrix|RowCol& |
|00004720| 6d 63 69 6e 29 0a 58 7b | 0a 58 20 20 20 52 45 50 |mcin).X{|.X   REP|
|00004730| 4f 52 54 0a 58 20 20 20 | 72 65 61 6c 2a 20 65 6c |ORT.X   |real* el|
|00004740| 20 3d 20 6d 63 69 6e 2e | 73 74 6f 72 65 3b 20 69 | = mcin.|store; i|
|00004750| 6e 74 20 69 20 3d 20 6d | 63 69 6e 2e 73 6b 69 70 |nt i = m|cin.skip|
|00004760| 3b 0a 58 20 20 20 77 68 | 69 6c 65 20 28 69 2d 2d |;.X   wh|ile (i--|
|00004770| 29 20 2a 65 6c 2b 2b 20 | 3d 20 30 2e 30 3b 0a 58 |) *el++ |= 0.0;.X|
|00004780| 20 20 20 65 6c 20 2b 3d | 20 6d 63 69 6e 2e 73 74 |   el +=| mcin.st|
|00004790| 6f 72 61 67 65 3b 20 69 | 20 3d 20 6e 72 6f 77 73 |orage; i| = nrows|
|000047a0| 20 2d 20 6d 63 69 6e 2e | 73 6b 69 70 20 2d 20 6d | - mcin.|skip - m|
|000047b0| 63 69 6e 2e 73 74 6f 72 | 61 67 65 3b 0a 58 20 20 |cin.stor|age;.X  |
|000047c0| 20 77 68 69 6c 65 20 28 | 69 2d 2d 29 20 2a 65 6c | while (|i--) *el|
|000047d0| 2b 2b 20 3d 20 30 2e 30 | 3b 0a 58 20 20 20 6c 75 |++ = 0.0|;.X   lu|
|000047e0| 62 6b 73 62 28 6d 63 69 | 6e 2e 73 74 6f 72 65 2c |bksb(mci|n.store,|
|000047f0| 20 6d 63 6f 75 74 2e 73 | 6b 69 70 29 3b 0a 58 7d | mcout.s|kip);.X}|
|00004800| 0a 58 0a 58 0a 58 2f 2f | 20 44 6f 20 77 65 20 6e |.X.X.X//| Do we n|
|00004810| 65 65 64 20 63 68 65 63 | 6b 20 66 6f 72 20 65 6e |eed chec|k for en|
|00004820| 74 69 72 65 6c 79 20 7a | 65 72 6f 20 6f 75 74 70 |tirely z|ero outp|
|00004830| 75 74 3f 0a 58 0a 58 76 | 6f 69 64 20 55 70 70 65 |ut?.X.Xv|oid Uppe|
|00004840| 72 54 72 69 61 6e 67 75 | 6c 61 72 4d 61 74 72 69 |rTriangu|larMatri|
|00004850| 78 3a 3a 53 6f 6c 76 65 | 72 28 4d 61 74 72 69 78 |x::Solve|r(Matrix|
|00004860| 52 6f 77 43 6f 6c 26 20 | 6d 63 6f 75 74 2c 0a 58 |RowCol& |mcout,.X|
|00004870| 20 20 20 63 6f 6e 73 74 | 20 4d 61 74 72 69 78 52 |   const| MatrixR|
|00004880| 6f 77 43 6f 6c 26 20 6d | 63 69 6e 29 0a 58 7b 0a |owCol& m|cin).X{.|
|00004890| 58 20 20 20 52 45 50 4f | 52 54 0a 58 20 20 20 72 |X   REPO|RT.X   r|
|000048a0| 65 61 6c 2a 20 65 6c 78 | 20 3d 20 6d 63 69 6e 2e |eal* elx| = mcin.|
|000048b0| 73 74 6f 72 65 2b 6d 63 | 6f 75 74 2e 73 6b 69 70 |store+mc|out.skip|
|000048c0| 3b 20 69 6e 74 20 69 20 | 3d 20 6d 63 69 6e 2e 73 |; int i |= mcin.s|
|000048d0| 6b 69 70 2d 6d 63 6f 75 | 74 2e 73 6b 69 70 3b 0a |kip-mcou|t.skip;.|
|000048e0| 58 20 20 20 77 68 69 6c | 65 20 28 69 2d 2d 20 3e |X   whil|e (i-- >|
|000048f0| 20 30 29 20 2a 65 6c 78 | 2b 2b 20 3d 20 30 2e 30 | 0) *elx|++ = 0.0|
|00004900| 3b 0a 58 20 20 20 69 6e | 74 20 6e 72 20 3d 20 6d |;.X   in|t nr = m|
|00004910| 63 69 6e 2e 73 6b 69 70 | 2b 6d 63 69 6e 2e 73 74 |cin.skip|+mcin.st|
|00004920| 6f 72 61 67 65 3b 20 65 | 6c 78 20 3d 20 6d 63 69 |orage; e|lx = mci|
|00004930| 6e 2e 73 74 6f 72 65 2b | 6e 72 3b 20 72 65 61 6c |n.store+|nr; real|
|00004940| 2a 20 65 6c 20 3d 20 65 | 6c 78 3b 0a 58 20 20 20 |* el = e|lx;.X   |
|00004950| 69 6e 74 20 6a 20 3d 20 | 6d 63 6f 75 74 2e 73 6b |int j = |mcout.sk|
|00004960| 69 70 2b 6d 63 6f 75 74 | 2e 73 74 6f 72 61 67 65 |ip+mcout|.storage|
|00004970| 2d 6e 72 3b 20 69 6e 74 | 20 6e 63 20 3d 20 6e 63 |-nr; int| nc = nc|
|00004980| 6f 6c 73 2d 6e 72 3b 20 | 69 20 3d 20 6e 72 2d 6d |ols-nr; |i = nr-m|
|00004990| 63 6f 75 74 2e 73 6b 69 | 70 3b 0a 58 20 20 20 77 |cout.ski|p;.X   w|
|000049a0| 68 69 6c 65 20 28 6a 2d | 2d 20 3e 20 30 29 20 2a |hile (j-|- > 0) *|
|000049b0| 65 6c 78 2b 2b 20 3d 20 | 30 2e 30 3b 0a 58 20 20 |elx++ = |0.0;.X  |
|000049c0| 20 72 65 61 6c 2a 20 41 | 65 6c 20 3d 20 73 74 6f | real* A|el = sto|
|000049d0| 72 65 20 2b 20 28 6e 72 | 2a 28 32 2a 6e 63 6f 6c |re + (nr|*(2*ncol|
|000049e0| 73 2d 6e 72 2b 31 29 29 | 2f 32 3b 20 6a 20 3d 20 |s-nr+1))|/2; j = |
|000049f0| 30 3b 0a 58 20 20 20 77 | 68 69 6c 65 20 28 69 2d |0;.X   w|hile (i-|
|00004a00| 2d 20 3e 20 30 29 0a 58 | 20 20 20 7b 0a 58 20 20 |- > 0).X|   {.X  |
|00004a10| 20 20 20 20 65 6c 78 20 | 3d 20 65 6c 3b 20 72 65 |    elx |= el; re|
|00004a20| 61 6c 20 73 75 6d 20 3d | 20 30 2e 30 3b 20 69 6e |al sum =| 0.0; in|
|00004a30| 74 20 6a 78 20 3d 20 6a | 2b 2b 3b 20 41 65 6c 20 |t jx = j|++; Ael |
|00004a40| 2d 3d 20 6e 63 3b 0a 58 | 20 20 20 20 20 20 77 68 |-= nc;.X|      wh|
|00004a50| 69 6c 65 20 28 6a 78 2d | 2d 29 20 73 75 6d 20 2b |ile (jx-|-) sum +|
|00004a60| 3d 20 2a 28 2d 2d 41 65 | 6c 29 20 2a 20 2a 28 2d |= *(--Ae|l) * *(-|
|00004a70| 2d 65 6c 78 29 3b 0a 58 | 20 20 20 20 20 20 65 6c |-elx);.X|      el|
|00004a80| 78 2d 2d 3b 20 2a 65 6c | 78 20 3d 20 28 2a 65 6c |x--; *el|x = (*el|
|00004a90| 78 20 2d 20 73 75 6d 29 | 20 2f 20 2a 28 2d 2d 41 |x - sum)| / *(--A|
|00004aa0| 65 6c 29 3b 0a 58 20 20 | 20 7d 0a 58 7d 0a 58 0a |el);.X  | }.X}.X.|
|00004ab0| 58 76 6f 69 64 20 4c 6f | 77 65 72 54 72 69 61 6e |Xvoid Lo|werTrian|
|00004ac0| 67 75 6c 61 72 4d 61 74 | 72 69 78 3a 3a 53 6f 6c |gularMat|rix::Sol|
|00004ad0| 76 65 72 28 4d 61 74 72 | 69 78 52 6f 77 43 6f 6c |ver(Matr|ixRowCol|
|00004ae0| 26 20 6d 63 6f 75 74 2c | 0a 58 20 20 20 63 6f 6e |& mcout,|.X   con|
|00004af0| 73 74 20 4d 61 74 72 69 | 78 52 6f 77 43 6f 6c 26 |st Matri|xRowCol&|
|00004b00| 20 6d 63 69 6e 29 0a 58 | 7b 0a 58 20 20 20 52 45 | mcin).X|{.X   RE|
|00004b10| 50 4f 52 54 0a 58 20 20 | 20 72 65 61 6c 2a 20 65 |PORT.X  | real* e|
|00004b20| 6c 78 20 3d 20 6d 63 69 | 6e 2e 73 74 6f 72 65 2b |lx = mci|n.store+|
|00004b30| 6d 63 6f 75 74 2e 73 6b | 69 70 3b 20 69 6e 74 20 |mcout.sk|ip; int |
|00004b40| 69 20 3d 20 6d 63 69 6e | 2e 73 6b 69 70 2d 6d 63 |i = mcin|.skip-mc|
|00004b50| 6f 75 74 2e 73 6b 69 70 | 3b 0a 58 20 20 20 77 68 |out.skip|;.X   wh|
|00004b60| 69 6c 65 20 28 69 2d 2d | 20 3e 20 30 29 20 2a 65 |ile (i--| > 0) *e|
|00004b70| 6c 78 2b 2b 20 3d 20 30 | 2e 30 3b 0a 58 20 20 20 |lx++ = 0|.0;.X   |
|00004b80| 69 6e 74 20 6e 63 20 3d | 20 6d 63 69 6e 2e 73 6b |int nc =| mcin.sk|
|00004b90| 69 70 3b 20 69 20 3d 20 | 6e 63 2b 6d 63 69 6e 2e |ip; i = |nc+mcin.|
|00004ba0| 73 74 6f 72 61 67 65 3b | 20 65 6c 78 20 3d 20 6d |storage;| elx = m|
|00004bb0| 63 69 6e 2e 73 74 6f 72 | 65 2b 69 3b 0a 58 20 20 |cin.stor|e+i;.X  |
|00004bc0| 20 69 6e 74 20 6e 72 20 | 3d 20 6d 63 6f 75 74 2e | int nr |= mcout.|
|00004bd0| 73 6b 69 70 2b 6d 63 6f | 75 74 2e 73 74 6f 72 61 |skip+mco|ut.stora|
|00004be0| 67 65 3b 20 69 6e 74 20 | 6a 20 3d 20 6e 72 2d 69 |ge; int |j = nr-i|
|00004bf0| 3b 20 69 20 3d 20 6e 72 | 2d 6e 63 3b 0a 58 20 20 |; i = nr|-nc;.X  |
|00004c00| 20 77 68 69 6c 65 20 28 | 6a 2d 2d 20 3e 20 30 29 | while (|j-- > 0)|
|00004c10| 20 2a 65 6c 78 2b 2b 20 | 3d 20 30 2e 30 3b 0a 58 | *elx++ |= 0.0;.X|
|00004c20| 20 20 20 72 65 61 6c 2a | 20 65 6c 20 3d 20 6d 63 |   real*| el = mc|
|00004c30| 69 6e 2e 73 74 6f 72 65 | 2b 6e 63 3b 20 72 65 61 |in.store|+nc; rea|
|00004c40| 6c 2a 20 41 65 6c 20 3d | 20 73 74 6f 72 65 20 2b |l* Ael =| store +|
|00004c50| 20 28 6e 63 2a 28 6e 63 | 2b 31 29 29 2f 32 3b 20 | (nc*(nc|+1))/2; |
|00004c60| 6a 20 3d 20 30 3b 0a 58 | 20 20 20 77 68 69 6c 65 |j = 0;.X|   while|
|00004c70| 20 28 69 2d 2d 20 3e 20 | 30 29 0a 58 20 20 20 7b | (i-- > |0).X   {|
|00004c80| 0a 58 20 20 20 20 20 20 | 65 6c 78 20 3d 20 65 6c |.X      |elx = el|
|00004c90| 3b 20 72 65 61 6c 20 73 | 75 6d 20 3d 20 30 2e 30 |; real s|um = 0.0|
|00004ca0| 3b 20 69 6e 74 20 6a 78 | 20 3d 20 6a 2b 2b 3b 20 |; int jx| = j++; |
|00004cb0| 41 65 6c 20 2b 3d 20 6e | 63 3b 0a 58 20 20 20 20 |Ael += n|c;.X    |
|00004cc0| 20 20 77 68 69 6c 65 20 | 28 6a 78 2d 2d 29 20 73 |  while |(jx--) s|
|00004cd0| 75 6d 20 2b 3d 20 2a 41 | 65 6c 2b 2b 20 2a 20 2a |um += *A|el++ * *|
|00004ce0| 65 6c 78 2b 2b 3b 0a 58 | 20 20 20 20 20 20 2a 65 |elx++;.X|      *e|
|00004cf0| 6c 78 20 3d 20 28 2a 65 | 6c 78 20 2d 20 73 75 6d |lx = (*e|lx - sum|
|00004d00| 29 20 2f 20 2a 41 65 6c | 2b 2b 3b 0a 58 20 20 20 |) / *Ael|++;.X   |
|00004d10| 7d 0a 58 7d 0a 58 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a |}.X}.X.X|/*******|
|00004d20| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 20 63 61 72 |********|**** car|
|00004d30| 72 79 20 6f 75 74 20 62 | 69 6e 61 72 79 20 6f 70 |ry out b|inary op|
|00004d40| 65 72 61 74 69 6f 6e 73 | 20 2a 2a 2a 2a 2a 2a 2a |erations| *******|
|00004d50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004d60| 2a 2a 2f 0a 58 0a 58 73 | 74 61 74 69 63 20 76 6f |**/.X.Xs|tatic vo|
|00004d70| 69 64 20 41 64 64 28 47 | 65 6e 65 72 61 6c 4d 61 |id Add(G|eneralMa|
|00004d80| 74 72 69 78 2a 2c 47 65 | 6e 65 72 61 6c 4d 61 74 |trix*,Ge|neralMat|
|00004d90| 72 69 78 2a 2c 20 47 65 | 6e 65 72 61 6c 4d 61 74 |rix*, Ge|neralMat|
|00004da0| 72 69 78 2a 29 3b 0a 58 | 73 74 61 74 69 63 20 76 |rix*);.X|static v|
|00004db0| 6f 69 64 20 41 64 64 28 | 47 65 6e 65 72 61 6c 4d |oid Add(|GeneralM|
|00004dc0| 61 74 72 69 78 2a 2c 47 | 65 6e 65 72 61 6c 4d 61 |atrix*,G|eneralMa|
|00004dd0| 74 72 69 78 2a 29 3b 0a | 58 0a 58 73 74 61 74 69 |trix*);.|X.Xstati|
|00004de0| 63 20 76 6f 69 64 20 53 | 75 62 74 72 61 63 74 28 |c void S|ubtract(|
|00004df0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 2c 47 |GeneralM|atrix*,G|
|00004e00| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 2c 20 47 |eneralMa|trix*, G|
|00004e10| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 29 3b 0a |eneralMa|trix*);.|
|00004e20| 58 73 74 61 74 69 63 20 | 76 6f 69 64 20 53 75 62 |Xstatic |void Sub|
|00004e30| 74 72 61 63 74 28 47 65 | 6e 65 72 61 6c 4d 61 74 |tract(Ge|neralMat|
|00004e40| 72 69 78 2a 2c 47 65 6e | 65 72 61 6c 4d 61 74 72 |rix*,Gen|eralMatr|
|00004e50| 69 78 2a 29 3b 0a 58 73 | 74 61 74 69 63 20 76 6f |ix*);.Xs|tatic vo|
|00004e60| 69 64 20 52 65 76 65 72 | 73 65 53 75 62 74 72 61 |id Rever|seSubtra|
|00004e70| 63 74 28 47 65 6e 65 72 | 61 6c 4d 61 74 72 69 78 |ct(Gener|alMatrix|
|00004e80| 2a 2c 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |*,Genera|lMatrix*|
|00004e90| 29 3b 0a 58 0a 58 73 74 | 61 74 69 63 20 47 65 6e |);.X.Xst|atic Gen|
|00004ea0| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 47 65 6e 65 |eralMatr|ix* Gene|
|00004eb0| 72 61 6c 41 64 64 28 47 | 65 6e 65 72 61 6c 4d 61 |ralAdd(G|eneralMa|
|00004ec0| 74 72 69 78 2a 2c 47 65 | 6e 65 72 61 6c 4d 61 74 |trix*,Ge|neralMat|
|00004ed0| 72 69 78 2a 2c 4d 61 74 | 72 69 78 54 79 70 65 29 |rix*,Mat|rixType)|
|00004ee0| 3b 0a 58 73 74 61 74 69 | 63 20 47 65 6e 65 72 61 |;.Xstati|c Genera|
|00004ef0| 6c 4d 61 74 72 69 78 2a | 20 47 65 6e 65 72 61 6c |lMatrix*| General|
|00004f00| 53 75 62 28 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |Sub(Gene|ralMatri|
|00004f10| 78 2a 2c 47 65 6e 65 72 | 61 6c 4d 61 74 72 69 78 |x*,Gener|alMatrix|
|00004f20| 2a 2c 4d 61 74 72 69 78 | 54 79 70 65 29 3b 0a 58 |*,Matrix|Type);.X|
|00004f30| 73 74 61 74 69 63 20 47 | 65 6e 65 72 61 6c 4d 61 |static G|eneralMa|
|00004f40| 74 72 69 78 2a 20 47 65 | 6e 65 72 61 6c 4d 75 6c |trix* Ge|neralMul|
|00004f50| 74 28 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |t(Genera|lMatrix*|
|00004f60| 2c 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 2c |,General|Matrix*,|
|00004f70| 4d 61 74 72 69 78 54 79 | 70 65 29 3b 0a 58 73 74 |MatrixTy|pe);.Xst|
|00004f80| 61 74 69 63 20 47 65 6e | 65 72 61 6c 4d 61 74 72 |atic Gen|eralMatr|
|00004f90| 69 78 2a 20 47 65 6e 65 | 72 61 6c 53 6f 6c 76 28 |ix* Gene|ralSolv(|
|00004fa0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 2c 47 |GeneralM|atrix*,G|
|00004fb0| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 2c 4d 61 |eneralMa|trix*,Ma|
|00004fc0| 74 72 69 78 54 79 70 65 | 29 3b 0a 58 0a 58 47 65 |trixType|);.X.XGe|
|00004fd0| 6e 65 72 61 6c 4d 61 74 | 72 69 78 2a 20 41 64 64 |neralMat|rix* Add|
|00004fe0| 65 64 4d 61 74 72 69 78 | 3a 3a 45 76 61 6c 75 61 |edMatrix|::Evalua|
|00004ff0| 74 65 28 4d 61 74 72 69 | 78 54 79 70 65 20 6d 74 |te(Matri|xType mt|
|00005000| 29 0a 58 7b 20 52 45 50 | 4f 52 54 20 72 65 74 75 |).X{ REP|ORT retu|
|00005010| 72 6e 20 47 65 6e 65 72 | 61 6c 41 64 64 28 62 6d |rn Gener|alAdd(bm|
|00005020| 31 2d 3e 45 76 61 6c 75 | 61 74 65 28 29 2c 20 62 |1->Evalu|ate(), b|
|00005030| 6d 32 2d 3e 45 76 61 6c | 75 61 74 65 28 29 2c 20 |m2->Eval|uate(), |
|00005040| 6d 74 29 3b 20 7d 0a 58 | 0a 58 47 65 6e 65 72 61 |mt); }.X|.XGenera|
|00005050| 6c 4d 61 74 72 69 78 2a | 20 53 75 62 74 72 61 63 |lMatrix*| Subtrac|
|00005060| 74 65 64 4d 61 74 72 69 | 78 3a 3a 45 76 61 6c 75 |tedMatri|x::Evalu|
|00005070| 61 74 65 28 4d 61 74 72 | 69 78 54 79 70 65 20 6d |ate(Matr|ixType m|
|00005080| 74 29 0a 58 7b 20 52 45 | 50 4f 52 54 20 72 65 74 |t).X{ RE|PORT ret|
|00005090| 75 72 6e 20 47 65 6e 65 | 72 61 6c 53 75 62 28 62 |urn Gene|ralSub(b|
|000050a0| 6d 31 2d 3e 45 76 61 6c | 75 61 74 65 28 29 2c 20 |m1->Eval|uate(), |
|000050b0| 62 6d 32 2d 3e 45 76 61 | 6c 75 61 74 65 28 29 2c |bm2->Eva|luate(),|
|000050c0| 20 6d 74 29 3b 20 7d 0a | 58 0a 58 47 65 6e 65 72 | mt); }.|X.XGener|
|000050d0| 61 6c 4d 61 74 72 69 78 | 2a 20 4d 75 6c 74 69 70 |alMatrix|* Multip|
|000050e0| 6c 69 65 64 4d 61 74 72 | 69 78 3a 3a 45 76 61 6c |liedMatr|ix::Eval|
|000050f0| 75 61 74 65 28 4d 61 74 | 72 69 78 54 79 70 65 20 |uate(Mat|rixType |
|00005100| 6d 74 29 0a 58 7b 0a 58 | 20 20 20 52 45 50 4f 52 |mt).X{.X|   REPOR|
|00005110| 54 0a 58 20 20 20 47 65 | 6e 65 72 61 6c 4d 61 74 |T.X   Ge|neralMat|
|00005120| 72 69 78 2a 20 67 6d 32 | 20 3d 20 62 6d 32 2d 3e |rix* gm2| = bm2->|
|00005130| 45 76 61 6c 75 61 74 65 | 28 29 3b 0a 58 20 20 20 |Evaluate|();.X   |
|00005140| 4d 61 74 72 69 78 54 79 | 70 65 20 6d 74 73 79 6d |MatrixTy|pe mtsym|
|00005150| 28 4d 61 74 72 69 78 54 | 79 70 65 3a 3a 53 79 6d |(MatrixT|ype::Sym|
|00005160| 29 3b 0a 58 20 20 20 69 | 66 20 28 67 6d 32 2d 3e |);.X   i|f (gm2->|
|00005170| 54 79 70 65 28 29 20 3d | 3d 20 6d 74 73 79 6d 29 |Type() =|= mtsym)|
|00005180| 20 67 6d 32 20 3d 20 67 | 6d 32 2d 3e 45 76 61 6c | gm2 = g|m2->Eval|
|00005190| 75 61 74 65 28 4d 61 74 | 72 69 78 54 79 70 65 3a |uate(Mat|rixType:|
|000051a0| 3a 52 65 63 74 29 3b 0a | 58 20 20 20 72 65 74 75 |:Rect);.|X   retu|
|000051b0| 72 6e 20 47 65 6e 65 72 | 61 6c 4d 75 6c 74 28 62 |rn Gener|alMult(b|
|000051c0| 6d 31 2d 3e 45 76 61 6c | 75 61 74 65 28 29 2c 20 |m1->Eval|uate(), |
|000051d0| 67 6d 32 2c 20 6d 74 29 | 3b 0a 58 7d 0a 58 0a 58 |gm2, mt)|;.X}.X.X|
|000051e0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 53 |GeneralM|atrix* S|
|000051f0| 6f 6c 76 65 64 4d 61 74 | 72 69 78 3a 3a 45 76 61 |olvedMat|rix::Eva|
|00005200| 6c 75 61 74 65 28 4d 61 | 74 72 69 78 54 79 70 65 |luate(Ma|trixType|
|00005210| 20 6d 74 29 0a 58 7b 20 | 52 45 50 4f 52 54 20 72 | mt).X{ |REPORT r|
|00005220| 65 74 75 72 6e 20 47 65 | 6e 65 72 61 6c 53 6f 6c |eturn Ge|neralSol|
|00005230| 76 28 62 6d 31 2d 3e 45 | 76 61 6c 75 61 74 65 28 |v(bm1->E|valuate(|
|00005240| 29 2c 20 62 6d 32 2d 3e | 45 76 61 6c 75 61 74 65 |), bm2->|Evaluate|
|00005250| 28 29 2c 20 6d 74 29 3b | 20 7d 0a 58 0a 58 2f 2f |(), mt);| }.X.X//|
|00005260| 20 72 6f 75 74 69 6e 65 | 73 20 66 6f 72 20 61 64 | routine|s for ad|
|00005270| 64 69 6e 67 20 6f 72 20 | 73 75 62 74 72 61 63 74 |ding or |subtract|
|00005280| 69 6e 67 20 6d 61 74 72 | 69 63 65 73 20 6f 66 20 |ing matr|ices of |
|00005290| 69 64 65 6e 74 69 63 61 | 6c 20 73 74 6f 72 61 67 |identica|l storag|
|000052a0| 65 20 73 74 72 75 63 74 | 75 72 65 0a 58 0a 58 73 |e struct|ure.X.Xs|
|000052b0| 74 61 74 69 63 20 76 6f | 69 64 20 41 64 64 28 47 |tatic vo|id Add(G|
|000052c0| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 20 67 6d |eneralMa|trix* gm|
|000052d0| 2c 20 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |, Genera|lMatrix*|
|000052e0| 20 67 6d 31 2c 20 47 65 | 6e 65 72 61 6c 4d 61 74 | gm1, Ge|neralMat|
|000052f0| 72 69 78 2a 20 67 6d 32 | 29 0a 58 7b 0a 58 20 20 |rix* gm2|).X{.X  |
|00005300| 20 52 45 50 4f 52 54 0a | 58 20 20 20 72 65 67 69 | REPORT.|X   regi|
|00005310| 73 74 65 72 20 72 65 61 | 6c 2a 20 73 31 3d 67 6d |ster rea|l* s1=gm|
|00005320| 31 2d 3e 53 74 6f 72 65 | 28 29 3b 20 72 65 67 69 |1->Store|(); regi|
|00005330| 73 74 65 72 20 72 65 61 | 6c 2a 20 73 32 3d 67 6d |ster rea|l* s2=gm|
|00005340| 32 2d 3e 53 74 6f 72 65 | 28 29 3b 0a 58 20 20 20 |2->Store|();.X   |
|00005350| 72 65 67 69 73 74 65 72 | 20 72 65 61 6c 2a 20 73 |register| real* s|
|00005360| 3d 67 6d 2d 3e 53 74 6f | 72 65 28 29 3b 20 72 65 |=gm->Sto|re(); re|
|00005370| 67 69 73 74 65 72 20 69 | 6e 74 20 69 3d 67 6d 2d |gister i|nt i=gm-|
|00005380| 3e 53 74 6f 72 61 67 65 | 28 29 3b 0a 58 20 20 20 |>Storage|();.X   |
|00005390| 77 68 69 6c 65 20 28 69 | 2d 2d 29 20 2a 73 2b 2b |while (i|--) *s++|
|000053a0| 20 3d 20 2a 73 31 2b 2b | 20 2b 20 2a 73 32 2b 2b | = *s1++| + *s2++|
|000053b0| 3b 0a 58 7d 0a 58 20 20 | 20 0a 58 73 74 61 74 69 |;.X}.X  | .Xstati|
|000053c0| 63 20 76 6f 69 64 20 41 | 64 64 28 47 65 6e 65 72 |c void A|dd(Gener|
|000053d0| 61 6c 4d 61 74 72 69 78 | 2a 20 67 6d 2c 20 47 65 |alMatrix|* gm, Ge|
|000053e0| 6e 65 72 61 6c 4d 61 74 | 72 69 78 2a 20 67 6d 32 |neralMat|rix* gm2|
|000053f0| 29 0a 58 7b 0a 58 20 20 | 20 52 45 50 4f 52 54 0a |).X{.X  | REPORT.|
|00005400| 58 20 20 20 72 65 67 69 | 73 74 65 72 20 72 65 61 |X   regi|ster rea|
|00005410| 6c 2a 20 73 32 3d 67 6d | 32 2d 3e 53 74 6f 72 65 |l* s2=gm|2->Store|
|00005420| 28 29 3b 0a 58 20 20 20 | 72 65 67 69 73 74 65 72 |();.X   |register|
|00005430| 20 72 65 61 6c 2a 20 73 | 3d 67 6d 2d 3e 53 74 6f | real* s|=gm->Sto|
|00005440| 72 65 28 29 3b 20 72 65 | 67 69 73 74 65 72 20 69 |re(); re|gister i|
|00005450| 6e 74 20 69 3d 67 6d 2d | 3e 53 74 6f 72 61 67 65 |nt i=gm-|>Storage|
|00005460| 28 29 3b 0a 58 20 20 20 | 77 68 69 6c 65 20 28 69 |();.X   |while (i|
|00005470| 2d 2d 29 20 2a 73 2b 2b | 20 2b 3d 20 2a 73 32 2b |--) *s++| += *s2+|
|00005480| 2b 3b 0a 58 7d 0a 58 0a | 58 73 74 61 74 69 63 20 |+;.X}.X.|Xstatic |
|00005490| 76 6f 69 64 20 53 75 62 | 74 72 61 63 74 28 47 65 |void Sub|tract(Ge|
|000054a0| 6e 65 72 61 6c 4d 61 74 | 72 69 78 2a 20 67 6d 2c |neralMat|rix* gm,|
|000054b0| 20 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 | General|Matrix* |
|000054c0| 67 6d 31 2c 20 47 65 6e | 65 72 61 6c 4d 61 74 72 |gm1, Gen|eralMatr|
|000054d0| 69 78 2a 20 67 6d 32 29 | 0a 58 7b 0a 58 20 20 20 |ix* gm2)|.X{.X   |
|000054e0| 52 45 50 4f 52 54 0a 58 | 20 20 20 72 65 67 69 73 |REPORT.X|   regis|
|000054f0| 74 65 72 20 72 65 61 6c | 2a 20 73 31 3d 67 6d 31 |ter real|* s1=gm1|
|00005500| 2d 3e 53 74 6f 72 65 28 | 29 3b 20 72 65 67 69 73 |->Store(|); regis|
|00005510| 74 65 72 20 72 65 61 6c | 2a 20 73 32 3d 67 6d 32 |ter real|* s2=gm2|
|00005520| 2d 3e 53 74 6f 72 65 28 | 29 3b 0a 58 20 20 20 72 |->Store(|);.X   r|
|00005530| 65 67 69 73 74 65 72 20 | 72 65 61 6c 2a 20 73 3d |egister |real* s=|
|00005540| 67 6d 2d 3e 53 74 6f 72 | 65 28 29 3b 20 72 65 67 |gm->Stor|e(); reg|
|00005550| 69 73 74 65 72 20 69 6e | 74 20 69 3d 67 6d 2d 3e |ister in|t i=gm->|
|00005560| 53 74 6f 72 61 67 65 28 | 29 3b 0a 58 20 20 20 77 |Storage(|);.X   w|
|00005570| 68 69 6c 65 20 28 69 2d | 2d 29 20 2a 73 2b 2b 20 |hile (i-|-) *s++ |
|00005580| 3d 20 2a 73 31 2b 2b 20 | 2d 20 2a 73 32 2b 2b 3b |= *s1++ |- *s2++;|
|00005590| 0a 58 7d 0a 58 0a 58 73 | 74 61 74 69 63 20 76 6f |.X}.X.Xs|tatic vo|
|000055a0| 69 64 20 53 75 62 74 72 | 61 63 74 28 47 65 6e 65 |id Subtr|act(Gene|
|000055b0| 72 61 6c 4d 61 74 72 69 | 78 2a 20 67 6d 2c 20 47 |ralMatri|x* gm, G|
|000055c0| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 20 67 6d |eneralMa|trix* gm|
|000055d0| 32 29 0a 58 7b 0a 58 20 | 20 20 52 45 50 4f 52 54 |2).X{.X |  REPORT|
|000055e0| 0a 58 20 20 20 72 65 67 | 69 73 74 65 72 20 72 65 |.X   reg|ister re|
|000055f0| 61 6c 2a 20 73 32 3d 67 | 6d 32 2d 3e 53 74 6f 72 |al* s2=g|m2->Stor|
|00005600| 65 28 29 3b 0a 58 20 20 | 20 72 65 67 69 73 74 65 |e();.X  | registe|
|00005610| 72 20 72 65 61 6c 2a 20 | 73 3d 67 6d 2d 3e 53 74 |r real* |s=gm->St|
|00005620| 6f 72 65 28 29 3b 20 72 | 65 67 69 73 74 65 72 20 |ore(); r|egister |
|00005630| 69 6e 74 20 69 3d 67 6d | 2d 3e 53 74 6f 72 61 67 |int i=gm|->Storag|
|00005640| 65 28 29 3b 0a 58 20 20 | 20 77 68 69 6c 65 20 28 |e();.X  | while (|
|00005650| 69 2d 2d 29 20 2a 73 2b | 2b 20 2d 3d 20 2a 73 32 |i--) *s+|+ -= *s2|
|00005660| 2b 2b 3b 0a 58 7d 0a 58 | 0a 58 73 74 61 74 69 63 |++;.X}.X|.Xstatic|
|00005670| 20 76 6f 69 64 20 52 65 | 76 65 72 73 65 53 75 62 | void Re|verseSub|
|00005680| 74 72 61 63 74 28 47 65 | 6e 65 72 61 6c 4d 61 74 |tract(Ge|neralMat|
|00005690| 72 69 78 2a 20 67 6d 2c | 20 47 65 6e 65 72 61 6c |rix* gm,| General|
|000056a0| 4d 61 74 72 69 78 2a 20 | 67 6d 32 29 0a 58 7b 0a |Matrix* |gm2).X{.|
|000056b0| 58 20 20 20 52 45 50 4f | 52 54 0a 58 20 20 20 72 |X   REPO|RT.X   r|
|000056c0| 65 67 69 73 74 65 72 20 | 72 65 61 6c 2a 20 73 32 |egister |real* s2|
|000056d0| 3d 67 6d 32 2d 3e 53 74 | 6f 72 65 28 29 3b 0a 58 |=gm2->St|ore();.X|
|000056e0| 20 20 20 72 65 67 69 73 | 74 65 72 20 72 65 61 6c |   regis|ter real|
|000056f0| 2a 20 73 3d 67 6d 2d 3e | 53 74 6f 72 65 28 29 3b |* s=gm->|Store();|
|00005700| 20 72 65 67 69 73 74 65 | 72 20 69 6e 74 20 69 3d | registe|r int i=|
|00005710| 67 6d 2d 3e 53 74 6f 72 | 61 67 65 28 29 3b 0a 58 |gm->Stor|age();.X|
|00005720| 20 20 20 77 68 69 6c 65 | 20 28 69 2d 2d 29 20 7b |   while| (i--) {|
|00005730| 20 2a 73 20 3d 20 2a 73 | 32 2b 2b 20 2d 20 2a 73 | *s = *s|2++ - *s|
|00005740| 3b 20 73 2b 2b 3b 20 7d | 0a 58 7d 0a 58 0a 58 73 |; s++; }|.X}.X.Xs|
|00005750| 74 61 74 69 63 20 47 65 | 6e 65 72 61 6c 4d 61 74 |tatic Ge|neralMat|
|00005760| 72 69 78 2a 20 47 65 6e | 65 72 61 6c 41 64 64 28 |rix* Gen|eralAdd(|
|00005770| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00005780| 6d 31 2c 20 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |m1, Gene|ralMatri|
|00005790| 78 2a 20 67 6d 32 2c 0a | 58 20 20 20 4d 61 74 72 |x* gm2,.|X   Matr|
|000057a0| 69 78 54 79 70 65 20 6d | 74 78 29 0a 58 7b 0a 58 |ixType m|tx).X{.X|
|000057b0| 20 20 20 69 6e 74 20 6e | 72 3d 67 6d 31 2d 3e 4e |   int n|r=gm1->N|
|000057c0| 72 6f 77 73 28 29 3b 20 | 69 6e 74 20 6e 63 3d 67 |rows(); |int nc=g|
|000057d0| 6d 31 2d 3e 4e 63 6f 6c | 73 28 29 3b 0a 58 20 20 |m1->Ncol|s();.X  |
|000057e0| 20 69 66 20 28 6e 72 21 | 3d 67 6d 32 2d 3e 4e 72 | if (nr!|=gm2->Nr|
|000057f0| 6f 77 73 28 29 20 7c 7c | 20 6e 63 21 3d 67 6d 32 |ows() ||| nc!=gm2|
|00005800| 2d 3e 4e 63 6f 6c 73 28 | 29 29 0a 58 20 20 20 20 |->Ncols(|)).X    |
|00005810| 20 20 4d 61 74 72 69 78 | 45 72 72 6f 72 28 22 49 |  Matrix|Error("I|
|00005820| 6e 63 6f 6d 70 61 74 69 | 62 6c 65 20 64 69 6d 65 |ncompati|ble dime|
|00005830| 6e 73 69 6f 6e 73 22 29 | 3b 0a 58 20 20 20 69 66 |nsions")|;.X   if|
|00005840| 20 28 21 6d 74 78 29 20 | 6d 74 78 20 3d 20 67 6d | (!mtx) |mtx = gm|
|00005850| 31 2d 3e 54 79 70 65 28 | 29 20 2b 20 67 6d 32 2d |1->Type(|) + gm2-|
|00005860| 3e 54 79 70 65 28 29 3b | 0a 58 20 20 20 69 66 20 |>Type();|.X   if |
|00005870| 28 6d 74 78 20 3d 3d 20 | 67 6d 31 2d 3e 54 79 70 |(mtx == |gm1->Typ|
|00005880| 65 28 29 20 26 26 20 6d | 74 78 20 3d 3d 20 67 6d |e() && m|tx == gm|
|00005890| 32 2d 3e 54 79 70 65 28 | 29 29 0a 58 09 09 09 09 |2->Type(|)).X....|
|000058a0| 09 20 20 20 20 20 2f 2f | 20 6d 6f 64 69 66 79 20 |.     //| modify |
|000058b0| 77 68 65 6e 20 62 61 6e | 64 20 6f 72 20 73 70 61 |when ban|d or spa|
|000058c0| 72 73 65 0a 58 09 09 09 | 09 09 20 20 20 20 20 2f |rse.X...|..     /|
|000058d0| 2f 20 6d 61 74 72 69 63 | 65 73 20 61 72 65 20 69 |/ matric|es are i|
|000058e0| 6e 63 6c 75 64 65 64 2e | 0a 58 20 20 20 7b 0a 58 |ncluded.|.X   {.X|
|000058f0| 20 20 20 20 20 20 69 66 | 20 28 67 6d 31 2d 3e 72 |      if| (gm1->r|
|00005900| 65 75 73 65 28 29 29 20 | 7b 20 52 45 50 4f 52 54 |euse()) |{ REPORT|
|00005910| 20 41 64 64 28 67 6d 31 | 2c 67 6d 32 29 3b 20 67 | Add(gm1|,gm2); g|
|00005920| 6d 32 2d 3e 74 44 65 6c | 65 74 65 28 29 3b 20 72 |m2->tDel|ete(); r|
|00005930| 65 74 75 72 6e 20 67 6d | 31 3b 20 7d 0a 58 20 20 |eturn gm|1; }.X  |
|00005940| 20 20 20 20 65 6c 73 65 | 20 69 66 20 28 67 6d 32 |    else| if (gm2|
|00005950| 2d 3e 72 65 75 73 65 28 | 29 29 20 7b 20 52 45 50 |->reuse(|)) { REP|
|00005960| 4f 52 54 20 41 64 64 28 | 67 6d 32 2c 67 6d 31 29 |ORT Add(|gm2,gm1)|
|00005970| 3b 20 72 65 74 75 72 6e | 20 67 6d 32 3b 20 7d 0a |; return| gm2; }.|
|00005980| 58 20 20 20 20 20 20 65 | 6c 73 65 0a 58 20 20 20 |X      e|lse.X   |
|00005990| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 20 20 20 52 |   {.X  |       R|
|000059a0| 45 50 4f 52 54 0a 58 20 | 20 20 20 20 20 20 20 20 |EPORT.X |        |
|000059b0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|000059c0| 6d 78 20 3d 20 67 6d 31 | 2d 3e 54 79 70 65 28 29 |mx = gm1|->Type()|
|000059d0| 2e 4e 65 77 28 6e 72 2c | 6e 63 29 3b 20 67 6d 78 |.New(nr,|nc); gmx|
|000059e0| 2d 3e 52 65 6c 65 61 73 | 65 41 6e 64 44 65 6c 65 |->Releas|eAndDele|
|000059f0| 74 65 28 29 3b 0a 58 20 | 20 20 20 20 20 20 20 20 |te();.X |        |
|00005a00| 67 6d 78 2d 3e 52 65 6c | 65 61 73 65 41 6e 64 44 |gmx->Rel|easeAndD|
|00005a10| 65 6c 65 74 65 28 29 3b | 20 41 64 64 28 67 6d 78 |elete();| Add(gmx|
|00005a20| 2c 67 6d 31 2c 67 6d 32 | 29 3b 20 72 65 74 75 72 |,gm1,gm2|); retur|
|00005a30| 6e 20 67 6d 78 3b 0a 58 | 20 20 20 20 20 20 7d 0a |n gmx;.X|      }.|
|00005a40| 58 20 20 20 7d 0a 58 20 | 20 20 65 6c 73 65 0a 58 |X   }.X |  else.X|
|00005a50| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 69 66 20 28 |   {.X  |    if (|
|00005a60| 67 6d 31 2d 3e 54 79 70 | 65 28 29 3d 3d 6d 74 78 |gm1->Typ|e()==mtx|
|00005a70| 20 26 26 20 67 6d 31 2d | 3e 72 65 75 73 65 28 29 | && gm1-|>reuse()|
|00005a80| 20 29 20 20 20 20 2f 2f | 20 6d 75 73 74 20 68 61 | )    //| must ha|
|00005a90| 76 65 20 74 79 70 65 20 | 74 65 73 74 20 66 69 72 |ve type |test fir|
|00005aa0| 73 74 0a 58 20 20 20 20 | 20 20 7b 0a 58 09 20 52 |st.X    |  {.X. R|
|00005ab0| 45 50 4f 52 54 0a 58 09 | 20 4d 61 74 72 69 78 52 |EPORT.X.| MatrixR|
|00005ac0| 6f 77 20 6d 72 31 28 67 | 6d 31 2c 20 53 74 6f 72 |ow mr1(g|m1, Stor|
|00005ad0| 65 4f 6e 45 78 69 74 2b | 4c 6f 61 64 4f 6e 45 6e |eOnExit+|LoadOnEn|
|00005ae0| 74 72 79 2b 44 69 72 65 | 63 74 50 61 72 74 29 3b |try+Dire|ctPart);|
|00005af0| 0a 58 09 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 |.X. Matr|ixRow mr|
|00005b00| 32 28 67 6d 32 2c 20 4c | 6f 61 64 4f 6e 45 6e 74 |2(gm2, L|oadOnEnt|
|00005b10| 72 79 29 3b 0a 58 09 20 | 77 68 69 6c 65 20 28 6e |ry);.X. |while (n|
|00005b20| 72 2d 2d 29 20 7b 20 6d | 72 31 2e 41 64 64 28 6d |r--) { m|r1.Add(m|
|00005b30| 72 32 29 3b 20 6d 72 31 | 2e 4e 65 78 74 28 29 3b |r2); mr1|.Next();|
|00005b40| 20 6d 72 32 2e 4e 65 78 | 74 28 29 3b 20 7d 0a 58 | mr2.Nex|t(); }.X|
|00005b50| 20 20 20 20 20 20 20 20 | 20 67 6d 32 2d 3e 74 44 |        | gm2->tD|
|00005b60| 65 6c 65 74 65 28 29 3b | 20 72 65 74 75 72 6e 20 |elete();| return |
|00005b70| 67 6d 31 3b 0a 58 20 20 | 20 20 20 20 7d 0a 58 20 |gm1;.X  |    }.X |
|00005b80| 20 20 20 20 20 65 6c 73 | 65 20 69 66 20 28 67 6d |     els|e if (gm|
|00005b90| 32 2d 3e 54 79 70 65 28 | 29 3d 3d 6d 74 78 20 26 |2->Type(|)==mtx &|
|00005ba0| 26 20 67 6d 32 2d 3e 72 | 65 75 73 65 28 29 20 29 |& gm2->r|euse() )|
|00005bb0| 0a 58 20 20 20 20 20 20 | 7b 0a 58 20 20 20 20 20 |.X      |{.X     |
|00005bc0| 20 20 20 20 52 45 50 4f | 52 54 0a 58 09 20 4d 61 |    REPO|RT.X. Ma|
|00005bd0| 74 72 69 78 52 6f 77 20 | 6d 72 31 28 67 6d 31 2c |trixRow |mr1(gm1,|
|00005be0| 20 4c 6f 61 64 4f 6e 45 | 6e 74 72 79 29 3b 0a 58 | LoadOnE|ntry);.X|
|00005bf0| 09 20 4d 61 74 72 69 78 | 52 6f 77 20 6d 72 32 28 |. Matrix|Row mr2(|
|00005c00| 67 6d 32 2c 20 53 74 6f | 72 65 4f 6e 45 78 69 74 |gm2, Sto|reOnExit|
|00005c10| 2b 4c 6f 61 64 4f 6e 45 | 6e 74 72 79 2b 44 69 72 |+LoadOnE|ntry+Dir|
|00005c20| 65 63 74 50 61 72 74 29 | 3b 0a 58 09 20 77 68 69 |ectPart)|;.X. whi|
|00005c30| 6c 65 20 28 6e 72 2d 2d | 29 20 7b 20 6d 72 32 2e |le (nr--|) { mr2.|
|00005c40| 41 64 64 28 6d 72 31 29 | 3b 20 6d 72 31 2e 4e 65 |Add(mr1)|; mr1.Ne|
|00005c50| 78 74 28 29 3b 20 6d 72 | 32 2e 4e 65 78 74 28 29 |xt(); mr|2.Next()|
|00005c60| 3b 20 7d 0a 58 09 20 69 | 66 20 28 67 6d 31 2d 3e |; }.X. i|f (gm1->|
|00005c70| 54 79 70 65 28 29 21 3d | 6d 74 78 29 20 67 6d 31 |Type()!=|mtx) gm1|
|00005c80| 2d 3e 74 44 65 6c 65 74 | 65 28 29 3b 0a 58 20 20 |->tDelet|e();.X  |
|00005c90| 20 20 20 20 20 20 20 72 | 65 74 75 72 6e 20 67 6d |       r|eturn gm|
|00005ca0| 32 3b 0a 58 20 20 20 20 | 20 20 7d 0a 58 20 20 20 |2;.X    |  }.X   |
|00005cb0| 20 20 20 65 6c 73 65 0a | 58 20 20 20 20 20 20 7b |   else.|X      {|
|00005cc0| 0a 58 20 20 20 20 20 20 | 20 20 20 52 45 50 4f 52 |.X      |   REPOR|
|00005cd0| 54 0a 58 09 20 47 65 6e | 65 72 61 6c 4d 61 74 72 |T.X. Gen|eralMatr|
|00005ce0| 69 78 2a 20 67 6d 78 20 | 3d 20 6d 74 78 2e 4e 65 |ix* gmx |= mtx.Ne|
|00005cf0| 77 28 6e 72 2c 6e 63 29 | 3b 0a 58 20 20 20 20 20 |w(nr,nc)|;.X     |
|00005d00| 20 20 20 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 |    Matr|ixRow mr|
|00005d10| 31 28 67 6d 31 2c 20 4c | 6f 61 64 4f 6e 45 6e 74 |1(gm1, L|oadOnEnt|
|00005d20| 72 79 29 3b 0a 58 09 20 | 4d 61 74 72 69 78 52 6f |ry);.X. |MatrixRo|
|00005d30| 77 20 6d 72 32 28 67 6d | 32 2c 20 4c 6f 61 64 4f |w mr2(gm|2, LoadO|
|00005d40| 6e 45 6e 74 72 79 29 3b | 0a 58 09 20 4d 61 74 72 |nEntry);|.X. Matr|
|00005d50| 69 78 52 6f 77 20 6d 72 | 78 28 67 6d 78 2c 20 53 |ixRow mr|x(gmx, S|
|00005d60| 74 6f 72 65 4f 6e 45 78 | 69 74 2b 44 69 72 65 63 |toreOnEx|it+Direc|
|00005d70| 74 50 61 72 74 29 3b 0a | 58 09 20 77 68 69 6c 65 |tPart);.|X. while|
|00005d80| 20 28 6e 72 2d 2d 29 0a | 58 09 20 7b 20 6d 72 78 | (nr--).|X. { mrx|
|00005d90| 2e 41 64 64 28 6d 72 31 | 2c 6d 72 32 29 3b 20 6d |.Add(mr1|,mr2); m|
|00005da0| 72 31 2e 4e 65 78 74 28 | 29 3b 20 6d 72 32 2e 4e |r1.Next(|); mr2.N|
|00005db0| 65 78 74 28 29 3b 20 6d | 72 78 2e 4e 65 78 74 28 |ext(); m|rx.Next(|
|00005dc0| 29 3b 20 7d 0a 58 09 20 | 69 66 20 28 67 6d 31 2d |); }.X. |if (gm1-|
|00005dd0| 3e 54 79 70 65 28 29 21 | 3d 6d 74 78 29 20 67 6d |>Type()!|=mtx) gm|
|00005de0| 31 2d 3e 74 44 65 6c 65 | 74 65 28 29 3b 0a 58 09 |1->tDele|te();.X.|
|00005df0| 20 69 66 20 28 67 6d 32 | 2d 3e 54 79 70 65 28 29 | if (gm2|->Type()|
|00005e00| 21 3d 6d 74 78 29 20 67 | 6d 32 2d 3e 74 44 65 6c |!=mtx) g|m2->tDel|
|00005e10| 65 74 65 28 29 3b 0a 58 | 20 20 20 20 20 20 20 20 |ete();.X|        |
|00005e20| 20 67 6d 78 2d 3e 52 65 | 6c 65 61 73 65 41 6e 64 | gmx->Re|leaseAnd|
|00005e30| 44 65 6c 65 74 65 28 29 | 3b 20 72 65 74 75 72 6e |Delete()|; return|
|00005e40| 20 67 6d 78 3b 0a 58 20 | 20 20 20 20 20 7d 0a 58 | gmx;.X |     }.X|
|00005e50| 20 20 20 7d 0a 58 7d 0a | 58 0a 58 0a 58 73 74 61 |   }.X}.|X.X.Xsta|
|00005e60| 74 69 63 20 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |tic Gene|ralMatri|
|00005e70| 78 2a 20 47 65 6e 65 72 | 61 6c 53 75 62 28 47 65 |x* Gener|alSub(Ge|
|00005e80| 6e 65 72 61 6c 4d 61 74 | 72 69 78 2a 20 67 6d 31 |neralMat|rix* gm1|
|00005e90| 2c 20 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |, Genera|lMatrix*|
|00005ea0| 20 67 6d 32 2c 0a 58 20 | 20 20 4d 61 74 72 69 78 | gm2,.X |  Matrix|
|00005eb0| 54 79 70 65 20 6d 74 78 | 29 0a 58 7b 0a 58 20 20 |Type mtx|).X{.X  |
|00005ec0| 20 69 66 20 28 21 6d 74 | 78 29 20 6d 74 78 20 3d | if (!mt|x) mtx =|
|00005ed0| 20 67 6d 31 2d 3e 54 79 | 70 65 28 29 20 2d 20 67 | gm1->Ty|pe() - g|
|00005ee0| 6d 32 2d 3e 54 79 70 65 | 28 29 3b 0a 58 20 20 20 |m2->Type|();.X   |
|00005ef0| 69 6e 74 20 6e 72 3d 67 | 6d 31 2d 3e 4e 72 6f 77 |int nr=g|m1->Nrow|
|00005f00| 73 28 29 3b 20 69 6e 74 | 20 6e 63 3d 67 6d 31 2d |s(); int| nc=gm1-|
|00005f10| 3e 4e 63 6f 6c 73 28 29 | 3b 0a 58 20 20 20 69 66 |>Ncols()|;.X   if|
|00005f20| 20 28 6e 72 21 3d 67 6d | 32 2d 3e 4e 72 6f 77 73 | (nr!=gm|2->Nrows|
|00005f30| 28 29 20 7c 7c 20 6e 63 | 21 3d 67 6d 32 2d 3e 4e |() || nc|!=gm2->N|
|00005f40| 63 6f 6c 73 28 29 29 0a | 58 20 20 20 20 20 20 4d |cols()).|X      M|
|00005f50| 61 74 72 69 78 45 72 72 | 6f 72 28 22 49 6e 63 6f |atrixErr|or("Inco|
|00005f60| 6d 70 61 74 69 62 6c 65 | 20 64 69 6d 65 6e 73 69 |mpatible| dimensi|
|00005f70| 6f 6e 73 22 29 3b 0a 58 | 20 20 20 69 66 20 28 6d |ons");.X|   if (m|
|00005f80| 74 78 20 3d 3d 20 67 6d | 31 2d 3e 54 79 70 65 28 |tx == gm|1->Type(|
|00005f90| 29 20 26 26 20 6d 74 78 | 20 3d 3d 20 67 6d 32 2d |) && mtx| == gm2-|
|00005fa0| 3e 54 79 70 65 28 29 29 | 0a 58 09 09 09 09 09 20 |>Type())|.X..... |
|00005fb0| 20 20 20 20 2f 2f 20 6d | 6f 64 69 66 79 20 77 68 |    // m|odify wh|
|00005fc0| 65 6e 20 62 61 6e 64 20 | 6f 72 20 73 70 61 72 73 |en band |or spars|
|00005fd0| 65 0a 58 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |e.X     |        |
|00005fe0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|00005ff0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|00006000| 2f 2f 20 6d 61 74 72 69 | 63 65 73 20 61 72 65 20 |// matri|ces are |
|00006010| 69 6e 63 6c 75 64 65 64 | 2e 0a 58 20 20 20 7b 0a |included|..X   {.|
|00006020| 58 20 20 20 20 20 20 69 | 66 20 28 67 6d 31 2d 3e |X      i|f (gm1->|
|00006030| 72 65 75 73 65 28 29 29 | 0a 58 20 20 20 20 20 20 |reuse())|.X      |
|00006040| 7b 20 52 45 50 4f 52 54 | 20 53 75 62 74 72 61 63 |{ REPORT| Subtrac|
|00006050| 74 28 67 6d 31 2c 67 6d | 32 29 3b 20 67 6d 32 2d |t(gm1,gm|2); gm2-|
|00006060| 3e 74 44 65 6c 65 74 65 | 28 29 3b 20 72 65 74 75 |>tDelete|(); retu|
|00006070| 72 6e 20 67 6d 31 3b 20 | 7d 0a 58 20 20 20 20 20 |rn gm1; |}.X     |
|00006080| 20 65 6c 73 65 20 69 66 | 20 28 67 6d 32 2d 3e 72 | else if| (gm2->r|
|00006090| 65 75 73 65 28 29 29 20 | 7b 20 52 45 50 4f 52 54 |euse()) |{ REPORT|
|000060a0| 20 52 65 76 65 72 73 65 | 53 75 62 74 72 61 63 74 | Reverse|Subtract|
|000060b0| 28 67 6d 32 2c 67 6d 31 | 29 3b 20 72 65 74 75 72 |(gm2,gm1|); retur|
|000060c0| 6e 20 67 6d 32 3b 20 7d | 0a 58 20 20 20 20 20 20 |n gm2; }|.X      |
|000060d0| 65 6c 73 65 0a 58 20 20 | 20 20 20 20 7b 0a 58 20 |else.X  |    {.X |
|000060e0| 20 20 20 20 20 20 20 20 | 52 45 50 4f 52 54 0a 58 |        |REPORT.X|
|000060f0| 09 20 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |. Genera|lMatrix*|
|00006100| 20 67 6d 78 20 3d 20 67 | 6d 31 2d 3e 54 79 70 65 | gmx = g|m1->Type|
|00006110| 28 29 2e 4e 65 77 28 6e | 72 2c 6e 63 29 3b 20 67 |().New(n|r,nc); g|
|00006120| 6d 78 2d 3e 52 65 6c 65 | 61 73 65 41 6e 64 44 65 |mx->Rele|aseAndDe|
|00006130| 6c 65 74 65 28 29 3b 0a | 58 20 20 20 20 20 20 20 |lete();.|X       |
|00006140| 20 20 67 6d 78 2d 3e 52 | 65 6c 65 61 73 65 41 6e |  gmx->R|eleaseAn|
|00006150| 64 44 65 6c 65 74 65 28 | 29 3b 20 53 75 62 74 72 |dDelete(|); Subtr|
|00006160| 61 63 74 28 67 6d 78 2c | 67 6d 31 2c 67 6d 32 29 |act(gmx,|gm1,gm2)|
|00006170| 3b 20 72 65 74 75 72 6e | 20 67 6d 78 3b 0a 58 20 |; return| gmx;.X |
|00006180| 20 20 20 20 20 7d 0a 58 | 20 20 20 7d 0a 58 20 20 |     }.X|   }.X  |
|00006190| 20 65 6c 73 65 0a 58 20 | 20 20 7b 0a 58 20 20 20 | else.X |  {.X   |
|000061a0| 20 20 20 69 66 20 28 67 | 6d 31 2d 3e 54 79 70 65 |   if (g|m1->Type|
|000061b0| 28 29 20 3d 3d 20 6d 74 | 78 20 26 26 20 67 6d 31 |() == mt|x && gm1|
|000061c0| 2d 3e 72 65 75 73 65 28 | 29 20 29 0a 58 20 20 20 |->reuse(|) ).X   |
|000061d0| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 20 20 20 52 |   {.X  |       R|
|000061e0| 45 50 4f 52 54 0a 58 09 | 20 4d 61 74 72 69 78 52 |EPORT.X.| MatrixR|
|000061f0| 6f 77 20 6d 72 31 28 67 | 6d 31 2c 20 4c 6f 61 64 |ow mr1(g|m1, Load|
|00006200| 4f 6e 45 6e 74 72 79 2b | 53 74 6f 72 65 4f 6e 45 |OnEntry+|StoreOnE|
|00006210| 78 69 74 2b 44 69 72 65 | 63 74 50 61 72 74 29 3b |xit+Dire|ctPart);|
|00006220| 0a 58 09 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 |.X. Matr|ixRow mr|
|00006230| 32 28 67 6d 32 2c 20 4c | 6f 61 64 4f 6e 45 6e 74 |2(gm2, L|oadOnEnt|
|00006240| 72 79 29 3b 0a 58 09 20 | 77 68 69 6c 65 20 28 6e |ry);.X. |while (n|
|00006250| 72 2d 2d 29 20 7b 20 6d | 72 31 2e 53 75 62 28 6d |r--) { m|r1.Sub(m|
|00006260| 72 32 29 3b 20 6d 72 31 | 2e 4e 65 78 74 28 29 3b |r2); mr1|.Next();|
|00006270| 20 6d 72 32 2e 4e 65 78 | 74 28 29 3b 20 7d 0a 58 | mr2.Nex|t(); }.X|
|00006280| 20 20 20 20 20 20 20 20 | 20 67 6d 32 2d 3e 74 44 |        | gm2->tD|
|00006290| 65 6c 65 74 65 28 29 3b | 20 72 65 74 75 72 6e 20 |elete();| return |
|000062a0| 67 6d 31 3b 0a 58 20 20 | 20 20 20 20 7d 0a 58 20 |gm1;.X  |    }.X |
|000062b0| 20 20 20 20 20 65 6c 73 | 65 20 69 66 20 28 67 6d |     els|e if (gm|
|000062c0| 32 2d 3e 54 79 70 65 28 | 29 20 3d 3d 20 6d 74 78 |2->Type(|) == mtx|
|000062d0| 20 26 26 20 67 6d 32 2d | 3e 72 65 75 73 65 28 29 | && gm2-|>reuse()|
|000062e0| 20 29 0a 58 20 20 20 20 | 20 20 7b 0a 58 20 20 20 | ).X    |  {.X   |
|000062f0| 20 20 20 20 20 20 52 45 | 50 4f 52 54 0a 58 20 20 |      RE|PORT.X  |
|00006300| 20 20 20 20 20 20 20 4d | 61 74 72 69 78 52 6f 77 |       M|atrixRow|
|00006310| 20 6d 72 31 28 67 6d 31 | 2c 20 4c 6f 61 64 4f 6e | mr1(gm1|, LoadOn|
|00006320| 45 6e 74 72 79 29 3b 0a | 58 09 20 4d 61 74 72 69 |Entry);.|X. Matri|
|00006330| 78 52 6f 77 20 6d 72 32 | 28 67 6d 32 2c 20 4c 6f |xRow mr2|(gm2, Lo|
|00006340| 61 64 4f 6e 45 6e 74 72 | 79 2b 53 74 6f 72 65 4f |adOnEntr|y+StoreO|
|00006350| 6e 45 78 69 74 2b 44 69 | 72 65 63 74 50 61 72 74 |nExit+Di|rectPart|
|00006360| 29 3b 0a 58 09 20 77 68 | 69 6c 65 20 28 6e 72 2d |);.X. wh|ile (nr-|
|00006370| 2d 29 20 7b 20 6d 72 32 | 2e 52 65 76 53 75 62 28 |-) { mr2|.RevSub(|
|00006380| 6d 72 31 29 3b 20 6d 72 | 31 2e 4e 65 78 74 28 29 |mr1); mr|1.Next()|
|00006390| 3b 20 6d 72 32 2e 4e 65 | 78 74 28 29 3b 20 7d 0a |; mr2.Ne|xt(); }.|
|000063a0| 58 09 20 69 66 20 28 67 | 6d 31 2d 3e 54 79 70 65 |X. if (g|m1->Type|
|000063b0| 28 29 21 3d 6d 74 78 29 | 20 67 6d 31 2d 3e 74 44 |()!=mtx)| gm1->tD|
|000063c0| 65 6c 65 74 65 28 29 3b | 0a 58 09 20 72 65 74 75 |elete();|.X. retu|
|000063d0| 72 6e 20 67 6d 32 3b 0a | 58 20 20 20 20 20 20 7d |rn gm2;.|X      }|
|000063e0| 0a 58 20 20 20 20 20 20 | 65 6c 73 65 0a 58 20 20 |.X      |else.X  |
|000063f0| 20 20 20 20 7b 0a 58 20 | 20 20 20 20 20 20 20 20 |    {.X |        |
+--------+-------------------------+-------------------------+--------+--------+