/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / sources / misc / 4255

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 6d 69 73 63 0a 50 61 |.sources|.misc.Pa|
|00000020| 74 68 3a 20 73 70 61 72 | 6b 79 21 6b 65 6e 74 0a |th: spar|ky!kent.|
|00000030| 46 72 6f 6d 3a 20 72 6f | 62 65 72 74 64 40 6b 61 |From: ro|bertd@ka|
|00000040| 75 72 69 2e 76 75 77 2e | 61 63 2e 6e 7a 20 28 52 |uri.vuw.|ac.nz (R|
|00000050| 6f 62 65 72 74 20 44 61 | 76 69 65 73 29 0a 53 75 |obert Da|vies).Su|
|00000060| 62 6a 65 63 74 3a 20 76 | 33 34 69 31 31 33 3a 20 |bject: v|34i113: |
|00000070| 20 6e 65 77 6d 61 74 30 | 37 20 2d 20 41 20 6d 61 | newmat0|7 - A ma|
|00000080| 74 72 69 78 20 70 61 63 | 6b 61 67 65 20 69 6e 20 |trix pac|kage in |
|00000090| 43 2b 2b 2c 20 50 61 72 | 74 30 37 2f 30 38 0a 4d |C++, Par|t07/08.M|
|000000a0| 65 73 73 61 67 65 2d 49 | 44 3a 20 3c 31 39 39 33 |essage-I|D: <1993|
|000000b0| 4a 61 6e 31 31 2e 31 35 | 33 33 32 38 2e 32 36 31 |Jan11.15|3328.261|
|000000c0| 32 40 73 70 61 72 6b 79 | 2e 69 6d 64 2e 73 74 65 |2@sparky|.imd.ste|
|000000d0| 72 6c 69 6e 67 2e 63 6f | 6d 3e 0a 46 6f 6c 6c 6f |rling.co|m>.Follo|
|000000e0| 77 75 70 2d 54 6f 3a 20 | 63 6f 6d 70 2e 73 6f 75 |wup-To: |comp.sou|
|000000f0| 72 63 65 73 2e 64 0a 58 | 2d 4d 64 34 2d 53 69 67 |rces.d.X|-Md4-Sig|
|00000100| 6e 61 74 75 72 65 3a 20 | 32 38 31 32 30 39 65 38 |nature: |281209e8|
|00000110| 30 62 35 65 37 62 35 37 | 62 38 34 63 34 66 61 62 |0b5e7b57|b84c4fab|
|00000120| 32 63 61 34 39 34 30 36 | 0a 53 65 6e 64 65 72 3a |2ca49406|.Sender:|
|00000130| 20 6b 65 6e 74 40 73 70 | 61 72 6b 79 2e 69 6d 64 | kent@sp|arky.imd|
|00000140| 2e 73 74 65 72 6c 69 6e | 67 2e 63 6f 6d 20 28 4b |.sterlin|g.com (K|
|00000150| 65 6e 74 20 4c 61 6e 64 | 66 69 65 6c 64 29 0a 4f |ent Land|field).O|
|00000160| 72 67 61 6e 69 7a 61 74 | 69 6f 6e 3a 20 53 74 65 |rganizat|ion: Ste|
|00000170| 72 6c 69 6e 67 20 53 6f | 66 74 77 61 72 65 0a 52 |rling So|ftware.R|
|00000180| 65 66 65 72 65 6e 63 65 | 73 3a 20 3c 63 73 6d 2d |eference|s: <csm-|
|00000190| 76 33 34 69 31 30 37 3d | 6e 65 77 6d 61 74 30 37 |v34i107=|newmat07|
|000001a0| 2e 30 39 32 37 30 36 40 | 73 70 61 72 6b 79 2e 49 |.092706@|sparky.I|
|000001b0| 4d 44 2e 53 74 65 72 6c | 69 6e 67 2e 43 4f 4d 3e |MD.Sterl|ing.COM>|
|000001c0| 0a 44 61 74 65 3a 20 4d | 6f 6e 2c 20 31 31 20 4a |.Date: M|on, 11 J|
|000001d0| 61 6e 20 31 39 39 33 20 | 31 35 3a 33 33 3a 32 38 |an 1993 |15:33:28|
|000001e0| 20 47 4d 54 0a 41 70 70 | 72 6f 76 65 64 3a 20 6b | GMT.App|roved: k|
|000001f0| 65 6e 74 40 73 70 61 72 | 6b 79 2e 69 6d 64 2e 73 |ent@spar|ky.imd.s|
|00000200| 74 65 72 6c 69 6e 67 2e | 63 6f 6d 0a 4c 69 6e 65 |terling.|com.Line|
|00000210| 73 3a 20 32 33 30 31 0a | 0a 53 75 62 6d 69 74 74 |s: 2301.|.Submitt|
|00000220| 65 64 2d 62 79 3a 20 72 | 6f 62 65 72 74 64 40 6b |ed-by: r|obertd@k|
|00000230| 61 75 72 69 2e 76 75 77 | 2e 61 63 2e 6e 7a 20 28 |auri.vuw|.ac.nz (|
|00000240| 52 6f 62 65 72 74 20 44 | 61 76 69 65 73 29 0a 50 |Robert D|avies).P|
|00000250| 6f 73 74 69 6e 67 2d 6e | 75 6d 62 65 72 3a 20 56 |osting-n|umber: V|
|00000260| 6f 6c 75 6d 65 20 33 34 | 2c 20 49 73 73 75 65 20 |olume 34|, Issue |
|00000270| 31 31 33 0a 41 72 63 68 | 69 76 65 2d 6e 61 6d 65 |113.Arch|ive-name|
|00000280| 3a 20 6e 65 77 6d 61 74 | 30 37 2f 70 61 72 74 30 |: newmat|07/part0|
|00000290| 37 0a 45 6e 76 69 72 6f | 6e 6d 65 6e 74 3a 20 43 |7.Enviro|nment: C|
|000002a0| 2b 2b 0a 53 75 70 65 72 | 73 65 64 65 73 3a 20 6e |++.Super|sedes: n|
|000002b0| 65 77 6d 61 74 30 36 3a | 20 56 6f 6c 75 6d 65 20 |ewmat06:| Volume |
|000002c0| 33 34 2c 20 49 73 73 75 | 65 20 37 2d 31 33 0a 0a |34, Issu|e 7-13..|
|000002d0| 23 21 20 2f 62 69 6e 2f | 73 68 0a 23 20 54 68 69 |#! /bin/|sh.# Thi|
|000002e0| 73 20 69 73 20 61 20 73 | 68 65 6c 6c 20 61 72 63 |s is a s|hell arc|
|000002f0| 68 69 76 65 2e 20 20 52 | 65 6d 6f 76 65 20 61 6e |hive.  R|emove an|
|00000300| 79 74 68 69 6e 67 20 62 | 65 66 6f 72 65 20 74 68 |ything b|efore th|
|00000310| 69 73 20 6c 69 6e 65 2c | 20 74 68 65 6e 20 75 6e |is line,| then un|
|00000320| 70 61 63 6b 0a 23 20 69 | 74 20 62 79 20 73 61 76 |pack.# i|t by sav|
|00000330| 69 6e 67 20 69 74 20 69 | 6e 74 6f 20 61 20 66 69 |ing it i|nto a fi|
|00000340| 6c 65 20 61 6e 64 20 74 | 79 70 69 6e 67 20 22 73 |le and t|yping "s|
|00000350| 68 20 66 69 6c 65 22 2e | 20 20 54 6f 20 6f 76 65 |h file".|  To ove|
|00000360| 72 77 72 69 74 65 20 65 | 78 69 73 74 69 6e 67 0a |rwrite e|xisting.|
|00000370| 23 20 66 69 6c 65 73 2c | 20 74 79 70 65 20 22 73 |# files,| type "s|
|00000380| 68 20 66 69 6c 65 20 2d | 63 22 2e 20 20 59 6f 75 |h file -|c".  You|
|00000390| 20 63 61 6e 20 61 6c 73 | 6f 20 66 65 65 64 20 74 | can als|o feed t|
|000003a0| 68 69 73 20 61 73 20 73 | 74 61 6e 64 61 72 64 20 |his as s|tandard |
|000003b0| 69 6e 70 75 74 20 76 69 | 61 0a 23 20 75 6e 73 68 |input vi|a.# unsh|
|000003c0| 61 72 2c 20 6f 72 20 62 | 79 20 74 79 70 69 6e 67 |ar, or b|y typing|
|000003d0| 20 22 73 68 20 3c 66 69 | 6c 65 22 2c 20 65 2e 67 | "sh <fi|le", e.g|
|000003e0| 2e 2e 20 20 49 66 20 74 | 68 69 73 20 61 72 63 68 |..  If t|his arch|
|000003f0| 69 76 65 20 69 73 20 63 | 6f 6d 70 6c 65 74 65 2c |ive is c|omplete,|
|00000400| 20 79 6f 75 0a 23 20 77 | 69 6c 6c 20 73 65 65 20 | you.# w|ill see |
|00000410| 74 68 65 20 66 6f 6c 6c | 6f 77 69 6e 67 20 6d 65 |the foll|owing me|
|00000420| 73 73 61 67 65 20 61 74 | 20 74 68 65 20 65 6e 64 |ssage at| the end|
|00000430| 3a 0a 23 09 09 22 45 6e | 64 20 6f 66 20 61 72 63 |:.#.."En|d of arc|
|00000440| 68 69 76 65 20 37 20 28 | 6f 66 20 38 29 2e 22 0a |hive 7 (|of 8).".|
|00000450| 23 20 43 6f 6e 74 65 6e | 74 73 3a 20 20 74 6d 74 |# Conten|ts:  tmt|
|00000460| 2e 63 78 78 20 74 6d 74 | 31 2e 63 78 78 20 74 6d |.cxx tmt|1.cxx tm|
|00000470| 74 32 2e 63 78 78 20 74 | 6d 74 33 2e 63 78 78 20 |t2.cxx t|mt3.cxx |
|00000480| 74 6d 74 34 2e 63 78 78 | 20 74 6d 74 35 2e 63 78 |tmt4.cxx| tmt5.cx|
|00000490| 78 0a 23 20 20 20 74 6d | 74 36 2e 63 78 78 20 74 |x.#   tm|t6.cxx t|
|000004a0| 6d 74 37 2e 63 78 78 20 | 74 6d 74 38 2e 63 78 78 |mt7.cxx |tmt8.cxx|
|000004b0| 20 74 6d 74 39 2e 63 78 | 78 20 74 6d 74 61 2e 63 | tmt9.cx|x tmta.c|
|000004c0| 78 78 20 74 6d 74 62 2e | 63 78 78 20 74 6d 74 63 |xx tmtb.|cxx tmtc|
|000004d0| 2e 63 78 78 0a 23 20 20 | 20 74 6d 74 64 2e 63 78 |.cxx.#  | tmtd.cx|
|000004e0| 78 20 74 6d 74 65 2e 63 | 78 78 20 74 6d 74 66 2e |x tmte.c|xx tmtf.|
|000004f0| 63 78 78 20 74 6d 74 67 | 2e 63 78 78 20 74 6d 74 |cxx tmtg|.cxx tmt|
|00000500| 68 2e 63 78 78 20 74 6d | 74 69 2e 63 78 78 0a 23 |h.cxx tm|ti.cxx.#|
|00000510| 20 57 72 61 70 70 65 64 | 20 62 79 20 72 6f 62 65 | Wrapped| by robe|
|00000520| 72 74 40 6b 65 61 20 6f | 6e 20 53 75 6e 20 4a 61 |rt@kea o|n Sun Ja|
|00000530| 6e 20 31 30 20 32 33 3a | 35 38 3a 33 30 20 31 39 |n 10 23:|58:30 19|
|00000540| 39 33 0a 50 41 54 48 3d | 2f 62 69 6e 3a 2f 75 73 |93.PATH=|/bin:/us|
|00000550| 72 2f 62 69 6e 3a 2f 75 | 73 72 2f 75 63 62 20 3b |r/bin:/u|sr/ucb ;|
|00000560| 20 65 78 70 6f 72 74 20 | 50 41 54 48 0a 69 66 20 | export |PATH.if |
|00000570| 74 65 73 74 20 2d 66 20 | 27 74 6d 74 2e 63 78 78 |test -f |'tmt.cxx|
|00000580| 27 20 2d 61 20 22 24 7b | 31 7d 22 20 21 3d 20 22 |' -a "${|1}" != "|
|00000590| 2d 63 22 20 3b 20 74 68 | 65 6e 20 0a 20 20 65 63 |-c" ; th|en .  ec|
|000005a0| 68 6f 20 73 68 61 72 3a | 20 57 69 6c 6c 20 6e 6f |ho shar:| Will no|
|000005b0| 74 20 63 6c 6f 62 62 65 | 72 20 65 78 69 73 74 69 |t clobbe|r existi|
|000005c0| 6e 67 20 66 69 6c 65 20 | 5c 22 27 74 6d 74 2e 63 |ng file |\"'tmt.c|
|000005d0| 78 78 27 5c 22 0a 65 6c | 73 65 0a 65 63 68 6f 20 |xx'\".el|se.echo |
|000005e0| 73 68 61 72 3a 20 45 78 | 74 72 61 63 74 69 6e 67 |shar: Ex|tracting|
|000005f0| 20 5c 22 27 74 6d 74 2e | 63 78 78 27 5c 22 20 5c | \"'tmt.|cxx'\" \|
|00000600| 28 34 36 33 33 20 63 68 | 61 72 61 63 74 65 72 73 |(4633 ch|aracters|
|00000610| 5c 29 0a 73 65 64 20 22 | 73 2f 5e 58 2f 2f 22 20 |\).sed "|s/^X//" |
|00000620| 3e 27 74 6d 74 2e 63 78 | 78 27 20 3c 3c 27 45 4e |>'tmt.cx|x' <<'EN|
|00000630| 44 5f 4f 46 5f 46 49 4c | 45 27 0a 58 23 64 65 66 |D_OF_FIL|E'.X#def|
|00000640| 69 6e 65 20 57 41 4e 54 | 5f 53 54 52 45 41 4d 0a |ine WANT|_STREAM.|
|00000650| 58 0a 58 23 69 6e 63 6c | 75 64 65 20 22 69 6e 63 |X.X#incl|ude "inc|
|00000660| 6c 75 64 65 2e 68 22 0a | 58 0a 58 23 69 6e 63 6c |lude.h".|X.X#incl|
|00000670| 75 64 65 20 22 6e 65 77 | 6d 61 74 2e 68 22 0a 58 |ude "new|mat.h".X|
|00000680| 0a 58 2f 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X/*****|********|
|00000690| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 20 |********|******* |
|000006a0| 74 65 73 74 20 70 72 6f | 67 72 61 6d 20 2a 2a 2a |test pro|gram ***|
|000006b0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000006c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2f 0a 58 0a 58 |********|***/.X.X|
|000006d0| 63 6c 61 73 73 20 50 72 | 69 6e 74 43 6f 75 6e 74 |class Pr|intCount|
|000006e0| 65 72 0a 58 7b 0a 58 20 | 20 20 69 6e 74 20 63 6f |er.X{.X |  int co|
|000006f0| 75 6e 74 3b 0a 58 20 20 | 20 63 68 61 72 2a 20 73 |unt;.X  | char* s|
|00000700| 3b 0a 58 70 75 62 6c 69 | 63 3a 0a 58 20 20 20 7e |;.Xpubli|c:.X   ~|
|00000710| 50 72 69 6e 74 43 6f 75 | 6e 74 65 72 28 29 3b 0a |PrintCou|nter();.|
|00000720| 58 20 20 20 50 72 69 6e | 74 43 6f 75 6e 74 65 72 |X   Prin|tCounter|
|00000730| 28 63 68 61 72 20 2a 20 | 73 78 29 20 3a 20 63 6f |(char * |sx) : co|
|00000740| 75 6e 74 28 30 29 2c 20 | 73 28 73 78 29 20 7b 7d |unt(0), |s(sx) {}|
|00000750| 0a 58 20 20 20 76 6f 69 | 64 20 6f 70 65 72 61 74 |.X   voi|d operat|
|00000760| 6f 72 2b 2b 28 29 20 7b | 20 63 6f 75 6e 74 2b 2b |or++() {| count++|
|00000770| 3b 20 7d 0a 58 7d 3b 0a | 58 0a 58 20 20 20 50 72 |; }.X};.|X.X   Pr|
|00000780| 69 6e 74 43 6f 75 6e 74 | 65 72 20 50 43 5a 28 22 |intCount|er PCZ("|
|00000790| 4e 75 6d 62 65 72 20 6f | 66 20 6e 6f 6e 2d 7a 65 |Number o|f non-ze|
|000007a0| 72 6f 20 6d 61 74 72 69 | 63 65 73 20 3d 20 22 29 |ro matri|ces = ")|
|000007b0| 3b 0a 58 20 20 20 50 72 | 69 6e 74 43 6f 75 6e 74 |;.X   Pr|intCount|
|000007c0| 65 72 20 50 43 4e 28 22 | 4e 75 6d 62 65 72 20 6f |er PCN("|Number o|
|000007d0| 66 20 6d 61 74 72 69 63 | 65 73 20 74 65 73 74 65 |f matric|es teste|
|000007e0| 64 20 20 20 3d 20 22 29 | 3b 0a 58 0a 58 50 72 69 |d   = ")|;.X.XPri|
|000007f0| 6e 74 43 6f 75 6e 74 65 | 72 3a 3a 7e 50 72 69 6e |ntCounte|r::~Prin|
|00000800| 74 43 6f 75 6e 74 65 72 | 28 29 0a 58 7b 20 63 6f |tCounter|().X{ co|
|00000810| 75 74 20 3c 3c 20 73 20 | 3c 3c 20 63 6f 75 6e 74 |ut << s |<< count|
|00000820| 20 3c 3c 20 22 5c 6e 22 | 3b 20 7d 0a 58 0a 58 0a | << "\n"|; }.X.X.|
|00000830| 58 76 6f 69 64 20 50 72 | 69 6e 74 28 63 6f 6e 73 |Xvoid Pr|int(cons|
|00000840| 74 20 4d 61 74 72 69 78 | 26 20 58 29 0a 58 7b 0a |t Matrix|& X).X{.|
|00000850| 58 20 20 20 2b 2b 50 43 | 4e 3b 0a 58 20 20 20 63 |X   ++PC|N;.X   c|
|00000860| 6f 75 74 20 3c 3c 20 22 | 5c 6e 4d 61 74 72 69 78 |out << "|\nMatrix|
|00000870| 20 74 79 70 65 3a 20 22 | 20 3c 3c 20 28 63 68 61 | type: "| << (cha|
|00000880| 72 2a 29 58 2e 54 79 70 | 65 28 29 20 3c 3c 20 22 |r*)X.Typ|e() << "|
|00000890| 20 28 22 3b 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c | (";.X  | cout <<|
|000008a0| 20 58 2e 4e 72 6f 77 73 | 28 29 20 3c 3c 20 22 2c | X.Nrows|() << ",|
|000008b0| 20 22 3b 0a 58 20 20 20 | 63 6f 75 74 20 3c 3c 20 | ";.X   |cout << |
|000008c0| 58 2e 4e 63 6f 6c 73 28 | 29 20 3c 3c 20 22 29 5c |X.Ncols(|) << ")\|
|000008d0| 6e 5c 6e 22 3b 0a 58 20 | 20 20 69 66 20 28 58 2e |n\n";.X |  if (X.|
|000008e0| 49 73 5a 65 72 6f 28 29 | 29 20 7b 20 63 6f 75 74 |IsZero()|) { cout|
|000008f0| 20 3c 3c 20 22 41 6c 6c | 20 65 6c 65 6d 65 6e 74 | << "All| element|
|00000900| 73 20 61 72 65 20 7a 65 | 72 6f 5c 6e 22 20 3c 3c |s are ze|ro\n" <<|
|00000910| 20 66 6c 75 73 68 3b 20 | 72 65 74 75 72 6e 3b 20 | flush; |return; |
|00000920| 7d 0a 58 20 20 20 69 6e | 74 20 6e 72 3d 58 2e 4e |}.X   in|t nr=X.N|
|00000930| 72 6f 77 73 28 29 3b 20 | 69 6e 74 20 6e 63 3d 58 |rows(); |int nc=X|
|00000940| 2e 4e 63 6f 6c 73 28 29 | 3b 0a 58 20 20 20 66 6f |.Ncols()|;.X   fo|
|00000950| 72 20 28 69 6e 74 20 69 | 3d 31 3b 20 69 3c 3d 6e |r (int i|=1; i<=n|
|00000960| 72 3b 20 69 2b 2b 29 0a | 58 20 20 20 7b 0a 58 20 |r; i++).|X   {.X |
|00000970| 20 20 20 20 20 66 6f 72 | 20 28 69 6e 74 20 6a 3d |     for| (int j=|
|00000980| 31 3b 20 6a 3c 3d 6e 63 | 3b 20 6a 2b 2b 29 20 20 |1; j<=nc|; j++)  |
|00000990| 63 6f 75 74 20 3c 3c 20 | 58 28 69 2c 6a 29 20 3c |cout << |X(i,j) <|
|000009a0| 3c 20 22 5c 74 22 3b 0a | 58 20 20 20 20 20 20 63 |< "\t";.|X      c|
|000009b0| 6f 75 74 20 3c 3c 20 22 | 5c 6e 22 3b 0a 58 20 20 |out << "|\n";.X  |
|000009c0| 20 7d 0a 58 20 20 20 63 | 6f 75 74 20 3c 3c 20 66 | }.X   c|out << f|
|000009d0| 6c 75 73 68 3b 20 2b 2b | 50 43 5a 3b 0a 58 7d 0a |lush; ++|PCZ;.X}.|
|000009e0| 58 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |X.Xvoid |Print(co|
|000009f0| 6e 73 74 20 55 70 70 65 | 72 54 72 69 61 6e 67 75 |nst Uppe|rTriangu|
|00000a00| 6c 61 72 4d 61 74 72 69 | 78 26 20 58 29 0a 58 7b |larMatri|x& X).X{|
|00000a10| 0a 58 20 20 20 2b 2b 50 | 43 4e 3b 0a 58 20 20 20 |.X   ++P|CN;.X   |
|00000a20| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 4d 61 74 72 69 |cout << |"\nMatri|
|00000a30| 78 20 74 79 70 65 3a 20 | 22 20 3c 3c 20 28 63 68 |x type: |" << (ch|
|00000a40| 61 72 2a 29 58 2e 54 79 | 70 65 28 29 20 3c 3c 20 |ar*)X.Ty|pe() << |
|00000a50| 22 20 28 22 3b 0a 58 20 | 20 20 63 6f 75 74 20 3c |" (";.X |  cout <|
|00000a60| 3c 20 58 2e 4e 72 6f 77 | 73 28 29 20 3c 3c 20 22 |< X.Nrow|s() << "|
|00000a70| 2c 20 22 3b 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c |, ";.X  | cout <<|
|00000a80| 20 58 2e 4e 63 6f 6c 73 | 28 29 20 3c 3c 20 22 29 | X.Ncols|() << ")|
|00000a90| 5c 6e 5c 6e 22 3b 0a 58 | 20 20 20 69 66 20 28 58 |\n\n";.X|   if (X|
|00000aa0| 2e 49 73 5a 65 72 6f 28 | 29 29 20 7b 20 63 6f 75 |.IsZero(|)) { cou|
|00000ab0| 74 20 3c 3c 20 22 41 6c | 6c 20 65 6c 65 6d 65 6e |t << "Al|l elemen|
|00000ac0| 74 73 20 61 72 65 20 7a | 65 72 6f 5c 6e 22 20 3c |ts are z|ero\n" <|
|00000ad0| 3c 20 66 6c 75 73 68 3b | 20 72 65 74 75 72 6e 3b |< flush;| return;|
|00000ae0| 20 7d 0a 58 20 20 20 69 | 6e 74 20 6e 72 3d 58 2e | }.X   i|nt nr=X.|
|00000af0| 4e 72 6f 77 73 28 29 3b | 20 69 6e 74 20 6e 63 3d |Nrows();| int nc=|
|00000b00| 58 2e 4e 63 6f 6c 73 28 | 29 3b 0a 58 20 20 20 66 |X.Ncols(|);.X   f|
|00000b10| 6f 72 20 28 69 6e 74 20 | 69 3d 31 3b 20 69 3c 3d |or (int |i=1; i<=|
|00000b20| 6e 72 3b 20 69 2b 2b 29 | 0a 58 20 20 20 7b 0a 58 |nr; i++)|.X   {.X|
|00000b30| 20 20 20 20 20 20 66 6f | 72 20 28 69 6e 74 20 6a |      fo|r (int j|
|00000b40| 3d 31 3b 20 6a 3c 69 3b | 20 6a 2b 2b 29 20 63 6f |=1; j<i;| j++) co|
|00000b50| 75 74 20 3c 3c 20 22 5c | 74 22 3b 0a 58 20 20 20 |ut << "\|t";.X   |
|00000b60| 20 20 20 66 6f 72 20 28 | 6a 3d 69 3b 20 6a 3c 3d |   for (|j=i; j<=|
|00000b70| 6e 63 3b 20 6a 2b 2b 29 | 20 20 63 6f 75 74 20 3c |nc; j++)|  cout <|
|00000b80| 3c 20 58 28 69 2c 6a 29 | 20 3c 3c 20 22 5c 74 22 |< X(i,j)| << "\t"|
|00000b90| 3b 0a 58 20 20 20 20 20 | 20 63 6f 75 74 20 3c 3c |;.X     | cout <<|
|00000ba0| 20 22 5c 6e 22 3b 0a 58 | 20 20 20 7d 0a 58 20 20 | "\n";.X|   }.X  |
|00000bb0| 20 63 6f 75 74 20 3c 3c | 20 66 6c 75 73 68 3b 20 | cout <<| flush; |
|00000bc0| 2b 2b 50 43 5a 3b 0a 58 | 7d 0a 58 0a 58 76 6f 69 |++PCZ;.X|}.X.Xvoi|
|00000bd0| 64 20 50 72 69 6e 74 28 | 63 6f 6e 73 74 20 44 69 |d Print(|const Di|
|00000be0| 61 67 6f 6e 61 6c 4d 61 | 74 72 69 78 26 20 58 29 |agonalMa|trix& X)|
|00000bf0| 0a 58 7b 0a 58 20 20 20 | 2b 2b 50 43 4e 3b 0a 58 |.X{.X   |++PCN;.X|
|00000c00| 20 20 20 63 6f 75 74 20 | 3c 3c 20 22 5c 6e 4d 61 |   cout |<< "\nMa|
|00000c10| 74 72 69 78 20 74 79 70 | 65 3a 20 22 20 3c 3c 20 |trix typ|e: " << |
|00000c20| 28 63 68 61 72 2a 29 58 | 2e 54 79 70 65 28 29 20 |(char*)X|.Type() |
|00000c30| 3c 3c 20 22 20 28 22 3b | 0a 58 20 20 20 63 6f 75 |<< " (";|.X   cou|
|00000c40| 74 20 3c 3c 20 58 2e 4e | 72 6f 77 73 28 29 20 3c |t << X.N|rows() <|
|00000c50| 3c 20 22 2c 20 22 3b 0a | 58 20 20 20 63 6f 75 74 |< ", ";.|X   cout|
|00000c60| 20 3c 3c 20 58 2e 4e 63 | 6f 6c 73 28 29 20 3c 3c | << X.Nc|ols() <<|
|00000c70| 20 22 29 5c 6e 5c 6e 22 | 3b 0a 58 20 20 20 69 66 | ")\n\n"|;.X   if|
|00000c80| 20 28 58 2e 49 73 5a 65 | 72 6f 28 29 29 20 7b 20 | (X.IsZe|ro()) { |
|00000c90| 63 6f 75 74 20 3c 3c 20 | 22 41 6c 6c 20 65 6c 65 |cout << |"All ele|
|00000ca0| 6d 65 6e 74 73 20 61 72 | 65 20 7a 65 72 6f 5c 6e |ments ar|e zero\n|
|00000cb0| 22 20 3c 3c 20 66 6c 75 | 73 68 3b 20 72 65 74 75 |" << flu|sh; retu|
|00000cc0| 72 6e 3b 20 7d 0a 58 20 | 20 20 69 6e 74 20 6e 72 |rn; }.X |  int nr|
|00000cd0| 3d 58 2e 4e 72 6f 77 73 | 28 29 3b 20 69 6e 74 20 |=X.Nrows|(); int |
|00000ce0| 6e 63 3d 58 2e 4e 63 6f | 6c 73 28 29 3b 0a 58 20 |nc=X.Nco|ls();.X |
|00000cf0| 20 20 66 6f 72 20 28 69 | 6e 74 20 69 3d 31 3b 20 |  for (i|nt i=1; |
|00000d00| 69 3c 3d 6e 72 3b 20 69 | 2b 2b 29 0a 58 20 20 20 |i<=nr; i|++).X   |
|00000d10| 7b 0a 58 20 20 20 20 20 | 20 66 6f 72 20 28 69 6e |{.X     | for (in|
|00000d20| 74 20 6a 3d 31 3b 20 6a | 3c 69 3b 20 6a 2b 2b 29 |t j=1; j|<i; j++)|
|00000d30| 20 63 6f 75 74 20 3c 3c | 20 22 5c 74 22 3b 0a 58 | cout <<| "\t";.X|
|00000d40| 20 20 20 20 20 20 69 66 | 20 28 69 3c 3d 6e 63 29 |      if| (i<=nc)|
|00000d50| 20 63 6f 75 74 20 3c 3c | 20 58 28 69 2c 69 29 20 | cout <<| X(i,i) |
|00000d60| 3c 3c 20 22 5c 74 22 3b | 0a 58 20 20 20 20 20 20 |<< "\t";|.X      |
|00000d70| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 22 3b 0a 58 20 |cout << |"\n";.X |
|00000d80| 20 20 7d 0a 58 20 20 20 | 63 6f 75 74 20 3c 3c 20 |  }.X   |cout << |
|00000d90| 66 6c 75 73 68 3b 20 2b | 2b 50 43 5a 3b 0a 58 7d |flush; +|+PCZ;.X}|
|00000da0| 0a 58 0a 58 76 6f 69 64 | 20 50 72 69 6e 74 28 63 |.X.Xvoid| Print(c|
|00000db0| 6f 6e 73 74 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |onst Sym|metricMa|
|00000dc0| 74 72 69 78 26 20 58 29 | 0a 58 7b 0a 58 20 20 20 |trix& X)|.X{.X   |
|00000dd0| 2b 2b 50 43 4e 3b 0a 58 | 20 20 20 63 6f 75 74 20 |++PCN;.X|   cout |
|00000de0| 3c 3c 20 22 5c 6e 4d 61 | 74 72 69 78 20 74 79 70 |<< "\nMa|trix typ|
|00000df0| 65 3a 20 22 20 3c 3c 20 | 28 63 68 61 72 2a 29 58 |e: " << |(char*)X|
|00000e00| 2e 54 79 70 65 28 29 20 | 3c 3c 20 22 20 28 22 3b |.Type() |<< " (";|
|00000e10| 0a 58 20 20 20 63 6f 75 | 74 20 3c 3c 20 58 2e 4e |.X   cou|t << X.N|
|00000e20| 72 6f 77 73 28 29 20 3c | 3c 20 22 2c 20 22 3b 0a |rows() <|< ", ";.|
|00000e30| 58 20 20 20 63 6f 75 74 | 20 3c 3c 20 58 2e 4e 63 |X   cout| << X.Nc|
|00000e40| 6f 6c 73 28 29 20 3c 3c | 20 22 29 5c 6e 5c 6e 22 |ols() <<| ")\n\n"|
|00000e50| 3b 0a 58 20 20 20 69 66 | 20 28 58 2e 49 73 5a 65 |;.X   if| (X.IsZe|
|00000e60| 72 6f 28 29 29 20 7b 20 | 63 6f 75 74 20 3c 3c 20 |ro()) { |cout << |
|00000e70| 22 41 6c 6c 20 65 6c 65 | 6d 65 6e 74 73 20 61 72 |"All ele|ments ar|
|00000e80| 65 20 7a 65 72 6f 5c 6e | 22 20 3c 3c 20 66 6c 75 |e zero\n|" << flu|
|00000e90| 73 68 3b 20 72 65 74 75 | 72 6e 3b 20 7d 0a 58 20 |sh; retu|rn; }.X |
|00000ea0| 20 20 69 6e 74 20 6e 72 | 3d 58 2e 4e 72 6f 77 73 |  int nr|=X.Nrows|
|00000eb0| 28 29 3b 20 69 6e 74 20 | 6e 63 3d 58 2e 4e 63 6f |(); int |nc=X.Nco|
|00000ec0| 6c 73 28 29 3b 0a 58 20 | 20 20 66 6f 72 20 28 69 |ls();.X |  for (i|
|00000ed0| 6e 74 20 69 3d 31 3b 20 | 69 3c 3d 6e 72 3b 20 69 |nt i=1; |i<=nr; i|
|00000ee0| 2b 2b 29 0a 58 20 20 20 | 7b 0a 58 20 20 20 20 20 |++).X   |{.X     |
|00000ef0| 20 66 6f 72 20 28 69 6e | 74 20 6a 3d 31 3b 20 6a | for (in|t j=1; j|
|00000f00| 3c 69 3b 20 6a 2b 2b 29 | 20 63 6f 75 74 20 3c 3c |<i; j++)| cout <<|
|00000f10| 20 58 28 6a 2c 69 29 20 | 3c 3c 20 22 5c 74 22 3b | X(j,i) |<< "\t";|
|00000f20| 0a 58 20 20 20 20 20 20 | 66 6f 72 20 28 6a 3d 69 |.X      |for (j=i|
|00000f30| 3b 20 6a 3c 3d 6e 63 3b | 20 6a 2b 2b 29 20 20 63 |; j<=nc;| j++)  c|
|00000f40| 6f 75 74 20 3c 3c 20 58 | 28 69 2c 6a 29 20 3c 3c |out << X|(i,j) <<|
|00000f50| 20 22 5c 74 22 3b 0a 58 | 20 20 20 20 20 20 63 6f | "\t";.X|      co|
|00000f60| 75 74 20 3c 3c 20 22 5c | 6e 22 3b 0a 58 20 20 20 |ut << "\|n";.X   |
|00000f70| 7d 0a 58 20 20 20 63 6f | 75 74 20 3c 3c 20 66 6c |}.X   co|ut << fl|
|00000f80| 75 73 68 3b 20 2b 2b 50 | 43 5a 3b 0a 58 7d 0a 58 |ush; ++P|CZ;.X}.X|
|00000f90| 0a 58 76 6f 69 64 20 50 | 72 69 6e 74 28 63 6f 6e |.Xvoid P|rint(con|
|00000fa0| 73 74 20 4c 6f 77 65 72 | 54 72 69 61 6e 67 75 6c |st Lower|Triangul|
|00000fb0| 61 72 4d 61 74 72 69 78 | 26 20 58 29 0a 58 7b 0a |arMatrix|& X).X{.|
|00000fc0| 58 20 20 20 2b 2b 50 43 | 4e 3b 0a 58 20 20 20 63 |X   ++PC|N;.X   c|
|00000fd0| 6f 75 74 20 3c 3c 20 22 | 5c 6e 4d 61 74 72 69 78 |out << "|\nMatrix|
|00000fe0| 20 74 79 70 65 3a 20 22 | 20 3c 3c 20 28 63 68 61 | type: "| << (cha|
|00000ff0| 72 2a 29 58 2e 54 79 70 | 65 28 29 20 3c 3c 20 22 |r*)X.Typ|e() << "|
|00001000| 20 28 22 3b 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c | (";.X  | cout <<|
|00001010| 20 58 2e 4e 72 6f 77 73 | 28 29 20 3c 3c 20 22 2c | X.Nrows|() << ",|
|00001020| 20 22 3b 0a 58 20 20 20 | 63 6f 75 74 20 3c 3c 20 | ";.X   |cout << |
|00001030| 58 2e 4e 63 6f 6c 73 28 | 29 20 3c 3c 20 22 29 5c |X.Ncols(|) << ")\|
|00001040| 6e 5c 6e 22 3b 0a 58 20 | 20 20 69 66 20 28 58 2e |n\n";.X |  if (X.|
|00001050| 49 73 5a 65 72 6f 28 29 | 29 20 7b 20 63 6f 75 74 |IsZero()|) { cout|
|00001060| 20 3c 3c 20 22 41 6c 6c | 20 65 6c 65 6d 65 6e 74 | << "All| element|
|00001070| 73 20 61 72 65 20 7a 65 | 72 6f 5c 6e 22 20 3c 3c |s are ze|ro\n" <<|
|00001080| 20 66 6c 75 73 68 3b 20 | 72 65 74 75 72 6e 3b 20 | flush; |return; |
|00001090| 7d 0a 58 20 20 20 69 6e | 74 20 6e 72 3d 58 2e 4e |}.X   in|t nr=X.N|
|000010a0| 72 6f 77 73 28 29 3b 0a | 58 20 20 20 66 6f 72 20 |rows();.|X   for |
|000010b0| 28 69 6e 74 20 69 3d 31 | 3b 20 69 3c 3d 6e 72 3b |(int i=1|; i<=nr;|
|000010c0| 20 69 2b 2b 29 0a 58 20 | 20 20 7b 0a 58 20 20 20 | i++).X |  {.X   |
|000010d0| 20 20 20 66 6f 72 20 28 | 69 6e 74 20 6a 3d 31 3b |   for (|int j=1;|
|000010e0| 20 6a 3c 3d 69 3b 20 6a | 2b 2b 29 20 63 6f 75 74 | j<=i; j|++) cout|
|000010f0| 20 3c 3c 20 58 28 69 2c | 6a 29 20 3c 3c 20 22 5c | << X(i,|j) << "\|
|00001100| 74 22 3b 0a 58 20 20 20 | 20 20 20 63 6f 75 74 20 |t";.X   |   cout |
|00001110| 3c 3c 20 22 5c 6e 22 3b | 0a 58 20 20 20 7d 0a 58 |<< "\n";|.X   }.X|
|00001120| 20 20 20 63 6f 75 74 20 | 3c 3c 20 66 6c 75 73 68 |   cout |<< flush|
|00001130| 3b 20 2b 2b 50 43 5a 3b | 0a 58 7d 0a 58 0a 58 0a |; ++PCZ;|.X}.X.X.|
|00001140| 58 76 6f 69 64 20 43 6c | 65 61 6e 28 4d 61 74 72 |Xvoid Cl|ean(Matr|
|00001150| 69 78 26 20 41 2c 20 52 | 65 61 6c 20 63 29 0a 58 |ix& A, R|eal c).X|
|00001160| 7b 0a 58 20 20 20 69 6e | 74 20 6e 72 20 3d 20 41 |{.X   in|t nr = A|
|00001170| 2e 4e 72 6f 77 73 28 29 | 3b 20 69 6e 74 20 6e 63 |.Nrows()|; int nc|
|00001180| 20 3d 20 41 2e 4e 63 6f | 6c 73 28 29 3b 0a 58 20 | = A.Nco|ls();.X |
|00001190| 20 20 66 6f 72 20 28 69 | 6e 74 20 69 3d 31 3b 20 |  for (i|nt i=1; |
|000011a0| 69 3c 3d 6e 72 3b 20 69 | 2b 2b 29 0a 58 20 20 20 |i<=nr; i|++).X   |
|000011b0| 7b 0a 58 20 20 20 20 20 | 20 66 6f 72 20 28 20 69 |{.X     | for ( i|
|000011c0| 6e 74 20 6a 3d 31 3b 20 | 6a 3c 3d 6e 63 3b 20 6a |nt j=1; |j<=nc; j|
|000011d0| 2b 2b 29 0a 58 20 20 20 | 20 20 20 7b 20 52 65 61 |++).X   |   { Rea|
|000011e0| 6c 20 61 20 3d 20 41 28 | 69 2c 6a 29 3b 20 69 66 |l a = A(|i,j); if|
|000011f0| 20 28 28 61 20 3c 20 63 | 29 20 26 26 20 28 61 20 | ((a < c|) && (a |
|00001200| 3e 20 2d 63 29 29 20 41 | 28 69 2c 6a 29 20 3d 20 |> -c)) A|(i,j) = |
|00001210| 30 2e 30 3b 20 7d 0a 58 | 20 20 20 7d 0a 58 7d 0a |0.0; }.X|   }.X}.|
|00001220| 58 0a 58 76 6f 69 64 20 | 43 6c 65 61 6e 28 44 69 |X.Xvoid |Clean(Di|
|00001230| 61 67 6f 6e 61 6c 4d 61 | 74 72 69 78 26 20 41 2c |agonalMa|trix& A,|
|00001240| 20 52 65 61 6c 20 63 29 | 0a 58 7b 0a 58 20 20 20 | Real c)|.X{.X   |
|00001250| 69 6e 74 20 6e 72 20 3d | 20 41 2e 4e 72 6f 77 73 |int nr =| A.Nrows|
|00001260| 28 29 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 6e 74 |();.X   |for (int|
|00001270| 20 69 3d 31 3b 20 69 3c | 3d 6e 72 3b 20 69 2b 2b | i=1; i<|=nr; i++|
|00001280| 29 0a 58 20 20 20 7b 20 | 52 65 61 6c 20 61 20 3d |).X   { |Real a =|
|00001290| 20 41 28 69 2c 69 29 3b | 20 69 66 20 28 28 61 20 | A(i,i);| if ((a |
|000012a0| 3c 20 63 29 20 26 26 20 | 28 61 20 3e 20 2d 63 29 |< c) && |(a > -c)|
|000012b0| 29 20 41 28 69 2c 69 29 | 20 3d 20 30 2e 30 3b 20 |) A(i,i)| = 0.0; |
|000012c0| 7d 0a 58 7d 0a 58 0a 58 | 0a 58 76 6f 69 64 20 74 |}.X}.X.X|.Xvoid t|
|000012d0| 72 79 6d 61 74 31 28 29 | 3b 20 76 6f 69 64 20 74 |rymat1()|; void t|
|000012e0| 72 79 6d 61 74 32 28 29 | 3b 20 76 6f 69 64 20 74 |rymat2()|; void t|
|000012f0| 72 79 6d 61 74 33 28 29 | 3b 0a 58 76 6f 69 64 20 |rymat3()|;.Xvoid |
|00001300| 74 72 79 6d 61 74 34 28 | 29 3b 20 76 6f 69 64 20 |trymat4(|); void |
|00001310| 74 72 79 6d 61 74 35 28 | 29 3b 20 76 6f 69 64 20 |trymat5(|); void |
|00001320| 74 72 79 6d 61 74 36 28 | 29 3b 0a 58 76 6f 69 64 |trymat6(|);.Xvoid|
|00001330| 20 74 72 79 6d 61 74 37 | 28 29 3b 20 76 6f 69 64 | trymat7|(); void|
|00001340| 20 74 72 79 6d 61 74 38 | 28 29 3b 20 76 6f 69 64 | trymat8|(); void|
|00001350| 20 74 72 79 6d 61 74 39 | 28 29 3b 0a 58 76 6f 69 | trymat9|();.Xvoi|
|00001360| 64 20 74 72 79 6d 61 74 | 61 28 29 3b 20 76 6f 69 |d trymat|a(); voi|
|00001370| 64 20 74 72 79 6d 61 74 | 62 28 29 3b 20 76 6f 69 |d trymat|b(); voi|
|00001380| 64 20 74 72 79 6d 61 74 | 63 28 29 3b 0a 58 76 6f |d trymat|c();.Xvo|
|00001390| 69 64 20 74 72 79 6d 61 | 74 64 28 29 3b 20 76 6f |id tryma|td(); vo|
|000013a0| 69 64 20 74 72 79 6d 61 | 74 65 28 29 3b 20 76 6f |id tryma|te(); vo|
|000013b0| 69 64 20 74 72 79 6d 61 | 74 66 28 29 3b 0a 58 76 |id tryma|tf();.Xv|
|000013c0| 6f 69 64 20 74 72 79 6d | 61 74 67 28 29 3b 20 76 |oid trym|atg(); v|
|000013d0| 6f 69 64 20 74 72 79 6d | 61 74 68 28 29 3b 20 76 |oid trym|ath(); v|
|000013e0| 6f 69 64 20 74 72 79 6d | 61 74 69 28 29 3b 0a 58 |oid trym|ati();.X|
|000013f0| 0a 58 6d 61 69 6e 28 29 | 0a 58 7b 0a 58 20 20 20 |.Xmain()|.X{.X   |
|00001400| 52 65 61 6c 2a 20 73 31 | 3b 20 52 65 61 6c 2a 20 |Real* s1|; Real* |
|00001410| 73 32 3b 20 52 65 61 6c | 2a 20 73 33 3b 20 52 65 |s2; Real|* s3; Re|
|00001420| 61 6c 2a 20 73 34 3b 0a | 58 20 20 20 63 6f 75 74 |al* s4;.|X   cout|
|00001430| 20 3c 3c 20 22 5c 6e 42 | 65 67 69 6e 20 74 65 73 | << "\nB|egin tes|
|00001440| 74 5c 6e 22 3b 20 20 20 | 2f 2f 20 46 6f 72 63 65 |t\n";   |// Force|
|00001450| 73 20 63 6f 75 74 20 74 | 6f 20 61 6c 6c 6f 63 61 |s cout t|o alloca|
|00001460| 74 65 20 6d 65 6d 6f 72 | 79 20 61 74 20 62 65 67 |te memor|y at beg|
|00001470| 69 6e 6e 69 6e 67 0a 58 | 20 20 20 7b 20 4d 61 74 |inning.X|   { Mat|
|00001480| 72 69 78 20 41 31 28 34 | 30 2c 32 30 30 29 3b 20 |rix A1(4|0,200); |
|00001490| 73 31 20 3d 20 41 31 2e | 53 74 6f 72 65 28 29 3b |s1 = A1.|Store();|
|000014a0| 20 7d 0a 58 20 20 20 7b | 20 4d 61 74 72 69 78 20 | }.X   {| Matrix |
|000014b0| 41 31 28 31 2c 31 29 3b | 20 73 33 20 3d 20 41 31 |A1(1,1);| s3 = A1|
|000014c0| 2e 53 74 6f 72 65 28 29 | 3b 20 7d 0a 58 20 20 20 |.Store()|; }.X   |
|000014d0| 7b 0a 58 20 20 20 20 20 | 20 54 72 61 63 65 72 20 |{.X     | Tracer |
|000014e0| 65 74 28 22 4d 61 74 72 | 69 78 20 74 65 73 74 20 |et("Matr|ix test |
|000014f0| 70 72 6f 67 72 61 6d 22 | 29 3b 0a 58 0a 58 20 20 |program"|);.X.X  |
|00001500| 20 20 20 20 4d 61 74 72 | 69 78 20 41 28 32 35 2c |    Matr|ix A(25,|
|00001510| 31 35 30 29 3b 0a 58 20 | 20 20 20 20 20 7b 0a 58 |150);.X |     {.X|
|00001520| 09 20 69 6e 74 20 69 3b | 0a 58 09 20 52 6f 77 56 |. int i;|.X. RowV|
|00001530| 65 63 74 6f 72 20 41 28 | 38 29 3b 0a 58 09 20 66 |ector A(|8);.X. f|
|00001540| 6f 72 20 28 69 3d 31 3b | 69 3c 3d 37 3b 69 2b 2b |or (i=1;|i<=7;i++|
|00001550| 29 20 41 28 69 29 3d 30 | 2e 30 3b 20 41 28 38 29 |) A(i)=0|.0; A(8)|
|00001560| 3d 31 2e 30 3b 0a 58 09 | 20 50 72 69 6e 74 28 41 |=1.0;.X.| Print(A|
|00001570| 29 3b 0a 58 20 20 20 20 | 20 20 7d 0a 58 20 20 20 |);.X    |  }.X   |
|00001580| 20 20 20 63 6f 75 74 20 | 3c 3c 20 22 5c 6e 22 3b |   cout |<< "\n";|
|00001590| 0a 58 0a 58 20 20 20 20 | 20 20 54 65 73 74 54 79 |.X.X    |  TestTy|
|000015a0| 70 65 41 64 64 28 29 3b | 20 54 65 73 74 54 79 70 |peAdd();| TestTyp|
|000015b0| 65 4d 75 6c 74 28 29 3b | 20 54 65 73 74 54 79 70 |eMult();| TestTyp|
|000015c0| 65 4f 72 64 65 72 28 29 | 3b 0a 58 0a 58 20 20 20 |eOrder()|;.X.X   |
|000015d0| 20 20 20 74 72 79 6d 61 | 74 31 28 29 3b 0a 58 20 |   tryma|t1();.X |
|000015e0| 20 20 20 20 20 74 72 79 | 6d 61 74 32 28 29 3b 0a |     try|mat2();.|
|000015f0| 58 20 20 20 20 20 20 74 | 72 79 6d 61 74 33 28 29 |X      t|rymat3()|
|00001600| 3b 0a 58 20 20 20 20 20 | 20 74 72 79 6d 61 74 34 |;.X     | trymat4|
|00001610| 28 29 3b 0a 58 20 20 20 | 20 20 20 74 72 79 6d 61 |();.X   |   tryma|
|00001620| 74 35 28 29 3b 0a 58 20 | 20 20 20 20 20 74 72 79 |t5();.X |     try|
|00001630| 6d 61 74 36 28 29 3b 0a | 58 20 20 20 20 20 20 74 |mat6();.|X      t|
|00001640| 72 79 6d 61 74 37 28 29 | 3b 0a 58 20 20 20 20 20 |rymat7()|;.X     |
|00001650| 20 74 72 79 6d 61 74 38 | 28 29 3b 0a 58 20 20 20 | trymat8|();.X   |
|00001660| 20 20 20 74 72 79 6d 61 | 74 39 28 29 3b 0a 58 20 |   tryma|t9();.X |
|00001670| 20 20 20 20 20 74 72 79 | 6d 61 74 61 28 29 3b 0a |     try|mata();.|
|00001680| 58 20 20 20 20 20 20 74 | 72 79 6d 61 74 62 28 29 |X      t|rymatb()|
|00001690| 3b 0a 58 20 20 20 20 20 | 20 74 72 79 6d 61 74 63 |;.X     | trymatc|
|000016a0| 28 29 3b 0a 58 20 20 20 | 20 20 20 74 72 79 6d 61 |();.X   |   tryma|
|000016b0| 74 64 28 29 3b 0a 58 20 | 20 20 20 20 20 74 72 79 |td();.X |     try|
|000016c0| 6d 61 74 65 28 29 3b 0a | 58 20 20 20 20 20 20 74 |mate();.|X      t|
|000016d0| 72 79 6d 61 74 66 28 29 | 3b 0a 58 20 20 20 20 20 |rymatf()|;.X     |
|000016e0| 20 74 72 79 6d 61 74 67 | 28 29 3b 0a 58 20 20 20 | trymatg|();.X   |
|000016f0| 20 20 20 74 72 79 6d 61 | 74 68 28 29 3b 0a 58 20 |   tryma|th();.X |
|00001700| 20 20 20 20 20 74 72 79 | 6d 61 74 69 28 29 3b 0a |     try|mati();.|
|00001710| 58 20 20 20 7d 0a 58 20 | 20 20 7b 20 4d 61 74 72 |X   }.X |  { Matr|
|00001720| 69 78 20 41 31 28 34 30 | 2c 32 30 30 29 3b 20 73 |ix A1(40|,200); s|
|00001730| 32 20 3d 20 41 31 2e 53 | 74 6f 72 65 28 29 3b 20 |2 = A1.S|tore(); |
|00001740| 7d 0a 58 20 20 20 63 6f | 75 74 20 3c 3c 20 22 5c |}.X   co|ut << "\|
|00001750| 6e 43 68 65 63 6b 69 6e | 67 20 66 6f 72 20 6c 6f |nCheckin|g for lo|
|00001760| 73 74 20 6d 65 6d 6f 72 | 79 3a 20 22 0a 58 20 20 |st memor|y: ".X  |
|00001770| 20 20 20 20 3c 3c 20 28 | 75 6e 73 69 67 6e 65 64 |    << (|unsigned|
|00001780| 20 6c 6f 6e 67 29 73 31 | 20 3c 3c 20 22 20 22 20 | long)s1| << " " |
|00001790| 3c 3c 20 28 75 6e 73 69 | 67 6e 65 64 20 6c 6f 6e |<< (unsi|gned lon|
|000017a0| 67 29 73 32 20 3c 3c 20 | 22 20 22 3b 0a 58 20 20 |g)s2 << |" ";.X  |
|000017b0| 20 69 66 20 28 73 31 20 | 21 3d 20 73 32 29 20 63 | if (s1 |!= s2) c|
|000017c0| 6f 75 74 20 3c 3c 20 22 | 20 2d 20 65 72 72 6f 72 |out << "| - error|
|000017d0| 5c 6e 22 3b 20 65 6c 73 | 65 20 63 6f 75 74 20 3c |\n"; els|e cout <|
|000017e0| 3c 20 22 20 2d 20 6f 6b | 5c 6e 5c 6e 22 3b 0a 58 |< " - ok|\n\n";.X|
|000017f0| 20 20 20 7b 20 4d 61 74 | 72 69 78 20 41 31 28 31 |   { Mat|rix A1(1|
|00001800| 2c 31 29 3b 20 73 34 20 | 3d 20 41 31 2e 53 74 6f |,1); s4 |= A1.Sto|
|00001810| 72 65 28 29 3b 20 7d 0a | 58 20 20 20 63 6f 75 74 |re(); }.|X   cout|
|00001820| 20 3c 3c 20 22 5c 6e 43 | 68 65 63 6b 69 6e 67 20 | << "\nC|hecking |
|00001830| 66 6f 72 20 6c 6f 73 74 | 20 6d 65 6d 6f 72 79 3a |for lost| memory:|
|00001840| 20 22 0a 58 20 20 20 20 | 20 20 3c 3c 20 28 75 6e | ".X    |  << (un|
|00001850| 73 69 67 6e 65 64 20 6c | 6f 6e 67 29 73 33 20 3c |signed l|ong)s3 <|
|00001860| 3c 20 22 20 22 20 3c 3c | 20 28 75 6e 73 69 67 6e |< " " <<| (unsign|
|00001870| 65 64 20 6c 6f 6e 67 29 | 73 34 20 3c 3c 20 22 20 |ed long)|s4 << " |
|00001880| 22 3b 0a 58 20 20 20 69 | 66 20 28 73 33 20 21 3d |";.X   i|f (s3 !=|
|00001890| 20 73 34 29 20 63 6f 75 | 74 20 3c 3c 20 22 20 2d | s4) cou|t << " -|
|000018a0| 20 65 72 72 6f 72 5c 6e | 22 3b 20 65 6c 73 65 20 | error\n|"; else |
|000018b0| 63 6f 75 74 20 3c 3c 20 | 22 20 2d 20 6f 6b 5c 6e |cout << |" - ok\n|
|000018c0| 5c 6e 22 3b 0a 58 23 69 | 66 64 65 66 20 44 4f 5f |\n";.X#i|fdef DO_|
|000018d0| 46 52 45 45 5f 43 48 45 | 43 4b 0a 58 20 20 20 46 |FREE_CHE|CK.X   F|
|000018e0| 72 65 65 43 68 65 63 6b | 3a 3a 53 74 61 74 75 73 |reeCheck|::Status|
|000018f0| 28 29 3b 0a 58 23 65 6e | 64 69 66 20 0a 58 20 20 |();.X#en|dif .X  |
|00001900| 20 72 65 74 75 72 6e 20 | 30 3b 0a 58 7d 0a 45 4e | return |0;.X}.EN|
|00001910| 44 5f 4f 46 5f 46 49 4c | 45 0a 69 66 20 74 65 73 |D_OF_FIL|E.if tes|
|00001920| 74 20 34 36 33 33 20 2d | 6e 65 20 60 77 63 20 2d |t 4633 -|ne `wc -|
|00001930| 63 20 3c 27 74 6d 74 2e | 63 78 78 27 60 3b 20 74 |c <'tmt.|cxx'`; t|
|00001940| 68 65 6e 0a 20 20 20 20 | 65 63 68 6f 20 73 68 61 |hen.    |echo sha|
|00001950| 72 3a 20 5c 22 27 74 6d | 74 2e 63 78 78 27 5c 22 |r: \"'tm|t.cxx'\"|
|00001960| 20 75 6e 70 61 63 6b 65 | 64 20 77 69 74 68 20 77 | unpacke|d with w|
|00001970| 72 6f 6e 67 20 73 69 7a | 65 21 0a 66 69 0a 23 20 |rong siz|e!.fi.# |
|00001980| 65 6e 64 20 6f 66 20 27 | 74 6d 74 2e 63 78 78 27 |end of '|tmt.cxx'|
|00001990| 0a 66 69 0a 69 66 20 74 | 65 73 74 20 2d 66 20 27 |.fi.if t|est -f '|
|000019a0| 74 6d 74 31 2e 63 78 78 | 27 20 2d 61 20 22 24 7b |tmt1.cxx|' -a "${|
|000019b0| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|000019c0| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en .  ec|ho shar:|
|000019d0| 20 57 69 6c 6c 20 6e 6f | 74 20 63 6c 6f 62 62 65 | Will no|t clobbe|
|000019e0| 72 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |r existi|ng file |
|000019f0| 5c 22 27 74 6d 74 31 2e | 63 78 78 27 5c 22 0a 65 |\"'tmt1.|cxx'\".e|
|00001a00| 6c 73 65 0a 65 63 68 6f | 20 73 68 61 72 3a 20 45 |lse.echo| shar: E|
|00001a10| 78 74 72 61 63 74 69 6e | 67 20 5c 22 27 74 6d 74 |xtractin|g \"'tmt|
|00001a20| 31 2e 63 78 78 27 5c 22 | 20 5c 28 32 33 30 34 20 |1.cxx'\"| \(2304 |
|00001a30| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00001a40| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 74 6d 74 31 | "s/^X//|" >'tmt1|
|00001a50| 2e 63 78 78 27 20 3c 3c | 27 45 4e 44 5f 4f 46 5f |.cxx' <<|'END_OF_|
|00001a60| 46 49 4c 45 27 0a 58 0a | 58 2f 2f 23 64 65 66 69 |FILE'.X.|X//#defi|
|00001a70| 6e 65 20 57 41 4e 54 5f | 53 54 52 45 41 4d 0a 58 |ne WANT_|STREAM.X|
|00001a80| 0a 58 0a 58 0a 58 23 69 | 6e 63 6c 75 64 65 20 22 |.X.X.X#i|nclude "|
|00001a90| 69 6e 63 6c 75 64 65 2e | 68 22 0a 58 0a 58 23 69 |include.|h".X.X#i|
|00001aa0| 6e 63 6c 75 64 65 20 22 | 6e 65 77 6d 61 74 2e 68 |nclude "|newmat.h|
|00001ab0| 22 0a 58 0a 58 2f 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |".X.X/**|********|
|00001ac0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001ad0| 2a 2a 20 74 65 73 74 20 | 70 72 6f 67 72 61 6d 20 |** test |program |
|00001ae0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001af0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|00001b00| 58 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |X.Xvoid |Print(co|
|00001b10| 6e 73 74 20 4d 61 74 72 | 69 78 26 20 58 29 3b 0a |nst Matr|ix& X);.|
|00001b20| 58 76 6f 69 64 20 50 72 | 69 6e 74 28 63 6f 6e 73 |Xvoid Pr|int(cons|
|00001b30| 74 20 55 70 70 65 72 54 | 72 69 61 6e 67 75 6c 61 |t UpperT|riangula|
|00001b40| 72 4d 61 74 72 69 78 26 | 20 58 29 3b 0a 58 76 6f |rMatrix&| X);.Xvo|
|00001b50| 69 64 20 50 72 69 6e 74 | 28 63 6f 6e 73 74 20 44 |id Print|(const D|
|00001b60| 69 61 67 6f 6e 61 6c 4d | 61 74 72 69 78 26 20 58 |iagonalM|atrix& X|
|00001b70| 29 3b 0a 58 76 6f 69 64 | 20 50 72 69 6e 74 28 63 |);.Xvoid| Print(c|
|00001b80| 6f 6e 73 74 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |onst Sym|metricMa|
|00001b90| 74 72 69 78 26 20 58 29 | 3b 0a 58 76 6f 69 64 20 |trix& X)|;.Xvoid |
|00001ba0| 50 72 69 6e 74 28 63 6f | 6e 73 74 20 4c 6f 77 65 |Print(co|nst Lowe|
|00001bb0| 72 54 72 69 61 6e 67 75 | 6c 61 72 4d 61 74 72 69 |rTriangu|larMatri|
|00001bc0| 78 26 20 58 29 3b 0a 58 | 0a 58 0a 58 76 6f 69 64 |x& X);.X|.X.Xvoid|
|00001bd0| 20 74 72 79 6d 61 74 31 | 28 29 0a 58 7b 0a 58 2f | trymat1|().X{.X/|
|00001be0| 2f 20 20 20 63 6f 75 74 | 20 3c 3c 20 22 5c 6e 46 |/   cout| << "\nF|
|00001bf0| 69 72 73 74 20 74 65 73 | 74 20 6f 66 20 4d 61 74 |irst tes|t of Mat|
|00001c00| 72 69 78 20 70 61 63 6b | 61 67 65 5c 6e 5c 6e 22 |rix pack|age\n\n"|
|00001c10| 3b 0a 58 20 20 20 54 72 | 61 63 65 72 20 65 74 28 |;.X   Tr|acer et(|
|00001c20| 22 46 69 72 73 74 20 74 | 65 73 74 20 6f 66 20 4d |"First t|est of M|
|00001c30| 61 74 72 69 78 20 70 61 | 63 6b 61 67 65 22 29 3b |atrix pa|ckage");|
|00001c40| 0a 58 20 20 20 45 78 63 | 65 70 74 69 6f 6e 3a 3a |.X   Exc|eption::|
|00001c50| 50 72 69 6e 74 54 72 61 | 63 65 28 54 52 55 45 29 |PrintTra|ce(TRUE)|
|00001c60| 3b 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 54 |;.X   {.|X      T|
|00001c70| 72 61 63 65 72 20 65 74 | 31 28 22 53 74 61 67 65 |racer et|1("Stage|
|00001c80| 20 31 22 29 3b 0a 58 20 | 20 20 20 20 20 69 6e 74 | 1");.X |     int|
|00001c90| 20 69 2c 6a 3b 0a 58 0a | 58 20 20 20 20 20 20 4c | i,j;.X.|X      L|
|00001ca0| 6f 77 65 72 54 72 69 61 | 6e 67 75 6c 61 72 4d 61 |owerTria|ngularMa|
|00001cb0| 74 72 69 78 20 4c 28 31 | 30 29 3b 0a 58 20 20 20 |trix L(1|0);.X   |
|00001cc0| 20 20 20 66 6f 72 20 28 | 69 3d 31 3b 69 3c 3d 31 |   for (|i=1;i<=1|
|00001cd0| 30 3b 69 2b 2b 29 20 66 | 6f 72 20 28 6a 3d 31 3b |0;i++) f|or (j=1;|
|00001ce0| 6a 3c 3d 69 3b 6a 2b 2b | 29 20 4c 28 69 2c 6a 29 |j<=i;j++|) L(i,j)|
|00001cf0| 3d 32 2e 30 2b 69 2a 69 | 2b 6a 3b 0a 58 20 20 20 |=2.0+i*i|+j;.X   |
|00001d00| 20 20 20 53 79 6d 6d 65 | 74 72 69 63 4d 61 74 72 |   Symme|tricMatr|
|00001d10| 69 78 20 53 28 31 30 29 | 3b 0a 58 20 20 20 20 20 |ix S(10)|;.X     |
|00001d20| 20 66 6f 72 20 28 69 3d | 31 3b 69 3c 3d 31 30 3b | for (i=|1;i<=10;|
|00001d30| 69 2b 2b 29 20 66 6f 72 | 20 28 6a 3d 31 3b 6a 3c |i++) for| (j=1;j<|
|00001d40| 3d 69 3b 6a 2b 2b 29 20 | 53 28 69 2c 6a 29 3d 69 |=i;j++) |S(i,j)=i|
|00001d50| 2a 6a 2b 31 2e 30 3b 0a | 58 20 20 20 20 20 20 53 |*j+1.0;.|X      S|
|00001d60| 79 6d 6d 65 74 72 69 63 | 4d 61 74 72 69 78 20 53 |ymmetric|Matrix S|
|00001d70| 31 20 3d 20 53 20 2f 20 | 32 2e 30 3b 0a 58 20 20 |1 = S / |2.0;.X  |
|00001d80| 20 20 20 20 53 20 3d 20 | 53 31 20 2a 20 32 2e 30 |    S = |S1 * 2.0|
|00001d90| 3b 0a 58 20 20 20 20 20 | 20 55 70 70 65 72 54 72 |;.X     | UpperTr|
|00001da0| 69 61 6e 67 75 6c 61 72 | 4d 61 74 72 69 78 20 55 |iangular|Matrix U|
|00001db0| 3d 4c 2e 74 28 29 2a 32 | 2e 30 3b 0a 58 20 20 20 |=L.t()*2|.0;.X   |
|00001dc0| 20 20 20 50 72 69 6e 74 | 28 4c 6f 77 65 72 54 72 |   Print|(LowerTr|
|00001dd0| 69 61 6e 67 75 6c 61 72 | 4d 61 74 72 69 78 28 4c |iangular|Matrix(L|
|00001de0| 2d 55 2e 74 28 29 2a 30 | 2e 35 29 29 3b 0a 58 20 |-U.t()*0|.5));.X |
|00001df0| 20 20 20 20 20 44 69 61 | 67 6f 6e 61 6c 4d 61 74 |     Dia|gonalMat|
|00001e00| 72 69 78 20 44 28 31 30 | 29 3b 0a 58 20 20 20 20 |rix D(10|);.X    |
|00001e10| 20 20 66 6f 72 20 28 69 | 3d 31 3b 69 3c 3d 31 30 |  for (i|=1;i<=10|
|00001e20| 3b 69 2b 2b 29 20 44 28 | 69 2c 69 29 3d 28 69 2d |;i++) D(|i,i)=(i-|
|00001e30| 34 29 2a 28 69 2d 35 29 | 2a 28 69 2d 36 29 3b 0a |4)*(i-5)|*(i-6);.|
|00001e40| 58 20 20 20 20 20 20 4d | 61 74 72 69 78 20 4d 3d |X      M|atrix M=|
|00001e50| 28 53 2b 55 2d 44 2b 4c | 29 2a 28 4c 2b 55 2d 44 |(S+U-D+L|)*(L+U-D|
|00001e60| 2b 53 29 3b 0a 58 20 20 | 20 20 20 20 44 69 61 67 |+S);.X  |    Diag|
|00001e70| 6f 6e 61 6c 4d 61 74 72 | 69 78 20 44 44 3d 44 2a |onalMatr|ix DD=D*|
|00001e80| 44 3b 0a 58 20 20 20 20 | 20 20 4c 6f 77 65 72 54 |D;.X    |  LowerT|
|00001e90| 72 69 61 6e 67 75 6c 61 | 72 4d 61 74 72 69 78 20 |riangula|rMatrix |
|00001ea0| 4c 44 3d 4c 2a 44 3b 0a | 58 20 20 20 20 20 20 2f |LD=L*D;.|X      /|
|00001eb0| 2f 20 65 78 70 72 65 73 | 73 69 6f 6e 73 20 73 70 |/ expres|sions sp|
|00001ec0| 6c 69 74 20 66 6f 72 20 | 54 75 72 62 6f 20 43 0a |lit for |Turbo C.|
|00001ed0| 58 20 20 20 20 20 20 4d | 61 74 72 69 78 20 4d 31 |X      M|atrix M1|
|00001ee0| 20 3d 20 53 2a 4c 20 2b | 20 55 2a 4c 20 2d 20 44 | = S*L +| U*L - D|
|00001ef0| 2a 4c 20 2b 20 4c 2a 4c | 20 2b 20 31 30 2e 30 3b |*L + L*L| + 10.0;|
|00001f00| 0a 58 20 20 20 20 20 20 | 7b 20 4d 31 20 3d 20 4d |.X      |{ M1 = M|
|00001f10| 31 20 2b 20 53 2a 55 20 | 2b 20 55 2a 55 20 2d 20 |1 + S*U |+ U*U - |
|00001f20| 44 2a 55 20 2b 20 4c 2a | 55 20 2d 20 53 2a 44 3b |D*U + L*|U - S*D;|
|00001f30| 20 7d 0a 58 20 20 20 20 | 20 20 7b 20 4d 31 20 3d | }.X    |  { M1 =|
|00001f40| 20 4d 31 20 2d 20 55 2a | 44 20 2b 20 44 44 20 2d | M1 - U*|D + DD -|
|00001f50| 20 4c 44 20 2b 20 53 2a | 53 3b 20 7d 0a 58 20 20 | LD + S*|S; }.X  |
|00001f60| 20 20 20 20 7b 20 4d 31 | 20 3d 20 4d 31 20 2b 20 |    { M1| = M1 + |
|00001f70| 55 2a 53 20 2d 20 44 2a | 53 20 2b 20 4c 2a 53 20 |U*S - D*|S + L*S |
|00001f80| 2d 20 31 30 2e 30 3b 20 | 7d 0a 58 20 20 20 20 20 |- 10.0; |}.X     |
|00001f90| 20 4d 3d 4d 31 2d 4d 3b | 0a 58 20 20 20 20 20 20 | M=M1-M;|.X      |
|00001fa0| 50 72 69 6e 74 28 4d 29 | 3b 0a 58 20 20 20 7d 0a |Print(M)|;.X   }.|
|00001fb0| 58 20 20 20 7b 0a 58 20 | 20 20 20 20 20 54 72 61 |X   {.X |     Tra|
|00001fc0| 63 65 72 20 65 74 31 28 | 22 53 74 61 67 65 20 32 |cer et1(|"Stage 2|
|00001fd0| 22 29 3b 0a 58 20 20 20 | 20 20 20 69 6e 74 20 69 |");.X   |   int i|
|00001fe0| 2c 6a 3b 0a 58 0a 58 20 | 20 20 20 20 20 4c 6f 77 |,j;.X.X |     Low|
|00001ff0| 65 72 54 72 69 61 6e 67 | 75 6c 61 72 4d 61 74 72 |erTriang|ularMatr|
|00002000| 69 78 20 4c 28 39 29 3b | 0a 58 20 20 20 20 20 20 |ix L(9);|.X      |
|00002010| 66 6f 72 20 28 69 3d 31 | 3b 69 3c 3d 39 3b 69 2b |for (i=1|;i<=9;i+|
|00002020| 2b 29 20 66 6f 72 20 28 | 6a 3d 31 3b 6a 3c 3d 69 |+) for (|j=1;j<=i|
|00002030| 3b 6a 2b 2b 29 20 4c 28 | 69 2c 6a 29 3d 31 2e 30 |;j++) L(|i,j)=1.0|
|00002040| 2b 6a 3b 0a 58 20 20 20 | 20 20 20 55 70 70 65 72 |+j;.X   |   Upper|
|00002050| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|00002060| 20 55 31 28 39 29 3b 0a | 58 20 20 20 20 20 20 66 | U1(9);.|X      f|
|00002070| 6f 72 20 28 6a 3d 31 3b | 6a 3c 3d 39 3b 6a 2b 2b |or (j=1;|j<=9;j++|
|00002080| 29 20 66 6f 72 20 28 69 | 3d 31 3b 69 3c 3d 6a 3b |) for (i|=1;i<=j;|
|00002090| 69 2b 2b 29 20 55 31 28 | 69 2c 6a 29 3d 31 2e 30 |i++) U1(|i,j)=1.0|
|000020a0| 2b 69 3b 0a 58 20 20 20 | 20 20 20 4c 6f 77 65 72 |+i;.X   |   Lower|
|000020b0| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|000020c0| 20 4c 58 28 39 29 3b 0a | 58 20 20 20 20 20 20 66 | LX(9);.|X      f|
|000020d0| 6f 72 20 28 69 3d 31 3b | 69 3c 3d 39 3b 69 2b 2b |or (i=1;|i<=9;i++|
|000020e0| 29 20 66 6f 72 20 28 6a | 3d 31 3b 6a 3c 3d 69 3b |) for (j|=1;j<=i;|
|000020f0| 6a 2b 2b 29 20 4c 58 28 | 69 2c 6a 29 3d 31 2e 30 |j++) LX(|i,j)=1.0|
|00002100| 2b 69 2a 69 3b 0a 58 20 | 20 20 20 20 20 55 70 70 |+i*i;.X |     Upp|
|00002110| 65 72 54 72 69 61 6e 67 | 75 6c 61 72 4d 61 74 72 |erTriang|ularMatr|
|00002120| 69 78 20 55 58 28 39 29 | 3b 0a 58 20 20 20 20 20 |ix UX(9)|;.X     |
|00002130| 20 66 6f 72 20 28 6a 3d | 31 3b 6a 3c 3d 39 3b 6a | for (j=|1;j<=9;j|
|00002140| 2b 2b 29 20 66 6f 72 20 | 28 69 3d 31 3b 69 3c 3d |++) for |(i=1;i<=|
|00002150| 6a 3b 69 2b 2b 29 20 55 | 58 28 69 2c 6a 29 3d 31 |j;i++) U|X(i,j)=1|
|00002160| 2e 30 2b 6a 2a 6a 3b 0a | 58 20 20 20 20 20 20 7b |.0+j*j;.|X      {|
|00002170| 0a 58 20 20 20 20 20 20 | 20 20 20 4c 3d 4c 2b 4c |.X      |   L=L+L|
|00002180| 58 2f 30 2e 35 3b 20 20 | 20 4c 3d 4c 2d 4c 58 2a |X/0.5;  | L=L-LX*|
|00002190| 33 2e 30 3b 20 20 20 4c | 3d 4c 58 2a 32 2e 30 2b |3.0;   L|=LX*2.0+|
|000021a0| 4c 3b 0a 58 20 20 20 20 | 20 20 20 20 20 55 31 3d |L;.X    |     U1=|
|000021b0| 55 31 2b 55 58 2a 32 2e | 30 3b 20 55 31 3d 55 31 |U1+UX*2.|0; U1=U1|
|000021c0| 2d 55 58 2a 33 2e 30 3b | 20 55 31 3d 55 58 2a 32 |-UX*3.0;| U1=UX*2|
|000021d0| 2e 30 2b 55 31 3b 0a 58 | 20 20 20 20 20 20 7d 0a |.0+U1;.X|      }.|
|000021e0| 58 0a 58 0a 58 20 20 20 | 20 20 20 53 79 6d 6d 65 |X.X.X   |   Symme|
|000021f0| 74 72 69 63 4d 61 74 72 | 69 78 20 53 28 39 29 3b |tricMatr|ix S(9);|
|00002200| 0a 58 20 20 20 20 20 20 | 66 6f 72 20 28 69 3d 31 |.X      |for (i=1|
|00002210| 3b 69 3c 3d 39 3b 69 2b | 2b 29 20 66 6f 72 20 28 |;i<=9;i+|+) for (|
|00002220| 6a 3d 31 3b 6a 3c 3d 69 | 3b 6a 2b 2b 29 20 53 28 |j=1;j<=i|;j++) S(|
|00002230| 69 2c 6a 29 3d 69 2a 69 | 2b 6a 3b 0a 58 20 20 20 |i,j)=i*i|+j;.X   |
|00002240| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 20 20 20 53 |   {.X  |       S|
|00002250| 79 6d 6d 65 74 72 69 63 | 4d 61 74 72 69 78 20 53 |ymmetric|Matrix S|
|00002260| 31 20 3d 20 53 3b 0a 58 | 20 20 20 20 20 20 20 20 |1 = S;.X|        |
|00002270| 20 53 3d 53 31 2b 35 2e | 30 3b 0a 58 20 20 20 20 | S=S1+5.|0;.X    |
|00002280| 20 20 20 20 20 53 3d 53 | 2d 33 2e 30 3b 0a 58 20 |     S=S|-3.0;.X |
|00002290| 20 20 20 20 20 7d 0a 58 | 0a 58 20 20 20 20 20 20 |     }.X|.X      |
|000022a0| 44 69 61 67 6f 6e 61 6c | 4d 61 74 72 69 78 20 44 |Diagonal|Matrix D|
|000022b0| 28 39 29 3b 0a 58 20 20 | 20 20 20 20 66 6f 72 20 |(9);.X  |    for |
|000022c0| 28 69 3d 31 3b 69 3c 3d | 39 3b 69 2b 2b 29 20 44 |(i=1;i<=|9;i++) D|
|000022d0| 28 69 2c 69 29 3d 53 28 | 69 2c 69 29 3b 0a 58 20 |(i,i)=S(|i,i);.X |
|000022e0| 20 20 20 20 20 55 70 70 | 65 72 54 72 69 61 6e 67 |     Upp|erTriang|
|000022f0| 75 6c 61 72 4d 61 74 72 | 69 78 20 55 3d 4c 2e 74 |ularMatr|ix U=L.t|
|00002300| 28 29 2a 32 2e 30 3b 0a | 58 20 20 20 20 20 20 7b |()*2.0;.|X      {|
|00002310| 0a 58 20 20 20 20 20 20 | 20 20 20 55 31 3d 55 31 |.X      |   U1=U1|
|00002320| 2a 32 2e 30 20 2d 20 55 | 3b 20 20 50 72 69 6e 74 |*2.0 - U|;  Print|
|00002330| 28 55 31 29 3b 0a 58 20 | 20 20 20 20 20 20 20 20 |(U1);.X |        |
|00002340| 4c 3d 4c 2a 32 2e 30 2d | 44 3b 20 55 3d 55 2d 44 |L=L*2.0-|D; U=U-D|
|00002350| 3b 0a 58 20 20 20 20 20 | 20 7d 0a 58 20 20 20 20 |;.X     | }.X    |
|00002360| 20 20 4d 61 74 72 69 78 | 20 4d 3d 55 2b 4c 3b 20 |  Matrix| M=U+L; |
|00002370| 53 3d 53 2a 32 2e 30 3b | 20 4d 3d 53 2d 4d 3b 20 |S=S*2.0;| M=S-M; |
|00002380| 50 72 69 6e 74 28 4d 29 | 3b 0a 58 20 20 20 7d 0a |Print(M)|;.X   }.|
|00002390| 58 2f 2f 20 20 20 63 6f | 75 74 20 3c 3c 20 22 5c |X//   co|ut << "\|
|000023a0| 6e 45 6e 64 20 6f 66 20 | 66 69 72 73 74 20 74 65 |nEnd of |first te|
|000023b0| 73 74 5c 6e 22 3b 0a 58 | 7d 0a 58 0a 45 4e 44 5f |st\n";.X|}.X.END_|
|000023c0| 4f 46 5f 46 49 4c 45 0a | 69 66 20 74 65 73 74 20 |OF_FILE.|if test |
|000023d0| 32 33 30 34 20 2d 6e 65 | 20 60 77 63 20 2d 63 20 |2304 -ne| `wc -c |
|000023e0| 3c 27 74 6d 74 31 2e 63 | 78 78 27 60 3b 20 74 68 |<'tmt1.c|xx'`; th|
|000023f0| 65 6e 0a 20 20 20 20 65 | 63 68 6f 20 73 68 61 72 |en.    e|cho shar|
|00002400| 3a 20 5c 22 27 74 6d 74 | 31 2e 63 78 78 27 5c 22 |: \"'tmt|1.cxx'\"|
|00002410| 20 75 6e 70 61 63 6b 65 | 64 20 77 69 74 68 20 77 | unpacke|d with w|
|00002420| 72 6f 6e 67 20 73 69 7a | 65 21 0a 66 69 0a 23 20 |rong siz|e!.fi.# |
|00002430| 65 6e 64 20 6f 66 20 27 | 74 6d 74 31 2e 63 78 78 |end of '|tmt1.cxx|
|00002440| 27 0a 66 69 0a 69 66 20 | 74 65 73 74 20 2d 66 20 |'.fi.if |test -f |
|00002450| 27 74 6d 74 32 2e 63 78 | 78 27 20 2d 61 20 22 24 |'tmt2.cx|x' -a "$|
|00002460| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|00002470| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen .  e|cho shar|
|00002480| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 63 6c 6f 62 62 |: Will n|ot clobb|
|00002490| 65 72 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |er exist|ing file|
|000024a0| 20 5c 22 27 74 6d 74 32 | 2e 63 78 78 27 5c 22 0a | \"'tmt2|.cxx'\".|
|000024b0| 65 6c 73 65 0a 65 63 68 | 6f 20 73 68 61 72 3a 20 |else.ech|o shar: |
|000024c0| 45 78 74 72 61 63 74 69 | 6e 67 20 5c 22 27 74 6d |Extracti|ng \"'tm|
|000024d0| 74 32 2e 63 78 78 27 5c | 22 20 5c 28 32 32 36 34 |t2.cxx'\|" \(2264|
|000024e0| 20 63 68 61 72 61 63 74 | 65 72 73 5c 29 0a 73 65 | charact|ers\).se|
|000024f0| 64 20 22 73 2f 5e 58 2f | 2f 22 20 3e 27 74 6d 74 |d "s/^X/|/" >'tmt|
|00002500| 32 2e 63 78 78 27 20 3c | 3c 27 45 4e 44 5f 4f 46 |2.cxx' <|<'END_OF|
|00002510| 5f 46 49 4c 45 27 0a 58 | 0a 58 2f 2f 23 64 65 66 |_FILE'.X|.X//#def|
|00002520| 69 6e 65 20 57 41 4e 54 | 5f 53 54 52 45 41 4d 0a |ine WANT|_STREAM.|
|00002530| 58 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 20 22 69 |X.X.X#in|clude "i|
|00002540| 6e 63 6c 75 64 65 2e 68 | 22 0a 58 0a 58 23 69 6e |nclude.h|".X.X#in|
|00002550| 63 6c 75 64 65 20 22 6e | 65 77 6d 61 74 2e 68 22 |clude "n|ewmat.h"|
|00002560| 0a 58 0a 58 2f 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X.X/***|********|
|00002570| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002580| 2a 20 74 65 73 74 20 70 | 72 6f 67 72 61 6d 20 2a |* test p|rogram *|
|00002590| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000025a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2f 0a 58 |********|*****/.X|
|000025b0| 0a 58 76 6f 69 64 20 50 | 72 69 6e 74 28 63 6f 6e |.Xvoid P|rint(con|
|000025c0| 73 74 20 4d 61 74 72 69 | 78 26 20 58 29 3b 0a 58 |st Matri|x& X);.X|
|000025d0| 76 6f 69 64 20 50 72 69 | 6e 74 28 63 6f 6e 73 74 |void Pri|nt(const|
|000025e0| 20 55 70 70 65 72 54 72 | 69 61 6e 67 75 6c 61 72 | UpperTr|iangular|
|000025f0| 4d 61 74 72 69 78 26 20 | 58 29 3b 0a 58 76 6f 69 |Matrix& |X);.Xvoi|
|00002600| 64 20 50 72 69 6e 74 28 | 63 6f 6e 73 74 20 44 69 |d Print(|const Di|
|00002610| 61 67 6f 6e 61 6c 4d 61 | 74 72 69 78 26 20 58 29 |agonalMa|trix& X)|
|00002620| 3b 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |;.Xvoid |Print(co|
|00002630| 6e 73 74 20 53 79 6d 6d | 65 74 72 69 63 4d 61 74 |nst Symm|etricMat|
|00002640| 72 69 78 26 20 58 29 3b | 0a 58 76 6f 69 64 20 50 |rix& X);|.Xvoid P|
|00002650| 72 69 6e 74 28 63 6f 6e | 73 74 20 4c 6f 77 65 72 |rint(con|st Lower|
|00002660| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|00002670| 26 20 58 29 3b 0a 58 0a | 58 0a 58 76 6f 69 64 20 |& X);.X.|X.Xvoid |
|00002680| 74 72 79 6d 61 74 32 28 | 29 0a 58 7b 0a 58 2f 2f |trymat2(|).X{.X//|
|00002690| 20 20 20 63 6f 75 74 20 | 3c 3c 20 22 5c 6e 53 65 |   cout |<< "\nSe|
|000026a0| 63 6f 6e 64 20 74 65 73 | 74 20 6f 66 20 4d 61 74 |cond tes|t of Mat|
|000026b0| 72 69 78 20 70 61 63 6b | 61 67 65 5c 6e 5c 6e 22 |rix pack|age\n\n"|
|000026c0| 3b 0a 58 20 20 20 54 72 | 61 63 65 72 20 65 74 28 |;.X   Tr|acer et(|
|000026d0| 22 53 65 63 6f 6e 64 20 | 74 65 73 74 20 6f 66 20 |"Second |test of |
|000026e0| 4d 61 74 72 69 78 20 70 | 61 63 6b 61 67 65 22 29 |Matrix p|ackage")|
|000026f0| 3b 0a 58 20 20 20 45 78 | 63 65 70 74 69 6f 6e 3a |;.X   Ex|ception:|
|00002700| 3a 50 72 69 6e 74 54 72 | 61 63 65 28 54 52 55 45 |:PrintTr|ace(TRUE|
|00002710| 29 3b 0a 58 0a 58 20 20 | 20 69 6e 74 20 69 2c 6a |);.X.X  | int i,j|
|00002720| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 20 4d 28 33 |;.X   Ma|trix M(3|
|00002730| 2c 35 29 3b 0a 58 20 20 | 20 66 6f 72 20 28 69 3d |,5);.X  | for (i=|
|00002740| 31 3b 20 69 3c 3d 33 3b | 20 69 2b 2b 29 20 66 6f |1; i<=3;| i++) fo|
|00002750| 72 20 28 6a 3d 31 3b 20 | 6a 3c 3d 35 3b 20 6a 2b |r (j=1; |j<=5; j+|
|00002760| 2b 29 20 4d 28 69 2c 6a | 29 20 3d 20 31 30 30 2a |+) M(i,j|) = 100*|
|00002770| 69 20 2b 20 6a 3b 0a 58 | 20 20 20 4d 61 74 72 69 |i + j;.X|   Matri|
|00002780| 78 20 58 28 38 2c 31 30 | 29 3b 0a 58 20 20 20 66 |x X(8,10|);.X   f|
|00002790| 6f 72 20 28 69 3d 31 3b | 20 69 3c 3d 38 3b 20 69 |or (i=1;| i<=8; i|
|000027a0| 2b 2b 29 20 66 6f 72 20 | 28 6a 3d 31 3b 20 6a 3c |++) for |(j=1; j<|
|000027b0| 3d 31 30 3b 20 6a 2b 2b | 29 20 58 28 69 2c 6a 29 |=10; j++|) X(i,j)|
|000027c0| 20 3d 20 31 30 30 30 2a | 69 20 2b 20 31 30 2a 6a | = 1000*|i + 10*j|
|000027d0| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 20 59 20 3d |;.X   Ma|trix Y =|
|000027e0| 20 58 3b 20 4d 61 74 72 | 69 78 20 5a 20 3d 20 58 | X; Matr|ix Z = X|
|000027f0| 3b 0a 58 20 20 20 7b 20 | 58 2e 53 75 62 4d 61 74 |;.X   { |X.SubMat|
|00002800| 72 69 78 28 32 2c 34 2c | 33 2c 37 29 20 3c 3c 20 |rix(2,4,|3,7) << |
|00002810| 4d 3b 20 7d 0a 58 20 20 | 20 66 6f 72 20 28 69 3d |M; }.X  | for (i=|
|00002820| 31 3b 20 69 3c 3d 33 3b | 20 69 2b 2b 29 20 66 6f |1; i<=3;| i++) fo|
|00002830| 72 20 28 6a 3d 31 3b 20 | 6a 3c 3d 35 3b 20 6a 2b |r (j=1; |j<=5; j+|
|00002840| 2b 29 20 59 28 69 2b 31 | 2c 6a 2b 32 29 20 3d 20 |+) Y(i+1|,j+2) = |
|00002850| 31 30 30 2a 69 20 2b 20 | 6a 3b 0a 58 20 20 20 50 |100*i + |j;.X   P|
|00002860| 72 69 6e 74 28 4d 61 74 | 72 69 78 28 58 2d 59 29 |rint(Mat|rix(X-Y)|
|00002870| 29 3b 0a 58 0a 58 0a 58 | 20 20 20 52 65 61 6c 20 |);.X.X.X|   Real |
|00002880| 61 5b 31 35 5d 3b 20 52 | 65 61 6c 2a 20 72 20 3d |a[15]; R|eal* r =|
|00002890| 20 61 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 3d 31 | a;.X   |for (i=1|
|000028a0| 3b 20 69 3c 3d 33 3b 20 | 69 2b 2b 29 20 66 6f 72 |; i<=3; |i++) for|
|000028b0| 20 28 6a 3d 31 3b 20 6a | 3c 3d 35 3b 20 6a 2b 2b | (j=1; j|<=5; j++|
|000028c0| 29 20 2a 72 2b 2b 20 3d | 20 31 30 30 2a 69 20 2b |) *r++ =| 100*i +|
|000028d0| 20 6a 3b 0a 58 20 20 20 | 7b 20 5a 2e 53 75 62 4d | j;.X   |{ Z.SubM|
|000028e0| 61 74 72 69 78 28 32 2c | 34 2c 33 2c 37 29 20 3c |atrix(2,|4,3,7) <|
|000028f0| 3c 20 61 3b 20 7d 0a 58 | 20 20 20 50 72 69 6e 74 |< a; }.X|   Print|
|00002900| 28 4d 61 74 72 69 78 28 | 5a 2d 59 29 29 3b 0a 58 |(Matrix(|Z-Y));.X|
|00002910| 0a 58 20 20 20 7b 20 4d | 3d 33 33 3b 20 58 2e 53 |.X   { M|=33; X.S|
|00002920| 75 62 4d 61 74 72 69 78 | 28 32 2c 34 2c 33 2c 37 |ubMatrix|(2,4,3,7|
|00002930| 29 20 3c 3c 20 4d 3b 20 | 7d 0a 58 20 20 20 7b 20 |) << M; |}.X   { |
|00002940| 5a 2e 53 75 62 4d 61 74 | 72 69 78 28 32 2c 34 2c |Z.SubMat|rix(2,4,|
|00002950| 33 2c 37 29 20 3d 20 33 | 33 3b 20 7d 0a 58 20 20 |3,7) = 3|3; }.X  |
|00002960| 20 50 72 69 6e 74 28 4d | 61 74 72 69 78 28 5a 2d | Print(M|atrix(Z-|
|00002970| 58 29 29 3b 0a 58 0a 58 | 20 20 20 66 6f 72 20 28 |X));.X.X|   for (|
|00002980| 69 3d 31 3b 20 69 3c 3d | 38 3b 20 69 2b 2b 29 20 |i=1; i<=|8; i++) |
|00002990| 66 6f 72 20 28 6a 3d 31 | 3b 20 6a 3c 3d 31 30 3b |for (j=1|; j<=10;|
|000029a0| 20 6a 2b 2b 29 20 58 28 | 69 2c 6a 29 20 3d 20 31 | j++) X(|i,j) = 1|
|000029b0| 30 30 30 2a 69 20 2b 20 | 31 30 2a 6a 3b 0a 58 20 |000*i + |10*j;.X |
|000029c0| 20 20 59 20 3d 20 58 3b | 0a 58 20 20 20 55 70 70 |  Y = X;|.X   Upp|
|000029d0| 65 72 54 72 69 61 6e 67 | 75 6c 61 72 4d 61 74 72 |erTriang|ularMatr|
|000029e0| 69 78 20 55 28 35 29 3b | 0a 58 20 20 20 66 6f 72 |ix U(5);|.X   for|
|000029f0| 20 28 69 3d 31 3b 20 69 | 3c 3d 35 3b 20 69 2b 2b | (i=1; i|<=5; i++|
|00002a00| 29 20 66 6f 72 20 28 6a | 3d 69 3b 20 6a 3c 3d 35 |) for (j|=i; j<=5|
|00002a10| 3b 20 6a 2b 2b 29 20 55 | 28 69 2c 6a 29 20 3d 20 |; j++) U|(i,j) = |
|00002a20| 31 30 30 2a 69 20 2b 20 | 6a 3b 0a 58 20 20 20 7b |100*i + |j;.X   {|
|00002a30| 20 58 2e 53 75 62 4d 61 | 74 72 69 78 28 33 2c 37 | X.SubMa|trix(3,7|
|00002a40| 2c 35 2c 39 29 20 3c 3c | 20 55 3b 20 7d 0a 58 20 |,5,9) <<| U; }.X |
|00002a50| 20 20 66 6f 72 20 28 69 | 3d 31 3b 20 69 3c 3d 35 |  for (i|=1; i<=5|
|00002a60| 3b 20 69 2b 2b 29 20 66 | 6f 72 20 28 6a 3d 69 3b |; i++) f|or (j=i;|
|00002a70| 20 6a 3c 3d 35 3b 20 6a | 2b 2b 29 20 59 28 69 2b | j<=5; j|++) Y(i+|
|00002a80| 32 2c 6a 2b 34 29 20 3d | 20 31 30 30 2a 69 20 2b |2,j+4) =| 100*i +|
|00002a90| 20 6a 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 3d 31 | j;.X   |for (i=1|
|00002aa0| 3b 20 69 3c 3d 35 3b 20 | 69 2b 2b 29 20 66 6f 72 |; i<=5; |i++) for|
|00002ab0| 20 28 6a 3d 31 3b 20 6a | 3c 69 3b 20 6a 2b 2b 29 | (j=1; j|<i; j++)|
|00002ac0| 20 59 28 69 2b 32 2c 6a | 2b 34 29 20 3d 20 30 2e | Y(i+2,j|+4) = 0.|
|00002ad0| 30 3b 0a 58 20 20 20 50 | 72 69 6e 74 28 4d 61 74 |0;.X   P|rint(Mat|
|00002ae0| 72 69 78 28 58 2d 59 29 | 29 3b 0a 58 20 20 20 66 |rix(X-Y)|);.X   f|
|00002af0| 6f 72 20 28 69 3d 31 3b | 20 69 3c 3d 38 3b 20 69 |or (i=1;| i<=8; i|
|00002b00| 2b 2b 29 20 66 6f 72 20 | 28 6a 3d 31 3b 20 6a 3c |++) for |(j=1; j<|
|00002b10| 3d 31 30 3b 20 6a 2b 2b | 29 20 58 28 69 2c 6a 29 |=10; j++|) X(i,j)|
|00002b20| 20 3d 20 31 30 30 30 2a | 69 20 2b 20 31 30 2a 6a | = 1000*|i + 10*j|
|00002b30| 3b 0a 58 20 20 20 59 20 | 3d 20 58 3b 0a 58 20 20 |;.X   Y |= X;.X  |
|00002b40| 20 66 6f 72 20 28 69 3d | 31 3b 20 69 3c 3d 35 3b | for (i=|1; i<=5;|
|00002b50| 20 69 2b 2b 29 20 66 6f | 72 20 28 6a 3d 69 3b 20 | i++) fo|r (j=i; |
|00002b60| 6a 3c 3d 35 3b 20 6a 2b | 2b 29 20 55 28 69 2c 6a |j<=5; j+|+) U(i,j|
|00002b70| 29 20 3d 20 31 30 30 2a | 69 20 2b 20 6a 3b 0a 58 |) = 100*|i + j;.X|
|00002b80| 20 20 20 7b 20 58 2e 53 | 75 62 4d 61 74 72 69 78 |   { X.S|ubMatrix|
|00002b90| 28 33 2c 37 2c 35 2c 39 | 29 2e 49 6e 6a 65 63 74 |(3,7,5,9|).Inject|
|00002ba0| 28 55 29 3b 20 7d 0a 58 | 20 20 20 66 6f 72 20 28 |(U); }.X|   for (|
|00002bb0| 69 3d 31 3b 20 69 3c 3d | 35 3b 20 69 2b 2b 29 20 |i=1; i<=|5; i++) |
|00002bc0| 66 6f 72 20 28 6a 3d 69 | 3b 20 6a 3c 3d 35 3b 20 |for (j=i|; j<=5; |
|00002bd0| 6a 2b 2b 29 20 59 28 69 | 2b 32 2c 6a 2b 34 29 20 |j++) Y(i|+2,j+4) |
|00002be0| 3d 20 31 30 30 2a 69 20 | 2b 20 6a 3b 0a 58 20 20 |= 100*i |+ j;.X  |
|00002bf0| 20 50 72 69 6e 74 28 4d | 61 74 72 69 78 28 58 2d | Print(M|atrix(X-|
|00002c00| 59 29 29 3b 0a 58 0a 58 | 0a 58 20 20 20 2f 2f 20 |Y));.X.X|.X   // |
|00002c10| 74 65 73 74 20 67 72 6f | 77 69 6e 67 20 61 6e 64 |test gro|wing and|
|00002c20| 20 73 68 72 69 6e 6b 69 | 6e 67 20 61 20 76 65 63 | shrinki|ng a vec|
|00002c30| 74 6f 72 0a 58 20 20 20 | 43 6f 6c 75 6d 6e 56 65 |tor.X   |ColumnVe|
|00002c40| 63 74 6f 72 20 56 28 31 | 30 30 29 3b 0a 58 20 20 |ctor V(1|00);.X  |
|00002c50| 20 66 6f 72 20 28 69 3d | 31 3b 69 3c 3d 31 30 30 | for (i=|1;i<=100|
|00002c60| 3b 69 2b 2b 29 20 56 28 | 69 29 20 3d 20 69 2a 69 |;i++) V(|i) = i*i|
|00002c70| 2b 69 3b 0a 58 20 20 20 | 56 20 3d 20 56 2e 52 6f |+i;.X   |V = V.Ro|
|00002c80| 77 73 28 31 2c 35 30 29 | 3b 20 20 20 20 20 20 20 |ws(1,50)|;       |
|00002c90| 20 20 20 20 20 20 20 20 | 2f 2f 20 74 6f 20 67 65 |        |// to ge|
|00002ca0| 74 20 66 69 72 73 74 20 | 35 30 20 76 6c 61 75 65 |t first |50 vlaue|
|00002cb0| 73 2e 0a 58 0a 58 20 20 | 20 7b 0a 58 20 20 20 20 |s..X.X  | {.X    |
|00002cc0| 20 20 56 2e 52 65 6c 65 | 61 73 65 28 29 3b 20 43 |  V.Rele|ase(); C|
|00002cd0| 6f 6c 75 6d 6e 56 65 63 | 74 6f 72 20 56 58 3d 56 |olumnVec|tor VX=V|
|00002ce0| 3b 0a 58 20 20 20 20 20 | 20 56 2e 52 65 44 69 6d |;.X     | V.ReDim|
|00002cf0| 65 6e 73 69 6f 6e 28 31 | 30 30 29 3b 20 56 20 3d |ension(1|00); V =|
|00002d00| 20 30 2e 30 3b 20 56 2e | 52 6f 77 73 28 31 2c 35 | 0.0; V.|Rows(1,5|
|00002d10| 30 29 3d 56 58 3b 0a 58 | 20 20 20 7d 20 20 20 20 |0)=VX;.X|   }    |
|00002d20| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|00002d30| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 56 20 |        |   // V |
|00002d40| 6e 6f 77 20 6c 65 6e 67 | 74 68 20 31 30 30 0a 58 |now leng|th 100.X|
|00002d50| 0a 58 23 69 66 20 5f 5f | 5a 54 43 5f 5f 0a 58 20 |.X#if __|ZTC__.X |
|00002d60| 20 20 43 6f 6c 75 6d 6e | 56 65 63 74 6f 72 20 56 |  Column|Vector V|
|00002d70| 4d 20 3d 20 56 3b 20 56 | 4d 20 3d 20 31 30 30 3b |M = V; V|M = 100;|
|00002d80| 0a 58 23 65 6c 73 65 0a | 58 20 20 20 4d 3d 56 3b |.X#else.|X   M=V;|
|00002d90| 20 4d 3d 31 30 30 3b 20 | 20 20 20 20 20 20 20 20 | M=100; |        |
|00002da0| 20 20 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 74 |        |    // t|
|00002db0| 6f 20 6d 61 6b 65 20 73 | 75 72 65 20 56 20 77 69 |o make s|ure V wi|
|00002dc0| 6c 6c 20 68 6f 6c 64 20 | 69 74 73 20 76 61 6c 75 |ll hold |its valu|
|00002dd0| 65 73 0a 58 23 65 6e 64 | 69 66 0a 58 20 20 20 66 |es.X#end|if.X   f|
|00002de0| 6f 72 20 28 69 3d 31 3b | 69 3c 3d 35 30 3b 69 2b |or (i=1;|i<=50;i+|
|00002df0| 2b 29 20 56 28 69 29 20 | 2d 3d 20 69 2a 69 2b 69 |+) V(i) |-= i*i+i|
|00002e00| 3b 0a 58 20 20 20 50 72 | 69 6e 74 28 56 29 3b 0a |;.X   Pr|int(V);.|
|00002e10| 58 0a 58 0a 58 2f 2f 20 | 20 20 63 6f 75 74 20 3c |X.X.X// |  cout <|
|00002e20| 3c 20 22 5c 6e 45 6e 64 | 20 6f 66 20 73 65 63 6f |< "\nEnd| of seco|
|00002e30| 6e 64 20 74 65 73 74 5c | 6e 22 3b 0a 58 7d 0a 45 |nd test\|n";.X}.E|
|00002e40| 4e 44 5f 4f 46 5f 46 49 | 4c 45 0a 69 66 20 74 65 |ND_OF_FI|LE.if te|
|00002e50| 73 74 20 32 32 36 34 20 | 2d 6e 65 20 60 77 63 20 |st 2264 |-ne `wc |
|00002e60| 2d 63 20 3c 27 74 6d 74 | 32 2e 63 78 78 27 60 3b |-c <'tmt|2.cxx'`;|
|00002e70| 20 74 68 65 6e 0a 20 20 | 20 20 65 63 68 6f 20 73 | then.  |  echo s|
|00002e80| 68 61 72 3a 20 5c 22 27 | 74 6d 74 32 2e 63 78 78 |har: \"'|tmt2.cxx|
|00002e90| 27 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |'\" unpa|cked wit|
|00002ea0| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 66 69 |h wrong |size!.fi|
|00002eb0| 0a 23 20 65 6e 64 20 6f | 66 20 27 74 6d 74 32 2e |.# end o|f 'tmt2.|
|00002ec0| 63 78 78 27 0a 66 69 0a | 69 66 20 74 65 73 74 20 |cxx'.fi.|if test |
|00002ed0| 2d 66 20 27 74 6d 74 33 | 2e 63 78 78 27 20 2d 61 |-f 'tmt3|.cxx' -a|
|00002ee0| 20 22 24 7b 31 7d 22 20 | 21 3d 20 22 2d 63 22 20 | "${1}" |!= "-c" |
|00002ef0| 3b 20 74 68 65 6e 20 0a | 20 20 65 63 68 6f 20 73 |; then .|  echo s|
|00002f00| 68 61 72 3a 20 57 69 6c | 6c 20 6e 6f 74 20 63 6c |har: Wil|l not cl|
|00002f10| 6f 62 62 65 72 20 65 78 | 69 73 74 69 6e 67 20 66 |obber ex|isting f|
|00002f20| 69 6c 65 20 5c 22 27 74 | 6d 74 33 2e 63 78 78 27 |ile \"'t|mt3.cxx'|
|00002f30| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|00002f40| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|00002f50| 27 74 6d 74 33 2e 63 78 | 78 27 5c 22 20 5c 28 32 |'tmt3.cx|x'\" \(2|
|00002f60| 31 39 32 20 63 68 61 72 | 61 63 74 65 72 73 5c 29 |192 char|acters\)|
|00002f70| 0a 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 27 |.sed "s/|^X//" >'|
|00002f80| 74 6d 74 33 2e 63 78 78 | 27 20 3c 3c 27 45 4e 44 |tmt3.cxx|' <<'END|
|00002f90| 5f 4f 46 5f 46 49 4c 45 | 27 0a 58 0a 58 2f 2f 23 |_OF_FILE|'.X.X//#|
|00002fa0| 64 65 66 69 6e 65 20 57 | 41 4e 54 5f 53 54 52 45 |define W|ANT_STRE|
|00002fb0| 41 4d 0a 58 0a 58 23 69 | 6e 63 6c 75 64 65 20 22 |AM.X.X#i|nclude "|
|00002fc0| 69 6e 63 6c 75 64 65 2e | 68 22 0a 58 0a 58 23 69 |include.|h".X.X#i|
|00002fd0| 6e 63 6c 75 64 65 20 22 | 6e 65 77 6d 61 74 2e 68 |nclude "|newmat.h|
|00002fe0| 22 0a 58 0a 58 0a 58 2f | 2a 2a 2a 2a 2a 2a 2a 2a |".X.X.X/|********|
|00002ff0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003000| 2a 2a 2a 2a 20 74 65 73 | 74 20 70 72 6f 67 72 61 |**** tes|t progra|
|00003010| 6d 20 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |m ******|********|
|00003020| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003030| 2f 0a 58 0a 58 76 6f 69 | 64 20 50 72 69 6e 74 28 |/.X.Xvoi|d Print(|
|00003040| 63 6f 6e 73 74 20 4d 61 | 74 72 69 78 26 20 58 29 |const Ma|trix& X)|
|00003050| 3b 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |;.Xvoid |Print(co|
|00003060| 6e 73 74 20 55 70 70 65 | 72 54 72 69 61 6e 67 75 |nst Uppe|rTriangu|
|00003070| 6c 61 72 4d 61 74 72 69 | 78 26 20 58 29 3b 0a 58 |larMatri|x& X);.X|
|00003080| 76 6f 69 64 20 50 72 69 | 6e 74 28 63 6f 6e 73 74 |void Pri|nt(const|
|00003090| 20 44 69 61 67 6f 6e 61 | 6c 4d 61 74 72 69 78 26 | Diagona|lMatrix&|
|000030a0| 20 58 29 3b 0a 58 76 6f | 69 64 20 50 72 69 6e 74 | X);.Xvo|id Print|
|000030b0| 28 63 6f 6e 73 74 20 53 | 79 6d 6d 65 74 72 69 63 |(const S|ymmetric|
|000030c0| 4d 61 74 72 69 78 26 20 | 58 29 3b 0a 58 76 6f 69 |Matrix& |X);.Xvoi|
|000030d0| 64 20 50 72 69 6e 74 28 | 63 6f 6e 73 74 20 4c 6f |d Print(|const Lo|
|000030e0| 77 65 72 54 72 69 61 6e | 67 75 6c 61 72 4d 61 74 |werTrian|gularMat|
|000030f0| 72 69 78 26 20 58 29 3b | 0a 58 0a 58 76 6f 69 64 |rix& X);|.X.Xvoid|
|00003100| 20 74 72 79 6d 61 74 33 | 28 29 0a 58 7b 0a 58 2f | trymat3|().X{.X/|
|00003110| 2f 20 20 20 63 6f 75 74 | 20 3c 3c 20 22 5c 6e 54 |/   cout| << "\nT|
|00003120| 68 69 72 64 20 74 65 73 | 74 20 6f 66 20 4d 61 74 |hird tes|t of Mat|
|00003130| 72 69 78 20 70 61 63 6b | 61 67 65 5c 6e 22 3b 0a |rix pack|age\n";.|
|00003140| 58 20 20 20 54 72 61 63 | 65 72 20 65 74 28 22 54 |X   Trac|er et("T|
|00003150| 68 69 72 64 20 74 65 73 | 74 20 6f 66 20 4d 61 74 |hird tes|t of Mat|
|00003160| 72 69 78 20 70 61 63 6b | 61 67 65 22 29 3b 0a 58 |rix pack|age");.X|
|00003170| 20 20 20 45 78 63 65 70 | 74 69 6f 6e 3a 3a 50 72 |   Excep|tion::Pr|
|00003180| 69 6e 74 54 72 61 63 65 | 28 54 52 55 45 29 3b 0a |intTrace|(TRUE);.|
|00003190| 58 0a 58 20 20 20 69 6e | 74 20 69 2c 6a 3b 0a 58 |X.X   in|t i,j;.X|
|000031a0| 0a 58 20 20 20 7b 0a 58 | 20 20 20 20 20 20 54 72 |.X   {.X|      Tr|
|000031b0| 61 63 65 72 20 65 74 31 | 28 22 53 74 61 67 65 20 |acer et1|("Stage |
|000031c0| 31 22 29 3b 0a 58 20 20 | 20 20 20 20 53 79 6d 6d |1");.X  |    Symm|
|000031d0| 65 74 72 69 63 4d 61 74 | 72 69 78 20 53 28 37 29 |etricMat|rix S(7)|
|000031e0| 3b 0a 58 20 20 20 20 20 | 20 66 6f 72 20 28 69 3d |;.X     | for (i=|
|000031f0| 31 3b 69 3c 3d 37 3b 69 | 2b 2b 29 20 66 6f 72 20 |1;i<=7;i|++) for |
|00003200| 28 6a 3d 31 3b 6a 3c 3d | 69 3b 6a 2b 2b 29 20 53 |(j=1;j<=|i;j++) S|
|00003210| 28 69 2c 6a 29 3d 69 2a | 69 2b 6a 3b 0a 58 20 20 |(i,j)=i*|i+j;.X  |
|00003220| 20 20 20 20 53 3d 2d 53 | 2b 32 2e 30 3b 0a 58 0a |    S=-S|+2.0;.X.|
|00003230| 58 20 20 20 20 20 20 44 | 69 61 67 6f 6e 61 6c 4d |X      D|iagonalM|
|00003240| 61 74 72 69 78 20 44 28 | 37 29 3b 0a 58 20 20 20 |atrix D(|7);.X   |
|00003250| 20 20 20 66 6f 72 20 28 | 69 3d 31 3b 69 3c 3d 37 |   for (|i=1;i<=7|
|00003260| 3b 69 2b 2b 29 20 44 28 | 69 2c 69 29 3d 53 28 69 |;i++) D(|i,i)=S(i|
|00003270| 2c 69 29 3b 0a 58 0a 58 | 20 20 20 20 20 20 4d 61 |,i);.X.X|      Ma|
|00003280| 74 72 69 78 20 4d 34 28 | 37 2c 37 29 3b 20 7b 20 |trix M4(|7,7); { |
|00003290| 4d 34 3d 44 2b 28 44 2b | 34 2e 30 29 3b 20 4d 34 |M4=D+(D+|4.0); M4|
|000032a0| 3d 4d 34 2d 44 2a 32 2e | 30 3b 20 20 4d 34 3d 4d |=M4-D*2.|0;  M4=M|
|000032b0| 34 2d 34 2e 30 3b 20 50 | 72 69 6e 74 28 4d 34 29 |4-4.0; P|rint(M4)|
|000032c0| 3b 20 7d 0a 58 20 20 20 | 20 20 20 53 79 6d 6d 65 |; }.X   |   Symme|
|000032d0| 74 72 69 63 4d 61 74 72 | 69 78 20 53 32 3d 44 3b |tricMatr|ix S2=D;|
|000032e0| 20 4d 61 74 72 69 78 20 | 4d 32 3d 53 32 3b 20 20 | Matrix |M2=S2;  |
|000032f0| 7b 20 4d 32 3d 2d 44 2b | 4d 32 3b 20 50 72 69 6e |{ M2=-D+|M2; Prin|
|00003300| 74 28 4d 32 29 3b 20 7d | 0a 58 20 20 20 20 20 20 |t(M2); }|.X      |
|00003310| 55 70 70 65 72 54 72 69 | 61 6e 67 75 6c 61 72 4d |UpperTri|angularM|
|00003320| 61 74 72 69 78 20 55 32 | 3d 44 3b 20 7b 20 4d 32 |atrix U2|=D; { M2|
|00003330| 3d 55 32 3b 20 4d 32 3d | 44 2d 4d 32 3b 20 50 72 |=U2; M2=|D-M2; Pr|
|00003340| 69 6e 74 28 4d 32 29 3b | 20 7d 0a 58 20 20 20 20 |int(M2);| }.X    |
|00003350| 20 20 4c 6f 77 65 72 54 | 72 69 61 6e 67 75 6c 61 |  LowerT|riangula|
|00003360| 72 4d 61 74 72 69 78 20 | 4c 32 3d 44 3b 20 7b 20 |rMatrix |L2=D; { |
|00003370| 4d 32 3d 4c 32 3b 20 4d | 32 3d 44 2d 4d 32 3b 20 |M2=L2; M|2=D-M2; |
|00003380| 50 72 69 6e 74 28 4d 32 | 29 3b 20 7d 0a 58 20 20 |Print(M2|); }.X  |
|00003390| 20 20 20 20 4d 32 3d 44 | 3b 20 4d 32 3d 4d 32 2d |    M2=D|; M2=M2-|
|000033a0| 44 3b 20 50 72 69 6e 74 | 28 4d 32 29 3b 0a 58 20 |D; Print|(M2);.X |
|000033b0| 20 20 20 20 20 66 6f 72 | 20 28 69 3d 31 3b 69 3c |     for| (i=1;i<|
|000033c0| 3d 37 3b 69 2b 2b 29 20 | 66 6f 72 20 28 6a 3d 31 |=7;i++) |for (j=1|
|000033d0| 3b 6a 3c 3d 69 3b 6a 2b | 2b 29 20 4c 32 28 69 2c |;j<=i;j+|+) L2(i,|
|000033e0| 6a 29 3d 32 2e 30 2d 69 | 2a 69 2d 6a 3b 0a 58 20 |j)=2.0-i|*i-j;.X |
|000033f0| 20 20 20 20 20 55 32 3d | 4c 32 2e 74 28 29 3b 20 |     U2=|L2.t(); |
|00003400| 44 3d 44 2e 74 28 29 3b | 20 53 3d 53 2e 74 28 29 |D=D.t();| S=S.t()|
|00003410| 3b 0a 58 20 20 20 20 20 | 20 4d 34 3d 28 4c 32 2d |;.X     | M4=(L2-|
|00003420| 31 2e 30 29 2b 28 55 32 | 2b 31 2e 30 29 2d 44 2d |1.0)+(U2|+1.0)-D-|
|00003430| 53 3b 20 50 72 69 6e 74 | 28 4d 34 29 3b 0a 58 20 |S; Print|(M4);.X |
|00003440| 20 20 20 20 20 4d 34 3d | 28 2d 4c 32 2b 31 2e 30 |     M4=|(-L2+1.0|
|00003450| 29 2b 28 2d 55 32 2d 31 | 2e 30 29 2b 44 2b 53 3b |)+(-U2-1|.0)+D+S;|
|00003460| 20 50 72 69 6e 74 28 4d | 34 29 3b 0a 58 20 20 20 | Print(M|4);.X   |
|00003470| 7d 0a 58 0a 58 20 20 20 | 7b 0a 58 20 20 20 20 20 |}.X.X   |{.X     |
|00003480| 20 54 72 61 63 65 72 20 | 65 74 31 28 22 53 74 61 | Tracer |et1("Sta|
|00003490| 67 65 20 32 22 29 3b 0a | 58 20 20 20 20 20 20 44 |ge 2");.|X      D|
|000034a0| 69 61 67 6f 6e 61 6c 4d | 61 74 72 69 78 20 44 28 |iagonalM|atrix D(|
|000034b0| 36 29 3b 0a 58 20 20 20 | 20 20 20 66 6f 72 20 28 |6);.X   |   for (|
|000034c0| 69 3d 31 3b 69 3c 3d 36 | 3b 69 2b 2b 29 20 44 28 |i=1;i<=6|;i++) D(|
|000034d0| 69 2c 69 29 3d 69 2a 33 | 2e 30 2b 69 2a 69 2b 32 |i,i)=i*3|.0+i*i+2|
|000034e0| 2e 30 3b 0a 58 20 20 20 | 20 20 20 55 70 70 65 72 |.0;.X   |   Upper|
|000034f0| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|00003500| 20 55 32 28 37 29 3b 20 | 4c 6f 77 65 72 54 72 69 | U2(7); |LowerTri|
|00003510| 61 6e 67 75 6c 61 72 4d | 61 74 72 69 78 20 4c 32 |angularM|atrix L2|
|00003520| 28 37 29 3b 0a 58 20 20 | 20 20 20 20 66 6f 72 20 |(7);.X  |    for |
|00003530| 28 69 3d 31 3b 69 3c 3d | 37 3b 69 2b 2b 29 20 66 |(i=1;i<=|7;i++) f|
|00003540| 6f 72 20 28 6a 3d 31 3b | 6a 3c 3d 69 3b 6a 2b 2b |or (j=1;|j<=i;j++|
|00003550| 29 20 4c 32 28 69 2c 6a | 29 3d 32 2e 30 2d 69 2a |) L2(i,j|)=2.0-i*|
|00003560| 69 2b 6a 3b 0a 58 20 20 | 20 20 20 20 7b 20 55 32 |i+j;.X  |    { U2|
|00003570| 3d 4c 32 2e 74 28 29 3b | 20 7d 0a 58 20 20 20 20 |=L2.t();| }.X    |
|00003580| 20 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |  Diagon|alMatrix|
|00003590| 20 44 31 28 37 29 3b 20 | 66 6f 72 20 28 69 3d 31 | D1(7); |for (i=1|
|000035a0| 3b 69 3c 3d 37 3b 69 2b | 2b 29 20 44 31 28 69 2c |;i<=7;i+|+) D1(i,|
|000035b0| 69 29 3d 28 69 2d 32 29 | 2a 28 69 2d 34 29 3b 0a |i)=(i-2)|*(i-4);.|
|000035c0| 58 20 20 20 20 20 20 4d | 61 74 72 69 78 20 4d 32 |X      M|atrix M2|
|000035d0| 28 36 2c 37 29 3b 0a 58 | 20 20 20 20 20 20 66 6f |(6,7);.X|      fo|
|000035e0| 72 20 28 69 3d 31 3b 69 | 3c 3d 36 3b 69 2b 2b 29 |r (i=1;i|<=6;i++)|
|000035f0| 20 66 6f 72 20 28 6a 3d | 31 3b 6a 3c 3d 37 3b 6a | for (j=|1;j<=7;j|
|00003600| 2b 2b 29 20 4d 32 28 69 | 2c 6a 29 3d 32 2e 30 2b |++) M2(i|,j)=2.0+|
|00003610| 69 2a 6a 2b 69 2a 69 2b | 32 2e 30 2a 6a 2a 6a 3b |i*j+i*i+|2.0*j*j;|
|00003620| 0a 58 20 20 20 20 20 20 | 4d 61 74 72 69 78 20 4d |.X      |Matrix M|
|00003630| 44 3d 44 3b 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |D=D; Sym|metricMa|
|00003640| 74 72 69 78 20 4d 44 31 | 28 31 29 3b 20 4d 44 31 |trix MD1|(1); MD1|
|00003650| 3d 44 31 3b 0a 58 20 20 | 20 20 20 20 4d 61 74 72 |=D1;.X  |    Matr|
|00003660| 69 78 20 4d 58 3d 4d 44 | 2a 4d 32 2a 4d 44 31 20 |ix MX=MD|*M2*MD1 |
|00003670| 2d 20 44 2a 28 4d 32 2a | 44 31 29 3b 20 50 72 69 |- D*(M2*|D1); Pri|
|00003680| 6e 74 28 4d 58 29 3b 0a | 58 20 20 20 20 20 20 4d |nt(MX);.|X      M|
|00003690| 58 3d 4d 44 2a 4d 32 2a | 4d 44 31 20 2d 20 28 44 |X=MD*M2*|MD1 - (D|
|000036a0| 2a 4d 32 29 2a 44 31 3b | 20 50 72 69 6e 74 28 4d |*M2)*D1;| Print(M|
|000036b0| 58 29 3b 0a 58 20 20 20 | 20 20 20 7b 0a 58 20 20 |X);.X   |   {.X  |
|000036c0| 20 20 20 20 20 20 20 44 | 2e 52 65 44 69 6d 65 6e |       D|.ReDimen|
|000036d0| 73 69 6f 6e 28 37 29 3b | 20 66 6f 72 20 28 69 3d |sion(7);| for (i=|
|000036e0| 31 3b 69 3c 3d 37 3b 69 | 2b 2b 29 20 44 28 69 2c |1;i<=7;i|++) D(i,|
|000036f0| 69 29 3d 69 2a 33 2e 30 | 2b 69 2a 69 2b 32 2e 30 |i)=i*3.0|+i*i+2.0|
|00003700| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 4c 6f 77 65 |;.X     |    Lowe|
|00003710| 72 54 72 69 61 6e 67 75 | 6c 61 72 4d 61 74 72 69 |rTriangu|larMatri|
|00003720| 78 20 4c 44 28 31 29 3b | 20 4c 44 3d 44 3b 0a 58 |x LD(1);| LD=D;.X|
|00003730| 20 20 20 20 20 20 20 20 | 20 55 70 70 65 72 54 72 |        | UpperTr|
|00003740| 69 61 6e 67 75 6c 61 72 | 4d 61 74 72 69 78 20 55 |iangular|Matrix U|
|00003750| 44 28 31 29 3b 20 55 44 | 3d 44 3b 0a 58 20 20 20 |D(1); UD|=D;.X   |
|00003760| 20 20 20 20 20 20 4d 32 | 3d 55 32 3b 20 4d 32 3d |      M2|=U2; M2=|
|00003770| 4c 44 2a 4d 32 2a 4d 44 | 31 20 2d 20 44 2a 28 55 |LD*M2*MD|1 - D*(U|
|00003780| 32 2a 44 31 29 3b 20 50 | 72 69 6e 74 28 4d 32 29 |2*D1); P|rint(M2)|
|00003790| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 4d 32 3d 55 |;.X     |    M2=U|
|000037a0| 32 3b 20 4d 32 3d 55 44 | 2a 4d 32 2a 4d 44 31 20 |2; M2=UD|*M2*MD1 |
|000037b0| 2d 20 28 44 2a 55 32 29 | 2a 44 31 3b 20 50 72 69 |- (D*U2)|*D1; Pri|
|000037c0| 6e 74 28 4d 32 29 3b 0a | 58 20 20 20 20 20 20 20 |nt(M2);.|X       |
|000037d0| 20 20 4d 32 3d 4c 32 3b | 20 4d 32 3d 4c 44 2a 4d |  M2=L2;| M2=LD*M|
|000037e0| 32 2a 4d 44 31 20 2d 20 | 44 2a 28 4c 32 2a 44 31 |2*MD1 - |D*(L2*D1|
|000037f0| 29 3b 20 50 72 69 6e 74 | 28 4d 32 29 3b 0a 58 20 |); Print|(M2);.X |
|00003800| 20 20 20 20 20 20 20 20 | 4d 32 3d 4c 32 3b 20 4d |        |M2=L2; M|
|00003810| 32 3d 55 44 2a 4d 32 2a | 4d 44 31 20 2d 20 28 44 |2=UD*M2*|MD1 - (D|
|00003820| 2a 4c 32 29 2a 44 31 3b | 20 50 72 69 6e 74 28 4d |*L2)*D1;| Print(M|
|00003830| 32 29 3b 0a 58 20 20 20 | 20 20 20 7d 0a 58 20 20 |2);.X   |   }.X  |
|00003840| 20 7d 0a 58 0a 58 2f 2f | 20 20 20 63 6f 75 74 20 | }.X.X//|   cout |
|00003850| 3c 3c 20 22 5c 6e 45 6e | 64 20 6f 66 20 74 68 69 |<< "\nEn|d of thi|
|00003860| 72 64 20 74 65 73 74 5c | 6e 22 3b 0a 58 7d 0a 58 |rd test\|n";.X}.X|
|00003870| 0a 45 4e 44 5f 4f 46 5f | 46 49 4c 45 0a 69 66 20 |.END_OF_|FILE.if |
|00003880| 74 65 73 74 20 32 31 39 | 32 20 2d 6e 65 20 60 77 |test 219|2 -ne `w|
|00003890| 63 20 2d 63 20 3c 27 74 | 6d 74 33 2e 63 78 78 27 |c -c <'t|mt3.cxx'|
|000038a0| 60 3b 20 74 68 65 6e 0a | 20 20 20 20 65 63 68 6f |`; then.|    echo|
|000038b0| 20 73 68 61 72 3a 20 5c | 22 27 74 6d 74 33 2e 63 | shar: \|"'tmt3.c|
|000038c0| 78 78 27 5c 22 20 75 6e | 70 61 63 6b 65 64 20 77 |xx'\" un|packed w|
|000038d0| 69 74 68 20 77 72 6f 6e | 67 20 73 69 7a 65 21 0a |ith wron|g size!.|
|000038e0| 66 69 0a 23 20 65 6e 64 | 20 6f 66 20 27 74 6d 74 |fi.# end| of 'tmt|
|000038f0| 33 2e 63 78 78 27 0a 66 | 69 0a 69 66 20 74 65 73 |3.cxx'.f|i.if tes|
|00003900| 74 20 2d 66 20 27 74 6d | 74 34 2e 63 78 78 27 20 |t -f 'tm|t4.cxx' |
|00003910| 2d 61 20 22 24 7b 31 7d | 22 20 21 3d 20 22 2d 63 |-a "${1}|" != "-c|
|00003920| 22 20 3b 20 74 68 65 6e | 20 0a 20 20 65 63 68 6f |" ; then| .  echo|
|00003930| 20 73 68 61 72 3a 20 57 | 69 6c 6c 20 6e 6f 74 20 | shar: W|ill not |
|00003940| 63 6c 6f 62 62 65 72 20 | 65 78 69 73 74 69 6e 67 |clobber |existing|
|00003950| 20 66 69 6c 65 20 5c 22 | 27 74 6d 74 34 2e 63 78 | file \"|'tmt4.cx|
|00003960| 78 27 5c 22 0a 65 6c 73 | 65 0a 65 63 68 6f 20 73 |x'\".els|e.echo s|
|00003970| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|00003980| 5c 22 27 74 6d 74 34 2e | 63 78 78 27 5c 22 20 5c |\"'tmt4.|cxx'\" \|
|00003990| 28 32 34 32 38 20 63 68 | 61 72 61 63 74 65 72 73 |(2428 ch|aracters|
|000039a0| 5c 29 0a 73 65 64 20 22 | 73 2f 5e 58 2f 2f 22 20 |\).sed "|s/^X//" |
|000039b0| 3e 27 74 6d 74 34 2e 63 | 78 78 27 20 3c 3c 27 45 |>'tmt4.c|xx' <<'E|
|000039c0| 4e 44 5f 4f 46 5f 46 49 | 4c 45 27 0a 58 0a 58 2f |ND_OF_FI|LE'.X.X/|
|000039d0| 2f 23 64 65 66 69 6e 65 | 20 57 41 4e 54 5f 53 54 |/#define| WANT_ST|
|000039e0| 52 45 41 4d 0a 58 0a 58 | 0a 58 23 69 6e 63 6c 75 |REAM.X.X|.X#inclu|
|000039f0| 64 65 20 22 69 6e 63 6c | 75 64 65 2e 68 22 0a 58 |de "incl|ude.h".X|
|00003a00| 0a 58 23 69 6e 63 6c 75 | 64 65 20 22 6e 65 77 6d |.X#inclu|de "newm|
|00003a10| 61 74 2e 68 22 0a 58 0a | 58 0a 58 2f 2a 2a 2a 2a |at.h".X.|X.X/****|
|00003a20| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003a30| 2a 2a 2a 2a 2a 2a 2a 2a | 20 74 65 73 74 20 70 72 |********| test pr|
|00003a40| 6f 67 72 61 6d 20 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |ogram **|********|
|00003a50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003a60| 2a 2a 2a 2a 2f 0a 58 0a | 58 76 6f 69 64 20 50 72 |****/.X.|Xvoid Pr|
|00003a70| 69 6e 74 28 63 6f 6e 73 | 74 20 4d 61 74 72 69 78 |int(cons|t Matrix|
|00003a80| 26 20 58 29 3b 0a 58 76 | 6f 69 64 20 50 72 69 6e |& X);.Xv|oid Prin|
|00003a90| 74 28 63 6f 6e 73 74 20 | 55 70 70 65 72 54 72 69 |t(const |UpperTri|
|00003aa0| 61 6e 67 75 6c 61 72 4d | 61 74 72 69 78 26 20 58 |angularM|atrix& X|
|00003ab0| 29 3b 0a 58 76 6f 69 64 | 20 50 72 69 6e 74 28 63 |);.Xvoid| Print(c|
|00003ac0| 6f 6e 73 74 20 44 69 61 | 67 6f 6e 61 6c 4d 61 74 |onst Dia|gonalMat|
|00003ad0| 72 69 78 26 20 58 29 3b | 0a 58 76 6f 69 64 20 50 |rix& X);|.Xvoid P|
|00003ae0| 72 69 6e 74 28 63 6f 6e | 73 74 20 53 79 6d 6d 65 |rint(con|st Symme|
|00003af0| 74 72 69 63 4d 61 74 72 | 69 78 26 20 58 29 3b 0a |tricMatr|ix& X);.|
|00003b00| 58 76 6f 69 64 20 50 72 | 69 6e 74 28 63 6f 6e 73 |Xvoid Pr|int(cons|
|00003b10| 74 20 4c 6f 77 65 72 54 | 72 69 61 6e 67 75 6c 61 |t LowerT|riangula|
|00003b20| 72 4d 61 74 72 69 78 26 | 20 58 29 3b 0a 58 0a 58 |rMatrix&| X);.X.X|
|00003b30| 76 6f 69 64 20 74 72 79 | 6d 61 74 34 28 29 0a 58 |void try|mat4().X|
|00003b40| 7b 0a 58 2f 2f 20 20 20 | 63 6f 75 74 20 3c 3c 20 |{.X//   |cout << |
|00003b50| 22 5c 6e 46 6f 75 72 74 | 68 20 74 65 73 74 20 6f |"\nFourt|h test o|
|00003b60| 66 20 4d 61 74 72 69 78 | 20 70 61 63 6b 61 67 65 |f Matrix| package|
|00003b70| 5c 6e 22 3b 0a 58 20 20 | 20 54 72 61 63 65 72 20 |\n";.X  | Tracer |
|00003b80| 65 74 28 22 46 6f 75 72 | 74 68 20 74 65 73 74 20 |et("Four|th test |
|00003b90| 6f 66 20 4d 61 74 72 69 | 78 20 70 61 63 6b 61 67 |of Matri|x packag|
|00003ba0| 65 22 29 3b 0a 58 20 20 | 20 45 78 63 65 70 74 69 |e");.X  | Excepti|
|00003bb0| 6f 6e 3a 3a 50 72 69 6e | 74 54 72 61 63 65 28 54 |on::Prin|tTrace(T|
|00003bc0| 52 55 45 29 3b 0a 58 0a | 58 20 20 20 69 6e 74 20 |RUE);.X.|X   int |
|00003bd0| 69 2c 6a 3b 0a 58 0a 58 | 20 20 20 7b 0a 58 20 20 |i,j;.X.X|   {.X  |
|00003be0| 20 20 20 20 54 72 61 63 | 65 72 20 65 74 31 28 22 |    Trac|er et1("|
|00003bf0| 53 74 61 67 65 20 31 22 | 29 3b 0a 58 20 20 20 20 |Stage 1"|);.X    |
|00003c00| 20 20 4d 61 74 72 69 78 | 20 4d 28 31 30 2c 31 30 |  Matrix| M(10,10|
|00003c10| 29 3b 0a 58 20 20 20 20 | 20 20 55 70 70 65 72 54 |);.X    |  UpperT|
|00003c20| 72 69 61 6e 67 75 6c 61 | 72 4d 61 74 72 69 78 20 |riangula|rMatrix |
|00003c30| 55 28 31 30 29 3b 0a 58 | 20 20 20 20 20 20 66 6f |U(10);.X|      fo|
|00003c40| 72 20 28 69 3d 31 3b 69 | 3c 3d 31 30 3b 69 2b 2b |r (i=1;i|<=10;i++|
|00003c50| 29 20 66 6f 72 20 28 6a | 3d 31 3b 6a 3c 3d 31 30 |) for (j|=1;j<=10|
|00003c60| 3b 6a 2b 2b 29 20 4d 28 | 69 2c 6a 29 20 3d 20 31 |;j++) M(|i,j) = 1|
|00003c70| 30 30 2a 69 2b 6a 3b 0a | 58 20 20 20 20 20 20 55 |00*i+j;.|X      U|
|00003c80| 20 3c 3c 20 2d 4d 3b 0a | 58 20 20 20 20 20 20 4d | << -M;.|X      M|
|00003c90| 61 74 72 69 78 20 58 31 | 20 3d 20 4d 2e 52 6f 77 |atrix X1| = M.Row|
|00003ca0| 73 28 32 2c 34 29 3b 0a | 58 20 20 20 20 20 20 4d |s(2,4);.|X      M|
|00003cb0| 61 74 72 69 78 20 59 31 | 20 3d 20 55 2e 74 28 29 |atrix Y1| = U.t()|
|00003cc0| 2e 52 6f 77 73 28 32 2c | 34 29 3b 0a 58 20 20 20 |.Rows(2,|4);.X   |
|00003cd0| 20 20 20 4d 61 74 72 69 | 78 20 58 20 3d 20 55 3b |   Matri|x X = U;|
|00003ce0| 20 7b 20 50 72 69 6e 74 | 28 4d 61 74 72 69 78 28 | { Print|(Matrix(|
|00003cf0| 58 2e 43 6f 6c 75 6d 6e | 73 28 32 2c 34 29 2e 74 |X.Column|s(2,4).t|
|00003d00| 28 29 2d 59 31 29 29 3b | 20 7d 0a 58 20 20 20 20 |()-Y1));| }.X    |
|00003d10| 20 20 52 6f 77 56 65 63 | 74 6f 72 20 52 56 20 3d |  RowVec|tor RV =|
|00003d20| 20 4d 2e 52 6f 77 28 35 | 29 3b 0a 58 20 20 20 20 | M.Row(5|);.X    |
|00003d30| 20 20 7b 0a 58 20 20 20 | 20 20 20 20 20 20 58 2e |  {.X   |      X.|
|00003d40| 52 65 44 69 6d 65 6e 73 | 69 6f 6e 28 33 2c 31 30 |ReDimens|ion(3,10|
|00003d50| 29 3b 0a 58 20 20 20 20 | 20 20 20 20 20 58 2e 52 |);.X    |     X.R|
|00003d60| 6f 77 28 31 29 20 3c 3c | 20 4d 2e 52 6f 77 28 32 |ow(1) <<| M.Row(2|
|00003d70| 29 3b 20 58 2e 52 6f 77 | 28 32 29 20 3c 3c 20 4d |); X.Row|(2) << M|
|00003d80| 2e 52 6f 77 28 33 29 3b | 20 58 2e 52 6f 77 28 33 |.Row(3);| X.Row(3|
|00003d90| 29 20 3c 3c 20 4d 2e 52 | 6f 77 28 34 29 3b 0a 58 |) << M.R|ow(4);.X|
|00003da0| 20 20 20 20 20 20 20 20 | 20 50 72 69 6e 74 28 4d |        | Print(M|
|00003db0| 61 74 72 69 78 28 58 2d | 58 31 29 29 3b 0a 58 20 |atrix(X-|X1));.X |
|00003dc0| 20 20 20 20 20 7d 0a 58 | 20 20 20 20 20 20 7b 0a |     }.X|      {.|
|00003dd0| 58 20 20 20 20 20 20 20 | 20 20 55 70 70 65 72 54 |X       |  UpperT|
|00003de0| 72 69 61 6e 67 75 6c 61 | 72 4d 61 74 72 69 78 20 |riangula|rMatrix |
|00003df0| 56 20 3d 20 55 2e 53 79 | 6d 53 75 62 4d 61 74 72 |V = U.Sy|mSubMatr|
|00003e00| 69 78 28 33 2c 35 29 3b | 0a 58 20 20 20 20 20 20 |ix(3,5);|.X      |
|00003e10| 20 20 20 4d 61 74 72 69 | 78 20 4d 56 20 3d 20 55 |   Matri|x MV = U|
|00003e20| 2e 53 75 62 4d 61 74 72 | 69 78 28 33 2c 35 2c 33 |.SubMatr|ix(3,5,3|
|00003e30| 2c 35 29 3b 20 7b 20 50 | 72 69 6e 74 28 4d 61 74 |,5); { P|rint(Mat|
|00003e40| 72 69 78 28 4d 56 2d 56 | 29 29 3b 20 7d 0a 58 20 |rix(MV-V|)); }.X |
|00003e50| 20 20 20 20 20 20 20 20 | 4d 61 74 72 69 78 20 58 |        |Matrix X|
|00003e60| 32 20 3d 20 4d 2e 74 28 | 29 2e 43 6f 6c 75 6d 6e |2 = M.t(|).Column|
|00003e70| 73 28 32 2c 34 29 3b 20 | 7b 20 50 72 69 6e 74 28 |s(2,4); |{ Print(|
|00003e80| 4d 61 74 72 69 78 28 58 | 32 2d 58 31 2e 74 28 29 |Matrix(X|2-X1.t()|
|00003e90| 29 29 3b 20 7d 0a 58 20 | 20 20 20 20 20 20 20 20 |)); }.X |        |
|00003ea0| 4d 61 74 72 69 78 20 59 | 32 20 3d 20 55 2e 43 6f |Matrix Y|2 = U.Co|
|00003eb0| 6c 75 6d 6e 73 28 32 2c | 34 29 3b 20 7b 20 50 72 |lumns(2,|4); { Pr|
|00003ec0| 69 6e 74 28 4d 61 74 72 | 69 78 28 59 32 2d 59 31 |int(Matr|ix(Y2-Y1|
|00003ed0| 2e 74 28 29 29 29 3b 20 | 7d 0a 58 20 20 20 20 20 |.t())); |}.X     |
|00003ee0| 20 20 20 20 43 6f 6c 75 | 6d 6e 56 65 63 74 6f 72 |    Colu|mnVector|
|00003ef0| 20 43 56 20 3d 20 4d 2e | 74 28 29 2e 43 6f 6c 75 | CV = M.|t().Colu|
|00003f00| 6d 6e 28 35 29 3b 20 7b | 20 50 72 69 6e 74 28 43 |mn(5); {| Print(C|
|00003f10| 6f 6c 75 6d 6e 56 65 63 | 74 6f 72 28 43 56 2d 52 |olumnVec|tor(CV-R|
|00003f20| 56 2e 74 28 29 29 29 3b | 20 7d 0a 58 20 20 20 20 |V.t()));| }.X    |
|00003f30| 20 20 20 20 20 58 2e 52 | 65 44 69 6d 65 6e 73 69 |     X.R|eDimensi|
|00003f40| 6f 6e 28 31 30 2c 33 29 | 3b 20 4d 20 3d 20 4d 2e |on(10,3)|; M = M.|
|00003f50| 74 28 29 3b 0a 58 20 20 | 20 20 20 20 20 20 20 58 |t();.X  |       X|
|00003f60| 2e 43 6f 6c 75 6d 6e 28 | 31 29 20 3c 3c 20 4d 2e |.Column(|1) << M.|
|00003f70| 43 6f 6c 75 6d 6e 28 32 | 29 3b 20 58 2e 43 6f 6c |Column(2|); X.Col|
|00003f80| 75 6d 6e 28 32 29 20 3c | 3c 20 4d 2e 43 6f 6c 75 |umn(2) <|< M.Colu|
|00003f90| 6d 6e 28 33 29 3b 0a 58 | 20 20 20 20 20 20 20 20 |mn(3);.X|        |
|00003fa0| 20 58 2e 43 6f 6c 75 6d | 6e 28 33 29 20 3c 3c 20 | X.Colum|n(3) << |
|00003fb0| 4d 2e 43 6f 6c 75 6d 6e | 28 34 29 3b 0a 58 20 20 |M.Column|(4);.X  |
|00003fc0| 20 20 20 20 20 20 20 50 | 72 69 6e 74 28 4d 61 74 |       P|rint(Mat|
|00003fd0| 72 69 78 28 58 2d 58 32 | 29 29 3b 0a 58 20 20 20 |rix(X-X2|));.X   |
|00003fe0| 20 20 20 7d 0a 58 20 20 | 20 7d 0a 58 0a 58 20 20 |   }.X  | }.X.X  |
|00003ff0| 20 7b 0a 58 20 20 20 20 | 20 20 54 72 61 63 65 72 | {.X    |  Tracer|
|00004000| 20 65 74 31 28 22 53 74 | 61 67 65 20 32 22 29 3b | et1("St|age 2");|
|00004010| 0a 58 20 20 20 20 20 20 | 4d 61 74 72 69 78 20 4d |.X      |Matrix M|
|00004020| 3b 20 4d 61 74 72 69 78 | 20 58 3b 20 4d 2e 52 65 |; Matrix| X; M.Re|
|00004030| 44 69 6d 65 6e 73 69 6f | 6e 28 35 2c 38 29 3b 0a |Dimensio|n(5,8);.|
|00004040| 58 20 20 20 20 20 20 66 | 6f 72 20 28 69 3d 31 3b |X      f|or (i=1;|
|00004050| 69 3c 3d 35 3b 69 2b 2b | 29 20 66 6f 72 20 28 6a |i<=5;i++|) for (j|
|00004060| 3d 31 3b 6a 3c 3d 38 3b | 6a 2b 2b 29 20 4d 28 69 |=1;j<=8;|j++) M(i|
|00004070| 2c 6a 29 20 3d 20 31 30 | 30 2a 69 2b 6a 3b 0a 58 |,j) = 10|0*i+j;.X|
|00004080| 20 20 20 20 20 20 7b 0a | 58 20 20 20 20 20 20 20 |      {.|X       |
|00004090| 20 20 58 20 3d 20 4d 2e | 43 6f 6c 75 6d 6e 73 28 |  X = M.|Columns(|
|000040a0| 35 2c 38 29 3b 20 4d 2e | 43 6f 6c 75 6d 6e 73 28 |5,8); M.|Columns(|
|000040b0| 35 2c 38 29 20 3c 3c 20 | 4d 2e 43 6f 6c 75 6d 6e |5,8) << |M.Column|
|000040c0| 73 28 31 2c 34 29 3b 0a | 58 20 20 20 20 20 20 20 |s(1,4);.|X       |
|000040d0| 20 20 20 20 20 20 4d 2e | 43 6f 6c 75 6d 6e 73 28 |      M.|Columns(|
|000040e0| 31 2c 34 29 20 3c 3c 20 | 58 3b 0a 58 20 20 20 20 |1,4) << |X;.X    |
|000040f0| 20 20 20 20 20 58 20 3d | 20 4d 2e 43 6f 6c 75 6d |     X =| M.Colum|
|00004100| 6e 73 28 33 2c 34 29 3b | 20 4d 2e 43 6f 6c 75 6d |ns(3,4);| M.Colum|
|00004110| 6e 73 28 33 2c 34 29 20 | 3c 3c 20 4d 2e 43 6f 6c |ns(3,4) |<< M.Col|
|00004120| 75 6d 6e 73 28 31 2c 32 | 29 3b 0a 58 20 20 20 20 |umns(1,2|);.X    |
|00004130| 20 20 20 20 20 20 20 20 | 20 4d 2e 43 6f 6c 75 6d |        | M.Colum|
|00004140| 6e 73 28 31 2c 32 29 20 | 3c 3c 20 58 3b 0a 58 20 |ns(1,2) |<< X;.X |
|00004150| 20 20 20 20 20 20 20 20 | 58 20 3d 20 4d 2e 43 6f |        |X = M.Co|
|00004160| 6c 75 6d 6e 73 28 37 2c | 38 29 3b 20 4d 2e 43 6f |lumns(7,|8); M.Co|
|00004170| 6c 75 6d 6e 73 28 37 2c | 38 29 20 3c 3c 20 4d 2e |lumns(7,|8) << M.|
|00004180| 43 6f 6c 75 6d 6e 73 28 | 35 2c 36 29 3b 0a 58 20 |Columns(|5,6);.X |
|00004190| 20 20 20 20 20 20 20 20 | 20 20 20 20 4d 2e 43 6f |        |    M.Co|
|000041a0| 6c 75 6d 6e 73 28 35 2c | 36 29 20 3c 3c 20 58 3b |lumns(5,|6) << X;|
|000041b0| 0a 58 20 20 20 20 20 20 | 7d 0a 58 20 20 20 20 20 |.X      |}.X     |
|000041c0| 20 7b 0a 58 20 20 20 20 | 20 20 20 20 20 58 20 3d | {.X    |     X =|
|000041d0| 20 4d 2e 43 6f 6c 75 6d | 6e 28 32 29 3b 20 4d 2e | M.Colum|n(2); M.|
|000041e0| 43 6f 6c 75 6d 6e 28 32 | 29 20 3d 20 4d 2e 43 6f |Column(2|) = M.Co|
|000041f0| 6c 75 6d 6e 28 31 29 3b | 20 4d 2e 43 6f 6c 75 6d |lumn(1);| M.Colum|
|00004200| 6e 28 31 29 20 3d 20 58 | 3b 0a 58 20 20 20 20 20 |n(1) = X|;.X     |
|00004210| 20 20 20 20 58 20 3d 20 | 4d 2e 43 6f 6c 75 6d 6e |    X = |M.Column|
|00004220| 28 34 29 3b 20 4d 2e 43 | 6f 6c 75 6d 6e 28 34 29 |(4); M.C|olumn(4)|
|00004230| 20 3d 20 4d 2e 43 6f 6c | 75 6d 6e 28 33 29 3b 20 | = M.Col|umn(3); |
|00004240| 4d 2e 43 6f 6c 75 6d 6e | 28 33 29 20 3d 20 58 3b |M.Column|(3) = X;|
|00004250| 0a 58 20 20 20 20 20 20 | 20 20 20 58 20 3d 20 4d |.X      |   X = M|
|00004260| 2e 43 6f 6c 75 6d 6e 28 | 36 29 3b 20 4d 2e 43 6f |.Column(|6); M.Co|
|00004270| 6c 75 6d 6e 28 36 29 20 | 3d 20 4d 2e 43 6f 6c 75 |lumn(6) |= M.Colu|
|00004280| 6d 6e 28 35 29 3b 20 4d | 2e 43 6f 6c 75 6d 6e 28 |mn(5); M|.Column(|
|00004290| 35 29 20 3d 20 58 3b 0a | 58 20 20 20 20 20 20 20 |5) = X;.|X       |
|000042a0| 20 20 58 20 3d 20 4d 2e | 43 6f 6c 75 6d 6e 28 38 |  X = M.|Column(8|
|000042b0| 29 3b 20 4d 2e 43 6f 6c | 75 6d 6e 28 38 29 20 3d |); M.Col|umn(8) =|
|000042c0| 20 4d 2e 43 6f 6c 75 6d | 6e 28 37 29 3b 20 4d 2e | M.Colum|n(7); M.|
|000042d0| 43 6f 6c 75 6d 6e 28 37 | 29 20 3d 20 58 3b 0a 58 |Column(7|) = X;.X|
|000042e0| 20 20 20 20 20 20 20 20 | 20 58 2e 52 65 44 69 6d |        | X.ReDim|
|000042f0| 65 6e 73 69 6f 6e 28 35 | 2c 38 29 3b 0a 58 20 20 |ension(5|,8);.X  |
|00004300| 20 20 20 20 7d 0a 58 20 | 20 20 20 20 20 66 6f 72 |    }.X |     for|
|00004310| 20 28 69 3d 31 3b 69 3c | 3d 35 3b 69 2b 2b 29 20 | (i=1;i<|=5;i++) |
|00004320| 66 6f 72 20 28 6a 3d 31 | 3b 6a 3c 3d 38 3b 6a 2b |for (j=1|;j<=8;j+|
|00004330| 2b 29 20 58 28 69 2c 39 | 2d 6a 29 20 3d 20 31 30 |+) X(i,9|-j) = 10|
|00004340| 30 2a 69 2b 6a 3b 0a 58 | 20 20 20 20 20 20 50 72 |0*i+j;.X|      Pr|
|00004350| 69 6e 74 28 4d 61 74 72 | 69 78 28 58 2d 4d 29 29 |int(Matr|ix(X-M))|
|00004360| 3b 0a 58 20 20 20 7d 0a | 58 0a 58 2f 2f 20 20 20 |;.X   }.|X.X//   |
|00004370| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 45 6e 64 20 6f |cout << |"\nEnd o|
|00004380| 66 20 66 6f 75 72 74 68 | 20 74 65 73 74 5c 6e 22 |f fourth| test\n"|
|00004390| 3b 0a 58 7d 0a 58 0a 45 | 4e 44 5f 4f 46 5f 46 49 |;.X}.X.E|ND_OF_FI|
|000043a0| 4c 45 0a 69 66 20 74 65 | 73 74 20 32 34 32 38 20 |LE.if te|st 2428 |
|000043b0| 2d 6e 65 20 60 77 63 20 | 2d 63 20 3c 27 74 6d 74 |-ne `wc |-c <'tmt|
|000043c0| 34 2e 63 78 78 27 60 3b | 20 74 68 65 6e 0a 20 20 |4.cxx'`;| then.  |
|000043d0| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 27 |  echo s|har: \"'|
|000043e0| 74 6d 74 34 2e 63 78 78 | 27 5c 22 20 75 6e 70 61 |tmt4.cxx|'\" unpa|
|000043f0| 63 6b 65 64 20 77 69 74 | 68 20 77 72 6f 6e 67 20 |cked wit|h wrong |
|00004400| 73 69 7a 65 21 0a 66 69 | 0a 23 20 65 6e 64 20 6f |size!.fi|.# end o|
|00004410| 66 20 27 74 6d 74 34 2e | 63 78 78 27 0a 66 69 0a |f 'tmt4.|cxx'.fi.|
|00004420| 69 66 20 74 65 73 74 20 | 2d 66 20 27 74 6d 74 35 |if test |-f 'tmt5|
|00004430| 2e 63 78 78 27 20 2d 61 | 20 22 24 7b 31 7d 22 20 |.cxx' -a| "${1}" |
|00004440| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|00004450| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c |  echo s|har: Wil|
|00004460| 6c 20 6e 6f 74 20 63 6c | 6f 62 62 65 72 20 65 78 |l not cl|obber ex|
|00004470| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 27 74 |isting f|ile \"'t|
|00004480| 6d 74 35 2e 63 78 78 27 | 5c 22 0a 65 6c 73 65 0a |mt5.cxx'|\".else.|
|00004490| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|000044a0| 63 74 69 6e 67 20 5c 22 | 27 74 6d 74 35 2e 63 78 |cting \"|'tmt5.cx|
|000044b0| 78 27 5c 22 20 5c 28 31 | 37 32 35 20 63 68 61 72 |x'\" \(1|725 char|
|000044c0| 61 63 74 65 72 73 5c 29 | 0a 73 65 64 20 22 73 2f |acters\)|.sed "s/|
|000044d0| 5e 58 2f 2f 22 20 3e 27 | 74 6d 74 35 2e 63 78 78 |^X//" >'|tmt5.cxx|
|000044e0| 27 20 3c 3c 27 45 4e 44 | 5f 4f 46 5f 46 49 4c 45 |' <<'END|_OF_FILE|
|000044f0| 27 0a 58 0a 58 2f 2f 23 | 64 65 66 69 6e 65 20 57 |'.X.X//#|define W|
|00004500| 41 4e 54 5f 53 54 52 45 | 41 4d 0a 58 0a 58 0a 58 |ANT_STRE|AM.X.X.X|
|00004510| 23 69 6e 63 6c 75 64 65 | 20 22 69 6e 63 6c 75 64 |#include| "includ|
|00004520| 65 2e 68 22 0a 58 0a 58 | 23 69 6e 63 6c 75 64 65 |e.h".X.X|#include|
|00004530| 20 22 6e 65 77 6d 61 74 | 2e 68 22 0a 58 0a 58 0a | "newmat|.h".X.X.|
|00004540| 58 2f 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |X/******|********|
|00004550| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 20 74 |********|****** t|
|00004560| 65 73 74 20 70 72 6f 67 | 72 61 6d 20 2a 2a 2a 2a |est prog|ram ****|
|00004570| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004580| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2f 0a 58 0a 58 76 |********|**/.X.Xv|
|00004590| 6f 69 64 20 50 72 69 6e | 74 28 63 6f 6e 73 74 20 |oid Prin|t(const |
|000045a0| 4d 61 74 72 69 78 26 20 | 58 29 3b 0a 58 76 6f 69 |Matrix& |X);.Xvoi|
|000045b0| 64 20 50 72 69 6e 74 28 | 63 6f 6e 73 74 20 55 70 |d Print(|const Up|
|000045c0| 70 65 72 54 72 69 61 6e | 67 75 6c 61 72 4d 61 74 |perTrian|gularMat|
|000045d0| 72 69 78 26 20 58 29 3b | 0a 58 76 6f 69 64 20 50 |rix& X);|.Xvoid P|
|000045e0| 72 69 6e 74 28 63 6f 6e | 73 74 20 44 69 61 67 6f |rint(con|st Diago|
|000045f0| 6e 61 6c 4d 61 74 72 69 | 78 26 20 58 29 3b 0a 58 |nalMatri|x& X);.X|
|00004600| 76 6f 69 64 20 50 72 69 | 6e 74 28 63 6f 6e 73 74 |void Pri|nt(const|
|00004610| 20 53 79 6d 6d 65 74 72 | 69 63 4d 61 74 72 69 78 | Symmetr|icMatrix|
|00004620| 26 20 58 29 3b 0a 58 76 | 6f 69 64 20 50 72 69 6e |& X);.Xv|oid Prin|
|00004630| 74 28 63 6f 6e 73 74 20 | 4c 6f 77 65 72 54 72 69 |t(const |LowerTri|
|00004640| 61 6e 67 75 6c 61 72 4d | 61 74 72 69 78 26 20 58 |angularM|atrix& X|
|00004650| 29 3b 0a 58 0a 58 76 6f | 69 64 20 74 72 79 6d 61 |);.X.Xvo|id tryma|
|00004660| 74 35 28 29 0a 58 7b 0a | 58 2f 2f 20 20 20 63 6f |t5().X{.|X//   co|
|00004670| 75 74 20 3c 3c 20 22 5c | 6e 46 69 66 74 68 20 74 |ut << "\|nFifth t|
|00004680| 65 73 74 20 6f 66 20 4d | 61 74 72 69 78 20 70 61 |est of M|atrix pa|
|00004690| 63 6b 61 67 65 5c 6e 22 | 3b 0a 58 20 20 20 54 72 |ckage\n"|;.X   Tr|
|000046a0| 61 63 65 72 20 65 74 28 | 22 46 69 66 74 68 20 74 |acer et(|"Fifth t|
|000046b0| 65 73 74 20 6f 66 20 4d | 61 74 72 69 78 20 70 61 |est of M|atrix pa|
|000046c0| 63 6b 61 67 65 22 29 3b | 0a 58 20 20 20 45 78 63 |ckage");|.X   Exc|
|000046d0| 65 70 74 69 6f 6e 3a 3a | 50 72 69 6e 74 54 72 61 |eption::|PrintTra|
|000046e0| 63 65 28 54 52 55 45 29 | 3b 0a 58 0a 58 20 20 20 |ce(TRUE)|;.X.X   |
|000046f0| 69 6e 74 20 69 2c 6a 3b | 0a 58 0a 58 20 20 20 4d |int i,j;|.X.X   M|
|00004700| 61 74 72 69 78 20 41 28 | 35 2c 36 29 3b 0a 58 20 |atrix A(|5,6);.X |
|00004710| 20 20 66 6f 72 20 28 69 | 3d 31 3b 69 3c 3d 35 3b |  for (i|=1;i<=5;|
|00004720| 69 2b 2b 29 20 66 6f 72 | 20 28 6a 3d 31 3b 6a 3c |i++) for| (j=1;j<|
|00004730| 3d 36 3b 6a 2b 2b 29 20 | 41 28 69 2c 6a 29 3d 31 |=6;j++) |A(i,j)=1|
|00004740| 2b 69 2a 6a 2b 69 2a 69 | 2b 6a 2a 6a 3b 0a 58 20 |+i*j+i*i|+j*j;.X |
|00004750| 20 20 43 6f 6c 75 6d 6e | 56 65 63 74 6f 72 20 43 |  Column|Vector C|
|00004760| 56 28 36 29 3b 0a 58 20 | 20 20 66 6f 72 20 28 69 |V(6);.X |  for (i|
|00004770| 3d 31 3b 69 3c 3d 36 3b | 69 2b 2b 29 20 43 56 28 |=1;i<=6;|i++) CV(|
|00004780| 69 29 3d 69 2a 69 2b 33 | 3b 0a 58 20 20 20 43 6f |i)=i*i+3|;.X   Co|
|00004790| 6c 75 6d 6e 56 65 63 74 | 6f 72 20 43 56 32 28 35 |lumnVect|or CV2(5|
|000047a0| 29 3b 20 66 6f 72 20 28 | 69 3d 31 3b 69 3c 3d 35 |); for (|i=1;i<=5|
|000047b0| 3b 69 2b 2b 29 20 43 56 | 32 28 69 29 3d 31 2e 30 |;i++) CV|2(i)=1.0|
|000047c0| 3b 0a 58 20 20 20 43 6f | 6c 75 6d 6e 56 65 63 74 |;.X   Co|lumnVect|
|000047d0| 6f 72 20 43 56 31 3d 43 | 56 3b 0a 58 0a 58 20 20 |or CV1=C|V;.X.X  |
|000047e0| 20 7b 0a 58 20 20 20 20 | 20 20 43 56 3d 41 2a 43 | {.X    |  CV=A*C|
|000047f0| 56 3b 0a 58 20 20 20 20 | 20 20 52 6f 77 56 65 63 |V;.X    |  RowVec|
|00004800| 74 6f 72 20 52 56 3d 43 | 56 2e 74 28 29 3b 20 2f |tor RV=C|V.t(); /|
|00004810| 2f 20 52 6f 77 56 65 63 | 74 6f 72 20 52 56 3b 20 |/ RowVec|tor RV; |
|00004820| 52 56 3d 43 56 2e 74 28 | 29 3b 0a 58 20 20 20 20 |RV=CV.t(|);.X    |
|00004830| 20 20 52 56 3d 52 56 2d | 31 2e 30 3b 0a 58 20 20 |  RV=RV-|1.0;.X  |
|00004840| 20 20 20 20 43 56 3d 28 | 52 56 2a 41 29 2e 74 28 |    CV=(|RV*A).t(|
|00004850| 29 2b 41 2e 74 28 29 2a | 43 56 32 3b 20 43 56 31 |)+A.t()*|CV2; CV1|
|00004860| 3d 28 41 2e 74 28 29 2a | 41 29 2a 43 56 31 20 2d |=(A.t()*|A)*CV1 -|
|00004870| 20 43 56 3b 0a 58 20 20 | 20 20 20 20 50 72 69 6e | CV;.X  |    Prin|
|00004880| 74 28 43 56 31 29 3b 0a | 58 20 20 20 7d 0a 58 0a |t(CV1);.|X   }.X.|
|00004890| 58 20 20 20 43 56 31 2e | 52 65 44 69 6d 65 6e 73 |X   CV1.|ReDimens|
|000048a0| 69 6f 6e 28 36 29 3b 0a | 58 20 20 20 43 56 32 2e |ion(6);.|X   CV2.|
|000048b0| 52 65 44 69 6d 65 6e 73 | 69 6f 6e 28 36 29 3b 0a |ReDimens|ion(6);.|
|000048c0| 58 20 20 20 43 56 2e 52 | 65 44 69 6d 65 6e 73 69 |X   CV.R|eDimensi|
|000048d0| 6f 6e 28 36 29 3b 0a 58 | 20 20 20 66 6f 72 20 28 |on(6);.X|   for (|
|000048e0| 69 3d 31 3b 69 3c 3d 36 | 3b 69 2b 2b 29 20 7b 20 |i=1;i<=6|;i++) { |
|000048f0| 43 56 31 28 69 29 3d 69 | 2a 33 2b 31 3b 20 43 56 |CV1(i)=i|*3+1; CV|
|00004900| 32 28 69 29 3d 31 30 2d | 69 3b 20 43 56 28 69 29 |2(i)=10-|i; CV(i)|
|00004910| 3d 31 31 2b 69 2a 32 3b | 20 7d 0a 58 20 20 20 43 |=11+i*2;| }.X   C|
|00004920| 6f 6c 75 6d 6e 56 65 63 | 74 6f 72 20 43 58 3d 43 |olumnVec|tor CX=C|
|00004930| 56 32 2d 43 56 3b 20 7b | 20 43 58 3d 43 58 2b 43 |V2-CV; {| CX=CX+C|
|00004940| 56 31 3b 20 50 72 69 6e | 74 28 43 58 29 3b 20 7d |V1; Prin|t(CX); }|
|00004950| 0a 58 20 20 20 50 72 69 | 6e 74 28 43 6f 6c 75 6d |.X   Pri|nt(Colum|
|00004960| 6e 56 65 63 74 6f 72 28 | 43 56 31 2b 43 56 32 2d |nVector(|CV1+CV2-|
|00004970| 43 56 29 29 3b 0a 58 20 | 20 20 52 6f 77 56 65 63 |CV));.X |  RowVec|
|00004980| 74 6f 72 20 52 56 3d 43 | 56 2e 74 28 29 3b 20 52 |tor RV=C|V.t(); R|
|00004990| 6f 77 56 65 63 74 6f 72 | 20 52 56 31 3d 43 56 31 |owVector| RV1=CV1|
|000049a0| 2e 74 28 29 3b 0a 58 20 | 20 20 52 6f 77 56 65 63 |.t();.X |  RowVec|
|000049b0| 74 6f 72 20 52 3d 52 56 | 2d 52 56 31 3b 20 50 72 |tor R=RV|-RV1; Pr|
|000049c0| 69 6e 74 28 52 6f 77 56 | 65 63 74 6f 72 28 52 2d |int(RowV|ector(R-|
|000049d0| 43 56 32 2e 74 28 29 29 | 29 3b 0a 58 0a 58 2f 2f |CV2.t())|);.X.X//|
|000049e0| 20 74 65 73 74 20 6c 6f | 61 64 69 6e 67 20 6f 66 | test lo|ading of|
|000049f0| 20 6c 69 73 74 0a 58 0a | 58 20 20 20 52 56 2e 52 | list.X.|X   RV.R|
|00004a00| 65 44 69 6d 65 6e 73 69 | 6f 6e 28 31 30 29 3b 0a |eDimensi|on(10);.|
|00004a10| 58 20 20 20 66 6f 72 20 | 28 69 3d 31 3b 69 3c 3d |X   for |(i=1;i<=|
|00004a20| 31 30 3b 69 2b 2b 29 20 | 52 56 28 69 29 20 3d 20 |10;i++) |RV(i) = |
|00004a30| 69 2a 69 3b 0a 58 20 20 | 20 52 56 31 2e 52 65 44 |i*i;.X  | RV1.ReD|
|00004a40| 69 6d 65 6e 73 69 6f 6e | 28 31 30 29 3b 0a 58 20 |imension|(10);.X |
|00004a50| 20 20 52 56 31 20 3c 3c | 20 31 20 3c 3c 20 34 20 |  RV1 <<| 1 << 4 |
|00004a60| 3c 3c 20 39 20 3c 3c 20 | 31 36 20 3c 3c 20 32 35 |<< 9 << |16 << 25|
|00004a70| 20 3c 3c 20 33 36 20 3c | 3c 20 34 39 20 3c 3c 20 | << 36 <|< 49 << |
|00004a80| 36 34 20 3c 3c 20 38 31 | 20 3c 3c 20 31 30 30 3b |64 << 81| << 100;|
|00004a90| 20 2f 2f 20 3c 3c 20 31 | 32 31 3b 0a 58 20 20 20 | // << 1|21;.X   |
|00004aa0| 50 72 69 6e 74 28 52 6f | 77 56 65 63 74 6f 72 28 |Print(Ro|wVector(|
|00004ab0| 52 56 2d 52 56 31 29 29 | 3b 0a 58 20 20 20 65 74 |RV-RV1))|;.X   et|
|00004ac0| 2e 52 65 4e 61 6d 65 28 | 22 46 69 66 74 68 20 74 |.ReName(|"Fifth t|
|00004ad0| 65 73 74 20 6f 66 20 4d | 61 74 72 69 78 20 70 61 |est of M|atrix pa|
|00004ae0| 63 6b 61 67 65 20 2d 20 | 61 74 20 65 6e 64 22 29 |ckage - |at end")|
|00004af0| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 20 58 28 32 |;.X   Ma|trix X(2|
|00004b00| 2c 33 29 3b 0a 58 20 20 | 20 58 20 3c 3c 20 31 31 |,3);.X  | X << 11|
|00004b10| 20 3c 3c 20 31 32 20 3c | 3c 20 31 33 0a 58 20 20 | << 12 <|< 13.X  |
|00004b20| 20 20 20 3c 3c 20 32 31 | 20 3c 3c 20 32 32 20 3c |   << 21| << 22 <|
|00004b30| 3c 20 32 33 3b 0a 58 20 | 20 20 58 28 31 2c 31 29 |< 23;.X |  X(1,1)|
|00004b40| 20 2d 3d 20 31 31 3b 20 | 58 28 31 2c 32 29 20 2d | -= 11; |X(1,2) -|
|00004b50| 3d 20 31 32 3b 20 58 28 | 31 2c 33 29 20 2d 3d 20 |= 12; X(|1,3) -= |
|00004b60| 31 33 3b 0a 58 20 20 20 | 58 28 32 2c 31 29 20 2d |13;.X   |X(2,1) -|
|00004b70| 3d 20 32 31 3b 20 58 28 | 32 2c 32 29 20 2d 3d 20 |= 21; X(|2,2) -= |
|00004b80| 32 32 3b 20 58 28 32 2c | 33 29 20 2d 3d 20 32 33 |22; X(2,|3) -= 23|
|00004b90| 3b 0a 58 20 20 20 50 72 | 69 6e 74 28 58 29 3b 0a |;.X   Pr|int(X);.|
|00004ba0| 58 2f 2f 20 20 20 42 61 | 6e 64 4d 61 74 72 69 78 |X//   Ba|ndMatrix|
|00004bb0| 20 42 4d 28 37 2c 32 2c | 32 29 3b 20 42 4d 20 3c | BM(7,2,|2); BM <|
|00004bc0| 3c 20 32 33 3b 0a 58 0a | 58 2f 2f 20 20 20 63 6f |< 23;.X.|X//   co|
|00004bd0| 75 74 20 3c 3c 20 22 5c | 6e 45 6e 64 20 6f 66 20 |ut << "\|nEnd of |
|00004be0| 66 69 66 74 68 20 74 65 | 73 74 5c 6e 22 3b 0a 58 |fifth te|st\n";.X|
|00004bf0| 7d 0a 45 4e 44 5f 4f 46 | 5f 46 49 4c 45 0a 69 66 |}.END_OF|_FILE.if|
|00004c00| 20 74 65 73 74 20 31 37 | 32 35 20 2d 6e 65 20 60 | test 17|25 -ne `|
|00004c10| 77 63 20 2d 63 20 3c 27 | 74 6d 74 35 2e 63 78 78 |wc -c <'|tmt5.cxx|
|00004c20| 27 60 3b 20 74 68 65 6e | 0a 20 20 20 20 65 63 68 |'`; then|.    ech|
|00004c30| 6f 20 73 68 61 72 3a 20 | 5c 22 27 74 6d 74 35 2e |o shar: |\"'tmt5.|
|00004c40| 63 78 78 27 5c 22 20 75 | 6e 70 61 63 6b 65 64 20 |cxx'\" u|npacked |
|00004c50| 77 69 74 68 20 77 72 6f | 6e 67 20 73 69 7a 65 21 |with wro|ng size!|
|00004c60| 0a 66 69 0a 23 20 65 6e | 64 20 6f 66 20 27 74 6d |.fi.# en|d of 'tm|
|00004c70| 74 35 2e 63 78 78 27 0a | 66 69 0a 69 66 20 74 65 |t5.cxx'.|fi.if te|
|00004c80| 73 74 20 2d 66 20 27 74 | 6d 74 36 2e 63 78 78 27 |st -f 't|mt6.cxx'|
|00004c90| 20 2d 61 20 22 24 7b 31 | 7d 22 20 21 3d 20 22 2d | -a "${1|}" != "-|
|00004ca0| 63 22 20 3b 20 74 68 65 | 6e 20 0a 20 20 65 63 68 |c" ; the|n .  ech|
|00004cb0| 6f 20 73 68 61 72 3a 20 | 57 69 6c 6c 20 6e 6f 74 |o shar: |Will not|
|00004cc0| 20 63 6c 6f 62 62 65 72 | 20 65 78 69 73 74 69 6e | clobber| existin|
|00004cd0| 67 20 66 69 6c 65 20 5c | 22 27 74 6d 74 36 2e 63 |g file \|"'tmt6.c|
|00004ce0| 78 78 27 5c 22 0a 65 6c | 73 65 0a 65 63 68 6f 20 |xx'\".el|se.echo |
|00004cf0| 73 68 61 72 3a 20 45 78 | 74 72 61 63 74 69 6e 67 |shar: Ex|tracting|
|00004d00| 20 5c 22 27 74 6d 74 36 | 2e 63 78 78 27 5c 22 20 | \"'tmt6|.cxx'\" |
|00004d10| 5c 28 31 33 32 36 20 63 | 68 61 72 61 63 74 65 72 |\(1326 c|haracter|
|00004d20| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|00004d30| 20 3e 27 74 6d 74 36 2e | 63 78 78 27 20 3c 3c 27 | >'tmt6.|cxx' <<'|
|00004d40| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 27 0a 58 0a 58 |END_OF_F|ILE'.X.X|
|00004d50| 2f 2f 23 64 65 66 69 6e | 65 20 57 41 4e 54 5f 53 |//#defin|e WANT_S|
|00004d60| 54 52 45 41 4d 0a 58 0a | 58 23 69 6e 63 6c 75 64 |TREAM.X.|X#includ|
|00004d70| 65 20 22 69 6e 63 6c 75 | 64 65 2e 68 22 0a 58 0a |e "inclu|de.h".X.|
|00004d80| 58 23 69 6e 63 6c 75 64 | 65 20 22 6e 65 77 6d 61 |X#includ|e "newma|
|00004d90| 74 2e 68 22 0a 58 0a 58 | 0a 58 2f 2a 2a 2a 2a 2a |t.h".X.X|.X/*****|
|00004da0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004db0| 2a 2a 2a 2a 2a 2a 2a 20 | 74 65 73 74 20 70 72 6f |******* |test pro|
|00004dc0| 67 72 61 6d 20 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |gram ***|********|
|00004dd0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004de0| 2a 2a 2a 2f 0a 58 0a 58 | 76 6f 69 64 20 50 72 69 |***/.X.X|void Pri|
|00004df0| 6e 74 28 63 6f 6e 73 74 | 20 4d 61 74 72 69 78 26 |nt(const| Matrix&|
|00004e00| 20 58 29 3b 0a 58 76 6f | 69 64 20 50 72 69 6e 74 | X);.Xvo|id Print|
|00004e10| 28 63 6f 6e 73 74 20 55 | 70 70 65 72 54 72 69 61 |(const U|pperTria|
|00004e20| 6e 67 75 6c 61 72 4d 61 | 74 72 69 78 26 20 58 29 |ngularMa|trix& X)|
|00004e30| 3b 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |;.Xvoid |Print(co|
|00004e40| 6e 73 74 20 44 69 61 67 | 6f 6e 61 6c 4d 61 74 72 |nst Diag|onalMatr|
|00004e50| 69 78 26 20 58 29 3b 0a | 58 76 6f 69 64 20 50 72 |ix& X);.|Xvoid Pr|
|00004e60| 69 6e 74 28 63 6f 6e 73 | 74 20 53 79 6d 6d 65 74 |int(cons|t Symmet|
|00004e70| 72 69 63 4d 61 74 72 69 | 78 26 20 58 29 3b 0a 58 |ricMatri|x& X);.X|
|00004e80| 76 6f 69 64 20 50 72 69 | 6e 74 28 63 6f 6e 73 74 |void Pri|nt(const|
|00004e90| 20 4c 6f 77 65 72 54 72 | 69 61 6e 67 75 6c 61 72 | LowerTr|iangular|
|00004ea0| 4d 61 74 72 69 78 26 20 | 58 29 3b 0a 58 0a 58 76 |Matrix& |X);.X.Xv|
|00004eb0| 6f 69 64 20 74 72 79 6d | 61 74 36 28 29 0a 58 7b |oid trym|at6().X{|
|00004ec0| 0a 58 2f 2f 20 20 20 63 | 6f 75 74 20 3c 3c 20 22 |.X//   c|out << "|
|00004ed0| 5c 6e 53 69 78 74 68 20 | 74 65 73 74 20 6f 66 20 |\nSixth |test of |
|00004ee0| 4d 61 74 72 69 78 20 70 | 61 63 6b 61 67 65 5c 6e |Matrix p|ackage\n|
|00004ef0| 22 3b 0a 58 20 20 20 54 | 72 61 63 65 72 20 65 74 |";.X   T|racer et|
|00004f00| 28 22 53 69 78 74 68 20 | 74 65 73 74 20 6f 66 20 |("Sixth |test of |
|00004f10| 4d 61 74 72 69 78 20 70 | 61 63 6b 61 67 65 22 29 |Matrix p|ackage")|
|00004f20| 3b 0a 58 20 20 20 45 78 | 63 65 70 74 69 6f 6e 3a |;.X   Ex|ception:|
|00004f30| 3a 50 72 69 6e 74 54 72 | 61 63 65 28 54 52 55 45 |:PrintTr|ace(TRUE|
|00004f40| 29 3b 0a 58 0a 58 20 20 | 20 69 6e 74 20 69 2c 6a |);.X.X  | int i,j|
|00004f50| 3b 0a 58 0a 58 0a 58 20 | 20 20 44 69 61 67 6f 6e |;.X.X.X |  Diagon|
|00004f60| 61 6c 4d 61 74 72 69 78 | 20 44 28 36 29 3b 0a 58 |alMatrix| D(6);.X|
|00004f70| 20 20 20 55 70 70 65 72 | 54 72 69 61 6e 67 75 6c |   Upper|Triangul|
|00004f80| 61 72 4d 61 74 72 69 78 | 20 55 28 36 29 3b 0a 58 |arMatrix| U(6);.X|
|00004f90| 20 20 20 66 6f 72 20 28 | 69 3d 31 3b 69 3c 3d 36 |   for (|i=1;i<=6|
|00004fa0| 3b 69 2b 2b 29 20 7b 20 | 66 6f 72 20 28 6a 3d 69 |;i++) { |for (j=i|
|00004fb0| 3b 6a 3c 3d 36 3b 6a 2b | 2b 29 20 55 28 69 2c 6a |;j<=6;j+|+) U(i,j|
|00004fc0| 29 3d 69 2a 69 2a 69 2d | 35 30 3b 20 44 28 69 2c |)=i*i*i-|50; D(i,|
|00004fd0| 69 29 3d 69 2a 69 2b 69 | 2d 31 30 3b 20 7d 0a 58 |i)=i*i+i|-10; }.X|
|00004fe0| 20 20 20 4c 6f 77 65 72 | 54 72 69 61 6e 67 75 6c |   Lower|Triangul|
|00004ff0| 61 72 4d 61 74 72 69 78 | 20 4c 3d 28 55 2a 33 2e |arMatrix| L=(U*3.|
|00005000| 30 29 2e 74 28 29 3b 0a | 58 20 20 20 53 79 6d 6d |0).t();.|X   Symm|
|00005010| 65 74 72 69 63 4d 61 74 | 72 69 78 20 53 28 36 29 |etricMat|rix S(6)|
|00005020| 3b 0a 58 20 20 20 66 6f | 72 20 28 69 3d 31 3b 69 |;.X   fo|r (i=1;i|
|00005030| 3c 3d 36 3b 69 2b 2b 29 | 20 66 6f 72 20 28 6a 3d |<=6;i++)| for (j=|
|00005040| 69 3b 6a 3c 3d 36 3b 6a | 2b 2b 29 20 53 28 69 2c |i;j<=6;j|++) S(i,|
|00005050| 6a 29 3d 69 2a 69 2b 32 | 2e 30 2b 6a 3b 0a 58 20 |j)=i*i+2|.0+j;.X |
|00005060| 20 20 4d 61 74 72 69 78 | 20 4d 44 3d 44 3b 20 4d |  Matrix| MD=D; M|
|00005070| 61 74 72 69 78 20 4d 4c | 3d 4c 3b 20 4d 61 74 72 |atrix ML|=L; Matr|
|00005080| 69 78 20 4d 55 3d 55 3b | 20 4d 61 74 72 69 78 20 |ix MU=U;| Matrix |
|00005090| 4d 53 3d 53 3b 0a 58 20 | 20 20 4d 61 74 72 69 78 |MS=S;.X |  Matrix|
|000050a0| 20 4d 28 36 2c 36 29 3b | 0a 58 20 20 20 66 6f 72 | M(6,6);|.X   for|
|000050b0| 20 28 69 3d 31 3b 69 3c | 3d 36 3b 69 2b 2b 29 20 | (i=1;i<|=6;i++) |
|000050c0| 66 6f 72 20 28 6a 3d 31 | 3b 6a 3c 3d 36 3b 6a 2b |for (j=1|;j<=6;j+|
|000050d0| 2b 29 20 4d 28 69 2c 6a | 29 3d 69 2a 6a 2b 69 2a |+) M(i,j|)=i*j+i*|
|000050e0| 69 2d 31 30 2e 30 3b 20 | 20 0a 58 20 20 20 7b 0a |i-10.0; | .X   {.|
|000050f0| 58 20 20 20 20 20 20 54 | 72 61 63 65 72 20 65 74 |X      T|racer et|
|00005100| 31 28 22 53 74 61 67 65 | 20 31 22 29 3b 0a 58 20 |1("Stage| 1");.X |
|00005110| 20 20 20 20 20 50 72 69 | 6e 74 28 4d 61 74 72 69 |     Pri|nt(Matri|
|00005120| 78 28 4d 53 2b 28 2d 4d | 53 29 29 29 3b 0a 58 20 |x(MS+(-M|S)));.X |
|00005130| 20 20 20 20 20 50 72 69 | 6e 74 28 4d 61 74 72 69 |     Pri|nt(Matri|
|00005140| 78 28 28 53 2b 4d 29 2d | 28 4d 53 2b 4d 29 29 29 |x((S+M)-|(MS+M)))|
|00005150| 3b 0a 58 20 20 20 20 20 | 20 50 72 69 6e 74 28 4d |;.X     | Print(M|
|00005160| 61 74 72 69 78 28 28 4d | 2b 55 29 2d 28 4d 2b 4d |atrix((M|+U)-(M+M|
|00005170| 55 29 29 29 3b 0a 58 20 | 20 20 20 20 20 50 72 69 |U)));.X |     Pri|
|00005180| 6e 74 28 4d 61 74 72 69 | 78 28 28 4d 2b 4c 29 2d |nt(Matri|x((M+L)-|
|00005190| 28 4d 2b 4d 4c 29 29 29 | 3b 0a 58 20 20 20 7d 0a |(M+ML)))|;.X   }.|
|000051a0| 58 20 20 20 7b 0a 58 20 | 20 20 20 20 20 54 72 61 |X   {.X |     Tra|
|000051b0| 63 65 72 20 65 74 31 28 | 22 53 74 61 67 65 20 32 |cer et1(|"Stage 2|
|000051c0| 22 29 3b 0a 58 20 20 20 | 20 20 20 50 72 69 6e 74 |");.X   |   Print|
|000051d0| 28 4d 61 74 72 69 78 28 | 28 4d 2b 44 29 2d 28 4d |(Matrix(|(M+D)-(M|
|000051e0| 2b 4d 44 29 29 29 3b 0a | 58 20 20 20 20 20 20 50 |+MD)));.|X      P|
|000051f0| 72 69 6e 74 28 4d 61 74 | 72 69 78 28 28 55 2b 44 |rint(Mat|rix((U+D|
|00005200| 29 2d 28 4d 55 2b 4d 44 | 29 29 29 3b 0a 58 20 20 |)-(MU+MD|)));.X  |
|00005210| 20 20 20 20 50 72 69 6e | 74 28 4d 61 74 72 69 78 |    Prin|t(Matrix|
|00005220| 28 28 44 2b 4c 29 2d 28 | 4d 4c 2b 4d 44 29 29 29 |((D+L)-(|ML+MD)))|
|00005230| 3b 0a 58 20 20 20 20 20 | 20 50 72 69 6e 74 28 4d |;.X     | Print(M|
|00005240| 61 74 72 69 78 28 28 2d | 55 2b 44 29 2b 4d 55 2d |atrix((-|U+D)+MU-|
|00005250| 4d 44 29 29 3b 0a 58 20 | 20 20 20 20 20 50 72 69 |MD));.X |     Pri|
|00005260| 6e 74 28 4d 61 74 72 69 | 78 28 28 2d 4c 2b 44 29 |nt(Matri|x((-L+D)|
|00005270| 2b 4d 4c 2d 4d 44 29 29 | 3b 0a 58 20 20 20 7d 0a |+ML-MD))|;.X   }.|
|00005280| 58 0a 58 0a 58 0a 58 2f | 2f 20 20 20 63 6f 75 74 |X.X.X.X/|/   cout|
|00005290| 20 3c 3c 20 22 5c 6e 45 | 6e 64 20 6f 66 20 73 69 | << "\nE|nd of si|
|000052a0| 78 74 68 20 74 65 73 74 | 5c 6e 22 3b 0a 58 7d 0a |xth test|\n";.X}.|
|000052b0| 58 0a 45 4e 44 5f 4f 46 | 5f 46 49 4c 45 0a 69 66 |X.END_OF|_FILE.if|
|000052c0| 20 74 65 73 74 20 31 33 | 32 36 20 2d 6e 65 20 60 | test 13|26 -ne `|
|000052d0| 77 63 20 2d 63 20 3c 27 | 74 6d 74 36 2e 63 78 78 |wc -c <'|tmt6.cxx|
|000052e0| 27 60 3b 20 74 68 65 6e | 0a 20 20 20 20 65 63 68 |'`; then|.    ech|
|000052f0| 6f 20 73 68 61 72 3a 20 | 5c 22 27 74 6d 74 36 2e |o shar: |\"'tmt6.|
|00005300| 63 78 78 27 5c 22 20 75 | 6e 70 61 63 6b 65 64 20 |cxx'\" u|npacked |
|00005310| 77 69 74 68 20 77 72 6f | 6e 67 20 73 69 7a 65 21 |with wro|ng size!|
|00005320| 0a 66 69 0a 23 20 65 6e | 64 20 6f 66 20 27 74 6d |.fi.# en|d of 'tm|
|00005330| 74 36 2e 63 78 78 27 0a | 66 69 0a 69 66 20 74 65 |t6.cxx'.|fi.if te|
|00005340| 73 74 20 2d 66 20 27 74 | 6d 74 37 2e 63 78 78 27 |st -f 't|mt7.cxx'|
|00005350| 20 2d 61 20 22 24 7b 31 | 7d 22 20 21 3d 20 22 2d | -a "${1|}" != "-|
|00005360| 63 22 20 3b 20 74 68 65 | 6e 20 0a 20 20 65 63 68 |c" ; the|n .  ech|
|00005370| 6f 20 73 68 61 72 3a 20 | 57 69 6c 6c 20 6e 6f 74 |o shar: |Will not|
|00005380| 20 63 6c 6f 62 62 65 72 | 20 65 78 69 73 74 69 6e | clobber| existin|
|00005390| 67 20 66 69 6c 65 20 5c | 22 27 74 6d 74 37 2e 63 |g file \|"'tmt7.c|
|000053a0| 78 78 27 5c 22 0a 65 6c | 73 65 0a 65 63 68 6f 20 |xx'\".el|se.echo |
|000053b0| 73 68 61 72 3a 20 45 78 | 74 72 61 63 74 69 6e 67 |shar: Ex|tracting|
|000053c0| 20 5c 22 27 74 6d 74 37 | 2e 63 78 78 27 5c 22 20 | \"'tmt7|.cxx'\" |
|000053d0| 5c 28 31 35 32 30 20 63 | 68 61 72 61 63 74 65 72 |\(1520 c|haracter|
|000053e0| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|000053f0| 20 3e 27 74 6d 74 37 2e | 63 78 78 27 20 3c 3c 27 | >'tmt7.|cxx' <<'|
|00005400| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 27 0a 58 0a 58 |END_OF_F|ILE'.X.X|
|00005410| 2f 2f 23 64 65 66 69 6e | 65 20 57 41 4e 54 5f 53 |//#defin|e WANT_S|
|00005420| 54 52 45 41 4d 0a 58 0a | 58 23 69 6e 63 6c 75 64 |TREAM.X.|X#includ|
|00005430| 65 20 22 69 6e 63 6c 75 | 64 65 2e 68 22 0a 58 0a |e "inclu|de.h".X.|
|00005440| 58 23 69 6e 63 6c 75 64 | 65 20 22 6e 65 77 6d 61 |X#includ|e "newma|
|00005450| 74 2e 68 22 0a 58 0a 58 | 0a 58 2f 2a 2a 2a 2a 2a |t.h".X.X|.X/*****|
|00005460| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005470| 2a 2a 2a 2a 2a 2a 2a 20 | 74 65 73 74 20 70 72 6f |******* |test pro|
|00005480| 67 72 61 6d 20 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |gram ***|********|
|00005490| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000054a0| 2a 2a 2a 2f 0a 58 0a 58 | 76 6f 69 64 20 50 72 69 |***/.X.X|void Pri|
|000054b0| 6e 74 28 63 6f 6e 73 74 | 20 4d 61 74 72 69 78 26 |nt(const| Matrix&|
|000054c0| 20 58 29 3b 0a 58 76 6f | 69 64 20 50 72 69 6e 74 | X);.Xvo|id Print|
|000054d0| 28 63 6f 6e 73 74 20 55 | 70 70 65 72 54 72 69 61 |(const U|pperTria|
|000054e0| 6e 67 75 6c 61 72 4d 61 | 74 72 69 78 26 20 58 29 |ngularMa|trix& X)|
|000054f0| 3b 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |;.Xvoid |Print(co|
|00005500| 6e 73 74 20 44 69 61 67 | 6f 6e 61 6c 4d 61 74 72 |nst Diag|onalMatr|
|00005510| 69 78 26 20 58 29 3b 0a | 58 76 6f 69 64 20 50 72 |ix& X);.|Xvoid Pr|
|00005520| 69 6e 74 28 63 6f 6e 73 | 74 20 53 79 6d 6d 65 74 |int(cons|t Symmet|
|00005530| 72 69 63 4d 61 74 72 69 | 78 26 20 58 29 3b 0a 58 |ricMatri|x& X);.X|
|00005540| 76 6f 69 64 20 50 72 69 | 6e 74 28 63 6f 6e 73 74 |void Pri|nt(const|
|00005550| 20 4c 6f 77 65 72 54 72 | 69 61 6e 67 75 6c 61 72 | LowerTr|iangular|
|00005560| 4d 61 74 72 69 78 26 20 | 58 29 3b 0a 58 0a 58 76 |Matrix& |X);.X.Xv|
|00005570| 6f 69 64 20 43 6c 65 61 | 6e 28 4d 61 74 72 69 78 |oid Clea|n(Matrix|
|00005580| 26 2c 20 52 65 61 6c 29 | 3b 0a 58 0a 58 76 6f 69 |&, Real)|;.X.Xvoi|
|00005590| 64 20 74 72 79 6d 61 74 | 37 28 29 0a 58 7b 0a 58 |d trymat|7().X{.X|
|000055a0| 2f 2f 20 20 20 63 6f 75 | 74 20 3c 3c 20 22 5c 6e |//   cou|t << "\n|
|000055b0| 53 65 76 65 6e 74 68 20 | 74 65 73 74 20 6f 66 20 |Seventh |test of |
|000055c0| 4d 61 74 72 69 78 20 70 | 61 63 6b 61 67 65 5c 6e |Matrix p|ackage\n|
|000055d0| 22 3b 0a 58 20 20 20 54 | 72 61 63 65 72 20 65 74 |";.X   T|racer et|
|000055e0| 28 22 53 65 76 65 6e 74 | 68 20 74 65 73 74 20 6f |("Sevent|h test o|
|000055f0| 66 20 4d 61 74 72 69 78 | 20 70 61 63 6b 61 67 65 |f Matrix| package|
|00005600| 22 29 3b 0a 58 20 20 20 | 45 78 63 65 70 74 69 6f |");.X   |Exceptio|
|00005610| 6e 3a 3a 50 72 69 6e 74 | 54 72 61 63 65 28 54 52 |n::Print|Trace(TR|
|00005620| 55 45 29 3b 0a 58 0a 58 | 20 20 20 69 6e 74 20 69 |UE);.X.X|   int i|
|00005630| 2c 6a 3b 0a 58 0a 58 0a | 58 20 20 20 44 69 61 67 |,j;.X.X.|X   Diag|
|00005640| 6f 6e 61 6c 4d 61 74 72 | 69 78 20 44 28 36 29 3b |onalMatr|ix D(6);|
|00005650| 0a 58 20 20 20 55 70 70 | 65 72 54 72 69 61 6e 67 |.X   Upp|erTriang|
|00005660| 75 6c 61 72 4d 61 74 72 | 69 78 20 55 28 36 29 3b |ularMatr|ix U(6);|
|00005670| 0a 58 20 20 20 66 6f 72 | 20 28 69 3d 31 3b 69 3c |.X   for| (i=1;i<|
|00005680| 3d 36 3b 69 2b 2b 29 20 | 7b 20 66 6f 72 20 28 6a |=6;i++) |{ for (j|
|00005690| 3d 69 3b 6a 3c 3d 36 3b | 6a 2b 2b 29 20 55 28 69 |=i;j<=6;|j++) U(i|
|000056a0| 2c 6a 29 3d 69 2a 69 2a | 6a 2d 35 30 3b 20 44 28 |,j)=i*i*|j-50; D(|
|000056b0| 69 2c 69 29 3d 69 2a 69 | 2b 69 2d 31 30 3b 20 7d |i,i)=i*i|+i-10; }|
|000056c0| 0a 58 20 20 20 4c 6f 77 | 65 72 54 72 69 61 6e 67 |.X   Low|erTriang|
|000056d0| 75 6c 61 72 4d 61 74 72 | 69 78 20 4c 3d 28 55 2a |ularMatr|ix L=(U*|
|000056e0| 33 2e 30 29 2e 74 28 29 | 3b 0a 58 20 20 20 53 79 |3.0).t()|;.X   Sy|
|000056f0| 6d 6d 65 74 72 69 63 4d | 61 74 72 69 78 20 53 28 |mmetricM|atrix S(|
|00005700| 36 29 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 3d 31 |6);.X   |for (i=1|
|00005710| 3b 69 3c 3d 36 3b 69 2b | 2b 29 20 66 6f 72 20 28 |;i<=6;i+|+) for (|
|00005720| 6a 3d 69 3b 6a 3c 3d 36 | 3b 6a 2b 2b 29 20 53 28 |j=i;j<=6|;j++) S(|
|00005730| 69 2c 6a 29 3d 69 2a 69 | 2b 32 2e 30 2b 6a 3b 0a |i,j)=i*i|+2.0+j;.|
|00005740| 58 20 20 20 4d 61 74 72 | 69 78 20 4d 44 3d 44 3b |X   Matr|ix MD=D;|
|00005750| 20 4d 61 74 72 69 78 20 | 4d 4c 3d 4c 3b 20 4d 61 | Matrix |ML=L; Ma|
|00005760| 74 72 69 78 20 4d 55 3d | 55 3b 0a 58 20 20 20 4d |trix MU=|U;.X   M|
|00005770| 61 74 72 69 78 20 4d 53 | 3d 53 3b 0a 58 20 20 20 |atrix MS|=S;.X   |
|00005780| 4d 61 74 72 69 78 20 4d | 28 36 2c 36 29 3b 0a 58 |Matrix M|(6,6);.X|
|00005790| 20 20 20 66 6f 72 20 28 | 69 3d 31 3b 69 3c 3d 36 |   for (|i=1;i<=6|
|000057a0| 3b 69 2b 2b 29 20 66 6f | 72 20 28 6a 3d 31 3b 6a |;i++) fo|r (j=1;j|
|000057b0| 3c 3d 36 3b 6a 2b 2b 29 | 20 4d 28 69 2c 6a 29 3d |<=6;j++)| M(i,j)=|
|000057c0| 69 2a 6a 2b 69 2a 69 2d | 31 30 2e 30 3b 20 20 0a |i*j+i*i-|10.0;  .|
|000057d0| 58 20 20 20 7b 0a 58 20 | 20 20 20 20 20 54 72 61 |X   {.X |     Tra|
|000057e0| 63 65 72 20 65 74 31 28 | 22 53 74 61 67 65 20 31 |cer et1(|"Stage 1|
|000057f0| 22 29 3b 0a 58 20 20 20 | 20 20 20 50 72 69 6e 74 |");.X   |   Print|
|00005800| 28 4d 61 74 72 69 78 28 | 28 53 2d 4d 29 2d 28 4d |(Matrix(|(S-M)-(M|
|00005810| 53 2d 4d 29 29 29 3b 0a | 58 20 20 20 20 20 20 50 |S-M)));.|X      P|
|00005820| 72 69 6e 74 28 4d 61 74 | 72 69 78 28 28 2d 4d 2d |rint(Mat|rix((-M-|
|00005830| 53 29 2b 28 4d 53 2b 4d | 29 29 29 3b 0a 58 20 20 |S)+(MS+M|)));.X  |
|00005840| 20 20 20 20 50 72 69 6e | 74 28 4d 61 74 72 69 78 |    Prin|t(Matrix|
|00005850| 28 28 55 2d 4d 29 2d 28 | 4d 55 2d 4d 29 29 29 3b |((U-M)-(|MU-M)));|
|00005860| 0a 58 20 20 20 7d 0a 58 | 20 20 20 7b 0a 58 20 20 |.X   }.X|   {.X  |
|00005870| 20 20 20 20 54 72 61 63 | 65 72 20 65 74 31 28 22 |    Trac|er et1("|
|00005880| 53 74 61 67 65 20 32 22 | 29 3b 0a 58 20 20 20 20 |Stage 2"|);.X    |
|00005890| 20 20 50 72 69 6e 74 28 | 4d 61 74 72 69 78 28 28 |  Print(|Matrix((|
|000058a0| 4c 2d 4d 29 2b 28 4d 2d | 4d 4c 29 29 29 3b 0a 58 |L-M)+(M-|ML)));.X|
|000058b0| 20 20 20 20 20 20 50 72 | 69 6e 74 28 4d 61 74 72 |      Pr|int(Matr|
|000058c0| 69 78 28 28 44 2d 4d 29 | 2b 28 4d 2d 4d 44 29 29 |ix((D-M)|+(M-MD))|
|000058d0| 29 3b 0a 58 20 20 20 20 | 20 20 50 72 69 6e 74 28 |);.X    |  Print(|
|000058e0| 4d 61 74 72 69 78 28 28 | 44 2d 53 29 2b 28 4d 53 |Matrix((|D-S)+(MS|
|000058f0| 2d 4d 44 29 29 29 3b 0a | 58 20 20 20 20 20 20 50 |-MD)));.|X      P|
|00005900| 72 69 6e 74 28 4d 61 74 | 72 69 78 28 28 44 2d 4c |rint(Mat|rix((D-L|
|00005910| 29 2b 28 4d 4c 2d 4d 44 | 29 29 29 3b 0a 58 20 20 |)+(ML-MD|)));.X  |
|00005920| 20 7d 0a 58 0a 58 20 20 | 20 7b 20 4d 3d 4d 55 2e | }.X.X  | { M=MU.|
|00005930| 74 28 29 3b 20 7d 0a 58 | 20 20 20 4c 6f 77 65 72 |t(); }.X|   Lower|
|00005940| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|00005950| 20 4c 59 3d 44 2e 69 28 | 29 2a 55 2e 74 28 29 3b | LY=D.i(|)*U.t();|
|00005960| 0a 58 20 20 20 7b 0a 58 | 20 20 20 20 20 20 54 72 |.X   {.X|      Tr|
|00005970| 61 63 65 72 20 65 74 31 | 28 22 53 74 61 67 65 20 |acer et1|("Stage |
|00005980| 33 22 29 3b 0a 58 20 20 | 20 20 20 20 4d 53 3d 44 |3");.X  |    MS=D|
|00005990| 2a 4c 59 2d 4d 3b 20 43 | 6c 65 61 6e 28 4d 53 2c |*LY-M; C|lean(MS,|
|000059a0| 30 2e 30 30 30 30 30 30 | 30 31 29 3b 20 50 72 69 |0.000000|01); Pri|
|000059b0| 6e 74 28 4d 53 29 3b 0a | 58 20 20 20 20 20 20 4c |nt(MS);.|X      L|
|000059c0| 3d 55 2e 74 28 29 3b 0a | 58 20 20 20 20 20 20 4c |=U.t();.|X      L|
|000059d0| 59 3d 44 2e 69 28 29 2a | 4c 3b 20 4d 53 3d 44 2a |Y=D.i()*|L; MS=D*|
|000059e0| 4c 59 2d 4d 3b 20 43 6c | 65 61 6e 28 4d 53 2c 30 |LY-M; Cl|ean(MS,0|
|000059f0| 2e 30 30 30 30 30 30 30 | 31 29 3b 20 50 72 69 6e |.0000000|1); Prin|
|00005a00| 74 28 4d 53 29 3b 0a 58 | 20 20 20 7d 0a 58 0a 58 |t(MS);.X|   }.X.X|
|00005a10| 2f 2f 20 20 20 63 6f 75 | 74 20 3c 3c 20 22 5c 6e |//   cou|t << "\n|
|00005a20| 45 6e 64 20 6f 66 20 73 | 65 76 65 6e 74 68 20 74 |End of s|eventh t|
|00005a30| 65 73 74 5c 6e 22 3b 0a | 58 7d 0a 45 4e 44 5f 4f |est\n";.|X}.END_O|
|00005a40| 46 5f 46 49 4c 45 0a 69 | 66 20 74 65 73 74 20 31 |F_FILE.i|f test 1|
|00005a50| 35 32 30 20 2d 6e 65 20 | 60 77 63 20 2d 63 20 3c |520 -ne |`wc -c <|
|00005a60| 27 74 6d 74 37 2e 63 78 | 78 27 60 3b 20 74 68 65 |'tmt7.cx|x'`; the|
|00005a70| 6e 0a 20 20 20 20 65 63 | 68 6f 20 73 68 61 72 3a |n.    ec|ho shar:|
|00005a80| 20 5c 22 27 74 6d 74 37 | 2e 63 78 78 27 5c 22 20 | \"'tmt7|.cxx'\" |
|00005a90| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|00005aa0| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|00005ab0| 6e 64 20 6f 66 20 27 74 | 6d 74 37 2e 63 78 78 27 |nd of 't|mt7.cxx'|
|00005ac0| 0a 66 69 0a 69 66 20 74 | 65 73 74 20 2d 66 20 27 |.fi.if t|est -f '|
|00005ad0| 74 6d 74 38 2e 63 78 78 | 27 20 2d 61 20 22 24 7b |tmt8.cxx|' -a "${|
|00005ae0| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|00005af0| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en .  ec|ho shar:|
|00005b00| 20 57 69 6c 6c 20 6e 6f | 74 20 63 6c 6f 62 62 65 | Will no|t clobbe|
|00005b10| 72 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |r existi|ng file |
|00005b20| 5c 22 27 74 6d 74 38 2e | 63 78 78 27 5c 22 0a 65 |\"'tmt8.|cxx'\".e|
|00005b30| 6c 73 65 0a 65 63 68 6f | 20 73 68 61 72 3a 20 45 |lse.echo| shar: E|
|00005b40| 78 74 72 61 63 74 69 6e | 67 20 5c 22 27 74 6d 74 |xtractin|g \"'tmt|
|00005b50| 38 2e 63 78 78 27 5c 22 | 20 5c 28 31 34 31 39 20 |8.cxx'\"| \(1419 |
|00005b60| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00005b70| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 74 6d 74 38 | "s/^X//|" >'tmt8|
|00005b80| 2e 63 78 78 27 20 3c 3c | 27 45 4e 44 5f 4f 46 5f |.cxx' <<|'END_OF_|
|00005b90| 46 49 4c 45 27 0a 58 0a | 58 2f 2f 23 64 65 66 69 |FILE'.X.|X//#defi|
|00005ba0| 6e 65 20 57 41 4e 54 5f | 53 54 52 45 41 4d 0a 58 |ne WANT_|STREAM.X|
|00005bb0| 0a 58 0a 58 23 69 6e 63 | 6c 75 64 65 20 22 69 6e |.X.X#inc|lude "in|
|00005bc0| 63 6c 75 64 65 2e 68 22 | 0a 58 0a 58 23 69 6e 63 |clude.h"|.X.X#inc|
|00005bd0| 6c 75 64 65 20 22 6e 65 | 77 6d 61 74 2e 68 22 0a |lude "ne|wmat.h".|
|00005be0| 58 0a 58 0a 58 2f 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |X.X.X/**|********|
|00005bf0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005c00| 2a 2a 20 74 65 73 74 20 | 70 72 6f 67 72 61 6d 20 |** test |program |
|00005c10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005c20| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|00005c30| 58 0a 58 76 6f 69 64 20 | 50 72 69 6e 74 28 63 6f |X.Xvoid |Print(co|
|00005c40| 6e 73 74 20 4d 61 74 72 | 69 78 26 20 58 29 3b 0a |nst Matr|ix& X);.|
|00005c50| 58 76 6f 69 64 20 50 72 | 69 6e 74 28 63 6f 6e 73 |Xvoid Pr|int(cons|
|00005c60| 74 20 55 70 70 65 72 54 | 72 69 61 6e 67 75 6c 61 |t UpperT|riangula|
|00005c70| 72 4d 61 74 72 69 78 26 | 20 58 29 3b 0a 58 76 6f |rMatrix&| X);.Xvo|
|00005c80| 69 64 20 50 72 69 6e 74 | 28 63 6f 6e 73 74 20 44 |id Print|(const D|
|00005c90| 69 61 67 6f 6e 61 6c 4d | 61 74 72 69 78 26 20 58 |iagonalM|atrix& X|
|00005ca0| 29 3b 0a 58 76 6f 69 64 | 20 50 72 69 6e 74 28 63 |);.Xvoid| Print(c|
|00005cb0| 6f 6e 73 74 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |onst Sym|metricMa|
|00005cc0| 74 72 69 78 26 20 58 29 | 3b 0a 58 76 6f 69 64 20 |trix& X)|;.Xvoid |
|00005cd0| 50 72 69 6e 74 28 63 6f | 6e 73 74 20 4c 6f 77 65 |Print(co|nst Lowe|
|00005ce0| 72 54 72 69 61 6e 67 75 | 6c 61 72 4d 61 74 72 69 |rTriangu|larMatri|
|00005cf0| 78 26 20 58 29 3b 0a 58 | 0a 58 76 6f 69 64 20 74 |x& X);.X|.Xvoid t|
|00005d00| 72 79 6d 61 74 38 28 29 | 0a 58 7b 0a 58 2f 2f 20 |rymat8()|.X{.X// |
|00005d10| 20 20 63 6f 75 74 20 3c | 3c 20 22 5c 6e 45 69 67 |  cout <|< "\nEig|
|00005d20| 68 74 68 20 74 65 73 74 | 20 6f 66 20 4d 61 74 72 |hth test| of Matr|
|00005d30| 69 78 20 70 61 63 6b 61 | 67 65 5c 6e 22 3b 0a 58 |ix packa|ge\n";.X|
|00005d40| 20 20 20 54 72 61 63 65 | 72 20 65 74 28 22 45 69 |   Trace|r et("Ei|
|00005d50| 67 68 74 68 20 74 65 73 | 74 20 6f 66 20 4d 61 74 |ghth tes|t of Mat|
|00005d60| 72 69 78 20 70 61 63 6b | 61 67 65 22 29 3b 0a 58 |rix pack|age");.X|
|00005d70| 20 20 20 45 78 63 65 70 | 74 69 6f 6e 3a 3a 50 72 |   Excep|tion::Pr|
|00005d80| 69 6e 74 54 72 61 63 65 | 28 54 52 55 45 29 3b 0a |intTrace|(TRUE);.|
|00005d90| 58 0a 58 20 20 20 69 6e | 74 20 69 3b 0a 58 0a 58 |X.X   in|t i;.X.X|
|00005da0| 0a 58 20 20 20 44 69 61 | 67 6f 6e 61 6c 4d 61 74 |.X   Dia|gonalMat|
|00005db0| 72 69 78 20 44 28 36 29 | 3b 0a 58 20 20 20 66 6f |rix D(6)|;.X   fo|
|00005dc0| 72 20 28 69 3d 31 3b 69 | 3c 3d 36 3b 69 2b 2b 29 |r (i=1;i|<=6;i++)|
|00005dd0| 20 20 44 28 69 2c 69 29 | 3d 69 2a 69 2b 69 2d 31 |  D(i,i)|=i*i+i-1|
|00005de0| 30 3b 0a 58 20 20 20 44 | 69 61 67 6f 6e 61 6c 4d |0;.X   D|iagonalM|
|00005df0| 61 74 72 69 78 20 44 32 | 3d 44 3b 0a 58 20 20 20 |atrix D2|=D;.X   |
|00005e00| 4d 61 74 72 69 78 20 4d | 44 3d 44 3b 0a 58 0a 58 |Matrix M|D=D;.X.X|
|00005e10| 20 20 20 44 69 61 67 6f | 6e 61 6c 4d 61 74 72 69 |   Diago|nalMatri|
|00005e20| 78 20 44 31 28 36 29 3b | 20 66 6f 72 20 28 69 3d |x D1(6);| for (i=|
|00005e30| 31 3b 69 3c 3d 36 3b 69 | 2b 2b 29 20 44 31 28 69 |1;i<=6;i|++) D1(i|
|00005e40| 2c 69 29 3d 2d 31 30 30 | 2b 69 2a 69 2a 69 3b 0a |,i)=-100|+i*i*i;.|
|00005e50| 58 20 20 20 4d 61 74 72 | 69 78 20 4d 44 31 3d 44 |X   Matr|ix MD1=D|
|00005e60| 31 3b 0a 58 20 20 20 50 | 72 69 6e 74 28 4d 61 74 |1;.X   P|rint(Mat|
|00005e70| 72 69 78 28 44 2a 44 31 | 2d 4d 44 2a 4d 44 31 29 |rix(D*D1|-MD*MD1)|
|00005e80| 29 3b 0a 58 20 20 20 50 | 72 69 6e 74 28 4d 61 74 |);.X   P|rint(Mat|
|00005e90| 72 69 78 28 28 2d 44 29 | 2a 44 31 2b 4d 44 2a 4d |rix((-D)|*D1+MD*M|
|00005ea0| 44 31 29 29 3b 0a 58 20 | 20 20 50 72 69 6e 74 28 |D1));.X |  Print(|
|00005eb0| 4d 61 74 72 69 78 28 44 | 2a 28 2d 44 31 29 2b 4d |Matrix(D|*(-D1)+M|
|00005ec0| 44 2a 4d 44 31 29 29 3b | 0a 58 20 20 20 44 69 61 |D*MD1));|.X   Dia|
|00005ed0| 67 6f 6e 61 6c 4d 61 74 | 72 69 78 20 44 58 3d 44 |gonalMat|rix DX=D|
|00005ee0| 3b 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 54 |;.X   {.|X      T|
|00005ef0| 72 61 63 65 72 20 65 74 | 31 28 22 53 74 61 67 65 |racer et|1("Stage|
|00005f00| 20 31 22 29 3b 0a 58 20 | 20 20 20 20 20 44 58 3d | 1");.X |     DX=|
|00005f10| 28 44 58 2b 44 31 29 2a | 44 58 3b 20 50 72 69 6e |(DX+D1)*|DX; Prin|
|00005f20| 74 28 4d 61 74 72 69 78 | 28 44 58 2d 28 4d 44 2b |t(Matrix|(DX-(MD+|
|00005f30| 4d 44 31 29 2a 4d 44 29 | 29 3b 0a 58 20 20 20 20 |MD1)*MD)|);.X    |
|00005f40| 20 20 44 58 3d 44 3b 0a | 58 20 20 20 20 20 20 44 |  DX=D;.|X      D|
|00005f50| 58 3d 2d 44 58 2a 44 58 | 2b 28 44 58 2d 28 2d 44 |X=-DX*DX|+(DX-(-D|
|00005f60| 31 29 29 2a 28 28 2d 44 | 31 29 2b 44 58 29 3b 0a |1))*((-D|1)+DX);.|
|00005f70| 58 23 69 66 64 65 66 20 | 5f 5f 47 4e 55 47 5f 5f |X#ifdef |__GNUG__|
|00005f80| 0a 58 20 20 20 20 20 20 | 4d 61 74 72 69 78 20 4d |.X      |Matrix M|
|00005f90| 58 20 3d 20 4d 61 74 72 | 69 78 28 4d 44 31 29 3b |X = Matr|ix(MD1);|
|00005fa0| 0a 58 20 20 20 20 20 20 | 4d 44 31 3d 44 58 2b 28 |.X      |MD1=DX+(|
|00005fb0| 4d 58 2e 74 28 29 29 2a | 28 4d 58 2e 74 28 29 29 |MX.t())*|(MX.t())|
|00005fc0| 3b 20 50 72 69 6e 74 28 | 4d 44 31 29 3b 0a 58 23 |; Print(|MD1);.X#|
|00005fd0| 65 6c 73 65 0a 58 20 20 | 20 20 20 20 4d 44 31 3d |else.X  |    MD1=|
|00005fe0| 44 58 2b 28 4d 61 74 72 | 69 78 28 4d 44 31 29 2e |DX+(Matr|ix(MD1).|
|00005ff0| 74 28 29 29 2a 28 4d 61 | 74 72 69 78 28 4d 44 31 |t())*(Ma|trix(MD1|
|00006000| 29 2e 74 28 29 29 3b 20 | 50 72 69 6e 74 28 4d 44 |).t()); |Print(MD|
|00006010| 31 29 3b 0a 58 23 65 6e | 64 69 66 0a 58 20 20 20 |1);.X#en|dif.X   |
|00006020| 20 20 20 44 58 3d 44 3b | 20 44 58 3d 44 58 3b 20 |   DX=D;| DX=DX; |
|00006030| 44 58 3d 44 32 2d 44 58 | 3b 20 50 72 69 6e 74 28 |DX=D2-DX|; Print(|
|00006040| 44 69 61 67 6f 6e 61 6c | 4d 61 74 72 69 78 28 44 |Diagonal|Matrix(D|
|00006050| 58 29 29 3b 0a 58 20 20 | 20 20 20 20 44 58 3d 44 |X));.X  |    DX=D|
|00006060| 3b 0a 58 20 20 20 7d 0a | 58 20 20 20 7b 0a 58 20 |;.X   }.|X   {.X |
|00006070| 20 20 20 20 20 54 72 61 | 63 65 72 20 65 74 31 28 |     Tra|cer et1(|
|00006080| 22 53 74 61 67 65 20 32 | 22 29 3b 0a 58 20 20 20 |"Stage 2|");.X   |
|00006090| 20 20 20 44 2e 52 65 6c | 65 61 73 65 28 32 29 3b |   D.Rel|ease(2);|
|000060a0| 0a 58 20 20 20 20 20 20 | 44 31 3d 44 3b 20 44 32 |.X      |D1=D; D2|
|000060b0| 3d 44 3b 0a 58 20 20 20 | 20 20 20 50 72 69 6e 74 |=D;.X   |   Print|
|000060c0| 28 44 69 61 67 6f 6e 61 | 6c 4d 61 74 72 69 78 28 |(Diagona|lMatrix(|
|000060d0| 44 31 2d 44 58 29 29 3b | 0a 58 20 20 20 20 20 20 |D1-DX));|.X      |
|000060e0| 50 72 69 6e 74 28 44 69 | 61 67 6f 6e 61 6c 4d 61 |Print(Di|agonalMa|
|000060f0| 74 72 69 78 28 44 32 2d | 44 58 29 29 3b 0a 58 20 |trix(D2-|DX));.X |
|00006100| 20 20 20 20 20 4d 44 31 | 3d 31 2e 30 3b 0a 58 20 |     MD1|=1.0;.X |
|00006110| 20 20 20 20 20 50 72 69 | 6e 74 28 4d 61 74 72 69 |     Pri|nt(Matri|
|00006120| 78 28 4d 44 31 2d 31 2e | 30 29 29 3b 0a 58 20 20 |x(MD1-1.|0));.X  |
|00006130| 20 7d 0a 58 0a 58 2f 2f | 20 20 20 63 6f 75 74 20 | }.X.X//|   cout |
|00006140| 3c 3c 20 22 5c 6e 45 6e | 64 20 6f 66 20 65 69 67 |<< "\nEn|d of eig|
|00006150| 68 74 68 20 74 65 73 74 | 5c 6e 22 3b 0a 58 7d 0a |hth test|\n";.X}.|
|00006160| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 0a 69 66 20 74 |END_OF_F|ILE.if t|
|00006170| 65 73 74 20 31 34 31 39 | 20 2d 6e 65 20 60 77 63 |est 1419| -ne `wc|
|00006180| 20 2d 63 20 3c 27 74 6d | 74 38 2e 63 78 78 27 60 | -c <'tm|t8.cxx'`|
|00006190| 3b 20 74 68 65 6e 0a 20 | 20 20 20 65 63 68 6f 20 |; then. |   echo |
|000061a0| 73 68 61 72 3a 20 5c 22 | 27 74 6d 74 38 2e 63 78 |shar: \"|'tmt8.cx|
|000061b0| 78 27 5c 22 20 75 6e 70 | 61 63 6b 65 64 20 77 69 |x'\" unp|acked wi|
|000061c0| 74 68 20 77 72 6f 6e 67 | 20 73 69 7a 65 21 0a 66 |th wrong| size!.f|
|000061d0| 69 0a 23 20 65 6e 64 20 | 6f 66 20 27 74 6d 74 38 |i.# end |of 'tmt8|
|000061e0| 2e 63 78 78 27 0a 66 69 | 0a 69 66 20 74 65 73 74 |.cxx'.fi|.if test|
|000061f0| 20 2d 66 20 27 74 6d 74 | 39 2e 63 78 78 27 20 2d | -f 'tmt|9.cxx' -|
|00006200| 61 20 22 24 7b 31 7d 22 | 20 21 3d 20 22 2d 63 22 |a "${1}"| != "-c"|
|00006210| 20 3b 20 74 68 65 6e 20 | 0a 20 20 65 63 68 6f 20 | ; then |.  echo |
|00006220| 73 68 61 72 3a 20 57 69 | 6c 6c 20 6e 6f 74 20 63 |shar: Wi|ll not c|
|00006230| 6c 6f 62 62 65 72 20 65 | 78 69 73 74 69 6e 67 20 |lobber e|xisting |
|00006240| 66 69 6c 65 20 5c 22 27 | 74 6d 74 39 2e 63 78 78 |file \"'|tmt9.cxx|
|00006250| 27 5c 22 0a 65 6c 73 65 | 0a 65 63 68 6f 20 73 68 |'\".else|.echo sh|
|00006260| 61 72 3a 20 45 78 74 72 | 61 63 74 69 6e 67 20 5c |ar: Extr|acting \|
|00006270| 22 27 74 6d 74 39 2e 63 | 78 78 27 5c 22 20 5c 28 |"'tmt9.c|xx'\" \(|
|00006280| 31 35 37 37 20 63 68 61 | 72 61 63 74 65 72 73 5c |1577 cha|racters\|
|00006290| 29 0a 73 65 64 20 22 73 | 2f 5e 58 2f 2f 22 20 3e |).sed "s|/^X//" >|
|000062a0| 27 74 6d 74 39 2e 63 78 | 78 27 20 3c 3c 27 45 4e |'tmt9.cx|x' <<'EN|
|000062b0| 44 5f 4f 46 5f 46 49 4c | 45 27 0a 58 0a 58 2f 2f |D_OF_FIL|E'.X.X//|
|000062c0| 23 64 65 66 69 6e 65 20 | 57 41 4e 54 5f 53 54 52 |#define |WANT_STR|
|000062d0| 45 41 4d 0a 58 0a 58 0a | 58 23 69 6e 63 6c 75 64 |EAM.X.X.|X#includ|
|000062e0| 65 20 22 69 6e 63 6c 75 | 64 65 2e 68 22 0a 58 23 |e "inclu|de.h".X#|
|000062f0| 69 6e 63 6c 75 64 65 20 | 22 6e 65 77 6d 61 74 2e |include |"newmat.|
|00006300| 68 22 0a 58 0a 58 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a |h".X.X.X|/*******|
|00006310| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00006320| 2a 2a 2a 2a 2a 20 74 65 | 73 74 20 70 72 6f 67 72 |***** te|st progr|
|00006330| 61 6d 20 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |am *****|********|
|00006340| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00006350| 2a 2f 0a 58 0a 58 76 6f | 69 64 20 50 72 69 6e 74 |*/.X.Xvo|id Print|
|00006360| 28 63 6f 6e 73 74 20 4d | 61 74 72 69 78 26 20 58 |(const M|atrix& X|
|00006370| 29 3b 0a 58 76 6f 69 64 | 20 50 72 69 6e 74 28 63 |);.Xvoid| Print(c|
|00006380| 6f 6e 73 74 20 55 70 70 | 65 72 54 72 69 61 6e 67 |onst Upp|erTriang|
|00006390| 75 6c 61 72 4d 61 74 72 | 69 78 26 20 58 29 3b 0a |ularMatr|ix& X);.|
|000063a0| 58 76 6f 69 64 20 50 72 | 69 6e 74 28 63 6f 6e 73 |Xvoid Pr|int(cons|
|000063b0| 74 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |t Diagon|alMatrix|
|000063c0| 26 20 58 29 3b 0a 58 76 | 6f 69 64 20 50 72 69 6e |& X);.Xv|oid Prin|
|000063d0| 74 28 63 6f 6e 73 74 20 | 53 79 6d 6d 65 74 72 69 |t(const |Symmetri|
|000063e0| 63 4d 61 74 72 69 78 26 | 20 58 29 3b 0a 58 76 6f |cMatrix&| X);.Xvo|
|000063f0| 69 64 20 50 72 69 6e 74 | 28 63 6f 6e 73 74 20 4c |id Print|(const L|
+--------+-------------------------+-------------------------+--------+--------+