home *** CD-ROM | disk | FTP | other *** search
| Unknown | 1993-08-10 | 46.9 KB |
open in: 
MacOS 8.1
|
Win98
|
DOS
view JSON data
|
view as text
This file was not able to be converted.
This format is not currently supported by dexvert.
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| file
| data
| default
| |
| 100%
| gt2
| Kopftext: 'holesky'
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 68 6f 6c 65 73 6b 79 91 | ff 2a aa 2c 96 03 55 55 |holesky.|.*.,..UU|
|00000010| f3 0a 0b a0 62 3e 00 0a | 00 00 00 0a 00 00 00 06 |....b>..|........|
|00000020| 63 6d 6d 69 31 30 b5 4c | 44 90 47 1d 4c 8d 9f fc |cmmi10.L|D.G.L...|
|00000030| 5e ff f3 09 30 65 97 72 | 00 07 00 00 00 07 00 00 |^...0e.r|........|
|00000040| 00 05 63 6d 6d 69 37 b4 | 54 8e 91 09 a2 0b b2 66 |..cmmi7.|T......f|
|00000050| 61 63 74 6f 72 93 61 6e | 64 93 73 6f 6c 76 90 b8 |actor.an|d.solv..|
|00000060| e3 65 8e 8d 92 01 ac 6f | 56 37 33 8e 8e a1 8d 8d |.e.....o|V73.....|
|00000070| 8d 8d 91 23 ff f2 da 4c | 44 4c 66 61 63 74 6f 72 |...#...L|DLfactor|
|00000080| 2c 91 05 3f fd 4c 44 4c | 73 6f 6c 76 65 8e 8e 8e |,..?.LDL|solve...|
|00000090| 8e a1 8d 8d 8d 8d 91 0f | 00 00 4c 55 66 61 63 74 |........|..LUfact|
|000000a0| 6f 72 2c 96 05 3f fd 4c | 55 73 6f 6c 76 65 2c 93 |or,..?.L|Usolve,.|
|000000b0| 4c 55 54 73 6f 6c 76 65 | 8e 8e 8e 8d 92 00 e2 2b |LUTsolve|.......+|
|000000c0| 56 b5 4c 55 91 04 6c 70 | b2 66 61 63 74 6f 72 96 |V.LU..lp|.factor.|
|000000d0| 03 55 55 28 47 61 75 73 | 73 69 61 6e 93 65 6c 69 |.UU(Gaus|sian.eli|
|000000e0| 6d 69 6e 61 74 69 6f 6e | 29 93 61 6e 64 93 73 6f |mination|).and.so|
|000000f0| 6c 76 90 b8 e3 65 8e 8d | 92 01 ac 6f 56 37 35 8e |lv...e..|...oV75.|
|00000100| 8e a1 8d 8d 8d 8d 91 23 | ff f2 da 4c 55 63 6f 6e |.......#|...LUcon|
|00000110| 64 65 73 74 2c 91 05 3f | fd 6d 8d 91 00 a1 45 89 |dest,..?|.m....E.|
|00000120| 00 00 66 66 00 03 26 66 | 8e 91 03 c7 ab 69 6e 76 |..ff..&f|.....inv|
|00000130| 65 72 73 65 8e 8e 8e 8e | a1 8d 8d 8d 8d 91 0f 00 |erse....|........|
|00000140| 00 51 52 66 61 63 74 6f | 72 2c 96 05 3f fd 51 52 |.QRfacto|r,..?.QR|
|00000150| 73 6f 6c 76 65 2c 93 51 | 52 54 73 6f 6c 76 65 8e |solve,.Q|RTsolve.|
|00000160| 8e 8e 8d 92 00 e2 2b 56 | b5 51 52 91 03 69 1c b2 |......+V|.QR..i..|
|00000170| 66 61 63 74 6f 72 96 03 | 55 55 61 6e 64 93 73 6f |factor..|UUand.so|
|00000180| 6c 76 90 b8 e3 65 8e 8d | 92 01 ac 6f 56 37 37 8e |lv...e..|...oV77.|
|00000190| 8e a1 8d 8d 8d 8d 91 23 | ff f2 da 51 52 43 50 66 |.......#|...QRCPf|
|000001a0| 61 63 74 6f 72 2c 91 05 | 3f fd 51 52 43 50 73 6f |actor,..|?.QRCPso|
|000001b0| 6c 76 65 8e 8e 8e 8d 92 | 00 f6 2b 58 b2 77 69 74 |lve.....|..+X.wit|
|000001c0| 68 2f 6f 75 74 96 03 55 | 55 63 6f 6c 75 6d 6e 93 |h/out..U|Ucolumn.|
|000001d0| 70 69 76 90 b8 e3 6f 74 | 69 6e 67 8e 8e a1 8d 8d |piv...ot|ing.....|
|000001e0| 8d 8d 91 23 ff f2 da 51 | 52 63 6f 6e 64 65 73 74 |...#...Q|Rcondest|
|000001f0| 8e 8e 8e 8e a1 8d 8d 8d | 8d 91 0f 00 00 6d 61 6b |........|.....mak|
|00000200| 65 51 2c 91 05 3f fd 6d | 61 6b 65 52 8e 8e 8e 8d |eQ,..?.m|akeR....|
|00000210| 92 00 e2 2b 56 b2 65 78 | 74 72 61 63 74 96 03 55 |...+V.ex|tract..U|
|00000220| 55 6d 61 74 72 69 63 65 | 73 93 66 72 6f 6d 93 63 |Umatrice|s.from.c|
|00000230| 6f 6d 70 61 63 74 93 66 | 6f 72 6d 8e 8d 92 01 ac |ompact.f|orm.....|
|00000240| 6f 56 37 39 8e 8e a1 8d | 8d 92 00 f6 2b 58 28 b5 |oV79....|....+X(.|
|00000250| 51 52 91 03 69 1c b2 6f | 6e 6c 79 29 8e 8e a1 8d |QR..i..o|nly)....|
|00000260| 8d 8d 8d 91 0f 00 00 da | 67 69 76 65 6e 73 2c 91 |........|givens,.|
|00000270| 05 3f fd 72 6f 74 8d 91 | 00 a1 45 89 00 00 66 66 |.?.rot..|..E...ff|
|00000280| 00 03 26 66 8e 91 03 c7 | ab 63 6f 6c 73 8e 8e 8e |..&f....|.cols...|
|00000290| 8d 92 00 e2 2b 56 b2 63 | 6f 6d 70 75 74 65 96 03 |....+V.c|ompute..|
|000002a0| 55 55 61 6e 64 93 61 70 | 70 6c 79 93 47 69 76 90 |UUand.ap|ply.Giv.|
|000002b0| b8 e3 65 6e 73 27 93 72 | 6f 74 61 74 69 6f 6e 73 |..ens'.r|otations|
|000002c0| 8e 8d 92 01 ac 6f 56 38 | 30 8e 8e a1 8d 8d 8d 8d |.....oV8|0.......|
|000002d0| 91 23 ff f2 da 72 6f 74 | 8d 91 00 a1 45 89 00 00 |.#...rot|....E...|
|000002e0| 66 66 00 03 26 66 8e 96 | 03 c7 ab 72 6f 77 73 2c |ff..&f..|...rows,|
|000002f0| 91 05 3f fd 72 6f 74 8d | 91 00 a1 45 89 00 00 66 |..?.rot.|...E...f|
|00000300| 66 00 03 26 66 8e 93 76 | 65 63 8e 8e 8e 8e a1 8d |f..&f..v|ec......|
|00000310| 8d 8d 8d 91 0f 00 00 68 | 68 76 65 63 2c 91 05 3f |.......h|hvec,..?|
|00000320| fd 68 68 74 72 76 65 63 | 8e 8e 8e 8d 92 00 e2 2b |.hhtrvec|.......+|
|00000330| 56 b2 48 6f 75 73 65 68 | 6f 6c 64 65 72 91 03 55 |V.Househ|older..U|
|00000340| 55 74 72 61 6e 73 66 6f | 72 6d 61 74 69 6f 6e 73 |Utransfo|rmations|
|00000350| 8e 8d 92 01 ac 6f 56 38 | 32 8e 8e a1 8d 8d 8d 8d |.....oV8|2.......|
|00000360| 91 23 ff f2 da 68 68 74 | 72 72 6f 77 73 2c 91 05 |.#...hht|rrows,..|
|00000370| 3f fd 68 68 74 72 63 6f | 6c 73 8e 8e 8e 8e a1 8d |?.hhtrco|ls......|
|00000380| 8d 8d 8d 91 0f 00 00 44 | 73 6f 6c 76 65 2c 96 05 |.......D|solve,..|
|00000390| 3f fd 4c 73 6f 6c 76 65 | 2c 93 4c 54 73 6f 6c 76 |?.Lsolve|,.LTsolv|
|000003a0| 65 8e 8e 8e 8d 92 00 e2 | 2b 56 b2 73 6f 6c 76 90 |e.......|+V.solv.|
|000003b0| b8 e3 65 96 03 55 55 66 | 6f 72 93 64 69 61 67 6f |..e..UUf|or.diago|
|000003c0| 6e 61 6c 93 61 6e 64 93 | 74 72 69 61 6e 67 75 6c |nal.and.|triangul|
|000003d0| 61 72 93 6d 61 74 72 69 | 63 65 73 8e 8d 92 01 ac |ar.matri|ces.....|
|000003e0| 6f 56 38 34 8e 8e a1 8d | 8d 8d 8d 91 23 ff f2 da |oV84....|....#...|
|000003f0| 55 73 6f 6c 76 65 2c 91 | 05 3f fd 55 54 73 6f 6c |Usolve,.|.?.UTsol|
|00000400| 76 65 8e 8e 8e 8e a6 8d | 8d 8d 8d 91 0f 00 00 4c |ve......|.......L|
|00000410| 44 4c 75 70 64 61 74 65 | 2c 91 05 3f fd 51 52 75 |DLupdate|,..?.QRu|
|00000420| 70 64 61 74 65 8e 8e 8e | 8d 92 00 e2 2b 56 b2 75 |pdate...|....+V.u|
|00000430| 70 9a 47 1d 64 61 74 65 | 96 03 55 55 72 6f 75 74 |p.G.date|..UUrout|
|00000440| 69 6e 65 73 93 66 6f 72 | 93 b5 4c 44 98 4c 8d 9f |ines.for|..LD.L..|
|00000450| fc 5e ff b4 54 8e 91 09 | a2 0b b2 61 6e 64 8e 8d |.^..T...|...and..|
|00000460| 92 01 ac 6f 56 38 36 8e | 8e a1 8d 8d 92 00 f6 2b |...oV86.|.......+|
|00000470| 58 b5 51 52 91 03 69 1c | b2 66 61 63 74 6f 72 69 |X.QR..i.|.factori|
|00000480| 73 61 74 69 6f 6e 73 8e | 8e a1 8d 8d 8d 8d 91 0f |sations.|........|
|00000490| 00 00 da 73 63 68 75 72 | 2c 96 05 3f fd 73 79 6d |...schur|,..?.sym|
|000004a0| 6d 65 69 67 2c 93 74 72 | 69 65 69 67 8e 8e 8e 8d |meig,.tr|ieig....|
|000004b0| 92 00 e2 2b 56 b2 65 69 | 67 65 6e 90 b8 e3 76 91 |...+V.ei|gen...v.|
|000004c0| ff 71 c7 61 6c 75 65 91 | 03 55 55 72 6f 75 74 69 |.q.alue.|.UUrouti|
|000004d0| 6e 65 73 8e 8d 92 01 ac | 6f 56 38 38 8e 8e a1 8d |nes.....|oV88....|
|000004e0| 8d 8d 8d 91 0f 00 00 da | 73 63 68 75 72 8d 91 00 |........|schur...|
|000004f0| a1 45 89 00 00 66 66 00 | 03 26 66 8e 96 03 c7 ab |.E...ff.|.&f.....|
|00000500| 65 76 61 6c 73 2c 91 05 | 3f fd 73 63 68 75 72 8d |evals,..|?.schur.|
|00000510| 91 00 a1 45 89 00 00 66 | 66 00 03 26 66 8e 93 76 |...E...f|f..&f..v|
|00000520| 65 63 73 8e 8e 8e 8d 92 | 00 e2 2b 56 b2 65 69 67 |ecs.....|..+V.eig|
|00000530| 65 6e 95 b8 e3 76 91 ff | 71 c7 61 6c 75 65 2f 76 |en...v..|q.alue/v|
|00000540| 93 65 63 74 6f 72 96 03 | 55 55 65 78 74 72 61 63 |.ector..|UUextrac|
|00000550| 74 69 6f 6e 93 72 6f 75 | 74 69 6e 65 73 8e 8d 92 |tion.rou|tines...|
|00000560| 01 ac 6f 56 39 30 8e 8e | a1 8d 8d 8d 8d 91 0f 00 |..oV90..|........|
|00000570| 00 da 73 76 64 2c 91 05 | 3f fd 62 69 73 76 64 8e |..svd,..|?.bisvd.|
|00000580| 8e 8e 8d 92 00 e2 2b 56 | b2 73 69 6e 67 75 6c 61 |......+V|.singula|
|00000590| 72 96 03 55 55 76 91 ff | 71 c7 61 6c 75 65 93 64 |r..UUv..|q.alue.d|
|000005a0| 65 63 6f 6d 70 90 47 1d | 6f 73 69 74 69 6f 6e 8e |ecomp.G.|osition.|
|000005b0| 8d 92 01 ac 6f 56 39 31 | 8e 8e a1 8d 8d 8d 8d 91 |....oV91|........|
|000005c0| 0f 00 00 da 66 66 74 2c | 91 05 3f fd 69 66 66 74 |....fft,|..?.ifft|
|000005d0| 8e 8e 8e 8d 92 00 e2 2b | 56 b2 46 96 ff 2a aa 61 |.......+|V.F..*.a|
|000005e0| 73 74 9b 03 55 55 46 93 | 6f 75 72 69 65 72 98 54 |st..UUF.|ourier.T|
|000005f0| 93 72 61 6e 73 66 6f 72 | 6d 8e 8d 92 01 ac 6f 56 |.ransfor|m.....oV|
|00000600| 39 33 8e 8e 8e 8e 8e a0 | 00 ab e8 d4 8d 91 0f 00 |93......|........|
|00000610| 00 54 91 ff 2a aa 6f 96 | 03 55 55 75 73 65 93 74 |.T..*.o.|.UUuse.t|
|00000620| 68 65 73 65 93 72 6f 75 | 74 69 6e 65 73 93 75 73 |hese.rou|tines.us|
|00000630| 65 93 74 68 65 93 69 6e | 63 6c 75 64 65 93 73 74 |e.the.in|clude.st|
|00000640| 61 74 65 6d 65 6e 90 b8 | e3 74 73 8e 9f 14 00 00 |atemen..|.ts.....|
|00000650| 8d da 23 69 6e 63 6c 75 | 64 65 91 05 3f fd 22 6d |..#inclu|de..?."m|
|00000660| 61 74 72 69 78 2e 68 22 | 8e 9f 0c 00 00 8d 23 69 |atrix.h"|......#i|
|00000670| 6e 63 6c 75 64 65 91 05 | 3f fd 22 6d 61 74 72 69 |nclude..|?."matri|
|00000680| 78 32 2e 68 22 8e 8e 9f | 19 4b 81 8d 92 00 e5 bc |x2.h"...|.K......|
|00000690| 37 b2 37 30 8e 8e 8c 8b | 00 00 00 47 00 00 00 00 |7.70....|...G....|
|000006a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000006c0| 00 00 00 2a a0 02 a9 52 | b3 8d a0 fd 62 ad 4d 8d |...*...R|....b.M.|
|000006d0| 8d 92 01 cb 78 6e f3 33 | 70 ae 30 4a 00 0a 00 00 |....xn.3|p.0J....|
|000006e0| 00 0a 00 00 00 06 63 6d | 73 6c 31 30 de 37 31 8e |......cm|sl10.71.|
|000006f0| 8e a0 02 84 07 32 8d a0 | fd 98 0a 15 8d 8d f3 30 |.....2..|.......0|
|00000700| 6d 23 bd 52 00 0a 00 00 | 00 0a 00 00 00 06 63 6d |m#.R....|......cm|
|00000710| 73 73 31 30 db 4e 41 4d | 45 8e 8e a9 10 44 95 8d |ss10.NAM|E....D..|
|00000720| 91 0f 00 00 da 42 4b 50 | 66 61 63 74 6f 72 2c 91 |.....BKP|factor,.|
|00000730| 05 3f fd 42 4b 50 73 6f | 6c 76 65 96 03 55 55 b2 |.?.BKPso|lve..UU.|
|00000740| 7b 93 42 75 6e 63 95 b8 | e3 68 7b 4b 61 75 66 6d |{.Bunc..|.h{Kaufm|
|00000750| 61 6e 7b 50 93 61 72 6c | 65 74 74 96 03 55 55 73 |an{P.arl|ett..UUs|
|00000760| 79 6d 6d 65 74 72 69 63 | 93 69 6e 64 65 0c 6e 69 |ymmetric|.inde.ni|
|00000770| 74 65 93 66 61 63 74 6f | 72 69 73 65 93 61 6e 64 |te.facto|rise.and|
|00000780| 93 73 6f 6c 76 90 b8 e3 | 65 8e 9f 13 44 95 8d 8d |.solv...|e...D...|
|00000790| db 53 59 4e 4f 50 53 49 | 53 8e 8e 9f 17 67 bf 8d |.SYNOPSI|S....g..|
|000007a0| da 23 69 6e 63 6c 75 64 | 65 91 05 3f fd 22 6d 61 |.#includ|e..?."ma|
|000007b0| 74 72 69 78 2e 68 22 8e | a4 0c 00 00 8d 23 69 6e |trix.h".|.....#in|
|000007c0| 63 6c 75 64 65 91 05 3f | fd 22 6d 61 74 72 69 78 |clude..?|."matrix|
|000007d0| 32 2e 68 22 8e a1 8d 4d | 41 54 91 1a 3f f1 2a 42 |2.h"...M|AT..?.*B|
|000007e0| 4b 50 66 61 63 74 6f 72 | 28 41 2c 70 69 76 6f 74 |KPfactor|(A,pivot|
|000007f0| 2c 62 6c 6f 63 6b 73 29 | 8e a1 8d 4d 41 54 91 1a |,blocks)|...MAT..|
|00000800| 3f f1 2a 41 3b 8e a1 8d | 50 45 52 4d 91 14 ff f4 |?.*A;...|PERM....|
|00000810| 2a 70 69 76 6f 74 2c 91 | 05 3f fd 2a 62 6c 6f 63 |*pivot,.|.?.*bloc|
|00000820| 6b 73 3b 8e 9f 19 30 89 | 8d 56 45 43 91 1a 3f f1 |ks;...0.|.VEC..?.|
|00000830| 2a 42 4b 50 73 6f 6c 76 | 65 28 41 2c 70 69 76 6f |*BKPsolv|e(A,pivo|
|00000840| 74 2c 62 6c 6f 63 6b 73 | 2c 62 2c 78 29 8e a1 8d |t,blocks|,b,x)...|
|00000850| 4d 41 54 91 1a 3f f1 2a | 41 3b 8e a1 8d 50 45 52 |MAT..?.*|A;...PER|
|00000860| 4d 91 14 ff f4 2a 70 69 | 76 6f 74 2c 91 05 3f fd |M....*pi|vot,..?.|
|00000870| 2a 62 6c 6f 63 6b 73 3b | 8e a1 8d 56 45 43 91 1a |*blocks;|...VEC..|
|00000880| 3f f1 2a 62 2c 91 05 3f | fd 2a 78 3b 8e 9f 1c 30 |?.*b,..?|.*x;...0|
|00000890| 89 8d 8d db 44 45 53 43 | 52 49 50 54 49 4f 4e 8e |....DESC|RIPTION.|
|000008a0| 8e a6 8d 91 0f 00 00 b2 | 54 68 65 96 03 55 55 72 |........|The..UUr|
|000008b0| 6f 75 74 69 6e 65 93 da | 42 4b 50 66 61 63 74 6f |outine..|BKPfacto|
|000008c0| 72 28 29 93 b2 66 6f 72 | 6d 73 93 f3 31 fd 00 27 |r()..for|ms..1..'|
|000008d0| 3a 00 0a 00 00 00 0a 00 | 00 00 06 63 6d 74 69 31 |:.......|...cmti1|
|000008e0| 30 dc 69 6e 91 03 93 e7 | 73 69 74 75 91 04 19 b8 |0.in....|situ....|
|000008f0| b2 61 93 73 79 6d 6d 65 | 74 72 69 63 93 69 6e 64 |.a.symme|tric.ind|
|00000900| 65 0c 6e 69 74 65 93 66 | 61 63 74 6f 72 69 73 61 |e.nite.f|actorisa|
|00000910| 74 69 6f 6e 93 6f 66 93 | 74 68 65 93 6d 61 74 72 |tion.of.|the.matr|
|00000920| 69 78 93 da 41 93 b2 6f | 66 93 74 68 65 93 66 6f |ix..A..o|f.the.fo|
|00000930| 72 6d 8e a4 14 eb f4 8d | 92 00 c3 23 ca b5 50 8d |rm......|...#..P.|
|00000940| 91 01 63 8f 9f fb de ff | b4 54 8e 91 07 b0 45 b5 |..c.....|.T....E.|
|00000950| 41 50 91 04 2a a7 b2 3d | 91 02 c7 18 b5 4d 96 01 |AP..*..=|.....M..|
|00000960| 17 1b 44 90 47 1d 4d 8d | 93 9f fb de ff b4 54 8e |..D.G.M.|......T.|
|00000970| 8e a1 8d b2 77 68 65 72 | 65 96 04 83 91 b5 50 91 |....wher|e.....P.|
|00000980| 05 e7 20 b2 69 73 93 61 | 93 70 90 47 1d 65 72 6d |.. .is.a|.p.G.erm|
|00000990| 9a b8 e3 75 74 61 74 69 | 6f 6e 93 6d 61 74 72 69 |...utati|on.matri|
|000009a0| 78 2c 91 04 cf 20 b5 4d | 91 05 9a ac b2 69 73 93 |x,... .M|.....is.|
|000009b0| 6c 6f 98 77 98 65 72 93 | 74 72 69 61 6e 67 75 6c |lo.w.er.|triangul|
|000009c0| 61 72 2c 91 04 cf 20 61 | 6e 64 93 b5 44 91 04 ca |ar,... a|nd..D...|
|000009d0| ae b2 69 73 93 62 6c 6f | 90 47 1d 63 98 6b 93 64 |..is.blo|.G.c.k.d|
|000009e0| 69 61 67 6f 6e 61 6c 2c | 91 04 cf 20 77 69 74 68 |iagonal,|... with|
|000009f0| 93 31 96 03 02 5c f3 0d | 21 22 2c 9a 00 0a 00 00 |.1...\..|!",.....|
|00000a00| 00 0a 00 00 00 06 63 6d | 73 79 31 30 b8 02 93 b2 |......cm|sy10....|
|00000a10| 31 96 04 83 91 6f 72 93 | 32 96 03 02 5c b8 02 93 |1....or.|2...\...|
|00000a20| b2 32 8e a4 0c 00 00 8d | 62 6c 6f 90 47 1d 63 9a |.2......|blo.G.c.|
|00000a30| b8 e3 6b 73 2e 91 07 ed | 0a 54 68 65 96 04 7e 6b |..ks....|.The..~k|
|00000a40| 6d 61 74 72 69 78 93 b5 | 50 91 05 e1 fa b2 69 73 |matrix..|P.....is|
|00000a50| 93 72 65 70 72 65 73 65 | 6e 98 74 65 64 93 62 98 |.represe|n.ted.b.|
|00000a60| 79 93 74 68 65 93 70 90 | 47 1d 65 72 6d 98 75 74 |y.the.p.|G.erm.ut|
|00000a70| 61 74 69 6f 6e 93 da 70 | 69 76 6f 74 93 b2 61 6e |ation..p|ivot..an|
|00000a80| 64 93 b5 44 8d 9f 01 7f | ff b4 69 69 8e 91 0a a7 |d..D....|..ii....|
|00000a90| 03 b2 69 73 93 61 93 31 | 96 02 fe ee b8 02 93 b2 |..is.a.1|........|
|00000aa0| 31 96 04 7e 6b 62 6c 6f | 90 47 1d 63 98 6b 93 69 |1..~kblo|.G.c.k.i|
|00000ab0| 66 93 61 6e 64 93 6f 6e | 6c 79 93 69 66 8e a1 8d |f.and.on|ly.if...|
|00000ac0| da 62 6c 6f 63 6b 73 2d | 3e 70 65 5b 69 5d 96 05 |.blocks-|>pe[i]..|
|00000ad0| 3f fd 3d 3d 93 69 b2 3b | 91 04 1b 11 6f 74 68 65 |?.==.i.;|....othe|
|00000ae0| 72 77 69 73 65 96 03 d9 | 28 da 62 6c 6f 63 6b 73 |rwise...|(.blocks|
|00000af0| 2d 3e 70 65 5b 69 5d 93 | b2 69 73 93 74 68 65 93 |->pe[i].|.is.the.|
|00000b00| 69 6e 64 65 78 93 6f 66 | 93 74 68 65 93 6f 74 68 |index.of|.the.oth|
|00000b10| 65 72 93 72 6f 9a b8 e3 | 77 2f 63 6f 6c 75 6d 6e |er.ro...|w/column|
|00000b20| 93 69 6e 93 74 68 65 93 | 32 96 02 90 c1 b8 02 93 |.in.the.|2.......|
|00000b30| b2 32 91 03 d9 28 62 6c | 6f 90 47 1d 63 98 6b 2e |.2...(bl|o.G.c.k.|
|00000b40| 8e a1 8d 41 66 74 65 72 | 96 03 e9 c8 74 68 65 93 |...After|....the.|
|00000b50| 72 6f 75 74 69 6e 65 93 | 74 68 65 93 b5 44 91 04 |routine.|the..D..|
|00000b60| 30 e5 b2 61 6e 64 93 b5 | 4d 91 05 00 e3 b2 66 61 |0..and..|M.....fa|
|00000b70| 63 74 6f 72 73 93 61 72 | 65 93 73 74 6f 72 65 64 |ctors.ar|e.stored|
|00000b80| 93 69 6e 93 da 41 91 03 | e9 a1 b2 69 6e 93 63 6f |.in..A..|...in.co|
|00000b90| 6d 70 61 63 74 93 66 6f | 72 6d 2e 91 06 2f 1f 54 |mpact.fo|rm.../.T|
|00000ba0| 68 69 73 93 61 95 b8 e3 | 76 93 6f 69 64 73 96 03 |his.a...|v.oids..|
|00000bb0| e9 c8 74 68 65 93 72 65 | 71 75 69 72 65 6d 65 6e |..the.re|quiremen|
|00000bc0| 90 b8 e3 74 93 66 6f 72 | 8e a1 8d 61 64 64 69 74 |...t.for|...addit|
|00000bd0| 69 6f 6e 61 6c 96 03 55 | 55 76 90 b8 e3 65 63 74 |ional..U|Uv...ect|
|00000be0| 6f 72 73 93 6f 72 93 6d | 61 74 72 69 63 65 73 93 |ors.or.m|atrices.|
|00000bf0| 66 6f 72 93 73 74 6f 72 | 61 67 65 2e 8e a6 8d 91 |for.stor|age.....|
|00000c00| 0f 00 00 4e 6f 74 65 96 | 02 38 e3 74 68 61 74 93 |...Note.|.8.that.|
|00000c10| da 70 69 76 6f 74 93 b2 | 61 6e 64 93 da 62 6c 6f |.pivot..|and..blo|
|00000c20| 63 6b 73 93 b2 6d 90 b8 | e3 75 73 74 93 62 9a 47 |cks..m..|.ust.b.G|
|00000c30| 1d 6f 74 68 93 62 98 65 | 93 6e 6f 6e 7b 4e 55 4c |.oth.b.e|.non{NUL|
|00000c40| 4c 91 02 38 9b 61 6e 64 | 93 64 69 0b 65 72 65 6e |L..8.and|.di.eren|
|00000c50| 90 b8 e3 74 93 66 6f 72 | 93 62 98 6f 74 68 93 da |...t.for|.b.oth..|
|00000c60| 42 4b 50 66 61 63 74 6f | 72 28 29 93 b2 61 6e 64 |BKPfacto|r()..and|
|00000c70| 93 da 42 4b 50 73 6f 6c | 76 65 28 29 b2 2e 8e a6 |..BKPsol|ve()....|
|00000c80| 8d 91 0f 00 00 54 68 65 | 96 03 f4 37 72 6f 75 74 |.....The|...7rout|
|00000c90| 69 6e 65 93 da 42 4b 50 | 73 6f 6c 76 65 28 29 93 |ine..BKP|solve().|
|00000ca0| b2 73 6f 6c 76 9a b8 e3 | 65 73 93 74 68 65 93 65 |.solv...|es.the.e|
|00000cb0| 71 75 61 74 69 6f 6e 93 | b5 41 78 96 03 cf e4 b2 |quation.|.Ax.....|
|00000cc0| 3d 93 b5 62 96 03 f4 37 | b2 66 6f 72 93 b5 78 b2 |=..b...7|.for..x.|
|00000cd0| 2e 91 06 4e 6c 54 68 65 | 93 73 6f 6c 76 98 65 93 |...NlThe|.solv.e.|
|00000ce0| 72 6f 75 74 69 6e 65 93 | da 42 4b 50 73 6f 6c 76 |routine.|.BKPsolv|
|00000cf0| 65 28 29 93 b2 69 73 93 | 64 65 73 69 67 6e 65 64 |e()..is.|designed|
|00000d00| 8e a1 8d 73 70 9a 47 1d | 65 63 69 0c 63 61 6c 6c |...sp.G.|eci.call|
|00000d10| 79 96 03 46 ec 74 6f 93 | 77 90 b8 e3 6f 72 6b 93 |y..F.to.|w...ork.|
|00000d20| 77 69 74 68 93 da 42 4b | 50 66 61 63 74 6f 72 28 |with..BK|Pfactor(|
|00000d30| 29 93 b2 61 73 93 74 68 | 65 79 93 6f 70 98 65 72 |)..as.th|ey.op.er|
|00000d40| 61 74 65 93 6f 6e 93 74 | 68 65 93 73 61 6d 65 93 |ate.on.t|he.same.|
|00000d50| 63 6f 6d 70 61 63 74 93 | 73 74 6f 72 61 67 65 93 |compact.|storage.|
|00000d60| 73 63 90 b8 e3 68 65 6d | 65 2e 91 04 6c f9 4e 6f |sc...hem|e...l.No|
|00000d70| 74 65 93 74 68 61 74 93 | 74 68 65 8e a1 8d 66 61 |te.that.|the...fa|
|00000d80| 63 74 6f 72 69 73 61 74 | 69 6f 6e 96 03 da de 6d |ctorisat|ion....m|
|00000d90| 61 9a b8 e3 79 93 73 75 | 63 63 65 65 64 93 77 68 |a...y.su|cceed.wh|
|00000da0| 65 6e 93 74 68 65 93 6d | 61 74 72 69 78 93 da 41 |en.the.m|atrix..A|
|00000db0| 91 03 da bc b2 70 61 73 | 73 65 64 93 69 73 93 73 |.....pas|sed.is.s|
|00000dc0| 69 6e 67 75 6c 61 72 2c | 91 03 fc 40 61 6e 64 93 |ingular,|...@and.|
|00000dd0| 74 68 61 74 93 74 68 65 | 93 73 6f 6c 76 98 65 93 |that.the|.solv.e.|
|00000de0| 72 6f 75 74 69 6e 65 93 | 6d 61 98 79 93 74 68 65 |routine.|ma.y.the|
|00000df0| 6e 93 66 61 69 6c 2c 8e | a1 8d 72 61 69 73 69 6e |n.fail,.|..raisin|
|00000e00| 67 96 03 82 84 61 6e 93 | da 45 5f 53 49 4e 47 91 |g....an.|.E_SING.|
|00000e10| 03 82 78 b2 65 72 72 6f | 72 2e 91 04 f9 54 54 68 |..x.erro|r....TTh|
|00000e20| 65 93 73 6f 6c 76 9a b8 | e3 65 93 72 6f 75 74 69 |e.solv..|.e.routi|
|00000e30| 6e 65 93 6d 61 98 79 93 | 62 9a 47 1d 65 93 75 73 |ne.ma.y.|b.G.e.us|
|00000e40| 65 64 93 dc 69 6e 91 03 | bd 79 73 69 74 75 91 04 |ed..in..|.ysitu..|
|00000e50| 46 e7 b2 77 69 74 68 93 | da 62 96 05 3f fd 3d 3d |F..with.|.b..?.==|
|00000e60| 93 78 b2 2e 91 04 f9 54 | 49 66 96 03 82 84 da 78 |.x.....T|If.....x|
|00000e70| 93 b2 69 73 93 4e 55 4c | 4c 91 03 82 78 6f 72 93 |..is.NUL|L...xor.|
|00000e80| 74 6f 98 6f 93 73 6d 61 | 6c 6c 93 74 6f 8e a1 8d |to.o.sma|ll.to...|
|00000e90| 68 6f 6c 64 96 03 33 da | 74 68 65 93 72 65 73 75 |hold..3.|the.resu|
|00000ea0| 6c 74 2c 91 03 3a 8c 74 | 68 65 6e 93 61 93 6e 65 |lt,..:.t|hen.a.ne|
|00000eb0| 77 93 76 90 b8 e3 65 63 | 74 6f 72 93 69 73 93 63 |w.v...ec|tor.is.c|
|00000ec0| 72 65 61 74 65 64 93 6f | 66 93 74 68 65 93 61 70 |reated.o|f.the.ap|
|00000ed0| 70 72 6f 70 72 69 61 74 | 65 93 73 69 7a 65 93 66 |propriat|e.size.f|
|00000ee0| 6f 72 93 73 74 6f 72 69 | 6e 67 93 74 68 65 93 72 |or.stori|ng.the.r|
|00000ef0| 65 73 75 6c 74 2e 91 04 | 66 9e 49 6e 93 65 69 74 |esult...|f.In.eit|
|00000f00| 68 65 72 93 63 61 73 65 | 93 74 68 65 8e a1 8d 72 |her.case|.the...r|
|00000f10| 65 73 75 6c 74 69 6e 67 | 96 03 55 55 73 6f 6c 75 |esulting|..UUsolu|
|00000f20| 74 69 6f 6e 93 76 90 b8 | e3 65 63 74 6f 72 93 69 |tion.v..|.ector.i|
|00000f30| 73 93 72 65 74 75 72 6e | 65 64 2e 8e a6 8d 91 0f |s.return|ed......|
|00000f40| 00 00 54 68 69 73 96 04 | 32 65 66 61 63 74 6f 72 |..This..|2efactor|
|00000f50| 69 73 61 74 69 6f 6e 93 | 72 6f 75 74 69 6e 65 2c |isation.|routine,|
|00000f60| 91 04 69 a9 61 6e 64 93 | 74 68 65 93 61 63 63 6f |..i.and.|the.acco|
|00000f70| 6d 70 61 6e 9a b8 e3 79 | 69 6e 67 93 73 6f 6c 76 |mpan...y|ing.solv|
|00000f80| 98 65 93 72 6f 75 74 69 | 6e 65 93 61 72 65 93 64 |.e.routi|ne.are.d|
|00000f90| 65 72 69 76 98 65 64 93 | 66 72 6f 6d 93 5c 44 65 |eriv.ed.|from.\De|
|00000fa0| 63 6f 6d 70 90 47 1d 6f | 73 69 74 69 6f 6e 93 6f |comp.G.o|sition.o|
|00000fb0| 66 93 61 8e a1 8d 53 79 | 6d 6d 65 74 72 69 63 96 |f.a...Sy|mmetric.|
|00000fc0| 03 55 55 4d 61 74 72 69 | 78 22 93 62 9a b8 e3 79 |.UUMatri|x".b...y|
|00000fd0| 93 4a 2e 93 42 75 6e 63 | 98 68 2c 93 4c 2e 93 4b |.J..Bunc|.h,.L..K|
|00000fe0| 61 75 66 6d 61 6e 93 61 | 6e 64 93 42 2e 93 50 98 |aufman.a|nd.B..P.|
|00000ff0| 61 72 6c 65 74 74 2c 93 | dc 4e 75 6d 65 72 69 73 |arlett,.|.Numeris|
|00001000| 63 68 65 91 03 93 e7 4d | 61 74 68 65 6d 61 74 69 |che....M|athemati|
|00001010| 6b 91 04 68 e2 f3 32 1a | f2 22 56 00 0a 00 00 00 |k..h..2.|."V.....|
|00001020| 0a 00 00 00 06 63 6d 62 | 78 31 30 dd 32 37 b2 2c |.....cmb|x10.27.,|
|00001030| 93 39 35 7b 31 30 39 93 | 28 31 39 37 36 29 2e 8e |.95{109.|(1976)..|
|00001040| a6 8d 91 0f 00 00 45 72 | 72 6f 72 73 96 03 70 97 |......Er|rors..p.|
|00001050| 77 69 6c 6c 93 62 90 47 | 1d 65 93 72 61 69 73 65 |will.b.G|.e.raise|
|00001060| 64 93 69 66 93 da 41 9b | 03 70 90 b2 6f 72 93 da |d.if..A.|.p..or..|
|00001070| 70 69 76 6f 74 93 b2 6f | 72 93 da 62 6c 6f 63 6b |pivot..o|r..block|
|00001080| 73 93 b2 61 72 65 93 4e | 55 4c 4c 2c 93 6f 72 93 |s..are.N|ULL,.or.|
|00001090| 69 66 93 da 41 98 b2 69 | 73 93 6e 6f 74 93 73 71 |if..A..i|s.not.sq|
|000010a0| 75 61 72 65 2c 91 03 77 | 68 6f 72 93 69 66 93 74 |uare,..w|hor.if.t|
|000010b0| 68 65 93 73 69 7a 65 73 | 93 6f 66 93 da 41 b2 2c |he.sizes|.of..A.,|
|000010c0| 93 da 70 69 76 6f 74 8e | a1 8d b2 6f 72 96 03 55 |..pivot.|...or..U|
|000010d0| 55 da 62 6c 6f 63 6b 73 | 93 b2 61 72 65 93 6e 6f |U.blocks|..are.no|
|000010e0| 74 93 63 6f 6d 70 61 74 | 69 62 6c 65 2e 8e 9f 13 |t.compat|ible....|
|000010f0| 44 95 8d 8d db 45 58 41 | 4d 50 4c 45 8e 8e 9f 19 |D....EXA|MPLE....|
|00001100| 30 89 8d da 4d 41 54 91 | 14 ff f4 2a 41 3b 8e a1 |0...MAT.|...*A;..|
|00001110| 8d 50 45 52 4d 91 0f bf | f7 2a 70 69 76 6f 74 2c |.PERM...|.*pivot,|
|00001120| 91 05 3f fd 2a 62 6c 6f | 63 6b 73 3b 8e a1 8d 56 |..?.*blo|cks;...V|
|00001130| 45 43 91 14 ff f4 2a 78 | 2c 91 05 3f fd 2a 62 3b |EC....*x|,..?.*b;|
|00001140| 8e a1 8d 91 0a 7f fa 2e | 2e 2e 2e 2e 2e 8e a1 8d |........|........|
|00001150| 41 96 05 3f fd 3d 93 69 | 6e 5f 6d 61 74 28 4d 4e |A..?.=.i|n_mat(MN|
|00001160| 55 4c 4c 29 3b 8e a1 8d | 62 96 05 3f fd 3d 93 69 |ULL);...|b..?.=.i|
|00001170| 6e 5f 76 65 63 28 56 4e | 55 4c 4c 29 3b 8e a1 8d |n_vec(VN|ULL);...|
|00001180| 70 69 76 6f 74 96 05 3f | fd 3d 93 67 65 74 5f 70 |pivot..?|.=.get_p|
|00001190| 65 72 6d 28 41 2d 3e 6d | 29 3b 8e a1 8d 62 6c 6f |erm(A->m|);...blo|
|000011a0| 63 6b 73 96 05 3f fd 3d | 93 67 65 74 5f 70 65 72 |cks..?.=|.get_per|
|000011b0| 6d 28 41 2d 3e 6d 29 3b | 8e a1 8d 2f 2a 96 05 3f |m(A->m);|.../*..?|
|000011c0| fd 61 73 73 75 6d 69 6e | 67 93 41 93 73 79 6d 6d |.assumin|g.A.symm|
|000011d0| 65 74 72 69 63 93 2a 2f | 8e a1 8d 42 4b 50 66 61 |etric.*/|...BKPfa|
|000011e0| 63 74 6f 72 28 41 2c 70 | 69 76 6f 74 2c 62 6c 6f |ctor(A,p|ivot,blo|
|000011f0| 63 6b 73 29 3b 8e a1 8d | 78 96 05 3f fd 3d 93 42 |cks);...|x..?.=.B|
|00001200| 4b 50 73 6f 6c 76 65 28 | 41 2c 70 69 76 6f 74 2c |KPsolve(|A,pivot,|
|00001210| 62 6c 6f 63 6b 73 2c 62 | 2c 56 4e 55 4c 4c 29 3b |blocks,b|,VNULL);|
|00001220| 8e 8e 8e 8c 8b 00 00 00 | 48 00 00 00 00 00 00 00 |........|H.......|
|00001230| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001240| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 08 |........|........|
|00001250| cb a0 02 a9 52 b3 8d a0 | fd 62 ad 4d 8d 8d b2 37 |....R...|.b.M...7|
|00001260| 32 92 00 99 14 c1 de 43 | 48 41 50 54 45 52 96 03 |2......C|HAPTER..|
|00001270| 55 55 35 2e 91 07 c7 1c | 44 45 4e 53 45 93 4d 41 |UU5.....|DENSE.MA|
|00001280| 9b ff 2a aa 54 52 49 58 | 93 46 91 fe e3 8d 41 90 |..*.TRIX|.F....A.|
|00001290| b8 e3 43 54 4f 52 49 53 | 41 98 54 49 4f 4e 93 4f |..CTORIS|A.TION.O|
|000012a0| 50 45 52 41 98 54 49 4f | 4e 53 8e 8e a0 02 84 07 |PERA.TIO|NS......|
|000012b0| 32 8d a0 fd 98 0a 15 8d | 8d db 53 45 45 91 03 55 |2.......|..SEE..U|
|000012c0| 55 41 4c 53 4f 8e 8e 9f | 10 44 95 8d 91 0f 00 00 |UALSO...|.D......|
|000012d0| da 43 48 66 61 63 74 6f | 72 28 29 96 03 55 55 b2 |.CHfacto|r()..UU.|
|000012e0| 61 6e 64 93 da 43 48 73 | 6f 6c 76 65 28 29 8e 9f |and..CHs|olve()..|
|000012f0| 13 44 95 8d 8d db 53 4f | 55 52 43 45 91 03 55 55 |.D....SO|URCE..UU|
|00001300| 46 49 4c 45 3a 91 14 00 | 02 da 42 4b 50 66 61 63 |FILE:...|..BKPfac|
|00001310| 74 6f 72 2e 63 8e 8e 8e | 8e 8c 8b 00 00 00 49 00 |tor.c...|......I.|
|00001320| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001330| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001340| 00 00 00 00 00 14 58 a0 | 02 a9 52 b3 8d a0 fd 62 |......X.|..R....b|
|00001350| ad 4d 8d 8d 92 01 cb 78 | 6e de 37 33 8e 8e a0 02 |.M.....x|n.73....|
|00001360| 84 07 32 8d a0 fd 98 0a | 15 8d 8d db 4e 41 4d 45 |..2.....|....NAME|
|00001370| 8e 8e 9f 10 44 95 8d 91 | 0f 00 00 da 43 48 66 61 |....D...|....CHfa|
|00001380| 63 74 6f 72 2c 96 05 3f | fd 4d 43 48 66 61 63 74 |ctor,..?|.MCHfact|
|00001390| 6f 72 2c 93 43 48 73 6f | 6c 76 65 2c 93 4c 44 4c |or,.CHso|lve,.LDL|
|000013a0| 66 61 63 74 6f 72 2c 93 | 4c 44 4c 73 6f 6c 76 65 |factor,.|LDLsolve|
|000013b0| 96 03 55 55 b2 7b 93 43 | 68 6f 6c 65 73 6b 79 93 |..UU.{.C|holesky.|
|000013c0| 66 61 63 74 6f 72 93 61 | 6e 64 93 73 6f 6c 76 90 |factor.a|nd.solv.|
|000013d0| b8 e3 65 8e 9f 13 44 95 | 8d 8d db 53 59 4e 4f 50 |..e...D.|...SYNOP|
|000013e0| 53 49 53 8e 8e 9f 18 9d | 9e 8d da 23 69 6e 63 6c |SIS.....|...#incl|
|000013f0| 75 64 65 91 05 3f fd 22 | 6d 61 74 72 69 78 2e 68 |ude..?."|matrix.h|
|00001400| 22 8e a4 0c 00 00 8d 23 | 69 6e 63 6c 75 64 65 91 |"......#|include.|
|00001410| 05 3f fd 22 6d 61 74 72 | 69 78 32 2e 68 22 8e a1 |.?."matr|ix2.h"..|
|00001420| 8d 4d 41 54 91 1a 3f f1 | 2a 43 48 66 61 63 74 6f |.MAT..?.|*CHfacto|
|00001430| 72 28 41 29 8e a1 8d 4d | 41 54 91 1a 3f f1 2a 41 |r(A)...M|AT..?.*A|
|00001440| 3b 8e a9 1a ca 23 8d 4d | 41 54 91 1a 3f f1 2a 4d |;....#.M|AT..?.*M|
|00001450| 43 48 66 61 63 74 6f 72 | 28 41 2c 74 6f 6c 29 8e |CHfactor|(A,tol).|
|00001460| a1 8d 4d 41 54 91 1a 3f | f1 2a 41 3b 8e a1 8d 64 |..MAT..?|.*A;...d|
|00001470| 6f 75 62 6c 65 91 0a 7f | fa 74 6f 6c 3b 8e a6 8d |ouble...|.tol;...|
|00001480| 56 45 43 91 1a 3f f1 2a | 43 48 73 6f 6c 76 65 28 |VEC..?.*|CHsolve(|
|00001490| 41 2c 62 2c 78 29 8e a1 | 8d 4d 41 54 91 1a 3f f1 |A,b,x)..|.MAT..?.|
|000014a0| 2a 41 3b 8e a1 8d 56 45 | 43 91 1a 3f f1 2a 62 2c |*A;...VE|C..?.*b,|
|000014b0| 91 05 3f fd 2a 78 3b 8e | a6 8d 4d 41 54 91 1a 3f |..?.*x;.|..MAT..?|
|000014c0| f1 2a 4c 44 4c 66 61 63 | 74 6f 72 28 41 29 8e a1 |.*LDLfac|tor(A)..|
|000014d0| 8d 4d 41 54 91 1a 3f f1 | 2a 41 3b 8e a6 8d 56 45 |.MAT..?.|*A;...VE|
|000014e0| 43 91 1a 3f f1 2a 4c 44 | 4c 73 6f 6c 76 65 28 41 |C..?.*LD|Lsolve(A|
|000014f0| 2c 62 2c 78 29 8e a1 8d | 4d 41 54 91 1a 3f f1 2a |,b,x)...|MAT..?.*|
|00001500| 41 3b 8e a1 8d 56 45 43 | 91 1a 3f f1 2a 62 2c 91 |A;...VEC|..?.*b,.|
|00001510| 05 3f fd 2a 78 3b 8e 9f | 1d ca 23 8d 8d db 44 45 |.?.*x;..|..#...DE|
|00001520| 53 43 52 49 50 54 49 4f | 4e 8e 8e a9 10 44 95 8d |SCRIPTIO|N....D..|
|00001530| 91 0f 00 00 b2 42 6f 74 | 68 96 04 a3 26 da 43 48 |.....Bot|h...&.CH|
|00001540| 66 61 63 74 6f 72 28 29 | 93 b2 61 6e 64 93 da 4c |factor()|..and..L|
|00001550| 44 4c 66 61 63 74 6f 72 | 28 29 93 b2 66 61 63 74 |DLfactor|()..fact|
|00001560| 6f 72 93 74 68 65 93 6d | 61 74 72 69 78 93 da 41 |or.the.m|atrix..A|
|00001570| 91 04 a2 d0 dc 69 6e 91 | 04 c7 04 73 69 74 75 91 |.....in.|...situ.|
|00001580| 05 67 89 b2 61 6e 64 93 | 72 65 74 75 72 6e 73 93 |.g..and.|returns.|
|00001590| 74 68 65 93 66 61 63 74 | 6f 72 65 64 93 6d 61 74 |the.fact|ored.mat|
|000015a0| 72 69 78 93 28 69 6e 8e | a1 8d 63 6f 6d 70 61 63 |rix.(in.|..compac|
|000015b0| 74 96 02 5d cd 66 6f 72 | 6d 29 2e 91 04 1f 45 54 |t..].for|m)....ET|
|000015c0| 68 65 93 43 68 6f 6c 65 | 73 6b 79 93 66 61 63 74 |he.Chole|sky.fact|
|000015d0| 6f 72 69 73 61 74 69 6f | 6e 93 72 6f 75 74 69 6e |orisatio|n.routin|
|000015e0| 65 93 61 6e 64 93 74 68 | 65 93 b5 4c 44 9a 47 1d |e.and.th|e..LD.G.|
|000015f0| 4c 8d 9f fc 5e ff b4 54 | 8e 91 08 aa 83 b2 72 6f |L...^..T|......ro|
|00001600| 75 74 69 6e 65 73 93 62 | 98 6f 74 68 93 75 73 65 |utines.b|.oth.use|
|00001610| 93 6f 6e 6c 79 93 74 68 | 65 93 6c 6f 95 b8 e3 77 |.only.th|e.lo...w|
|00001620| 93 65 72 91 02 5d cd 74 | 72 69 61 6e 67 75 6c 61 |.er..].t|riangula|
|00001630| 72 8e a1 8d 70 61 72 74 | 96 03 55 55 6f 66 93 da |r...part|..UUof..|
|00001640| 41 b2 2c 93 62 75 74 93 | 74 68 65 93 43 68 6f 6c |A.,.but.|the.Chol|
|00001650| 65 73 6b 79 93 66 61 63 | 74 6f 72 69 73 61 74 69 |esky.fac|torisati|
|00001660| 6f 6e 93 72 6f 75 74 69 | 6e 65 93 0c 6c 6c 73 93 |on.routi|ne..lls.|
|00001670| 74 68 65 93 75 70 70 90 | 47 1d 65 72 93 74 72 69 |the.upp.|G.er.tri|
|00001680| 61 6e 67 75 6c 61 72 93 | 70 61 72 74 93 6f 66 93 |angular.|part.of.|
|00001690| da 41 93 b2 61 6c 73 6f | 2e 8e a6 8d 91 0f 00 00 |.A..also|........|
|000016a0| 54 68 65 73 65 96 04 34 | d1 72 6f 75 74 69 6e 65 |These..4|.routine|
|000016b0| 73 93 72 65 71 75 69 72 | 65 93 74 68 61 74 93 da |s.requir|e.that..|
|000016c0| 41 91 04 34 97 b2 69 73 | 93 73 71 75 61 72 65 2e |A..4..is|.square.|
|000016d0| 91 07 10 3a 54 68 65 93 | 43 68 6f 6c 65 73 6b 79 |...:The.|Cholesky|
|000016e0| 93 66 61 63 74 6f 72 69 | 73 61 74 69 6f 6e 2c 9b |.factori|sation,.|
|000016f0| 04 6c af 69 6e 93 70 61 | 72 74 69 63 75 6c 61 72 |.l.in.pa|rticular|
|00001700| 2c 98 72 65 71 75 69 72 | 65 73 93 74 68 61 74 93 |,.requir|es.that.|
|00001710| da 41 91 04 34 97 b2 62 | 90 47 1d 65 8e a1 8d 73 |.A..4..b|.G.e...s|
|00001720| 75 0e 63 69 65 6e 95 b8 | e3 74 6c 79 9b 03 54 6c |u.cien..|.tly..Tl|
|00001730| 70 90 47 1d 6f 73 69 74 | 69 76 93 65 98 64 65 0c |p.G.osit|iv.e.de.|
|00001740| 6e 69 74 65 98 28 65 2e | 67 2e 98 6c 6f 93 77 93 |nite.(e.|g..lo.w.|
|00001750| 65 73 74 98 65 69 67 65 | 6e 93 76 91 ff 71 c7 61 |est.eige|n.v..q.a|
|00001760| 6c 75 65 98 6f 66 98 da | 41 98 b2 69 73 98 61 74 |lue.of..|A..is.at|
|00001770| 98 6c 65 61 73 74 98 6d | 61 63 93 68 69 6e 65 98 |.least.m|ac.hine.|
|00001780| 65 70 73 69 6c 6f 6e 98 | 61 93 77 93 61 93 79 98 |epsilon.|a.w.a.y.|
|00001790| 66 72 6f 6d 98 7a 65 72 | 6f 29 2e 91 04 71 79 49 |from.zer|o)...qyI|
|000017a0| 66 98 6e 6f 6e 7b 8e a1 | 8d 70 9a 47 1d 6f 73 69 |f.non{..|.p.G.osi|
|000017b0| 74 69 76 90 b8 e3 65 96 | 03 93 43 64 65 0c 6e 69 |tiv...e.|..Cde.ni|
|000017c0| 74 65 6e 65 73 73 93 69 | 73 93 64 65 74 65 63 74 |teness.i|s.detect|
|000017d0| 65 64 93 64 75 72 69 6e | 67 93 66 61 63 74 6f 72 |ed.durin|g.factor|
|000017e0| 69 73 61 74 69 6f 6e 2c | 91 03 a2 be 74 68 65 6e |isation,|....then|
|000017f0| 93 61 6e 93 da 45 5f 50 | 4f 53 44 45 46 91 03 93 |.an..E_P|OSDEF...|
|00001800| 33 b2 65 72 72 6f 72 93 | 77 69 6c 6c 93 62 98 65 |3.error.|will.b.e|
|00001810| 93 72 61 69 73 65 64 2e | 91 05 2b 90 49 66 93 79 |.raised.|..+.If.y|
|00001820| 90 b8 e3 6f 75 93 77 69 | 73 68 93 74 6f 8e a1 8d |...ou.wi|sh.to...|
|00001830| 63 61 74 63 95 b8 e3 68 | 9b 03 cc fb 73 75 63 93 |catc...h|....suc.|
|00001840| 68 98 61 6e 98 65 72 72 | 6f 72 2c 91 03 ea e4 73 |h.an.err|or,....s|
|00001850| 65 65 98 69 6e 66 6f 72 | 6d 61 74 69 6f 6e 98 6f |ee.infor|mation.o|
|00001860| 6e 98 74 68 65 98 da 63 | 61 74 63 68 28 29 98 b2 |n.the..c|atch()..|
|00001870| 6d 61 63 72 6f 2e 91 05 | d8 b8 49 66 98 79 93 6f |macro...|..If.y.o|
|00001880| 75 72 98 6d 61 74 72 69 | 78 98 69 73 98 69 6e 64 |ur.matri|x.is.ind|
|00001890| 65 0c 6e 69 74 65 2c 91 | 03 ea e4 74 68 65 6e 98 |e.nite,.|...then.|
|000018a0| 69 74 98 77 93 6f 75 6c | 64 98 62 90 47 1d 65 8e |it.w.oul|d.b.G.e.|
|000018b0| a1 8d 62 90 47 1d 65 73 | 74 96 03 55 55 74 6f 93 |..b.G.es|t..UUto.|
|000018c0| 75 73 65 93 74 68 65 93 | da 42 4b 50 66 61 63 74 |use.the.|.BKPfact|
|000018d0| 6f 72 28 29 93 b2 61 6e | 64 93 da 42 4b 50 73 6f |or()..an|d..BKPso|
|000018e0| 6c 76 65 28 29 93 b2 72 | 6f 75 74 69 6e 65 73 2e |lve()..r|outines.|
|000018f0| 8e a6 8d 91 0f 00 00 54 | 68 65 96 04 c3 32 72 6f |.......T|he...2ro|
|00001900| 75 74 69 6e 65 93 da 4d | 43 48 66 61 63 74 6f 72 |utine..M|CHfactor|
|00001910| 28 29 93 b2 63 6f 6d 70 | 75 74 65 73 93 61 93 dc |()..comp|utes.a..|
|00001920| 6d 6f 96 ff 7d 27 64 69 | 0c 65 93 64 91 05 cb bb |mo..}'di|.e.d....|
|00001930| b2 43 68 6f 6c 65 73 6b | 79 96 04 c3 32 66 61 63 |.Cholesk|y...2fac|
|00001940| 74 6f 72 69 73 61 74 69 | 6f 6e 2e 91 08 bb 5e 54 |torisati|on....^T|
|00001950| 68 69 73 93 69 73 93 6e | 6f 74 93 61 93 74 72 75 |his.is.n|ot.a.tru|
|00001960| 65 93 43 68 6f 6c 65 73 | 6b 79 8e a1 8d 66 61 63 |e.Choles|ky...fac|
|00001970| 74 6f 72 69 73 61 74 69 | 6f 6e 2c 91 02 e8 a4 62 |torisati|on,....b|
|00001980| 75 74 96 02 cd 77 72 61 | 74 68 65 72 93 74 68 65 |ut...wra|ther.the|
|00001990| 93 43 68 6f 6c 65 73 6b | 79 93 66 61 63 74 6f 72 |.Cholesk|y.factor|
|000019a0| 69 73 61 74 69 6f 6e 93 | 6f 66 93 b5 41 96 01 29 |isation.|of..A..)|
|000019b0| 27 b2 2b 93 b5 44 96 03 | 14 94 b2 77 68 65 72 65 |'.+..D..|...where|
|000019c0| 9b 02 cd 77 b5 44 93 b2 | 69 73 98 61 98 64 69 61 |...w.D..|is.a.dia|
|000019d0| 67 6f 6e 61 6c 98 6d 61 | 74 72 69 78 98 77 69 74 |gonal.ma|trix.wit|
|000019e0| 68 98 6e 6f 6e 2d 6e 65 | 67 61 74 69 76 90 b8 e3 |h.non-ne|gativ...|
|000019f0| 65 8e a1 8d 65 6e 9a b8 | e3 74 72 69 65 73 2e 91 |e...en..|.tries..|
|00001a00| 06 53 0d 57 68 65 74 68 | 65 72 96 03 f5 c2 74 68 |.S.Wheth|er....th|
|00001a10| 65 93 b5 41 93 b2 6d 61 | 74 72 69 78 93 69 73 93 |e..A..ma|trix.is.|
|00001a20| 6d 6f 90 47 1d 64 69 0c | 65 64 93 69 6e 93 74 68 |mo.G.di.|ed.in.th|
|00001a30| 69 73 93 77 98 61 98 79 | 93 69 73 93 64 65 74 65 |is.w.a.y|.is.dete|
|00001a40| 72 6d 69 6e 65 64 93 62 | 98 79 93 74 68 65 93 da |rmined.b|.y.the..|
|00001a50| 74 6f 6c 93 b2 70 61 72 | 61 6d 65 74 65 72 3b 91 |tol..par|ameter;.|
|00001a60| 04 45 f8 74 68 65 93 64 | 69 61 67 6f 6e 61 6c 8e |.E.the.d|iagonal.|
|00001a70| a1 8d 65 6e 9a b8 e3 74 | 72 79 96 03 fc e9 6f 66 |..en...t|ry....of|
|00001a80| 93 74 68 65 93 43 68 6f | 6c 65 73 6b 79 93 66 61 |.the.Cho|lesky.fa|
|00001a90| 63 74 6f 72 69 73 61 74 | 69 6f 6e 93 69 73 93 65 |ctorisat|ion.is.e|
|00001aa0| 6e 73 75 72 65 64 93 74 | 6f 93 62 90 47 1d 65 91 |nsured.t|o.b.G.e.|
|00001ab0| 07 db 4c b8 15 8d 8d 91 | 03 de 63 9f f7 e4 fb 70 |..L.....|..c....p|
|00001ac0| 8e 8d 91 0c 33 ba 9f f7 | e4 fb 89 00 00 66 65 00 |....3...|.....fe.|
|00001ad0| 0f bf f7 9f 08 1b 05 8d | da 74 6f 6c 8e 8e 8e 91 |........|.tol....|
|00001ae0| 1b f3 b1 b2 2e 91 06 68 | 84 54 68 65 93 b5 44 91 |.......h|.The..D.|
|00001af0| 04 44 06 b2 6d 61 74 72 | 69 78 93 69 73 93 67 75 |.D..matr|ix.is.gu|
|00001b00| 61 72 61 6e 98 74 65 65 | 64 93 74 6f 93 62 90 47 |aran.tee|d.to.b.G|
|00001b10| 1d 65 93 7a 65 72 6f 93 | 69 6e 8e a1 8d 65 78 61 |.e.zero.|in...exa|
|00001b20| 63 74 96 03 55 55 61 72 | 69 74 68 6d 65 74 69 63 |ct..UUar|ithmetic|
|00001b30| 93 69 66 93 b5 75 8d 9f | fc 5e ff b4 54 8e 9b 06 |.if..u..|.^..T...|
|00001b40| 4c b6 b5 41 75 96 02 c7 | 18 b8 15 8d 93 da 74 6f |L..Au...|......to|
|00001b50| 6c 8e 91 12 87 0f b5 75 | 8d 9f fc 5e ff b4 54 8e |l......u|...^..T.|
|00001b60| 98 b5 75 96 03 55 55 b2 | 66 6f 72 93 61 6c 6c 93 |..u..UU.|for.all.|
|00001b70| b5 75 b2 2e 8e 9f 13 44 | 95 8d 8d db 45 58 41 4d |.u.....D|....EXAM|
|00001b80| 50 4c 45 8e 8e 9f 1a ca | 23 8d da 4d 41 54 91 14 |PLE.....|#..MAT..|
|00001b90| ff f4 2a 41 2c 96 05 3f | fd 2a 4c 4c 54 2c 93 2a |..*A,..?|.*LLT,.*|
|00001ba0| 4c 44 4c 3b 8e a1 8d 56 | 45 43 91 14 ff f4 2a 62 |LDL;...V|EC....*b|
|00001bb0| 2c 91 05 3f fd 2a 78 3b | 8e a1 8d 64 6f 75 62 6c |,..?.*x;|...doubl|
|00001bc0| 65 91 05 3f fd 74 6f 6c | 3b 8e a1 8d 91 0a 7f fa |e..?.tol|;.......|
|00001bd0| 2e 2e 2e 2e 2e 2e 8e a1 | 8d 41 96 05 3f fd 3d 93 |........|.A..?.=.|
|00001be0| 69 6e 5f 6d 61 74 28 4d | 4e 55 4c 4c 29 3b 8e a1 |in_mat(M|NULL);..|
|00001bf0| 8d 62 96 05 3f fd 3d 93 | 69 6e 5f 76 65 63 28 56 |.b..?.=.|in_vec(V|
|00001c00| 4e 55 4c 4c 29 3b 8e a1 | 8d 69 6e 70 75 74 28 22 |NULL);..|.input("|
|00001c10| 49 6e 70 75 74 96 05 3f | fd 74 6f 6c 93 66 6f 72 |Input..?|.tol.for|
|00001c20| 93 6d 6f 64 69 66 69 65 | 64 93 43 68 6f 6c 65 73 |.modifie|d.Choles|
|00001c30| 6b 79 3a 93 22 2c 93 22 | 25 6c 66 22 2c 93 26 74 |ky:.",."|%lf",.&t|
|00001c40| 6f 6c 29 3b 8e a1 8d 4c | 4c 54 96 05 3f fd 3d 93 |ol);...L|LT..?.=.|
|00001c50| 63 70 5f 6d 61 74 28 41 | 2c 4d 4e 55 4c 4c 29 3b |cp_mat(A|,MNULL);|
|00001c60| 8e 8e 8e 8c 8b 00 00 00 | 4a 00 00 00 00 00 00 00 |........|J.......|
|00001c70| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001c80| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 15 |........|........|
|00001c90| 4e a0 02 a9 52 b3 8d a0 | fd 62 ad 4d 8d 8d b2 37 |N...R...|.b.M...7|
|00001ca0| 34 92 00 99 14 c1 de 43 | 48 41 50 54 45 52 96 03 |4......C|HAPTER..|
|00001cb0| 55 55 35 2e 91 07 c7 1c | 44 45 4e 53 45 93 4d 41 |UU5.....|DENSE.MA|
|00001cc0| 9b ff 2a aa 54 52 49 58 | 93 46 91 fe e3 8d 41 90 |..*.TRIX|.F....A.|
|00001cd0| b8 e3 43 54 4f 52 49 53 | 41 98 54 49 4f 4e 93 4f |..CTORIS|A.TION.O|
|00001ce0| 50 45 52 41 98 54 49 4f | 4e 53 8e 8e a0 02 84 07 |PERA.TIO|NS......|
|00001cf0| 32 8d a0 fd 98 0a 15 8d | da 2f 2a 96 05 3f fd 49 |2.......|./*..?.I|
|00001d00| 66 93 41 93 70 6f 73 69 | 74 69 76 65 93 64 65 66 |f.A.posi|tive.def|
|00001d10| 69 6e 69 74 65 2e 2e 2e | 93 2a 2f 8e a4 0c 00 00 |inite...|.*/.....|
|00001d20| 8d 43 48 66 61 63 74 6f | 72 28 4c 4c 54 29 3b 8e |.CHfacto|r(LLT);.|
|00001d30| a1 8d 78 96 05 3f fd 3d | 93 43 48 73 6f 6c 76 65 |..x..?.=|.CHsolve|
|00001d40| 28 4c 4c 54 2c 62 2c 56 | 4e 55 4c 4c 29 3b 8e a1 |(LLT,b,V|NULL);..|
|00001d50| 8d 2f 2a 96 05 3f fd 2e | 2e 2e 6f 74 68 65 72 77 |./*..?..|..otherw|
|00001d60| 69 73 65 2c 93 67 65 74 | 93 61 70 70 72 6f 78 69 |ise,.get|.approxi|
|00001d70| 6d 61 74 65 93 73 6f 6c | 75 74 69 6f 6e 2e 2e 2e |mate.sol|ution...|
|00001d80| 93 2a 2f 8e a1 8d 4c 4c | 54 96 05 3f fd 3d 93 63 |.*/...LL|T..?.=.c|
|00001d90| 70 5f 6d 61 74 28 41 2c | 4d 4e 55 4c 4c 29 3b 8e |p_mat(A,|MNULL);.|
|00001da0| a1 8d 4d 43 48 66 61 63 | 74 6f 72 28 4c 4c 54 2c |..MCHfac|tor(LLT,|
|00001db0| 74 6f 6c 29 3b 91 1f 7f | ee 2f 2a 96 05 3f fd 4c |tol);...|./*..?.L|
|00001dc0| 4c 54 93 6e 6f 77 93 68 | 61 73 93 66 61 63 74 6f |LT.now.h|as.facto|
|00001dd0| 72 73 93 6f 66 93 41 93 | 2b 93 44 93 2a 2f 8e a1 |rs.of.A.|+.D.*/..|
|00001de0| 8d 4d 43 48 73 6f 6c 76 | 65 28 4c 4c 54 2c 62 2c |.MCHsolv|e(LLT,b,|
|00001df0| 78 29 3b 8e a1 8d 2f 2a | 96 05 3f fd 2e 2e 2e 6f |x);.../*|..?....o|
|00001e00| 72 93 75 73 65 93 4c 44 | 4c 93 66 61 63 74 6f 72 |r.use.LD|L.factor|
|00001e10| 69 73 61 74 69 6f 6e 93 | 2a 2f 8e a1 8d 4c 44 4c |isation.|*/...LDL|
|00001e20| 96 05 3f fd 3d 93 63 70 | 5f 6d 61 74 28 41 2c 4d |..?.=.cp|_mat(A,M|
|00001e30| 4e 55 4c 4c 29 3b 8e a1 | 8d 4c 44 4c 66 61 63 74 |NULL);..|.LDLfact|
|00001e40| 6f 72 28 4c 44 4c 29 3b | 8e a1 8d 4c 44 4c 73 6f |or(LDL);|...LDLso|
|00001e50| 6c 76 65 28 4c 44 4c 2c | 62 2c 78 29 3b 8e 9f 1d |lve(LDL,|b,x);...|
|00001e60| 44 95 8d 8d db 53 45 45 | 91 03 55 55 41 4c 53 4f |D....SEE|..UUALSO|
|00001e70| 8e 8e 9f 10 44 95 8d 91 | 0f 00 00 da 63 61 74 63 |....D...|....catc|
|00001e80| 68 28 29 96 03 55 55 b2 | 61 6e 64 93 da 42 4b 50 |h()..UU.|and..BKP|
|00001e90| 66 61 63 74 6f 72 28 29 | 8e 9f 13 44 95 8d 8d db |factor()|...D....|
|00001ea0| 53 4f 55 52 43 45 91 03 | 55 55 46 49 4c 45 3a 91 |SOURCE..|UUFILE:.|
|00001eb0| 14 00 02 da 43 48 66 61 | 63 74 6f 72 2e 63 8e 8e |....CHfa|ctor.c..|
|00001ec0| 8e 8e 8c 8b 00 00 00 4b | 00 00 00 00 00 00 00 00 |.......K|........|
|00001ed0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001ee0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 1e 98 |........|........|
|00001ef0| a0 02 a9 52 b3 8d a0 fd | 62 ad 4d 8d 8d 92 01 cb |...R....|b.M.....|
|00001f00| 78 6e de 37 35 8e 8e a0 | 02 84 07 32 8d a0 fd 98 |xn.75...|...2....|
|00001f10| 0a 15 8d 8d db 4e 41 4d | 45 8e 8e 9f 10 44 95 8d |.....NAM|E....D..|
|00001f20| 91 0f 00 00 da 4c 55 66 | 61 63 74 6f 72 2c 96 05 |.....LUf|actor,..|
|00001f30| 3f fd 4c 55 73 6f 6c 76 | 65 2c 93 4c 55 54 73 6f |?.LUsolv|e,.LUTso|
|00001f40| 6c 76 65 2c 93 4c 55 63 | 6f 6e 64 65 73 74 2c 93 |lve,.LUc|ondest,.|
|00001f50| 6d 8d 91 00 a1 45 89 00 | 00 66 66 00 03 26 66 8e |m....E..|.ff..&f.|
|00001f60| 91 03 c7 ab 69 6e 76 65 | 72 73 65 96 02 38 e3 b2 |....inve|rse..8..|
|00001f70| 7b 93 b5 4c 55 91 03 4f | fe b2 66 61 63 74 6f 72 |{..LU..O|..factor|
|00001f80| 69 73 61 74 69 6f 6e 93 | 28 47 61 75 73 73 69 61 |isation.|(Gaussia|
|00001f90| 6e 93 65 6c 69 6d 69 6e | 61 74 69 6f 6e 29 93 61 |n.elimin|ation).a|
|00001fa0| 6e 64 8e a4 0c 00 00 8d | 73 6f 6c 76 90 b8 e3 65 |nd......|solv...e|
|00001fb0| 8e 9f 13 44 95 8d 8d db | 53 59 4e 4f 50 53 49 53 |...D....|SYNOPSIS|
|00001fc0| 8e 8e 9f 16 a5 d4 8d da | 23 69 6e 63 6c 75 64 65 |........|#include|
|00001fd0| 91 05 3f fd 22 6d 61 74 | 72 69 78 2e 68 22 8e a1 |..?."mat|rix.h"..|
|00001fe0| 8d 23 69 6e 63 6c 75 64 | 65 91 05 3f fd 22 6d 61 |.#includ|e..?."ma|
|00001ff0| 74 72 69 78 32 2e 68 22 | 8e a1 8d 4d 41 54 91 1a |trix2.h"|...MAT..|
|00002000| 3f f1 2a 4c 55 66 61 63 | 74 6f 72 28 41 2c 70 69 |?.*LUfac|tor(A,pi|
|00002010| 76 6f 74 29 8e a1 8d 4d | 41 54 91 1a 3f f1 2a 41 |vot)...M|AT..?.*A|
|00002020| 3b 8e a1 8d 50 45 52 4d | 91 14 ff f4 2a 70 69 76 |;...PERM|....*piv|
|00002030| 6f 74 3b 8e a9 18 3e 23 | 8d 56 45 43 91 1a 3f f1 |ot;...>#|.VEC..?.|
|00002040| 2a 4c 55 73 6f 6c 76 65 | 28 41 2c 70 69 76 6f 74 |*LUsolve|(A,pivot|
|00002050| 2c 62 2c 78 29 8e a1 8d | 4d 41 54 91 1a 3f f1 2a |,b,x)...|MAT..?.*|
|00002060| 41 3b 8e a1 8d 50 45 52 | 4d 91 14 ff f4 2a 70 69 |A;...PER|M....*pi|
|00002070| 76 6f 74 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 62 |vot;...V|EC..?.*b|
|00002080| 2c 91 05 3f fd 2a 78 3b | 8e a6 8d 56 45 43 91 1a |,..?.*x;|...VEC..|
|00002090| 3f f1 2a 4c 55 54 73 6f | 6c 76 65 28 41 2c 70 69 |?.*LUTso|lve(A,pi|
|000020a0| 76 6f 74 2c 62 2c 78 29 | 8e a1 8d 4d 41 54 91 1a |vot,b,x)|...MAT..|
|000020b0| 3f f1 2a 41 3b 8e a1 8d | 50 45 52 4d 91 14 ff f4 |?.*A;...|PERM....|
|000020c0| 2a 70 69 76 6f 74 3b 8e | a1 8d 56 45 43 91 1a 3f |*pivot;.|..VEC..?|
|000020d0| f1 2a 62 2c 91 05 3f fd | 2a 78 3b 8e a6 8d 64 6f |.*b,..?.|*x;...do|
|000020e0| 75 62 6c 65 91 0a 7f fa | 4c 55 63 6f 6e 64 65 73 |uble....|LUcondes|
|000020f0| 74 28 4c 55 2c 70 69 76 | 6f 74 29 8e a1 8d 4d 41 |t(LU,piv|ot)...MA|
|00002100| 54 91 1a 3f f1 2a 4c 55 | 3b 8e a1 8d 50 45 52 4d |T..?.*LU|;...PERM|
|00002110| 91 14 ff f4 2a 70 69 76 | 6f 74 3b 8e a6 8d 4d 41 |....*piv|ot;...MA|
|00002120| 54 91 1a 3f f1 2a 6d 5f | 69 6e 76 65 72 73 65 28 |T..?.*m_|inverse(|
|00002130| 41 2c 6f 75 74 29 8e a1 | 8d 4d 41 54 91 1a 3f f1 |A,out)..|.MAT..?.|
|00002140| 2a 41 2c 91 05 3f fd 2a | 6f 75 74 3b 8e 9f 1b 3e |*A,..?.*|out;...>|
|00002150| 23 8d 8d db 44 45 53 43 | 52 49 50 54 49 4f 4e 8e |#...DESC|RIPTION.|
|00002160| 8e 9f 10 44 95 8d 91 0f | 00 00 b2 54 68 65 96 03 |...D....|...The..|
|00002170| b1 da 72 6f 75 74 69 6e | 65 93 da 4c 55 66 61 63 |..routin|e..LUfac|
|00002180| 74 6f 72 28 29 93 b2 70 | 90 47 1d 65 72 66 6f 72 |tor()..p|.G.erfor|
|00002190| 6d 73 93 b5 4c 55 91 04 | c8 f5 b2 66 61 63 74 6f |ms..LU..|...facto|
|000021a0| 72 69 73 61 74 69 6f 6e | 2c 91 03 c8 fb 77 68 69 |risation|,....whi|
|000021b0| 63 9a b8 e3 68 93 69 73 | 93 6f 74 68 65 72 77 69 |c...h.is|.otherwi|
|000021c0| 73 65 93 6b 6e 6f 98 77 | 6e 93 61 73 93 47 61 75 |se.kno.w|n.as.Gau|
|000021d0| 73 73 69 61 6e 93 65 6c | 69 6d 69 6e 61 74 69 6f |ssian.el|iminatio|
|000021e0| 6e 8e a1 8d 77 69 74 68 | 96 04 0b 1a 69 6d 70 6c |n...with|....impl|
|000021f0| 69 63 69 74 93 73 63 61 | 6c 65 64 93 70 61 72 74 |icit.sca|led.part|
|00002200| 69 61 6c 93 70 69 76 90 | b8 e3 6f 74 69 6e 67 2e |ial.piv.|..oting.|
|00002210| 91 06 93 17 54 68 65 93 | b5 4c 55 91 05 22 35 b2 |....The.|.LU.."5.|
|00002220| 66 61 63 74 6f 72 73 93 | 6f 66 93 da 41 9b 04 0a |factors.|of..A...|
|00002230| eb b2 61 72 65 93 73 74 | 6f 72 65 64 93 69 6e 93 |..are.st|ored.in.|
|00002240| da 41 98 b2 69 6e 93 63 | 6f 6d 70 61 63 74 93 66 |.A..in.c|ompact.f|
|00002250| 6f 72 6d 2e 91 06 93 17 | 4f 6e 63 65 93 74 68 69 |orm.....|Once.thi|
|00002260| 73 93 69 73 8e a1 8d 64 | 6f 6e 65 2c 91 03 85 b8 |s.is...d|one,....|
|00002270| 74 68 65 96 03 7c 0b 72 | 6f 75 74 69 6e 65 93 da |the..|.r|outine..|
|00002280| 4c 55 73 6f 6c 76 65 28 | 29 93 b2 63 61 6e 93 62 |LUsolve(|)..can.b|
|00002290| 90 47 1d 65 93 75 73 65 | 64 93 74 6f 93 73 6f 6c |.G.e.use|d.to.sol|
|000022a0| 76 9a b8 e3 65 93 65 71 | 75 61 74 69 6f 6e 73 93 |v...e.eq|uations.|
|000022b0| 6f 66 93 74 68 65 93 66 | 6f 72 6d 93 b5 41 78 96 |of.the.f|orm..Ax.|
|000022c0| 03 07 9c b2 3d 93 b5 62 | 96 03 7c 0b b2 66 6f 72 |....=..b|..|..for|
|000022d0| 93 b5 78 93 b2 62 98 79 | 93 66 6f 72 77 98 61 72 |..x..b.y|.forw.ar|
|000022e0| 64 93 61 6e 64 93 62 61 | 63 98 6b 8e a1 8d 73 75 |d.and.ba|c.k...su|
|000022f0| 62 73 74 69 74 75 74 69 | 6f 6e 2e 91 06 26 b7 54 |bstituti|on...&.T|
|00002300| 68 65 96 03 e6 fa 73 79 | 73 74 65 6d 93 b5 41 8d |he....sy|stem..A.|
|00002310| 9f fc 5e ff b4 54 8e 91 | 06 4c b6 b5 78 96 03 b9 |..^..T..|.L..x...|
|00002320| d5 b2 3d 93 b5 62 96 03 | e6 fa b2 63 61 6e 93 62 |..=..b..|...can.b|
|00002330| 90 47 1d 65 93 73 6f 6c | 76 9a b8 e3 65 64 93 62 |.G.e.sol|v...ed.b|
|00002340| 98 79 93 75 73 69 6e 67 | 93 da 4c 55 54 73 6f 6c |.y.using|..LUTsol|
|00002350| 76 65 28 29 b2 2e 91 06 | 26 b7 54 68 65 93 63 6f |ve()....|&.The.co|
|00002360| 90 47 1d 64 65 93 66 6f | 72 93 61 93 66 75 6c 6c |.G.de.fo|r.a.full|
|00002370| 93 66 61 63 74 6f 72 69 | 73 65 93 61 6e 64 8e a1 |.factori|se.and..|
|00002380| 8d 73 6f 6c 76 69 6e 67 | 9b 03 55 55 b5 41 78 96 |.solving|..UU.Ax.|
|00002390| 02 c7 18 b2 3d 93 b5 62 | 98 b2 61 6e 64 98 b5 41 |....=..b|..and..A|
|000023a0| 8d 9f fc 5e ff b4 54 8e | 91 06 4c b6 b5 79 91 03 |...^..T.|..L..y..|
|000023b0| 22 f1 b2 3d 93 b5 62 98 | b2 69 73 3a 8e a9 16 a5 |"..=..b.|.is:....|
|000023c0| d4 8d da 2f 2a 96 05 3f | fd 73 65 74 93 75 70 93 |.../*..?|.set.up.|
|000023d0| 41 93 61 6e 64 93 62 93 | 2a 2f 8e a1 8d 91 0a 7f |A.and.b.|*/......|
|000023e0| fa 2e 2e 2e 2e 2e 2e 8e | a1 8d 70 69 76 6f 74 96 |........|..pivot.|
|000023f0| 05 3f fd 3d 93 67 65 74 | 5f 70 65 72 6d 28 41 2d |.?.=.get|_perm(A-|
|00002400| 3e 6d 29 3b 8e a1 8d 78 | 96 05 3f fd 3d 93 67 65 |>m);...x|..?.=.ge|
|00002410| 74 5f 76 65 63 28 41 2d | 3e 6e 29 3b 8e a1 8d 79 |t_vec(A-|>n);...y|
|00002420| 96 05 3f fd 3d 93 67 65 | 74 5f 76 65 63 28 41 2d |..?.=.ge|t_vec(A-|
|00002430| 3e 6d 29 3b 8e a1 8d 4c | 55 96 05 3f fd 3d 93 63 |>m);...L|U..?.=.c|
|00002440| 70 5f 6d 61 74 28 41 2c | 4d 4e 55 4c 4c 29 3b 8e |p_mat(A,|MNULL);.|
|00002450| a1 8d 4c 55 66 61 63 74 | 6f 72 28 4c 55 2c 70 69 |..LUfact|or(LU,pi|
|00002460| 76 6f 74 29 3b 8e a1 8d | 78 96 05 3f fd 3d 93 4c |vot);...|x..?.=.L|
|00002470| 55 73 6f 6c 76 65 28 4c | 55 2c 70 69 76 6f 74 2c |Usolve(L|U,pivot,|
|00002480| 62 2c 78 29 3b 8e a1 8d | 79 96 05 3f fd 3d 93 4c |b,x);...|y..?.=.L|
|00002490| 55 54 73 6f 6c 76 65 28 | 4c 55 2c 70 69 76 6f 74 |UTsolve(|LU,pivot|
|000024a0| 2c 62 2c 79 29 3b 8e a1 | 8d 63 6f 6e 64 69 74 69 |,b,y);..|.conditi|
|000024b0| 6f 6e 96 05 3f fd 3d 93 | 4c 55 63 6f 6e 64 65 73 |on..?.=.|LUcondes|
|000024c0| 74 28 4c 55 2c 70 69 76 | 6f 74 29 3b 8e a6 8d 91 |t(LU,piv|ot);....|
|000024d0| 0f 00 00 b2 41 91 02 ce | c3 66 75 6c 6c 96 02 ce |....A...|.full...|
|000024e0| e5 64 65 73 63 72 69 70 | 74 69 6f 6e 93 6f 66 93 |.descrip|tion.of.|
|000024f0| 47 61 75 73 73 69 61 6e | 93 65 6c 69 6d 69 6e 61 |Gaussian|.elimina|
|00002500| 74 69 6f 6e 93 77 69 74 | 68 93 70 61 72 74 69 61 |tion.wit|h.partia|
|00002510| 6c 93 70 69 76 9a b8 e3 | 6f 74 69 6e 67 93 61 6e |l.piv...|oting.an|
|00002520| 64 93 69 74 73 93 6e 98 | 75 6d 65 72 69 63 61 6c |d.its.n.|umerical|
|00002530| 93 62 90 47 1d 65 68 61 | 98 76 69 6f 75 72 93 63 |.b.G.eha|.viour.c|
|00002540| 61 6e 93 62 90 47 1d 65 | 93 66 6f 75 6e 64 8e a1 |an.b.G.e|.found..|
|00002550| 8d 69 6e 96 03 ad 60 61 | 93 6e 95 b8 e3 75 6d 93 |.in...`a|.n...um.|
|00002560| 62 9a 47 1d 65 72 96 03 | ad 60 6f 66 93 62 98 6f |b.G.er..|.`of.b.o|
|00002570| 98 6f 6b 73 2c 91 03 c3 | 63 74 68 6f 75 67 68 93 |.oks,...|cthough.|
|00002580| 77 90 b8 e3 65 93 72 65 | 66 65 72 93 74 68 65 93 |w...e.re|fer.the.|
|00002590| 72 65 61 64 65 72 93 73 | 70 98 65 63 69 0c 63 61 |reader.s|p.eci.ca|
|000025a0| 6c 6c 79 93 74 6f 93 dc | 4d 61 74 72 69 78 91 03 |lly.to..|Matrix..|
|000025b0| e4 e7 43 6f 6d 70 75 74 | 61 74 69 6f 6e 73 91 04 |..Comput|ations..|
|000025c0| 7f 81 b2 62 90 b8 e3 79 | 93 47 2e 48 2e 93 47 6f |...b...y|.G.H..Go|
|000025d0| 6c 75 62 93 61 6e 64 8e | a1 8d 43 2e 96 03 8f 4e |lub.and.|..C....N|
|000025e0| 76 91 ff 71 c7 61 6e 93 | 4c 6f 61 6e 2c 9b 03 9d |v..q.an.|Loan,...|
|000025f0| cd 4e 6f 72 74 68 93 4f | 78 66 6f 72 64 93 41 63 |.North.O|xford.Ac|
|00002600| 61 64 65 6d 69 63 2c 8d | 98 b8 78 8e 8d 91 08 0f |ademic,.|..x.....|
|00002610| 96 78 8e 91 0c 81 5f b2 | 33 2e 32 7b 33 2e 34 2c |.x...._.|3.2{3.4,|
|00002620| 98 70 70 2e 93 39 32 7b | 31 32 32 2c 98 32 6e 64 |.pp..92{|122,.2nd|
|00002630| 93 45 64 69 74 69 6f 6e | 93 28 31 39 38 39 29 2e |.Edition|.(1989).|
|00002640| 91 05 1f b3 54 68 65 93 | 76 91 ff 71 c7 61 72 69 |....The.|v..q.ari|
|00002650| 61 6e 90 b8 e3 74 93 68 | 65 72 65 93 69 73 93 74 |an...t.h|ere.is.t|
|00002660| 68 61 74 8e a1 8d 73 63 | 61 6c 69 6e 67 96 04 62 |hat...sc|aling..b|
|00002670| 70 69 73 93 75 73 65 64 | 93 dc 69 6d 70 6c 69 63 |pis.used|..implic|
|00002680| 69 74 6c 79 b2 2e 91 07 | 99 19 54 68 61 74 93 69 |itly....|..That.i|
|00002690| 73 2c 91 04 a5 b7 73 63 | 61 6c 69 6e 67 93 69 73 |s,....sc|aling.is|
|000026a0| 93 6f 6e 6c 79 93 75 73 | 65 64 93 74 6f 93 64 65 |.only.us|ed.to.de|
|000026b0| 63 69 64 65 93 77 68 69 | 63 9a b8 e3 68 93 72 6f |cide.whi|c...h.ro|
|000026c0| 98 77 73 93 74 6f 93 73 | 77 98 61 70 93 64 75 72 |.ws.to.s|w.ap.dur|
|000026d0| 69 6e 67 93 74 68 65 93 | 70 61 72 74 69 61 6c 8e |ing.the.|partial.|
|000026e0| a1 8d 70 69 76 90 b8 e3 | 6f 74 69 6e 67 91 03 55 |..piv...|oting..U|
|000026f0| 55 70 72 6f 90 47 1d 63 | 65 73 73 2e 8e 8e 8e 8c |Upro.G.c|ess.....|
|00002700| 8b 00 00 00 4c 00 00 00 | 00 00 00 00 00 00 00 00 |....L...|........|
|00002710| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002720| 00 00 00 00 00 00 00 00 | 00 00 00 20 f7 a0 02 a9 |........|... ....|
|00002730| 52 b3 8d a0 fd 62 ad 4d | 8d 8d b2 37 36 92 00 99 |R....b.M|...76...|
|00002740| 14 c1 de 43 48 41 50 54 | 45 52 96 03 55 55 35 2e |...CHAPT|ER..UU5.|
|00002750| 91 07 c7 1c 44 45 4e 53 | 45 93 4d 41 9b ff 2a aa |....DENS|E.MA..*.|
|00002760| 54 52 49 58 93 46 91 fe | e3 8d 41 90 b8 e3 43 54 |TRIX.F..|..A...CT|
|00002770| 4f 52 49 53 41 98 54 49 | 4f 4e 93 4f 50 45 52 41 |ORISA.TI|ON.OPERA|
|00002780| 98 54 49 4f 4e 53 8e 8e | a0 02 84 07 32 8d a0 fd |.TIONS..|....2...|
|00002790| 98 0a 15 8d 91 0f 00 00 | b2 4e 6f 74 65 96 03 5e |........|.Note..^|
|000027a0| 73 74 68 61 74 93 74 68 | 65 93 66 61 63 74 6f 72 |sthat.th|e.factor|
|000027b0| 69 73 61 74 69 6f 6e 93 | 72 6f 75 74 69 6e 65 93 |isation.|routine.|
|000027c0| da 4c 55 66 61 63 74 6f | 72 28 29 93 b2 6d 61 9a |.LUfacto|r()..ma.|
|000027d0| b8 e3 79 93 73 75 63 63 | 65 65 64 93 77 68 65 72 |..y.succ|eed.wher|
|000027e0| 65 93 74 68 65 93 73 6f | 6c 76 98 65 93 72 6f 75 |e.the.so|lv.e.rou|
|000027f0| 74 69 6e 65 93 da 4c 55 | 73 6f 6c 76 65 28 29 93 |tine..LU|solve().|
|00002800| b2 66 61 69 6c 73 93 69 | 66 2c 8e a4 0c 00 00 8d |.fails.i|f,......|
|00002810| 66 6f 72 96 03 4d 4c 65 | 78 61 6d 70 6c 65 2c 91 |for..MLe|xample,.|
|00002820| 03 4e e7 da 41 9b 03 4d | 4a b2 69 73 93 73 69 6e |.N..A..M|J.is.sin|
|00002830| 67 75 6c 61 72 2e 91 04 | 6f 19 41 6c 73 6f 93 6e |gular...|o.Also.n|
|00002840| 6f 74 65 93 74 68 61 74 | 93 b5 4c 55 91 04 64 67 |ote.that|..LU..dg|
|00002850| b2 66 61 63 74 6f 72 69 | 73 61 74 69 6f 6e 93 61 |.factori|sation.a|
|00002860| 6c 73 6f 93 73 75 63 63 | 65 65 64 73 93 77 68 65 |lso.succ|eeds.whe|
|00002870| 6e 93 da 41 98 b2 69 73 | 93 6e 6f 74 93 65 76 90 |n..A..is|.not.ev.|
|00002880| b8 e3 65 6e 93 73 71 75 | 61 72 65 2c 91 03 4e e7 |..en.squ|are,..N.|
|00002890| 74 68 6f 75 67 68 8e a1 | 8d 74 68 69 73 96 03 30 |though..|.this..0|
|000028a0| dc 69 73 93 61 93 72 65 | 71 75 69 72 65 6d 65 6e |.is.a.re|quiremen|
|000028b0| 9a b8 e3 74 93 66 6f 72 | 93 74 68 65 93 73 75 63 |...t.for|.the.suc|
|000028c0| 63 65 73 73 93 6f 66 93 | da 4c 55 73 6f 6c 76 65 |cess.of.|.LUsolve|
|000028d0| 28 29 b2 2e 91 04 65 9f | 45 72 72 6f 72 73 93 61 |()....e.|Errors.a|
|000028e0| 72 65 93 72 61 69 73 65 | 64 93 62 98 79 93 da 4c |re.raise|d.b.y..L|
|000028f0| 55 66 61 63 74 6f 72 28 | 29 93 b2 69 66 93 da 41 |Ufactor(|)..if..A|
|00002900| 91 03 30 d3 b2 6f 72 93 | da 70 69 76 6f 74 93 b2 |..0..or.|.pivot..|
|00002910| 69 73 93 4e 55 4c 4c 2c | 8e a1 8d 69 66 96 03 d7 |is.NULL,|...if...|
|00002920| 58 74 68 65 93 73 69 7a | 65 93 6f 66 93 da 70 69 |Xthe.siz|e.of..pi|
|00002930| 76 6f 74 93 b2 69 73 93 | 6c 65 73 73 93 74 68 61 |vot..is.|less.tha|
|00002940| 6e 93 74 68 65 93 6e 95 | b8 e3 75 6d 93 62 90 47 |n.the.n.|..um.b.G|
|00002950| 1d 65 72 96 03 d7 58 6f | 66 93 72 6f 9a b8 e3 77 |.er...Xo|f.ro...w|
|00002960| 73 93 6f 66 93 da 41 b2 | 2e 93 45 72 72 6f 72 73 |s.of..A.|..Errors|
|00002970| 93 61 72 65 93 72 61 69 | 73 65 64 93 62 98 79 93 |.are.rai|sed.b.y.|
|00002980| da 4c 55 73 6f 6c 76 65 | 93 b2 69 66 93 74 68 65 |.LUsolve|..if.the|
|00002990| 73 65 93 63 6f 6e 64 69 | 74 69 6f 6e 73 8e a1 8d |se.condi|tions...|
|000029a0| 6f 90 47 1d 63 63 75 72 | 2c 91 03 3d 09 69 66 96 |o.G.ccur|,..=.if.|
|000029b0| 03 36 f6 da 62 93 b2 69 | 73 93 4e 55 4c 4c 2c 93 |.6..b..i|s.NULL,.|
|000029c0| 69 66 93 da 41 9b 03 36 | ee b2 69 73 93 6e 6f 74 |if..A..6|..is.not|
|000029d0| 93 73 71 75 61 72 65 2e | 91 04 67 a7 54 68 65 6e |.square.|..g.Then|
|000029e0| 93 69 66 93 da 78 93 b2 | 69 73 93 4e 55 4c 4c 98 |.if..x..|is.NULL.|
|000029f0| 6f 72 93 74 6f 90 47 1d | 6f 93 73 6d 61 6c 6c 93 |or.to.G.|o.small.|
|00002a00| 74 6f 93 63 6f 6e 9a b8 | e3 74 61 69 6e 93 74 68 |to.con..|.tain.th|
|00002a10| 65 93 72 65 73 75 6c 74 | 93 61 93 6e 65 77 93 76 |e.result|.a.new.v|
|00002a20| 98 65 63 74 6f 72 93 6f | 66 8e a1 8d 74 68 65 96 |.ector.o|f...the.|
|00002a30| 02 dd 03 61 70 70 72 6f | 70 72 69 61 74 65 93 73 |...appro|priate.s|
|00002a40| 69 7a 65 93 69 73 93 63 | 72 65 61 74 65 64 2e 9b |ize.is.c|reated..|
|00002a50| 04 49 ac 49 6e 93 65 69 | 74 68 65 72 93 63 61 73 |.I.In.ei|ther.cas|
|00002a60| 65 93 74 68 65 93 73 6f | 6c 75 74 69 6f 6e 93 6f |e.the.so|lution.o|
|00002a70| 66 93 b5 41 78 96 02 c7 | 18 b2 3d 93 b5 62 b2 2c |f..Ax...|..=..b.,|
|00002a80| 96 02 f5 14 da 78 b2 2c | 93 69 73 96 02 dd 03 72 |.....x.,|.is....r|
|00002a90| 65 74 75 72 6e 65 64 2e | 98 54 68 65 93 72 6f 75 |eturned.|.The.rou|
|00002aa0| 74 69 6e 65 73 93 da 4c | 55 73 6f 6c 76 65 28 29 |tines..L|Usolve()|
|00002ab0| 8e a1 8d b2 61 6e 64 96 | 03 98 fe da 4c 55 54 73 |....and.|....LUTs|
|00002ac0| 6f 6c 76 65 28 29 93 b2 | 6d 61 90 b8 e3 79 93 dc |olve()..|ma...y..|
|00002ad0| 6e 6f 74 91 04 8b d6 b2 | 62 9a 47 1d 65 93 75 73 |not.....|b.G.e.us|
|00002ae0| 65 64 93 dc 69 6e 91 03 | d2 27 73 69 74 75 91 00 |ed..in..|.'situ..|
|00002af0| c4 63 b2 3b 91 03 ba d3 | 74 68 61 74 93 69 73 2c |.c.;....|that.is,|
|00002b00| 91 03 a9 e8 77 69 74 68 | 93 da 62 96 05 3f fd 3d |....with|..b..?.=|
|00002b10| 3d 93 78 b2 2e 91 05 3c | c3 54 68 69 73 96 03 98 |=.x....<|.This...|
|00002b20| fe 69 73 93 62 98 65 63 | 61 75 73 65 93 6f 66 93 |.is.b.ec|ause.of.|
|00002b30| 74 68 65 93 70 98 65 72 | 6d 90 b8 e3 75 74 61 74 |the.p.er|m...utat|
|00002b40| 69 6f 6e 93 6f 66 93 74 | 68 65 8e a1 8d 76 90 b8 |ion.of.t|he...v..|
|00002b50| e3 65 63 74 6f 72 2e 8e | a9 10 44 95 8d 91 0f 00 |.ector..|..D.....|
|00002b60| 00 54 68 65 96 03 aa 07 | 63 6f 6e 64 69 74 69 6f |.The....|conditio|
|00002b70| 6e 93 6e 95 b8 e3 75 6d | 93 62 90 47 1d 65 72 9b |n.n...um|.b.G.er.|
|00002b80| 03 aa 07 28 72 65 6c 61 | 74 69 76 93 65 98 74 6f |...(rela|tiv.e.to|
|00002b90| 98 74 68 65 98 69 6e 0c | 6e 69 74 93 79 98 6e 6f |.the.in.|nit.y.no|
|00002ba0| 72 6d 29 98 63 61 6e 98 | 62 90 47 1d 65 98 dc 65 |rm).can.|b.G.e..e|
|00002bb0| 73 74 69 6d 61 74 65 91 | ff 7d 27 64 91 04 b2 90 |stimate.|.}'d....|
|00002bc0| b2 75 73 69 6e 67 98 74 | 68 65 98 72 6f 75 74 69 |.using.t|he.routi|
|00002bd0| 6e 65 98 da 4c 55 63 6f | 6e 64 65 73 74 28 29 b2 |ne..LUco|ndest().|
|00002be0| 2e 8e a1 8d 54 68 69 73 | 96 02 8c d4 65 73 74 69 |....This|....esti|
|00002bf0| 6d 61 74 65 93 69 73 93 | 6e 6f 74 93 67 75 61 72 |mate.is.|not.guar|
|00002c00| 61 6e 9a b8 e3 74 65 65 | 64 93 74 6f 93 75 6e 64 |an...tee|d.to.und|
|00002c10| 65 72 2d 93 6f 72 93 6f | 98 76 98 65 72 2d 65 73 |er-.or.o|.v.er-es|
|00002c20| 74 69 6d 61 74 65 93 74 | 68 65 93 74 72 75 65 93 |timate.t|he.true.|
|00002c30| 63 6f 6e 64 69 74 69 6f | 6e 93 6e 98 75 6d 98 62 |conditio|n.n.um.b|
|00002c40| 90 47 1d 65 72 3b 91 02 | cf aa 68 6f 98 77 98 65 |.G.er;..|..ho.w.e|
|00002c50| 76 98 65 72 2c 91 02 b4 | ee 69 74 93 63 61 6e 93 |v.er,...|.it.can.|
|00002c60| 75 73 75 61 6c 6c 79 8e | a1 8d 62 90 47 1d 65 96 |usually.|..b.G.e.|
|00002c70| 02 d9 c1 72 65 6c 69 65 | 64 93 6f 6e 93 74 6f 93 |...relie|d.on.to.|
|00002c80| 67 69 76 9a b8 e3 65 93 | 61 6e 93 65 73 74 69 6d |giv...e.|an.estim|
|00002c90| 61 74 65 93 63 6f 72 72 | 65 63 74 93 74 6f 93 77 |ate.corr|ect.to.w|
|00002ca0| 69 74 68 69 6e 93 61 6e | 93 6f 72 64 65 72 93 6f |ithin.an|.order.o|
|00002cb0| 66 93 6d 61 67 6e 69 74 | 75 64 65 2c 91 02 f2 78 |f.magnit|ude,...x|
|00002cc0| 77 68 69 63 98 68 93 69 | 73 93 75 73 75 61 6c 6c |whic.h.i|s.usuall|
|00002cd0| 79 93 61 6c 6c 93 74 68 | 61 74 93 69 73 93 72 65 |y.all.th|at.is.re|
|00002ce0| 71 75 69 72 65 64 2e 8e | a6 8d 91 0f 00 00 54 68 |quired..|......Th|
|00002cf0| 65 96 03 5f 44 72 6f 75 | 74 69 6e 65 93 da 6d 5f |e.._Drou|tine..m_|
|00002d00| 69 6e 76 65 72 73 65 28 | 29 93 b2 63 6f 6d 70 75 |inverse(|)..compu|
|00002d10| 74 65 73 93 74 68 65 93 | 69 6e 95 b8 e3 76 93 65 |tes.the.|in...v.e|
|00002d20| 72 73 65 96 03 5f 44 6f | 66 93 da 41 91 03 5f 42 |rse.._Do|f..A.._B|
|00002d30| b2 61 6e 64 93 72 65 74 | 75 72 6e 73 93 69 74 93 |.and.ret|urns.it.|
|00002d40| 69 6e 93 da 6f 75 74 b2 | 2e 91 04 8f 95 54 68 69 |in..out.|.....Thi|
|00002d50| 73 93 69 73 93 63 61 72 | 72 69 65 64 93 6f 75 74 |s.is.car|ried.out|
|00002d60| 93 75 73 69 6e 67 93 74 | 68 65 8e a1 8d 4c 55 91 |.using.t|he...LU.|
|00002d70| 03 44 aa 66 61 63 74 6f | 72 69 73 61 74 69 6f 6e |.D.facto|risation|
|00002d80| 96 03 44 af 72 6f 75 74 | 69 6e 65 73 2e 91 04 6c |..D.rout|ines...l|
|00002d90| 3a 41 73 93 69 73 93 75 | 73 75 61 6c 6c 79 93 6e |:As.is.u|sually.n|
|00002da0| 6f 74 65 64 93 69 6e 93 | 6e 9a b8 e3 75 6d 65 72 |oted.in.|n...umer|
|00002db0| 69 63 61 6c 93 61 6e 61 | 6c 79 73 69 73 93 74 65 |ical.ana|lysis.te|
|00002dc0| 78 74 73 2c 91 03 48 03 | 69 6e 98 76 98 65 72 73 |xts,..H.|in.v.ers|
|00002dd0| 65 93 6d 61 74 72 69 63 | 65 73 93 73 68 6f 75 6c |e.matric|es.shoul|
|00002de0| 64 93 72 61 72 65 6c 79 | 93 62 90 47 1d 65 8e a1 |d.rarely|.b.G.e..|
|00002df0| 8d 63 6f 6d 70 75 74 65 | 64 2e 91 04 3b 69 49 66 |.compute|d...;iIf|
|00002e00| 96 02 b2 39 61 93 73 79 | 73 74 65 6d 93 6f 66 93 |...9a.sy|stem.of.|
|00002e10| 65 71 75 61 74 69 6f 6e | 73 93 6e 65 65 64 93 74 |equation|s.need.t|
|00002e20| 6f 93 62 9a 47 1d 65 93 | 73 6f 6c 76 90 b8 e3 65 |o.b.G.e.|solv...e|
|00002e30| 64 2c 91 02 d2 d8 75 73 | 65 93 74 68 65 93 61 62 |d,....us|e.the.ab|
|00002e40| 98 6f 95 b8 e3 76 93 65 | 96 02 b2 39 63 6f 98 64 |.o...v.e|...9co.d|
|00002e50| 65 93 63 61 6c 6c 69 6e | 67 93 da 4c 55 66 61 63 |e.callin|g..LUfac|
|00002e60| 74 6f 72 28 29 93 b2 61 | 6e 64 93 da 4c 55 73 6f |tor()..a|nd..LUso|
|00002e70| 6c 76 65 28 29 8e a1 8d | b2 64 69 72 65 63 74 6c |lve()...|.directl|
|00002e80| 79 91 ff 2a aa 2e 8e 9f | 13 44 95 8d 8d db 53 4f |y..*....|.D....SO|
|00002e90| 55 52 43 45 91 03 55 55 | 46 49 4c 45 3a 91 14 00 |URCE..UU|FILE:...|
|00002ea0| 02 da 4c 55 66 61 63 74 | 6f 72 2e 63 8e 8e 8e 8e |..LUfact|or.c....|
|00002eb0| 8c 8b 00 00 00 4d 00 00 | 00 00 00 00 00 00 00 00 |.....M..|........|
|00002ec0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002ed0| 00 00 00 00 00 00 00 00 | 00 00 00 00 29 34 a0 02 |........|....)4..|
|00002ee0| a9 52 b3 8d a0 fd 62 ad | 4d 8d 8d 92 01 cb 78 6e |.R....b.|M.....xn|
|00002ef0| de 37 37 8e 8e a0 02 84 | 07 32 8d a0 fd 98 0a 15 |.77.....|.2......|
|00002f00| 8d 8d db 4e 41 4d 45 8e | 8e 9f 10 44 95 8d 91 0f |...NAME.|...D....|
|00002f10| 00 00 da 51 52 66 61 63 | 74 6f 72 2c 96 05 3f fd |...QRfac|tor,..?.|
|00002f20| 51 52 43 50 66 61 63 74 | 6f 72 2c 93 51 52 73 6f |QRCPfact|or,.QRso|
|00002f30| 6c 76 65 2c 93 51 52 54 | 73 6f 6c 76 65 2c 93 51 |lve,.QRT|solve,.Q|
|00002f40| 52 63 6f 6e 64 65 73 74 | 96 03 55 55 b2 7b 93 b5 |Rcondest|..UU.{..|
|00002f50| 51 52 91 03 69 1c b2 66 | 61 63 74 6f 72 69 73 61 |QR..i..f|actorisa|
|00002f60| 74 69 6f 6e 93 61 6e 64 | 93 73 6f 6c 76 90 b8 e3 |tion.and|.solv...|
|00002f70| 65 8e 9f 13 44 95 8d 8d | db 53 59 4e 4f 50 53 49 |e...D...|.SYNOPSI|
|00002f80| 53 8e 8e 9f 17 fd af 8d | da 23 69 6e 63 6c 75 64 |S.......|.#includ|
|00002f90| 65 91 05 3f fd 22 6d 61 | 74 72 69 78 2e 68 22 8e |e..?."ma|trix.h".|
|00002fa0| a4 0c 00 00 8d 23 69 6e | 63 6c 75 64 65 91 05 3f |.....#in|clude..?|
|00002fb0| fd 22 6d 61 74 72 69 78 | 32 2e 68 22 8e a1 8d 4d |."matrix|2.h"...M|
|00002fc0| 41 54 91 1a 3f f1 2a 51 | 52 66 61 63 74 6f 72 28 |AT..?.*Q|Rfactor(|
|00002fd0| 41 2c 64 69 61 67 29 8e | a1 8d 4d 41 54 91 1a 3f |A,diag).|..MAT..?|
|00002fe0| f1 2a 41 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 64 |.*A;...V|EC..?.*d|
|00002ff0| 69 61 67 3b 8e a9 19 eb | f5 8d 4d 41 54 91 1a 3f |iag;....|..MAT..?|
|00003000| f1 2a 51 52 43 50 66 61 | 63 74 6f 72 28 41 2c 64 |.*QRCPfa|ctor(A,d|
|00003010| 69 61 67 2c 70 69 76 6f | 74 29 8e a1 8d 4d 41 54 |iag,pivo|t)...MAT|
|00003020| 91 1a 3f f1 2a 41 3b 8e | a1 8d 56 45 43 91 1a 3f |..?.*A;.|..VEC..?|
|00003030| f1 2a 64 69 61 67 3b 8e | a1 8d 50 45 52 4d 91 14 |.*diag;.|..PERM..|
|00003040| ff f4 2a 70 69 76 6f 74 | 3b 8e a6 8d 56 45 43 91 |..*pivot|;...VEC.|
|00003050| 1a 3f f1 2a 51 52 73 6f | 6c 76 65 28 41 2c 64 69 |.?.*QRso|lve(A,di|
|00003060| 61 67 2c 62 2c 78 29 8e | a1 8d 4d 41 54 91 1a 3f |ag,b,x).|..MAT..?|
|00003070| f1 2a 41 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 64 |.*A;...V|EC..?.*d|
|00003080| 69 61 67 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 62 |iag;...V|EC..?.*b|
|00003090| 2c 91 05 3f fd 2a 78 3b | 8e a6 8d 56 45 43 91 1a |,..?.*x;|...VEC..|
|000030a0| 3f f1 2a 51 52 54 73 6f | 6c 76 65 28 41 2c 64 69 |?.*QRTso|lve(A,di|
|000030b0| 61 67 2c 62 2c 78 29 8e | a1 8d 4d 41 54 91 1a 3f |ag,b,x).|..MAT..?|
|000030c0| f1 2a 41 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 64 |.*A;...V|EC..?.*d|
|000030d0| 69 61 67 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 62 |iag;...V|EC..?.*b|
|000030e0| 2c 91 05 3f fd 2a 78 3b | 8e a6 8d 56 45 43 91 1a |,..?.*x;|...VEC..|
|000030f0| 3f f1 2a 51 52 43 50 73 | 6f 6c 76 65 28 41 2c 64 |?.*QRCPs|olve(A,d|
|00003100| 69 61 67 2c 70 69 76 6f | 74 2c 62 2c 78 29 8e a1 |iag,pivo|t,b,x)..|
|00003110| 8d 4d 41 54 91 1a 3f f1 | 2a 41 3b 8e a1 8d 56 45 |.MAT..?.|*A;...VE|
|00003120| 43 91 1a 3f f1 2a 64 69 | 61 67 3b 8e a1 8d 56 45 |C..?.*di|ag;...VE|
|00003130| 43 91 1a 3f f1 2a 62 2c | 91 05 3f fd 2a 78 3b 8e |C..?.*b,|..?.*x;.|
|00003140| a6 8d 64 6f 75 62 6c 65 | 91 0a 7f fa 51 52 63 6f |..double|....QRco|
|00003150| 6e 64 65 73 74 28 51 52 | 29 8e a1 8d 4d 41 54 91 |ndest(QR|)...MAT.|
|00003160| 1a 3f f1 2a 51 52 3b 8e | 9f 1c eb f5 8d 8d db 44 |.?.*QR;.|.......D|
|00003170| 45 53 43 52 49 50 54 49 | 4f 4e 8e 8e a9 10 44 95 |ESCRIPTI|ON....D.|
|00003180| 8d 91 0f 00 00 b2 54 68 | 65 96 02 f9 80 72 6f 75 |......Th|e....rou|
|00003190| 74 69 6e 65 93 da 51 52 | 66 61 63 74 6f 72 28 29 |tine..QR|factor()|
|000031a0| 93 b2 70 90 47 1d 65 72 | 66 6f 72 6d 73 93 61 93 |..p.G.er|forms.a.|
|000031b0| 73 74 72 61 69 67 68 95 | b8 e3 74 66 6f 72 77 93 |straigh.|..tforw.|
|000031c0| 61 72 64 96 02 f9 80 b5 | 51 52 91 03 0d 47 b2 66 |ard.....|QR...G.f|
|000031d0| 61 63 74 6f 72 69 73 61 | 74 69 6f 6e 93 6f 66 93 |actorisa|tion.of.|
|000031e0| da 41 b2 2e 93 46 91 ff | 2a aa 6f 72 93 74 68 6f |.A...F..|*.or.tho|
|000031f0| 73 65 93 75 6e 66 61 6d | 69 6c 69 61 72 93 77 69 |se.unfam|iliar.wi|
|00003200| 74 68 93 74 68 65 8e a1 | 8d 74 65 72 6d 69 6e 6f |th.the..|.termino|
|00003210| 6c 6f 67 79 91 ff 2a aa | 2c 96 03 55 55 74 68 65 |logy..*.|,..UUthe|
|00003220| 93 b5 51 52 91 03 69 1c | b2 66 61 63 74 6f 72 69 |..QR..i.|.factori|
|00003230| 73 61 74 69 6f 6e 93 6f | 66 93 b5 41 93 b2 69 73 |sation.o|f..A..is|
|00003240| 93 61 93 66 61 63 74 6f | 72 69 73 61 74 69 6f 6e |.a.facto|risation|
|00003250| 93 6f 66 93 74 68 65 93 | 66 6f 72 6d 8e a4 15 a7 |.of.the.|form....|
|00003260| 60 8d 92 00 d7 b2 8c b5 | 41 96 02 c7 18 b2 3d 93 |`.......|A.....=.|
|00003270| b5 51 91 01 aa a8 52 8e | a1 8d b2 77 68 65 72 65 |.Q....R.|...where|
|00003280| 96 03 9a 60 b5 51 93 b2 | 69 73 93 6f 72 74 68 6f |...`.Q..|is.ortho|
|00003290| 67 6f 6e 61 6c 93 61 6e | 64 93 b5 52 91 03 ae 27 |gonal.an|d..R...'|
|000032a0| b2 69 73 93 75 70 70 90 | 47 1d 65 72 93 74 72 69 |.is.upp.|G.er.tri|
|000032b0| 61 6e 67 75 6c 61 72 2e | 91 05 40 e8 54 68 69 73 |angular.|..@.This|
|000032c0| 93 66 61 63 74 6f 72 69 | 73 61 74 69 6f 6e 93 65 |.factori|sation.e|
|000032d0| 78 69 73 74 73 93 77 68 | 65 74 68 65 72 93 6f 72 |xists.wh|ether.or|
|000032e0| 93 6e 6f 74 93 b5 41 93 | b2 69 73 93 73 69 6e 67 |.not..A.|.is.sing|
|000032f0| 75 6c 61 72 93 6f 72 8e | a4 0c 00 00 8d 65 76 90 |ular.or.|.....ev.|
|00003300| b8 e3 65 6e 96 03 64 f8 | 73 71 75 61 72 65 2e 9b |..en..d.|square..|
|00003310| 04 a0 b0 54 68 65 93 b5 | 51 52 91 03 78 bf b2 66 |...The..|QR..x..f|
|00003320| 61 63 74 6f 72 69 73 61 | 74 69 6f 6e 93 69 73 93 |actorisa|tion.is.|
|00003330| 70 90 47 1d 65 72 66 6f | 72 6d 65 64 93 75 73 69 |p.G.erfo|rmed.usi|
|00003340| 6e 67 93 48 6f 75 73 65 | 68 6f 6c 64 65 72 93 74 |ng.House|holder.t|
|00003350| 72 61 6e 73 66 6f 72 6d | 61 74 69 6f 6e 73 2e 98 |ransform|ations..|
|00003360| 28 54 68 65 73 65 93 61 | 72 65 93 6f 72 74 68 6f |(These.a|re.ortho|
|00003370| 67 6f 6e 61 6c 8e a1 8d | 6d 61 74 72 69 63 65 73 |gonal...|matrices|
|00003380| 96 03 55 55 6f 66 93 74 | 68 65 93 66 6f 72 6d 93 |..UUof.t|he.form.|
|00003390| b5 50 8d 9f 01 7f ff b4 | 69 8e 9b 06 1b 64 b2 3d |.P......|i....d.=|
|000033a0| 91 02 c7 18 b5 49 91 03 | 01 c2 b8 00 91 02 38 e0 |.....I..|......8.|
|000033b0| b5 0b 8d 9f 01 7f ff b4 | 69 8e 96 03 54 4c b5 76 |........|i...TL.v|
|000033c0| 8d 9f 01 7f ff b4 69 8e | 93 b5 76 8d 9f fc 5e ff |......i.|..v...^.|
|000033d0| 8d 90 5b d9 b4 54 8e 9f | 06 3b 5a 8d 69 8e 8e 91 |..[..T..|.;Z.i...|
|000033e0| 09 fd e4 b2 77 68 65 72 | 65 91 03 55 55 b5 0b 8d |....wher|e..UU...|
|000033f0| 9f 01 7f ff b4 69 8e 98 | b2 3d 91 02 c7 18 32 b5 |.....i..|.=....2.|
|00003400| 3d 76 8d 9f fc 5e ff 8d | 90 5b d9 b4 54 8e 9f 06 |=v...^..|.[..T...|
|00003410| 3b 5a 8d 69 8e 8e 91 06 | a8 8f b5 76 8d 9f 01 7f |;Z.i....|...v....|
|00003420| ff b4 69 8e 93 b2 2e 29 | 8e a6 8d 91 0f 00 00 54 |..i....)|.......T|
|00003430| 68 65 96 03 c5 9c 72 6f | 75 74 69 6e 65 93 da 51 |he....ro|utine..Q|
|00003440| 52 43 50 66 61 63 74 6f | 72 28 29 93 b2 70 90 47 |RCPfacto|r()..p.G|
|00003450| 1d 65 72 66 6f 72 6d 73 | 93 61 93 b5 51 52 91 03 |.erforms|.a..QR..|
|00003460| d9 63 b2 66 61 63 74 6f | 72 69 73 61 74 69 6f 6e |.c.facto|risation|
|00003470| 93 77 69 74 68 93 63 6f | 6c 75 6d 6e 93 70 69 76 |.with.co|lumn.piv|
|00003480| 95 b8 e3 6f 74 69 6e 67 | 2c 91 03 e1 ae 77 68 69 |...oting|,....whi|
|00003490| 63 93 68 96 03 c5 9c 69 | 73 93 61 93 66 61 63 74 |c.h....i|s.a.fact|
|000034a0| 6f 72 69 73 61 74 69 6f | 6e 8e a1 8d 6f 66 96 03 |orisatio|n...of..|
|000034b0| 55 55 74 68 65 93 66 6f | 72 6d 8e a1 8d 92 00 cf |UUthe.fo|rm......|
|000034c0| f6 dd b5 41 96 01 aa a8 | b2 05 8d 9f fb de ff b4 |...A....|........|
|000034d0| 54 8e 91 09 13 ce b2 3d | 91 02 c7 18 b5 51 93 52 |T......=|.....Q.R|
|000034e0| 8e 9f 11 ca d3 8d b2 77 | 68 65 72 65 96 03 35 1c |.......w|here..5.|
|000034f0| 61 64 64 69 74 69 6f 6e | 61 6c 6c 79 91 ff 2a aa |addition|ally..*.|
|00003500| 2c 91 03 3b 8e 05 93 69 | 73 93 61 93 70 90 47 1d |,..;...i|s.a.p.G.|
|00003510| 65 72 6d 9a b8 e3 75 74 | 61 74 69 6f 6e 93 6d 61 |erm...ut|ation.ma|
|00003520| 74 72 69 78 2e 91 04 67 | 09 54 68 65 93 05 93 6d |trix...g|.The...m|
|00003530| 61 74 72 69 78 93 69 73 | 93 72 65 70 72 65 73 65 |atrix.is|.represe|
|00003540| 6e 98 74 65 64 93 62 98 | 79 93 da 70 69 76 6f 74 |n.ted.b.|y..pivot|
|00003550| b2 2e 91 04 67 09 54 68 | 69 73 93 69 73 93 64 6f |....g.Th|is.is.do|
|00003560| 6e 65 93 65 78 61 63 74 | 6c 79 8e a1 8d 61 73 96 |ne.exact|ly...as.|
|00003570| 03 55 55 66 6f 72 93 da | 51 52 66 61 63 74 6f 72 |.UUfor..|QRfactor|
|00003580| 28 29 93 b2 65 78 63 65 | 70 74 93 66 6f 72 93 74 |()..exce|pt.for.t|
|00003590| 68 65 93 70 69 76 90 b8 | e3 6f 74 69 6e 67 2e 8e |he.piv..|.oting..|
|000035a0| a6 8d 91 0f 00 00 42 6f | 74 68 96 03 f9 f1 6f 66 |......Bo|th....of|
|000035b0| 93 74 68 65 73 65 93 66 | 61 63 74 6f 72 69 73 61 |.these.f|actorisa|
|000035c0| 74 69 6f 6e 73 93 61 72 | 65 93 70 90 47 1d 65 72 |tions.ar|e.p.G.er|
|000035d0| 66 6f 72 6d 65 64 93 dc | 69 6e 91 04 2b 58 73 69 |formed..|in..+Xsi|
|000035e0| 74 75 91 00 c4 63 b2 2c | 91 04 23 18 61 6e 64 93 |tu...c.,|..#.and.|
|000035f0| 73 74 6f 72 65 93 74 68 | 65 93 b5 51 93 b2 61 6e |store.th|e..Q..an|
|00003600| 64 93 b5 52 91 04 0d b8 | b2 66 61 63 74 6f 72 73 |d..R....|.factors|
|00003610| 93 63 6f 6d 70 61 63 74 | 6c 79 93 69 6e 93 da 41 |.compact|ly.in..A|
|00003620| 91 03 f9 c6 b2 61 6e 64 | 8e a1 8d da 64 69 61 67 |.....and|....diag|
|00003630| b2 2e 91 07 bf 6a 54 68 | 69 73 96 04 6f 36 63 6f |.....jTh|is..o6co|
|00003640| 6d 70 61 63 74 93 66 6f | 72 6d 93 69 73 93 75 73 |mpact.fo|rm.is.us|
|00003650| 65 64 93 63 6f 6e 73 69 | 73 74 65 6e 9a b8 e3 74 |ed.consi|sten...t|
|00003660| 6c 79 93 77 69 74 68 69 | 6e 93 74 68 69 73 93 70 |ly.withi|n.this.p|
|00003670| 61 63 98 6b 91 ff 71 c7 | 61 67 65 2c 91 04 b5 ae |ac.k..q.|age,....|
|00003680| 61 6e 64 93 69 73 93 65 | 73 73 65 6e 98 74 69 61 |and.is.e|ssen.tia|
|00003690| 6c 6c 79 93 74 68 61 74 | 93 6f 66 93 47 6f 6c 75 |lly.that|.of.Golu|
|000036a0| 62 93 61 6e 64 8e 8e 8e | 8c 8b 00 00 00 4e 00 00 |b.and...|.....N..|
|000036b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000036c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000036d0| 00 00 00 00 30 e5 a0 02 | a9 52 b3 8d a0 fd 62 ad |....0...|.R....b.|
|000036e0| 4d 8d 8d b2 37 38 92 00 | 99 14 c1 de 43 48 41 50 |M...78..|....CHAP|
|000036f0| 54 45 52 96 03 55 55 35 | 2e 91 07 c7 1c 44 45 4e |TER..UU5|.....DEN|
|00003700| 53 45 93 4d 41 9b ff 2a | aa 54 52 49 58 93 46 91 |SE.MA..*|.TRIX.F.|
|00003710| fe e3 8d 41 90 b8 e3 43 | 54 4f 52 49 53 41 98 54 |...A...C|TORISA.T|
|00003720| 49 4f 4e 93 4f 50 45 52 | 41 98 54 49 4f 4e 53 8e |ION.OPER|A.TIONS.|
|00003730| 8e a0 02 84 07 32 8d a0 | fd 98 0a 15 8d b2 76 91 |.....2..|......v.|
|00003740| ff 71 c7 61 6e 96 03 30 | 26 4c 6f 61 6e 27 73 93 |.q.an..0|&Loan's.|
|00003750| dc 4d 61 74 72 69 78 91 | 03 71 b2 43 6f 6d 70 75 |.Matrix.|.q.Compu|
|00003760| 74 61 74 69 6f 6e 73 b2 | 2c 8d 91 03 37 96 b8 78 |tations.|,...7..x|
|00003770| 8e 91 07 a9 5f b2 35 2e | 32 2c 9b 03 37 96 70 2e |...._.5.|2,..7.p.|
|00003780| 93 32 31 32 2c 98 32 6e | 64 93 65 64 69 74 69 6f |.212,.2n|d.editio|
|00003790| 6e 2c 98 28 31 39 38 39 | 29 93 65 78 63 65 70 74 |n,.(1989|).except|
|000037a0| 93 74 68 61 74 93 74 68 | 65 93 b5 76 90 5b d9 b2 |.that.th|e..v.[..|
|000037b0| 27 73 93 61 72 65 93 6e | 6f 74 93 6e 6f 72 6d 61 |'s.are.n|ot.norma|
|000037c0| 6c 69 73 65 64 93 69 6e | 8e a4 0c 00 00 8d 74 68 |lised.in|......th|
|000037d0| 69 73 96 03 27 b1 70 61 | 63 9a b8 e3 6b 91 ff 71 |is..'.pa|c...k..q|
|000037e0| c7 61 67 65 2e 91 04 62 | 90 54 68 65 93 64 69 6d |.age...b|.The.dim|
|000037f0| 65 6e 73 69 6f 6e 73 93 | 6f 66 93 62 90 47 1d 6f |ensions.|of.b.G.o|
|00003800| 74 68 93 da 64 69 61 67 | 93 b2 6d 98 75 73 74 93 |th..diag|..m.ust.|
|00003810| 62 90 47 1d 65 93 61 74 | 93 6c 65 61 73 74 93 61 |b.G.e.at|.least.a|
|00003820| 73 93 6c 61 72 67 65 93 | 61 73 93 74 68 65 93 6d |s.large.|as.the.m|
|00003830| 69 6e 69 6d 98 75 6d 93 | 6f 66 93 74 68 65 93 6e |inim.um.|of.the.n|
|00003840| 98 75 6d 98 62 90 47 1d | 65 72 93 6f 66 93 72 6f |.um.b.G.|er.of.ro|
|00003850| 98 77 73 8e a1 8d 61 6e | 64 96 03 55 55 63 6f 6c |.ws...an|d..UUcol|
|00003860| 75 6d 6e 73 93 6f 66 93 | da 41 b2 2e 8e 9f 10 44 |umns.of.|.A.....D|
|00003870| 95 8d 91 0f 00 00 4f 6e | 63 65 96 04 62 9d da 41 |......On|ce..b..A|
|00003880| b2 2c 93 da 64 69 61 67 | 93 b2 63 6f 6e 9a b8 e3 |.,..diag|..con...|
|00003890| 74 61 69 6e 93 74 68 69 | 73 93 63 6f 6d 70 61 63 |tain.thi|s.compac|
|000038a0| 74 93 72 65 70 72 65 73 | 65 6e 98 74 61 74 69 6f |t.repres|en.tatio|
|000038b0| 6e 93 6f 66 93 74 68 65 | 93 b5 51 52 91 04 76 64 |n.of.the|..QR..vd|
|000038c0| b2 66 61 63 74 6f 72 73 | 93 6f 66 93 b5 41 b2 2c |.factors|.of..A.,|
|000038d0| 91 04 a5 ef 77 98 65 93 | 63 61 6e 93 75 73 65 93 |....w.e.|can.use.|
|000038e0| da 51 52 73 6f 6c 76 65 | 28 29 93 b2 74 6f 8e a1 |.QRsolve|()..to..|
|000038f0| 8d 73 6f 6c 76 90 b8 e3 | 65 96 03 38 10 73 79 73 |.solv...|e..8.sys|
|00003900| 74 65 6d 73 93 6f 66 93 | 6c 69 6e 65 61 72 93 65 |tems.of.|linear.e|
|00003910| 71 75 61 74 69 6f 6e 73 | 2c 9b 03 3d eb 61 6e 64 |quations|,..=.and|
|00003920| 93 69 6e 64 65 65 64 2c | 98 0c 6e 64 93 6c 65 61 |.indeed,|..nd.lea|
|00003930| 73 74 93 73 71 75 61 72 | 65 93 65 72 72 6f 72 93 |st.squar|e.error.|
|00003940| 73 6f 6c 75 74 69 6f 6e | 73 93 74 6f 93 6f 95 b8 |solution|s.to.o..|
|00003950| e3 76 93 65 72 64 65 74 | 65 72 6d 69 6e 65 64 96 |.v.erdet|ermined.|
|00003960| 03 38 10 73 79 73 74 65 | 6d 73 93 6f 66 8e a1 8d |.8.syste|ms.of...|
|00003970| 65 71 75 61 74 69 6f 6e | 73 2e 9b 04 71 c7 53 65 |equation|s...q.Se|
|00003980| 65 96 03 55 55 dc 4d 61 | 74 72 69 78 91 03 93 e7 |e..UU.Ma|trix....|
|00003990| 43 6f 6d 70 75 74 61 74 | 69 6f 6e 73 b2 2c 8d 93 |Computat|ions.,..|
|000039a0| b8 78 8e 91 07 c7 1e b2 | 31 2e 34 2c 93 70 2e 93 |.x......|1.4,.p..|
|000039b0| 31 31 93 66 6f 72 93 61 | 6e 93 65 78 61 6d 70 6c |11.for.a|n.exampl|
|000039c0| 65 2e 98 49 6e 64 65 65 | 64 2c 93 74 68 65 93 63 |e..Indee|d,.the.c|
|000039d0| 6f 90 47 1d 64 65 8e a9 | 15 2d 43 8d da 4d 41 54 |o.G.de..|.-C..MAT|
|000039e0| 91 14 ff f4 2a 51 52 3b | 8e a1 8d 91 0a 7f fa 2e |....*QR;|........|
|000039f0| 2e 2e 2e 2e 2e 8e a1 8d | 51 52 96 05 3f fd 3d 93 |........|QR..?.=.|
|00003a00| 63 70 5f 6d 61 74 28 41 | 2c 4d 4e 55 4c 4c 29 3b |cp_mat(A|,MNULL);|
|00003a10| 8e a1 8d 51 52 66 61 63 | 74 6f 72 28 51 52 2c 64 |...QRfac|tor(QR,d|
|00003a20| 69 61 67 29 3b 8e a1 8d | 51 52 73 6f 6c 76 65 28 |iag);...|QRsolve(|
|00003a30| 51 52 2c 64 69 61 67 2c | 62 2c 78 29 3b 8e a6 8d |QR,diag,|b,x);...|
|00003a40| b2 0c 6e 64 73 96 03 55 | 55 74 68 65 93 6c 65 61 |..nds..U|Uthe.lea|
|00003a50| 73 74 93 73 71 75 61 72 | 65 73 93 73 6f 6c 75 74 |st.squar|es.solut|
|00003a60| 69 6f 6e 93 b5 78 93 b2 | 74 6f 8e a1 8d 92 00 d9 |ion..x..|to......|
|00003a70| 17 cb b5 41 91 01 aa a8 | 78 96 02 c7 18 b8 19 93 |...A....|x.......|
|00003a80| b5 62 3a 8e 9f 0f ae 83 | 8d b2 53 69 6d 69 6c 61 |.b:.....|..Simila|
|00003a90| 72 6c 79 91 ff 2a aa 2c | 9b 02 d3 f1 69 66 96 02 |rly..*.,|....if..|
|00003aa0| b3 98 da 51 52 43 50 66 | 61 63 74 6f 72 28 29 93 |...QRCPf|actor().|
|00003ab0| b2 69 73 93 74 6f 93 62 | 90 47 1d 65 93 75 73 65 |.is.to.b|.G.e.use|
|00003ac0| 64 93 74 6f 93 66 61 63 | 74 6f 72 93 b5 41 b2 2c |d.to.fac|tor..A.,|
|00003ad0| 98 74 68 65 6e 93 da 51 | 52 43 50 73 6f 6c 76 65 |.then..Q|RCPsolve|
|00003ae0| 28 29 93 b2 63 61 6e 93 | 62 90 47 1d 65 93 75 73 |()..can.|b.G.e.us|
|00003af0| 65 64 93 74 6f 93 73 6f | 6c 76 90 b8 e3 65 93 74 |ed.to.so|lv...e.t|
|00003b00| 68 65 93 6c 65 61 73 74 | 93 73 71 75 61 72 65 73 |he.least|.squares|
|00003b10| 8e a1 8d 70 72 6f 62 6c | 65 6d 9b 03 55 55 b5 41 |...probl|em..UU.A|
|00003b20| 78 96 02 c7 18 b8 19 93 | b5 62 b2 2e 91 04 71 c7 |x.......|.b....q.|
|00003b30| 54 68 65 98 63 6f 90 47 | 1d 64 65 98 74 6f 98 64 |The.co.G|.de.to.d|
|00003b40| 6f 98 74 68 69 73 98 69 | 73 3a 8e a6 8d da 51 52 |o.this.i|s:....QR|
|00003b50| 96 05 3f fd 3d 93 63 70 | 5f 6d 61 74 28 41 2c 4d |..?.=.cp|_mat(A,M|
|00003b60| 4e 55 4c 4c 29 3b 8e a1 | 8d 51 52 43 50 66 61 63 |NULL);..|.QRCPfac|
|00003b70| 74 6f 72 28 51 52 2c 64 | 69 61 67 2c 70 69 76 6f |tor(QR,d|iag,pivo|
|00003b80| 74 29 3b 8e a1 8d 51 52 | 43 50 73 6f 6c 76 65 28 |t);...QR|CPsolve(|
|00003b90| 51 52 2c 64 69 61 67 2c | 70 69 76 6f 74 2c 62 2c |QR,diag,|pivot,b,|
|00003ba0| 78 29 3b 8e a6 8d 91 0f | 00 00 b2 4e 6f 74 65 96 |x);.....|...Note.|
|00003bb0| 04 12 d3 74 68 61 74 93 | da 51 52 54 73 6f 6c 76 |...that.|.QRTsolv|
|00003bc0| 65 28 51 52 2c 64 69 61 | 67 2c 62 2c 78 29 93 b2 |e(QR,dia|g,b,x)..|
|00003bd0| 73 6f 6c 76 90 b8 e3 65 | 73 93 74 68 65 93 dc 75 |solv...e|s.the..u|
|00003be0| 6e 64 65 72 96 ff 7d 27 | 64 65 74 65 72 6d 69 6e |nder..}'|determin|
|00003bf0| 65 93 64 91 05 1b 5c b2 | 70 72 6f 62 6c 65 6d 9b |e.d...\.|problem.|
|00003c00| 04 12 d3 b5 41 78 96 04 | 02 e9 b2 3d 93 b5 62 b2 |....Ax..|...=..b.|
|00003c10| 3b 91 04 71 92 74 68 61 | 74 98 69 73 2c 91 04 42 |;..q.tha|t.is,..B|
|00003c20| 33 69 74 98 63 6f 6d 70 | 75 74 65 73 8e a1 8d 74 |3it.comp|utes...t|
|00003c30| 68 65 96 03 55 55 6d 69 | 6e 69 6d 90 b8 e3 75 6d |he..UUmi|nim...um|
|00003c40| 93 32 2d 6e 6f 72 6d 93 | b5 78 93 b2 74 68 61 74 |.2-norm.|.x..that|
|00003c50| 93 73 61 74 69 73 0c 65 | 73 93 b5 41 78 96 02 c7 |.satis.e|s..Ax...|
|00003c60| 18 b2 3d 93 b5 62 96 03 | 55 55 b2 66 6f 72 93 b5 |..=..b..|UU.for..|
|00003c70| 6d 96 02 c7 18 b8 14 93 | b5 6e b2 2e 8e 9f 10 44 |m.......|.n.....D|
|00003c80| 95 8d 91 0f 00 00 54 68 | 65 96 03 4d 2e 63 6f 6e |......Th|e..M.con|
|00003c90| 64 69 74 69 6f 6e 93 6e | 95 b8 e3 75 6d 93 62 9a |dition.n|...um.b.|
|00003ca0| 47 1d 65 72 96 03 4d 2e | 6f 66 93 61 93 6d 61 74 |G.er..M.|of.a.mat|
|00003cb0| 72 69 78 93 66 61 63 74 | 6f 72 65 64 93 75 73 69 |rix.fact|ored.usi|
|00003cc0| 6e 67 93 65 69 74 68 65 | 72 93 da 51 52 66 61 63 |ng.eithe|r..QRfac|
|00003cd0| 74 6f 72 28 29 93 b2 6f | 72 93 da 51 52 43 50 66 |tor()..o|r..QRCPf|
|00003ce0| 61 63 74 6f 72 28 29 93 | b2 63 61 6e 93 62 98 65 |actor().|.can.b.e|
|00003cf0| 93 65 73 74 69 6d 61 74 | 65 64 8e a1 8d 75 73 69 |.estimat|ed...usi|
|00003d00| 6e 67 91 03 55 55 da 51 | 52 63 6f 6e 64 65 73 74 |ng..UU.Q|Rcondest|
|00003d10| 28 29 b2 3a 8e a6 8d da | 70 72 69 6e 74 66 28 22 |().:....|printf("|
|00003d20| 32 2d 6e 6f 72 6d 96 05 | 3f fd 63 6f 6e 64 69 74 |2-norm..|?.condit|
|00003d30| 69 6f 6e 93 6e 75 6d 62 | 65 72 93 61 70 70 72 6f |ion.numb|er.appro|
|00003d40| 78 2e 93 3d 93 25 67 5c | 6e 22 2c 93 51 52 63 6f |x..=.%g\|n",.QRco|
|00003d50| 6e 64 65 73 74 28 51 52 | 29 29 3b 8e a6 8d b2 54 |ndest(QR|));....T|
|00003d60| 68 65 96 04 42 f2 66 75 | 6e 63 74 69 6f 6e 93 da |he..B.fu|nction..|
|00003d70| 51 52 63 6f 6e 64 65 73 | 74 28 29 93 b2 72 65 74 |QRcondes|t()..ret|
|00003d80| 75 72 6e 73 93 61 93 dc | 6c 6f 77 65 72 91 04 6e |urns.a..|lower..n|
|00003d90| 82 62 91 ff 7d 27 6f 75 | 6e 64 91 05 4b 7b b2 66 |.b..}'ou|nd..K{.f|
|00003da0| 6f 72 93 74 68 65 93 6c | 65 61 73 74 93 73 71 75 |or.the.l|east.squ|
|00003db0| 61 72 65 73 93 63 6f 6e | 64 69 74 69 6f 6e 93 6e |ares.con|dition.n|
|00003dc0| 95 b8 e3 75 6d 93 62 90 | 47 1d 65 72 96 04 42 f2 |...um.b.|G.er..B.|
|00003dd0| 6f 66 93 74 68 65 93 66 | 61 63 74 6f 72 65 64 8e |of.the.f|actored.|
|00003de0| a1 8d 6d 61 74 72 69 78 | 91 03 55 55 b5 41 8e a1 |..matrix|..UU.A..|
|00003df0| 8d 92 00 ba 9a 8b 14 8d | 9f 01 7f ff b4 4c 53 8e |........|.....LS.|
|00003e00| 91 0b 47 7b b2 28 b5 41 | b2 29 96 02 c7 18 3d 93 |..G{.(.A|.)....=.|
|00003e10| b8 6b b5 41 b8 6b 8d 9f | 01 7f ff f3 06 d9 93 a0 |.k.A.k..|........|
|00003e20| 52 00 07 00 00 00 07 00 | 00 00 04 63 6d 72 37 b1 |R.......|...cmr7.|
|00003e30| 32 8e 91 04 7c 73 b8 6b | b5 41 8d 9f fb de ff b1 |2...|s.k|.A......|
|00003e40| 2b 8e 91 06 a3 90 b8 6b | 8d 9f 01 7f ff b1 32 8e |+......k|......2.|
|00003e50| 8e 9f 0f ae 83 8d b2 70 | 72 6f 9a b8 e3 76 69 64 |.......p|ro...vid|
|00003e60| 65 64 96 03 55 55 b5 41 | 93 b2 68 61 73 93 66 75 |ed..UU.A|..has.fu|
|00003e70| 6c 6c 93 72 61 6e 6b 2e | 91 04 71 c7 49 66 93 b5 |ll.rank.|..q.If..|
|00003e80| 41 93 b2 69 73 93 73 71 | 75 61 72 65 2c 93 74 68 |A..is.sq|uare,.th|
|00003e90| 65 6e 93 74 68 69 73 93 | 69 73 93 65 78 61 63 74 |en.this.|is.exact|
|00003ea0| 6c 79 93 65 71 75 61 6c | 93 74 6f 93 74 68 65 93 |ly.equal|.to.the.|
|00003eb0| 32 2d 6e 6f 72 6d 93 63 | 6f 6e 64 69 74 69 6f 6e |2-norm.c|ondition|
|00003ec0| 93 6e 98 75 6d 98 62 90 | 47 1d 65 72 8e a4 12 22 |.n.um.b.|G.er..."|
|00003ed0| da 8d 92 00 ba 90 0e b5 | 14 8d 9f 01 7f ff b1 32 |........|.......2|
|00003ee0| 8e 9b 04 7c 73 b2 28 b5 | 41 b2 29 96 02 c7 18 3d |...|s.(.|A.)....=|
|00003ef0| 93 b8 6b b5 41 b8 6b 8d | 9f 01 7f ff b1 32 8e 98 |..k.A.k.|.....2..|
|00003f00| b8 6b b5 41 8d 9f fb de | ff f3 0c 4f 21 e2 85 00 |.k.A....|...O!...|
|00003f10| 07 00 00 00 07 00 00 00 | 05 63 6d 73 79 37 b7 00 |........|.cmsy7..|
|00003f20| b1 31 8e 91 0a bc 74 b8 | 6b 8d 9f 01 7f ff b1 32 |.1....t.|k......2|
|00003f30| 8e 98 b5 3a 8e a1 8d b2 | 49 66 96 03 55 55 74 68 |...:....|If..UUth|
|00003f40| 65 93 b5 51 52 91 03 69 | 1c b2 66 61 63 74 6f 72 |e..QR..i|..factor|
|00003f50| 73 93 61 72 65 93 65 78 | 61 63 74 6c 79 93 73 69 |s.are.ex|actly.si|
|00003f60| 6e 67 75 6c 61 72 2c 93 | 74 68 65 6e 93 da 51 52 |ngular,.|then..QR|
|00003f70| 63 6f 6e 64 65 73 74 28 | 29 93 b2 77 69 6c 6c 93 |condest(|)..will.|
|00003f80| 72 65 74 75 72 6e 93 da | 48 55 47 45 93 b2 28 da |return..|HUGE..(.|
|00003f90| 48 55 47 45 5f 56 41 4c | 93 b2 66 6f 72 93 41 4e |HUGE_VAL|..for.AN|
|00003fa0| 53 49 93 43 29 2e 8e a9 | 10 44 95 8d 91 0f 00 00 |SI.C)...|.D......|
|00003fb0| 54 68 65 96 03 2e 96 65 | 73 74 69 6d 61 74 65 93 |The....e|stimate.|
|00003fc0| 69 73 93 6f 62 74 61 69 | 6e 65 64 93 62 90 b8 e3 |is.obtai|ned.b...|
|00003fd0| 79 93 6f 62 74 61 69 6e | 69 6e 67 93 65 73 74 69 |y.obtain|ing.esti|
|00003fe0| 6d 61 74 65 73 93 66 6f | 72 93 b8 6b b5 52 9a 13 |mates.fo|r..k.R..|
|00003ff0| c7 b8 6b 8d 9f 01 7f ff | b1 32 8e 91 07 ab 09 b2 |..k.....|.2......|
|00004000| 61 6e 64 93 b8 6b b5 52 | 8d 98 9f fc 5e ff b7 00 |and..k.R|....^...|
|00004010| b1 31 8e 91 0a d0 3b b8 | 6b 8d 9f 01 7f ff b1 32 |.1....;.|k......2|
|00004020| 8e 91 04 7c 73 b2 2e 91 | 04 64 dd 4e 6f 74 65 93 |...|s...|.d.Note.|
|00004030| 74 68 61 74 93 b5 51 93 | b2 61 6e 64 93 05 93 64 |that..Q.|.and...d|
|00004040| 6f 93 6e 6f 74 93 61 0b | 65 63 74 8e a4 0c 00 00 |o.not.a.|ect.....|
|00004050| 8d 74 68 65 96 03 fd f0 | 32 2d 6e 6f 72 6d 93 6f |.the....|2-norm.o|
|00004060| 72 93 6c 65 61 73 74 93 | 73 71 75 61 72 65 73 93 |r.least.|squares.|
|00004070| 63 6f 6e 64 69 74 69 6f | 6e 93 6e 95 b8 e3 75 6d |conditio|n.n...um|
|00004080| 93 62 90 47 1d 65 72 73 | 2e 91 06 6b 98 54 68 65 |.b.G.ers|...k.The|
|00004090| 96 03 fd f0 65 73 74 69 | 6d 61 74 65 93 6f 66 93 |....esti|mate.of.|
|000040a0| b8 6b b5 52 8d 90 13 c7 | 9f fc 5e ff b7 00 b1 31 |.k.R....|..^....1|
|000040b0| 8e 91 0a d0 3b b8 6b 8d | 9f 01 7f ff b1 32 8e 91 |....;.k.|.....2..|
|000040c0| 08 7a 63 b2 69 73 93 66 | 6f 75 6e 64 93 75 73 69 |.zc.is.f|ound.usi|
|000040d0| 6e 67 93 74 68 65 93 74 | 65 63 90 b8 e3 68 6e 69 |ng.the.t|ec...hni|
|000040e0| 71 75 65 73 93 6f 66 8e | a1 8d da 4c 55 63 6f 6e |ques.of.|...LUcon|
|000040f0| 64 65 73 74 28 29 96 03 | a4 1b b2 74 6f 93 6f 62 |dest()..|...to.ob|
|00004100| 74 61 69 6e 93 61 93 76 | 9a b8 e3 65 63 74 6f 72 |tain.a.v|...ector|
|00004110| 93 b5 79 91 03 ff f4 b2 | 77 69 74 68 93 75 6e 69 |..y.....|with.uni|
|00004120| 74 93 b8 31 b2 2d 6e 6f | 72 6d 93 73 75 63 98 68 |t..1.-no|rm.suc.h|
|00004130| 93 74 68 61 74 93 b8 6b | b5 52 8d 90 13 c7 9f fc |.that..k|.R......|
|00004140| 5e ff b7 00 b1 31 8e 91 | 0a d0 3b b5 79 90 5b d9 |^....1..|..;.y.[.|
|00004150| b8 6b 8d 9f 01 7f ff b7 | 31 8e 91 0c 1d 00 b2 69 |.k......|1......i|
|00004160| 73 93 71 75 69 74 65 93 | 73 6d 61 6c 6c 2e 91 05 |s.quite.|small...|
|00004170| 5e 19 54 68 69 73 93 69 | 73 93 64 65 73 63 72 69 |^.This.i|s.descri|
|00004180| 62 90 47 1d 65 64 8e a1 | 8d 69 6e 96 03 5d cc 47 |b.G.ed..|.in..].G|
|00004190| 6f 6c 75 62 93 61 6e 64 | 93 76 91 ff 71 c7 61 6e |olub.and|.v..q.an|
|000041a0| 93 4c 6f 61 6e 2c 9b 03 | 5f ea 32 6e 64 93 45 64 |.Loan,..|_.2nd.Ed|
|000041b0| 69 74 69 6f 6e 93 70 70 | 2e 93 31 32 38 7b 31 33 |ition.pp|..128{13|
|000041c0| 30 2c 98 28 31 39 38 39 | 29 2e 91 04 8b 2b 54 68 |0,.(1989|)....+Th|
|000041d0| 65 6e 93 74 68 65 93 70 | 9a 47 1d 6f 95 b8 e3 77 |en.the.p|.G.o...w|
|000041e0| 93 65 72 96 03 5d cc 6d | 65 74 68 6f 98 64 93 69 |.er..].m|etho.d.i|
|000041f0| 73 93 61 70 70 6c 69 65 | 64 93 74 6f 93 74 68 65 |s.applie|d.to.the|
|00004200| 93 6d 61 74 72 69 78 8e | a1 8d 28 b5 52 8d 90 13 |.matrix.|..(.R...|
|00004210| c7 9f fc 5e ff b4 54 8e | 91 06 60 7d b5 52 9a 13 |...^..T.|..`}.R..|
|00004220| c7 b2 29 8d 9f fc 5e ff | b7 00 b1 31 8e 91 0e b8 |..)...^.|...1....|
|00004230| d5 b2 61 96 03 fc 61 74 | 6f 74 61 6c 93 6f 66 93 |..a...at|otal.of.|
|00004240| 74 68 72 65 65 93 74 69 | 6d 65 73 93 77 69 74 68 |three.ti|mes.with|
|00004250| 93 69 6e 69 74 69 61 6c | 93 76 90 b8 e3 65 63 74 |.initial|.v...ect|
|00004260| 6f 72 93 b5 79 90 5b d9 | b2 2e 91 06 66 eb 54 68 |or..y.[.|....f.Th|
|00004270| 65 93 63 6f 72 72 65 73 | 70 90 47 1d 6f 6e 64 69 |e.corres|p.G.ondi|
|00004280| 6e 67 93 65 73 74 69 6d | 61 74 65 93 6f 66 93 b8 |ng.estim|ate.of..|
|00004290| 6b b5 52 98 b8 6b 8d 9f | 01 7f ff b1 32 8e 91 08 |k.R..k..|....2...|
|000042a0| 78 d4 b2 69 73 93 6f 62 | 74 61 69 6e 65 64 93 62 |x..is.ob|tained.b|
|000042b0| 90 b8 e3 79 8e a1 8d 61 | 96 04 1d d9 72 65 6c 61 |...y...a|....rela|
|000042c0| 74 65 64 93 6d 65 74 68 | 6f 9a 47 1d 64 93 6f 66 |ted.meth|o.G.d.of|
|000042d0| 93 0c 6e 64 69 6e 67 93 | 61 93 76 90 b8 e3 65 63 |..nding.|a.v...ec|
|000042e0| 74 6f 72 93 b5 79 91 04 | 79 b2 b2 77 69 74 68 93 |tor..y..|y..with.|
|000042f0| 75 6e 69 74 93 b8 31 b2 | 2d 6e 6f 72 6d 93 61 6e |unit..1.|-norm.an|
|00004300| 64 93 b8 6b b5 52 90 13 | c7 79 90 5b d9 b8 6b 8d |d..k.R..|.y.[..k.|
|00004310| 9f 01 7f ff b7 31 8e 91 | 0c 96 be b2 71 75 69 74 |.....1..|....quit|
|00004320| 65 93 6c 61 72 67 65 2e | 91 06 cb 52 54 68 65 93 |e.large.|...RThe.|
|00004330| 70 98 6f 95 b8 e3 77 93 | 65 72 96 04 1d d9 6d 65 |p.o...w.|er....me|
|00004340| 74 68 6f 98 64 93 69 73 | 8e a1 8d 61 70 70 6c 69 |tho.d.is|...appli|
|00004350| 65 64 96 04 26 e0 74 6f | 93 74 68 65 93 6d 61 74 |ed..&.to|.the.mat|
|00004360| 72 69 78 93 b5 52 8d 90 | 13 c7 9f fc 5e ff b4 54 |rix..R..|....^..T|
|00004370| 8e 91 06 60 7d b5 52 90 | 13 c7 b2 2e 91 06 e6 67 |...`}.R.|.......g|
|00004380| 54 91 ff 2a aa 61 6b 69 | 6e 67 93 73 71 75 61 72 |T..*.aki|ng.squar|
|00004390| 65 93 72 6f 90 47 1d 6f | 74 93 6f 66 93 74 68 65 |e.ro.G.o|t.of.the|
|000043a0| 93 65 73 74 69 6d 61 74 | 65 64 93 65 69 67 65 6e |.estimat|ed.eigen|
|000043b0| 9a b8 e3 76 91 ff 71 c7 | 61 6c 75 65 73 93 67 69 |...v..q.|alues.gi|
|000043c0| 76 98 65 73 93 61 93 6c | 6f 98 77 98 65 72 93 62 |v.es.a.l|o.w.er.b|
|000043d0| 90 47 1d 6f 75 6e 64 93 | 74 6f 93 74 68 65 8e a1 |.G.ound.|to.the..|
|000043e0| 8d 32 2d 6e 6f 72 6d 96 | 03 55 55 63 6f 6e 64 69 |.2-norm.|.UUcondi|
|000043f0| 74 69 6f 6e 93 6e 95 b8 | e3 75 6d 93 62 90 47 1d |tion.n..|.um.b.G.|
|00004400| 65 72 96 03 55 55 6f 66 | 93 b5 52 90 13 c7 b2 2e |er..UUof|..R.....|
|00004410| 8e a6 8d 91 0f 00 00 41 | 91 02 38 9b 73 69 6d 70 |.......A|..8.simp|
|00004420| 6c 65 2c 9b 02 71 c7 61 | 6e 64 96 02 38 e3 75 73 |le,..q.a|nd..8.us|
|00004430| 75 61 6c 6c 79 93 72 65 | 6c 69 61 62 6c 65 2c 98 |ually.re|liable,.|
|00004440| 65 73 74 69 6d 61 74 65 | 93 6f 66 93 74 68 65 93 |estimate|.of.the.|
|00004450| 72 61 6e 6b 93 6f 66 93 | 61 93 6d 61 74 72 69 78 |rank.of.|a.matrix|
|00004460| 93 69 73 93 74 6f 93 66 | 61 63 74 6f 72 93 74 68 |.is.to.f|actor.th|
|00004470| 65 93 6d 61 74 72 69 78 | 93 da 41 91 02 38 9b b2 |e.matrix|..A..8..|
|00004480| 75 73 69 6e 67 93 da 51 | 52 43 50 66 61 63 74 6f |using..Q|RCPfacto|
|00004490| 72 28 29 b2 2c 8e a1 8d | 61 6e 64 96 02 da 5a 74 |r().,...|and...Zt|
|000044a0| 68 65 6e 93 74 6f 93 63 | 6f 75 6e 9a b8 e3 74 93 |hen.to.c|oun...t.|
|000044b0| 74 68 65 93 6e 98 75 6d | 98 62 90 47 1d 65 72 93 |the.n.um|.b.G.er.|
|000044c0| 6f 66 93 64 69 61 67 6f | 6e 61 6c 93 65 6e 98 74 |of.diago|nal.en.t|
|000044d0| 72 69 65 73 93 6f 66 93 | da 41 9b 02 da 3b b2 67 |ries.of.|.A...;.g|
|000044e0| 72 65 61 74 65 72 93 74 | 68 61 6e 93 61 93 63 65 |reater.t|han.a.ce|
|000044f0| 72 74 61 69 6e 93 74 6f | 6c 65 72 61 6e 63 65 93 |rtain.to|lerance.|
|00004500| 69 6e 93 6d 61 67 6e 69 | 74 75 64 65 2e 91 04 48 |in.magni|tude...H|
|00004510| c9 41 98 6d 6f 72 65 8e | a1 8d 72 65 6c 69 61 62 |.A.more.|..reliab|
|00004520| 6c 65 96 03 55 55 61 70 | 70 72 6f 61 63 90 b8 e3 |le..UUap|proac...|
|00004530| 68 93 69 73 93 74 6f 93 | 75 73 65 93 74 68 65 93 |h.is.to.|use.the.|
|00004540| 53 69 6e 67 75 6c 61 72 | 93 56 91 ff 2a aa 61 6c |Singular|.V..*.al|
|00004550| 75 65 93 44 65 63 6f 6d | 70 90 47 1d 6f 73 69 74 |ue.Decom|p.G.osit|
|00004560| 69 6f 6e 2e 91 04 71 c7 | 53 65 65 93 da 73 76 64 |ion...q.|See..svd|
|00004570| 28 29 b2 2e 8e a4 13 44 | 95 8d 8d db 53 45 45 91 |().....D|....SEE.|
|00004580| 03 55 55 41 4c 53 4f 8e | 8e a6 8d 91 0f 00 00 b2 |.UUALSO.|........|
|00004590| 48 6f 75 73 65 68 6f 6c | 64 65 72 96 03 55 55 72 |Househol|der..UUr|
|000045a0| 6f 75 74 69 6e 65 73 93 | da 68 68 76 65 63 28 29 |outines.|.hhvec()|
|000045b0| b2 2c 93 da 68 68 74 72 | 76 65 63 28 29 b2 2c 93 |.,..hhtr|vec().,.|
|000045c0| da 68 68 74 72 72 6f 77 | 73 28 29 93 b2 61 6e 64 |.hhtrrow|s()..and|
|000045d0| 93 da 68 68 74 72 63 6f | 6c 73 28 29 b2 3b 93 da |..hhtrco|ls().;..|
|000045e0| 73 76 64 28 29 b2 2e 8e | a1 8d 8d db 53 4f 55 52 |svd()...|....SOUR|
|000045f0| 43 45 91 03 55 55 46 49 | 4c 45 3a 91 14 00 02 da |CE..UUFI|LE:.....|
|00004600| 51 52 66 61 63 74 6f 72 | 2e 63 8e 8e 8e 8e 8c 8b |QRfactor|.c......|
|00004610| 00 00 00 4f 00 00 00 00 | 00 00 00 00 00 00 00 00 |...O....|........|
|00004620| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00004630| 00 00 00 00 00 00 00 00 | 00 00 38 dd a0 02 a9 52 |........|..8....R|
|00004640| b3 8d a0 fd 62 ad 4d 8d | 8d 92 01 cb 78 6e de 37 |....b.M.|....xn.7|
|00004650| 39 8e 8e a0 02 84 07 32 | 8d a0 fd 98 0a 15 8d 8d |9......2|........|
|00004660| db 4e 41 4d 45 8e 8e a9 | 10 44 95 8d 91 0f 00 00 |.NAME...|.D......|
|00004670| da 6d 61 6b 65 51 2c 91 | 05 3f fd 6d 61 6b 65 52 |.makeQ,.|.?.makeR|
|00004680| 96 03 55 55 b2 7b 93 65 | 78 70 6c 69 63 69 74 6c |..UU.{.e|xplicitl|
|00004690| 79 93 66 6f 72 6d 93 b5 | 51 93 b2 61 6e 64 93 b5 |y.form..|Q..and..|
|000046a0| 52 91 03 69 1c b2 66 61 | 63 74 6f 72 73 8e 9f 13 |R..i..fa|ctors...|
|000046b0| 44 95 8d 8d db 53 59 4e | 4f 50 53 49 53 8e 8e 9f |D....SYN|OPSIS...|
|000046c0| 18 44 95 8d da 23 69 6e | 63 6c 75 64 65 91 05 3f |.D...#in|clude..?|
|000046d0| fd 22 6d 61 74 72 69 78 | 2e 68 22 8e a4 0c 00 00 |."matrix|.h".....|
|000046e0| 8d 23 69 6e 63 6c 75 64 | 65 91 05 3f fd 22 6d 61 |.#includ|e..?."ma|
|000046f0| 74 72 69 78 32 2e 68 22 | 8e a1 8d 4d 41 54 91 1a |trix2.h"|...MAT..|
|00004700| 3f f1 2a 6d 61 6b 65 51 | 28 51 52 2c 64 69 61 67 |?.*makeQ|(QR,diag|
|00004710| 2c 51 6f 75 74 29 8e a1 | 8d 4d 41 54 91 1a 3f f1 |,Qout)..|.MAT..?.|
|00004720| 2a 51 52 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 64 |*QR;...V|EC..?.*d|
|00004730| 69 61 67 3b 8e a1 8d 4d | 41 54 91 1a 3f f1 2a 51 |iag;...M|AT..?.*Q|
|00004740| 6f 75 74 3b 8e 9f 1a 44 | 95 8d 4d 41 54 91 1a 3f |out;...D|..MAT..?|
|00004750| f1 2a 6d 61 6b 65 52 28 | 51 52 2c 52 6f 75 74 29 |.*makeR(|QR,Rout)|
|00004760| 8e a1 8d 4d 41 54 91 1a | 3f f1 2a 51 52 2c 91 05 |...MAT..|?.*QR,..|
|00004770| 3f fd 2a 52 6f 75 74 3b | 8e 9f 1d 44 95 8d 8d db |?.*Rout;|...D....|
|00004780| 44 45 53 43 52 49 50 54 | 49 4f 4e 8e 8e a6 8d 91 |DESCRIPT|ION.....|
|00004790| 0f 00 00 b2 54 68 65 96 | 02 fc 7d 72 6f 75 74 69 |....The.|..}routi|
|000047a0| 6e 65 93 da 6d 61 6b 65 | 51 28 29 93 b2 65 78 70 |ne..make|Q()..exp|
|000047b0| 6c 69 63 69 74 6c 79 93 | 66 6f 72 6d 73 93 74 68 |licitly.|forms.th|
|000047c0| 65 93 6f 72 74 68 6f 67 | 6f 6e 61 6c 93 b5 51 93 |e.orthog|onal..Q.|
|000047d0| b2 6d 61 74 72 69 78 93 | 6f 66 93 74 68 65 93 b5 |.matrix.|of.the..|
|000047e0| 51 52 91 03 10 44 b2 66 | 61 63 74 6f 72 69 73 61 |QR...D.f|actorisa|
|000047f0| 74 69 6f 6e 93 66 72 6f | 6d 93 74 68 65 93 63 6f |tion.fro|m.the.co|
|00004800| 6d 70 61 63 74 8e a1 8d | 72 65 70 72 65 73 65 6e |mpact...|represen|
|00004810| 90 b8 e3 74 61 74 69 6f | 6e 96 03 0e 88 69 6e 93 |...tatio|n....in.|
|00004820| da 51 52 b2 2c 93 da 64 | 69 61 67 b2 2e 9b 04 5a |.QR.,..d|iag....Z|
|00004830| 2e 54 68 65 93 72 65 73 | 75 6c 74 93 69 73 93 73 |.The.res|ult.is.s|
|00004840| 74 6f 72 65 64 93 69 6e | 93 da 51 6f 75 74 b2 2e |tored.in|..Qout..|
|00004850| 98 54 68 69 73 93 72 6f | 75 74 69 6e 65 93 6d 61 |.This.ro|utine.ma|
|00004860| 90 b8 e3 79 93 6e 6f 74 | 93 62 90 47 1d 65 93 75 |...y.not|.b.G.e.u|
|00004870| 73 65 64 93 74 6f 93 66 | 6f 72 6d 93 da 51 6f 75 |sed.to.f|orm..Qou|
|00004880| 74 93 dc 69 6e 91 03 52 | c4 73 69 74 75 b2 2e 8e |t..in..R|.situ...|
|00004890| a6 8d 91 0f 00 00 54 68 | 65 96 03 36 1c 72 6f 75 |......Th|e..6.rou|
|000048a0| 74 69 6e 65 93 da 6d 61 | 6b 65 52 28 29 93 b2 65 |tine..ma|keR()..e|
|000048b0| 78 70 6c 69 63 69 74 6c | 79 93 66 6f 72 6d 73 93 |xplicitl|y.forms.|
|000048c0| 74 68 65 93 75 70 70 90 | 47 1d 65 72 93 74 72 69 |the.upp.|G.er.tri|
|000048d0| 61 6e 67 75 6c 61 72 93 | b5 52 9b 03 49 e3 b2 6d |angular.|.R..I..m|
|000048e0| 61 74 72 69 78 93 6f 66 | 93 74 68 65 93 b5 51 52 |atrix.of|.the..QR|
|000048f0| 98 b2 66 61 63 74 6f 72 | 69 73 61 74 69 6f 6e 2e |..factor|isation.|
|00004900| 91 04 67 5f 54 68 65 93 | 72 65 73 75 6c 74 8e a1 |..g_The.|result..|
|00004910| 8d 69 73 96 03 c0 16 73 | 74 6f 72 65 64 93 69 6e |.is....s|tored.in|
|00004920| 93 da 52 6f 75 74 b2 2e | 9b 05 b2 0a 54 68 69 73 |..Rout..|....This|
|00004930| 93 72 6f 75 74 69 6e 65 | 93 6d 61 90 b8 e3 79 93 |.routine|.ma...y.|
|00004940| 62 90 47 1d 65 93 75 73 | 65 64 93 dc 69 6e 91 03 |b.G.e.us|ed..in..|
|00004950| f6 1e 73 69 74 75 91 00 | c4 63 b2 3b 91 03 f5 76 |..situ..|.c.;...v|
|00004960| 74 68 61 74 93 69 73 2c | 91 03 da c6 77 69 74 68 |that.is,|....with|
|00004970| 93 da 51 52 96 05 3f fd | 3d 3d 93 52 6f 75 74 b2 |..QR..?.|==.Rout.|
|00004980| 2e 98 28 41 63 74 75 61 | 6c 6c 79 96 03 c0 16 74 |..(Actua|lly....t|
|00004990| 68 65 93 72 6f 75 74 69 | 6e 65 93 6a 75 73 74 8e |he.routi|ne.just.|
|000049a0| a1 8d 7a 65 72 6f 73 96 | 03 55 55 74 68 65 93 73 |..zeros.|.UUthe.s|
|000049b0| 74 72 69 63 74 6c 79 93 | 6c 6f 95 b8 e3 77 93 65 |trictly.|lo...w.e|
|000049c0| 72 96 03 55 55 74 72 69 | 61 6e 67 75 6c 61 72 93 |r..UUtri|angular.|
|000049d0| 68 61 6c 66 93 6f 66 93 | da 51 52 b2 2e 29 8e a6 |half.of.|.QR..)..|
|000049e0| 8d 91 0f 00 00 49 66 96 | 03 55 55 da 51 6f 75 74 |.....If.|.UU.Qout|
|000049f0| 93 b2 6f 72 93 da 52 6f | 75 74 93 b2 69 73 93 4e |..or..Ro|ut..is.N|
|00004a00| 55 4c 4c 93 6f 72 93 74 | 6f 90 47 1d 6f 93 73 6d |ULL.or.t|o.G.o.sm|
|00004a10| 61 6c 6c 93 74 6f 93 63 | 6f 6e 90 b8 e3 74 61 69 |all.to.c|on...tai|
|00004a20| 6e 93 74 68 65 93 72 65 | 73 75 6c 74 93 74 68 65 |n.the.re|sult.the|
|00004a30| 6e 93 61 93 6e 65 77 93 | 6d 61 74 72 69 78 93 69 |n.a.new.|matrix.i|
|00004a40| 73 93 63 72 65 61 74 65 | 64 93 61 6e 64 93 72 65 |s.create|d.and.re|
|00004a50| 74 75 72 6e 65 64 2e 8e | 9f 13 44 95 8d 8d db 45 |turned..|..D....E|
|00004a60| 58 41 4d 50 4c 45 8e 8e | 9f 1a 44 95 8d da 4d 41 |XAMPLE..|..D...MA|
|00004a70| 54 91 14 ff f4 2a 41 2c | 96 05 3f fd 2a 51 52 2c |T....*A,|..?.*QR,|
|00004a80| 93 2a 51 2c 93 2a 52 3b | 8e a1 8d 56 45 43 91 14 |.*Q,.*R;|...VEC..|
|00004a90| ff f4 2a 64 69 61 67 3b | 8e a1 8d 91 0a 7f fa 2e |..*diag;|........|
|00004aa0| 2e 2e 2e 2e 2e 8e a1 8d | 64 69 61 67 96 05 3f fd |........|diag..?.|
|00004ab0| 3d 93 67 65 74 5f 76 65 | 63 28 41 2d 3e 6d 29 3b |=.get_ve|c(A->m);|
|00004ac0| 8e a1 8d 51 52 96 05 3f | fd 3d 93 63 70 5f 6d 61 |...QR..?|.=.cp_ma|
|00004ad0| 74 28 41 2c 4d 4e 55 4c | 4c 29 3b 8e a1 8d 51 52 |t(A,MNUL|L);...QR|
|00004ae0| 66 61 63 74 6f 72 28 51 | 52 2c 64 69 61 67 29 3b |factor(Q|R,diag);|
|00004af0| 8e a1 8d 51 96 05 3f fd | 3d 93 6d 61 6b 65 51 28 |...Q..?.|=.makeQ(|
|00004b00| 51 52 2c 64 69 61 67 2c | 4d 4e 55 4c 4c 29 3b 8e |QR,diag,|MNULL);.|
|00004b10| a1 8d 52 96 05 3f fd 3d | 93 6d 61 6b 65 52 28 51 |..R..?.=|.makeR(Q|
|00004b20| 52 2c 4d 4e 55 4c 4c 29 | 3b 8e a1 8d 2f 2a 96 05 |R,MNULL)|;.../*..|
|00004b30| 3f fd 6d 61 6b 65 52 28 | 51 52 2c 51 52 29 3b 93 |?.makeR(|QR,QR);.|
|00004b40| 72 65 70 6c 61 63 65 73 | 93 51 52 93 77 69 74 68 |replaces|.QR.with|
|00004b50| 93 74 68 65 93 52 93 6d | 61 74 72 69 78 93 2a 2f |.the.R.m|atrix.*/|
|00004b60| 8e 9f 1d 44 95 8d 8d db | 53 4f 55 52 43 45 91 03 |...D....|SOURCE..|
|00004b70| 55 55 46 49 4c 45 3a 91 | 14 00 02 da 51 52 66 61 |UUFILE:.|....QRfa|
|00004b80| 63 74 6f 72 2e 63 8e 8e | 8e 8e 8c 8b 00 00 00 50 |ctor.c..|.......P|
|00004b90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00004ba0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00004bb0| 00 00 00 00 00 00 48 43 | a0 02 a9 52 b3 8d a0 fd |......HC|...R....|
|00004bc0| 62 ad 4d 8d 8d b2 38 30 | 92 00 99 14 c1 de 43 48 |b.M...80|......CH|
|00004bd0| 41 50 54 45 52 96 03 55 | 55 35 2e 91 07 c7 1c 44 |APTER..U|U5.....D|
|00004be0| 45 4e 53 45 93 4d 41 9b | ff 2a aa 54 52 49 58 93 |ENSE.MA.|.*.TRIX.|
|00004bf0| 46 91 fe e3 8d 41 90 b8 | e3 43 54 4f 52 49 53 41 |F....A..|.CTORISA|
|00004c00| 98 54 49 4f 4e 93 4f 50 | 45 52 41 98 54 49 4f 4e |.TION.OP|ERA.TION|
|00004c10| 53 8e 8e a0 02 84 07 32 | 8d a0 fd 98 0a 15 8d 8d |S......2|........|
|00004c20| db 4e 41 4d 45 8e 8e 9f | 10 44 95 8d 91 0f 00 00 |.NAME...|.D......|
|00004c30| da 67 69 76 65 6e 73 2c | 96 05 3f fd 72 6f 74 8d |.givens,|..?.rot.|
|00004c40| 91 00 a1 45 89 00 00 66 | 66 00 03 26 66 8e 9b 03 |...E...f|f..&f...|
|00004c50| c7 ab 63 6f 6c 73 2c 93 | 72 6f 74 8d 91 00 a1 45 |..cols,.|rot....E|
|00004c60| 89 00 00 66 66 00 03 26 | 66 8e 98 72 6f 77 73 2c |...ff..&|f..rows,|
|00004c70| 93 72 6f 74 8d 91 00 a1 | 45 89 00 00 66 66 00 03 |.rot....|E...ff..|
|00004c80| 26 66 8e 98 76 65 63 96 | 03 55 55 b2 7b 93 47 69 |&f..vec.|.UU.{.Gi|
|00004c90| 76 90 b8 e3 65 6e 73 27 | 93 72 6f 74 61 74 69 6f |v...ens'|.rotatio|
|00004ca0| 6e 73 93 72 6f 75 74 69 | 6e 65 73 8e 9f 13 44 95 |ns.routi|nes...D.|
|00004cb0| 8d 8d db 53 59 4e 4f 50 | 53 49 53 8e 8e 9f 16 9c |...SYNOP|SIS.....|
|00004cc0| 9a 8d da 23 69 6e 63 6c | 75 64 65 91 05 3f fd 22 |...#incl|ude..?."|
|00004cd0| 6d 61 74 72 69 78 2e 68 | 22 8e a4 0c 00 00 8d 23 |matrix.h|"......#|
|00004ce0| 69 6e 63 6c 75 64 65 91 | 05 3f fd 22 6d 61 74 72 |include.|.?."matr|
|00004cf0| 69 78 32 2e 68 22 8e a1 | 8d 76 6f 69 64 91 14 ff |ix2.h"..|.void...|
|00004d00| f4 67 69 76 65 6e 73 28 | 78 2c 79 2c 63 2c 73 29 |.givens(|x,y,c,s)|
|00004d10| 8e a1 8d 64 6f 75 62 6c | 65 91 0a 7f fa 78 2c 91 |...doubl|e....x,.|
|00004d20| 05 3f fd 79 3b 8e a1 8d | 64 6f 75 62 6c 65 91 0a |.?.y;...|double..|
|00004d30| 7f fa 2a 63 2c 91 05 3f | fd 2a 73 3b 8e a9 18 32 |..*c,..?|.*s;...2|
|00004d40| 9c 8d 4d 41 54 91 1a 3f | f1 2a 72 6f 74 5f 63 6f |..MAT..?|.*rot_co|
|00004d50| 6c 73 28 41 2c 69 2c 6b | 2c 63 2c 73 2c 6f 75 74 |ls(A,i,k|,c,s,out|
|00004d60| 29 8e a1 8d 4d 41 54 91 | 1a 3f f1 2a 41 2c 91 05 |)...MAT.|.?.*A,..|
|00004d70| 3f fd 2a 6f 75 74 3b 8e | a1 8d 69 6e 74 91 1a 3f |?.*out;.|..int..?|
|00004d80| f1 69 2c 91 05 3f fd 6b | 3b 8e a1 8d 64 6f 75 62 |.i,..?.k|;...doub|
|00004d90| 6c 65 91 0a 7f fa 63 2c | 91 05 3f fd 73 3b 8e a6 |le....c,|..?.s;..|
|00004da0| 8d 4d 41 54 91 1a 3f f1 | 2a 72 6f 74 5f 72 6f 77 |.MAT..?.|*rot_row|
|00004db0| 73 28 41 2c 69 2c 6b 2c | 63 2c 73 2c 6f 75 74 29 |s(A,i,k,|c,s,out)|
|00004dc0| 8e a1 8d 4d 41 54 91 1a | 3f f1 2a 41 2c 91 05 3f |...MAT..|?.*A,..?|
|00004dd0| fd 2a 6f 75 74 3b 8e a1 | 8d 69 6e 74 91 1a 3f f1 |.*out;..|.int..?.|
|00004de0| 69 2c 91 05 3f fd 6b 3b | 8e a1 8d 64 6f 75 62 6c |i,..?.k;|...doubl|
|00004df0| 65 91 0a 7f fa 63 2c 91 | 05 95 8d 91 0f 00 00 b2 |e....c,.|........|
|00004e00| 55 70 90 47 1d 64 61 74 | 69 6e 67 91 03 55 55 b5 |Up.G.dat|ing..UU.|
|00004e10| 51 52 91 03 69 1c b2 66 | 61 63 74 6f 72 69 73 61 |QR..i..f|actorisa|
|00004e20| 74 69 6f 6e 3a 8e a1 8d | da 4d 41 54 91 14 ff f4 |tion:...|.MAT....|
|00004e30| 2a 41 2c 96 05 3f fd 2a | 51 52 2c 93 2a 51 2c 93 |*A,..?.*|QR,.*Q,.|
|00004e40| 2a 52 3b 8e a6 8d 56 45 | 43 91 14 ff f4 2a 64 69 |*R;...VE|C....*di|
|00004e50| 61 67 2c 96 05 3f fd 2a | 62 65 74 61 2c 93 2a 75 |ag,..?.*|beta,.*u|
|00004e60| 2c 93 2a 76 2c 93 2a 77 | 3b 8e a6 8d 91 0a 7f fa |,.*v,.*w|;.......|
|00004e70| 2e 2e 2e 2e 2e 2e 8e a6 | 8d 51 52 96 05 3f fd 3d |........|.QR..?.=|
|00004e80| 93 63 70 5f 6d 61 74 28 | 41 2c 4d 4e 55 4c 4c 29 |.cp_mat(|A,MNULL)|
|00004e90| 3b 8e a6 8d 51 52 66 61 | 63 74 6f 72 28 51 52 2c |;...QRfa|ctor(QR,|
|00004ea0| 64 69 61 67 2c 62 65 74 | 61 29 3b 8e a6 8d 51 96 |diag,bet|a);...Q.|
|00004eb0| 05 3f fd 3d 93 6d 61 6b | 65 51 28 51 52 2c 64 69 |.?.=.mak|eQ(QR,di|
|00004ec0| 61 67 2c 62 65 74 61 2c | 4d 4e 55 4c 4c 29 3b 8e |ag,beta,|MNULL);.|
|00004ed0| a6 8d 52 96 05 3f fd 3d | 93 6d 61 6b 65 52 28 51 |..R..?.=|.makeR(Q|
|00004ee0| 52 2c 4d 4e 55 4c 4c 29 | 3b 8e a6 8d 91 0a 7f fa |R,MNULL)|;.......|
|00004ef0| 2e 2e 2e 2e 2e 2e 8e a6 | 8d 2f 2a 96 05 3f fd 41 |........|./*..?.A|
|00004f00| 93 3c 2d 93 41 93 2b 93 | 77 2e 76 5e 54 93 2a 2f |.<-.A.+.|w.v^T.*/|
|00004f10| 8e a6 8d 75 96 05 3f fd | 3d 93 67 65 74 5f 76 65 |...u..?.|=.get_ve|
|00004f20| 63 28 51 2d 3e 6d 29 3b | 8e a6 8d 75 96 05 3f fd |c(Q->m);|...u..?.|
|00004f30| 3d 93 76 6d 5f 6d 6c 74 | 28 51 2c 77 2c 75 29 3b |=.vm_mlt|(Q,w,u);|
|00004f40| 8e a6 8d 51 52 75 70 64 | 61 74 65 28 51 2c 52 2c |...QRupd|ate(Q,R,|
|00004f50| 75 2c 76 29 3b 8e 9f 1b | 44 95 8d 8d db 53 4f 55 |u,v);...|D....SOU|
|00004f60| 52 43 45 91 03 55 55 46 | 49 4c 45 3a 91 14 00 02 |RCE..UUF|ILE:....|
|00004f70| da 75 70 64 61 74 65 2e | 63 8e 8e 8e 8e 8c 8b 00 |.update.|c.......|
|00004f80| 00 00 58 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |..X.....|........|
|00004f90| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00004fa0| 00 00 00 00 00 00 00 00 | 00 7f 9e a0 02 a9 52 b3 |........|......R.|
|00004fb0| 8d a0 fd 62 ad 4d 8d 8d | b2 38 38 92 00 99 14 c1 |...b.M..|.88.....|
|00004fc0| de 43 48 41 50 54 45 52 | 96 03 55 55 35 2e 91 07 |.CHAPTER|..UU5...|
|00004fd0| c7 1c 44 45 4e 53 45 93 | 4d 41 9b ff 2a aa 54 52 |..DENSE.|MA..*.TR|
|00004fe0| 49 58 93 46 91 fe e3 8d | 41 90 b8 e3 43 54 4f 52 |IX.F....|A...CTOR|
|00004ff0| 49 53 41 98 54 49 4f 4e | 93 4f 50 45 52 41 98 54 |ISA.TION|.OPERA.T|
|00005000| 49 4f 4e 53 8e 8e a0 02 | 84 07 32 8d a0 fd 98 0a |IONS....|..2.....|
|00005010| 15 8d 8d db 4e 41 4d 45 | 8e 8e 9f 10 44 95 8d 91 |....NAME|....D...|
|00005020| 0f 00 00 da 73 63 68 75 | 72 2c 96 05 3f fd 73 79 |....schu|r,..?.sy|
|00005030| 6d 6d 65 69 67 2c 93 74 | 72 69 65 69 67 96 03 55 |mmeig,.t|rieig..U|
|00005040| 55 b2 7b 93 65 69 67 65 | 6e 90 b8 e3 76 91 ff 71 |U.{.eige|n...v..q|
|00005050| c7 61 6c 75 65 93 72 6f | 75 74 69 6e 65 73 8e 9f |.alue.ro|utines..|
|00005060| 13 44 95 8d 8d db 53 59 | 4e 4f 50 53 49 53 8e 8e |.D....SY|NOPSIS..|
|00005070| 9f 18 20 e3 8d da 23 69 | 6e 63 6c 75 64 65 91 05 |.. ...#i|nclude..|
|00005080| 3f fd 22 6d 61 74 72 69 | 78 2e 68 22 8e a4 0c 00 |?."matri|x.h"....|
|00005090| 00 8d 23 69 6e 63 6c 75 | 64 65 91 05 3f fd 22 6d |..#inclu|de..?."m|
|000050a0| 61 74 72 69 78 32 2e 68 | 22 8e a1 8d 4d 41 54 91 |atrix2.h|"...MAT.|
|000050b0| 1a 3f f1 2a 73 63 68 75 | 72 28 41 2c 51 29 8e a1 |.?.*schu|r(A,Q)..|
|000050c0| 8d 4d 41 54 91 1a 3f f1 | 2a 41 2c 91 05 3f fd 2a |.MAT..?.|*A,..?.*|
|000050d0| 51 3b 8e a9 1a 17 f6 8d | 56 45 43 91 1a 3f f1 2a |Q;......|VEC..?.*|
|000050e0| 73 79 6d 6d 65 69 67 28 | 41 2c 51 2c 6f 75 74 29 |symmeig(|A,Q,out)|
|000050f0| 8e a1 8d 4d 41 54 91 1a | 3f f1 2a 41 2c 91 05 3f |...MAT..|?.*A,..?|
|00005100| fd 2a 51 3b 8e a1 8d 56 | 45 43 91 1a 3f f1 2a 6f |.*Q;...V|EC..?.*o|
|00005110| 75 74 3b 8e a6 8d 56 45 | 43 91 1a 3f f1 2a 74 72 |ut;...VE|C..?.*tr|
|00005120| 69 65 69 67 28 61 2c 62 | 2c 51 29 8e a1 8d 56 45 |ieig(a,b|,Q)...VE|
|00005130| 43 91 1a 3f f1 2a 61 2c | 91 05 3f fd 2a 62 3b 8e |C..?.*a,|..?.*b;.|
|00005140| a1 8d 4d 41 54 91 1a 3f | f1 2a 51 3b 8e 9f 1d 17 |..MAT..?|.*Q;....|
|00005150| f6 8d 8d db 44 45 53 43 | 52 49 50 54 49 4f 4e 8e |....DESC|RIPTION.|
|00005160| 8e a9 10 44 95 8d 91 0f | 00 00 b2 54 68 65 96 04 |...D....|...The..|
|00005170| 05 a8 72 6f 75 74 69 6e | 65 93 da 73 63 68 75 72 |..routin|e..schur|
|00005180| 28 29 93 b2 63 6f 6d 70 | 75 74 65 73 93 74 68 65 |()..comp|utes.the|
|00005190| 93 52 65 61 6c 93 53 63 | 95 b8 e3 68 93 75 72 96 |.Real.Sc|...h.ur.|
|000051a0| 04 05 a8 64 65 63 6f 6d | 70 90 47 1d 6f 73 69 74 |...decom|p.G.osit|
|000051b0| 69 6f 6e 93 6f 66 93 74 | 68 65 93 6d 61 74 72 69 |ion.of.t|he.matri|
|000051c0| 78 93 da 41 b2 2e 93 54 | 68 61 74 93 69 73 2c 91 |x..A...T|hat.is,.|
|000051d0| 04 31 bc 69 74 93 63 6f | 6d 70 75 74 65 73 93 61 |.1.it.co|mputes.a|
|000051e0| 8e a1 8d 62 6c 6f 95 47 | 1d 63 90 b8 e3 6b 9b 03 |...blo.G|.c...k..|
|000051f0| 55 55 75 70 70 93 65 72 | 98 74 72 69 61 6e 67 75 |UUupp.er|.triangu|
|00005200| 6c 61 72 98 6d 61 74 72 | 69 78 98 b5 54 91 04 b8 |lar.matr|ix..T...|
|00005210| e4 b2 61 6e 64 98 61 6e | 98 6f 72 74 68 6f 67 6f |..and.an|.orthogo|
|00005220| 6e 61 6c 98 6d 61 74 72 | 69 78 98 b5 51 98 b2 73 |nal.matr|ix..Q..s|
|00005230| 75 63 90 b8 e3 68 98 74 | 68 61 74 8e a4 15 d3 61 |uc...h.t|hat....a|
|00005240| 8d 92 00 d0 89 2a b5 51 | 8d 9f fb de ff b4 54 8e |.....*.Q|......T.|
|00005250| 91 06 4c b6 b5 41 51 96 | 02 c7 18 b2 3d 93 b5 54 |..L..AQ.|....=..T|
|00005260| 91 00 d5 56 3a 8e a1 8d | b2 54 68 65 96 04 11 e6 |...V:...|.The....|
|00005270| 6d 61 74 72 69 78 93 b5 | 54 91 05 75 75 b2 68 61 |matrix..|T..uu.ha|
|00005280| 73 93 64 69 61 67 6f 6e | 61 6c 93 62 6c 6f 90 47 |s.diagon|al.blo.G|
|00005290| 1d 63 9a b8 e3 6b 73 93 | 6f 66 93 73 69 7a 65 73 |.c...ks.|of.sizes|
|000052a0| 93 31 96 02 b6 96 b8 02 | 93 b2 31 96 04 11 e6 61 |.1......|..1....a|
|000052b0| 6e 64 93 32 96 02 b6 96 | b8 02 93 b2 32 2e 91 06 |nd.2....|....2...|
|000052c0| a7 7b 54 68 65 96 04 11 | e6 65 69 67 65 6e 98 76 |.{The...|.eigen.v|
|000052d0| 91 ff 71 c7 61 6c 75 65 | 73 93 6f 66 93 74 68 65 |..q.alue|s.of.the|
|000052e0| 73 65 93 64 69 61 67 6f | 6e 61 6c 93 62 6c 6f 90 |se.diago|nal.blo.|
|000052f0| 47 1d 63 98 6b 73 93 61 | 72 65 8e a4 0c 00 00 8d |G.c.ks.a|re......|
|00005300| 74 68 65 96 03 42 66 65 | 69 67 65 6e 9a b8 e3 76 |the..Bfe|igen...v|
|00005310| 91 ff 71 c7 61 6c 75 65 | 73 93 6f 66 93 74 68 65 |..q.alue|s.of.the|
|00005320| 93 6f 72 69 67 69 6e 61 | 6c 93 b5 41 93 b2 6d 61 |.origina|l..A..ma|
|00005330| 74 72 69 78 2e 91 04 6b | 77 54 68 65 93 61 6c 67 |trix...k|wThe.alg|
|00005340| 6f 72 69 74 68 6d 93 75 | 73 65 64 93 74 6f 93 0c |orithm.u|sed.to..|
|00005350| 6e 64 93 74 68 65 93 65 | 69 67 65 6e 98 76 91 ff |nd.the.e|igen.v..|
|00005360| 71 c7 61 6c 75 65 73 93 | 6f 66 93 b5 41 93 b2 69 |q.alues.|of..A..i|
|00005370| 73 93 74 68 65 93 46 91 | ff 2a aa 72 61 6e 63 69 |s.the.F.|.*.ranci|
|00005380| 73 93 b5 51 52 8e a1 8d | b2 61 6c 67 6f 72 69 74 |s..QR...|.algorit|
|00005390| 68 6d 2e 91 04 59 c7 54 | 68 69 73 96 03 0d 54 61 |hm...Y.T|his...Ta|
|000053a0| 6c 67 6f 72 69 74 68 6d | 93 69 73 93 64 65 73 63 |lgorithm|.is.desc|
|000053b0| 72 69 62 90 47 1d 65 64 | 93 69 6e 93 dc 4d 61 74 |rib.G.ed|.in..Mat|
|000053c0| 72 69 78 91 03 51 a9 43 | 6f 6d 70 75 74 61 74 69 |rix..Q.C|omputati|
|000053d0| 6f 6e 73 91 03 df 75 b2 | 62 90 b8 e3 79 93 47 2e |ons...u.|b...y.G.|
|000053e0| 93 47 6f 6c 75 62 93 61 | 6e 64 93 43 2e 93 76 91 |.Golub.a|nd.C..v.|
|000053f0| ff 71 c7 61 6e 93 4c 6f | 61 6e 2c 91 03 1b bb 70 |.q.an.Lo|an,....p|
|00005400| 70 2e 93 32 33 31 7b 32 | 33 36 2c 8e a1 8d 31 73 |p..231{2|36,...1s|
|00005410| 74 96 03 55 55 45 64 69 | 74 69 6f 6e 93 28 31 39 |t..UUEdi|tion.(19|
|00005420| 38 33 29 2c 93 70 70 2e | 93 33 37 37 7b 33 38 31 |83),.pp.|.377{381|
|00005430| 2c 93 32 6e 64 93 45 64 | 69 74 69 6f 6e 93 28 31 |,.2nd.Ed|ition.(1|
|00005440| 39 38 39 29 2e 8e a6 8d | 91 0f 00 00 54 68 65 96 |989)....|....The.|
|00005450| 04 0c 93 6d 61 74 72 69 | 78 93 da 41 9b 04 0c 64 |...matri|x..A...d|
|00005460| b2 69 73 93 6f 95 b8 e3 | 76 93 65 72 77 72 69 74 |.is.o...|v.erwrit|
|00005470| 74 65 6e 96 04 0c 93 77 | 69 74 68 93 b5 54 91 01 |ten....w|ith..T..|
|00005480| 63 8f b2 2c 91 04 3a 62 | 61 6e 64 93 69 66 93 da |c..,..:b|and.if..|
|00005490| 51 98 b2 69 73 93 6e 6f | 74 93 4e 55 4c 4c 98 61 |Q..is.no|t.NULL.a|
|000054a0| 6e 64 93 74 68 65 93 63 | 6f 72 72 65 63 74 93 73 |nd.the.c|orrect.s|
|000054b0| 69 7a 65 2c 91 04 3a 62 | 74 68 65 6e 93 74 68 65 |ize,..:b|then.the|
|000054c0| 93 b5 51 93 b2 6d 61 74 | 72 69 78 93 69 73 8e a1 |..Q..mat|rix.is..|
|000054d0| 8d 73 74 6f 72 65 64 96 | 03 55 55 69 6e 93 69 74 |.stored.|.UUin.it|
|000054e0| 2e 8e a6 8d 91 0f 00 00 | 54 68 65 96 02 e7 e3 72 |........|The....r|
|000054f0| 6f 75 74 69 6e 65 93 da | 73 79 6d 6d 65 69 67 28 |outine..|symmeig(|
|00005500| 29 93 b2 63 6f 6d 70 75 | 74 65 73 93 74 68 65 93 |)..compu|tes.the.|
|00005510| 65 69 67 65 6e 90 b8 e3 | 76 91 ff 71 c7 61 6c 75 |eigen...|v..q.alu|
|00005520| 65 73 93 6f 66 93 61 93 | dc 73 79 6d 6d 65 74 72 |es.of.a.|.symmetr|
|00005530| 69 63 91 03 78 99 b2 6d | 61 74 72 69 78 2e 91 04 |ic..x..m|atrix...|
|00005540| 4d 4c 49 74 93 61 6c 73 | 6f 93 63 6f 6d 70 75 74 |MLIt.als|o.comput|
|00005550| 65 73 93 61 6e 93 6f 72 | 74 68 6f 67 6f 6e 61 6c |es.an.or|thogonal|
|00005560| 8e a1 8d 6d 61 74 72 69 | 78 96 03 55 55 b5 51 93 |...matri|x..UU.Q.|
|00005570| b2 73 75 63 90 b8 e3 68 | 93 74 68 61 74 8e a1 8d |.suc...h|.that...|
|00005580| 92 00 d1 ca 80 b5 51 8d | 9f fb de ff b4 54 8e 91 |......Q.|.....T..|
|00005590| 06 4c b6 b5 41 51 96 02 | c7 18 b2 3d 93 03 8e 9f |.L..AQ..|...=....|
|000055a0| 11 e5 3a 8d 77 68 65 72 | 65 96 03 a4 f2 03 93 69 |..:.wher|e......i|
|000055b0| 73 93 74 68 65 93 64 69 | 61 67 6f 6e 61 6c 93 6d |s.the.di|agonal.m|
|000055c0| 61 74 72 69 78 93 6f 66 | 93 65 69 67 65 6e 9a b8 |atrix.of|.eigen..|
|000055d0| e3 76 91 ff 71 c7 61 6c | 75 65 73 2e 91 05 60 9d |.v..q.al|ues...`.|
|000055e0| 54 68 65 93 61 6c 67 6f | 72 69 74 68 6d 93 75 73 |The.algo|rithm.us|
|000055f0| 65 64 93 74 6f 93 0c 6e | 64 93 74 68 65 93 65 69 |ed.to..n|d.the.ei|
|00005600| 67 65 6e 98 76 91 ff 71 | c7 61 6c 75 65 73 93 6f |gen.v..q|.alues.o|
|00005610| 66 93 b5 41 93 b2 63 6f | 6e 73 69 73 74 73 93 6f |f..A..co|nsists.o|
|00005620| 66 8e a1 8d 63 6f 6e 95 | b8 e3 76 93 65 72 73 69 |f...con.|..v.ersi|
|00005630| 6f 6e 96 03 d0 ee 74 6f | 93 73 79 6d 6d 65 74 72 |on....to|.symmetr|
|00005640| 69 63 93 48 65 73 73 65 | 6e 90 b8 e3 62 90 47 1d |ic.Hesse|n...b.G.|
|00005650| 65 72 67 93 28 73 79 6d | 6d 65 74 72 69 63 93 74 |erg.(sym|metric.t|
|00005660| 72 69 64 69 61 67 6f 6e | 61 6c 29 93 66 6f 72 6d |ridiagon|al).form|
|00005670| 2c 91 03 ef d5 61 6e 64 | 93 74 68 65 6e 93 61 70 |,....and|.then.ap|
|00005680| 70 6c 79 69 6e 67 93 da | 74 72 69 65 69 67 28 29 |plying..|trieig()|
|00005690| 93 b2 74 6f 93 6f 62 74 | 61 69 6e 8e a1 8d 74 68 |..to.obt|ain...th|
|000056a0| 65 96 03 55 55 65 69 67 | 65 6e 90 b8 e3 76 91 ff |e..UUeig|en...v..|
|000056b0| 71 c7 61 6c 75 65 73 93 | 6f 66 93 74 68 65 93 74 |q.alues.|of.the.t|
|000056c0| 72 69 64 69 61 67 6f 6e | 61 6c 93 6d 61 74 72 69 |ridiagon|al.matri|
|000056d0| 78 2e 8e a6 8d 91 0f 00 | 00 54 68 65 96 04 32 8c |x.......|.The..2.|
|000056e0| 65 69 67 65 6e 9a b8 e3 | 76 91 ff 71 c7 61 6c 75 |eigen...|v..q.alu|
|000056f0| 65 73 93 61 72 65 93 73 | 74 6f 72 65 64 93 69 6e |es.are.s|tored.in|
|00005700| 93 da 6f 75 74 93 b2 70 | 72 6f 98 76 69 64 65 64 |..out..p|ro.vided|
|00005710| 93 69 74 93 69 73 93 6e | 6f 74 93 4e 55 4c 4c 91 |.it.is.n|ot.NULL.|
|00005720| 04 32 53 61 6e 64 93 69 | 73 93 73 75 0e 63 69 65 |.2Sand.i|s.su.cie|
|00005730| 6e 98 74 6c 79 93 6c 61 | 72 67 65 93 74 6f 93 63 |n.tly.la|rge.to.c|
|00005740| 6f 6e 98 74 61 69 6e 93 | 61 6c 6c 93 74 68 65 8e |on.tain.|all.the.|
|00005750| a1 8d 65 69 67 65 6e 95 | b8 e3 76 91 ff 71 c7 61 |..eigen.|..v..q.a|
|00005760| 6c 75 65 73 2e 91 04 71 | c7 54 68 65 9b 03 55 55 |lues...q|.The..UU|
|00005770| 76 93 65 63 74 6f 72 98 | 63 6f 6e 93 74 61 69 6e |v.ector.|con.tain|
|00005780| 69 6e 67 98 74 68 65 98 | 65 69 67 65 6e 93 76 91 |ing.the.|eigen.v.|
|00005790| ff 71 c7 61 6c 75 65 73 | 98 69 73 98 72 65 74 75 |.q.alues|.is.retu|
|000057a0| 72 6e 65 64 2e 91 04 71 | c7 54 68 65 98 6d 61 74 |rned...q|.The.mat|
|000057b0| 72 69 78 98 da 41 98 b2 | 69 73 98 dc 6e 6f 74 91 |rix..A..|is..not.|
|000057c0| 04 48 2d b2 6f 93 76 93 | 65 72 77 72 69 74 74 65 |.H-.o.v.|erwritte|
|000057d0| 6e 2e 8e a6 8d 91 0f 00 | 00 54 68 65 96 03 55 55 |n.......|.The..UU|
|000057e0| 72 6f 75 74 69 6e 65 93 | da 74 72 69 65 69 67 28 |routine.|.trieig(|
|000057f0| 29 93 b2 63 6f 6d 70 75 | 74 65 73 93 74 68 65 93 |)..compu|tes.the.|
|00005800| 65 69 67 65 6e 90 b8 e3 | 76 91 ff 71 c7 61 6c 75 |eigen...|v..q.alu|
|00005810| 65 73 93 6f 66 93 74 68 | 65 93 73 79 6d 6d 65 74 |es.of.th|e.symmet|
|00005820| 72 69 63 93 74 72 69 64 | 69 61 67 6f 6e 61 6c 93 |ric.trid|iagonal.|
|00005830| 6d 61 74 72 69 78 8e 9f | 34 34 19 8d 8d 8d 8d 92 |matrix..|44......|
|00005840| 00 9f 7d 66 b5 54 91 04 | 2a a7 b2 3d 8d 8d 9f d9 |..}f.T..|*..=....|
|00005850| e6 4e 8d 91 02 c7 18 ab | 32 8e 9f 11 99 a6 8d 91 |.N......|2.......|
|00005860| 02 c7 18 36 8e a4 06 00 | 04 8d 91 02 c7 18 36 8e |...6....|......6.|
|00005870| a1 8d 91 02 c7 18 36 8e | a1 8d 91 02 c7 18 36 8e |......6.|......6.|
|00005880| a1 8d 91 02 c7 18 36 8e | a1 8d 91 02 c7 18 36 8e |......6.|......6.|
|00005890| 9f 06 66 69 8d 91 02 c7 | 18 34 8e 8e 8d 9f e1 3e |..fi....|.4.....>|
|000058a0| 90 8d 8d 91 09 71 c4 b5 | 61 8d 9f 01 7f ff b1 30 |.....q..|a......0|
|000058b0| 8e 8e 8d 91 1d b6 aa b5 | 62 8d 9f 01 7f ff b1 30 |........|b......0|
|000058c0| 8e 8e 8e a1 8d 8d 91 09 | f1 07 b5 62 8d 9f 01 7f |........|...b....|
|000058d0| ff b1 30 8e 8e 8d 91 1d | 37 67 b5 61 8d 9f 01 7f |..0.....|7g.a....|
|000058e0| ff b1 31 8e 8e 8d 91 32 | 6e cf b5 62 8d 9f 01 7f |..1....2|n..b....|
|000058f0| ff b1 31 8e 8e 8e a4 12 | a7 d4 8d 8d 91 1d b6 aa |..1.....|........|
|00005900| b5 62 8d 9f 01 7f ff b1 | 31 8e 8e 8d 91 31 ef 8c |.b......|1....1..|
|00005910| b5 61 8d 9f 01 7f ff b1 | 32 8e 8e 8d 8d 8d 9f f9 |.a......|2.......|
|00005920| 00 00 8d 91 4b 5c e4 b2 | 2e 8e 8e 8d 91 4f 40 71 |....K\..|.....O@q|
|00005930| 9f fc 00 00 2e 8e 8d 91 | 53 23 fe 9f ff 00 00 2e |........|S#......|
|00005940| 8e 8e 8e 8e a1 8d 8d 8d | 8d 9f f9 00 00 8d 91 31 |........|.......1|
|00005950| 8b 42 2e 8e 8e 8d 91 35 | 6e cf 9f fc 00 00 2e 8e |.B.....5|n.......|
|00005960| 8d 91 39 52 5c 9f ff 00 | 00 2e 8e 8e 8e 8d 8d 8d |..9R\...|........|
|00005970| 9f f9 00 00 8d 91 4b 5c | e4 2e 8e 8e 8d 91 4f 40 |......K\|......O@|
|00005980| 71 9f fc 00 00 2e 8e 8d | 91 53 23 fe 9f ff 00 00 |q.......|.S#.....|
|00005990| 2e 8e 8e 8e 8d 91 65 1f | 90 b5 62 8d 9f 01 7f ff |......e.|..b.....|
|000059a0| b4 6e b7 00 b1 32 8e 8e | 8e 9f 0c 00 00 8d 8d 91 |.n...2..|........|
|000059b0| 46 a7 b1 b5 62 8d 9f 01 | 7f ff b4 6e b7 00 b1 32 |F...b...|...n...2|
|000059c0| 8e 8e 8d 91 64 a0 4d b5 | 61 8d 9f 01 7f ff b4 6e |....d.M.|a......n|
|000059d0| b7 00 b1 31 8e 8e 8e 8e | 8d 9f d9 e6 4e 8d 91 79 |...1....|....N..y|
|000059e0| 97 6f ab 33 8e 9f 11 99 | a6 8d 91 79 97 6f 37 8e |.o.3....|...y.o7.|
|000059f0| a4 06 00 04 8d 91 79 97 | 6f 37 8e a1 8d 91 79 97 |......y.|o7....y.|
|00005a00| 6f 37 8e a1 8d 91 79 97 | 6f 37 8e a1 8d 91 79 97 |o7....y.|o7....y.|
|00005a10| 6f 37 8e a1 8d 91 79 97 | 6f 37 8e 9f 06 66 69 8d |o7....y.|o7...fi.|
|00005a20| 91 79 97 6f 35 8e 8e 8e | 92 00 81 ec c3 b5 3a 8e |.y.o5...|......:.|
|00005a30| 8e 8e 8d 8d 8d 8d 8f 01 | b2 28 35 2e 34 29 8e 8e |........|.(5.4)..|
|00005a40| 8e 8e 8e 9f 33 ec fc 8d | 54 68 65 96 03 cb 42 61 |....3...|The...Ba|
|00005a50| 6c 67 6f 72 69 74 68 6d | 93 75 73 65 64 93 69 73 |lgorithm|.used.is|
|00005a60| 93 61 93 5c 63 9a b8 e3 | 68 61 73 69 6e 67 22 93 |.a.\c...|hasing".|
|00005a70| 74 65 63 98 68 6e 69 71 | 75 65 93 64 65 73 63 72 |tec.hniq|ue.descr|
|00005a80| 69 62 90 47 1d 65 64 93 | 69 6e 93 dc 4d 61 74 72 |ib.G.ed.|in..Matr|
|00005a90| 69 78 91 04 00 65 43 6f | 6d 70 75 74 61 74 69 6f |ix...eCo|mputatio|
|00005aa0| 6e 73 b2 2c 9b 03 e8 bd | 70 70 2e 93 32 37 38 7b |ns.,....|pp..278{|
|00005ab0| 32 38 31 2c 98 31 73 74 | 93 45 64 69 74 69 6f 6e |281,.1st|.Edition|
|00005ac0| 2c 8e a1 8d 70 70 2e 96 | 03 82 d3 34 32 31 7b 34 |,...pp..|...421{4|
|00005ad0| 32 34 2c 91 03 8e 33 32 | 6e 64 93 45 64 69 74 69 |24,...32|nd.Editi|
|00005ae0| 6f 6e 2e 91 04 fa 41 49 | 74 93 61 6c 73 6f 93 61 |on....AI|t.also.a|
|00005af0| 63 63 75 6d 9a b8 e3 75 | 6c 61 74 65 73 93 74 68 |ccum...u|lates.th|
|00005b00| 65 93 6d 61 74 72 69 78 | 93 b5 51 93 b2 73 75 63 |e.matrix|..Q..suc|
|00005b10| 98 68 93 74 68 61 74 93 | b5 51 8d 9f fc 5e ff b4 |.h.that.|.Q...^..|
|00005b20| 54 8e 91 06 4c b6 b5 54 | 91 01 63 8f 51 93 b2 69 |T...L..T|..c.Q..i|
|00005b30| 73 93 64 69 61 67 6f 6e | 61 6c 2e 91 04 fa 41 54 |s.diagon|al....AT|
|00005b40| 91 ff 2a aa 6f 93 63 6f | 6d 70 75 74 65 93 74 68 |..*.o.co|mpute.th|
|00005b50| 65 8e a1 8d 63 6f 72 72 | 65 63 74 96 03 55 55 b5 |e...corr|ect..UU.|
|00005b60| 51 93 b2 6d 61 74 72 69 | 78 2c 93 da 51 93 b2 73 |Q..matri|x,..Q..s|
|00005b70| 68 6f 75 6c 64 93 62 90 | 47 1d 65 93 69 6e 69 74 |hould.b.|G.e.init|
|00005b80| 69 61 6c 69 73 65 64 93 | 74 6f 93 74 68 65 93 69 |ialised.|to.the.i|
|00005b90| 64 65 6e 95 b8 e3 74 69 | 74 93 79 96 03 55 55 6d |den...ti|t.y..UUm|
|00005ba0| 61 74 72 69 78 93 6f 6e | 93 65 6e 90 b8 e3 74 72 |atrix.on|.en...tr|
|00005bb0| 79 93 74 6f 93 da 74 72 | 69 65 69 67 28 29 b2 2e |y.to..tr|ieig()..|
|00005bc0| 91 04 71 c7 28 53 65 65 | 93 da 69 64 5f 6d 61 74 |..q.(See|..id_mat|
|00005bd0| 28 29 b2 2e 29 8e 8e 8e | 8c 8b 00 00 00 59 00 00 |()..)...|.....Y..|
|00005be0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005bf0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00005c00| 00 00 00 00 81 b2 a0 02 | a9 52 b3 8d a0 fd 62 ad |........|.R....b.|
|00005c10| 4d 8d 8d 92 01 cb 78 6e | de 38 39 8e 8e a0 02 84 |M.....xn|.89.....|
|00005c20| 07 32 8d a0 fd 98 0a 15 | 8d 91 0f 00 00 b2 54 68 |.2......|......Th|
|00005c30| 65 96 03 3e 2a 76 91 ff | 71 c7 61 6c 75 65 73 93 |e..>*v..|q.alues.|
|00005c40| 69 6e 93 74 68 65 93 da | 61 93 b2 61 6e 64 93 da |in.the..|a..and..|
|00005c50| 62 93 b2 76 9a b8 e3 65 | 63 74 6f 72 73 93 61 72 |b..v...e|ctors.ar|
|00005c60| 65 93 6f 98 76 98 65 72 | 77 72 69 74 74 65 6e 2e |e.o.v.er|written.|
|00005c70| 91 04 6a 0e 41 98 74 93 | 74 68 65 93 65 6e 64 93 |..j.A.t.|the.end.|
|00005c80| 6f 66 93 74 68 65 93 72 | 6f 75 74 69 6e 65 2c 91 |of.the.r|outine,.|
|00005c90| 03 42 cc da 61 93 b2 63 | 6f 6e 98 74 61 69 6e 73 |.B..a..c|on.tains|
|00005ca0| 93 74 68 65 93 65 69 67 | 65 6e 98 76 91 ff 71 c7 |.the.eig|en.v..q.|
|00005cb0| 61 6c 75 65 73 2c 8e a4 | 0c 00 00 8d 61 6e 64 96 |alues,..|....and.|
|00005cc0| 03 55 55 74 68 65 93 da | 62 93 b2 76 90 b8 e3 65 |.UUthe..|b..v...e|
|00005cd0| 63 74 6f 72 93 69 73 93 | 7a 65 72 6f 2e 8e a9 10 |ctor.is.|zero....|
|00005ce0| 44 95 8d 91 0f 00 00 49 | 6e 96 03 99 3f 61 6c 6c |D......I|n...?all|
|00005cf0| 93 6f 66 93 74 68 65 93 | 61 62 90 47 1d 6f 95 b8 |.of.the.|ab.G.o..|
|00005d00| e3 76 93 65 96 03 99 3f | 72 6f 75 74 69 6e 65 73 |.v.e...?|routines|
|00005d10| 2c 91 03 aa 39 69 66 93 | 74 68 65 93 6d 61 74 72 |,...9if.|the.matr|
|00005d20| 69 78 93 da 51 9b 03 99 | 2d b2 69 73 93 4e 55 4c |ix..Q...|-.is.NUL|
|00005d30| 4c 98 6f 6e 93 65 6e 90 | b8 e3 74 72 79 91 ff 2a |L.on.en.|..try..*|
|00005d40| aa 2c 91 03 aa 39 74 68 | 65 6e 93 6e 6f 93 63 61 |.,...9th|en.no.ca|
|00005d50| 6c 63 75 6c 61 74 69 6f | 6e 93 6f 66 93 74 68 65 |lculatio|n.of.the|
|00005d60| 93 b5 51 93 b2 6d 61 74 | 72 69 63 65 73 93 69 73 |..Q..mat|rices.is|
|00005d70| 8e a1 8d 70 9a 47 1d 65 | 72 66 6f 72 6d 65 64 2e |...p.G.e|rformed.|
|00005d80| 91 04 71 c7 54 68 69 73 | 96 03 55 55 73 68 6f 75 |..q.This|..UUshou|
|00005d90| 6c 64 93 73 70 98 65 65 | 64 93 75 70 93 74 68 65 |ld.sp.ee|d.up.the|
|00005da0| 93 72 6f 75 74 69 6e 65 | 73 93 73 6f 6d 65 77 68 |.routine|s.somewh|
|00005db0| 61 74 93 69 66 93 6f 6e | 6c 79 93 74 68 65 93 65 |at.if.on|ly.the.e|
|00005dc0| 69 67 65 6e 90 b8 e3 76 | 91 ff 71 c7 61 6c 75 65 |igen...v|..q.alue|
|00005dd0| 73 93 61 72 65 93 6e 65 | 65 64 65 64 2e 8e 9f 13 |s.are.ne|eded....|
|00005de0| 44 95 8d 8d db 45 58 41 | 4d 50 4c 45 8e 8e a6 8d |D....EXA|MPLE....|
|00005df0| 91 0f 00 00 b2 43 6f 6d | 70 75 74 69 6e 67 96 03 |.....Com|puting..|
|00005e00| 55 55 72 65 61 6c 93 53 | 63 95 b8 e3 68 93 75 72 |UUreal.S|c...h.ur|
|00005e10| 96 03 55 55 64 65 63 6f | 6d 70 9a 47 1d 6f 73 69 |..UUdeco|mp.G.osi|
|00005e20| 74 69 6f 6e 93 6f 66 93 | 28 70 98 6f 73 73 69 62 |tion.of.|(p.ossib|
|00005e30| 6c 79 29 93 6e 6f 6e 73 | 79 6d 6d 65 74 72 69 63 |ly).nons|ymmetric|
|00005e40| 93 b5 41 b2 3a 8e a9 18 | 44 95 8d da 4d 41 54 91 |..A.:...|D...MAT.|
|00005e50| 0f bf f7 2a 41 2c 96 05 | 3f fd 2a 53 2c 93 2a 51 |...*A,..|?.*S,.*Q|
|00005e60| 2c 93 2a 58 5f 72 65 2c | 93 2a 58 5f 69 6d 3b 8e |,.*X_re,|.*X_im;.|
|00005e70| a1 8d 56 45 43 91 0f bf | f7 2a 65 76 61 6c 73 5f |..VEC...|.*evals_|
|00005e80| 72 65 2c 91 05 3f fd 2a | 65 76 61 6c 73 5f 69 6d |re,..?.*|evals_im|
|00005e90| 3b 8e a1 8d 91 0a 7f fa | 2e 2e 2e 2e 2e 2e 8e a1 |;.......|........|
|00005ea0| 8d 53 96 05 3f fd 3d 93 | 63 70 5f 6d 61 74 28 41 |.S..?.=.|cp_mat(A|
|00005eb0| 2c 4d 4e 55 4c 4c 29 3b | 8e a1 8d 51 96 05 3f fd |,MNULL);|...Q..?.|
|00005ec0| 3d 93 67 65 74 5f 6d 61 | 74 28 41 2d 3e 6d 2c 41 |=.get_ma|t(A->m,A|
|00005ed0| 2d 3e 6d 29 3b 8e a1 8d | 73 63 68 75 72 28 53 2c |->m);...|schur(S,|
|00005ee0| 51 29 3b 8e a1 8d 2f 2a | 96 05 3f fd 67 65 74 93 |Q);.../*|..?.get.|
|00005ef0| 65 69 67 65 6e 76 61 6c | 75 65 73 93 28 72 65 61 |eigenval|ues.(rea|
|00005f00| 6c 2c 93 69 6d 61 67 69 | 6e 61 72 79 93 70 61 72 |l,.imagi|nary.par|
|00005f10| 74 73 29 93 2a 2f 8e a1 | 8d 65 76 61 6c 73 5f 72 |ts).*/..|.evals_r|
|00005f20| 65 96 05 3f fd 3d 93 67 | 65 74 5f 76 65 63 28 41 |e..?.=.g|et_vec(A|
|00005f30| 2d 3e 6d 29 3b 8e a1 8d | 65 76 61 6c 73 5f 69 6d |->m);...|evals_im|
|00005f40| 96 05 3f fd 3d 93 67 65 | 74 5f 76 65 63 28 41 2d |..?.=.ge|t_vec(A-|
|00005f50| 3e 6d 29 3b 8e a1 8d 73 | 63 68 75 72 5f 65 76 61 |>m);...s|chur_eva|
|00005f60| 6c 73 28 53 2c 65 76 61 | 6c 73 5f 72 65 2c 65 76 |ls(S,eva|ls_re,ev|
|00005f70| 61 6c 73 5f 69 6d 29 3b | 8e a1 8d 2f 2a 96 05 3f |als_im);|.../*..?|
|00005f80| fd 67 65 74 93 65 69 67 | 65 6e 76 65 63 74 6f 72 |.get.eig|envector|
|00005f90| 73 93 28 72 65 61 6c 2c | 93 69 6d 61 67 69 6e 61 |s.(real,|.imagina|
|00005fa0| 72 79 93 70 61 72 74 73 | 29 93 2a 2f 8e a1 8d 58 |ry.parts|).*/...X|
|00005fb0| 5f 72 65 96 05 3f fd 3d | 93 67 65 74 5f 6d 61 74 |_re..?.=|.get_mat|
|00005fc0| 28 41 2d 3e 6d 2c 41 2d | 3e 6d 29 3b 8e a1 8d 58 |(A->m,A-|>m);...X|
|00005fd0| 5f 69 6d 96 05 3f fd 3d | 93 67 65 74 5f 6d 61 74 |_im..?.=|.get_mat|
|00005fe0| 28 41 2d 3e 6d 2c 41 2d | 3e 6d 29 3b 8e a1 8d 73 |(A->m,A-|>m);...s|
|00005ff0| 63 68 75 72 5f 65 76 65 | 63 73 28 53 2c 51 2c 58 |chur_eve|cs(S,Q,X|
|00006000| 5f 72 65 2c 58 5f 69 6d | 29 3b 8e a6 8d 91 0f 00 |_re,X_im|);......|
|00006010| 00 b2 43 6f 6d 70 75 74 | 69 6e 67 96 03 55 55 65 |..Comput|ing..UUe|
|00006020| 69 67 65 6e 9a b8 e3 76 | 91 ff 71 c7 61 6c 75 65 |igen...v|..q.alue|
|00006030| 73 93 61 6e 64 93 65 69 | 67 65 6e 98 76 98 65 63 |s.and.ei|gen.v.ec|
|00006040| 74 6f 72 73 93 6f 66 93 | 61 93 72 65 61 6c 93 73 |tors.of.|a.real.s|
|00006050| 79 6d 6d 65 74 72 69 63 | 93 6d 61 74 72 69 78 3a |ymmetric|.matrix:|
|00006060| 8e a6 8d da 4d 41 54 91 | 0f bf f7 2a 41 2c 91 05 |....MAT.|...*A,..|
|00006070| 3f fd 2a 51 3b 8e a1 8d | 56 45 43 91 0f bf f7 2a |?.*Q;...|VEC....*|
|00006080| 65 76 61 6c 73 3b 8e a1 | 8d 91 0a 7f fa 2e 2e 2e |evals;..|........|
|00006090| 2e 2e 2e 8e a1 8d 65 76 | 61 6c 73 96 05 3f fd 3d |......ev|als..?.=|
|000060a0| 93 67 65 74 5f 76 65 63 | 28 41 2d 3e 6d 29 3b 8e |.get_vec|(A->m);.|
|000060b0| a1 8d 65 76 61 6c 73 96 | 05 3f fd 3d 93 73 79 6d |..evals.|.?.=.sym|
|000060c0| 6d 65 69 67 28 41 2c 51 | 2c 65 76 61 6c 73 29 3b |meig(A,Q|,evals);|
|000060d0| 8e a6 8d b2 54 68 65 96 | 03 55 55 b5 51 93 b2 6d |....The.|.UU.Q..m|
|000060e0| 61 74 72 69 78 93 63 6f | 6e 9a b8 e3 74 61 69 6e |atrix.co|n...tain|
|000060f0| 73 93 74 68 65 93 65 69 | 67 65 6e 98 76 98 65 63 |s.the.ei|gen.v.ec|
|00006100| 74 6f 72 73 2e 8e 9f 10 | 44 95 8d 91 0f 00 00 43 |tors....|D......C|
|00006110| 6f 6d 70 75 74 69 6e 67 | 96 03 72 99 74 68 65 93 |omputing|..r.the.|
|00006120| 65 69 67 65 6e 9a b8 e3 | 76 91 ff 71 c7 61 6c 75 |eigen...|v..q.alu|
|00006130| 65 73 93 61 6e 64 93 65 | 69 67 65 6e 98 76 98 65 |es.and.e|igen.v.e|
|00006140| 63 74 6f 72 73 93 6f 66 | 93 61 93 73 79 6d 6d 65 |ctors.of|.a.symme|
|00006150| 74 72 69 63 93 74 72 69 | 64 69 61 67 6f 6e 61 6c |tric.tri|diagonal|
|00006160| 93 6d 61 74 72 69 78 93 | 64 65 0c 6e 65 64 93 62 |.matrix.|de.ned.b|
|00006170| 98 79 93 74 68 65 93 76 | 98 65 63 74 6f 72 73 93 |.y.the.v|.ectors.|
|00006180| b5 61 8e a1 8d b2 28 74 | 68 65 96 03 55 55 64 69 |.a....(t|he..UUdi|
|00006190| 61 67 6f 6e 61 6c 93 65 | 6e 9a b8 e3 74 72 69 65 |agonal.e|n...trie|
|000061a0| 73 29 93 61 6e 64 93 b5 | 62 93 b2 28 74 68 65 93 |s).and..|b..(the.|
|000061b0| 6f 0b 2d 64 69 61 67 6f | 6e 61 6c 93 65 6e 98 74 |o.-diago|nal.en.t|
|000061c0| 72 69 65 73 29 3a 8e a6 | 8d da 4d 41 54 91 14 ff |ries):..|..MAT...|
|000061d0| f4 2a 51 3b 8e a1 8d 56 | 45 43 91 14 ff f4 2a 61 |.*Q;...V|EC....*a|
|000061e0| 2c 91 05 3f fd 2a 62 3b | 8e a1 8d 91 0a 7f fa 2e |,..?.*b;|........|
|000061f0| 2e 2e 2e 2e 2e 8e a1 8d | 51 96 05 3f fd 3d 93 67 |........|Q..?.=.g|
|00006200| 65 74 5f 6d 61 74 28 61 | 2d 3e 64 69 6d 2c 61 2d |et_mat(a|->dim,a-|
|00006210| 3e 64 69 6d 29 3b 8e a1 | 8d 69 64 5f 6d 61 74 28 |>dim);..|.id_mat(|
|00006220| 51 29 3b 91 14 ff f4 2f | 2a 96 05 3f fd 6d 75 73 |Q);..../|*..?.mus|
|00006230| 74 93 69 6e 69 74 69 61 | 6c 69 73 65 93 51 93 2a |t.initia|lise.Q.*|
|00006240| 2f 8e a1 8d 74 72 69 65 | 69 67 28 61 2c 62 2c 51 |/...trie|ig(a,b,Q|
|00006250| 29 3b 8e a1 8d 2f 2a 96 | 05 3f fd 61 93 69 73 93 |);.../*.|.?.a.is.|
|00006260| 6e 6f 77 93 74 68 65 93 | 76 65 63 74 6f 72 93 6f |now.the.|vector.o|
|00006270| 66 93 65 69 67 65 6e 76 | 61 6c 75 65 73 93 2a 2f |f.eigenv|alues.*/|
|00006280| 8e 9f 1b 44 95 8d 8d db | 53 45 45 91 03 55 55 41 |...D....|SEE..UUA|
|00006290| 4c 53 4f 8e 8e a4 10 44 | 95 8d 91 0f 00 00 b2 54 |LSO....D|.......T|
|000062a0| 68 65 96 03 55 55 48 65 | 73 73 65 6e 90 b8 e3 62 |he..UUHe|ssen...b|
|000062b0| 90 47 1d 65 72 67 93 72 | 6f 75 74 69 6e 65 73 93 |.G.erg.r|outines.|
|000062c0| 69 6e 93 da 68 65 73 73 | 65 6e 2e 63 b2 2e 8e a9 |in..hess|en.c....|
|000062d0| 13 44 95 8d 8d db 42 55 | 47 53 8e 8e a1 8d 91 0f |.D....BU|GS......|
|000062e0| 00 00 b2 49 74 96 03 9b | 22 69 73 93 75 70 93 74 |...It...|"is.up.t|
|000062f0| 6f 93 74 68 65 93 63 61 | 6c 6c 65 72 93 6f 66 93 |o.the.ca|ller.of.|
|00006300| da 73 79 6d 6d 65 69 67 | 28 29 93 b2 74 6f 93 65 |.symmeig|()..to.e|
|00006310| 6e 73 75 72 65 93 74 68 | 61 74 93 74 68 65 93 da |nsure.th|at.the..|
|00006320| 41 9b 03 9b 10 b2 6d 61 | 74 72 69 78 93 69 73 93 |A.....ma|trix.is.|
|00006330| 73 79 6d 6d 65 74 72 69 | 63 2e 91 05 43 2f 53 79 |symmetri|c...C/Sy|
|00006340| 6d 6d 65 74 72 79 93 6f | 66 93 da 41 98 b2 69 73 |mmetry.o|f..A..is|
|00006350| 93 6e 65 69 74 68 65 72 | 8e 9f 0c 00 00 8d 63 95 |.neither|......c.|
|00006360| b8 e3 68 65 63 93 6b 93 | 65 64 96 03 55 55 6e 6f |..hec.k.|ed..UUno|
|00006370| 72 93 65 6e 66 6f 72 63 | 65 64 93 69 6e 93 da 73 |r.enforc|ed.in..s|
|00006380| 79 6d 6d 65 69 67 28 29 | b2 2e 8e a6 8d 8d db 53 |ymmeig()|.......S|
|00006390| 4f 55 52 43 45 91 03 55 | 55 46 49 4c 45 3a 91 14 |OURCE..U|UFILE:..|
|000063a0| 00 02 da 73 79 6d 6d 65 | 69 67 2e 63 2c 91 05 3f |...symme|ig.c,..?|
|000063b0| fd 73 63 68 75 72 2e 63 | 8e 8e 8e 8e 8c 8b 00 00 |.schur.c|........|
|000063c0| 00 5a 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.Z......|........|
|000063d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000063e0| 00 00 00 00 00 00 00 00 | 8e 0d a0 02 a9 52 b3 8d |........|.....R..|
|000063f0| a0 fd 62 ad 4d 8d 8d b2 | 39 30 92 00 99 14 c1 de |..b.M...|90......|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
