home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
view JSON data
|
view as text
This file was processed as: ArcFS Archive
(archive/arcFS).
You can browse this item here: meschach
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| ArcFS Archive (archive/arcFS)
| magic
| Supported |
| 100%
| file
| RISC OS archive (ArcFS format)
| default
| |
| 99%
| file
| data
| default
| |
| 100%
| TrID
| Acorn ArcFS Archive
| default
| |
| 100%
| gt2
| Kopftext: 'Archive'
| default (weak)
| |
| 100%
| detectItEasy
| Archive: ArcFS
| default
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 41 72 63 68 69 76 65 00 | 94 11 00 00 34 14 00 00 |Archive.|....4...|
|00000010| 28 00 00 00 64 00 00 00 | 00 00 00 00 00 00 00 00 |(...d...|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000050| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000060| 82 21 4d 65 73 63 68 61 | 63 68 00 00 ff ff ff ff |.!Mescha|ch......|
|00000070| 45 ff ff ff 77 ef f4 a0 | 3b 00 00 00 ff ff ff ff |E...w...|;.......|
|00000080| 54 0f 00 80 82 21 42 6f | 6f 74 00 00 08 9b db 8c |T....!Bo|ot......|
|00000090| c0 00 00 00 45 eb ff ff | 2e 5c 08 80 3b 00 58 66 |....E...|.\..;.Xf|
|000000a0| c0 00 00 00 24 00 00 00 | 82 21 48 65 6c 70 00 00 |....$...|.!Help..|
|000000b0| 08 9b db 8c 06 05 00 00 | 45 ff ff ff f9 61 af 81 |........|E....a..|
|000000c0| 3b 00 61 36 06 05 00 00 | 08 01 00 00 82 21 52 75 |;.a6....|.....!Ru|
|000000d0| 6e 00 85 00 08 9b db 8c | 7c 00 00 00 45 eb ff ff |n.......||...E...|
|000000e0| f0 7a 08 80 3b 00 5c f3 | 7c 00 00 00 32 06 00 00 |.z..;.\.||...2...|
|000000f0| 82 21 53 70 72 69 74 65 | 73 00 db 8c e0 02 00 00 |.!Sprite|s.......|
|00000100| 45 f9 ff ff 0b bc 0a 80 | 3b 00 95 7d e0 02 00 00 |E.......|;..}....|
|00000110| d2 06 00 00 82 21 53 70 | 72 69 74 65 73 32 32 00 |.....!Sp|rites22.|
|00000120| 88 05 00 00 45 f9 ff ff | 84 20 0a 80 3b 00 8c 55 |....E...|. ..;..U|
|00000130| 88 05 00 00 d6 09 00 00 | 82 21 53 70 72 69 74 65 |........|.!Sprite|
|00000140| 73 32 33 00 58 02 00 00 | 45 f9 ff ff 70 62 0a 80 |s23.X...|E...pb..|
|00000150| 3b 00 2b 18 58 02 00 00 | 82 0f 00 00 82 63 00 86 |;.+.X...|.....c..|
|00000160| 0e 0f 85 00 08 9b db 8c | ff ff ff ff 45 ff ff ff |........|....E...|
|00000170| a7 f2 f4 a0 3b 00 00 00 | ff ff ff ff ac 0b 00 80 |....;...|........|
|00000180| 82 61 72 6e 6f 6c 64 69 | 00 00 00 00 ec 12 00 00 |.arnoldi|........|
|00000190| 45 ff ff ff 78 ef 40 17 | 3b 00 7d cf ec 12 00 00 |E...x.@.|;.}.....|
|000001a0| fe 11 00 00 82 62 64 66 | 61 63 74 6f 72 00 00 00 |.....bdf|actor...|
|000001b0| 44 34 00 00 45 ff ff ff | d8 1c 7d 57 3b 00 be 45 |D4..E...|..}W;..E|
|000001c0| 44 34 00 00 0e 25 00 00 | 82 62 6b 70 66 61 63 74 |D4...%..|.bkpfact|
|000001d0| 6f 00 00 00 02 21 00 00 | 45 ff ff ff 6c f1 40 17 |o....!..|E...l.@.|
|000001e0| 3b 00 85 7d 02 21 00 00 | 76 59 00 00 82 63 68 66 |;..}.!..|vY...chf|
|000001f0| 61 63 74 6f 72 00 00 00 | 9e 13 00 00 45 ff ff ff |actor...|....E...|
|00000200| 94 1a 40 17 3b 00 8f 93 | 9e 13 00 00 9c 7a 00 00 |..@.;...|.....z..|
|00000210| 82 63 6f 6e 6a 67 72 61 | 64 00 00 00 3a 20 00 00 |.conjgra|d...: ..|
|00000220| 45 ff ff ff 44 1f 40 17 | 3b 00 de 4c 3a 20 00 00 |E...D.@.|;..L: ..|
|00000230| 5e 8e 00 00 82 63 6f 70 | 79 00 00 00 00 00 00 00 |^....cop|y.......|
|00000240| 9a 15 00 00 45 ff ff ff | 78 27 40 17 3b 00 85 e4 |....E...|x'@.;...|
|00000250| 9a 15 00 00 bc ae 00 00 | 82 64 6d 61 63 68 65 70 |........|.dmachep|
|00000260| 73 00 00 00 7a 05 00 00 | 45 ff ff ff 1c 2e 40 17 |s...z...|E.....@.|
|00000270| 3b 00 b4 09 7a 05 00 00 | 7a c4 00 00 82 65 72 72 |;...z...|z....err|
|00000280| 00 00 00 00 00 00 00 00 | 1c 26 00 00 45 ff ff ff |........|.&..E...|
|00000290| 0c 39 40 17 3b 00 00 69 | 1c 26 00 00 18 ca 00 00 |.9@.;..i|.&......|
|000002a0| 82 65 78 74 72 61 73 00 | 00 00 00 00 af 2a 00 00 |.extras.|.....*..|
|000002b0| 45 ff ff ff 68 58 01 58 | 3b 00 a9 d8 af 2a 00 00 |E...hX.X|;....*..|
|000002c0| 58 f0 00 00 82 66 66 74 | 00 00 00 00 00 00 00 00 |X....fft|........|
|000002d0| 34 0f 00 00 45 ff ff ff | 00 ea 40 17 3b 00 3e 9f |4...E...|..@.;.>.|
|000002e0| 34 0f 00 00 2b 1b 01 00 | 82 66 6d 61 63 68 65 70 |4...+...|.fmachep|
|000002f0| 73 00 00 00 78 05 00 00 | 45 ff ff ff b8 5f 40 17 |s...x...|E...._@.|
|00000300| 3b 00 8f c4 78 05 00 00 | 83 2a 01 00 82 67 69 76 |;...x...|.*...giv|
|00000310| 65 6e 73 00 00 00 00 00 | 13 0e 00 00 45 ff ff ff |ens.....|....E...|
|00000320| ac 61 40 17 3b 00 b0 a4 | 13 0e 00 00 1f 30 01 00 |.a@.;...|.....0..|
|00000330| 82 68 65 73 73 65 6e 00 | 00 00 00 00 df 0f 00 00 |.hessen.|........|
|00000340| 45 ff ff ff 54 14 40 17 | 3b 00 50 87 df 0f 00 00 |E...T.@.|;.P.....|
|00000350| 56 3e 01 00 82 68 73 65 | 68 6c 64 72 00 00 00 00 |V>...hse|hldr....|
|00000360| 7f 12 00 00 45 ff ff ff | e4 15 40 17 3b 00 ef b2 |....E...|..@.;...|
|00000370| 7f 12 00 00 59 4e 01 00 | 82 69 6e 69 74 00 00 00 |....YN..|.init...|
|00000380| 00 00 00 00 01 17 00 00 | 45 ff ff ff e8 25 40 17 |........|E....%@.|
|00000390| 3b 00 d0 eb 01 17 00 00 | fc 60 01 00 82 69 6f 74 |;.......|.`...iot|
|000003a0| 6f 72 74 00 00 00 00 00 | 8c 0d 00 00 45 ff ff ff |ort.....|....E...|
|000003b0| 45 41 f7 7f 3b 00 83 90 | 8c 0d 00 00 21 78 01 00 |EA..;...|....!x..|
|000003c0| 82 69 74 65 72 30 00 00 | 00 00 00 00 2f 23 00 00 |.iter0..|..../#..|
|000003d0| 45 ff ff ff 23 1b 87 7f | 3b 00 75 fc 2f 23 00 00 |E...#...|;.u./#..|
|000003e0| d1 85 01 00 82 69 74 65 | 72 6e 73 79 6d 00 00 00 |.....ite|rnsym...|
|000003f0| 98 75 00 00 45 ff ff ff | d0 39 8f 21 3b 00 99 22 |.u..E...|.9.!;.."|
|00000400| 98 75 00 00 24 a9 01 00 | 82 69 74 65 72 73 79 6d |.u..$...|.itersym|
|00000410| 00 00 00 00 02 36 00 00 | 45 ff ff ff 38 53 40 17 |.....6..|E...8S@.|
|00000420| 3b 00 56 b4 02 36 00 00 | e0 1e 02 00 82 69 74 65 |;.V..6..|.....ite|
|00000430| 72 74 6f 72 74 00 00 00 | 32 3f 00 00 45 ff ff ff |rtort...|2?..E...|
|00000440| b4 63 a9 17 3b 00 0b 40 | 32 3f 00 00 06 55 02 00 |.c..;..@|2?...U..|
|00000450| 82 69 76 65 63 6f 70 00 | 00 00 00 00 48 26 00 00 |.ivecop.|....H&..|
|00000460| 45 ff ff ff d4 e8 40 17 | 3b 00 f1 6a 48 26 00 00 |E.....@.|;..jH&..|
|00000470| 5c 94 02 00 82 6c 61 6e | 63 7a 6f 73 00 00 00 00 |\....lan|czos....|
|00000480| 35 1e 00 00 45 ff ff ff | e8 5d 3f 17 3b 00 60 7a |5...E...|.]?.;.`z|
|00000490| 35 1e 00 00 c8 ba 02 00 | 82 6c 75 66 61 63 74 6f |5.......|.lufacto|
|000004a0| 72 00 00 00 c2 19 00 00 | 45 ff ff ff 98 c6 02 36 |r.......|E......6|
|000004b0| 3b 00 cf 89 c2 19 00 00 | 21 d9 02 00 82 6d 61 63 |;.......|!....mac|
|000004c0| 68 69 6e 65 00 00 00 00 | 2b 10 00 00 45 ff ff ff |hine....|+...E...|
|000004d0| 04 06 9d 7f 3b 00 d9 af | 2b 10 00 00 07 f3 02 00 |....;...|+.......|
|000004e0| 82 6d 61 74 6c 61 62 00 | 00 00 00 00 c7 13 00 00 |.matlab.|........|
|000004f0| 45 ff ff ff f8 84 3f 17 | 3b 00 5d 8b c7 13 00 00 |E.....?.|;.].....|
|00000500| 56 03 03 00 82 6d 61 74 | 6f 70 00 00 00 00 00 00 |V....mat|op......|
|00000510| ee 2e 00 00 45 ff ff ff | 64 8c 3f 17 3b 00 83 2e |....E...|d.?.;...|
|00000520| ee 2e 00 00 41 17 03 00 | 82 6d 61 74 72 69 78 69 |....A...|.matrixi|
|00000530| 6f 00 00 00 61 34 00 00 | 45 ff ff ff 74 9a 3f 17 |o...a4..|E...t.?.|
|00000540| 3b 00 81 6a 61 34 00 00 | 53 46 03 00 82 6d 61 78 |;..ja4..|SF...max|
|00000550| 69 6e 74 00 00 00 00 00 | e9 04 00 00 45 ff ff ff |int.....|....E...|
|00000560| 50 a0 3f 17 3b 00 f6 73 | e9 04 00 00 d8 7a 03 00 |P.?.;..s|.....z..|
|00000570| 82 6d 65 6d 69 6e 66 6f | 00 00 00 00 c3 23 00 00 |.meminfo|.....#..|
|00000580| 45 ff ff ff a0 b4 3f 17 | 3b 00 16 bf c3 23 00 00 |E.....?.|;....#..|
|00000590| e5 7f 03 00 82 6d 65 6d | 6f 72 79 00 00 00 00 00 |.....mem|ory.....|
|000005a0| ce 4d 00 00 45 ff ff ff | 30 ff 67 41 3b 00 55 9a |.M..E...|0.gA;.U.|
|000005b0| ce 4d 00 00 cc a3 03 00 | 82 6d 65 6d 73 74 61 74 |.M......|.memstat|
|000005c0| 00 00 00 00 c4 21 00 00 | 45 ff ff ff 54 c9 3f 17 |.....!..|E...T.?.|
|000005d0| 3b 00 11 33 c4 21 00 00 | be f1 03 00 82 6d 65 6d |;..3.!..|.....mem|
|000005e0| 74 6f 72 74 00 00 00 00 | c1 43 00 00 45 ff ff ff |tort....|.C..E...|
|000005f0| b4 63 a9 17 3b 00 fa 75 | c1 43 00 00 a6 13 04 00 |.c..;..u|.C......|
|00000600| 82 6d 66 75 6e 63 00 00 | 00 00 00 00 3a 23 00 00 |.mfunc..|....:#..|
|00000610| 45 ff ff ff fc df 3f 17 | 3b 00 5e 1b 3a 23 00 00 |E.....?.|;.^.:#..|
|00000620| 8b 57 04 00 82 6d 66 75 | 6e 74 6f 72 74 00 00 00 |.W...mfu|ntort...|
|00000630| b5 11 00 00 45 ff ff ff | b0 60 ab 17 3b 00 15 f8 |....E...|.`..;...|
|00000640| b5 11 00 00 e9 7a 04 00 | 82 6e 6f 72 6d 00 00 00 |.....z..|.norm...|
|00000650| 00 00 00 00 5b 10 00 00 | 45 ff ff ff c0 e9 3f 17 |....[...|E.....?.|
|00000660| 3b 00 10 e9 5b 10 00 00 | c2 8c 04 00 82 6f 74 68 |;...[...|.....oth|
|00000670| 65 72 69 6f 00 00 00 00 | 82 10 00 00 45 ff ff ff |erio....|....E...|
|00000680| 90 f1 3f 17 3b 00 f1 e3 | 82 10 00 00 41 9d 04 00 |..?.;...|....A...|
|00000690| 82 70 78 6f 70 00 00 00 | 00 00 00 00 49 1d 00 00 |.pxop...|....I...|
|000006a0| 45 ff ff ff e0 c6 2d 3b | 3b 00 e8 ba 49 1d 00 00 |E.....-;|;...I...|
|000006b0| e7 ad 04 00 82 71 72 66 | 61 63 74 6f 72 00 00 00 |.....qrf|actor...|
|000006c0| 45 34 00 00 45 ff ff ff | d8 fe 3f 17 3b 00 b4 21 |E4..E...|..?.;..!|
|000006d0| 45 34 00 00 54 cb 04 00 | 82 73 63 68 75 72 00 00 |E4..T...|.schur..|
|000006e0| 00 00 00 00 07 48 00 00 | 45 ff ff ff e0 e1 b1 37 |.....H..|E......7|
|000006f0| 3b 00 a0 00 07 48 00 00 | bd ff 04 00 82 73 6f 6c |;....H..|.....sol|
|00000700| 76 65 00 00 00 00 00 00 | 99 19 00 00 45 ff ff ff |ve......|....E...|
|00000710| e4 7f 3f 17 3b 00 91 43 | 99 19 00 00 e8 47 05 00 |..?.;..C|.....G..|
|00000720| 82 73 70 61 72 73 65 00 | 00 00 00 00 d5 5c 00 00 |.sparse.|.....\..|
|00000730| 45 ff ff ff fc 2f 10 33 | 3b 00 36 35 d5 5c 00 00 |E..../.3|;.65.\..|
|00000740| a5 61 05 00 82 73 70 61 | 72 73 65 69 6f 00 00 00 |.a...spa|rseio...|
|00000750| 21 20 00 00 45 ff ff ff | d8 e5 3f 17 3b 00 d1 cb |! ..E...|..?.;...|
|00000760| 21 20 00 00 9e be 05 00 | 82 73 70 62 6b 70 00 00 |! ......|.spbkp..|
|00000770| 00 00 00 00 14 8b 00 00 | 45 ff ff ff 4c d5 40 17 |........|E...L.@.|
|00000780| 3b 00 24 72 14 8b 00 00 | e3 de 05 00 82 73 70 63 |;.$r....|.....spc|
|00000790| 68 66 63 74 72 00 00 00 | 29 3e 00 00 45 ff ff ff |hfctr...|)>..E...|
|000007a0| 0c a3 3f 17 3b 00 2b f0 | 29 3e 00 00 1b 6a 06 00 |..?.;.+.|)>...j..|
|000007b0| 82 73 70 6c 75 66 63 74 | 72 00 00 00 38 2a 00 00 |.splufct|r...8*..|
|000007c0| 45 ff ff ff a4 e7 5d 53 | 3b 00 cb dd 38 2a 00 00 |E.....]S|;...8*..|
|000007d0| 68 a8 06 00 82 73 70 72 | 6f 77 00 00 00 00 00 00 |h....spr|ow......|
|000007e0| 03 45 00 00 45 ff ff ff | 7c 05 40 17 3b 00 cc 79 |.E..E...||.@.;..y|
|000007f0| 03 45 00 00 c4 d2 06 00 | 82 73 70 73 77 61 70 00 |.E......|.spswap.|
|00000800| 00 00 00 00 f0 1c 00 00 | 45 ff ff ff dc d6 40 17 |........|E.....@.|
|00000810| 3b 00 5e ad f0 1c 00 00 | eb 17 07 00 82 73 70 74 |;.^.....|.....spt|
|00000820| 6f 72 74 00 00 00 00 00 | 16 2c 00 00 45 ff ff ff |ort.....|.,..E...|
|00000830| 0c 23 52 2f 3b 00 39 92 | 16 2c 00 00 ff 34 07 00 |.#R/;.9.|.,...4..|
|00000840| 82 73 75 62 6d 61 74 00 | 00 00 00 00 b5 11 00 00 |.submat.|........|
|00000850| 45 ff ff ff e0 56 3f 17 | 3b 00 f7 1c b5 11 00 00 |E....V?.|;.......|
|00000860| 39 61 07 00 82 73 76 64 | 00 00 00 00 00 00 00 00 |9a...svd|........|
|00000870| bd 26 00 00 45 ff ff ff | 2c d2 40 17 3b 00 96 ff |.&..E...|,.@.;...|
|00000880| bd 26 00 00 12 73 07 00 | 82 73 79 6d 6d 65 69 67 |.&...s..|.symmeig|
|00000890| 00 00 00 00 4f 17 00 00 | 45 ff ff ff fc 7a b6 28 |....O...|E....z.(|
|000008a0| 3b 00 34 bf 4f 17 00 00 | f3 99 07 00 82 74 6f 72 |;.4.O...|.....tor|
|000008b0| 74 75 72 65 00 00 00 00 | 25 6e 00 00 45 ff ff ff |ture....|%n..E...|
|000008c0| f3 38 f6 7f 3b 00 aa dc | 25 6e 00 00 66 b1 07 00 |.8..;...|%n..f...|
|000008d0| 82 74 75 74 61 64 76 00 | 00 00 00 00 9c 11 00 00 |.tutadv.|........|
|000008e0| 45 ff ff ff 40 f2 20 58 | 3b 00 1f c2 9c 11 00 00 |E...@. X|;.......|
|000008f0| af 1f 08 00 82 74 75 74 | 6f 72 69 61 6c 00 00 00 |.....tut|orial...|
|00000900| 57 1e 00 00 45 ff ff ff | 74 96 20 58 3b 00 86 40 |W...E...|t. X;..@|
|00000910| 57 1e 00 00 6f 31 08 00 | 82 75 70 64 61 74 65 00 |W...o1..|.update.|
|00000920| 00 00 00 00 71 0d 00 00 | 45 ff ff ff c8 28 3f 17 |....q...|E....(?.|
|00000930| 3b 00 4a 64 71 0d 00 00 | ea 4f 08 00 82 76 65 63 |;.Jdq...|.O...vec|
|00000940| 6f 70 00 00 00 00 00 00 | 76 34 00 00 45 ff ff ff |op......|v4..E...|
|00000950| 00 8b 10 33 3b 00 6c 1b | 76 34 00 00 7f 5d 08 00 |...3;.l.|v4...]..|
|00000960| 82 76 65 72 73 69 6f 6e | 00 00 00 00 02 0a 00 00 |.version|........|
|00000970| 45 ff ff ff 2c 7a 2e 3b | 3b 00 df 29 02 0a 00 00 |E...,z.;|;..)....|
|00000980| 19 92 08 00 82 7a 63 6f | 70 79 00 00 00 00 00 00 |.....zco|py......|
|00000990| 54 14 00 00 45 ff ff ff | f4 e3 39 17 3b 00 a7 74 |T...E...|..9.;..t|
|000009a0| 54 14 00 00 3f 9c 08 00 | 82 7a 66 75 6e 63 00 00 |T...?...|.zfunc..|
|000009b0| 00 00 00 00 a5 11 00 00 | 45 ff ff ff cb 66 85 7f |........|E....f..|
|000009c0| 3b 00 f6 0e a5 11 00 00 | b7 b0 08 00 82 7a 67 69 |;.......|.....zgi|
|000009d0| 76 65 6e 73 00 00 00 00 | c1 12 00 00 45 ff ff ff |vens....|....E...|
|000009e0| b8 43 09 33 3b 00 8d a9 | c1 12 00 00 80 c2 08 00 |.C.3;...|........|
|000009f0| 82 7a 68 65 73 73 65 6e | 00 00 00 00 6c 0f 00 00 |.zhessen|....l...|
|00000a00| 45 ff ff ff 1c bc 39 17 | 3b 00 b0 7b 6c 0f 00 00 |E.....9.|;..{l...|
|00000a10| 65 d5 08 00 82 7a 68 73 | 65 68 6c 64 72 00 00 00 |e....zhs|ehldr...|
|00000a20| 9c 15 00 00 45 ff ff ff | 4c 14 6d 42 3b 00 54 46 |....E...|L.mB;.TF|
|00000a30| 9c 15 00 00 f5 e4 08 00 | 82 7a 6c 75 66 63 74 72 |........|.zlufctr|
|00000a40| 00 00 00 00 ec 1a 00 00 | 45 ff ff ff 70 ae 39 17 |........|E...p.9.|
|00000a50| 3b 00 d9 7a ec 1a 00 00 | b5 fa 08 00 82 7a 6d 61 |;..z....|.....zma|
|00000a60| 63 68 69 6e 65 00 00 00 | 9f 10 00 00 45 ff ff ff |chine...|....E...|
|00000a70| 30 a8 39 17 3b 00 50 ff | 9f 10 00 00 c5 15 09 00 |0.9.;.P.|........|
|00000a80| 82 7a 6d 61 74 69 6f 00 | 00 00 00 00 99 29 00 00 |.zmatio.|.....)..|
|00000a90| 45 ff ff ff f4 98 39 17 | 3b 00 7b 1a 99 29 00 00 |E.....9.|;.{..)..|
|00000aa0| 88 26 09 00 82 7a 6d 61 | 74 6c 61 62 00 00 00 00 |.&...zma|tlab....|
|00000ab0| d5 16 00 00 45 ff ff ff | f8 8f 39 17 3b 00 6d 5c |....E...|..9.;.m\|
|00000ac0| d5 16 00 00 45 50 09 00 | 82 7a 6d 61 74 6f 70 00 |....EP..|.zmatop.|
|00000ad0| 00 00 00 00 95 3b 00 00 | 45 ff ff ff c4 87 39 17 |.....;..|E.....9.|
|00000ae0| 3b 00 0e 0b 95 3b 00 00 | 3e 67 09 00 82 7a 6d 65 |;....;..|>g...zme|
|00000af0| 6d 6f 72 79 00 00 00 00 | ab 3b 00 00 45 ff ff ff |mory....|.;..E...|
|00000b00| 60 10 68 41 3b 00 a8 b0 | ab 3b 00 00 f7 a2 09 00 |`.hA;...|.;......|
|00000b10| 82 7a 6e 6f 72 6d 00 00 | 00 00 00 00 10 12 00 00 |.znorm..|........|
|00000b20| 45 ff ff ff ac 40 39 17 | 3b 00 7b 0b 10 12 00 00 |E....@9.|;.{.....|
|00000b30| c6 de 09 00 82 7a 71 72 | 66 63 74 72 00 00 00 00 |.....zqr|fctr....|
|00000b40| dc 35 00 00 45 ff ff ff | ac de f8 7f 3b 00 ff 3d |.5..E...|....;..=|
|00000b50| dc 35 00 00 fa f0 09 00 | 82 7a 73 63 68 75 72 00 |.5......|.zschur.|
|00000b60| 00 00 00 00 e1 2a 00 00 | 45 ff ff ff c8 33 39 17 |.....*..|E....39.|
|00000b70| 3b 00 38 37 e1 2a 00 00 | fa 26 0a 00 82 7a 73 6f |;.87.*..|.&...zso|
|00000b80| 6c 76 65 00 00 00 00 00 | 95 1d 00 00 45 ff ff ff |lve.....|....E...|
|00000b90| 80 26 39 17 3b 00 d0 ee | 95 1d 00 00 ff 51 0a 00 |.&9.;...|.....Q..|
|00000ba0| 82 7a 74 6f 72 74 75 72 | 65 00 00 00 8b 4e 00 00 |.ztortur|e....N..|
|00000bb0| 45 ff ff ff ae 50 f9 7f | 3b 00 04 52 8b 4e 00 00 |E....P..|;..R.N..|
|00000bc0| b8 6f 0a 00 82 7a 76 65 | 63 6f 70 00 00 00 00 00 |.o...zve|cop.....|
|00000bd0| d9 2b 00 00 45 ff ff ff | 28 95 10 33 3b 00 bc 31 |.+..E...|(..3;..1|
|00000be0| d9 2b 00 00 67 be 0a 00 | 00 00 00 00 00 00 00 00 |.+..g...|........|
|00000bf0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000c00| 00 00 00 00 00 00 00 00 | 00 00 00 00 82 63 6f 70 |........|.....cop|
|00000c10| 79 72 69 67 68 74 00 8c | 64 04 00 00 45 ff ff ff |yright..|d...E...|
|00000c20| 60 bb 41 17 3b 00 ab 51 | 64 04 00 00 64 ea 0a 00 |`.A.;..Q|d...d...|
|00000c30| 82 44 4f 43 00 0f 85 00 | 08 9b db 8c ff ff ff ff |.DOC....|........|
|00000c40| 45 ff ff ff 95 21 f5 a0 | 3b 00 00 00 ff ff ff ff |E....!..|;.......|
|00000c50| 60 0c 00 80 82 66 6e 69 | 6e 64 65 78 00 00 00 00 |`....fni|ndex....|
|00000c60| 22 43 00 00 45 ff ff ff | 9c c5 aa 17 3b 00 72 59 |"C..E...|....;.rY|
|00000c70| 22 43 00 00 ec ee 0a 00 | 82 74 75 74 6f 72 69 61 |"C......|.tutoria|
|00000c80| 6c 00 00 00 9c b3 00 00 | 45 ff ff ff 00 c6 aa 17 |l.......|E.......|
|00000c90| 3b 00 7c 4e 9c b3 00 00 | 32 32 0b 00 00 00 00 00 |;.|N....|22......|
|00000ca0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000cb0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000cc0| 82 68 00 86 0e 0f 85 00 | 08 9b db 8c ff ff ff ff |.h......|........|
|00000cd0| 45 ff ff ff 0f 24 f5 a0 | 3b 00 00 00 ff ff ff ff |E....$..|;.......|
|00000ce0| 7c 0e 00 80 82 65 72 72 | 00 00 00 00 00 00 00 00 ||....err|........|
|00000cf0| 2e 16 00 00 45 ff ff ff | 90 3c 40 17 3b 00 2d 6e |....E...|.<@.;.-n|
|00000d00| 2e 16 00 00 f2 e5 0b 00 | 82 69 74 65 72 00 00 00 |........|.iter...|
|00000d10| 00 00 00 00 ff 1a 00 00 | 45 ff ff ff b0 5d 10 33 |........|E....].3|
|00000d20| 3b 00 72 0c ff 1a 00 00 | 44 fc 0b 00 82 6d 61 63 |;.r.....|D....mac|
|00000d30| 68 69 6e 65 00 00 00 00 | 04 0f 00 00 45 ff ff ff |hine....|....E...|
|00000d40| b8 56 8b 7f 3b 00 72 f7 | 04 0f 00 00 67 17 0c 00 |.V..;.r.|....g...|
|00000d50| 82 6d 61 74 6c 61 62 00 | 00 00 00 00 b6 0b 00 00 |.matlab.|........|
|00000d60| 45 ff ff ff 44 e7 40 17 | 3b 00 3b b4 b6 0b 00 00 |E...D.@.|;.;.....|
|00000d70| 8f 26 0c 00 82 6d 61 74 | 72 69 78 00 00 00 00 00 |.&...mat|rix.....|
|00000d80| 20 4b 00 00 45 ff ff ff | 44 3a 09 47 3b 00 dd ad | K..E...|D:.G;...|
|00000d90| 20 4b 00 00 69 32 0c 00 | 82 6d 61 74 72 69 78 32 | K..i2..|.matrix2|
|00000da0| 00 00 00 00 63 20 00 00 | 45 ff ff ff b8 97 3f 17 |....c ..|E.....?.|
|00000db0| 3b 00 aa 2a 63 20 00 00 | ad 7d 0c 00 82 6d 65 6d |;..*c ..|.}...mem|
|00000dc0| 69 6e 66 6f 00 00 00 00 | 34 10 00 00 45 ff ff ff |info....|4...E...|
|00000dd0| c0 b7 3f 17 3b 00 4f d4 | 34 10 00 00 34 9e 0c 00 |..?.;.O.|4...4...|
|00000de0| 82 6f 6c 64 6e 61 6d 65 | 73 00 00 00 0d 0f 00 00 |.oldname|s.......|
|00000df0| 45 ff ff ff 70 ee 3f 17 | 3b 00 03 03 0d 0f 00 00 |E...p.?.|;.......|
|00000e00| 8c ae 0c 00 82 4f 6d 61 | 63 68 69 6e 65 00 00 00 |.....Oma|chine...|
|00000e10| ae 11 00 00 45 ff ff ff | 5c 70 0b 33 3b 00 8f 52 |....E...|\p.3;..R|
|00000e20| ae 11 00 00 bd bd 0c 00 | 82 73 70 61 72 73 65 00 |........|.sparse.|
|00000e30| 00 00 00 00 53 19 00 00 | 45 ff ff ff b4 d2 3f 17 |....S...|E.....?.|
|00000e40| 3b 00 38 58 53 19 00 00 | 8f cf 0c 00 82 73 70 61 |;.8XS...|.....spa|
|00000e50| 72 73 65 32 00 00 00 00 | 58 0c 00 00 45 ff ff ff |rse2....|X...E...|
|00000e60| 58 d9 3f 17 3b 00 bf 1b | 58 0c 00 00 06 e9 0c 00 |X.?.;...|X.......|
|00000e70| 82 7a 6d 61 74 72 69 78 | 00 00 00 00 4e 22 00 00 |.zmatrix|....N"..|
|00000e80| 45 ff ff ff 48 7f 10 33 | 3b 00 8b 57 4e 22 00 00 |E...H..3|;..WN"..|
|00000e90| 82 f5 0c 00 82 7a 6d 61 | 74 72 69 78 32 00 00 00 |.....zma|trix2...|
|00000ea0| 37 10 00 00 45 ff ff ff | 64 7e 39 17 3b 00 2f 28 |7...E...|d~9.;./(|
|00000eb0| 37 10 00 00 f4 17 0d 00 | 00 00 00 00 00 00 00 00 |7.......|........|
|00000ec0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000ed0| 00 00 00 00 00 00 00 00 | 00 00 00 00 82 4d 61 6b |........|.....Mak|
|00000ee0| 65 4c 69 62 00 9b db 8c | a2 47 00 00 45 e1 ff ff |eLib....|.G..E...|
|00000ef0| ad 1e f9 7f 3b 00 5e a2 | a2 47 00 00 4f 28 0d 00 |....;.^.|.G..O(..|
|00000f00| 82 6f 00 86 0e 0f 85 00 | 08 9b db 8c ff ff ff ff |.o......|........|
|00000f10| 45 ff ff ff 5d 2c f5 a0 | 3b 00 00 00 ff ff ff ff |E...],..|;.......|
|00000f20| 0c 0f 00 80 82 6d 65 73 | 63 68 61 63 68 00 00 00 |.....mes|chach...|
|00000f30| b8 e0 04 00 45 fd ff ff | c4 e4 f8 7f 3b 00 aa a9 |....E...|....;...|
|00000f40| b8 e0 04 00 15 70 0d 00 | 00 00 00 00 00 00 00 00 |.....p..|........|
|00000f50| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000f60| 00 00 00 00 00 00 00 00 | 00 00 00 00 82 52 65 61 |........|.....Rea|
|00000f70| 64 4d 65 00 08 9b db 8c | 56 46 00 00 45 ff ff ff |dMe.....|VF..E...|
|00000f80| 60 ac 89 7e 3b 00 90 f9 | 56 46 00 00 f1 50 12 00 |`..~;...|VF...P..|
|00000f90| 00 de d9 86 0e 0f 85 00 | 08 9b db 8c 0d 0c 8c 85 |........|........|
|00000fa0| 8c de d9 0c 82 98 db 8d | dd 8c 40 86 8c 93 da 8d |........|..@.....|
|00000fb0| 89 aa 02 09 82 42 6f 6f | 6b 00 00 00 00 00 00 00 |.....Boo|k.......|
|00000fc0| ff ff ff ff 45 ff ff ff | 00 30 f5 a0 3b 00 00 00 |....E...|.0..;...|
|00000fd0| ff ff ff ff 4c 11 00 80 | 82 63 74 73 00 b3 00 00 |....L...|.cts....|
|00000fe0| 0a 0a 09 09 08 75 00 00 | 44 e4 fc ff 40 da ed c6 |.....u..|D...@...|
|00000ff0| 3b 00 ae 4c 08 75 00 00 | 6b 97 12 00 82 64 61 74 |;..L.u..|k....dat|
|00001000| 61 73 74 72 75 63 74 00 | 48 78 00 00 44 e4 fc ff |astruct.|Hx..D...|
|00001010| 84 73 ed c6 3b 00 85 60 | 48 78 00 00 97 0c 13 00 |.s..;..`|Hx......|
|00001020| 82 64 65 73 69 67 6e 6c | 69 62 00 09 dc 2e 01 00 |.designl|ib......|
|00001030| 44 e4 fc ff 4c 7c 95 f3 | 3b 00 75 c4 dc 2e 01 00 |D...L|..|;.u.....|
|00001040| 03 85 13 00 82 66 6e 69 | 6e 64 65 78 00 0a 09 09 |.....fni|ndex....|
|00001050| 34 79 00 00 44 e4 fc ff | 00 d4 ed c6 3b 00 44 6c |4y..D...|....;.Dl|
|00001060| 34 79 00 00 03 b4 14 00 | 82 69 6e 64 65 78 00 00 |4y......|.index..|
|00001070| 0a 0a 09 09 7c 30 00 00 | 44 e4 fc ff 84 d7 ed c6 |....|0..|D.......|
|00001080| 3b 00 5d 24 7c 30 00 00 | 5b 2d 15 00 82 69 6e 73 |;.]$|0..|[-...ins|
|00001090| 74 61 6c 6c 00 0a 09 09 | ec 48 00 00 44 e4 fc ff |tall....|.H..D...|
|000010a0| 90 bc ed c6 3b 00 40 04 | ec 48 00 00 fb 5d 15 00 |....;.@.|.H...]..|
|000010b0| 82 6d 61 6e 75 61 6c 00 | 0a 0a 09 09 70 1c 01 00 |.manual.|....p...|
|000010c0| 44 e4 fc ff ac 96 ed c6 | 3b 00 fc 25 70 1c 01 00 |D.......|;..%p...|
|000010d0| 0b a7 15 00 82 6d 61 6e | 75 61 6c 32 00 0a 09 09 |.....man|ual2....|
|000010e0| 10 b2 00 00 44 e4 fc ff | e8 a5 ed c6 3b 00 19 9f |....D...|....;...|
|000010f0| 10 b2 00 00 9f c3 16 00 | 82 6d 61 6e 75 61 6c 33 |........|.manual3|
|00001100| 00 0a 09 09 b0 bb 00 00 | 44 e4 fc ff 88 b5 ed c6 |........|D.......|
|00001110| 3b 00 4a f5 b0 bb 00 00 | d3 75 17 00 82 6e 75 6d |;.J.....|.u...num|
|00001120| 6c 69 6e 61 6c 67 00 09 | f0 a6 00 00 44 e4 fc ff |linalg..|....D...|
|00001130| 30 81 ed c6 3b 00 bc 9d | f0 a6 00 00 a7 31 18 00 |0...;...|.....1..|
|00001140| 82 72 65 61 64 69 6e 67 | 00 0a 09 09 30 0a 00 00 |.reading|....0...|
|00001150| 44 e4 fc ff 18 d0 ed c6 | 3b 00 ae 05 30 0a 00 00 |D.......|;...0...|
|00001160| bb d8 18 00 82 74 75 74 | 6f 72 69 61 6c 00 09 09 |.....tut|orial...|
|00001170| e8 bc 00 00 44 e4 fc ff | 30 68 ed c6 3b 00 14 e3 |....D...|0h..;...|
|00001180| e8 bc 00 00 0f e3 18 00 | 00 74 20 28 e0 b3 00 00 |........|.t (....|
|00001190| 0a 0a 09 09 09 20 20 20 | 4d 45 53 43 48 41 43 48 |..... |MESCHACH|
|000011a0| 20 56 45 52 53 49 4f 4e | bc 1f 2e 32 82 43 41 55 | VERSION|...2.CAU|
|000011b0| 47 00 00 00 00 00 00 00 | e3 1a 00 00 45 ff ff ff |G.......|....E...|
|000011c0| ad e9 af 81 3b 00 85 be | e3 1a 00 00 1b a0 19 00 |....;...|........|
|000011d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011f0| 00 00 00 00 68 65 20 4d | 65 73 63 68 61 63 68 3a |....he M|eschach:|
|00001200| 20 70 61 74 68 0a 0a 53 | 65 74 20 4d 65 73 63 68 | path..S|et Mesch|
|00001210| 61 63 68 24 44 69 72 20 | 3c 4f 62 65 79 24 44 69 |ach$Dir |<Obey$Di|
|00001220| 72 3e 0a 53 65 74 20 4d | 65 73 63 68 61 63 68 24 |r>.Set M|eschach$|
|00001230| 50 61 74 68 20 3c 4d 65 | 73 63 68 61 63 68 24 44 |Path <Me|schach$D|
|00001240| 69 72 3e 2e 0a 0a 0a 0a | 46 69 6c 65 44 61 74 61 |ir>.....|FileData|
|00001250| 82 21 48 65 6c 70 00 00 | 08 9b db 8c 06 05 00 00 |.!Help..|........|
|00001260| 45 ff ff ff f9 61 af 81 | 03 00 61 36 21 4d 65 73 |E....a..|..a6!Mes|
|00001270| 63 68 61 63 68 20 56 31 | 2e 32 62 0a 2d 2d 2d 2d |chach V1|.2b.----|
|00001280| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 0a 0a 43 6f 70 |--------|---..Cop|
|00001290| 79 72 69 67 68 74 20 28 | 43 29 20 31 39 39 33 20 |yright (|C) 1993 |
|000012a0| 44 61 76 69 64 20 45 2e | 20 53 74 65 77 61 72 64 |David E.| Steward|
|000012b0| 20 26 20 5a 62 69 67 6e | 69 65 77 20 4c 65 79 6b | & Zbign|iew Leyk|
|000012c0| 2c 20 61 6c 6c 20 72 69 | 67 68 74 73 20 72 65 73 |, all ri|ghts res|
|000012d0| 65 72 76 65 64 2e 0a 0a | 50 6f 72 74 20 66 6f 72 |erved...|Port for|
|000012e0| 20 41 63 6f 72 6e 20 52 | 49 53 43 20 4d 61 63 68 | Acorn R|ISC Mach|
|000012f0| 69 6e 65 73 20 62 79 20 | 50 61 75 6c 20 46 69 65 |ines by |Paul Fie|
|00001300| 6c 64 0a 0a 0a 4d 65 73 | 63 68 61 63 68 20 69 73 |ld...Mes|chach is|
|00001310| 20 61 20 6c 69 62 72 61 | 72 79 20 6f 66 20 43 20 | a libra|ry of C |
|00001320| 72 6f 75 74 69 6e 65 73 | 20 66 6f 72 20 70 65 72 |routines| for per|
|00001330| 66 6f 72 6d 69 6e 67 20 | 63 61 6c 63 75 6c 61 74 |forming |calculat|
|00001340| 69 6f 6e 20 6f 6e 20 6d | 61 74 72 69 63 65 73 20 |ion on m|atrices |
|00001350| 61 6e 64 20 76 65 63 74 | 6f 72 73 2e 20 46 6f 72 |and vect|ors. For|
|00001360| 20 6d 6f 72 65 20 69 6e | 66 6f 72 6d 61 74 69 6f | more in|formatio|
|00001370| 6e 20 72 65 61 64 20 74 | 68 65 20 52 65 61 64 4d |n read t|he ReadM|
|00001380| 65 20 66 69 6c 65 20 69 | 6e 20 74 68 65 20 4d 65 |e file i|n the Me|
|00001390| 73 63 68 61 63 68 20 64 | 69 72 65 63 74 6f 72 79 |schach d|irectory|
|000013a0| 20 61 6e 64 20 74 68 65 | 20 66 69 6c 65 73 20 69 | and the| files i|
|000013b0| 6e 20 74 68 65 20 44 4f | 43 20 64 69 72 65 63 74 |n the DO|C direct|
|000013c0| 6f 72 79 2e 0a 54 6f 20 | 75 73 65 20 74 68 65 20 |ory..To |use the |
|000013d0| 6c 69 62 72 61 72 79 20 | 6a 75 73 74 20 6c 69 6e |library |just lin|
|000013e0| 6b 20 77 69 74 68 20 27 | 4d 65 73 61 63 68 61 72 |k with '|Mesachar|
|000013f0| 63 68 3a 6f 2e 6d 65 73 | 63 68 61 72 63 68 27 2e |ch:o.mes|charch'.|
|00001400| 0a 0a 41 20 6d 61 6e 75 | 61 6c 20 66 6f 72 20 61 |..A manu|al for a|
|00001410| 6e 20 6f 6c 64 65 72 20 | 76 65 72 73 69 6f 6e 20 |n older |version |
|00001420| 6f 66 20 4d 65 73 63 68 | 61 63 68 20 68 61 73 20 |of Mesch|ach has |
|00001430| 62 65 65 6e 46 69 6c 65 | 44 61 74 61 82 21 42 6f |beenFile|Data.!Bo|
|00001440| 6f 74 00 00 08 9b db 8c | c0 00 00 00 45 eb ff ff |ot......|....E...|
|00001450| 2e 5c 08 80 03 00 58 66 | 49 63 6f 6e 53 70 72 69 |.\....Xf|IconSpri|
|00001460| 74 65 73 20 3c 4f 62 65 | 79 24 44 69 72 3e 2e 21 |tes <Obe|y$Dir>.!|
|00001470| 53 70 72 69 74 65 73 0a | 0a 49 46 20 22 3c 4d 65 |Sprites.|.IF "<Me|
|00001480| 73 63 68 61 63 68 24 44 | 69 72 3e 22 20 3d 20 22 |schach$D|ir>" = "|
|00001490| 22 20 54 48 45 4e 20 53 | 65 74 20 43 24 50 61 74 |" THEN S|et C$Pat|
|000014a0| 68 20 3c 43 24 50 61 74 | 68 3e 2e 2c 3c 4f 62 65 |h <C$Pat|h>.,<Obe|
|000014b0| 79 24 44 69 72 3e 2e 0a | 0a 0a 7c 20 53 65 74 20 |y$Dir>..|..| Set |
|000014c0| 75 70 20 74 68 65 20 4d | 65 73 63 68 61 63 68 3a |up the M|eschach:|
|000014d0| 20 70 61 74 68 0a 0a 53 | 65 74 20 4d 65 73 63 68 | path..S|et Mesch|
|000014e0| 61 63 68 24 44 69 72 20 | 3c 4f 62 65 79 24 44 69 |ach$Dir |<Obey$Di|
|000014f0| 72 3e 0a 53 65 74 20 4d | 65 73 63 68 61 63 68 24 |r>.Set M|eschach$|
|00001500| 50 61 74 68 20 3c 4d 65 | 73 63 68 61 63 68 24 44 |Path <Me|schach$D|
|00001510| 69 72 3e 2e 0a 0a 0a 0a | 46 69 6c 65 44 61 74 61 |ir>.....|FileData|
|00001520| 82 21 48 65 6c 70 00 00 | 08 9b db 8c 06 05 00 00 |.!Help..|........|
|00001530| 45 ff ff ff f9 61 af 81 | 03 00 61 36 21 4d 65 73 |E....a..|..a6!Mes|
|00001540| 63 68 61 63 68 20 56 31 | 2e 32 62 0a 2d 2d 2d 2d |chach V1|.2b.----|
|00001550| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 0a 0a 43 6f 70 |--------|---..Cop|
|00001560| 79 72 69 67 68 74 20 28 | 43 29 20 31 39 39 33 20 |yright (|C) 1993 |
|00001570| 44 61 76 69 64 20 45 2e | 20 53 74 65 77 61 72 64 |David E.| Steward|
|00001580| 20 26 20 5a 62 69 67 6e | 69 65 77 20 4c 65 79 6b | & Zbign|iew Leyk|
|00001590| 2c 20 61 6c 6c 20 72 69 | 67 68 74 73 20 72 65 73 |, all ri|ghts res|
|000015a0| 65 72 76 65 64 2e 0a 0a | 50 6f 72 74 20 66 6f 72 |erved...|Port for|
|000015b0| 20 41 63 6f 72 6e 20 52 | 49 53 43 20 4d 61 63 68 | Acorn R|ISC Mach|
|000015c0| 69 6e 65 73 20 62 79 20 | 50 61 75 6c 20 46 69 65 |ines by |Paul Fie|
|000015d0| 6c 64 0a 0a 0a 4d 65 73 | 63 68 61 63 68 20 69 73 |ld...Mes|chach is|
|000015e0| 20 61 20 6c 69 62 72 61 | 72 79 20 6f 66 20 43 20 | a libra|ry of C |
|000015f0| 72 6f 75 74 69 6e 65 73 | 20 66 6f 72 20 70 65 72 |routines| for per|
|00001600| 66 6f 72 6d 69 6e 67 20 | 63 61 6c 63 75 6c 61 74 |forming |calculat|
|00001610| 69 6f 6e 20 6f 6e 20 6d | 61 74 72 69 63 65 73 20 |ion on m|atrices |
|00001620| 61 6e 64 20 76 65 63 74 | 6f 72 73 2e 20 46 6f 72 |and vect|ors. For|
|00001630| 20 6d 6f 72 65 20 69 6e | 66 6f 72 6d 61 74 69 6f | more in|formatio|
|00001640| 6e 20 72 65 61 64 20 74 | 68 65 20 52 65 61 64 4d |n read t|he ReadM|
|00001650| 65 20 66 69 6c 65 20 69 | 6e 20 74 68 65 20 4d 65 |e file i|n the Me|
|00001660| 73 63 68 61 63 68 20 64 | 69 72 65 63 74 6f 72 79 |schach d|irectory|
|00001670| 20 61 6e 64 20 74 68 65 | 20 66 69 6c 65 73 20 69 | and the| files i|
|00001680| 6e 20 74 68 65 20 44 4f | 43 20 64 69 72 65 63 74 |n the DO|C direct|
|00001690| 6f 72 79 2e 0a 54 6f 20 | 75 73 65 20 74 68 65 20 |ory..To |use the |
|000016a0| 6c 69 62 72 61 72 79 20 | 6a 75 73 74 20 6c 69 6e |library |just lin|
|000016b0| 6b 20 77 69 74 68 20 27 | 4d 65 73 61 63 68 61 72 |k with '|Mesachar|
|000016c0| 63 68 3a 6f 2e 6d 65 73 | 63 68 61 72 63 68 27 2e |ch:o.mes|charch'.|
|000016d0| 0a 0a 41 20 6d 61 6e 75 | 61 6c 20 66 6f 72 20 61 |..A manu|al for a|
|000016e0| 6e 20 6f 6c 64 65 72 20 | 76 65 72 73 69 6f 6e 20 |n older |version |
|000016f0| 6f 66 20 4d 65 73 63 68 | 61 63 68 20 68 61 73 20 |of Mesch|ach has |
|00001700| 62 65 65 6e 20 69 6e 63 | 6c 75 64 65 64 20 69 6e |been inc|luded in|
|00001710| 20 74 68 65 20 42 6f 6f | 6b 20 64 69 72 65 63 74 | the Boo|k direct|
|00001720| 6f 72 79 20 62 75 74 20 | 69 74 20 69 73 20 69 6e |ory but |it is in|
|00001730| 20 44 56 49 20 66 6f 72 | 6d 61 74 20 61 6e 64 20 | DVI for|mat and |
|00001740| 6e 65 65 64 73 20 61 20 | 70 72 6f 67 72 61 6d 20 |needs a |program |
|00001750| 73 75 63 68 20 61 73 20 | 21 44 56 49 76 69 65 77 |such as |!DVIview|
|00001760| 20 74 6f 20 64 69 73 70 | 6c 61 79 20 69 74 2e 20 | to disp|lay it. |
|00001770| 21 44 56 49 76 69 65 77 | 20 69 73 20 61 76 61 69 |!DVIview| is avai|
|00001780| 6c 61 62 6c 65 20 66 72 | 6f 6d 20 76 61 72 69 6f |lable fr|om vario|
|00001790| 75 73 20 66 74 70 20 73 | 69 74 65 73 20 6f 72 20 |us ftp s|ites or |
|000017a0| 66 72 6f 6d 20 6d 65 20 | 28 73 65 65 20 61 64 64 |from me |(see add|
|000017b0| 72 65 73 73 20 62 65 6c | 6f 77 29 20 2d 20 6a 75 |ress bel|ow) - ju|
|000017c0| 73 74 20 73 65 6e 64 20 | 6d 65 20 33 20 66 6f 72 |st send |me 3 for|
|000017d0| 6d 61 74 74 65 64 20 64 | 69 73 63 73 20 61 6e 64 |matted d|iscs and|
|000017e0| 20 61 6e 20 53 41 45 2e | 0a 0a 4e 2e 42 2e 0a 2d | an SAE.|..N.B..-|
|000017f0| 2d 2d 2d 0a 54 68 65 20 | 41 63 6f 72 6e 20 76 65 |---.The |Acorn ve|
|00001800| 72 73 69 6f 6e 20 6f 66 | 20 4d 65 73 63 68 61 63 |rsion of| Meschac|
|00001810| 68 20 63 6f 75 6c 64 20 | 62 65 20 6d 61 64 65 20 |h could |be made |
|00001820| 72 61 74 68 65 72 20 6d | 6f 72 65 20 65 66 66 69 |rather m|ore effi|
|00001830| 63 69 65 6e 74 20 62 79 | 20 70 6c 61 79 69 6e 67 |cient by| playing|
|00001840| 20 77 69 74 68 20 74 68 | 65 20 74 69 6d 65 2d 63 | with th|e time-c|
|00001850| 72 69 74 69 63 61 6c 20 | 72 6f 75 74 69 6e 65 73 |ritical |routines|
|00001860| 20 69 6e 20 27 6d 61 63 | 68 69 6e 65 2e 63 27 2e | in 'mac|hine.c'.|
|00001870| 20 54 6f 20 62 65 67 69 | 6e 20 77 69 74 68 20 49 | To begi|n with I|
|00001880| 20 68 61 76 65 6e 27 74 | 20 63 68 65 63 6b 65 64 | haven't| checked|
|00001890| 20 77 68 65 74 68 65 72 | 20 75 6e 72 6f 6c 6c 69 | whether| unrolli|
|000018a0| 6e 67 20 6c 6f 6f 70 73 | 20 67 69 76 65 73 20 61 |ng loops| gives a|
|000018b0| 20 73 70 65 65 64 20 69 | 6e 63 72 65 61 73 65 2e | speed i|ncrease.|
|000018c0| 20 42 65 79 6f 6e 64 20 | 74 68 61 74 2c 20 74 68 | Beyond |that, th|
|000018d0| 65 20 72 6f 75 74 69 6e | 65 73 20 63 6f 75 6c 64 |e routin|es could|
|000018e0| 20 62 65 20 63 6f 6e 76 | 65 72 74 65 64 20 69 6e | be conv|erted in|
|000018f0| 74 6f 20 61 73 73 65 6d | 62 6c 65 72 3a 20 75 73 |to assem|bler: us|
|00001900| 69 6e 67 20 4c 46 4d 20 | 61 6e 64 20 53 46 4d 20 |ing LFM |and SFM |
|00001910| 74 6f 20 72 65 61 64 20 | 6d 75 6c 74 69 70 6c 65 |to read |multiple|
|00001920| 20 72 65 67 69 73 74 65 | 72 73 20 69 6e 20 6f 6e | registe|rs in on|
|00001930| 65 20 69 6e 73 74 72 75 | 63 74 69 6f 6e 20 69 73 |e instru|ction is|
|00001940| 20 6f 6e 65 20 77 61 79 | 20 74 6f 20 69 6d 70 72 | one way| to impr|
|00001950| 6f 76 65 20 73 70 65 65 | 64 2e 0a 55 6e 66 6f 72 |ove spee|d..Unfor|
|00001960| 74 75 6e 61 74 65 6c 79 | 2c 20 49 20 68 61 76 65 |tunately|, I have|
|00001970| 6e 27 74 20 67 6f 74 20 | 74 68 65 20 74 69 6d 65 |n't got |the time|
|00001980| 20 74 6f 20 64 6f 20 74 | 68 69 73 20 6d 79 73 65 | to do t|his myse|
|00001990| 6c 66 2e 0a 0a 0a 59 6f | 75 20 63 61 6e 20 63 6f |lf....Yo|u can co|
|000019a0| 6e 74 61 63 74 20 6d 65 | 3a 0a 0a 50 61 75 6c 20 |ntact me|:..Paul |
|000019b0| 46 69 65 6c 64 0a 33 34 | 36 20 43 68 61 64 77 65 |Field.34|6 Chadwe|
|000019c0| 6c 6c 20 48 65 61 74 68 | 20 4c 61 6e 65 0a 43 68 |ll Heath| Lane.Ch|
|000019d0| 61 64 77 65 6c 6c 20 48 | 65 61 74 68 0a 52 6f 6d |adwell H|eath.Rom|
|000019e0| 66 6f 72 64 0a 45 73 73 | 65 78 0a 52 4d 36 20 34 |ford.Ess|ex.RM6 4|
|000019f0| 59 48 0a 45 6e 67 6c 61 | 6e 64 0a 0a 54 65 6c 3a |YH.Engla|nd..Tel:|
|00001a00| 20 30 38 31 20 35 39 38 | 20 39 36 37 36 0a 0a 65 | 081 598| 9676..e|
|00001a10| 6d 61 69 6c 3a 20 70 61 | 75 6c 66 40 64 63 73 2e |mail: pa|ulf@dcs.|
|00001a20| 71 6d 77 2e 61 63 2e 75 | 6b 20 20 28 55 6e 74 69 |qmw.ac.u|k (Unti|
|00001a30| 6c 20 74 68 65 20 65 6e | 64 20 6f 66 20 31 39 39 |l the en|d of 199|
|00001a40| 34 29 46 69 6c 65 44 61 | 74 61 82 21 52 75 6e 00 |4)FileDa|ta.!Run.|
|00001a50| 85 00 08 9b db 8c 7c 00 | 00 00 45 eb ff ff f0 7a |......|.|..E....z|
|00001a60| 08 80 03 00 5c f3 49 63 | 6f 6e 53 70 72 69 74 65 |....\.Ic|onSprite|
|00001a70| 73 20 3c 4f 62 65 79 24 | 44 69 72 3e 2e 21 53 70 |s <Obey$|Dir>.!Sp|
|00001a80| 72 69 74 65 73 0a 0a 7c | 20 53 65 74 20 75 70 20 |rites..|| Set up |
|00001a90| 74 68 65 20 4d 65 73 63 | 68 61 63 68 3a 20 70 61 |the Mesc|hach: pa|
|00001aa0| 74 68 0a 0a 53 65 74 20 | 4d 65 73 63 68 61 63 68 |th..Set |Meschach|
|00001ab0| 24 44 69 72 20 3c 4f 62 | 65 79 24 44 69 72 3e 0a |$Dir <Ob|ey$Dir>.|
|00001ac0| 53 65 74 20 4d 65 73 63 | 68 61 63 68 24 50 61 74 |Set Mesc|hach$Pat|
|00001ad0| 68 20 3c 4d 65 73 63 68 | 61 63 68 24 44 69 72 3e |h <Mesch|ach$Dir>|
|00001ae0| 2e 0a 46 69 6c 65 44 61 | 74 61 82 21 53 70 72 69 |..FileDa|ta.!Spri|
|00001af0| 74 65 73 00 db 8c e0 02 | 00 00 45 f9 ff ff 0b bc |tes.....|..E.....|
|00001b00| 0a 80 03 00 95 7d 01 00 | 00 00 10 00 00 00 e4 02 |.....}..|........|
|00001b10| 00 00 d4 02 00 00 21 6d | 65 73 63 68 61 63 68 00 |......!m|eschach.|
|00001b20| 00 00 04 00 00 00 10 00 | 00 00 00 00 00 00 07 00 |........|........|
|00001b30| 00 00 2c 00 00 00 80 01 | 00 00 0c 00 00 00 77 77 |..,.....|......ww|
|00001b40| 00 00 00 03 00 30 00 00 | 00 00 00 00 00 77 77 00 |.....0..|.....ww.|
|00001b50| 00 00 77 00 03 33 33 02 | 30 33 00 30 33 33 55 03 |..w..33.|03.033U.|
|00001b60| ee ee 77 00 00 00 77 30 | 30 43 76 05 03 30 00 30 |..w...w0|0Cv..0.0|
|00001b70| 22 76 44 77 ee eb 77 00 | 00 00 77 02 30 33 76 05 |"vDw..w.|..w.03v.|
|00001b80| 30 33 00 30 44 57 54 75 | e4 eb 77 00 00 00 77 12 |03.0DWTu|..w...w.|
|00001b90| 00 53 75 05 00 00 00 00 | 55 57 5e 75 e5 eb 77 00 |.Su.....|UW^u..w.|
|00001ba0| 00 00 77 00 00 50 75 05 | 00 00 00 00 33 57 5b 75 |..w..Pu.|....3W[u|
|00001bb0| e3 bb 77 00 00 00 77 00 | 00 40 74 04 00 00 00 ee |..w...w.|.@t.....|
|00001bc0| 3e 73 77 47 b3 bb 77 00 | 00 00 77 00 00 00 74 04 |>swG..w.|..w...t.|
|00001bd0| 00 00 ee ee eb 73 77 47 | bb 1b 77 00 00 00 77 12 |.....swG|..w...w.|
|00001be0| 00 00 00 00 00 ee ee eb | eb eb bb bb 1b 11 77 00 |........|......w.|
|00001bf0| 00 00 77 00 00 00 00 00 | ee ee eb eb eb bb bb 1b |..w.....|........|
|00001c00| 11 00 77 00 00 00 77 00 | 00 33 55 e3 ee eb eb eb |..w...w.|.3U.....|
|00001c10| bb bb 33 12 00 00 77 00 | 00 00 77 00 22 76 44 77 |..3...w.|..w."vDw|
|00001c20| 9c eb eb bb 3b 43 76 05 | 00 00 77 00 00 00 77 00 |....;Cv.|..w...w.|
|00001c30| 44 57 54 75 e4 ee be bb | 3b 33 76 05 00 00 77 00 |DWTu....|;3v...w.|
|00001c40| 00 00 77 00 55 57 59 75 | 25 2e be 1b 11 53 75 35 |..w.UWYu|%....Su5|
|00001c50| 00 00 77 00 00 00 77 00 | 33 57 57 75 33 2e 42 14 |..w...w.|3WWu3.B.|
|00001c60| 30 53 75 35 03 00 77 00 | 00 00 77 00 30 73 77 47 |0Su5..w.|..w.0swG|
|00001c70| 43 11 24 32 30 40 74 34 | 00 00 77 00 00 00 77 77 |C.$20@t4|..w...ww|
|00001c80| 30 73 77 47 13 00 30 03 | 30 00 74 24 03 77 77 00 |0swG..0.|0.t$.ww.|
|00001c90| 00 00 ff ff 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001ca0| 00 ff ff ff ff ff ff 00 | 00 00 f0 0f 00 00 00 00 |........|........|
|00001cb0| 00 f0 ff 0f 00 00 ff ff | ff ff ff 00 00 ff ff 0f |........|........|
|00001cc0| 00 00 00 00 f0 ff ff ff | 00 00 ff ff ff ff ff 00 |........|........|
|00001cd0| 00 ff ff 0f 00 00 00 00 | f0 ff 00 ff 0f 00 ff ff |........|........|
|00001ce0| ff ff ff 00 00 00 ff 0f | 00 00 00 00 f0 ff 00 ff |........|........|
|00001cf0| 0f 00 ff ff ff ff ff 00 | 00 00 ff 0f 00 00 00 00 |........|........|
|00001d00| f0 ff 00 ff 0f 00 ff ff | ff ff ff 00 00 00 ff 0f |........|........|
|00001d10| 00 00 00 00 00 ff ff ff | 00 00 ff ff ff ff ff 00 |........|........|
|00001d20| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 ff ff |........|........|
|00001d30| ff ff ff 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001d40| 00 00 ff ff ff ff ff 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001d50| 00 00 00 00 00 00 ff ff | ff ff ff 00 00 f0 ff 0f |........|........|
|00001d60| 00 00 00 00 00 00 f0 0f | 00 00 ff ff ff ff ff 00 |........|........|
|00001d70| f0 ff ff ff 00 00 00 00 | 00 ff ff 0f 00 00 ff ff |........|........|
|00001d80| ff ff ff 00 f0 ff 00 ff | 0f 00 00 00 00 ff ff 0f |........|........|
|00001d90| 00 00 ff ff ff ff ff 00 | f0 ff 00 ff 0f 00 00 00 |........|........|
|00001da0| 00 00 ff 0f 00 00 ff ff | ff ff ff 00 f0 ff 00 ff |........|........|
|00001db0| 0f 00 00 00 00 00 ff 0f | 00 00 ff ff ff ff ff 00 |........|........|
|00001dc0| 00 ff ff ff 00 00 00 00 | 00 00 ff 0f 00 00 ff ff |........|........|
|00001dd0| ff ff ff ff 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001de0| 00 ff ff ff ff ff 46 69 | 6c 65 44 61 74 61 82 21 |......Fi|leData.!|
|00001df0| 53 70 72 69 74 65 73 32 | 32 00 88 05 00 00 45 f9 |Sprites2|2.....E.|
|00001e00| ff ff 84 20 0a 80 03 00 | 8c 55 01 00 00 00 10 00 |... ....|.U......|
|00001e10| 00 00 8c 05 00 00 7c 05 | 00 00 21 6d 65 73 63 68 |......|.|..!mesch|
|00001e20| 61 63 68 00 00 00 04 00 | 00 00 21 00 00 00 00 00 |ach.....|..!.....|
|00001e30| 00 00 07 00 00 00 2c 00 | 00 00 d4 02 00 00 1f 00 |......,.|........|
|00001e40| 00 00 77 77 00 00 00 30 | 00 00 03 00 00 00 00 00 |..ww...0|........|
|00001e50| 00 77 77 00 00 00 77 77 | 00 00 00 00 00 00 03 00 |.ww...ww|........|
|00001e60| 00 00 00 00 e0 77 77 00 | 00 00 77 00 03 33 60 34 |.....ww.|..w..3`4|
|00001e70| 00 33 03 00 30 53 77 25 | ee ee 77 00 00 00 77 30 |.3..0Sw%|..w...w0|
|00001e80| 30 20 75 35 30 00 03 00 | 30 75 77 67 ee be 77 00 |0 u50...|0uwg..w.|
|00001e90| 00 00 77 30 20 76 77 35 | 30 00 03 00 20 77 22 76 |..w0 vw5|0... w"v|
|00001ea0| e2 eb 77 00 00 00 77 03 | 30 55 77 35 30 22 33 23 |..w...w.|0Uw50"3#|
|00001eb0| 30 57 e3 75 b3 be 77 00 | 00 00 77 03 00 33 75 35 |0W.u..w.|..w..3u5|
|00001ec0| 42 44 44 44 54 57 ee 75 | e5 eb 77 00 00 00 77 00 |BDDDTW.u|..w...w.|
|00001ed0| 00 00 75 35 44 44 44 44 | 54 57 ee 75 b5 be 77 00 |..u5DDDD|TW.u..w.|
|00001ee0| 00 00 77 03 00 00 75 45 | 44 03 00 20 54 57 e2 75 |..w...uE|D.. TW.u|
|00001ef0| e5 eb 77 00 00 00 77 00 | 00 00 75 45 24 00 00 00 |..w...w.|..uE$...|
|00001f00| 30 57 b4 75 b3 be 77 00 | 00 00 77 00 00 00 75 45 |0W.u..w.|..w...uE|
|00001f10| 02 00 00 00 2e 77 22 76 | e2 bb 77 00 00 00 77 00 |.....w"v|..w...w.|
|00001f20| 00 00 75 35 00 00 00 e0 | ee 75 77 67 be bb 77 00 |..u5....|.uwg..w.|
|00001f30| 00 00 77 00 00 00 74 24 | 00 00 00 ee ee 5b 77 25 |..w...t$|.....[w%|
|00001f40| bb bb 77 00 00 00 77 00 | 00 20 44 04 00 00 e0 ee |..w...w.|. D.....|
|00001f50| be be be be bb bb 77 00 | 00 00 77 00 00 20 44 03 |......w.|..w.. D.|
|00001f60| 00 00 ee ee eb eb eb bb | bb 1b 77 00 00 00 77 00 |........|..w...w.|
|00001f70| 00 30 44 02 00 e0 ee be | be be be bb bb 11 77 00 |.0D.....|......w.|
|00001f80| 00 00 77 03 00 30 44 02 | 00 ee ee eb eb eb bb bb |..w..0D.|........|
|00001f90| 1b 01 77 00 00 00 77 00 | 00 30 44 02 e0 ee be be |..w...w.|.0D.....|
|00001fa0| be be bb bb 11 00 77 00 | 00 00 77 00 00 30 44 02 |......w.|..w..0D.|
|00001fb0| ee ee eb eb eb bb bb 1b | 01 00 77 00 00 00 77 00 |........|..w...w.|
|00001fc0| 30 30 44 e3 ee be be be | be 2b 22 15 00 00 77 00 |00D.....|.+"...w.|
|00001fd0| 00 00 77 00 00 23 44 e4 | ee eb eb eb bb 4b 44 51 |..w..#D.|.....KDQ|
|00001fe0| 00 00 77 00 00 00 77 00 | 00 53 77 25 99 be be be |..w...w.|.Sw%....|
|00001ff0| bb 4b 64 34 00 00 77 00 | 00 00 77 00 00 75 77 67 |.Kd4..w.|..w..uwg|
|00002000| 9c eb eb bb bb 23 75 05 | 00 00 77 00 00 00 77 00 |.....#u.|..w...w.|
|00002010| 20 77 22 76 92 be be bb | 2b 76 77 05 00 00 77 00 | w"v....|+vw...w.|
|00002020| 00 00 77 00 30 57 30 75 | 23 ee be bb 42 55 77 05 |..w.0W0u|#...BUw.|
|00002030| 00 00 77 00 00 00 77 00 | 50 57 9e 75 45 24 22 42 |..w...w.|PW.uE$"B|
|00002040| 44 34 75 05 00 00 77 00 | 00 00 77 00 50 57 99 75 |D4u...w.|..w.PW.u|
|00002050| 45 44 44 44 44 03 75 35 | 00 00 77 00 00 00 77 00 |EDDDD.u5|..w...w.|
|00002060| 50 57 e9 75 45 44 44 44 | 23 00 75 35 00 00 77 00 |PW.uEDDD|#.u5..w.|
|00002070| 00 00 77 00 30 57 44 75 | 33 2e 22 05 30 33 75 35 |..w.0WDu|3.".03u5|
|00002080| 03 00 77 00 00 00 77 00 | 20 77 22 76 52 11 05 30 |..w...w.| w"vR..0|
|00002090| 30 00 75 35 00 00 77 00 | 00 00 77 00 30 75 77 67 |0.u5..w.|..w.0uwg|
|000020a0| 51 00 33 33 30 00 75 35 | 00 00 77 00 00 00 77 00 |Q.330.u5|..w...w.|
|000020b0| 30 50 77 25 03 00 03 00 | 30 00 74 34 00 00 77 00 |0Pw%....|0.t4..w.|
|000020c0| 00 00 77 77 30 33 03 00 | 03 00 30 03 30 00 03 00 |..ww03..|..0.0...|
|000020d0| 03 77 77 00 00 00 77 77 | 00 00 00 00 00 00 00 00 |.ww...ww|........|
|000020e0| 00 00 00 00 00 77 77 00 | 00 00 ff ff 00 00 00 00 |.....ww.|........|
|000020f0| 00 00 00 00 00 00 00 00 | 00 ff ff ff ff ff ff ff |........|........|
|00002100| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 ff ff ff |........|........|
|00002110| ff ff ff 00 00 00 f0 0f | 00 00 00 00 00 f0 ff ff |........|........|
|00002120| 00 00 ff ff ff ff ff 00 | 00 f0 ff 0f 00 00 00 00 |........|........|
|00002130| 00 ff ff ff 00 00 ff ff | ff ff ff 00 f0 ff ff 0f |........|........|
|00002140| 00 00 00 00 f0 ff ff ff | 0f 00 ff ff ff ff ff 00 |........|........|
|00002150| 00 ff ff 0f 00 00 00 00 | f0 ff 00 ff 0f 00 ff ff |........|........|
|00002160| ff ff ff 00 00 00 ff 0f | 00 00 00 00 f0 ff 00 ff |........|........|
|00002170| 0f 00 ff ff ff ff ff 00 | 00 00 ff 0f 00 00 00 00 |........|........|
|00002180| f0 ff 00 ff 0f 00 ff ff | ff ff ff 00 00 00 ff 0f |........|........|
|00002190| 00 00 00 00 f0 ff 0f ff | 0f 00 ff ff ff ff ff 00 |........|........|
|000021a0| 00 00 ff 0f 00 00 00 00 | f0 ff 0f ff 0f 00 ff ff |........|........|
|000021b0| ff ff ff 00 00 00 ff 0f | 00 00 00 00 f0 ff ff ff |........|........|
|000021c0| 0f 00 ff ff ff ff ff 00 | 00 00 ff 0f 00 00 00 00 |........|........|
|000021d0| 00 ff ff ff 00 00 ff ff | ff ff ff 00 00 00 ff 0f |........|........|
|000021e0| 00 00 00 00 00 f0 ff ff | 00 00 ff ff ff ff ff 00 |........|........|
|000021f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 ff ff |........|........|
|00002200| ff ff ff 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002210| 00 00 ff ff ff ff ff 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002220| 00 00 00 00 00 00 ff ff | ff ff ff 00 00 f0 00 00 |........|........|
|00002230| 00 00 00 00 00 00 00 00 | 00 00 ff ff ff ff ff 00 |........|........|
|00002240| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 ff ff |........|........|
|00002250| ff ff ff 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002260| 00 00 ff ff ff ff ff 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002270| 00 00 00 00 00 00 ff ff | ff ff ff 00 00 00 00 00 |........|........|
|00002280| 00 00 00 00 00 00 00 00 | 00 00 ff ff ff ff ff 00 |........|........|
|00002290| 00 f0 ff ff 00 00 00 00 | 00 00 f0 0f 00 00 ff ff |........|........|
|000022a0| ff ff ff 00 00 ff ff ff | 00 00 00 00 00 ff ff 0f |........|........|
|000022b0| 00 00 ff ff ff ff ff 00 | f0 ff ff ff 0f 00 00 00 |........|........|
|000022c0| f0 ff ff 0f 00 00 ff ff | ff ff ff 00 f0 ff f0 ff |........|........|
|000022d0| 0f 00 00 00 00 ff ff 0f | 00 00 ff ff ff ff ff 00 |........|........|
|000022e0| f0 ff 00 ff 0f 00 00 00 | 00 f0 ff 0f 00 00 ff ff |........|........|
|000022f0| ff ff ff 00 f0 ff 00 ff | 0f 00 00 00 00 00 ff 0f |........|........|
|00002300| 00 00 ff ff ff ff ff 00 | f0 ff 00 ff 0f 00 00 00 |........|........|
|00002310| 00 00 ff 0f 00 00 ff ff | ff ff ff 00 f0 ff 00 ff |........|........|
|00002320| 0f 00 00 00 00 00 ff 0f | 00 00 ff ff ff ff ff 00 |........|........|
|00002330| f0 ff ff ff 0f 00 00 00 | 00 00 ff 0f 00 00 ff ff |........|........|
|00002340| ff ff ff 00 00 ff ff ff | 00 00 00 00 00 00 ff 0f |........|........|
|00002350| 00 00 ff ff ff ff ff 00 | 00 f0 ff ff 00 00 00 00 |........|........|
|00002360| 00 00 ff 0f 00 00 ff ff | ff ff ff ff 00 00 00 00 |........|........|
|00002370| 00 00 00 00 00 00 00 00 | 00 ff ff ff ff ff ff ff |........|........|
|00002380| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 ff ff ff |........|........|
|00002390| ff ff 46 69 6c 65 44 61 | 74 61 82 21 53 70 72 69 |..FileDa|ta.!Spri|
|000023a0| 74 65 73 32 33 00 58 02 | 00 00 45 f9 ff ff 70 62 |tes23.X.|..E...pb|
|000023b0| 0a 80 03 00 2b 18 01 00 | 00 00 10 00 00 00 5c 02 |....+...|......\.|
|000023c0| 00 00 4c 02 00 00 21 6d | 65 73 63 68 61 63 68 00 |..L...!m|eschach.|
|000023d0| 00 00 01 00 00 00 21 00 | 00 00 00 00 00 00 01 00 |......!.|........|
|000023e0| 00 00 2c 00 00 00 3c 01 | 00 00 17 00 00 00 0f 08 |..,...<.|........|
|000023f0| 01 c0 03 00 00 00 0f 00 | 01 e0 03 00 00 00 d3 ce |........|........|
|00002400| e1 17 03 00 00 00 2b 2f | e1 8f 03 00 00 00 cb 2f |......+/|......./|
|00002410| c1 4c 03 00 00 00 e7 2f | e1 bf 03 00 00 00 c7 cf |.L...../|........|
|00002420| e1 5d 03 00 00 00 03 07 | e0 bc 03 00 00 00 07 07 |.]......|........|
|00002430| e0 5c 03 00 00 00 03 07 | e0 bf 03 00 00 00 03 07 |.\......|........|
|00002440| d0 cc 03 00 00 00 03 07 | c8 ef 03 00 00 00 03 07 |........|........|
|00002450| c4 f7 03 00 00 00 03 00 | a2 fa 03 00 00 00 03 00 |........|........|
|00002460| 51 7d 03 00 00 00 03 80 | a8 3e 03 00 00 00 07 40 |Q}......|.>.....@|
|00002470| 54 1f 03 00 00 00 03 20 | aa 0f 03 00 00 00 03 10 |T...... |........|
|00002480| d5 07 03 00 00 00 23 88 | ea 07 03 00 00 00 43 55 |......#.|......CU|
|00002490| f5 09 03 00 00 00 c3 87 | fa 0f 03 00 00 00 c3 4f |........|.......O|
|000024a0| 3d 07 03 00 00 00 c3 ac | df 07 03 00 00 00 e3 7c |=.......|.......||
|000024b0| df 07 03 00 00 00 e3 bf | 0f 07 03 00 00 00 e3 ff |........|........|
|000024c0| 07 0f 03 00 00 00 e3 ff | 01 0f 03 00 00 00 e3 3f |........|.......?|
|000024d0| e6 1f 03 00 00 00 c3 2c | 29 0f 03 00 00 00 e3 2f |.......,|)....../|
|000024e0| 2f 0f 03 00 00 00 a3 17 | 21 0f 03 00 00 00 ef 11 |/.......|!.......|
|000024f0| 26 d1 03 00 00 00 0f 00 | 00 c0 03 00 00 00 0f 00 |&.......|........|
|00002500| 00 c0 ff ff ff ff 0f 00 | 00 c0 ff ff ff ff 03 06 |........|........|
|00002510| 80 07 ff ff ff ff 03 07 | c0 0f ff ff ff ff c3 07 |........|........|
|00002520| c0 0c ff ff ff ff c3 07 | e0 1c ff ff ff ff 03 07 |........|........|
|00002530| 60 18 ff ff ff ff 03 07 | 60 18 ff ff ff ff 03 07 |`.......|`.......|
|00002540| 60 18 ff ff ff ff 03 07 | e0 1c ff ff ff ff 03 07 |`.......|........|
|00002550| c0 0c ff ff ff ff 03 07 | c0 0f ff ff ff ff 03 07 |........|........|
|00002560| 80 07 ff ff ff ff 03 00 | 00 00 ff ff ff ff 03 00 |........|........|
|00002570| 00 00 ff ff ff ff 03 00 | 00 00 ff ff ff ff 03 00 |........|........|
|00002580| 00 00 ff ff ff ff 03 00 | 00 00 ff ff ff ff 03 00 |........|........|
|00002590| 00 00 ff ff ff ff 03 00 | 00 00 ff ff ff ff 03 00 |........|........|
|000025a0| 00 00 ff ff ff ff 83 07 | 00 06 ff ff ff ff c3 0f |........|........|
|000025b0| 00 07 ff ff ff ff c3 0c | c0 07 ff ff ff ff e3 1c |........|........|
|000025c0| c0 07 ff ff ff ff 63 18 | 00 07 ff ff ff ff 63 18 |......c.|......c.|
|000025d0| 00 07 ff ff ff ff 63 18 | 00 07 ff ff ff ff e3 1c |......c.|........|
|000025e0| 00 07 ff ff ff ff c3 0c | 00 07 ff ff ff ff c3 0f |........|........|
|000025f0| 00 07 ff ff ff ff 83 07 | 00 07 ff ff ff ff 0f 00 |........|........|
|00002600| 00 c0 ff ff ff ff 0f 00 | 00 c0 ff ff ff ff 46 69 |........|......Fi|
|00002610| 6c 65 44 61 74 61 82 61 | 72 6e 6f 6c 64 69 00 00 |leData.a|rnoldi..|
|00002620| 00 00 ec 12 00 00 45 ff | ff ff 78 ef 40 17 03 00 |......E.|..x.@...|
|00002630| 7d cf 0a 2f 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |}../****|********|
|00002640| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002650| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002660| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002670| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 0a 2a |********|******.*|
|00002680| 2a 0a 2a 2a 20 43 6f 70 | 79 72 69 67 68 74 20 28 |*.** Cop|yright (|
|00002690| 43 29 20 31 39 39 33 20 | 44 61 76 69 64 20 45 2e |C) 1993 |David E.|
|000026a0| 20 53 74 65 77 61 72 64 | 20 26 20 5a 62 69 67 6e | Steward| & Zbign|
|000026b0| 69 65 77 20 4c 65 79 6b | 2c 20 61 6c 6c 20 72 69 |iew Leyk|, all ri|
|000026c0| 67 68 74 73 20 72 65 73 | 65 72 76 65 64 2e 0a 2a |ghts res|erved..*|
|000026d0| 2a 0a 2a 2a 09 09 09 20 | 20 20 20 20 4d 65 73 63 |*.**... | Mesc|
|000026e0| 68 61 63 68 20 4c 69 62 | 72 61 72 79 0a 2a 2a 20 |hach Lib|rary.** |
|000026f0| 0a 2a 2a 20 54 68 69 73 | 20 4d 65 73 63 68 61 63 |.** This| Meschac|
|00002700| 68 20 4c 69 62 72 61 72 | 79 20 69 73 20 70 72 6f |h Librar|y is pro|
|00002710| 76 69 64 65 64 20 22 61 | 73 20 69 73 22 20 77 69 |vided "a|s is" wi|
|00002720| 74 68 6f 75 74 20 61 6e | 79 20 65 78 70 72 65 73 |thout an|y expres|
|00002730| 73 20 0a 2a 2a 20 6f 72 | 20 69 6d 70 6c 69 65 64 |s .** or| implied|
|00002740| 20 77 61 72 72 61 6e 74 | 79 20 6f 66 20 61 6e 79 | warrant|y of any|
|00002750| 20 6b 69 6e 64 20 77 69 | 74 68 20 72 65 73 70 65 | kind wi|th respe|
|00002760| 63 74 20 74 6f 20 74 68 | 69 73 20 73 6f 66 74 77 |ct to th|is softw|
|00002770| 61 72 65 2e 20 0a 2a 2a | 20 49 6e 20 70 61 72 74 |are. .**| In part|
|00002780| 69 63 75 6c 61 72 20 74 | 68 65 20 61 75 74 68 6f |icular t|he autho|
|00002790| 72 73 20 73 68 61 6c 6c | 20 6e 6f 74 20 62 65 20 |rs shall| not be |
|000027a0| 6c 69 61 62 6c 65 20 66 | 6f 72 20 61 6e 79 20 64 |liable f|or any d|
|000027b0| 69 72 65 63 74 2c 20 0a | 2a 2a 20 69 6e 64 69 72 |irect, .|** indir|
|000027c0| 65 63 74 2c 20 73 70 65 | 63 69 61 6c 2c 20 69 6e |ect, spe|cial, in|
|000027d0| 63 69 64 65 6e 74 61 6c | 20 6f 72 20 63 6f 6e 73 |cidental| or cons|
|000027e0| 65 71 75 65 6e 74 69 61 | 6c 20 64 61 6d 61 67 65 |equentia|l damage|
|000027f0| 73 20 61 72 69 73 69 6e | 67 20 0a 2a 2a 20 69 6e |s arisin|g .** in|
|00002800| 20 61 6e 79 20 77 61 79 | 20 66 72 6f 6d 20 75 73 | any way| from us|
|00002810| 65 20 6f 66 20 74 68 65 | 20 73 6f 66 74 77 61 72 |e of the| softwar|
|00002820| 65 2e 0a 2a 2a 20 0a 2a | 2a 20 45 76 65 72 79 6f |e..** .*|* Everyo|
|00002830| 6e 65 20 69 73 20 67 72 | 61 6e 74 65 64 20 70 65 |ne is gr|anted pe|
|00002840| 72 6d 69 73 73 69 6f 6e | 20 74 6f 20 63 6f 70 79 |rmission| to copy|
|00002850| 2c 20 6d 6f 64 69 66 79 | 20 61 6e 64 20 72 65 64 |, modify| and red|
|00002860| 69 73 74 72 69 62 75 74 | 65 20 74 68 69 73 0a 2a |istribut|e this.*|
|00002870| 2a 20 4d 65 73 63 68 61 | 63 68 20 4c 69 62 72 61 |* Mescha|ch Libra|
|00002880| 72 79 2c 20 70 72 6f 76 | 69 64 65 64 3a 0a 2a 2a |ry, prov|ided:.**|
|00002890| 20 20 31 2e 20 20 41 6c | 6c 20 63 6f 70 69 65 73 | 1. Al|l copies|
|000028a0| 20 63 6f 6e 74 61 69 6e | 20 74 68 69 73 20 63 6f | contain| this co|
|000028b0| 70 79 72 69 67 68 74 20 | 6e 6f 74 69 63 65 2e 0a |pyright |notice..|
|000028c0| 2a 2a 20 20 32 2e 20 20 | 41 6c 6c 20 6d 6f 64 69 |** 2. |All modi|
|000028d0| 66 69 65 64 20 63 6f 70 | 69 65 73 20 73 68 61 6c |fied cop|ies shal|
|000028e0| 6c 20 63 61 72 72 79 20 | 61 20 6e 6f 74 69 63 65 |l carry |a notice|
|000028f0| 20 73 74 61 74 69 6e 67 | 20 77 68 6f 0a 2a 2a 20 | stating| who.** |
|00002900| 20 20 20 20 20 6d 61 64 | 65 20 74 68 65 20 6c 61 | mad|e the la|
|00002910| 73 74 20 6d 6f 64 69 66 | 69 63 61 74 69 6f 6e 20 |st modif|ication |
|00002920| 61 6e 64 20 74 68 65 20 | 64 61 74 65 20 6f 66 20 |and the |date of |
|00002930| 73 75 63 68 20 6d 6f 64 | 69 66 69 63 61 74 69 6f |such mod|ificatio|
|00002940| 6e 2e 0a 2a 2a 20 20 33 | 2e 20 20 4e 6f 20 63 68 |n..** 3|. No ch|
|00002950| 61 72 67 65 20 69 73 20 | 6d 61 64 65 20 66 6f 72 |arge is |made for|
|00002960| 20 74 68 69 73 20 73 6f | 66 74 77 61 72 65 20 6f | this so|ftware o|
|00002970| 72 20 77 6f 72 6b 73 20 | 64 65 72 69 76 65 64 20 |r works |derived |
|00002980| 66 72 6f 6d 20 69 74 2e | 20 20 0a 2a 2a 20 20 20 |from it.| .** |
|00002990| 20 20 20 54 68 69 73 20 | 63 6c 61 75 73 65 20 73 | This |clause s|
|000029a0| 68 61 6c 6c 20 6e 6f 74 | 20 62 65 20 63 6f 6e 73 |hall not| be cons|
|000029b0| 74 72 75 65 64 20 61 73 | 20 63 6f 6e 73 74 72 61 |trued as| constra|
|000029c0| 69 6e 69 6e 67 20 6f 74 | 68 65 72 20 73 6f 66 74 |ining ot|her soft|
|000029d0| 77 61 72 65 0a 2a 2a 20 | 20 20 20 20 20 64 69 73 |ware.** | dis|
|000029e0| 74 72 69 62 75 74 65 64 | 20 6f 6e 20 74 68 65 20 |tributed| on the |
|000029f0| 73 61 6d 65 20 6d 65 64 | 69 75 6d 20 61 73 20 74 |same med|ium as t|
|00002a00| 68 69 73 20 73 6f 66 74 | 77 61 72 65 2c 20 6e 6f |his soft|ware, no|
|00002a10| 72 20 69 73 20 61 0a 2a | 2a 20 20 20 20 20 20 64 |r is a.*|* d|
|00002a20| 69 73 74 72 69 62 75 74 | 69 6f 6e 20 66 65 65 20 |istribut|ion fee |
|00002a30| 63 6f 6e 73 69 64 65 72 | 65 64 20 61 20 63 68 61 |consider|ed a cha|
|00002a40| 72 67 65 2e 0a 2a 2a 0a | 2a 2a 2a 2a 2a 2a 2a 2a |rge..**.|********|
|00002a50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a60| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a70| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a80| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a90| 2a 2a 2a 2f 0a 0a 2f 2a | 0a 09 41 72 6e 6f 6c 64 |***/../*|..Arnold|
|00002aa0| 69 20 6d 65 74 68 6f 64 | 20 66 6f 72 20 66 69 6e |i method| for fin|
|00002ab0| 64 69 6e 67 20 65 69 67 | 65 6e 76 61 6c 75 65 73 |ding eig|envalues|
|00002ac0| 20 6f 66 20 6c 61 72 67 | 65 20 6e 6f 6e 2d 73 79 | of larg|e non-sy|
|00002ad0| 6d 6d 65 74 72 69 63 0a | 09 09 6d 61 74 72 69 63 |mmetric.|..matric|
|00002ae0| 65 73 0a 2a 2f 0a 23 69 | 6e 63 6c 75 64 65 09 3c |es.*/.#i|nclude.<|
|00002af0| 73 74 64 69 6f 2e 68 3e | 0a 23 69 6e 63 6c 75 64 |stdio.h>|.#includ|
|00002b00| 65 09 3c 6d 61 74 68 2e | 68 3e 0a 23 69 6e 63 6c |e.<math.|h>.#incl|
|00002b10| 75 64 65 09 22 6d 61 74 | 72 69 78 2e 68 22 0a 23 |ude."mat|rix.h".#|
|00002b20| 69 6e 63 6c 75 64 65 09 | 22 6d 61 74 72 69 78 32 |include.|"matrix2|
|00002b30| 2e 68 22 0a 23 69 6e 63 | 6c 75 64 65 09 22 73 70 |.h".#inc|lude."sp|
|00002b40| 61 72 73 65 2e 68 22 0a | 0a 73 74 61 74 69 63 20 |arse.h".|.static |
|00002b50| 63 68 61 72 20 72 63 73 | 69 64 5b 5d 20 3d 20 22 |char rcs|id[] = "|
|00002b60| 24 49 64 3a 20 61 72 6e | 6f 6c 64 69 2e 63 2c 76 |$Id: arn|oldi.c,v|
|00002b70| 20 31 2e 33 20 31 39 39 | 34 2f 30 31 2f 31 33 20 | 1.3 199|4/01/13 |
|00002b80| 30 35 3a 34 35 3a 34 30 | 20 64 65 73 20 45 78 70 |05:45:40| des Exp|
|00002b90| 20 24 22 3b 0a 0a 2f 2a | 20 61 72 6e 6f 6c 64 69 | $";../*| arnoldi|
|00002ba0| 20 2d 2d 20 61 6e 20 69 | 6d 70 6c 65 6d 65 6e 74 | -- an i|mplement|
|00002bb0| 61 74 69 6f 6e 20 6f 66 | 20 74 68 65 20 41 72 6e |ation of| the Arn|
|00002bc0| 6f 6c 64 69 20 6d 65 74 | 68 6f 64 20 2a 2f 0a 4d |oldi met|hod */.M|
|00002bd0| 41 54 09 2a 61 72 6e 6f | 6c 64 69 28 41 2c 41 5f |AT.*arno|ldi(A,A_|
|00002be0| 70 61 72 61 6d 2c 78 30 | 2c 6d 2c 68 5f 72 65 6d |param,x0|,m,h_rem|
|00002bf0| 2c 51 2c 48 29 0a 56 45 | 43 09 2a 28 2a 41 29 28 |,Q,H).VE|C.*(*A)(|
|00002c00| 29 3b 0a 76 6f 69 64 09 | 2a 41 5f 70 61 72 61 6d |);.void.|*A_param|
|00002c10| 3b 0a 56 45 43 09 2a 78 | 30 3b 0a 69 6e 74 09 6d |;.VEC.*x|0;.int.m|
|00002c20| 3b 0a 52 65 61 6c 09 2a | 68 5f 72 65 6d 3b 0a 4d |;.Real.*|h_rem;.M|
|00002c30| 41 54 09 2a 51 2c 20 2a | 48 3b 0a 7b 0a 09 73 74 |AT.*Q, *|H;.{..st|
|00002c40| 61 74 69 63 20 56 45 43 | 09 2a 76 3d 56 4e 55 4c |atic VEC|.*v=VNUL|
|00002c50| 4c 2c 20 2a 75 3d 56 4e | 55 4c 4c 2c 20 2a 72 3d |L, *u=VN|ULL, *r=|
|00002c60| 56 4e 55 4c 4c 2c 20 2a | 73 3d 56 4e 55 4c 4c 2c |VNULL, *|s=VNULL,|
|00002c70| 20 2a 74 6d 70 3d 56 4e | 55 4c 4c 3b 0a 09 69 6e | *tmp=VN|ULL;..in|
|00002c80| 74 09 69 3b 0a 09 52 65 | 61 6c 09 68 5f 76 61 6c |t.i;..Re|al.h_val|
|00002c90| 3b 0a 0a 09 69 66 20 28 | 20 21 20 41 20 7c 7c 20 |;...if (| ! A || |
|00002ca0| 21 20 51 20 7c 7c 20 21 | 20 78 30 20 29 0a 09 20 |! Q || !| x0 ).. |
|00002cb0| 20 20 20 65 72 72 6f 72 | 28 45 5f 4e 55 4c 4c 2c | error|(E_NULL,|
|00002cc0| 22 61 72 6e 6f 6c 64 69 | 22 29 3b 0a 09 69 66 20 |"arnoldi|");..if |
|00002cd0| 28 20 6d 20 3c 3d 20 30 | 20 29 0a 09 20 20 20 20 |( m <= 0| ).. |
|00002ce0| 65 72 72 6f 72 28 45 5f | 42 4f 55 4e 44 53 2c 22 |error(E_|BOUNDS,"|
|00002cf0| 61 72 6e 6f 6c 64 69 22 | 29 3b 0a 09 69 66 20 28 |arnoldi"|);..if (|
|00002d00| 20 51 2d 3e 6e 20 21 3d | 20 78 30 2d 3e 64 69 6d | Q->n !=| x0->dim|
|00002d10| 20 7c 7c 09 51 2d 3e 6d | 20 21 3d 20 6d 20 29 0a | ||.Q->m| != m ).|
|00002d20| 09 20 20 20 20 65 72 72 | 6f 72 28 45 5f 53 49 5a |. err|or(E_SIZ|
|00002d30| 45 53 2c 22 61 72 6e 6f | 6c 64 69 22 29 3b 0a 0a |ES,"arno|ldi");..|
|00002d40| 09 6d 5f 7a 65 72 6f 28 | 51 29 3b 0a 09 48 20 3d |.m_zero(|Q);..H =|
|00002d50| 20 6d 5f 72 65 73 69 7a | 65 28 48 2c 6d 2c 6d 29 | m_resiz|e(H,m,m)|
|00002d60| 3b 0a 09 6d 5f 7a 65 72 | 6f 28 48 29 3b 0a 09 75 |;..m_zer|o(H);..u|
|00002d70| 20 3d 20 76 5f 72 65 73 | 69 7a 65 28 75 2c 78 30 | = v_res|ize(u,x0|
|00002d80| 2d 3e 64 69 6d 29 3b 0a | 09 76 20 3d 20 76 5f 72 |->dim);.|.v = v_r|
|00002d90| 65 73 69 7a 65 28 76 2c | 78 30 2d 3e 64 69 6d 29 |esize(v,|x0->dim)|
|00002da0| 3b 0a 09 72 20 3d 20 76 | 5f 72 65 73 69 7a 65 28 |;..r = v|_resize(|
|00002db0| 72 2c 6d 29 3b 0a 09 73 | 20 3d 20 76 5f 72 65 73 |r,m);..s| = v_res|
|00002dc0| 69 7a 65 28 73 2c 6d 29 | 3b 0a 09 74 6d 70 20 3d |ize(s,m)|;..tmp =|
|00002dd0| 20 76 5f 72 65 73 69 7a | 65 28 74 6d 70 2c 78 30 | v_resiz|e(tmp,x0|
|00002de0| 2d 3e 64 69 6d 29 3b 0a | 09 4d 45 4d 5f 53 54 41 |->dim);.|.MEM_STA|
|00002df0| 54 5f 52 45 47 28 75 2c | 54 59 50 45 5f 56 45 43 |T_REG(u,|TYPE_VEC|
|00002e00| 29 3b 0a 09 4d 45 4d 5f | 53 54 41 54 5f 52 45 47 |);..MEM_|STAT_REG|
|00002e10| 28 76 2c 54 59 50 45 5f | 56 45 43 29 3b 0a 09 4d |(v,TYPE_|VEC);..M|
|00002e20| 45 4d 5f 53 54 41 54 5f | 52 45 47 28 72 2c 54 59 |EM_STAT_|REG(r,TY|
|00002e30| 50 45 5f 56 45 43 29 3b | 0a 09 4d 45 4d 5f 53 54 |PE_VEC);|..MEM_ST|
|00002e40| 41 54 5f 52 45 47 28 73 | 2c 54 59 50 45 5f 56 45 |AT_REG(s|,TYPE_VE|
|00002e50| 43 29 3b 0a 09 4d 45 4d | 5f 53 54 41 54 5f 52 45 |C);..MEM|_STAT_RE|
|00002e60| 47 28 74 6d 70 2c 54 59 | 50 45 5f 56 45 43 29 3b |G(tmp,TY|PE_VEC);|
|00002e70| 0a 09 73 76 5f 6d 6c 74 | 28 31 2e 30 2f 76 5f 6e |..sv_mlt|(1.0/v_n|
|00002e80| 6f 72 6d 32 28 78 30 29 | 2c 78 30 2c 76 29 3b 0a |orm2(x0)|,x0,v);.|
|00002e90| 09 66 6f 72 20 28 20 69 | 20 3d 20 30 3b 20 69 20 |.for ( i| = 0; i |
|00002ea0| 3c 20 6d 3b 20 69 2b 2b | 20 29 0a 09 7b 0a 09 20 |< m; i++| )..{.. |
|00002eb0| 20 20 20 73 65 74 5f 72 | 6f 77 28 51 2c 69 2c 76 | set_r|ow(Q,i,v|
|00002ec0| 29 3b 0a 09 20 20 20 20 | 75 20 3d 20 28 2a 41 29 |);.. |u = (*A)|
|00002ed0| 28 41 5f 70 61 72 61 6d | 2c 76 2c 75 29 3b 0a 09 |(A_param|,v,u);..|
|00002ee0| 20 20 20 20 72 20 3d 20 | 6d 76 5f 6d 6c 74 28 51 | r = |mv_mlt(Q|
|00002ef0| 2c 75 2c 72 29 3b 0a 09 | 20 20 20 20 74 6d 70 20 |,u,r);..| tmp |
|00002f00| 3d 20 76 6d 5f 6d 6c 74 | 28 51 2c 72 2c 74 6d 70 |= vm_mlt|(Q,r,tmp|
|00002f10| 29 3b 0a 09 20 20 20 20 | 76 5f 73 75 62 28 75 2c |);.. |v_sub(u,|
|00002f20| 74 6d 70 2c 75 29 3b 0a | 09 20 20 20 20 68 5f 76 |tmp,u);.|. h_v|
|00002f30| 61 6c 20 3d 20 76 5f 6e | 6f 72 6d 32 28 75 29 3b |al = v_n|orm2(u);|
|00002f40| 0a 09 20 20 20 20 2f 2a | 20 69 66 20 75 20 3d 3d |.. /*| if u ==|
|00002f50| 20 30 20 74 68 65 6e 20 | 77 65 20 68 61 76 65 20 | 0 then |we have |
|00002f60| 61 6e 20 65 78 61 63 74 | 20 73 75 62 73 70 61 63 |an exact| subspac|
|00002f70| 65 20 2a 2f 0a 09 20 20 | 20 20 69 66 20 28 20 68 |e */.. | if ( h|
|00002f80| 5f 76 61 6c 20 3d 3d 20 | 30 2e 30 20 29 0a 09 20 |_val == |0.0 ).. |
|00002f90| 20 20 20 7b 0a 09 09 2a | 68 5f 72 65 6d 20 3d 20 | {...*|h_rem = |
|00002fa0| 68 5f 76 61 6c 3b 0a 09 | 09 72 65 74 75 72 6e 20 |h_val;..|.return |
|00002fb0| 48 3b 0a 09 20 20 20 20 | 7d 0a 09 20 20 20 20 2f |H;.. |}.. /|
|00002fc0| 2a 20 69 74 65 72 61 74 | 69 76 65 20 72 65 66 69 |* iterat|ive refi|
|00002fd0| 6e 65 6d 65 6e 74 20 2d | 2d 20 65 6e 73 75 72 65 |nement -|- ensure|
|00002fe0| 73 20 6e 65 61 72 20 6f | 72 74 68 6f 67 6f 6e 61 |s near o|rthogona|
|00002ff0| 6c 69 74 79 20 2a 2f 0a | 09 20 20 20 20 64 6f 20 |lity */.|. do |
|00003000| 7b 0a 09 09 73 20 3d 20 | 6d 76 5f 6d 6c 74 28 51 |{...s = |mv_mlt(Q|
|00003010| 2c 75 2c 73 29 3b 0a 09 | 09 74 6d 70 20 3d 20 76 |,u,s);..|.tmp = v|
|00003020| 6d 5f 6d 6c 74 28 51 2c | 73 2c 74 6d 70 29 3b 0a |m_mlt(Q,|s,tmp);.|
|00003030| 09 09 76 5f 73 75 62 28 | 75 2c 74 6d 70 2c 75 29 |..v_sub(|u,tmp,u)|
|00003040| 3b 0a 09 09 76 5f 61 64 | 64 28 72 2c 73 2c 72 29 |;...v_ad|d(r,s,r)|
|00003050| 3b 0a 09 20 20 20 20 7d | 20 77 68 69 6c 65 20 28 |;.. }| while (|
|00003060| 20 76 5f 6e 6f 72 6d 32 | 28 73 29 20 3e 20 30 2e | v_norm2|(s) > 0.|
|00003070| 31 2a 28 68 5f 76 61 6c | 20 3d 20 76 5f 6e 6f 72 |1*(h_val| = v_nor|
|00003080| 6d 32 28 75 29 29 20 29 | 3b 0a 09 20 20 20 20 2f |m2(u)) )|;.. /|
|00003090| 2a 20 6e 6f 77 20 74 68 | 61 74 20 75 20 69 73 20 |* now th|at u is |
|000030a0| 6e 65 61 72 6c 79 20 6f | 72 74 68 6f 67 6f 6e 61 |nearly o|rthogona|
|000030b0| 6c 20 74 6f 20 51 2c 20 | 75 70 64 61 74 65 20 48 |l to Q, |update H|
|000030c0| 20 2a 2f 0a 09 20 20 20 | 20 73 65 74 5f 63 6f 6c | */.. | set_col|
|000030d0| 28 48 2c 69 2c 72 29 3b | 0a 09 20 20 20 20 69 66 |(H,i,r);|.. if|
|000030e0| 20 28 20 69 20 3d 3d 20 | 6d 2d 31 20 29 0a 09 20 | ( i == |m-1 ).. |
|000030f0| 20 20 20 7b 0a 09 09 2a | 68 5f 72 65 6d 20 3d 20 | {...*|h_rem = |
|00003100| 68 5f 76 61 6c 3b 0a 09 | 09 63 6f 6e 74 69 6e 75 |h_val;..|.continu|
|00003110| 65 3b 0a 09 20 20 20 20 | 7d 0a 09 20 20 20 20 2f |e;.. |}.. /|
|00003120| 2a 20 48 2d 3e 6d 65 5b | 69 2b 31 5d 5b 69 5d 20 |* H->me[|i+1][i] |
|00003130| 3d 20 68 5f 76 61 6c 3b | 20 2a 2f 0a 09 20 20 20 |= h_val;| */.. |
|00003140| 20 6d 5f 73 65 74 5f 76 | 61 6c 28 48 2c 69 2b 31 | m_set_v|al(H,i+1|
|00003150| 2c 69 2c 68 5f 76 61 6c | 29 3b 0a 09 20 20 20 20 |,i,h_val|);.. |
|00003160| 73 76 5f 6d 6c 74 28 31 | 2e 30 2f 68 5f 76 61 6c |sv_mlt(1|.0/h_val|
|00003170| 2c 75 2c 76 29 3b 0a 09 | 7d 0a 0a 09 72 65 74 75 |,u,v);..|}...retu|
|00003180| 72 6e 20 48 3b 0a 7d 0a | 0a 2f 2a 20 73 70 5f 61 |rn H;.}.|./* sp_a|
|00003190| 72 6e 6f 6c 64 69 20 2d | 2d 20 75 73 65 73 20 61 |rnoldi -|- uses a|
|000031a0| 72 6e 6f 6c 64 69 28 29 | 20 77 69 74 68 20 61 6e |rnoldi()| with an|
|000031b0| 20 65 78 70 6c 69 63 69 | 74 20 72 65 70 72 65 73 | explici|t repres|
|000031c0| 65 6e 74 61 74 69 6f 6e | 20 6f 66 20 41 20 2a 2f |entation| of A */|
|000031d0| 0a 4d 41 54 09 2a 73 70 | 5f 61 72 6e 6f 6c 64 69 |.MAT.*sp|_arnoldi|
|000031e0| 28 41 2c 78 30 2c 6d 2c | 68 5f 72 65 6d 2c 51 2c |(A,x0,m,|h_rem,Q,|
|000031f0| 48 29 0a 53 50 4d 41 54 | 09 2a 41 3b 0a 56 45 43 |H).SPMAT|.*A;.VEC|
|00003200| 09 2a 78 30 3b 0a 69 6e | 74 09 6d 3b 0a 52 65 61 |.*x0;.in|t.m;.Rea|
|00003210| 6c 09 2a 68 5f 72 65 6d | 3b 0a 4d 41 54 09 2a 51 |l.*h_rem|;.MAT.*Q|
|00003220| 2c 20 2a 48 3b 0a 7b 09 | 72 65 74 75 72 6e 20 61 |, *H;.{.|return a|
|00003230| 72 6e 6f 6c 64 69 28 73 | 70 5f 6d 76 5f 6d 6c 74 |rnoldi(s|p_mv_mlt|
|00003240| 2c 41 2c 78 30 2c 6d 2c | 68 5f 72 65 6d 2c 51 2c |,A,x0,m,|h_rem,Q,|
|00003250| 48 29 3b 09 7d 0a 0a 2f | 2a 20 67 6d 72 65 73 20 |H);.}../|* gmres |
|00003260| 2d 2d 20 67 65 6e 65 72 | 61 6c 69 73 65 64 20 6d |-- gener|alised m|
|00003270| 69 6e 69 6d 75 6d 20 72 | 65 73 69 64 75 61 6c 20 |inimum r|esidual |
|00003280| 61 6c 67 6f 72 69 74 68 | 6d 20 6f 66 20 53 61 61 |algorith|m of Saa|
|00003290| 64 20 26 20 53 63 68 75 | 6c 74 7a 0a 09 09 53 49 |d & Schu|ltz...SI|
|000032a0| 41 4d 20 4a 2e 20 53 63 | 69 2e 20 53 74 61 74 2e |AM J. Sc|i. Stat.|
|000032b0| 20 43 6f 6d 70 2e 20 76 | 2e 37 2c 20 70 70 2e 38 | Comp. v|.7, pp.8|
|000032c0| 35 36 2d 2d 38 36 39 20 | 28 31 39 38 36 29 0a 09 |56--869 |(1986)..|
|000032d0| 2d 2d 20 79 20 69 73 20 | 6f 76 65 72 77 72 69 74 |-- y is |overwrit|
|000032e0| 74 65 6e 20 77 69 74 68 | 20 74 68 65 20 73 6f 6c |ten with| the sol|
|000032f0| 75 74 69 6f 6e 20 2a 2f | 0a 56 45 43 09 2a 67 6d |ution */|.VEC.*gm|
|00003300| 72 65 73 28 41 2c 41 5f | 70 61 72 61 6d 2c 6d 2c |res(A,A_|param,m,|
|00003310| 51 2c 52 2c 62 2c 74 6f | 6c 2c 78 29 0a 56 45 43 |Q,R,b,to|l,x).VEC|
|00003320| 09 2a 28 2a 41 29 28 29 | 3b 0a 76 6f 69 64 09 2a |.*(*A)()|;.void.*|
|00003330| 41 5f 70 61 72 61 6d 3b | 0a 56 45 43 09 2a 62 2c |A_param;|.VEC.*b,|
|00003340| 20 2a 78 3b 0a 69 6e 74 | 09 6d 3b 0a 4d 41 54 09 | *x;.int|.m;.MAT.|
|00003350| 2a 51 2c 20 2a 52 3b 0a | 64 6f 75 62 6c 65 09 74 |*Q, *R;.|double.t|
|00003360| 6f 6c 3b 0a 7b 0a 20 20 | 20 20 73 74 61 74 69 63 |ol;.{. | static|
|00003370| 20 56 45 43 09 2a 76 3d | 56 4e 55 4c 4c 2c 20 2a | VEC.*v=|VNULL, *|
|00003380| 75 3d 56 4e 55 4c 4c 2c | 20 2a 72 3d 56 4e 55 4c |u=VNULL,| *r=VNUL|
|00003390| 4c 2c 20 2a 74 6d 70 3d | 56 4e 55 4c 4c 2c 20 2a |L, *tmp=|VNULL, *|
|000033a0| 72 68 73 3d 56 4e 55 4c | 4c 3b 0a 20 20 20 20 73 |rhs=VNUL|L;. s|
|000033b0| 74 61 74 69 63 20 56 45 | 43 09 2a 64 69 61 67 3d |tatic VE|C.*diag=|
|000033c0| 56 4e 55 4c 4c 2c 20 2a | 62 65 74 61 3d 56 4e 55 |VNULL, *|beta=VNU|
|000033d0| 4c 4c 3b 0a 20 20 20 20 | 69 6e 74 09 69 3b 0a 20 |LL;. |int.i;. |
|000033e0| 20 20 20 52 65 61 6c 09 | 68 5f 76 61 6c 2c 20 6e | Real.|h_val, n|
|000033f0| 6f 72 6d 5f 62 3b 0a 20 | 20 20 20 0a 20 20 20 20 |orm_b;. | . |
|00003400| 69 66 20 28 20 21 20 41 | 20 7c 7c 20 21 20 51 20 |if ( ! A| || ! Q |
|00003410| 7c 7c 20 21 20 62 20 7c | 7c 20 21 20 52 20 29 0a ||| ! b ||| ! R ).|
|00003420| 09 65 72 72 6f 72 28 45 | 5f 4e 55 4c 4c 2c 22 67 |.error(E|_NULL,"g|
|00003430| 6d 72 65 73 22 29 3b 0a | 20 20 20 20 69 66 20 28 |mres");.| if (|
|00003440| 20 6d 20 3c 3d 20 30 20 | 29 0a 09 65 72 72 6f 72 | m <= 0 |)..error|
|00003450| 28 45 5f 42 4f 55 4e 44 | 53 2c 22 67 6d 72 65 73 |(E_BOUND|S,"gmres|
|00003460| 22 29 3b 0a 20 20 20 20 | 69 66 20 28 20 51 2d 3e |");. |if ( Q->|
|00003470| 6e 20 21 3d 20 62 2d 3e | 64 69 6d 20 7c 7c 20 51 |n != b->|dim || Q|
|00003480| 2d 3e 6d 20 21 3d 20 6d | 20 29 0a 09 65 72 72 6f |->m != m| )..erro|
|00003490| 72 28 45 5f 53 49 5a 45 | 53 2c 22 67 6d 72 65 73 |r(E_SIZE|S,"gmres|
|000034a0| 22 29 3b 0a 20 20 20 20 | 0a 20 20 20 20 78 20 3d |");. |. x =|
|000034b0| 20 76 5f 63 6f 70 79 28 | 62 2c 78 29 3b 0a 20 20 | v_copy(|b,x);. |
|000034c0| 20 20 6d 5f 7a 65 72 6f | 28 51 29 3b 0a 20 20 20 | m_zero|(Q);. |
|000034d0| 20 52 20 3d 20 6d 5f 72 | 65 73 69 7a 65 28 52 2c | R = m_r|esize(R,|
|000034e0| 6d 2b 31 2c 6d 29 3b 0a | 20 20 20 20 6d 5f 7a 65 |m+1,m);.| m_ze|
|000034f0| 72 6f 28 52 29 3b 0a 20 | 20 20 20 75 20 3d 20 76 |ro(R);. | u = v|
|00003500| 5f 72 65 73 69 7a 65 28 | 75 2c 78 2d 3e 64 69 6d |_resize(|u,x->dim|
|00003510| 29 3b 0a 20 20 20 20 76 | 20 3d 20 76 5f 72 65 73 |);. v| = v_res|
|00003520| 69 7a 65 28 76 2c 78 2d | 3e 64 69 6d 29 3b 0a 20 |ize(v,x-|>dim);. |
|00003530| 20 20 20 74 6d 70 20 3d | 20 76 5f 72 65 73 69 7a | tmp =| v_resiz|
|00003540| 65 28 74 6d 70 2c 78 2d | 3e 64 69 6d 29 3b 0a 20 |e(tmp,x-|>dim);. |
|00003550| 20 20 20 72 68 73 20 3d | 20 76 5f 72 65 73 69 7a | rhs =| v_resiz|
|00003560| 65 28 72 68 73 2c 6d 2b | 31 29 3b 0a 20 20 20 20 |e(rhs,m+|1);. |
|00003570| 4d 45 4d 5f 53 54 41 54 | 5f 52 45 47 28 75 2c 54 |MEM_STAT|_REG(u,T|
|00003580| 59 50 45 5f 56 45 43 29 | 3b 0a 20 20 20 20 4d 45 |YPE_VEC)|;. ME|
|00003590| 4d 5f 53 54 41 54 5f 52 | 45 47 28 76 2c 54 59 50 |M_STAT_R|EG(v,TYP|
|000035a0| 45 5f 56 45 43 29 3b 0a | 20 20 20 20 4d 45 4d 5f |E_VEC);.| MEM_|
|000035b0| 53 54 41 54 5f 52 45 47 | 28 72 2c 54 59 50 45 5f |STAT_REG|(r,TYPE_|
|000035c0| 56 45 43 29 3b 0a 20 20 | 20 20 4d 45 4d 5f 53 54 |VEC);. | MEM_ST|
|000035d0| 41 54 5f 52 45 47 28 74 | 6d 70 2c 54 59 50 45 5f |AT_REG(t|mp,TYPE_|
|000035e0| 56 45 43 29 3b 0a 20 20 | 20 20 4d 45 4d 5f 53 54 |VEC);. | MEM_ST|
|000035f0| 41 54 5f 52 45 47 28 72 | 68 73 2c 54 59 50 45 5f |AT_REG(r|hs,TYPE_|
|00003600| 56 45 43 29 3b 0a 20 20 | 20 20 6e 6f 72 6d 5f 62 |VEC);. | norm_b|
|00003610| 20 3d 20 76 5f 6e 6f 72 | 6d 32 28 78 29 3b 0a 20 | = v_nor|m2(x);. |
|00003620| 20 20 20 69 66 20 28 20 | 6e 6f 72 6d 5f 62 20 3d | if ( |norm_b =|
|00003630| 3d 20 30 2e 30 20 29 0a | 09 65 72 72 6f 72 28 45 |= 0.0 ).|.error(E|
|00003640| 5f 52 41 4e 47 45 2c 22 | 67 6d 72 65 73 22 29 3b |_RANGE,"|gmres");|
|00003650| 0a 20 20 20 20 73 76 5f | 6d 6c 74 28 31 2e 30 2f |. sv_|mlt(1.0/|
|00003660| 6e 6f 72 6d 5f 62 2c 78 | 2c 76 29 3b 0a 20 20 20 |norm_b,x|,v);. |
|00003670| 20 0a 20 20 20 20 66 6f | 72 20 28 20 69 20 3d 20 | . fo|r ( i = |
|00003680| 30 3b 20 69 20 3c 20 6d | 3b 20 69 2b 2b 20 29 0a |0; i < m|; i++ ).|
|00003690| 20 20 20 20 7b 0a 09 73 | 65 74 5f 72 6f 77 28 51 | {..s|et_row(Q|
|000036a0| 2c 69 2c 76 29 3b 0a 09 | 74 72 61 63 65 63 61 74 |,i,v);..|tracecat|
|000036b0| 63 68 28 75 20 3d 20 28 | 2a 41 29 28 41 5f 70 61 |ch(u = (|*A)(A_pa|
|000036c0| 72 61 6d 2c 76 2c 75 29 | 2c 22 67 6d 72 65 73 22 |ram,v,u)|,"gmres"|
|000036d0| 29 3b 0a 09 72 20 3d 20 | 6d 76 5f 6d 6c 74 28 51 |);..r = |mv_mlt(Q|
|000036e0| 2c 75 2c 72 29 3b 0a 09 | 74 6d 70 20 3d 20 76 6d |,u,r);..|tmp = vm|
|000036f0| 5f 6d 6c 74 28 51 2c 72 | 2c 74 6d 70 29 3b 0a 09 |_mlt(Q,r|,tmp);..|
|00003700| 76 5f 73 75 62 28 75 2c | 74 6d 70 2c 75 29 3b 0a |v_sub(u,|tmp,u);.|
|00003710| 09 68 5f 76 61 6c 20 3d | 20 76 5f 6e 6f 72 6d 32 |.h_val =| v_norm2|
|00003720| 28 75 29 3b 0a 09 73 65 | 74 5f 63 6f 6c 28 52 2c |(u);..se|t_col(R,|
|00003730| 69 2c 72 29 3b 0a 09 52 | 2d 3e 6d 65 5b 69 2b 31 |i,r);..R|->me[i+1|
|00003740| 5d 5b 69 5d 20 3d 20 68 | 5f 76 61 6c 3b 0a 09 73 |][i] = h|_val;..s|
|00003750| 76 5f 6d 6c 74 28 31 2e | 30 2f 68 5f 76 61 6c 2c |v_mlt(1.|0/h_val,|
|00003760| 75 2c 76 29 3b 0a 20 20 | 20 20 7d 0a 20 20 20 20 |u,v);. | }. |
|00003770| 0a 20 20 20 20 2f 2a 20 | 75 73 65 20 69 20 78 20 |. /* |use i x |
|00003780| 69 20 73 75 62 6d 61 74 | 72 69 78 20 6f 66 20 52 |i submat|rix of R|
|00003790| 20 2a 2f 0a 20 20 20 20 | 52 20 3d 20 6d 5f 72 65 | */. |R = m_re|
|000037a0| 73 69 7a 65 28 52 2c 69 | 2b 31 2c 69 29 3b 0a 20 |size(R,i|+1,i);. |
|000037b0| 20 20 20 72 68 73 20 3d | 20 76 5f 72 65 73 69 7a | rhs =| v_resiz|
|000037c0| 65 28 72 68 73 2c 69 2b | 31 29 3b 0a 20 20 20 20 |e(rhs,i+|1);. |
|000037d0| 76 5f 7a 65 72 6f 28 72 | 68 73 29 3b 0a 20 20 20 |v_zero(r|hs);. |
|000037e0| 20 72 68 73 2d 3e 76 65 | 5b 30 5d 20 3d 20 6e 6f | rhs->ve|[0] = no|
|000037f0| 72 6d 5f 62 3b 0a 20 20 | 20 20 74 6d 70 20 3d 20 |rm_b;. | tmp = |
|00003800| 76 5f 72 65 73 69 7a 65 | 28 74 6d 70 2c 69 29 3b |v_resize|(tmp,i);|
|00003810| 0a 20 20 20 20 64 69 61 | 67 20 3d 20 76 5f 72 65 |. dia|g = v_re|
|00003820| 73 69 7a 65 28 64 69 61 | 67 2c 69 2b 31 29 3b 0a |size(dia|g,i+1);.|
|00003830| 20 20 20 20 62 65 74 61 | 20 3d 20 76 5f 72 65 73 | beta| = v_res|
|00003840| 69 7a 65 28 62 65 74 61 | 2c 69 2b 31 29 3b 0a 20 |ize(beta|,i+1);. |
|00003850| 20 20 20 4d 45 4d 5f 53 | 54 41 54 5f 52 45 47 28 | MEM_S|TAT_REG(|
|00003860| 62 65 74 61 2c 54 59 50 | 45 5f 56 45 43 29 3b 0a |beta,TYP|E_VEC);.|
|00003870| 20 20 20 20 4d 45 4d 5f | 53 54 41 54 5f 52 45 47 | MEM_|STAT_REG|
|00003880| 28 64 69 61 67 2c 54 59 | 50 45 5f 56 45 43 29 3b |(diag,TY|PE_VEC);|
|00003890| 0a 20 20 20 20 51 52 66 | 61 63 74 6f 72 28 52 2c |. QRf|actor(R,|
|000038a0| 64 69 61 67 20 2f 2a 20 | 2c 62 65 74 61 20 2a 2f |diag /* |,beta */|
|000038b0| 29 3b 0a 20 20 20 20 74 | 6d 70 20 3d 20 51 52 73 |);. t|mp = QRs|
|000038c0| 6f 6c 76 65 28 52 2c 64 | 69 61 67 2c 20 2f 2a 20 |olve(R,d|iag, /* |
|000038d0| 62 65 74 61 2c 20 2a 2f | 20 72 68 73 2c 74 6d 70 |beta, */| rhs,tmp|
|000038e0| 29 3b 0a 20 20 20 20 76 | 5f 72 65 73 69 7a 65 28 |);. v|_resize(|
|000038f0| 74 6d 70 2c 6d 29 3b 0a | 20 20 20 20 76 6d 5f 6d |tmp,m);.| vm_m|
|00003900| 6c 74 28 51 2c 74 6d 70 | 2c 78 29 3b 0a 0a 20 20 |lt(Q,tmp|,x);.. |
|00003910| 20 20 72 65 74 75 72 6e | 20 78 3b 0a 7d 0a 46 69 | return| x;.}.Fi|
|00003920| 6c 65 44 61 74 61 82 62 | 64 66 61 63 74 6f 72 00 |leData.b|dfactor.|
|00003930| 00 00 44 34 00 00 45 ff | ff ff d8 1c 7d 57 03 00 |..D4..E.|....}W..|
|00003940| be 45 0a 2f 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.E./****|********|
|00003950| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003960| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003970| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003980| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 0a 2a |********|******.*|
|00003990| 2a 0a 2a 2a 20 43 6f 70 | 79 72 69 67 68 74 20 28 |*.** Cop|yright (|
|000039a0| 43 29 20 31 39 39 33 20 | 44 61 76 69 64 20 45 2e |C) 1993 |David E.|
|000039b0| 20 53 74 65 77 61 72 64 | 20 26 20 5a 62 69 67 6e | Steward| & Zbign|
|000039c0| 69 65 77 20 4c 65 79 6b | 2c 20 61 6c 6c 20 72 69 |iew Leyk|, all ri|
|000039d0| 67 68 74 73 20 72 65 73 | 65 72 76 65 64 2e 0a 2a |ghts res|erved..*|
|000039e0| 2a 0a 2a 2a 09 09 09 20 | 20 20 20 20 4d 65 73 63 |*.**... | Mesc|
|000039f0| 68 61 63 68 20 4c 69 62 | 72 61 72 79 0a 2a 2a 20 |hach Lib|rary.** |
|00003a00| 0a 2a 2a 20 54 68 69 73 | 20 4d 65 73 63 68 61 63 |.** This| Meschac|
|00003a10| 68 20 4c 69 62 72 61 72 | 79 20 69 73 20 70 72 6f |h Librar|y is pro|
|00003a20| 76 69 64 65 64 20 22 61 | 73 20 69 73 22 20 77 69 |vided "a|s is" wi|
|00003a30| 74 68 6f 75 74 20 61 6e | 79 20 65 78 70 72 65 73 |thout an|y expres|
|00003a40| 73 20 0a 2a 2a 20 6f 72 | 20 69 6d 70 6c 69 65 64 |s .** or| implied|
|00003a50| 20 77 61 72 72 61 6e 74 | 79 20 6f 66 20 61 6e 79 | warrant|y of any|
|00003a60| 20 6b 69 6e 64 20 77 69 | 74 68 20 72 65 73 70 65 | kind wi|th respe|
|00003a70| 63 74 20 74 6f 20 74 68 | 69 73 20 73 6f 66 74 77 |ct to th|is softw|
|00003a80| 61 72 65 2e 20 0a 2a 2a | 20 49 6e 20 70 61 72 74 |are. .**| In part|
|00003a90| 69 63 75 6c 61 72 20 74 | 68 65 20 61 75 74 68 6f |icular t|he autho|
|00003aa0| 72 73 20 73 68 61 6c 6c | 20 6e 6f 74 20 62 65 20 |rs shall| not be |
|00003ab0| 6c 69 61 62 6c 65 20 66 | 6f 72 20 61 6e 79 20 64 |liable f|or any d|
|00003ac0| 69 72 65 63 74 2c 20 0a | 2a 2a 20 69 6e 64 69 72 |irect, .|** indir|
|00003ad0| 65 63 74 2c 20 73 70 65 | 63 69 61 6c 2c 20 69 6e |ect, spe|cial, in|
|00003ae0| 63 69 64 65 6e 74 61 6c | 20 6f 72 20 63 6f 6e 73 |cidental| or cons|
|00003af0| 65 71 75 65 6e 74 69 61 | 6c 20 64 61 6d 61 67 65 |equentia|l damage|
|00003b00| 73 20 61 72 69 73 69 6e | 67 20 0a 2a 2a 20 69 6e |s arisin|g .** in|
|00003b10| 20 61 6e 79 20 77 61 79 | 20 66 72 6f 6d 20 75 73 | any way| from us|
|00003b20| 65 20 6f 66 20 74 68 65 | 20 73 6f 66 74 77 61 72 |e of the| softwar|
|00003b30| 65 2e 0a 2a 2a 20 0a 2a | 2a 20 45 76 65 72 79 6f |e..** .*|* Everyo|
|00003b40| 6e 65 20 69 73 20 67 72 | 61 6e 74 65 64 20 70 65 |ne is gr|anted pe|
|00003b50| 72 6d 69 73 73 69 6f 6e | 20 74 6f 20 63 6f 70 79 |rmission| to copy|
|00003b60| 2c 20 6d 6f 64 69 66 79 | 20 61 6e 64 20 72 65 64 |, modify| and red|
|00003b70| 69 73 74 72 69 62 75 74 | 65 20 74 68 69 73 0a 2a |istribut|e this.*|
|00003b80| 2a 20 4d 65 73 63 68 61 | 63 68 20 4c 69 62 72 61 |* Mescha|ch Libra|
|00003b90| 72 79 2c 20 70 72 6f 76 | 69 64 65 64 3a 0a 2a 2a |ry, prov|ided:.**|
|00003ba0| 20 20 31 2e 20 20 41 6c | 6c 20 63 6f 70 69 65 73 | 1. Al|l copies|
|00003bb0| 20 63 6f 6e 74 61 69 6e | 20 74 68 69 73 20 63 6f | contain| this co|
|00003bc0| 70 79 72 69 67 68 74 20 | 6e 6f 74 69 63 65 2e 0a |pyright |notice..|
|00003bd0| 2a 2a 20 20 32 2e 20 20 | 41 6c 6c 20 6d 6f 64 69 |** 2. |All modi|
|00003be0| 66 69 65 64 20 63 6f 70 | 69 65 73 20 73 68 61 6c |fied cop|ies shal|
|00003bf0| 6c 20 63 61 72 72 79 20 | 61 20 6e 6f 74 69 63 65 |l carry |a notice|
|00003c00| 20 73 74 61 74 69 6e 67 | 20 77 68 6f 0a 2a 2a 20 | stating| who.** |
|00003c10| 20 20 20 20 20 6d 61 64 | 65 20 74 68 65 20 6c 61 | mad|e the la|
|00003c20| 73 74 20 6d 6f 64 69 66 | 69 63 61 74 69 6f 6e 20 |st modif|ication |
|00003c30| 61 6e 64 20 74 68 65 20 | 64 61 74 65 20 6f 66 20 |and the |date of |
|00003c40| 73 75 63 68 20 6d 6f 64 | 69 66 69 63 61 74 69 6f |such mod|ificatio|
|00003c50| 6e 2e 0a 2a 2a 20 20 33 | 2e 20 20 4e 6f 20 63 68 |n..** 3|. No ch|
|00003c60| 61 72 67 65 20 69 73 20 | 6d 61 64 65 20 66 6f 72 |arge is |made for|
|00003c70| 20 74 68 69 73 20 73 6f | 66 74 77 61 72 65 20 6f | this so|ftware o|
|00003c80| 72 20 77 6f 72 6b 73 20 | 64 65 72 69 76 65 64 20 |r works |derived |
|00003c90| 66 72 6f 6d 20 69 74 2e | 20 20 0a 2a 2a 20 20 20 |from it.| .** |
|00003ca0| 20 20 20 54 68 69 73 20 | 63 6c 61 75 73 65 20 73 | This |clause s|
|00003cb0| 68 61 6c 6c 20 6e 6f 74 | 20 62 65 20 63 6f 6e 73 |hall not| be cons|
|00003cc0| 74 72 75 65 64 20 61 73 | 20 63 6f 6e 73 74 72 61 |trued as| constra|
|00003cd0| 69 6e 69 6e 67 20 6f 74 | 68 65 72 20 73 6f 66 74 |ining ot|her soft|
|00003ce0| 77 61 72 65 0a 2a 2a 20 | 20 20 20 20 20 64 69 73 |ware.** | dis|
|00003cf0| 74 72 69 62 75 74 65 64 | 20 6f 6e 20 74 68 65 20 |tributed| on the |
|00003d00| 73 61 6d 65 20 6d 65 64 | 69 75 6d 20 61 73 20 74 |same med|ium as t|
|00003d10| 68 69 73 20 73 6f 66 74 | 77 61 72 65 2c 20 6e 6f |his soft|ware, no|
|00003d20| 72 20 69 73 20 61 0a 2a | 2a 20 20 20 20 20 20 64 |r is a.*|* d|
|00003d30| 69 73 74 72 69 62 75 74 | 69 6f 6e 20 66 65 65 20 |istribut|ion fee |
|00003d40| 63 6f 6e 73 69 64 65 72 | 65 64 20 61 20 63 68 61 |consider|ed a cha|
|00003d50| 72 67 65 2e 0a 2a 2a 0a | 2a 2a 2a 2a 2a 2a 2a 2a |rge..**.|********|
|00003d60| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003d70| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003d80| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003d90| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003da0| 2a 2a 2a 2f 0a 0a 0a 2f | 2a 0a 20 20 42 61 6e 64 |***/.../|*. Band|
|00003db0| 20 6d 61 74 72 69 78 20 | 66 61 63 74 6f 72 69 73 | matrix |factoris|
|00003dc0| 61 74 69 6f 6e 20 72 6f | 75 74 69 6e 65 73 0a 20 |ation ro|utines. |
|00003dd0| 20 2a 2f 0a 0a 2f 2a 20 | 62 64 66 61 63 74 6f 72 | */../* |bdfactor|
|00003de0| 2e 63 20 20 31 38 2f 31 | 31 2f 39 33 20 2a 2f 0a |.c 18/1|1/93 */.|
|00003df0| 73 74 61 74 69 63 09 63 | 68 61 72 09 72 63 73 69 |static.c|har.rcsi|
|00003e00| 64 5b 5d 20 3d 20 22 24 | 49 64 3a 20 62 64 66 61 |d[] = "$|Id: bdfa|
|00003e10| 63 74 6f 72 2e 63 2c 76 | 20 31 2e 35 20 31 39 39 |ctor.c,v| 1.5 199|
|00003e20| 34 2f 30 35 2f 31 37 20 | 32 33 3a 31 38 3a 33 32 |4/05/17 |23:18:32|
|00003e30| 20 64 65 73 20 45 78 70 | 20 24 22 3b 0a 0a 23 69 | des Exp| $";..#i|
|00003e40| 6e 63 6c 75 64 65 09 3c | 73 74 64 69 6f 2e 68 3e |nclude.<|stdio.h>|
|00003e50| 0a 23 69 6e 63 6c 75 64 | 65 09 3c 6d 61 74 68 2e |.#includ|e.<math.|
|00003e60| 68 3e 0a 23 69 6e 63 6c | 75 64 65 20 20 20 20 20 |h>.#incl|ude |
|00003e70| 20 20 20 22 6d 61 74 72 | 69 78 32 2e 68 22 0a 0a | "matr|ix2.h"..|
|00003e80| 0a 2f 2a 20 67 65 6e 65 | 72 61 74 65 20 62 61 6e |./* gene|rate ban|
|00003e90| 64 20 6d 61 74 72 69 78 | 20 0a 20 20 20 66 6f 72 |d matrix| . for|
|00003ea0| 20 61 20 6d 61 74 72 69 | 78 20 20 77 69 74 68 20 | a matri|x with |
|00003eb0| 6e 20 63 6f 6c 75 6d 6e | 73 2c 0a 20 20 20 6c 62 |n column|s,. lb|
|00003ec0| 20 73 75 62 64 69 61 67 | 6f 6e 61 6c 73 20 61 6e | subdiag|onals an|
|00003ed0| 64 20 75 62 20 73 75 70 | 65 72 64 69 61 67 6f 6e |d ub sup|erdiagon|
|00003ee0| 61 6c 73 3b 0a 0a 20 20 | 20 57 61 79 20 6f 66 20 |als;.. | Way of |
|00003ef0| 73 61 76 69 6e 67 20 61 | 20 62 61 6e 64 20 6f 66 |saving a| band of|
|00003f00| 20 61 20 6d 61 74 72 69 | 78 3a 0a 20 20 20 66 69 | a matri|x:. fi|
|00003f10| 72 73 74 20 77 65 20 73 | 61 76 65 20 73 75 62 64 |rst we s|ave subd|
|00003f20| 69 61 67 6f 6e 61 6c 73 | 20 28 66 72 6f 6d 20 30 |iagonals| (from 0|
|00003f30| 20 74 6f 20 6c 62 2d 31 | 29 3b 0a 20 20 20 74 68 | to lb-1|);. th|
|00003f40| 65 6e 20 6d 61 69 6e 20 | 64 69 61 67 6f 6e 61 6c |en main |diagonal|
|00003f50| 20 28 69 6e 20 74 68 65 | 20 6c 62 20 72 6f 77 29 | (in the| lb row)|
|00003f60| 0a 20 20 20 61 6e 64 20 | 74 68 65 6e 20 73 75 70 |. and |then sup|
|00003f70| 65 72 64 69 61 67 6f 6e | 61 6c 73 20 28 66 72 6f |erdiagon|als (fro|
|00003f80| 6d 20 6c 62 2b 31 20 74 | 6f 20 6c 62 2b 75 62 29 |m lb+1 t|o lb+ub)|
|00003f90| 0a 20 20 20 69 6e 20 73 | 75 63 68 20 61 20 77 61 |. in s|uch a wa|
|00003fa0| 79 20 74 68 61 74 20 74 | 68 65 20 65 6c 65 6d 65 |y that t|he eleme|
|00003fb0| 6e 74 73 20 77 68 69 63 | 68 20 77 65 72 65 20 70 |nts whic|h were p|
|00003fc0| 72 65 76 69 6f 75 73 6c | 79 0a 20 20 20 69 6e 20 |reviousl|y. in |
|00003fd0| 6f 6e 65 20 63 6f 6c 75 | 6d 6e 20 61 72 65 20 6e |one colu|mn are n|
|00003fe0| 6f 77 20 61 6c 73 6f 20 | 69 6e 20 6f 6e 65 20 63 |ow also |in one c|
|00003ff0| 6f 6c 75 6d 6e 0a 2a 2f | 0a 0a 42 41 4e 44 20 2a |olumn.*/|..BAND *|
|00004000| 62 64 5f 67 65 74 28 6c | 62 2c 75 62 2c 6e 29 0a |bd_get(l|b,ub,n).|
|00004010| 69 6e 74 20 6c 62 2c 20 | 75 62 2c 20 6e 3b 0a 7b |int lb, |ub, n;.{|
|00004020| 0a 20 20 20 42 41 4e 44 | 20 2a 41 3b 0a 0a 20 20 |. BAND| *A;.. |
|00004030| 20 69 66 20 28 6c 62 20 | 3c 20 30 20 7c 7c 20 75 | if (lb |< 0 || u|
|00004040| 62 20 3c 20 30 20 7c 7c | 20 6e 20 3c 3d 20 30 29 |b < 0 ||| n <= 0)|
|00004050| 0a 20 20 20 20 20 65 72 | 72 6f 72 28 45 5f 4e 45 |. er|ror(E_NE|
|00004060| 47 2c 22 62 64 5f 67 65 | 74 22 29 3b 0a 0a 20 20 |G,"bd_ge|t");.. |
|00004070| 20 69 66 20 28 28 41 20 | 3d 20 4e 45 57 28 42 41 | if ((A |= NEW(BA|
|00004080| 4e 44 29 29 20 3d 3d 20 | 28 42 41 4e 44 20 2a 29 |ND)) == |(BAND *)|
|00004090| 4e 55 4c 4c 29 0a 20 20 | 20 20 20 65 72 72 6f 72 |NULL). | error|
|000040a0| 28 45 5f 4d 45 4d 2c 22 | 62 64 5f 67 65 74 22 29 |(E_MEM,"|bd_get")|
|000040b0| 3b 0a 20 20 20 65 6c 73 | 65 20 69 66 20 28 6d 65 |;. els|e if (me|
|000040c0| 6d 5f 69 6e 66 6f 5f 69 | 73 5f 6f 6e 28 29 29 20 |m_info_i|s_on()) |
|000040d0| 7b 0a 20 20 20 20 20 20 | 6d 65 6d 5f 62 79 74 65 |{. |mem_byte|
|000040e0| 73 28 54 59 50 45 5f 42 | 41 4e 44 2c 30 2c 73 69 |s(TYPE_B|AND,0,si|
|000040f0| 7a 65 6f 66 28 42 41 4e | 44 29 29 3b 0a 20 20 20 |zeof(BAN|D));. |
|00004100| 20 20 20 6d 65 6d 5f 6e | 75 6d 76 61 72 28 54 59 | mem_n|umvar(TY|
|00004110| 50 45 5f 42 41 4e 44 2c | 31 29 3b 0a 20 20 20 7d |PE_BAND,|1);. }|
|00004120| 0a 0a 20 20 20 6c 62 20 | 3d 20 41 2d 3e 6c 62 20 |.. lb |= A->lb |
|00004130| 3d 20 6d 69 6e 28 6e 2d | 31 2c 6c 62 29 3b 0a 20 |= min(n-|1,lb);. |
|00004140| 20 20 75 62 20 3d 20 41 | 2d 3e 75 62 20 3d 20 6d | ub = A|->ub = m|
|00004150| 69 6e 28 6e 2d 31 2c 75 | 62 29 3b 0a 20 20 20 41 |in(n-1,u|b);. A|
|00004160| 2d 3e 6d 61 74 20 3d 20 | 6d 5f 67 65 74 28 6c 62 |->mat = |m_get(lb|
|00004170| 2b 75 62 2b 31 2c 6e 29 | 3b 0a 20 20 20 72 65 74 |+ub+1,n)|;. ret|
|00004180| 75 72 6e 20 41 3b 0a 7d | 0a 0a 69 6e 74 20 62 64 |urn A;.}|..int bd|
|00004190| 5f 66 72 65 65 28 41 29 | 0a 42 41 4e 44 20 2a 41 |_free(A)|.BAND *A|
|000041a0| 3b 0a 7b 0a 20 20 20 69 | 66 20 28 20 41 20 3d 3d |;.{. i|f ( A ==|
|000041b0| 20 28 42 41 4e 44 20 2a | 29 4e 55 4c 4c 20 7c 7c | (BAND *|)NULL |||
|000041c0| 20 41 2d 3e 6c 62 20 3c | 20 30 20 7c 7c 20 41 2d | A->lb <| 0 || A-|
|000041d0| 3e 75 62 20 3c 20 30 20 | 29 0a 20 20 20 20 20 2f |>ub < 0 |). /|
|000041e0| 2a 20 64 6f 6e 27 74 20 | 74 72 75 73 74 20 69 74 |* don't |trust it|
|000041f0| 20 2a 2f 0a 20 20 20 20 | 20 72 65 74 75 72 6e 20 | */. | return |
|00004200| 28 2d 31 29 3b 0a 0a 20 | 20 20 69 66 20 28 41 2d |(-1);.. | if (A-|
|00004210| 3e 6d 61 74 29 20 6d 5f | 66 72 65 65 28 41 2d 3e |>mat) m_|free(A->|
|00004220| 6d 61 74 29 3b 0a 0a 20 | 20 20 69 66 20 28 6d 65 |mat);.. | if (me|
|00004230| 6d 5f 69 6e 66 6f 5f 69 | 73 5f 6f 6e 28 29 29 20 |m_info_i|s_on()) |
|00004240| 7b 0a 20 20 20 20 20 20 | 6d 65 6d 5f 62 79 74 65 |{. |mem_byte|
|00004250| 73 28 54 59 50 45 5f 42 | 41 4e 44 2c 73 69 7a 65 |s(TYPE_B|AND,size|
|00004260| 6f 66 28 42 41 4e 44 29 | 2c 30 29 3b 0a 20 20 20 |of(BAND)|,0);. |
|00004270| 20 20 20 6d 65 6d 5f 6e | 75 6d 76 61 72 28 54 59 | mem_n|umvar(TY|
|00004280| 50 45 5f 42 41 4e 44 2c | 2d 31 29 3b 0a 20 20 20 |PE_BAND,|-1);. |
|00004290| 7d 0a 0a 20 20 20 66 72 | 65 65 28 28 63 68 61 72 |}.. fr|ee((char|
|000042a0| 20 2a 29 41 29 3b 0a 20 | 20 20 72 65 74 75 72 6e | *)A);. | return|
|000042b0| 20 30 3b 0a 7d 0a 0a 0a | 2f 2a 20 72 65 73 69 7a | 0;.}...|/* resiz|
|000042c0| 65 20 62 61 6e 64 20 6d | 61 74 72 69 78 20 2a 2f |e band m|atrix */|
|000042d0| 0a 0a 42 41 4e 44 20 2a | 62 64 5f 72 65 73 69 7a |..BAND *|bd_resiz|
|000042e0| 65 28 41 2c 6e 65 77 5f | 6c 62 2c 6e 65 77 5f 75 |e(A,new_|lb,new_u|
|000042f0| 62 2c 6e 65 77 5f 6e 29 | 0a 42 41 4e 44 20 2a 41 |b,new_n)|.BAND *A|
|00004300| 3b 0a 69 6e 74 20 6e 65 | 77 5f 6c 62 2c 6e 65 77 |;.int ne|w_lb,new|
|00004310| 5f 75 62 2c 6e 65 77 5f | 6e 3b 0a 7b 0a 20 20 20 |_ub,new_|n;.{. |
|00004320| 69 6e 74 20 6c 62 2c 75 | 62 2c 69 2c 6a 2c 6c 2c |int lb,u|b,i,j,l,|
|00004330| 73 68 69 66 74 2c 75 6d | 69 6e 3b 0a 20 20 20 52 |shift,um|in;. R|
|00004340| 65 61 6c 20 2a 2a 41 76 | 3b 0a 0a 20 20 20 69 66 |eal **Av|;.. if|
|00004350| 20 28 6e 65 77 5f 6c 62 | 20 3c 20 30 20 7c 7c 20 | (new_lb| < 0 || |
|00004360| 6e 65 77 5f 75 62 20 3c | 20 30 20 7c 7c 20 6e 65 |new_ub <| 0 || ne|
|00004370| 77 5f 6e 20 3c 3d 20 30 | 29 0a 20 20 20 20 20 65 |w_n <= 0|). e|
|00004380| 72 72 6f 72 28 45 5f 4e | 45 47 2c 22 62 64 5f 72 |rror(E_N|EG,"bd_r|
|00004390| 65 73 69 7a 65 22 29 3b | 0a 20 20 20 69 66 20 28 |esize");|. if (|
|000043a0| 20 21 20 41 20 29 0a 20 | 20 20 20 20 72 65 74 75 | ! A ). | retu|
|000043b0| 72 6e 20 62 64 5f 67 65 | 74 28 6e 65 77 5f 6c 62 |rn bd_ge|t(new_lb|
|000043c0| 2c 6e 65 77 5f 75 62 2c | 6e 65 77 5f 6e 29 3b 0a |,new_ub,|new_n);.|
|000043d0| 20 20 20 20 69 66 20 28 | 20 41 2d 3e 6c 62 2b 41 | if (| A->lb+A|
|000043e0| 2d 3e 75 62 2b 31 20 3e | 20 41 2d 3e 6d 61 74 2d |->ub+1 >| A->mat-|
|000043f0| 3e 6d 20 29 0a 09 65 72 | 72 6f 72 28 45 5f 49 4e |>m )..er|ror(E_IN|
|00004400| 54 45 52 4e 2c 22 62 64 | 5f 72 65 73 69 7a 65 22 |TERN,"bd|_resize"|
|00004410| 29 3b 0a 0a 20 20 20 69 | 66 20 28 20 41 2d 3e 6c |);.. i|f ( A->l|
|00004420| 62 20 3d 3d 20 6e 65 77 | 5f 6c 62 20 26 26 20 41 |b == new|_lb && A|
|00004430| 2d 3e 75 62 20 3d 3d 20 | 6e 65 77 5f 75 62 20 26 |->ub == |new_ub &|
|00004440| 26 20 41 2d 3e 6d 61 74 | 2d 3e 6e 20 3d 3d 20 6e |& A->mat|->n == n|
|00004450| 65 77 5f 6e 20 29 0a 09 | 72 65 74 75 72 6e 20 41 |ew_n )..|return A|
|00004460| 3b 0a 0a 20 20 20 6c 62 | 20 3d 20 41 2d 3e 6c 62 |;.. lb| = A->lb|
|00004470| 3b 0a 20 20 20 75 62 20 | 3d 20 41 2d 3e 75 62 3b |;. ub |= A->ub;|
|00004480| 0a 20 20 20 41 76 20 3d | 20 41 2d 3e 6d 61 74 2d |. Av =| A->mat-|
|00004490| 3e 6d 65 3b 0a 20 20 20 | 75 6d 69 6e 20 3d 20 6d |>me;. |umin = m|
|000044a0| 69 6e 28 75 62 2c 6e 65 | 77 5f 75 62 29 3b 0a 0a |in(ub,ne|w_ub);..|
|000044b0| 20 20 20 20 2f 2a 20 65 | 6e 73 75 72 65 20 74 68 | /* e|nsure th|
|000044c0| 61 74 20 75 6e 75 73 65 | 64 20 74 72 69 61 6e 67 |at unuse|d triang|
|000044d0| 6c 65 73 20 61 74 20 65 | 64 67 65 73 20 61 72 65 |les at e|dges are|
|000044e0| 20 7a 65 72 6f 27 64 20 | 2a 2f 0a 0a 20 20 20 66 | zero'd |*/.. f|
|000044f0| 6f 72 20 28 20 69 20 3d | 20 30 3b 20 69 20 3c 20 |or ( i =| 0; i < |
|00004500| 6c 62 3b 20 69 2b 2b 20 | 29 0a 20 20 20 20 20 20 |lb; i++ |). |
|00004510| 66 6f 72 20 28 20 6a 20 | 3d 20 41 2d 3e 6d 61 74 |for ( j |= A->mat|
|00004520| 2d 3e 6e 20 2d 20 6c 62 | 20 2b 20 69 3b 20 6a 20 |->n - lb| + i; j |
|00004530| 3c 20 41 2d 3e 6d 61 74 | 2d 3e 6e 3b 20 6a 2b 2b |< A->mat|->n; j++|
|00004540| 20 29 0a 09 41 76 5b 69 | 5d 5b 6a 5d 20 3d 20 30 | )..Av[i|][j] = 0|
|00004550| 2e 30 3b 20 20 0a 20 20 | 20 20 66 6f 72 20 28 20 |.0; . | for ( |
|00004560| 69 20 3d 20 6c 62 2b 31 | 2c 6c 3d 31 3b 20 6c 20 |i = lb+1|,l=1; l |
|00004570| 3c 3d 20 75 6d 69 6e 3b | 20 69 2b 2b 2c 6c 2b 2b |<= umin;| i++,l++|
|00004580| 20 29 0a 20 20 20 20 20 | 20 66 6f 72 20 28 20 6a | ). | for ( j|
|00004590| 20 3d 20 30 3b 20 6a 20 | 3c 20 6c 3b 20 6a 2b 2b | = 0; j |< l; j++|
|000045a0| 20 29 0a 09 41 76 5b 69 | 5d 5b 6a 5d 20 3d 20 30 | )..Av[i|][j] = 0|
|000045b0| 2e 30 3b 20 0a 0a 20 20 | 20 6e 65 77 5f 6c 62 20 |.0; .. | new_lb |
|000045c0| 3d 20 41 2d 3e 6c 62 20 | 3d 20 6d 69 6e 28 6e 65 |= A->lb |= min(ne|
|000045d0| 77 5f 6c 62 2c 6e 65 77 | 5f 6e 2d 31 29 3b 0a 20 |w_lb,new|_n-1);. |
|000045e0| 20 20 6e 65 77 5f 75 62 | 20 3d 20 41 2d 3e 75 62 | new_ub| = A->ub|
|000045f0| 20 3d 20 6d 69 6e 28 6e | 65 77 5f 75 62 2c 6e 65 | = min(n|ew_ub,ne|
|00004600| 77 5f 6e 2d 31 29 3b 0a | 20 20 20 41 2d 3e 6d 61 |w_n-1);.| A->ma|
|00004610| 74 20 3d 20 6d 5f 72 65 | 73 69 7a 65 28 41 2d 3e |t = m_re|size(A->|
|00004620| 6d 61 74 2c 6e 65 77 5f | 6c 62 2b 6e 65 77 5f 75 |mat,new_|lb+new_u|
|00004630| 62 2b 31 2c 6e 65 77 5f | 6e 29 3b 0a 20 20 20 41 |b+1,new_|n);. A|
|00004640| 76 20 3d 20 41 2d 3e 6d | 61 74 2d 3e 6d 65 3b 0a |v = A->m|at->me;.|
|00004650| 0a 20 20 20 2f 2a 20 69 | 66 20 6e 65 77 5f 6c 62 |. /* i|f new_lb|
|00004660| 20 21 3d 20 6c 62 20 74 | 68 65 6e 20 6d 6f 76 65 | != lb t|hen move|
|00004670| 20 74 68 65 20 72 6f 77 | 73 20 74 6f 20 67 65 74 | the row|s to get|
|00004680| 20 74 68 65 20 6d 61 69 | 6e 20 64 69 61 67 20 0a | the mai|n diag .|
|00004690| 20 20 20 20 20 20 69 6e | 20 74 68 65 20 6e 65 77 | in| the new|
|000046a0| 5f 6c 62 20 72 6f 77 20 | 2a 2f 0a 0a 20 20 20 69 |_lb row |*/.. i|
|000046b0| 66 20 28 6e 65 77 5f 6c | 62 20 3e 20 6c 62 29 20 |f (new_l|b > lb) |
|000046c0| 7b 0a 20 20 20 20 20 20 | 73 68 69 66 74 20 3d 20 |{. |shift = |
|000046d0| 6e 65 77 5f 6c 62 2d 6c | 62 3b 0a 0a 20 20 20 20 |new_lb-l|b;.. |
|000046e0| 20 20 66 6f 72 20 28 69 | 3d 6c 62 2b 75 6d 69 6e | for (i|=lb+umin|
|000046f0| 2c 20 6c 3d 69 2b 73 68 | 69 66 74 3b 20 69 20 3e |, l=i+sh|ift; i >|
|00004700| 3d 20 30 3b 20 69 2d 2d | 2c 6c 2d 2d 29 0a 09 4d |= 0; i--|,l--)..M|
|00004710| 45 4d 5f 43 4f 50 59 28 | 41 76 5b 69 5d 2c 41 76 |EM_COPY(|Av[i],Av|
|00004720| 5b 6c 5d 2c 6e 65 77 5f | 6e 2a 73 69 7a 65 6f 66 |[l],new_|n*sizeof|
|00004730| 28 52 65 61 6c 29 29 3b | 0a 20 20 20 20 20 20 66 |(Real));|. f|
|00004740| 6f 72 20 28 6c 3d 73 68 | 69 66 74 2d 31 3b 20 6c |or (l=sh|ift-1; l|
|00004750| 20 3e 3d 20 30 3b 20 6c | 2d 2d 29 0a 09 5f 5f 7a | >= 0; l|--)..__z|
|00004760| 65 72 6f 5f 5f 28 41 76 | 5b 6c 5d 2c 6e 65 77 5f |ero__(Av|[l],new_|
|00004770| 6e 29 3b 0a 20 20 20 7d | 0a 20 20 20 65 6c 73 65 |n);. }|. else|
|00004780| 20 7b 20 2f 2a 20 6e 65 | 77 5f 6c 62 20 3c 20 6c | { /* ne|w_lb < l|
|00004790| 62 20 2a 2f 0a 20 20 20 | 20 20 20 73 68 69 66 74 |b */. | shift|
|000047a0| 20 3d 20 6c 62 20 2d 20 | 6e 65 77 5f 6c 62 3b 0a | = lb - |new_lb;.|
|000047b0| 0a 20 20 20 20 20 20 66 | 6f 72 20 28 69 3d 73 68 |. f|or (i=sh|
|000047c0| 69 66 74 2c 20 6c 3d 30 | 3b 20 69 20 3c 3d 20 6c |ift, l=0|; i <= l|
|000047d0| 62 2b 75 6d 69 6e 3b 20 | 69 2b 2b 2c 6c 2b 2b 29 |b+umin; |i++,l++)|
|000047e0| 0a 09 4d 45 4d 5f 43 4f | 50 59 28 41 76 5b 69 5d |..MEM_CO|PY(Av[i]|
|000047f0| 2c 41 76 5b 6c 5d 2c 6e | 65 77 5f 6e 2a 73 69 7a |,Av[l],n|ew_n*siz|
|00004800| 65 6f 66 28 52 65 61 6c | 29 29 3b 0a 20 20 20 20 |eof(Real|));. |
|00004810| 20 20 66 6f 72 20 28 69 | 3d 6c 62 2b 75 6d 69 6e | for (i|=lb+umin|
|00004820| 2b 31 3b 20 69 20 3c 3d | 20 6e 65 77 5f 6c 62 2b |+1; i <=| new_lb+|
|00004830| 6e 65 77 5f 75 62 3b 20 | 69 2b 2b 29 0a 09 5f 5f |new_ub; |i++)..__|
|00004840| 7a 65 72 6f 5f 5f 28 41 | 76 5b 69 5d 2c 6e 65 77 |zero__(A|v[i],new|
|00004850| 5f 6e 29 3b 0a 20 20 20 | 7d 0a 0a 20 20 20 72 65 |_n);. |}.. re|
|00004860| 74 75 72 6e 20 41 3b 0a | 7d 0a 0a 0a 0a 42 41 4e |turn A;.|}....BAN|
|00004870| 44 20 2a 62 64 5f 63 6f | 70 79 28 41 2c 42 29 0a |D *bd_co|py(A,B).|
|00004880| 42 41 4e 44 20 2a 41 2c | 2a 42 3b 0a 7b 0a 20 20 |BAND *A,|*B;.{. |
|00004890| 20 69 6e 74 20 6c 62 2c | 75 62 2c 69 2c 6a 2c 6e | int lb,|ub,i,j,n|
|000048a0| 3b 0a 20 20 20 0a 20 20 | 20 69 66 20 28 20 21 41 |;. . | if ( !A|
|000048b0| 20 29 0a 20 20 20 20 20 | 65 72 72 6f 72 28 45 5f | ). |error(E_|
|000048c0| 4e 55 4c 4c 2c 22 62 64 | 5f 63 6f 70 79 22 29 3b |NULL,"bd|_copy");|
|000048d0| 0a 0a 20 20 20 69 66 20 | 28 41 20 3d 3d 20 42 29 |.. if |(A == B)|
|000048e0| 20 72 65 74 75 72 6e 20 | 42 3b 0a 20 20 20 0a 20 | return |B;. . |
|000048f0| 20 20 6e 20 3d 20 41 2d | 3e 6d 61 74 2d 3e 6e 3b | n = A-|>mat->n;|
|00004900| 0a 20 20 20 69 66 20 28 | 20 21 42 20 29 0a 20 20 |. if (| !B ). |
|00004910| 20 20 20 42 20 3d 20 62 | 64 5f 67 65 74 28 41 2d | B = b|d_get(A-|
|00004920| 3e 6c 62 2c 41 2d 3e 75 | 62 2c 6e 29 3b 0a 20 20 |>lb,A->u|b,n);. |
|00004930| 20 65 6c 73 65 20 69 66 | 20 28 42 2d 3e 6c 62 20 | else if| (B->lb |
|00004940| 21 3d 20 41 2d 3e 6c 62 | 20 7c 7c 20 42 2d 3e 75 |!= A->lb| || B->u|
|00004950| 62 20 21 3d 20 41 2d 3e | 75 62 20 7c 7c 20 42 2d |b != A->|ub || B-|
|00004960| 3e 6d 61 74 2d 3e 6e 20 | 21 3d 20 6e 20 29 0a 20 |>mat->n |!= n ). |
|00004970| 20 20 20 20 42 20 3d 20 | 62 64 5f 72 65 73 69 7a | B = |bd_resiz|
|00004980| 65 28 42 2c 41 2d 3e 6c | 62 2c 41 2d 3e 75 62 2c |e(B,A->l|b,A->ub,|
|00004990| 6e 29 3b 0a 20 20 20 0a | 20 20 20 69 66 20 28 41 |n);. .| if (A|
|000049a0| 2d 3e 6d 61 74 20 3d 3d | 20 42 2d 3e 6d 61 74 29 |->mat ==| B->mat)|
|000049b0| 20 72 65 74 75 72 6e 20 | 42 3b 0a 20 20 20 75 62 | return |B;. ub|
|000049c0| 20 3d 20 42 2d 3e 75 62 | 20 3d 20 41 2d 3e 75 62 | = B->ub| = A->ub|
|000049d0| 3b 0a 20 20 20 6c 62 20 | 3d 20 42 2d 3e 6c 62 20 |;. lb |= B->lb |
|000049e0| 3d 20 41 2d 3e 6c 62 3b | 0a 0a 20 20 20 66 6f 72 |= A->lb;|.. for|
|000049f0| 20 28 20 69 3d 30 2c 20 | 6a 3d 6e 2d 6c 62 3b 20 | ( i=0, |j=n-lb; |
|00004a00| 69 20 3c 3d 20 6c 62 3b | 20 69 2b 2b 2c 20 6a 2b |i <= lb;| i++, j+|
|00004a10| 2b 20 29 0a 20 20 20 20 | 20 4d 45 4d 5f 43 4f 50 |+ ). | MEM_COP|
|00004a20| 59 28 41 2d 3e 6d 61 74 | 2d 3e 6d 65 5b 69 5d 2c |Y(A->mat|->me[i],|
|00004a30| 42 2d 3e 6d 61 74 2d 3e | 6d 65 5b 69 5d 2c 6a 2a |B->mat->|me[i],j*|
|00004a40| 73 69 7a 65 6f 66 28 52 | 65 61 6c 29 29 3b 20 20 |sizeof(R|eal)); |
|00004a50| 20 0a 0a 20 20 20 66 6f | 72 20 28 20 69 3d 6c 62 | .. fo|r ( i=lb|
|00004a60| 2b 31 2c 20 6a 3d 31 3b | 20 69 20 3c 3d 20 6c 62 |+1, j=1;| i <= lb|
|00004a70| 2b 75 62 3b 20 69 2b 2b | 2c 20 6a 2b 2b 20 29 0a |+ub; i++|, j++ ).|
|00004a80| 20 20 20 20 20 4d 45 4d | 5f 43 4f 50 59 28 41 2d | MEM|_COPY(A-|
|00004a90| 3e 6d 61 74 2d 3e 6d 65 | 5b 69 5d 2b 6a 2c 42 2d |>mat->me|[i]+j,B-|
|00004aa0| 3e 6d 61 74 2d 3e 6d 65 | 5b 69 5d 2b 6a 2c 28 6e |>mat->me|[i]+j,(n|
|00004ab0| 20 2d 20 6a 29 2a 73 69 | 7a 65 6f 66 28 52 65 61 | - j)*si|zeof(Rea|
|00004ac0| 6c 29 29 3b 20 20 20 20 | 20 0a 0a 20 20 20 72 65 |l)); | .. re|
|00004ad0| 74 75 72 6e 20 42 3b 0a | 7d 0a 0a 0a 2f 2a 20 63 |turn B;.|}.../* c|
|00004ae0| 6f 70 79 20 62 61 6e 64 | 20 6d 61 74 72 69 78 20 |opy band| matrix |
|00004af0| 74 6f 20 61 20 73 71 75 | 61 72 65 20 6d 61 74 72 |to a squ|are matr|
|00004b00| 69 78 20 2a 2f 0a 4d 41 | 54 20 2a 62 61 6e 64 32 |ix */.MA|T *band2|
|00004b10| 6d 61 74 28 62 41 2c 41 | 29 0a 42 41 4e 44 20 2a |mat(bA,A|).BAND *|
|00004b20| 62 41 3b 0a 4d 41 54 20 | 2a 41 3b 0a 7b 0a 20 20 |bA;.MAT |*A;.{. |
|00004b30| 20 69 6e 74 20 69 2c 6a | 2c 6c 2c 6e 2c 6e 31 3b | int i,j|,l,n,n1;|
|00004b40| 0a 20 20 20 69 6e 74 20 | 6c 62 2c 20 75 62 3b 0a |. int |lb, ub;.|
|00004b50| 20 20 20 52 65 61 6c 20 | 2a 2a 62 6d 61 74 3b 0a | Real |**bmat;.|
|00004b60| 0a 20 20 20 69 66 20 28 | 20 21 62 41 20 7c 7c 20 |. if (| !bA || |
|00004b70| 21 41 29 0a 20 20 20 20 | 20 65 72 72 6f 72 28 45 |!A). | error(E|
|00004b80| 5f 4e 55 4c 4c 2c 22 62 | 61 6e 64 32 6d 61 74 22 |_NULL,"b|and2mat"|
|00004b90| 29 3b 0a 20 20 20 69 66 | 20 28 20 62 41 2d 3e 6d |);. if| ( bA->m|
|00004ba0| 61 74 20 3d 3d 20 41 20 | 29 0a 20 20 20 20 20 65 |at == A |). e|
|00004bb0| 72 72 6f 72 28 45 5f 49 | 4e 53 49 54 55 2c 22 62 |rror(E_I|NSITU,"b|
|00004bc0| 61 6e 64 32 6d 61 74 22 | 29 3b 0a 0a 20 20 20 75 |and2mat"|);.. u|
|00004bd0| 62 20 3d 20 62 41 2d 3e | 75 62 3b 0a 20 20 20 6c |b = bA->|ub;. l|
|00004be0| 62 20 3d 20 62 41 2d 3e | 6c 62 3b 0a 20 20 20 6e |b = bA->|lb;. n|
|00004bf0| 20 3d 20 62 41 2d 3e 6d | 61 74 2d 3e 6e 3b 0a 20 | = bA->m|at->n;. |
|00004c00| 20 20 6e 31 20 3d 20 6e | 2d 31 3b 0a 20 20 20 62 | n1 = n|-1;. b|
|00004c10| 6d 61 74 20 3d 20 62 41 | 2d 3e 6d 61 74 2d 3e 6d |mat = bA|->mat->m|
|00004c20| 65 3b 0a 0a 20 20 20 41 | 20 3d 20 6d 5f 72 65 73 |e;.. A| = m_res|
|00004c30| 69 7a 65 28 41 2c 6e 2c | 6e 29 3b 0a 20 20 20 6d |ize(A,n,|n);. m|
|00004c40| 5f 7a 65 72 6f 28 41 29 | 3b 0a 0a 20 20 20 66 6f |_zero(A)|;.. fo|
|00004c50| 72 20 28 6a 3d 30 3b 20 | 6a 20 3c 20 6e 3b 20 6a |r (j=0; |j < n; j|
|00004c60| 2b 2b 29 0a 20 20 20 20 | 20 66 6f 72 20 28 69 3d |++). | for (i=|
|00004c70| 6d 69 6e 28 6e 31 2c 6a | 2b 6c 62 29 2c 6c 3d 6c |min(n1,j|+lb),l=l|
|00004c80| 62 2b 6a 2d 69 3b 20 69 | 20 3e 3d 20 6d 61 78 28 |b+j-i; i| >= max(|
|00004c90| 30 2c 6a 2d 75 62 29 3b | 20 69 2d 2d 2c 6c 2b 2b |0,j-ub);| i--,l++|
|00004ca0| 29 0a 20 20 20 20 20 20 | 20 41 2d 3e 6d 65 5b 69 |). | A->me[i|
|00004cb0| 5d 5b 6a 5d 20 3d 20 62 | 6d 61 74 5b 6c 5d 5b 6a |][j] = b|mat[l][j|
|00004cc0| 5d 3b 0a 0a 20 20 20 72 | 65 74 75 72 6e 20 41 3b |];.. r|eturn A;|
|00004cd0| 0a 7d 0a 0a 2f 2a 20 63 | 6f 70 79 20 61 20 73 71 |.}../* c|opy a sq|
|00004ce0| 75 61 72 65 20 6d 61 74 | 72 69 78 20 74 6f 20 61 |uare mat|rix to a|
|00004cf0| 20 62 61 6e 64 20 6d 61 | 74 72 69 78 20 77 69 74 | band ma|trix wit|
|00004d00| 68 20 0a 20 20 20 6c 62 | 20 73 75 62 64 69 61 67 |h . lb| subdiag|
|00004d10| 6f 6e 61 6c 73 20 61 6e | 64 20 75 62 20 73 75 70 |onals an|d ub sup|
|00004d20| 65 72 64 69 61 67 6f 6e | 61 6c 73 20 2a 2f 0a 42 |erdiagon|als */.B|
|00004d30| 41 4e 44 20 2a 6d 61 74 | 32 62 61 6e 64 28 41 2c |AND *mat|2band(A,|
|00004d40| 6c 62 2c 75 62 2c 62 41 | 29 0a 42 41 4e 44 20 2a |lb,ub,bA|).BAND *|
|00004d50| 62 41 3b 0a 4d 41 54 20 | 2a 41 3b 0a 69 6e 74 20 |bA;.MAT |*A;.int |
|00004d60| 6c 62 2c 20 75 62 3b 0a | 7b 0a 20 20 20 69 6e 74 |lb, ub;.|{. int|
|00004d70| 20 69 2c 20 6a 2c 20 6c | 2c 20 6e 31 3b 0a 20 20 | i, j, l|, n1;. |
|00004d80| 20 52 65 61 6c 20 2a 2a | 62 6d 61 74 3b 0a 20 20 | Real **|bmat;. |
|00004d90| 20 0a 20 20 20 69 66 20 | 28 21 20 41 20 7c 7c 20 | . if |(! A || |
|00004da0| 21 20 62 41 29 0a 20 20 | 20 20 20 65 72 72 6f 72 |! bA). | error|
|00004db0| 28 45 5f 4e 55 4c 4c 2c | 22 6d 61 74 32 62 61 6e |(E_NULL,|"mat2ban|
|00004dc0| 64 22 29 3b 0a 20 20 20 | 69 66 20 28 75 62 20 3c |d");. |if (ub <|
|00004dd0| 20 30 20 7c 7c 20 6c 62 | 20 3c 20 30 29 0a 20 20 | 0 || lb| < 0). |
|00004de0| 20 20 20 65 72 72 6f 72 | 28 45 5f 53 49 5a 45 53 | error|(E_SIZES|
|00004df0| 2c 22 6d 61 74 32 62 61 | 6e 64 22 29 3b 0a 20 20 |,"mat2ba|nd");. |
|00004e00| 20 69 66 20 28 62 41 2d | 3e 6d 61 74 20 3d 3d 20 | if (bA-|>mat == |
|00004e10| 41 29 0a 20 20 20 20 20 | 65 72 72 6f 72 28 45 5f |A). |error(E_|
|00004e20| 49 4e 53 49 54 55 2c 22 | 6d 61 74 32 62 61 6e 64 |INSITU,"|mat2band|
|00004e30| 22 29 3b 0a 0a 20 20 20 | 6e 31 20 3d 20 41 2d 3e |");.. |n1 = A->|
|00004e40| 6e 2d 31 3b 0a 20 20 20 | 6c 62 20 3d 20 6d 69 6e |n-1;. |lb = min|
|00004e50| 28 6e 31 2c 6c 62 29 3b | 0a 20 20 20 75 62 20 3d |(n1,lb);|. ub =|
|00004e60| 20 6d 69 6e 28 6e 31 2c | 75 62 29 3b 0a 20 20 20 | min(n1,|ub);. |
|00004e70| 62 41 20 3d 20 62 64 5f | 72 65 73 69 7a 65 28 62 |bA = bd_|resize(b|
|00004e80| 41 2c 6c 62 2c 75 62 2c | 6e 31 2b 31 29 3b 0a 20 |A,lb,ub,|n1+1);. |
|00004e90| 20 20 62 6d 61 74 20 3d | 20 62 41 2d 3e 6d 61 74 | bmat =| bA->mat|
|00004ea0| 2d 3e 6d 65 3b 0a 0a 20 | 20 20 66 6f 72 20 28 6a |->me;.. | for (j|
|00004eb0| 3d 30 3b 20 6a 20 3c 3d | 20 6e 31 3b 20 6a 2b 2b |=0; j <=| n1; j++|
|00004ec0| 29 0a 20 20 20 20 20 66 | 6f 72 20 28 69 3d 6d 69 |). f|or (i=mi|
|00004ed0| 6e 28 6e 31 2c 6a 2b 6c | 62 29 2c 6c 3d 6c 62 2b |n(n1,j+l|b),l=lb+|
|00004ee0| 6a 2d 69 3b 20 69 20 3e | 3d 20 6d 61 78 28 30 2c |j-i; i >|= max(0,|
|00004ef0| 6a 2d 75 62 29 3b 20 69 | 2d 2d 2c 6c 2b 2b 29 0a |j-ub); i|--,l++).|
|00004f00| 20 20 20 20 20 20 20 62 | 6d 61 74 5b 6c 5d 5b 6a | b|mat[l][j|
|00004f10| 5d 20 3d 20 41 2d 3e 6d | 65 5b 69 5d 5b 6a 5d 3b |] = A->m|e[i][j];|
|00004f20| 0a 0a 20 20 20 72 65 74 | 75 72 6e 20 62 41 3b 0a |.. ret|urn bA;.|
|00004f30| 7d 0a 0a 0a 0a 2f 2a 20 | 74 72 61 6e 73 70 6f 73 |}..../* |transpos|
|00004f40| 69 74 69 6f 6e 20 6f 66 | 20 6d 61 74 72 69 78 20 |ition of| matrix |
|00004f50| 69 6e 3b 0a 20 20 20 6f | 75 74 20 2d 20 6d 61 74 |in;. o|ut - mat|
|00004f60| 72 69 78 20 61 66 74 65 | 72 20 74 72 61 6e 73 70 |rix afte|r transp|
|00004f70| 6f 73 69 74 69 6f 6e 3b | 0a 20 20 20 63 61 6e 20 |osition;|. can |
|00004f80| 62 65 20 64 6f 6e 65 20 | 69 6e 20 73 69 74 75 0a |be done |in situ.|
|00004f90| 2a 2f 0a 0a 42 41 4e 44 | 20 2a 62 64 5f 74 72 61 |*/..BAND| *bd_tra|
|00004fa0| 6e 73 70 28 69 6e 2c 6f | 75 74 29 0a 42 41 4e 44 |nsp(in,o|ut).BAND|
|00004fb0| 20 2a 69 6e 2c 20 2a 6f | 75 74 3b 0a 7b 0a 20 20 | *in, *o|ut;.{. |
|00004fc0| 20 69 6e 74 20 69 2c 20 | 6a 2c 20 6a 6a 2c 20 6c | int i, |j, jj, l|
|00004fd0| 2c 20 6b 2c 20 6c 62 2c | 20 75 62 2c 20 6c 75 62 |, k, lb,| ub, lub|
|00004fe0| 2c 20 6e 2c 20 6e 31 3b | 0a 20 20 20 69 6e 74 20 |, n, n1;|. int |
|00004ff0| 69 6e 5f 73 69 74 75 3b | 0a 20 20 20 52 65 61 6c |in_situ;|. Real|
|00005000| 20 20 2a 2a 69 6e 5f 76 | 2c 20 2a 2a 6f 75 74 5f | **in_v|, **out_|
|00005010| 76 3b 0a 20 20 20 0a 20 | 20 20 69 66 20 28 20 69 |v;. . | if ( i|
|00005020| 6e 20 3d 3d 20 28 42 41 | 4e 44 20 2a 29 4e 55 4c |n == (BA|ND *)NUL|
|00005030| 4c 20 7c 7c 20 69 6e 2d | 3e 6d 61 74 20 3d 3d 20 |L || in-|>mat == |
|00005040| 28 4d 41 54 20 2a 29 4e | 55 4c 4c 20 29 0a 20 20 |(MAT *)N|ULL ). |
|00005050| 20 20 20 65 72 72 6f 72 | 28 45 5f 4e 55 4c 4c 2c | error|(E_NULL,|
|00005060| 22 62 64 5f 74 72 61 6e | 73 70 22 29 3b 0a 0a 20 |"bd_tran|sp");.. |
|00005070| 20 20 6c 62 20 3d 20 69 | 6e 2d 3e 6c 62 3b 0a 20 | lb = i|n->lb;. |
|00005080| 20 20 75 62 20 3d 20 69 | 6e 2d 3e 75 62 3b 0a 20 | ub = i|n->ub;. |
|00005090| 20 20 6c 75 62 20 3d 20 | 6c 62 2b 75 62 3b 0a 20 | lub = |lb+ub;. |
|000050a0| 20 20 6e 20 3d 20 69 6e | 2d 3e 6d 61 74 2d 3e 6e | n = in|->mat->n|
|000050b0| 3b 0a 20 20 20 6e 31 20 | 3d 20 6e 2d 31 3b 0a 0a |;. n1 |= n-1;..|
|000050c0| 20 20 20 69 6e 5f 73 69 | 74 75 20 3d 20 28 20 69 | in_si|tu = ( i|
|000050d0| 6e 20 3d 3d 20 6f 75 74 | 20 29 3b 0a 20 20 20 69 |n == out| );. i|
|000050e0| 66 20 28 20 21 20 69 6e | 5f 73 69 74 75 20 29 0a |f ( ! in|_situ ).|
|000050f0| 20 20 20 20 20 20 20 6f | 75 74 20 3d 20 62 64 5f | o|ut = bd_|
|00005100| 72 65 73 69 7a 65 28 6f | 75 74 2c 75 62 2c 6c 62 |resize(o|ut,ub,lb|
|00005110| 2c 6e 29 3b 0a 20 20 20 | 65 6c 73 65 0a 20 20 20 |,n);. |else. |
|00005120| 7b 20 20 20 2f 2a 20 6f | 6e 6c 79 20 6e 65 65 64 |{ /* o|nly need|
|00005130| 20 74 6f 20 73 77 61 70 | 20 6c 62 20 61 6e 64 20 | to swap| lb and |
|00005140| 75 62 20 66 69 65 6c 64 | 73 20 2a 2f 0a 20 20 20 |ub field|s */. |
|00005150| 20 20 20 20 6f 75 74 2d | 3e 6c 62 20 3d 20 75 62 | out-|>lb = ub|
|00005160| 3b 0a 20 20 20 20 20 20 | 20 6f 75 74 2d 3e 75 62 |;. | out->ub|
|00005170| 20 3d 20 6c 62 3b 0a 20 | 20 20 7d 0a 0a 20 20 20 | = lb;. | }.. |
|00005180| 69 6e 5f 76 20 3d 20 69 | 6e 2d 3e 6d 61 74 2d 3e |in_v = i|n->mat->|
|00005190| 6d 65 3b 0a 20 20 20 0a | 20 20 20 69 66 20 28 21 |me;. .| if (!|
|000051a0| 20 69 6e 5f 73 69 74 75 | 29 20 7b 0a 20 20 20 20 | in_situ|) {. |
|000051b0| 20 20 69 6e 74 20 73 68 | 5f 69 6e 2c 73 68 5f 6f | int sh|_in,sh_o|
|000051c0| 75 74 3b 20 0a 0a 20 20 | 20 20 20 20 6f 75 74 5f |ut; .. | out_|
|000051d0| 76 20 3d 20 6f 75 74 2d | 3e 6d 61 74 2d 3e 6d 65 |v = out-|>mat->me|
|000051e0| 3b 0a 20 20 20 20 20 20 | 66 6f 72 20 28 69 3d 30 |;. |for (i=0|
|000051f0| 2c 20 6c 3d 6c 75 62 2c | 20 6b 3d 6c 62 2d 69 3b |, l=lub,| k=lb-i;|
|00005200| 20 69 20 3c 3d 20 6c 75 | 62 3b 20 69 2b 2b 2c 6c | i <= lu|b; i++,l|
|00005210| 2d 2d 2c 6b 2d 2d 29 20 | 7b 0a 09 20 73 68 5f 69 |--,k--) |{.. sh_i|
|00005220| 6e 20 3d 20 6d 61 78 28 | 2d 6b 2c 30 29 3b 0a 09 |n = max(|-k,0);..|
|00005230| 20 73 68 5f 6f 75 74 20 | 3d 20 6d 61 78 28 6b 2c | sh_out |= max(k,|
|00005240| 30 29 3b 0a 09 20 4d 45 | 4d 5f 43 4f 50 59 28 26 |0);.. ME|M_COPY(&|
|00005250| 28 69 6e 5f 76 5b 69 5d | 5b 73 68 5f 69 6e 5d 29 |(in_v[i]|[sh_in])|
|00005260| 2c 26 28 6f 75 74 5f 76 | 5b 6c 5d 5b 73 68 5f 6f |,&(out_v|[l][sh_o|
|00005270| 75 74 5d 29 2c 0a 09 09 | 20 20 28 6e 2d 73 68 5f |ut]),...| (n-sh_|
|00005280| 69 6e 2d 73 68 5f 6f 75 | 74 29 2a 73 69 7a 65 6f |in-sh_ou|t)*sizeo|
|00005290| 66 28 52 65 61 6c 29 29 | 3b 0a 09 20 2f 2a 2a 2a |f(Real))|;.. /***|
|000052a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000052b0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 0a |********|*******.|
|000052c0| 09 20 66 6f 72 20 28 6a | 3d 6e 31 2d 73 68 5f 6f |. for (j|=n1-sh_o|
|000052d0| 75 74 2c 20 6a 6a 3d 6e | 31 2d 73 68 5f 69 6e 3b |ut, jj=n|1-sh_in;|
|000052e0| 20 6a 20 3e 3d 20 73 68 | 5f 69 6e 3b 20 6a 2d 2d | j >= sh|_in; j--|
|000052f0| 2c 6a 6a 2d 2d 29 20 7b | 0a 09 20 20 20 20 6f 75 |,jj--) {|.. ou|
|00005300| 74 5f 76 5b 6c 5d 5b 6a | 6a 5d 20 3d 20 69 6e 5f |t_v[l][j|j] = in_|
|00005310| 76 5b 69 5d 5b 6a 5d 3b | 0a 09 20 7d 0a 09 20 2a |v[i][j];|.. }.. *|
|00005320| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005330| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005340| 2a 2f 0a 20 20 20 20 20 | 20 7d 0a 20 20 20 7d 0a |*/. | }. }.|
|00005350| 20 20 20 65 6c 73 65 20 | 69 66 20 28 75 62 20 3d | else |if (ub =|
|00005360| 3d 20 6c 62 29 20 7b 0a | 20 20 20 20 20 20 52 65 |= lb) {.| Re|
|00005370| 61 6c 20 74 6d 70 3b 0a | 0a 20 20 20 20 20 20 66 |al tmp;.|. f|
|00005380| 6f 72 20 28 69 3d 30 2c | 20 6c 3d 6c 75 62 2c 20 |or (i=0,| l=lub, |
|00005390| 6b 3d 6c 62 2d 69 3b 20 | 69 20 3c 20 6c 62 3b 20 |k=lb-i; |i < lb; |
|000053a0| 69 2b 2b 2c 6c 2d 2d 2c | 6b 2d 2d 29 20 7b 0a 09 |i++,l--,|k--) {..|
|000053b0| 20 66 6f 72 20 28 6a 3d | 6e 31 2d 6b 2c 20 6a 6a | for (j=|n1-k, jj|
|000053c0| 3d 6e 31 3b 20 6a 20 3e | 3d 20 30 3b 20 6a 2d 2d |=n1; j >|= 0; j--|
|000053d0| 2c 6a 6a 2d 2d 29 20 7b | 0a 09 20 20 20 20 74 6d |,jj--) {|.. tm|
|000053e0| 70 20 3d 20 69 6e 5f 76 | 5b 6c 5d 5b 6a 6a 5d 3b |p = in_v|[l][jj];|
|000053f0| 0a 09 20 20 20 20 69 6e | 5f 76 5b 6c 5d 5b 6a 6a |.. in|_v[l][jj|
|00005400| 5d 20 3d 20 69 6e 5f 76 | 5b 69 5d 5b 6a 5d 3b 0a |] = in_v|[i][j];.|
|00005410| 09 20 20 20 20 69 6e 5f | 76 5b 69 5d 5b 6a 5d 20 |. in_|v[i][j] |
|00005420| 3d 20 74 6d 70 3b 0a 09 | 20 7d 0a 20 20 20 20 20 |= tmp;..| }. |
|00005430| 20 7d 0a 20 20 20 7d 0a | 20 20 20 65 6c 73 65 20 | }. }.| else |
|00005440| 69 66 20 28 75 62 20 3e | 20 6c 62 29 20 7b 20 20 |if (ub >| lb) { |
|00005450| 2f 2a 20 68 65 6e 63 65 | 20 69 2d 75 62 20 3c 3d |/* hence| i-ub <=|
|00005460| 20 30 20 26 20 6c 2d 6c | 62 20 3e 3d 20 30 20 2a | 0 & l-l|b >= 0 *|
|00005470| 2f 0a 20 20 20 20 20 20 | 69 6e 74 20 70 2c 70 70 |/. |int p,pp|
|00005480| 2c 6c 62 69 3b 0a 20 20 | 20 20 20 20 0a 20 20 20 |,lbi;. | . |
|00005490| 20 20 20 66 6f 72 20 28 | 69 3d 30 2c 20 6c 3d 6c | for (|i=0, l=l|
|000054a0| 75 62 3b 20 69 20 3c 20 | 28 6c 75 62 2b 31 29 2f |ub; i < |(lub+1)/|
|000054b0| 32 3b 20 69 2b 2b 2c 6c | 2d 2d 29 20 7b 0a 09 20 |2; i++,l|--) {.. |
|000054c0| 6c 62 69 20 3d 20 6c 62 | 2d 69 3b 0a 09 20 66 6f |lbi = lb|-i;.. fo|
|000054d0| 72 20 28 6a 3d 6c 2d 6c | 62 2c 20 6a 6a 3d 30 2c |r (j=l-l|b, jj=0,|
|000054e0| 20 70 3d 6d 61 78 28 2d | 6c 62 69 2c 30 29 2c 20 | p=max(-|lbi,0), |
|000054f0| 70 70 20 3d 20 6d 61 78 | 28 6c 2d 75 62 2c 30 29 |pp = max|(l-ub,0)|
|00005500| 3b 20 6a 20 3c 3d 20 6e | 31 3b 20 0a 09 20 20 20 |; j <= n|1; .. |
|00005510| 20 20 20 6a 2b 2b 2c 6a | 6a 2b 2b 2c 70 2b 2b 2c | j++,j|j++,p++,|
|00005520| 70 70 2b 2b 29 20 7b 0a | 09 20 20 20 20 69 6e 5f |pp++) {.|. in_|
|00005530| 76 5b 6c 5d 5b 70 70 5d | 20 3d 20 69 6e 5f 76 5b |v[l][pp]| = in_v[|
|00005540| 69 5d 5b 70 5d 3b 0a 09 | 20 20 20 20 69 6e 5f 76 |i][p];..| in_v|
|00005550| 5b 69 5d 5b 6a 6a 5d 20 | 3d 20 69 6e 5f 76 5b 6c |[i][jj] |= in_v[l|
|00005560| 5d 5b 6a 5d 3b 0a 09 20 | 7d 0a 09 20 66 6f 72 20 |][j];.. |}.. for |
|00005570| 28 20 20 3b 20 70 20 3c | 3d 20 6e 31 2d 6d 61 78 |( ; p <|= n1-max|
|00005580| 28 6c 62 69 2c 30 29 3b | 20 70 2b 2b 2c 70 70 2b |(lbi,0);| p++,pp+|
|00005590| 2b 29 0a 09 20 20 20 69 | 6e 5f 76 5b 6c 5d 5b 70 |+).. i|n_v[l][p|
|000055a0| 70 5d 20 3d 20 69 6e 5f | 76 5b 69 5d 5b 70 5d 3b |p] = in_|v[i][p];|
|000055b0| 0a 20 20 20 20 20 20 7d | 0a 20 20 20 20 20 20 0a |. }|. .|
|000055c0| 20 20 20 20 20 20 69 66 | 20 28 6c 75 62 25 32 20 | if| (lub%2 |
|000055d0| 3d 3d 20 30 29 20 7b 20 | 2f 2a 20 73 68 69 66 74 |== 0) { |/* shift|
|000055e0| 20 6f 6e 6c 79 20 2a 2f | 0a 09 20 69 20 3d 20 6c | only */|.. i = l|
|000055f0| 75 62 2f 32 3b 0a 09 20 | 66 6f 72 20 28 6a 3d 6d |ub/2;.. |for (j=m|
|00005600| 61 78 28 69 2d 6c 62 2c | 30 29 2c 20 6a 6a 3d 30 |ax(i-lb,|0), jj=0|
|00005610| 3b 20 6a 6a 20 3c 3d 20 | 6e 31 2d 75 62 2b 69 3b |; jj <= |n1-ub+i;|
|00005620| 20 6a 2b 2b 2c 6a 6a 2b | 2b 29 20 0a 09 20 20 20 | j++,jj+|+) .. |
|00005630| 69 6e 5f 76 5b 69 5d 5b | 6a 6a 5d 20 3d 20 69 6e |in_v[i][|jj] = in|
|00005640| 5f 76 5b 69 5d 5b 6a 5d | 3b 0a 20 20 20 20 20 20 |_v[i][j]|;. |
|00005650| 7d 0a 20 20 20 7d 0a 20 | 20 20 65 6c 73 65 20 7b |}. }. | else {|
|00005660| 20 20 20 20 20 20 2f 2a | 20 75 62 20 3c 20 6c 62 | /*| ub < lb|
|00005670| 2c 20 68 65 6e 63 65 20 | 75 62 2d 6c 20 3c 3d 20 |, hence |ub-l <= |
|00005680| 30 20 26 20 6c 62 2d 69 | 20 3e 3d 20 30 20 2a 2f |0 & lb-i| >= 0 */|
|00005690| 0a 20 20 20 20 20 20 69 | 6e 74 20 70 2c 70 70 2c |. i|nt p,pp,|
|000056a0| 75 62 69 3b 0a 0a 20 20 | 20 20 20 20 66 6f 72 20 |ubi;.. | for |
|000056b0| 28 69 3d 30 2c 20 6c 3d | 6c 75 62 3b 20 69 20 3c |(i=0, l=|lub; i <|
|000056c0| 20 28 6c 75 62 2b 31 29 | 2f 32 3b 20 69 2b 2b 2c | (lub+1)|/2; i++,|
|000056d0| 6c 2d 2d 29 20 7b 0a 09 | 20 75 62 69 20 3d 20 69 |l--) {..| ubi = i|
|000056e0| 2d 75 62 3b 0a 09 20 66 | 6f 72 20 28 6a 3d 6e 31 |-ub;.. f|or (j=n1|
|000056f0| 2d 6d 61 78 28 6c 62 2d | 6c 2c 30 29 2c 20 6a 6a |-max(lb-|l,0), jj|
|00005700| 3d 6e 31 2d 6d 61 78 28 | 2d 75 62 69 2c 30 29 2c |=n1-max(|-ubi,0),|
|00005710| 20 70 3d 6e 31 2d 6c 62 | 2b 69 2c 20 70 70 3d 6e | p=n1-lb|+i, pp=n|
|00005720| 31 3b 0a 09 20 20 20 20 | 20 20 70 20 3e 3d 20 30 |1;.. | p >= 0|
|00005730| 3b 20 6a 2d 2d 2c 20 6a | 6a 2d 2d 2c 20 70 70 2d |; j--, j|j--, pp-|
|00005740| 2d 2c 20 70 2d 2d 29 20 | 7b 0a 09 20 20 20 20 69 |-, p--) |{.. i|
|00005750| 6e 5f 76 5b 69 5d 5b 6a | 6a 5d 20 3d 20 69 6e 5f |n_v[i][j|j] = in_|
|00005760| 76 5b 6c 5d 5b 6a 5d 3b | 0a 09 20 20 20 20 69 6e |v[l][j];|.. in|
|00005770| 5f 76 5b 6c 5d 5b 70 70 | 5d 20 3d 20 69 6e 5f 76 |_v[l][pp|] = in_v|
|00005780| 5b 69 5d 5b 70 5d 3b 0a | 09 20 7d 0a 09 20 66 6f |[i][p];.|. }.. fo|
|00005790| 72 20 28 20 20 3b 20 6a | 6a 20 3e 3d 20 6d 61 78 |r ( ; j|j >= max|
|000057a0| 28 75 62 69 2c 30 29 3b | 20 6a 2d 2d 2c 20 6a 6a |(ubi,0);| j--, jj|
|000057b0| 2d 2d 29 0a 09 20 20 20 | 69 6e 5f 76 5b 69 5d 5b |--).. |in_v[i][|
|000057c0| 6a 6a 5d 20 3d 20 69 6e | 5f 76 5b 6c 5d 5b 6a 5d |jj] = in|_v[l][j]|
|000057d0| 3b 0a 20 20 20 20 20 20 | 7d 0a 0a 20 20 20 20 20 |;. |}.. |
|000057e0| 20 69 66 20 28 6c 75 62 | 25 32 20 3d 3d 20 30 29 | if (lub|%2 == 0)|
|000057f0| 20 7b 20 20 2f 2a 20 73 | 68 69 66 74 20 6f 6e 6c | { /* s|hift onl|
|00005800| 79 20 2a 2f 0a 09 20 69 | 20 3d 20 6c 75 62 2f 32 |y */.. i| = lub/2|
|00005810| 3b 0a 09 20 66 6f 72 20 | 28 6a 3d 6e 31 2d 6c 62 |;.. for |(j=n1-lb|
|00005820| 2b 69 2c 20 6a 6a 3d 6e | 31 2d 6d 61 78 28 75 62 |+i, jj=n|1-max(ub|
|00005830| 2d 69 2c 30 29 3b 20 6a | 20 3e 3d 20 30 3b 20 6a |-i,0); j| >= 0; j|
|00005840| 2d 2d 2c 20 6a 6a 2d 2d | 29 20 0a 09 20 20 20 20 |--, jj--|) .. |
|00005850| 69 6e 5f 76 5b 69 5d 5b | 6a 6a 5d 20 3d 20 69 6e |in_v[i][|jj] = in|
|00005860| 5f 76 5b 69 5d 5b 6a 5d | 3b 0a 20 20 20 20 20 20 |_v[i][j]|;. |
|00005870| 7d 0a 20 20 20 7d 0a 0a | 20 20 20 72 65 74 75 72 |}. }..| retur|
|00005880| 6e 20 6f 75 74 3b 0a 7d | 0a 0a 0a 0a 2f 2a 20 62 |n out;.}|..../* b|
|00005890| 64 4c 55 66 61 63 74 6f | 72 20 2d 2d 20 67 61 75 |dLUfacto|r -- gau|
|000058a0| 73 73 69 61 6e 20 65 6c | 69 6d 69 6e 61 74 69 6f |ssian el|iminatio|
|000058b0| 6e 20 77 69 74 68 20 70 | 61 72 74 69 61 6c 20 70 |n with p|artial p|
|000058c0| 69 76 6f 74 69 6e 67 0a | 20 20 20 2d 2d 20 6f 6e |ivoting.| -- on|
|000058d0| 20 65 6e 74 72 79 2c 20 | 74 68 65 20 6d 61 74 72 | entry, |the matr|
|000058e0| 69 78 20 41 20 69 6e 20 | 62 61 6e 64 20 73 74 6f |ix A in |band sto|
|000058f0| 72 61 67 65 20 77 69 74 | 68 20 65 6c 65 6d 65 6e |rage wit|h elemen|
|00005900| 74 73 20 0a 20 20 20 20 | 20 20 69 6e 20 72 6f 77 |ts . | in row|
|00005910| 73 20 30 20 74 6f 20 6c | 62 2b 75 62 3b 20 0a 20 |s 0 to l|b+ub; . |
|00005920| 20 20 20 20 20 54 68 65 | 20 6a 74 68 20 63 6f 6c | The| jth col|
|00005930| 75 6d 6e 20 6f 66 20 41 | 20 69 73 20 73 74 6f 72 |umn of A| is stor|
|00005940| 65 64 20 69 6e 20 74 68 | 65 20 6a 74 68 20 63 6f |ed in th|e jth co|
|00005950| 6c 75 6d 6e 20 6f 66 20 | 0a 20 20 20 20 20 20 62 |lumn of |. b|
|00005960| 61 6e 64 20 41 20 28 62 | 41 29 20 61 73 20 66 6f |and A (b|A) as fo|
|00005970| 6c 6c 6f 77 73 3a 0a 20 | 20 20 20 20 20 62 41 2d |llows:. | bA-|
|00005980| 3e 6d 61 74 2d 3e 6d 65 | 5b 6c 62 2b 6a 2d 69 5d |>mat->me|[lb+j-i]|
|00005990| 5b 6a 5d 20 3d 20 41 2d | 3e 6d 65 5b 69 5d 5b 6a |[j] = A-|>me[i][j|
|000059a0| 5d 20 66 6f 72 20 0a 20 | 20 20 20 20 20 6d 61 78 |] for . | max|
|000059b0| 28 30 2c 6a 2d 6c 62 29 | 20 3c 3d 20 69 20 3c 3d |(0,j-lb)| <= i <=|
|000059c0| 20 6d 69 6e 28 41 2d 3e | 6e 2d 31 2c 6a 2b 75 62 | min(A->|n-1,j+ub|
|000059d0| 29 3b 0a 20 20 20 2d 2d | 20 6f 6e 20 65 78 69 74 |);. --| on exit|
|000059e0| 3a 20 55 20 69 73 20 73 | 74 6f 72 65 64 20 61 73 |: U is s|tored as|
|000059f0| 20 61 6e 20 75 70 70 65 | 72 20 74 72 69 61 6e 67 | an uppe|r triang|
|00005a00| 75 6c 61 72 20 6d 61 74 | 72 69 78 0a 20 20 20 20 |ular mat|rix. |
|00005a10| 20 20 77 69 74 68 20 6c | 62 2b 75 62 20 73 75 70 | with l|b+ub sup|
|00005a20| 65 72 64 69 61 67 6f 6e | 61 6c 73 20 69 6e 20 72 |erdiagon|als in r|
|00005a30| 6f 77 73 20 6c 62 20 74 | 6f 20 32 2a 6c 62 2b 75 |ows lb t|o 2*lb+u|
|00005a40| 62 2c 20 0a 20 20 20 20 | 20 20 61 6e 64 20 74 68 |b, . | and th|
|00005a50| 65 20 6d 61 74 72 69 78 | 20 4c 20 69 73 20 73 74 |e matrix| L is st|
|00005a60| 6f 72 65 64 20 69 6e 20 | 72 6f 77 73 20 30 20 74 |ored in |rows 0 t|
|00005a70| 6f 20 6c 62 2d 31 2e 0a | 20 20 20 20 20 20 4d 61 |o lb-1..| Ma|
|00005a80| 74 72 69 78 20 55 20 69 | 73 20 70 65 72 6d 75 74 |trix U i|s permut|
|00005a90| 65 64 2c 20 77 68 65 72 | 65 61 73 20 4c 20 69 73 |ed, wher|eas L is|
|00005aa0| 20 6e 6f 74 20 70 65 72 | 6d 75 74 65 64 20 21 21 | not per|muted !!|
|00005ab0| 21 0a 20 20 20 20 20 20 | 54 68 65 72 65 66 6f 72 |!. |Therefor|
|00005ac0| 65 20 77 65 20 73 61 76 | 65 20 73 6f 6d 65 20 6d |e we sav|e some m|
|00005ad0| 65 6d 6f 72 79 2e 0a 20 | 20 20 2a 2f 0a 42 41 4e |emory.. | */.BAN|
|00005ae0| 44 09 2a 62 64 4c 55 66 | 61 63 74 6f 72 28 62 41 |D.*bdLUf|actor(bA|
|00005af0| 2c 70 69 76 6f 74 29 0a | 42 41 4e 44 09 2a 62 41 |,pivot).|BAND.*bA|
|00005b00| 3b 0a 50 45 52 4d 09 2a | 70 69 76 6f 74 3b 0a 7b |;.PERM.*|pivot;.{|
|00005b10| 0a 20 20 20 69 6e 74 09 | 69 2c 20 6a 2c 20 6b 2c |. int.|i, j, k,|
|00005b20| 20 6c 2c 20 6e 2c 20 6e | 31 2c 20 6c 62 2c 20 75 | l, n, n|1, lb, u|
|00005b30| 62 2c 20 6c 75 62 2c 20 | 6b 5f 65 6e 64 2c 20 6b |b, lub, |k_end, k|
|00005b40| 5f 6c 75 62 3b 0a 20 20 | 20 69 6e 74 09 69 5f 6d |_lub;. | int.i_m|
|00005b50| 61 78 2c 20 73 68 69 66 | 74 3b 0a 20 20 20 52 65 |ax, shif|t;. Re|
|00005b60| 61 6c 09 2a 2a 62 41 5f | 76 3b 0a 20 20 20 52 65 |al.**bA_|v;. Re|
|00005b70| 61 6c 20 6d 61 78 31 2c | 20 74 65 6d 70 3b 0a 20 |al max1,| temp;. |
|00005b80| 20 20 0a 20 20 20 69 66 | 20 28 20 62 41 3d 3d 28 | . if| ( bA==(|
|00005b90| 42 41 4e 44 20 2a 29 4e | 55 4c 4c 20 7c 7c 20 70 |BAND *)N|ULL || p|
|00005ba0| 69 76 6f 74 3d 3d 28 50 | 45 52 4d 20 2a 29 4e 55 |ivot==(P|ERM *)NU|
|00005bb0| 4c 4c 20 29 0a 20 20 20 | 20 20 65 72 72 6f 72 28 |LL ). | error(|
|00005bc0| 45 5f 4e 55 4c 4c 2c 22 | 62 64 4c 55 66 61 63 74 |E_NULL,"|bdLUfact|
|00005bd0| 6f 72 22 29 3b 0a 0a 20 | 20 20 6c 62 20 3d 20 62 |or");.. | lb = b|
|00005be0| 41 2d 3e 6c 62 3b 0a 20 | 20 20 75 62 20 3d 20 62 |A->lb;. | ub = b|
|00005bf0| 41 2d 3e 75 62 3b 0a 20 | 20 20 6c 75 62 20 3d 20 |A->ub;. | lub = |
|00005c00| 6c 62 2b 75 62 3b 0a 20 | 20 20 6e 20 3d 20 62 41 |lb+ub;. | n = bA|
|00005c10| 2d 3e 6d 61 74 2d 3e 6e | 3b 0a 20 20 20 6e 31 20 |->mat->n|;. n1 |
|00005c20| 3d 20 6e 2d 31 3b 0a 20 | 20 20 6c 75 62 20 3d 20 |= n-1;. | lub = |
|00005c30| 6c 62 2b 75 62 3b 0a 0a | 20 20 20 69 66 20 28 20 |lb+ub;..| if ( |
|00005c40| 70 69 76 6f 74 2d 3e 73 | 69 7a 65 20 21 3d 20 6e |pivot->s|ize != n|
|00005c50| 20 29 0a 20 20 20 20 20 | 65 72 72 6f 72 28 45 5f | ). |error(E_|
|00005c60| 53 49 5a 45 53 2c 22 62 | 64 4c 55 66 61 63 74 6f |SIZES,"b|dLUfacto|
|00005c70| 72 22 29 3b 0a 0a 20 20 | 20 0a 20 20 20 2f 2a 20 |r");.. | . /* |
|00005c80| 69 6e 69 74 69 61 6c 69 | 73 65 20 70 69 76 6f 74 |initiali|se pivot|
|00005c90| 20 77 69 74 68 20 69 64 | 65 6e 74 69 74 79 20 70 | with id|entity p|
|00005ca0| 65 72 6d 75 74 61 74 69 | 6f 6e 20 2a 2f 0a 20 20 |ermutati|on */. |
|00005cb0| 20 66 6f 72 20 28 20 69 | 3d 30 3b 20 69 20 3c 20 | for ( i|=0; i < |
|00005cc0| 6e 3b 20 69 2b 2b 20 29 | 0a 20 20 20 20 20 70 69 |n; i++ )|. pi|
|00005cd0| 76 6f 74 2d 3e 70 65 5b | 69 5d 20 3d 20 69 3b 0a |vot->pe[|i] = i;.|
|00005ce0| 0a 20 20 20 2f 2a 20 65 | 78 74 65 6e 64 20 62 61 |. /* e|xtend ba|
|00005cf0| 6e 64 20 6d 61 74 72 69 | 78 20 2a 2f 0a 20 20 20 |nd matri|x */. |
|00005d00| 2f 2a 20 65 78 74 65 6e | 64 65 64 20 70 61 72 74 |/* exten|ded part|
|00005d10| 20 69 73 20 66 69 6c 6c | 65 64 20 77 69 74 68 20 | is fill|ed with |
|00005d20| 7a 65 72 6f 73 20 2a 2f | 0a 20 20 20 62 41 20 3d |zeros */|. bA =|
|00005d30| 20 62 64 5f 72 65 73 69 | 7a 65 28 62 41 2c 6c 62 | bd_resi|ze(bA,lb|
|00005d40| 2c 6d 69 6e 28 6e 31 2c | 6c 75 62 29 2c 6e 29 3b |,min(n1,|lub),n);|
|00005d50| 0a 20 20 20 62 41 5f 76 | 20 3d 20 62 41 2d 3e 6d |. bA_v| = bA->m|
|00005d60| 61 74 2d 3e 6d 65 3b 0a | 0a 0a 20 20 20 2f 2a 20 |at->me;.|.. /* |
|00005d70| 6d 61 69 6e 20 6c 6f 6f | 70 20 2a 2f 0a 0a 20 20 |main loo|p */.. |
|00005d80| 20 66 6f 72 20 28 20 6b | 3d 30 3b 20 6b 20 3c 20 | for ( k|=0; k < |
|00005d90| 6e 31 3b 20 6b 2b 2b 20 | 29 0a 20 20 20 7b 0a 20 |n1; k++ |). {. |
|00005da0| 20 20 20 20 20 6b 5f 65 | 6e 64 20 3d 20 6d 61 78 | k_e|nd = max|
|00005db0| 28 30 2c 6c 62 2b 6b 2d | 6e 31 29 3b 0a 20 20 20 |(0,lb+k-|n1);. |
|00005dc0| 20 20 20 6b 5f 6c 75 62 | 20 3d 20 6d 69 6e 28 6b | k_lub| = min(k|
|00005dd0| 2b 6c 75 62 2c 6e 31 29 | 3b 0a 0a 20 20 20 20 20 |+lub,n1)|;.. |
|00005de0| 20 2f 2a 20 66 69 6e 64 | 20 74 68 65 20 62 65 73 | /* find| the bes|
|00005df0| 74 20 70 69 76 6f 74 20 | 72 6f 77 20 2a 2f 0a 20 |t pivot |row */. |
|00005e00| 20 20 20 20 20 0a 20 20 | 20 20 20 20 6d 61 78 31 | . | max1|
|00005e10| 20 3d 20 30 2e 30 3b 09 | 0a 20 20 20 20 20 20 69 | = 0.0;.|. i|
|00005e20| 5f 6d 61 78 20 3d 20 2d | 31 3b 0a 20 20 20 20 20 |_max = -|1;. |
|00005e30| 20 66 6f 72 20 28 20 69 | 3d 6c 62 3b 20 69 20 3e | for ( i|=lb; i >|
|00005e40| 3d 20 6b 5f 65 6e 64 3b | 20 69 2d 2d 20 29 20 7b |= k_end;| i-- ) {|
|00005e50| 0a 09 20 74 65 6d 70 20 | 3d 20 66 61 62 73 28 62 |.. temp |= fabs(b|
|00005e60| 41 5f 76 5b 69 5d 5b 6b | 5d 29 3b 0a 09 20 69 66 |A_v[i][k|]);.. if|
|00005e70| 20 28 20 74 65 6d 70 20 | 3e 20 6d 61 78 31 20 29 | ( temp |> max1 )|
|00005e80| 0a 09 20 7b 20 6d 61 78 | 31 20 3d 20 74 65 6d 70 |.. { max|1 = temp|
|00005e90| 3b 09 69 5f 6d 61 78 20 | 3d 20 69 3b 20 7d 0a 20 |;.i_max |= i; }. |
|00005ea0| 20 20 20 20 20 7d 0a 20 | 20 20 20 20 20 0a 20 20 | }. | . |
|00005eb0| 20 20 20 20 2f 2a 20 69 | 66 20 6e 6f 20 70 69 76 | /* i|f no piv|
|00005ec0| 6f 74 20 74 68 65 6e 20 | 69 67 6e 6f 72 65 20 63 |ot then |ignore c|
|00005ed0| 6f 6c 75 6d 6e 20 6b 2e | 2e 2e 20 2a 2f 0a 20 20 |olumn k.|.. */. |
|00005ee0| 20 20 20 20 69 66 20 28 | 20 69 5f 6d 61 78 20 3d | if (| i_max =|
|00005ef0| 3d 20 2d 31 20 29 0a 09 | 63 6f 6e 74 69 6e 75 65 |= -1 )..|continue|
|00005f00| 3b 0a 20 20 20 20 20 20 | 0a 20 20 20 20 20 20 2f |;. |. /|
|00005f10| 2a 20 64 6f 20 77 65 20 | 70 69 76 6f 74 20 3f 20 |* do we |pivot ? |
|00005f20| 2a 2f 0a 20 20 20 20 20 | 20 69 66 20 28 20 69 5f |*/. | if ( i_|
|00005f30| 6d 61 78 20 21 3d 20 6c | 62 20 29 09 2f 2a 20 79 |max != l|b )./* y|
|00005f40| 65 73 20 77 65 20 64 6f | 2e 2e 2e 20 2a 2f 0a 20 |es we do|... */. |
|00005f50| 20 20 20 20 20 7b 0a 09 | 20 2f 2a 20 73 61 76 65 | {..| /* save|
|00005f60| 20 74 72 61 6e 73 70 6f | 73 69 74 69 6f 6e 20 75 | transpo|sition u|
|00005f70| 73 69 6e 67 20 6e 6f 6e | 2d 73 68 69 66 74 65 64 |sing non|-shifted|
|00005f80| 20 69 6e 64 69 63 65 73 | 20 2a 2f 0a 09 20 73 68 | indices| */.. sh|
|00005f90| 69 66 74 20 3d 20 6c 62 | 2d 69 5f 6d 61 78 3b 0a |ift = lb|-i_max;.|
|00005fa0| 09 20 70 78 5f 74 72 61 | 6e 73 70 28 70 69 76 6f |. px_tra|nsp(pivo|
|00005fb0| 74 2c 6b 2b 73 68 69 66 | 74 2c 6b 29 3b 0a 09 20 |t,k+shif|t,k);.. |
|00005fc0| 66 6f 72 20 28 20 69 3d | 6c 62 2c 20 6a 3d 6b 3b |for ( i=|lb, j=k;|
|00005fd0| 20 6a 20 3c 3d 20 6b 5f | 6c 75 62 3b 20 69 2b 2b | j <= k_|lub; i++|
|00005fe0| 2c 6a 2b 2b 20 29 0a 09 | 20 7b 0a 09 20 20 20 20 |,j++ )..| {.. |
|00005ff0| 74 65 6d 70 20 3d 20 62 | 41 5f 76 5b 69 5d 5b 6a |temp = b|A_v[i][j|
|00006000| 5d 3b 0a 09 20 20 20 20 | 62 41 5f 76 5b 69 5d 5b |];.. |bA_v[i][|
|00006010| 6a 5d 20 3d 20 62 41 5f | 76 5b 69 2d 73 68 69 66 |j] = bA_|v[i-shif|
|00006020| 74 5d 5b 6a 5d 3b 0a 09 | 20 20 20 20 62 41 5f 76 |t][j];..| bA_v|
|00006030| 5b 69 2d 73 68 69 66 74 | 5d 5b 6a 5d 20 3d 20 74 |[i-shift|][j] = t|
|00006040| 65 6d 70 3b 0a 09 20 7d | 0a 20 20 20 20 20 20 7d |emp;.. }|. }|
|00006050| 0a 20 20 20 20 20 20 0a | 20 20 20 20 20 20 2f 2a |. .| /*|
|00006060| 20 72 6f 77 20 6f 70 65 | 72 61 74 69 6f 6e 73 20 | row ope|rations |
|00006070| 2a 2f 0a 20 20 20 20 20 | 20 66 6f 72 20 28 20 69 |*/. | for ( i|
|00006080| 3d 6c 62 2d 31 3b 20 69 | 20 3e 3d 20 6b 5f 65 6e |=lb-1; i| >= k_en|
|00006090| 64 3b 20 69 2d 2d 20 29 | 20 7b 0a 09 20 74 65 6d |d; i-- )| {.. tem|
|000060a0| 70 20 3d 20 62 41 5f 76 | 5b 69 5d 5b 6b 5d 20 2f |p = bA_v|[i][k] /|
|000060b0| 3d 20 62 41 5f 76 5b 6c | 62 5d 5b 6b 5d 3b 0a 09 |= bA_v[l|b][k];..|
|000060c0| 20 73 68 69 66 74 20 3d | 20 6c 62 2d 69 3b 0a 09 | shift =| lb-i;..|
|000060d0| 20 66 6f 72 20 28 20 6a | 3d 6b 2b 31 2c 6c 3d 69 | for ( j|=k+1,l=i|
|000060e0| 2b 31 3b 20 6a 20 3c 3d | 20 6b 5f 6c 75 62 3b 20 |+1; j <=| k_lub; |
|000060f0| 6c 2b 2b 2c 6a 2b 2b 20 | 29 0a 09 20 20 20 62 41 |l++,j++ |).. bA|
|00006100| 5f 76 5b 6c 5d 5b 6a 5d | 20 2d 3d 20 74 65 6d 70 |_v[l][j]| -= temp|
|00006110| 2a 62 41 5f 76 5b 6c 2b | 73 68 69 66 74 5d 5b 6a |*bA_v[l+|shift][j|
|00006120| 5d 3b 0a 20 20 20 20 20 | 20 7d 0a 20 20 20 7d 0a |];. | }. }.|
|00006130| 20 20 20 0a 20 20 20 72 | 65 74 75 72 6e 20 62 41 | . r|eturn bA|
|00006140| 3b 0a 7d 0a 0a 0a 2f 2a | 20 62 64 4c 55 73 6f 6c |;.}.../*| bdLUsol|
|00006150| 76 65 20 2d 2d 20 67 69 | 76 65 6e 20 61 6e 20 4c |ve -- gi|ven an L|
|00006160| 55 20 66 61 63 74 6f 72 | 69 73 61 74 69 6f 6e 20 |U factor|isation |
|00006170| 69 6e 20 62 41 2c 20 73 | 6f 6c 76 65 20 62 41 2a |in bA, s|olve bA*|
|00006180| 78 3d 62 20 2a 2f 0a 2f | 2a 20 70 69 76 6f 74 20 |x=b */./|* pivot |
|00006190| 69 73 20 63 68 61 6e 67 | 65 64 20 75 70 6f 6e 20 |is chang|ed upon |
|000061a0| 72 65 74 75 72 6e 20 20 | 2a 2f 0a 56 45 43 09 2a |return |*/.VEC.*|
|000061b0| 62 64 4c 55 73 6f 6c 76 | 65 28 62 41 2c 70 69 76 |bdLUsolv|e(bA,piv|
|000061c0| 6f 74 2c 62 2c 78 29 0a | 42 41 4e 44 09 2a 62 41 |ot,b,x).|BAND.*bA|
|000061d0| 3b 0a 50 45 52 4d 09 2a | 70 69 76 6f 74 3b 0a 56 |;.PERM.*|pivot;.V|
|000061e0| 45 43 09 2a 62 2c 2a 78 | 3b 0a 7b 0a 20 20 20 69 |EC.*b,*x|;.{. i|
|000061f0| 6e 74 20 69 2c 6a 2c 6c | 2c 6e 2c 6e 31 2c 70 69 |nt i,j,l|,n,n1,pi|
|00006200| 2c 6c 62 2c 75 62 2c 6a | 6d 69 6e 2c 20 6d 61 78 |,lb,ub,j|min, max|
|00006210| 6a 3b 0a 20 20 20 52 65 | 61 6c 20 63 3b 0a 20 20 |j;. Re|al c;. |
|00006220| 20 52 65 61 6c 20 2a 2a | 62 41 5f 76 3b 0a 0a 20 | Real **|bA_v;.. |
|00006230| 20 20 69 66 20 28 20 62 | 41 3d 3d 28 42 41 4e 44 | if ( b|A==(BAND|
|00006240| 20 2a 29 4e 55 4c 4c 20 | 7c 7c 20 62 3d 3d 28 56 | *)NULL ||| b==(V|
|00006250| 45 43 20 2a 29 4e 55 4c | 4c 20 7c 7c 20 70 69 76 |EC *)NUL|L || piv|
|00006260| 6f 74 3d 3d 28 50 45 52 | 4d 20 2a 29 4e 55 4c 4c |ot==(PER|M *)NULL|
|00006270| 20 29 0a 20 20 20 20 20 | 65 72 72 6f 72 28 45 5f | ). |error(E_|
|00006280| 4e 55 4c 4c 2c 22 62 64 | 4c 55 73 6f 6c 76 65 22 |NULL,"bd|LUsolve"|
|00006290| 29 3b 0a 20 20 20 69 66 | 20 28 20 62 41 2d 3e 6d |);. if| ( bA->m|
|000062a0| 61 74 2d 3e 6e 20 21 3d | 20 62 2d 3e 64 69 6d 20 |at->n !=| b->dim |
|000062b0| 7c 7c 20 62 41 2d 3e 6d | 61 74 2d 3e 6e 20 21 3d ||| bA->m|at->n !=|
|000062c0| 20 70 69 76 6f 74 2d 3e | 73 69 7a 65 29 0a 20 20 | pivot->|size). |
|000062d0| 20 20 20 65 72 72 6f 72 | 28 45 5f 53 49 5a 45 53 | error|(E_SIZES|
|000062e0| 2c 22 62 64 4c 55 73 6f | 6c 76 65 22 29 3b 0a 20 |,"bdLUso|lve");. |
|000062f0| 0a 20 20 20 6c 62 20 3d | 20 62 41 2d 3e 6c 62 3b |. lb =| bA->lb;|
|00006300| 0a 20 20 20 75 62 20 3d | 20 62 41 2d 3e 75 62 3b |. ub =| bA->ub;|
|00006310| 0a 20 20 20 6e 20 3d 20 | 62 2d 3e 64 69 6d 3b 0a |. n = |b->dim;.|
|00006320| 20 20 20 6e 31 20 3d 20 | 6e 2d 31 3b 0a 20 20 20 | n1 = |n-1;. |
|00006330| 62 41 5f 76 20 3d 20 62 | 41 2d 3e 6d 61 74 2d 3e |bA_v = b|A->mat->|
|00006340| 6d 65 3b 0a 0a 20 20 20 | 78 20 3d 20 76 5f 72 65 |me;.. |x = v_re|
|00006350| 73 69 7a 65 28 78 2c 62 | 2d 3e 64 69 6d 29 3b 0a |size(x,b|->dim);.|
|00006360| 20 20 20 70 78 5f 76 65 | 63 28 70 69 76 6f 74 2c | px_ve|c(pivot,|
|00006370| 62 2c 78 29 3b 0a 0a 20 | 20 20 2f 2a 20 73 6f 6c |b,x);.. | /* sol|
|00006380| 76 65 20 4c 78 20 3d 20 | 62 3b 20 69 6d 70 6c 69 |ve Lx = |b; impli|
|00006390| 63 69 74 20 64 69 61 67 | 6f 6e 61 6c 20 3d 20 31 |cit diag|onal = 1|
|000063a0| 20 0a 20 20 20 20 20 20 | 4c 20 69 73 20 6e 6f 74 | . |L is not|
|000063b0| 20 70 65 72 6d 75 74 65 | 64 2c 20 74 68 65 72 65 | permute|d, there|
|000063c0| 66 6f 72 65 20 69 74 20 | 6d 75 73 74 20 62 65 20 |fore it |must be |
|000063d0| 70 65 72 6d 75 74 65 64 | 20 6e 6f 77 0a 20 20 20 |permuted| now. |
|000063e0| 20 2a 2f 0a 20 20 20 0a | 20 20 20 70 78 5f 69 6e | */. .| px_in|
|000063f0| 76 28 70 69 76 6f 74 2c | 70 69 76 6f 74 29 3b 0a |v(pivot,|pivot);.|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
