home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: SHell self-extracting ARchive
(archive/shar).
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| SHell self-extracting ARchive (archive/shar)
| magic
| Supported |
| 1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
| 100%
| file
| ASCII text
| default
| |
| 100%
| checkBytes
| Printable ASCII
| default
| |
| 100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
| 100%
| siegfried
| fmt/329 Shell Archive Format
| default
| |
| 100%
| detectItEasy
| Format: plain text[LF]
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 6d 69 73 63 0a 50 61 |.sources|.misc.Pa|
|00000020| 74 68 3a 20 73 70 61 72 | 6b 79 21 6b 65 6e 74 0a |th: spar|ky!kent.|
|00000030| 46 72 6f 6d 3a 20 72 6f | 62 65 72 74 64 40 6b 61 |From: ro|bertd@ka|
|00000040| 75 72 69 2e 76 75 77 2e | 61 63 2e 6e 7a 20 28 52 |uri.vuw.|ac.nz (R|
|00000050| 6f 62 65 72 74 20 44 61 | 76 69 65 73 29 0a 53 75 |obert Da|vies).Su|
|00000060| 62 6a 65 63 74 3a 20 76 | 33 34 69 31 31 31 3a 20 |bject: v|34i111: |
|00000070| 20 6e 65 77 6d 61 74 30 | 37 20 2d 20 41 20 6d 61 | newmat0|7 - A ma|
|00000080| 74 72 69 78 20 70 61 63 | 6b 61 67 65 20 69 6e 20 |trix pac|kage in |
|00000090| 43 2b 2b 2c 20 50 61 72 | 74 30 35 2f 30 38 0a 4d |C++, Par|t05/08.M|
|000000a0| 65 73 73 61 67 65 2d 49 | 44 3a 20 3c 31 39 39 33 |essage-I|D: <1993|
|000000b0| 4a 61 6e 31 31 2e 31 35 | 33 32 31 35 2e 32 34 36 |Jan11.15|3215.246|
|000000c0| 30 40 73 70 61 72 6b 79 | 2e 69 6d 64 2e 73 74 65 |0@sparky|.imd.ste|
|000000d0| 72 6c 69 6e 67 2e 63 6f | 6d 3e 0a 46 6f 6c 6c 6f |rling.co|m>.Follo|
|000000e0| 77 75 70 2d 54 6f 3a 20 | 63 6f 6d 70 2e 73 6f 75 |wup-To: |comp.sou|
|000000f0| 72 63 65 73 2e 64 0a 58 | 2d 4d 64 34 2d 53 69 67 |rces.d.X|-Md4-Sig|
|00000100| 6e 61 74 75 72 65 3a 20 | 39 35 63 33 62 66 65 33 |nature: |95c3bfe3|
|00000110| 65 36 33 34 64 63 30 33 | 61 62 31 62 35 34 36 63 |e634dc03|ab1b546c|
|00000120| 63 65 65 38 31 33 37 61 | 0a 53 65 6e 64 65 72 3a |cee8137a|.Sender:|
|00000130| 20 6b 65 6e 74 40 73 70 | 61 72 6b 79 2e 69 6d 64 | kent@sp|arky.imd|
|00000140| 2e 73 74 65 72 6c 69 6e | 67 2e 63 6f 6d 20 28 4b |.sterlin|g.com (K|
|00000150| 65 6e 74 20 4c 61 6e 64 | 66 69 65 6c 64 29 0a 4f |ent Land|field).O|
|00000160| 72 67 61 6e 69 7a 61 74 | 69 6f 6e 3a 20 53 74 65 |rganizat|ion: Ste|
|00000170| 72 6c 69 6e 67 20 53 6f | 66 74 77 61 72 65 0a 52 |rling So|ftware.R|
|00000180| 65 66 65 72 65 6e 63 65 | 73 3a 20 3c 63 73 6d 2d |eference|s: <csm-|
|00000190| 76 33 34 69 31 30 37 3d | 6e 65 77 6d 61 74 30 37 |v34i107=|newmat07|
|000001a0| 2e 30 39 32 37 30 36 40 | 73 70 61 72 6b 79 2e 49 |.092706@|sparky.I|
|000001b0| 4d 44 2e 53 74 65 72 6c | 69 6e 67 2e 43 4f 4d 3e |MD.Sterl|ing.COM>|
|000001c0| 0a 44 61 74 65 3a 20 4d | 6f 6e 2c 20 31 31 20 4a |.Date: M|on, 11 J|
|000001d0| 61 6e 20 31 39 39 33 20 | 31 35 3a 33 32 3a 31 35 |an 1993 |15:32:15|
|000001e0| 20 47 4d 54 0a 41 70 70 | 72 6f 76 65 64 3a 20 6b | GMT.App|roved: k|
|000001f0| 65 6e 74 40 73 70 61 72 | 6b 79 2e 69 6d 64 2e 73 |ent@spar|ky.imd.s|
|00000200| 74 65 72 6c 69 6e 67 2e | 63 6f 6d 0a 4c 69 6e 65 |terling.|com.Line|
|00000210| 73 3a 20 32 31 36 34 0a | 0a 53 75 62 6d 69 74 74 |s: 2164.|.Submitt|
|00000220| 65 64 2d 62 79 3a 20 72 | 6f 62 65 72 74 64 40 6b |ed-by: r|obertd@k|
|00000230| 61 75 72 69 2e 76 75 77 | 2e 61 63 2e 6e 7a 20 28 |auri.vuw|.ac.nz (|
|00000240| 52 6f 62 65 72 74 20 44 | 61 76 69 65 73 29 0a 50 |Robert D|avies).P|
|00000250| 6f 73 74 69 6e 67 2d 6e | 75 6d 62 65 72 3a 20 56 |osting-n|umber: V|
|00000260| 6f 6c 75 6d 65 20 33 34 | 2c 20 49 73 73 75 65 20 |olume 34|, Issue |
|00000270| 31 31 31 0a 41 72 63 68 | 69 76 65 2d 6e 61 6d 65 |111.Arch|ive-name|
|00000280| 3a 20 6e 65 77 6d 61 74 | 30 37 2f 70 61 72 74 30 |: newmat|07/part0|
|00000290| 35 0a 45 6e 76 69 72 6f | 6e 6d 65 6e 74 3a 20 43 |5.Enviro|nment: C|
|000002a0| 2b 2b 0a 53 75 70 65 72 | 73 65 64 65 73 3a 20 6e |++.Super|sedes: n|
|000002b0| 65 77 6d 61 74 30 36 3a | 20 56 6f 6c 75 6d 65 20 |ewmat06:| Volume |
|000002c0| 33 34 2c 20 49 73 73 75 | 65 20 37 2d 31 33 0a 0a |34, Issu|e 7-13..|
|000002d0| 23 21 20 2f 62 69 6e 2f | 73 68 0a 23 20 54 68 69 |#! /bin/|sh.# Thi|
|000002e0| 73 20 69 73 20 61 20 73 | 68 65 6c 6c 20 61 72 63 |s is a s|hell arc|
|000002f0| 68 69 76 65 2e 20 20 52 | 65 6d 6f 76 65 20 61 6e |hive. R|emove an|
|00000300| 79 74 68 69 6e 67 20 62 | 65 66 6f 72 65 20 74 68 |ything b|efore th|
|00000310| 69 73 20 6c 69 6e 65 2c | 20 74 68 65 6e 20 75 6e |is line,| then un|
|00000320| 70 61 63 6b 0a 23 20 69 | 74 20 62 79 20 73 61 76 |pack.# i|t by sav|
|00000330| 69 6e 67 20 69 74 20 69 | 6e 74 6f 20 61 20 66 69 |ing it i|nto a fi|
|00000340| 6c 65 20 61 6e 64 20 74 | 79 70 69 6e 67 20 22 73 |le and t|yping "s|
|00000350| 68 20 66 69 6c 65 22 2e | 20 20 54 6f 20 6f 76 65 |h file".| To ove|
|00000360| 72 77 72 69 74 65 20 65 | 78 69 73 74 69 6e 67 0a |rwrite e|xisting.|
|00000370| 23 20 66 69 6c 65 73 2c | 20 74 79 70 65 20 22 73 |# files,| type "s|
|00000380| 68 20 66 69 6c 65 20 2d | 63 22 2e 20 20 59 6f 75 |h file -|c". You|
|00000390| 20 63 61 6e 20 61 6c 73 | 6f 20 66 65 65 64 20 74 | can als|o feed t|
|000003a0| 68 69 73 20 61 73 20 73 | 74 61 6e 64 61 72 64 20 |his as s|tandard |
|000003b0| 69 6e 70 75 74 20 76 69 | 61 0a 23 20 75 6e 73 68 |input vi|a.# unsh|
|000003c0| 61 72 2c 20 6f 72 20 62 | 79 20 74 79 70 69 6e 67 |ar, or b|y typing|
|000003d0| 20 22 73 68 20 3c 66 69 | 6c 65 22 2c 20 65 2e 67 | "sh <fi|le", e.g|
|000003e0| 2e 2e 20 20 49 66 20 74 | 68 69 73 20 61 72 63 68 |.. If t|his arch|
|000003f0| 69 76 65 20 69 73 20 63 | 6f 6d 70 6c 65 74 65 2c |ive is c|omplete,|
|00000400| 20 79 6f 75 0a 23 20 77 | 69 6c 6c 20 73 65 65 20 | you.# w|ill see |
|00000410| 74 68 65 20 66 6f 6c 6c | 6f 77 69 6e 67 20 6d 65 |the foll|owing me|
|00000420| 73 73 61 67 65 20 61 74 | 20 74 68 65 20 65 6e 64 |ssage at| the end|
|00000430| 3a 0a 23 09 09 22 45 6e | 64 20 6f 66 20 61 72 63 |:.#.."En|d of arc|
|00000440| 68 69 76 65 20 35 20 28 | 6f 66 20 38 29 2e 22 0a |hive 5 (|of 8).".|
|00000450| 23 20 43 6f 6e 74 65 6e | 74 73 3a 20 20 6e 65 77 |# Conten|ts: new|
|00000460| 6d 61 74 35 2e 63 78 78 | 20 6e 65 77 6d 61 74 36 |mat5.cxx| newmat6|
|00000470| 2e 63 78 78 20 6e 65 77 | 6d 61 74 37 2e 63 78 78 |.cxx new|mat7.cxx|
|00000480| 20 6e 65 77 6d 61 74 38 | 2e 63 78 78 0a 23 20 20 | newmat8|.cxx.# |
|00000490| 20 6e 65 77 6d 61 74 39 | 2e 63 78 78 20 6e 65 77 | newmat9|.cxx new|
|000004a0| 6d 61 74 72 6d 2e 63 78 | 78 0a 23 20 57 72 61 70 |matrm.cx|x.# Wrap|
|000004b0| 70 65 64 20 62 79 20 72 | 6f 62 65 72 74 40 6b 65 |ped by r|obert@ke|
|000004c0| 61 20 6f 6e 20 53 75 6e | 20 4a 61 6e 20 31 30 20 |a on Sun| Jan 10 |
|000004d0| 32 33 3a 35 38 3a 30 37 | 20 31 39 39 33 0a 50 41 |23:58:07| 1993.PA|
|000004e0| 54 48 3d 2f 62 69 6e 3a | 2f 75 73 72 2f 62 69 6e |TH=/bin:|/usr/bin|
|000004f0| 3a 2f 75 73 72 2f 75 63 | 62 20 3b 20 65 78 70 6f |:/usr/uc|b ; expo|
|00000500| 72 74 20 50 41 54 48 0a | 69 66 20 74 65 73 74 20 |rt PATH.|if test |
|00000510| 2d 66 20 27 6e 65 77 6d | 61 74 35 2e 63 78 78 27 |-f 'newm|at5.cxx'|
|00000520| 20 2d 61 20 22 24 7b 31 | 7d 22 20 21 3d 20 22 2d | -a "${1|}" != "-|
|00000530| 63 22 20 3b 20 74 68 65 | 6e 20 0a 20 20 65 63 68 |c" ; the|n . ech|
|00000540| 6f 20 73 68 61 72 3a 20 | 57 69 6c 6c 20 6e 6f 74 |o shar: |Will not|
|00000550| 20 63 6c 6f 62 62 65 72 | 20 65 78 69 73 74 69 6e | clobber| existin|
|00000560| 67 20 66 69 6c 65 20 5c | 22 27 6e 65 77 6d 61 74 |g file \|"'newmat|
|00000570| 35 2e 63 78 78 27 5c 22 | 0a 65 6c 73 65 0a 65 63 |5.cxx'\"|.else.ec|
|00000580| 68 6f 20 73 68 61 72 3a | 20 45 78 74 72 61 63 74 |ho shar:| Extract|
|00000590| 69 6e 67 20 5c 22 27 6e | 65 77 6d 61 74 35 2e 63 |ing \"'n|ewmat5.c|
|000005a0| 78 78 27 5c 22 20 5c 28 | 31 31 35 36 37 20 63 68 |xx'\" \(|11567 ch|
|000005b0| 61 72 61 63 74 65 72 73 | 5c 29 0a 73 65 64 20 22 |aracters|\).sed "|
|000005c0| 73 2f 5e 58 2f 2f 22 20 | 3e 27 6e 65 77 6d 61 74 |s/^X//" |>'newmat|
|000005d0| 35 2e 63 78 78 27 20 3c | 3c 27 45 4e 44 5f 4f 46 |5.cxx' <|<'END_OF|
|000005e0| 5f 46 49 4c 45 27 0a 58 | 2f 2f 24 24 20 6e 65 77 |_FILE'.X|//$$ new|
|000005f0| 6d 61 74 35 2e 63 78 78 | 20 20 20 20 20 20 20 20 |mat5.cxx| |
|00000600| 20 54 72 61 6e 73 70 6f | 73 65 2c 20 65 76 61 6c | Transpo|se, eval|
|00000610| 75 61 74 65 20 65 74 63 | 0a 58 0a 58 2f 2f 20 43 |uate etc|.X.X// C|
|00000620| 6f 70 79 72 69 67 68 74 | 20 28 43 29 20 31 39 39 |opyright| (C) 199|
|00000630| 31 2c 32 2c 33 3a 20 52 | 20 42 20 44 61 76 69 65 |1,2,3: R| B Davie|
|00000640| 73 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 20 22 69 |s.X.X#in|clude "i|
|00000650| 6e 63 6c 75 64 65 2e 68 | 22 0a 58 0a 58 23 69 6e |nclude.h|".X.X#in|
|00000660| 63 6c 75 64 65 20 22 6e | 65 77 6d 61 74 2e 68 22 |clude "n|ewmat.h"|
|00000670| 0a 58 23 69 6e 63 6c 75 | 64 65 20 22 6e 65 77 6d |.X#inclu|de "newm|
|00000680| 61 74 72 63 2e 68 22 0a | 58 0a 58 2f 2f 23 64 65 |atrc.h".|X.X//#de|
|00000690| 66 69 6e 65 20 52 45 50 | 4f 52 54 20 7b 20 73 74 |fine REP|ORT { st|
|000006a0| 61 74 69 63 20 45 78 65 | 43 6f 75 6e 74 65 72 20 |atic Exe|Counter |
|000006b0| 45 78 65 43 6f 75 6e 74 | 28 5f 5f 4c 49 4e 45 5f |ExeCount|(__LINE_|
|000006c0| 5f 2c 35 29 3b 20 45 78 | 65 43 6f 75 6e 74 2b 2b |_,5); Ex|eCount++|
|000006d0| 3b 20 7d 0a 58 0a 58 23 | 64 65 66 69 6e 65 20 52 |; }.X.X#|define R|
|000006e0| 45 50 4f 52 54 20 7b 7d | 0a 58 0a 58 0a 58 2f 2a |EPORT {}|.X.X.X/*|
|000006f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000700| 2a 2a 2a 2a 2a 2a 2a 20 | 63 61 72 72 79 20 6f 75 |******* |carry ou|
|00000710| 74 20 6f 70 65 72 61 74 | 69 6f 6e 73 20 2a 2a 2a |t operat|ions ***|
|00000720| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000730| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2f 0a 58 0a 58 |********|***/.X.X|
|00000740| 0a 58 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |.XGenera|lMatrix*|
|00000750| 20 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 3a 3a | General|Matrix::|
|00000760| 54 72 61 6e 73 70 6f 73 | 65 28 54 72 61 6e 73 70 |Transpos|e(Transp|
|00000770| 6f 73 65 64 4d 61 74 72 | 69 78 2a 20 74 6d 2c 20 |osedMatr|ix* tm, |
|00000780| 4d 61 74 72 69 78 54 79 | 70 65 20 6d 74 29 0a 58 |MatrixTy|pe mt).X|
|00000790| 7b 0a 58 20 20 20 47 65 | 6e 65 72 61 6c 4d 61 74 |{.X Ge|neralMat|
|000007a0| 72 69 78 2a 20 67 6d 31 | 3b 0a 58 0a 58 20 20 20 |rix* gm1|;.X.X |
|000007b0| 69 66 20 28 43 6f 6d 70 | 61 72 65 28 54 79 70 65 |if (Comp|are(Type|
|000007c0| 28 29 2e 74 28 29 2c 6d | 74 29 29 0a 58 20 20 20 |().t(),m|t)).X |
|000007d0| 7b 0a 58 20 20 20 20 20 | 20 52 45 50 4f 52 54 0a |{.X | REPORT.|
|000007e0| 58 20 20 20 20 20 20 67 | 6d 31 20 3d 20 6d 74 2e |X g|m1 = mt.|
|000007f0| 4e 65 77 28 6e 63 6f 6c | 73 2c 6e 72 6f 77 73 2c |New(ncol|s,nrows,|
|00000800| 74 6d 29 3b 0a 58 20 20 | 20 20 20 20 66 6f 72 20 |tm);.X | for |
|00000810| 28 69 6e 74 20 69 3d 30 | 3b 20 69 3c 6e 63 6f 6c |(int i=0|; i<ncol|
|00000820| 73 3b 20 69 2b 2b 29 0a | 58 20 20 20 20 20 20 7b |s; i++).|X {|
|00000830| 0a 58 09 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 |.X. Matr|ixRow mr|
|00000840| 28 67 6d 31 2c 20 53 74 | 6f 72 65 4f 6e 45 78 69 |(gm1, St|oreOnExi|
|00000850| 74 2b 44 69 72 65 63 74 | 50 61 72 74 2c 20 69 29 |t+Direct|Part, i)|
|00000860| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 4d 61 74 72 |;.X | Matr|
|00000870| 69 78 43 6f 6c 20 6d 63 | 28 74 68 69 73 2c 20 6d |ixCol mc|(this, m|
|00000880| 72 2e 53 74 6f 72 65 28 | 29 2c 20 4c 6f 61 64 4f |r.Store(|), LoadO|
|00000890| 6e 45 6e 74 72 79 2c 20 | 69 29 3b 0a 58 20 20 20 |nEntry, |i);.X |
|000008a0| 20 20 20 7d 0a 58 20 20 | 20 7d 0a 58 20 20 20 65 | }.X | }.X e|
|000008b0| 6c 73 65 0a 58 20 20 20 | 7b 0a 58 20 20 20 20 20 |lse.X |{.X |
|000008c0| 20 52 45 50 4f 52 54 0a | 58 20 20 20 20 20 20 67 | REPORT.|X g|
|000008d0| 6d 31 20 3d 20 6d 74 2e | 4e 65 77 28 6e 63 6f 6c |m1 = mt.|New(ncol|
|000008e0| 73 2c 6e 72 6f 77 73 2c | 74 6d 29 3b 0a 58 20 20 |s,nrows,|tm);.X |
|000008f0| 20 20 20 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 | Matr|ixRow mr|
|00000900| 28 74 68 69 73 2c 20 4c | 6f 61 64 4f 6e 45 6e 74 |(this, L|oadOnEnt|
|00000910| 72 79 29 3b 0a 58 20 20 | 20 20 20 20 4d 61 74 72 |ry);.X | Matr|
|00000920| 69 78 43 6f 6c 20 6d 63 | 28 67 6d 31 2c 20 53 74 |ixCol mc|(gm1, St|
|00000930| 6f 72 65 4f 6e 45 78 69 | 74 2b 44 69 72 65 63 74 |oreOnExi|t+Direct|
|00000940| 50 61 72 74 29 3b 0a 58 | 20 20 20 20 20 20 69 6e |Part);.X| in|
|00000950| 74 20 69 20 3d 20 6e 72 | 6f 77 73 3b 0a 58 20 20 |t i = nr|ows;.X |
|00000960| 20 20 20 20 77 68 69 6c | 65 20 28 69 2d 2d 29 20 | whil|e (i--) |
|00000970| 7b 20 6d 63 2e 43 6f 70 | 79 28 6d 72 29 3b 20 6d |{ mc.Cop|y(mr); m|
|00000980| 72 2e 4e 65 78 74 28 29 | 3b 20 6d 63 2e 4e 65 78 |r.Next()|; mc.Nex|
|00000990| 74 28 29 3b 20 7d 0a 58 | 20 20 20 7d 0a 58 20 20 |t(); }.X| }.X |
|000009a0| 20 74 44 65 6c 65 74 65 | 28 29 3b 20 67 6d 31 2d | tDelete|(); gm1-|
|000009b0| 3e 52 65 6c 65 61 73 65 | 41 6e 64 44 65 6c 65 74 |>Release|AndDelet|
|000009c0| 65 28 29 3b 20 72 65 74 | 75 72 6e 20 67 6d 31 3b |e(); ret|urn gm1;|
|000009d0| 0a 58 7d 0a 58 0a 58 47 | 65 6e 65 72 61 6c 4d 61 |.X}.X.XG|eneralMa|
|000009e0| 74 72 69 78 2a 20 53 79 | 6d 6d 65 74 72 69 63 4d |trix* Sy|mmetricM|
|000009f0| 61 74 72 69 78 3a 3a 54 | 72 61 6e 73 70 6f 73 65 |atrix::T|ranspose|
|00000a00| 28 54 72 61 6e 73 70 6f | 73 65 64 4d 61 74 72 69 |(Transpo|sedMatri|
|00000a10| 78 2a 2c 20 4d 61 74 72 | 69 78 54 79 70 65 20 6d |x*, Matr|ixType m|
|00000a20| 74 29 0a 58 7b 20 52 45 | 50 4f 52 54 20 20 72 65 |t).X{ RE|PORT re|
|00000a30| 74 75 72 6e 20 45 76 61 | 6c 75 61 74 65 28 6d 74 |turn Eva|luate(mt|
|00000a40| 29 3b 20 7d 0a 58 0a 58 | 0a 58 47 65 6e 65 72 61 |); }.X.X|.XGenera|
|00000a50| 6c 4d 61 74 72 69 78 2a | 20 44 69 61 67 6f 6e 61 |lMatrix*| Diagona|
|00000a60| 6c 4d 61 74 72 69 78 3a | 3a 54 72 61 6e 73 70 6f |lMatrix:|:Transpo|
|00000a70| 73 65 28 54 72 61 6e 73 | 70 6f 73 65 64 4d 61 74 |se(Trans|posedMat|
|00000a80| 72 69 78 2a 2c 20 4d 61 | 74 72 69 78 54 79 70 65 |rix*, Ma|trixType|
|00000a90| 20 6d 74 29 0a 58 7b 20 | 52 45 50 4f 52 54 20 72 | mt).X{ |REPORT r|
|00000aa0| 65 74 75 72 6e 20 45 76 | 61 6c 75 61 74 65 28 6d |eturn Ev|aluate(m|
|00000ab0| 74 29 3b 20 7d 0a 58 0a | 58 42 6f 6f 6c 65 61 6e |t); }.X.|XBoolean|
|00000ac0| 20 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 3a 3a | General|Matrix::|
|00000ad0| 49 73 5a 65 72 6f 28 29 | 20 63 6f 6e 73 74 0a 58 |IsZero()| const.X|
|00000ae0| 7b 0a 58 20 20 20 52 45 | 50 4f 52 54 0a 58 20 20 |{.X RE|PORT.X |
|00000af0| 20 52 65 61 6c 2a 20 73 | 3d 73 74 6f 72 65 3b 20 | Real* s|=store; |
|00000b00| 69 6e 74 20 69 3d 73 74 | 6f 72 61 67 65 3b 0a 58 |int i=st|orage;.X|
|00000b10| 20 20 20 77 68 69 6c 65 | 20 28 69 2d 2d 29 20 7b | while| (i--) {|
|00000b20| 20 69 66 20 28 2a 73 2b | 2b 29 20 72 65 74 75 72 | if (*s+|+) retur|
|00000b30| 6e 20 46 41 4c 53 45 3b | 20 7d 0a 58 20 20 20 72 |n FALSE;| }.X r|
|00000b40| 65 74 75 72 6e 20 54 52 | 55 45 3b 0a 58 7d 0a 58 |eturn TR|UE;.X}.X|
|00000b50| 0a 58 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |.XGenera|lMatrix*|
|00000b60| 20 43 6f 6c 75 6d 6e 56 | 65 63 74 6f 72 3a 3a 54 | ColumnV|ector::T|
|00000b70| 72 61 6e 73 70 6f 73 65 | 28 54 72 61 6e 73 70 6f |ranspose|(Transpo|
|00000b80| 73 65 64 4d 61 74 72 69 | 78 2a 2c 20 4d 61 74 72 |sedMatri|x*, Matr|
|00000b90| 69 78 54 79 70 65 20 6d | 74 29 0a 58 7b 0a 58 20 |ixType m|t).X{.X |
|00000ba0| 20 20 52 45 50 4f 52 54 | 0a 58 20 20 20 47 65 6e | REPORT|.X Gen|
|00000bb0| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 67 6d 78 20 |eralMatr|ix* gmx |
|00000bc0| 3d 20 6e 65 77 20 52 6f | 77 56 65 63 74 6f 72 3b |= new Ro|wVector;|
|00000bd0| 20 4d 61 74 72 69 78 45 | 72 72 6f 72 4e 6f 53 70 | MatrixE|rrorNoSp|
|00000be0| 61 63 65 28 67 6d 78 29 | 3b 0a 58 20 20 20 67 6d |ace(gmx)|;.X gm|
|00000bf0| 78 2d 3e 6e 72 6f 77 73 | 20 3d 20 31 3b 20 67 6d |x->nrows| = 1; gm|
|00000c00| 78 2d 3e 6e 63 6f 6c 73 | 20 3d 20 67 6d 78 2d 3e |x->ncols| = gmx->|
|00000c10| 73 74 6f 72 61 67 65 20 | 3d 20 73 74 6f 72 61 67 |storage |= storag|
|00000c20| 65 3b 0a 58 20 20 20 72 | 65 74 75 72 6e 20 42 6f |e;.X r|eturn Bo|
|00000c30| 72 72 6f 77 53 74 6f 72 | 65 28 67 6d 78 2c 6d 74 |rrowStor|e(gmx,mt|
|00000c40| 29 3b 0a 58 7d 0a 58 0a | 58 47 65 6e 65 72 61 6c |);.X}.X.|XGeneral|
|00000c50| 4d 61 74 72 69 78 2a 20 | 52 6f 77 56 65 63 74 6f |Matrix* |RowVecto|
|00000c60| 72 3a 3a 54 72 61 6e 73 | 70 6f 73 65 28 54 72 61 |r::Trans|pose(Tra|
|00000c70| 6e 73 70 6f 73 65 64 4d | 61 74 72 69 78 2a 2c 20 |nsposedM|atrix*, |
|00000c80| 4d 61 74 72 69 78 54 79 | 70 65 20 6d 74 29 0a 58 |MatrixTy|pe mt).X|
|00000c90| 7b 0a 58 20 20 20 52 45 | 50 4f 52 54 0a 58 20 20 |{.X RE|PORT.X |
|00000ca0| 20 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 | General|Matrix* |
|00000cb0| 67 6d 78 20 3d 20 6e 65 | 77 20 43 6f 6c 75 6d 6e |gmx = ne|w Column|
|00000cc0| 56 65 63 74 6f 72 3b 20 | 4d 61 74 72 69 78 45 72 |Vector; |MatrixEr|
|00000cd0| 72 6f 72 4e 6f 53 70 61 | 63 65 28 67 6d 78 29 3b |rorNoSpa|ce(gmx);|
|00000ce0| 0a 58 20 20 20 67 6d 78 | 2d 3e 6e 63 6f 6c 73 20 |.X gmx|->ncols |
|00000cf0| 3d 20 31 3b 20 67 6d 78 | 2d 3e 6e 72 6f 77 73 20 |= 1; gmx|->nrows |
|00000d00| 3d 20 67 6d 78 2d 3e 73 | 74 6f 72 61 67 65 20 3d |= gmx->s|torage =|
|00000d10| 20 73 74 6f 72 61 67 65 | 3b 0a 58 20 20 20 72 65 | storage|;.X re|
|00000d20| 74 75 72 6e 20 42 6f 72 | 72 6f 77 53 74 6f 72 65 |turn Bor|rowStore|
|00000d30| 28 67 6d 78 2c 6d 74 29 | 3b 0a 58 7d 0a 58 0a 58 |(gmx,mt)|;.X}.X.X|
|00000d40| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 47 |GeneralM|atrix* G|
|00000d50| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 3a 3a 45 76 |eneralMa|trix::Ev|
|00000d60| 61 6c 75 61 74 65 28 4d | 61 74 72 69 78 54 79 70 |aluate(M|atrixTyp|
|00000d70| 65 20 6d 74 29 0a 58 7b | 0a 58 20 20 20 69 66 20 |e mt).X{|.X if |
|00000d80| 28 43 6f 6d 70 61 72 65 | 28 74 68 69 73 2d 3e 54 |(Compare|(this->T|
|00000d90| 79 70 65 28 29 2c 6d 74 | 29 29 20 7b 20 52 45 50 |ype(),mt|)) { REP|
|00000da0| 4f 52 54 20 72 65 74 75 | 72 6e 20 74 68 69 73 3b |ORT retu|rn this;|
|00000db0| 20 7d 0a 58 20 20 20 52 | 45 50 4f 52 54 0a 58 20 | }.X R|EPORT.X |
|00000dc0| 20 20 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a | Genera|lMatrix*|
|00000dd0| 20 67 6d 78 20 3d 20 6d | 74 2e 4e 65 77 28 6e 72 | gmx = m|t.New(nr|
|00000de0| 6f 77 73 2c 6e 63 6f 6c | 73 2c 74 68 69 73 29 3b |ows,ncol|s,this);|
|00000df0| 0a 58 20 20 20 4d 61 74 | 72 69 78 52 6f 77 20 6d |.X Mat|rixRow m|
|00000e00| 72 28 74 68 69 73 2c 20 | 4c 6f 61 64 4f 6e 45 6e |r(this, |LoadOnEn|
|00000e10| 74 72 79 29 3b 0a 58 20 | 20 20 4d 61 74 72 69 78 |try);.X | Matrix|
|00000e20| 52 6f 77 20 6d 72 78 28 | 67 6d 78 2c 20 53 74 6f |Row mrx(|gmx, Sto|
|00000e30| 72 65 4f 6e 45 78 69 74 | 2b 44 69 72 65 63 74 50 |reOnExit|+DirectP|
|00000e40| 61 72 74 29 3b 0a 58 20 | 20 20 69 6e 74 20 69 3d |art);.X | int i=|
|00000e50| 6e 72 6f 77 73 3b 0a 58 | 20 20 20 77 68 69 6c 65 |nrows;.X| while|
|00000e60| 20 28 69 2d 2d 29 20 7b | 20 6d 72 78 2e 43 6f 70 | (i--) {| mrx.Cop|
|00000e70| 79 28 6d 72 29 3b 20 6d | 72 78 2e 4e 65 78 74 28 |y(mr); m|rx.Next(|
|00000e80| 29 3b 20 6d 72 2e 4e 65 | 78 74 28 29 3b 20 7d 0a |); mr.Ne|xt(); }.|
|00000e90| 58 20 20 20 74 44 65 6c | 65 74 65 28 29 3b 20 67 |X tDel|ete(); g|
|00000ea0| 6d 78 2d 3e 52 65 6c 65 | 61 73 65 41 6e 64 44 65 |mx->Rele|aseAndDe|
|00000eb0| 6c 65 74 65 28 29 3b 20 | 72 65 74 75 72 6e 20 67 |lete(); |return g|
|00000ec0| 6d 78 3b 0a 58 7d 0a 58 | 0a 58 47 65 6e 65 72 61 |mx;.X}.X|.XGenera|
|00000ed0| 6c 4d 61 74 72 69 78 2a | 20 43 6f 6e 73 74 4d 61 |lMatrix*| ConstMa|
|00000ee0| 74 72 69 78 3a 3a 45 76 | 61 6c 75 61 74 65 28 4d |trix::Ev|aluate(M|
|00000ef0| 61 74 72 69 78 54 79 70 | 65 20 6d 74 29 0a 58 7b |atrixTyp|e mt).X{|
|00000f00| 0a 58 20 20 20 69 66 20 | 28 43 6f 6d 70 61 72 65 |.X if |(Compare|
|00000f10| 28 63 67 6d 2d 3e 54 79 | 70 65 28 29 2c 6d 74 29 |(cgm->Ty|pe(),mt)|
|00000f20| 29 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 52 |).X {.|X R|
|00000f30| 45 50 4f 52 54 0a 58 23 | 69 66 64 65 66 20 54 45 |EPORT.X#|ifdef TE|
|00000f40| 4d 50 53 5f 44 45 53 54 | 52 4f 59 45 44 5f 51 55 |MPS_DEST|ROYED_QU|
|00000f50| 49 43 4b 4c 59 0a 58 20 | 20 20 20 20 20 47 65 6e |ICKLY.X | Gen|
|00000f60| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 67 6d 78 20 |eralMatr|ix* gmx |
|00000f70| 3d 20 28 47 65 6e 65 72 | 61 6c 4d 61 74 72 69 78 |= (Gener|alMatrix|
|00000f80| 2a 29 63 67 6d 3b 20 64 | 65 6c 65 74 65 20 74 68 |*)cgm; d|elete th|
|00000f90| 69 73 3b 20 72 65 74 75 | 72 6e 20 67 6d 78 3b 0a |is; retu|rn gmx;.|
|00000fa0| 58 23 65 6c 73 65 0a 58 | 20 20 20 20 20 20 72 65 |X#else.X| re|
|00000fb0| 74 75 72 6e 20 28 47 65 | 6e 65 72 61 6c 4d 61 74 |turn (Ge|neralMat|
|00000fc0| 72 69 78 2a 29 63 67 6d | 3b 0a 58 23 65 6e 64 69 |rix*)cgm|;.X#endi|
|00000fd0| 66 0a 58 20 20 20 7d 0a | 58 20 20 20 52 45 50 4f |f.X }.|X REPO|
|00000fe0| 52 54 0a 58 20 20 20 47 | 65 6e 65 72 61 6c 4d 61 |RT.X G|eneralMa|
|00000ff0| 74 72 69 78 2a 20 67 6d | 78 20 3d 20 63 67 6d 2d |trix* gm|x = cgm-|
|00001000| 3e 54 79 70 65 28 29 2e | 4e 65 77 28 63 67 6d 2d |>Type().|New(cgm-|
|00001010| 3e 4e 72 6f 77 73 28 29 | 2c 63 67 6d 2d 3e 4e 63 |>Nrows()|,cgm->Nc|
|00001020| 6f 6c 73 28 29 2c 74 68 | 69 73 29 3b 0a 58 20 20 |ols(),th|is);.X |
|00001030| 20 4d 61 74 72 69 78 52 | 6f 77 20 6d 72 28 28 47 | MatrixR|ow mr((G|
|00001040| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 29 63 67 |eneralMa|trix*)cg|
|00001050| 6d 2c 20 4c 6f 61 64 4f | 6e 45 6e 74 72 79 29 3b |m, LoadO|nEntry);|
|00001060| 2f 2f 61 73 73 75 6d 65 | 20 77 6f 6e 27 74 20 63 |//assume| won't c|
|00001070| 68 61 6e 67 65 20 74 68 | 69 73 0a 58 20 20 20 4d |hange th|is.X M|
|00001080| 61 74 72 69 78 52 6f 77 | 20 6d 72 78 28 67 6d 78 |atrixRow| mrx(gmx|
|00001090| 2c 20 53 74 6f 72 65 4f | 6e 45 78 69 74 2b 44 69 |, StoreO|nExit+Di|
|000010a0| 72 65 63 74 50 61 72 74 | 29 3b 0a 58 20 20 20 69 |rectPart|);.X i|
|000010b0| 6e 74 20 69 3d 63 67 6d | 2d 3e 4e 72 6f 77 73 28 |nt i=cgm|->Nrows(|
|000010c0| 29 3b 0a 58 20 20 20 77 | 68 69 6c 65 20 28 69 2d |);.X w|hile (i-|
|000010d0| 2d 29 20 7b 20 6d 72 78 | 2e 43 6f 70 79 28 6d 72 |-) { mrx|.Copy(mr|
|000010e0| 29 3b 20 6d 72 78 2e 4e | 65 78 74 28 29 3b 20 6d |); mrx.N|ext(); m|
|000010f0| 72 2e 4e 65 78 74 28 29 | 3b 20 7d 0a 58 20 20 20 |r.Next()|; }.X |
|00001100| 67 6d 78 2d 3e 52 65 6c | 65 61 73 65 41 6e 64 44 |gmx->Rel|easeAndD|
|00001110| 65 6c 65 74 65 28 29 3b | 0a 58 23 69 66 64 65 66 |elete();|.X#ifdef|
|00001120| 20 54 45 4d 50 53 5f 44 | 45 53 54 52 4f 59 45 44 | TEMPS_D|ESTROYED|
|00001130| 5f 51 55 49 43 4b 4c 59 | 0a 58 20 20 20 64 65 6c |_QUICKLY|.X del|
|00001140| 65 74 65 20 74 68 69 73 | 3b 0a 58 23 65 6e 64 69 |ete this|;.X#endi|
|00001150| 66 0a 58 20 20 20 72 65 | 74 75 72 6e 20 67 6d 78 |f.X re|turn gmx|
|00001160| 3b 20 2f 2f 20 6e 6f 20 | 74 44 65 6c 65 74 65 0a |; // no |tDelete.|
|00001170| 58 7d 0a 58 0a 58 47 65 | 6e 65 72 61 6c 4d 61 74 |X}.X.XGe|neralMat|
|00001180| 72 69 78 2a 20 53 68 69 | 66 74 65 64 4d 61 74 72 |rix* Shi|ftedMatr|
|00001190| 69 78 3a 3a 45 76 61 6c | 75 61 74 65 28 4d 61 74 |ix::Eval|uate(Mat|
|000011a0| 72 69 78 54 79 70 65 20 | 6d 74 29 0a 58 7b 0a 58 |rixType |mt).X{.X|
|000011b0| 20 20 20 67 6d 3d 28 28 | 42 61 73 65 4d 61 74 72 | gm=((|BaseMatr|
|000011c0| 69 78 2a 26 29 62 6d 29 | 2d 3e 45 76 61 6c 75 61 |ix*&)bm)|->Evalua|
|000011d0| 74 65 28 29 3b 0a 58 20 | 20 20 69 6e 74 20 6e 72 |te();.X | int nr|
|000011e0| 3d 67 6d 2d 3e 4e 72 6f | 77 73 28 29 3b 20 69 6e |=gm->Nro|ws(); in|
|000011f0| 74 20 6e 63 3d 67 6d 2d | 3e 4e 63 6f 6c 73 28 29 |t nc=gm-|>Ncols()|
|00001200| 3b 0a 58 20 20 20 43 6f | 6d 70 61 72 65 28 67 6d |;.X Co|mpare(gm|
|00001210| 2d 3e 54 79 70 65 28 29 | 2e 41 64 64 45 71 75 61 |->Type()|.AddEqua|
|00001220| 6c 45 6c 28 29 2c 6d 74 | 29 3b 0a 58 20 20 20 69 |lEl(),mt|);.X i|
|00001230| 66 20 28 21 28 6d 74 3d | 3d 67 6d 2d 3e 54 79 70 |f (!(mt=|=gm->Typ|
|00001240| 65 28 29 29 29 0a 58 20 | 20 20 7b 0a 58 20 20 20 |e())).X | {.X |
|00001250| 20 20 20 52 45 50 4f 52 | 54 0a 58 20 20 20 20 20 | REPOR|T.X |
|00001260| 20 47 65 6e 65 72 61 6c | 4d 61 74 72 69 78 2a 20 | General|Matrix* |
|00001270| 67 6d 78 20 3d 20 6d 74 | 2e 4e 65 77 28 6e 72 2c |gmx = mt|.New(nr,|
|00001280| 6e 63 2c 74 68 69 73 29 | 3b 0a 58 20 20 20 20 20 |nc,this)|;.X |
|00001290| 20 4d 61 74 72 69 78 52 | 6f 77 20 6d 72 28 67 6d | MatrixR|ow mr(gm|
|000012a0| 2c 20 4c 6f 61 64 4f 6e | 45 6e 74 72 79 29 3b 20 |, LoadOn|Entry); |
|000012b0| 0a 58 20 20 20 20 20 20 | 4d 61 74 72 69 78 52 6f |.X |MatrixRo|
|000012c0| 77 20 6d 72 78 28 67 6d | 78 2c 20 53 74 6f 72 65 |w mrx(gm|x, Store|
|000012d0| 4f 6e 45 78 69 74 2b 44 | 69 72 65 63 74 50 61 72 |OnExit+D|irectPar|
|000012e0| 74 29 3b 0a 58 20 20 20 | 20 20 20 77 68 69 6c 65 |t);.X | while|
|000012f0| 20 28 6e 72 2d 2d 29 20 | 7b 20 6d 72 78 2e 41 64 | (nr--) |{ mrx.Ad|
|00001300| 64 28 6d 72 2c 66 29 3b | 20 6d 72 78 2e 4e 65 78 |d(mr,f);| mrx.Nex|
|00001310| 74 28 29 3b 20 6d 72 2e | 4e 65 78 74 28 29 3b 20 |t(); mr.|Next(); |
|00001320| 7d 0a 58 20 20 20 20 20 | 20 67 6d 78 2d 3e 52 65 |}.X | gmx->Re|
|00001330| 6c 65 61 73 65 41 6e 64 | 44 65 6c 65 74 65 28 29 |leaseAnd|Delete()|
|00001340| 3b 20 67 6d 2d 3e 74 44 | 65 6c 65 74 65 28 29 3b |; gm->tD|elete();|
|00001350| 0a 58 23 69 66 64 65 66 | 20 54 45 4d 50 53 5f 44 |.X#ifdef| TEMPS_D|
|00001360| 45 53 54 52 4f 59 45 44 | 5f 51 55 49 43 4b 4c 59 |ESTROYED|_QUICKLY|
|00001370| 0a 58 20 20 20 20 20 20 | 64 65 6c 65 74 65 20 74 |.X |delete t|
|00001380| 68 69 73 3b 0a 58 23 65 | 6e 64 69 66 0a 58 20 20 |his;.X#e|ndif.X |
|00001390| 20 20 20 20 72 65 74 75 | 72 6e 20 67 6d 78 3b 0a | retu|rn gmx;.|
|000013a0| 58 20 20 20 7d 0a 58 20 | 20 20 65 6c 73 65 20 69 |X }.X | else i|
|000013b0| 66 20 28 67 6d 2d 3e 72 | 65 75 73 65 28 29 29 0a |f (gm->r|euse()).|
|000013c0| 58 20 20 20 7b 0a 58 20 | 20 20 20 20 20 52 45 50 |X {.X | REP|
|000013d0| 4f 52 54 20 67 6d 2d 3e | 41 64 64 28 66 29 3b 0a |ORT gm->|Add(f);.|
|000013e0| 58 23 69 66 64 65 66 20 | 54 45 4d 50 53 5f 44 45 |X#ifdef |TEMPS_DE|
|000013f0| 53 54 52 4f 59 45 44 5f | 51 55 49 43 4b 4c 59 0a |STROYED_|QUICKLY.|
|00001400| 58 20 20 20 20 20 20 47 | 65 6e 65 72 61 6c 4d 61 |X G|eneralMa|
|00001410| 74 72 69 78 2a 20 67 6d | 78 20 3d 20 67 6d 3b 20 |trix* gm|x = gm; |
|00001420| 64 65 6c 65 74 65 20 74 | 68 69 73 3b 20 72 65 74 |delete t|his; ret|
|00001430| 75 72 6e 20 67 6d 78 3b | 0a 58 23 65 6c 73 65 0a |urn gmx;|.X#else.|
|00001440| 58 20 20 20 20 20 20 72 | 65 74 75 72 6e 20 67 6d |X r|eturn gm|
|00001450| 3b 0a 58 23 65 6e 64 69 | 66 0a 58 20 20 20 7d 0a |;.X#endi|f.X }.|
|00001460| 58 20 20 20 65 6c 73 65 | 0a 58 20 20 20 7b 0a 58 |X else|.X {.X|
|00001470| 20 20 20 20 20 20 52 45 | 50 4f 52 54 20 47 65 6e | RE|PORT Gen|
|00001480| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 67 6d 79 20 |eralMatr|ix* gmy |
|00001490| 3d 20 67 6d 2d 3e 54 79 | 70 65 28 29 2e 4e 65 77 |= gm->Ty|pe().New|
|000014a0| 28 6e 72 2c 6e 63 2c 74 | 68 69 73 29 3b 0a 58 20 |(nr,nc,t|his);.X |
|000014b0| 20 20 20 20 20 67 6d 79 | 2d 3e 52 65 6c 65 61 73 | gmy|->Releas|
|000014c0| 65 41 6e 64 44 65 6c 65 | 74 65 28 29 3b 20 67 6d |eAndDele|te(); gm|
|000014d0| 79 2d 3e 41 64 64 28 67 | 6d 2c 66 29 3b 0a 58 23 |y->Add(g|m,f);.X#|
|000014e0| 69 66 64 65 66 20 54 45 | 4d 50 53 5f 44 45 53 54 |ifdef TE|MPS_DEST|
|000014f0| 52 4f 59 45 44 5f 51 55 | 49 43 4b 4c 59 0a 58 20 |ROYED_QU|ICKLY.X |
|00001500| 20 20 20 20 20 64 65 6c | 65 74 65 20 74 68 69 73 | del|ete this|
|00001510| 3b 0a 58 23 65 6e 64 69 | 66 0a 58 20 20 20 20 20 |;.X#endi|f.X |
|00001520| 20 72 65 74 75 72 6e 20 | 67 6d 79 3b 0a 58 20 20 | return |gmy;.X |
|00001530| 20 7d 0a 58 7d 0a 58 0a | 58 47 65 6e 65 72 61 6c | }.X}.X.|XGeneral|
|00001540| 4d 61 74 72 69 78 2a 20 | 53 63 61 6c 65 64 4d 61 |Matrix* |ScaledMa|
|00001550| 74 72 69 78 3a 3a 45 76 | 61 6c 75 61 74 65 28 4d |trix::Ev|aluate(M|
|00001560| 61 74 72 69 78 54 79 70 | 65 20 6d 74 29 0a 58 7b |atrixTyp|e mt).X{|
|00001570| 0a 58 20 20 20 67 6d 3d | 28 28 42 61 73 65 4d 61 |.X gm=|((BaseMa|
|00001580| 74 72 69 78 2a 26 29 62 | 6d 29 2d 3e 45 76 61 6c |trix*&)b|m)->Eval|
|00001590| 75 61 74 65 28 29 3b 0a | 58 20 20 20 69 6e 74 20 |uate();.|X int |
|000015a0| 6e 72 3d 67 6d 2d 3e 4e | 72 6f 77 73 28 29 3b 20 |nr=gm->N|rows(); |
|000015b0| 69 6e 74 20 6e 63 3d 67 | 6d 2d 3e 4e 63 6f 6c 73 |int nc=g|m->Ncols|
|000015c0| 28 29 3b 0a 58 20 20 20 | 69 66 20 28 43 6f 6d 70 |();.X |if (Comp|
|000015d0| 61 72 65 28 67 6d 2d 3e | 54 79 70 65 28 29 2c 6d |are(gm->|Type(),m|
|000015e0| 74 29 29 0a 58 20 20 20 | 7b 0a 58 20 20 20 20 20 |t)).X |{.X |
|000015f0| 20 69 66 20 28 67 6d 2d | 3e 72 65 75 73 65 28 29 | if (gm-|>reuse()|
|00001600| 29 0a 58 20 20 20 20 20 | 20 7b 0a 58 20 20 20 20 |).X | {.X |
|00001610| 20 20 20 20 20 52 45 50 | 4f 52 54 20 67 6d 2d 3e | REP|ORT gm->|
|00001620| 4d 75 6c 74 69 70 6c 79 | 28 66 29 3b 0a 58 23 69 |Multiply|(f);.X#i|
|00001630| 66 64 65 66 20 54 45 4d | 50 53 5f 44 45 53 54 52 |fdef TEM|PS_DESTR|
|00001640| 4f 59 45 44 5f 51 55 49 | 43 4b 4c 59 0a 58 20 20 |OYED_QUI|CKLY.X |
|00001650| 20 20 20 20 20 20 20 47 | 65 6e 65 72 61 6c 4d 61 | G|eneralMa|
|00001660| 74 72 69 78 2a 20 67 6d | 78 20 3d 20 67 6d 3b 20 |trix* gm|x = gm; |
|00001670| 64 65 6c 65 74 65 20 74 | 68 69 73 3b 20 72 65 74 |delete t|his; ret|
|00001680| 75 72 6e 20 67 6d 78 3b | 0a 58 23 65 6c 73 65 0a |urn gmx;|.X#else.|
|00001690| 58 20 20 20 20 20 20 20 | 20 20 72 65 74 75 72 6e |X | return|
|000016a0| 20 67 6d 3b 0a 58 23 65 | 6e 64 69 66 0a 58 20 20 | gm;.X#e|ndif.X |
|000016b0| 20 20 20 20 7d 0a 58 20 | 20 20 20 20 20 65 6c 73 | }.X | els|
|000016c0| 65 0a 58 20 20 20 20 20 | 20 7b 0a 58 20 20 20 20 |e.X | {.X |
|000016d0| 20 20 20 20 20 52 45 50 | 4f 52 54 20 47 65 6e 65 | REP|ORT Gene|
|000016e0| 72 61 6c 4d 61 74 72 69 | 78 2a 20 67 6d 78 20 3d |ralMatri|x* gmx =|
|000016f0| 20 67 6d 2d 3e 54 79 70 | 65 28 29 2e 4e 65 77 28 | gm->Typ|e().New(|
|00001700| 6e 72 2c 6e 63 2c 74 68 | 69 73 29 3b 0a 58 20 20 |nr,nc,th|is);.X |
|00001710| 20 20 20 20 20 20 20 67 | 6d 78 2d 3e 52 65 6c 65 | g|mx->Rele|
|00001720| 61 73 65 41 6e 64 44 65 | 6c 65 74 65 28 29 3b 20 |aseAndDe|lete(); |
|00001730| 67 6d 78 2d 3e 4d 75 6c | 74 69 70 6c 79 28 67 6d |gmx->Mul|tiply(gm|
|00001740| 2c 66 29 3b 0a 58 23 69 | 66 64 65 66 20 54 45 4d |,f);.X#i|fdef TEM|
|00001750| 50 53 5f 44 45 53 54 52 | 4f 59 45 44 5f 51 55 49 |PS_DESTR|OYED_QUI|
|00001760| 43 4b 4c 59 0a 58 20 20 | 20 20 20 20 20 20 20 64 |CKLY.X | d|
|00001770| 65 6c 65 74 65 20 74 68 | 69 73 3b 0a 58 23 65 6e |elete th|is;.X#en|
|00001780| 64 69 66 0a 58 20 20 20 | 20 20 20 20 20 20 72 65 |dif.X | re|
|00001790| 74 75 72 6e 20 67 6d 78 | 3b 0a 58 20 20 20 20 20 |turn gmx|;.X |
|000017a0| 20 7d 0a 58 20 20 20 7d | 0a 58 20 20 20 65 6c 73 | }.X }|.X els|
|000017b0| 65 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 52 |e.X {.|X R|
|000017c0| 45 50 4f 52 54 0a 58 20 | 20 20 20 20 20 47 65 6e |EPORT.X | Gen|
|000017d0| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 67 6d 78 20 |eralMatr|ix* gmx |
|000017e0| 3d 20 6d 74 2e 4e 65 77 | 28 6e 72 2c 6e 63 2c 74 |= mt.New|(nr,nc,t|
|000017f0| 68 69 73 29 3b 0a 58 20 | 20 20 20 20 20 4d 61 74 |his);.X | Mat|
|00001800| 72 69 78 52 6f 77 20 6d | 72 28 67 6d 2c 20 4c 6f |rixRow m|r(gm, Lo|
|00001810| 61 64 4f 6e 45 6e 74 72 | 79 29 3b 20 0a 58 20 20 |adOnEntr|y); .X |
|00001820| 20 20 20 20 4d 61 74 72 | 69 78 52 6f 77 20 6d 72 | Matr|ixRow mr|
|00001830| 78 28 67 6d 78 2c 20 53 | 74 6f 72 65 4f 6e 45 78 |x(gmx, S|toreOnEx|
|00001840| 69 74 2b 44 69 72 65 63 | 74 50 61 72 74 29 3b 0a |it+Direc|tPart);.|
|00001850| 58 20 20 20 20 20 20 77 | 68 69 6c 65 20 28 6e 72 |X w|hile (nr|
|00001860| 2d 2d 29 20 7b 20 6d 72 | 78 2e 4d 75 6c 74 69 70 |--) { mr|x.Multip|
|00001870| 6c 79 28 6d 72 2c 66 29 | 3b 20 6d 72 78 2e 4e 65 |ly(mr,f)|; mrx.Ne|
|00001880| 78 74 28 29 3b 20 6d 72 | 2e 4e 65 78 74 28 29 3b |xt(); mr|.Next();|
|00001890| 20 7d 0a 58 20 20 20 20 | 20 20 67 6d 78 2d 3e 52 | }.X | gmx->R|
|000018a0| 65 6c 65 61 73 65 41 6e | 64 44 65 6c 65 74 65 28 |eleaseAn|dDelete(|
|000018b0| 29 3b 20 67 6d 2d 3e 74 | 44 65 6c 65 74 65 28 29 |); gm->t|Delete()|
|000018c0| 3b 0a 58 23 69 66 64 65 | 66 20 54 45 4d 50 53 5f |;.X#ifde|f TEMPS_|
|000018d0| 44 45 53 54 52 4f 59 45 | 44 5f 51 55 49 43 4b 4c |DESTROYE|D_QUICKL|
|000018e0| 59 0a 58 20 20 20 20 20 | 20 64 65 6c 65 74 65 20 |Y.X | delete |
|000018f0| 74 68 69 73 3b 0a 58 23 | 65 6e 64 69 66 0a 58 20 |this;.X#|endif.X |
|00001900| 20 20 20 20 20 72 65 74 | 75 72 6e 20 67 6d 78 3b | ret|urn gmx;|
|00001910| 0a 58 20 20 20 7d 0a 58 | 7d 0a 58 0a 58 47 65 6e |.X }.X|}.X.XGen|
|00001920| 65 72 61 6c 4d 61 74 72 | 69 78 2a 20 4e 65 67 61 |eralMatr|ix* Nega|
|00001930| 74 65 64 4d 61 74 72 69 | 78 3a 3a 45 76 61 6c 75 |tedMatri|x::Evalu|
|00001940| 61 74 65 28 4d 61 74 72 | 69 78 54 79 70 65 20 6d |ate(Matr|ixType m|
|00001950| 74 29 0a 58 7b 0a 58 20 | 20 20 67 6d 3d 28 28 42 |t).X{.X | gm=((B|
|00001960| 61 73 65 4d 61 74 72 69 | 78 2a 26 29 62 6d 29 2d |aseMatri|x*&)bm)-|
|00001970| 3e 45 76 61 6c 75 61 74 | 65 28 29 3b 0a 58 20 20 |>Evaluat|e();.X |
|00001980| 20 69 6e 74 20 6e 72 3d | 67 6d 2d 3e 4e 72 6f 77 | int nr=|gm->Nrow|
|00001990| 73 28 29 3b 20 69 6e 74 | 20 6e 63 3d 67 6d 2d 3e |s(); int| nc=gm->|
|000019a0| 4e 63 6f 6c 73 28 29 3b | 0a 58 20 20 20 69 66 20 |Ncols();|.X if |
|000019b0| 28 43 6f 6d 70 61 72 65 | 28 67 6d 2d 3e 54 79 70 |(Compare|(gm->Typ|
|000019c0| 65 28 29 2c 6d 74 29 29 | 0a 58 20 20 20 7b 0a 58 |e(),mt))|.X {.X|
|000019d0| 20 20 20 20 20 20 69 66 | 20 28 67 6d 2d 3e 72 65 | if| (gm->re|
|000019e0| 75 73 65 28 29 29 0a 58 | 20 20 20 20 20 20 7b 0a |use()).X| {.|
|000019f0| 58 20 20 20 20 20 20 20 | 20 20 52 45 50 4f 52 54 |X | REPORT|
|00001a00| 20 67 6d 2d 3e 4e 65 67 | 61 74 65 28 29 3b 0a 58 | gm->Neg|ate();.X|
|00001a10| 23 69 66 64 65 66 20 54 | 45 4d 50 53 5f 44 45 53 |#ifdef T|EMPS_DES|
|00001a20| 54 52 4f 59 45 44 5f 51 | 55 49 43 4b 4c 59 0a 58 |TROYED_Q|UICKLY.X|
|00001a30| 20 20 20 20 20 20 20 20 | 20 47 65 6e 65 72 61 6c | | General|
|00001a40| 4d 61 74 72 69 78 2a 20 | 67 6d 78 20 3d 20 67 6d |Matrix* |gmx = gm|
|00001a50| 3b 20 64 65 6c 65 74 65 | 20 74 68 69 73 3b 20 72 |; delete| this; r|
|00001a60| 65 74 75 72 6e 20 67 6d | 78 3b 0a 58 23 65 6c 73 |eturn gm|x;.X#els|
|00001a70| 65 0a 58 20 20 20 20 20 | 20 20 20 20 72 65 74 75 |e.X | retu|
|00001a80| 72 6e 20 67 6d 3b 0a 58 | 23 65 6e 64 69 66 0a 58 |rn gm;.X|#endif.X|
|00001a90| 20 20 20 20 20 20 7d 0a | 58 20 20 20 20 20 20 65 | }.|X e|
|00001aa0| 6c 73 65 0a 58 20 20 20 | 20 20 20 7b 0a 58 20 20 |lse.X | {.X |
|00001ab0| 20 20 20 20 20 20 20 52 | 45 50 4f 52 54 0a 58 20 | R|EPORT.X |
|00001ac0| 20 20 20 20 20 20 20 20 | 47 65 6e 65 72 61 6c 4d | |GeneralM|
|00001ad0| 61 74 72 69 78 2a 20 67 | 6d 78 20 3d 20 67 6d 2d |atrix* g|mx = gm-|
|00001ae0| 3e 54 79 70 65 28 29 2e | 4e 65 77 28 6e 72 2c 6e |>Type().|New(nr,n|
|00001af0| 63 2c 74 68 69 73 29 3b | 0a 58 20 20 20 20 20 20 |c,this);|.X |
|00001b00| 20 20 20 67 6d 78 2d 3e | 52 65 6c 65 61 73 65 41 | gmx->|ReleaseA|
|00001b10| 6e 64 44 65 6c 65 74 65 | 28 29 3b 20 67 6d 78 2d |ndDelete|(); gmx-|
|00001b20| 3e 4e 65 67 61 74 65 28 | 67 6d 29 3b 0a 58 23 69 |>Negate(|gm);.X#i|
|00001b30| 66 64 65 66 20 54 45 4d | 50 53 5f 44 45 53 54 52 |fdef TEM|PS_DESTR|
|00001b40| 4f 59 45 44 5f 51 55 49 | 43 4b 4c 59 0a 58 20 20 |OYED_QUI|CKLY.X |
|00001b50| 20 20 20 20 20 20 20 64 | 65 6c 65 74 65 20 74 68 | d|elete th|
|00001b60| 69 73 3b 0a 58 23 65 6e | 64 69 66 0a 58 20 20 20 |is;.X#en|dif.X |
|00001b70| 20 20 20 20 20 20 72 65 | 74 75 72 6e 20 67 6d 78 | re|turn gmx|
|00001b80| 3b 0a 58 20 20 20 20 20 | 20 7d 0a 58 20 20 20 7d |;.X | }.X }|
|00001b90| 0a 58 20 20 20 65 6c 73 | 65 0a 58 20 20 20 7b 0a |.X els|e.X {.|
|00001ba0| 58 20 20 20 20 20 20 52 | 45 50 4f 52 54 0a 58 20 |X R|EPORT.X |
|00001bb0| 20 20 20 20 20 47 65 6e | 65 72 61 6c 4d 61 74 72 | Gen|eralMatr|
|00001bc0| 69 78 2a 20 67 6d 78 20 | 3d 20 6d 74 2e 4e 65 77 |ix* gmx |= mt.New|
|00001bd0| 28 6e 72 2c 6e 63 2c 74 | 68 69 73 29 3b 0a 58 20 |(nr,nc,t|his);.X |
|00001be0| 20 20 20 20 20 4d 61 74 | 72 69 78 52 6f 77 20 6d | Mat|rixRow m|
|00001bf0| 72 28 67 6d 2c 20 4c 6f | 61 64 4f 6e 45 6e 74 72 |r(gm, Lo|adOnEntr|
|00001c00| 79 29 3b 20 0a 58 20 20 | 20 20 20 20 4d 61 74 72 |y); .X | Matr|
|00001c10| 69 78 52 6f 77 20 6d 72 | 78 28 67 6d 78 2c 20 53 |ixRow mr|x(gmx, S|
|00001c20| 74 6f 72 65 4f 6e 45 78 | 69 74 2b 44 69 72 65 63 |toreOnEx|it+Direc|
|00001c30| 74 50 61 72 74 29 3b 0a | 58 20 20 20 20 20 20 77 |tPart);.|X w|
|00001c40| 68 69 6c 65 20 28 6e 72 | 2d 2d 29 20 7b 20 6d 72 |hile (nr|--) { mr|
|00001c50| 78 2e 4e 65 67 61 74 65 | 28 6d 72 29 3b 20 6d 72 |x.Negate|(mr); mr|
|00001c60| 78 2e 4e 65 78 74 28 29 | 3b 20 6d 72 2e 4e 65 78 |x.Next()|; mr.Nex|
|00001c70| 74 28 29 3b 20 7d 0a 58 | 20 20 20 20 20 20 67 6d |t(); }.X| gm|
|00001c80| 78 2d 3e 52 65 6c 65 61 | 73 65 41 6e 64 44 65 6c |x->Relea|seAndDel|
|00001c90| 65 74 65 28 29 3b 20 67 | 6d 2d 3e 74 44 65 6c 65 |ete(); g|m->tDele|
|00001ca0| 74 65 28 29 3b 0a 58 23 | 69 66 64 65 66 20 54 45 |te();.X#|ifdef TE|
|00001cb0| 4d 50 53 5f 44 45 53 54 | 52 4f 59 45 44 5f 51 55 |MPS_DEST|ROYED_QU|
|00001cc0| 49 43 4b 4c 59 0a 58 20 | 20 20 20 20 20 64 65 6c |ICKLY.X | del|
|00001cd0| 65 74 65 20 74 68 69 73 | 3b 0a 58 23 65 6e 64 69 |ete this|;.X#endi|
|00001ce0| 66 0a 58 20 20 20 20 20 | 20 72 65 74 75 72 6e 20 |f.X | return |
|00001cf0| 67 6d 78 3b 0a 58 20 20 | 20 7d 0a 58 7d 20 20 20 |gmx;.X | }.X} |
|00001d00| 0a 58 0a 58 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |.X.XGene|ralMatri|
|00001d10| 78 2a 20 54 72 61 6e 73 | 70 6f 73 65 64 4d 61 74 |x* Trans|posedMat|
|00001d20| 72 69 78 3a 3a 45 76 61 | 6c 75 61 74 65 28 4d 61 |rix::Eva|luate(Ma|
|00001d30| 74 72 69 78 54 79 70 65 | 20 6d 74 29 0a 58 7b 0a |trixType| mt).X{.|
|00001d40| 58 20 20 20 52 45 50 4f | 52 54 0a 58 20 20 20 67 |X REPO|RT.X g|
|00001d50| 6d 3d 28 28 42 61 73 65 | 4d 61 74 72 69 78 2a 26 |m=((Base|Matrix*&|
|00001d60| 29 62 6d 29 2d 3e 45 76 | 61 6c 75 61 74 65 28 29 |)bm)->Ev|aluate()|
|00001d70| 3b 0a 58 20 20 20 43 6f | 6d 70 61 72 65 28 67 6d |;.X Co|mpare(gm|
|00001d80| 2d 3e 54 79 70 65 28 29 | 2e 74 28 29 2c 6d 74 29 |->Type()|.t(),mt)|
|00001d90| 3b 0a 58 20 20 20 47 65 | 6e 65 72 61 6c 4d 61 74 |;.X Ge|neralMat|
|00001da0| 72 69 78 2a 20 67 6d 78 | 3d 67 6d 2d 3e 54 72 61 |rix* gmx|=gm->Tra|
|00001db0| 6e 73 70 6f 73 65 28 74 | 68 69 73 2c 20 6d 74 29 |nspose(t|his, mt)|
|00001dc0| 3b 0a 58 23 69 66 64 65 | 66 20 54 45 4d 50 53 5f |;.X#ifde|f TEMPS_|
|00001dd0| 44 45 53 54 52 4f 59 45 | 44 5f 51 55 49 43 4b 4c |DESTROYE|D_QUICKL|
|00001de0| 59 0a 58 20 20 20 64 65 | 6c 65 74 65 20 74 68 69 |Y.X de|lete thi|
|00001df0| 73 3b 0a 58 23 65 6e 64 | 69 66 0a 58 20 20 20 72 |s;.X#end|if.X r|
|00001e00| 65 74 75 72 6e 20 67 6d | 78 3b 0a 58 7d 0a 58 20 |eturn gm|x;.X}.X |
|00001e10| 20 20 0a 58 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 | .XGene|ralMatri|
|00001e20| 78 2a 20 52 6f 77 65 64 | 4d 61 74 72 69 78 3a 3a |x* Rowed|Matrix::|
|00001e30| 45 76 61 6c 75 61 74 65 | 28 4d 61 74 72 69 78 54 |Evaluate|(MatrixT|
|00001e40| 79 70 65 20 6d 74 29 0a | 58 7b 0a 58 20 20 20 67 |ype mt).|X{.X g|
|00001e50| 6d 20 3d 20 28 28 42 61 | 73 65 4d 61 74 72 69 78 |m = ((Ba|seMatrix|
|00001e60| 2a 26 29 62 6d 29 2d 3e | 45 76 61 6c 75 61 74 65 |*&)bm)->|Evaluate|
|00001e70| 28 29 3b 0a 58 20 20 20 | 47 65 6e 65 72 61 6c 4d |();.X |GeneralM|
|00001e80| 61 74 72 69 78 2a 20 67 | 6d 78 20 3d 20 6e 65 77 |atrix* g|mx = new|
|00001e90| 20 52 6f 77 56 65 63 74 | 6f 72 3b 20 4d 61 74 72 | RowVect|or; Matr|
|00001ea0| 69 78 45 72 72 6f 72 4e | 6f 53 70 61 63 65 28 67 |ixErrorN|oSpace(g|
|00001eb0| 6d 78 29 3b 0a 58 20 20 | 20 67 6d 78 2d 3e 6e 72 |mx);.X | gmx->nr|
|00001ec0| 6f 77 73 20 3d 20 31 3b | 20 67 6d 78 2d 3e 6e 63 |ows = 1;| gmx->nc|
|00001ed0| 6f 6c 73 20 3d 20 67 6d | 78 2d 3e 73 74 6f 72 61 |ols = gm|x->stora|
|00001ee0| 67 65 20 3d 20 67 6d 2d | 3e 73 74 6f 72 61 67 65 |ge = gm-|>storage|
|00001ef0| 3b 0a 58 23 69 66 64 65 | 66 20 54 45 4d 50 53 5f |;.X#ifde|f TEMPS_|
|00001f00| 44 45 53 54 52 4f 59 45 | 44 5f 51 55 49 43 4b 4c |DESTROYE|D_QUICKL|
|00001f10| 59 0a 58 20 20 20 47 65 | 6e 65 72 61 6c 4d 61 74 |Y.X Ge|neralMat|
|00001f20| 72 69 78 2a 20 67 6d 79 | 20 3d 20 67 6d 3b 20 64 |rix* gmy| = gm; d|
|00001f30| 65 6c 65 74 65 20 74 68 | 69 73 3b 20 72 65 74 75 |elete th|is; retu|
|00001f40| 72 6e 20 67 6d 79 2d 3e | 42 6f 72 72 6f 77 53 74 |rn gmy->|BorrowSt|
|00001f50| 6f 72 65 28 67 6d 78 2c | 6d 74 29 3b 0a 58 23 65 |ore(gmx,|mt);.X#e|
|00001f60| 6c 73 65 0a 58 20 20 20 | 72 65 74 75 72 6e 20 67 |lse.X |return g|
|00001f70| 6d 2d 3e 42 6f 72 72 6f | 77 53 74 6f 72 65 28 67 |m->Borro|wStore(g|
|00001f80| 6d 78 2c 6d 74 29 3b 0a | 58 23 65 6e 64 69 66 0a |mx,mt);.|X#endif.|
|00001f90| 58 7d 0a 58 0a 58 47 65 | 6e 65 72 61 6c 4d 61 74 |X}.X.XGe|neralMat|
|00001fa0| 72 69 78 2a 20 43 6f 6c | 65 64 4d 61 74 72 69 78 |rix* Col|edMatrix|
|00001fb0| 3a 3a 45 76 61 6c 75 61 | 74 65 28 4d 61 74 72 69 |::Evalua|te(Matri|
|00001fc0| 78 54 79 70 65 20 6d 74 | 29 0a 58 7b 0a 58 20 20 |xType mt|).X{.X |
|00001fd0| 20 67 6d 20 3d 20 28 28 | 42 61 73 65 4d 61 74 72 | gm = ((|BaseMatr|
|00001fe0| 69 78 2a 26 29 62 6d 29 | 2d 3e 45 76 61 6c 75 61 |ix*&)bm)|->Evalua|
|00001ff0| 74 65 28 29 3b 0a 58 20 | 20 20 47 65 6e 65 72 61 |te();.X | Genera|
|00002000| 6c 4d 61 74 72 69 78 2a | 20 67 6d 78 20 3d 20 6e |lMatrix*| gmx = n|
|00002010| 65 77 20 43 6f 6c 75 6d | 6e 56 65 63 74 6f 72 3b |ew Colum|nVector;|
|00002020| 20 4d 61 74 72 69 78 45 | 72 72 6f 72 4e 6f 53 70 | MatrixE|rrorNoSp|
|00002030| 61 63 65 28 67 6d 78 29 | 3b 0a 58 20 20 20 67 6d |ace(gmx)|;.X gm|
|00002040| 78 2d 3e 6e 63 6f 6c 73 | 20 3d 20 31 3b 20 67 6d |x->ncols| = 1; gm|
|00002050| 78 2d 3e 6e 72 6f 77 73 | 20 3d 20 67 6d 78 2d 3e |x->nrows| = gmx->|
|00002060| 73 74 6f 72 61 67 65 20 | 3d 20 67 6d 2d 3e 73 74 |storage |= gm->st|
|00002070| 6f 72 61 67 65 3b 0a 58 | 23 69 66 64 65 66 20 54 |orage;.X|#ifdef T|
|00002080| 45 4d 50 53 5f 44 45 53 | 54 52 4f 59 45 44 5f 51 |EMPS_DES|TROYED_Q|
|00002090| 55 49 43 4b 4c 59 0a 58 | 20 20 20 47 65 6e 65 72 |UICKLY.X| Gener|
|000020a0| 61 6c 4d 61 74 72 69 78 | 2a 20 67 6d 79 20 3d 20 |alMatrix|* gmy = |
|000020b0| 67 6d 3b 20 64 65 6c 65 | 74 65 20 74 68 69 73 3b |gm; dele|te this;|
|000020c0| 20 72 65 74 75 72 6e 20 | 67 6d 79 2d 3e 42 6f 72 | return |gmy->Bor|
|000020d0| 72 6f 77 53 74 6f 72 65 | 28 67 6d 78 2c 6d 74 29 |rowStore|(gmx,mt)|
|000020e0| 3b 0a 58 23 65 6c 73 65 | 0a 58 20 20 20 72 65 74 |;.X#else|.X ret|
|000020f0| 75 72 6e 20 67 6d 2d 3e | 42 6f 72 72 6f 77 53 74 |urn gm->|BorrowSt|
|00002100| 6f 72 65 28 67 6d 78 2c | 6d 74 29 3b 0a 58 23 65 |ore(gmx,|mt);.X#e|
|00002110| 6e 64 69 66 0a 58 7d 0a | 58 0a 58 47 65 6e 65 72 |ndif.X}.|X.XGener|
|00002120| 61 6c 4d 61 74 72 69 78 | 2a 20 44 69 61 67 65 64 |alMatrix|* Diaged|
|00002130| 4d 61 74 72 69 78 3a 3a | 45 76 61 6c 75 61 74 65 |Matrix::|Evaluate|
|00002140| 28 4d 61 74 72 69 78 54 | 79 70 65 20 6d 74 29 0a |(MatrixT|ype mt).|
|00002150| 58 7b 0a 58 20 20 20 67 | 6d 20 3d 20 28 28 42 61 |X{.X g|m = ((Ba|
|00002160| 73 65 4d 61 74 72 69 78 | 2a 26 29 62 6d 29 2d 3e |seMatrix|*&)bm)->|
|00002170| 45 76 61 6c 75 61 74 65 | 28 29 3b 0a 58 20 20 20 |Evaluate|();.X |
|00002180| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00002190| 6d 78 20 3d 20 6e 65 77 | 20 44 69 61 67 6f 6e 61 |mx = new| Diagona|
|000021a0| 6c 4d 61 74 72 69 78 3b | 20 4d 61 74 72 69 78 45 |lMatrix;| MatrixE|
|000021b0| 72 72 6f 72 4e 6f 53 70 | 61 63 65 28 67 6d 78 29 |rrorNoSp|ace(gmx)|
|000021c0| 3b 0a 58 20 20 20 67 6d | 78 2d 3e 6e 72 6f 77 73 |;.X gm|x->nrows|
|000021d0| 20 3d 20 67 6d 78 2d 3e | 6e 63 6f 6c 73 20 3d 20 | = gmx->|ncols = |
|000021e0| 67 6d 78 2d 3e 73 74 6f | 72 61 67 65 20 3d 20 67 |gmx->sto|rage = g|
|000021f0| 6d 2d 3e 73 74 6f 72 61 | 67 65 3b 0a 58 23 69 66 |m->stora|ge;.X#if|
|00002200| 64 65 66 20 54 45 4d 50 | 53 5f 44 45 53 54 52 4f |def TEMP|S_DESTRO|
|00002210| 59 45 44 5f 51 55 49 43 | 4b 4c 59 0a 58 20 20 20 |YED_QUIC|KLY.X |
|00002220| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 67 |GeneralM|atrix* g|
|00002230| 6d 79 20 3d 20 67 6d 3b | 20 64 65 6c 65 74 65 20 |my = gm;| delete |
|00002240| 74 68 69 73 3b 20 72 65 | 74 75 72 6e 20 67 6d 79 |this; re|turn gmy|
|00002250| 2d 3e 42 6f 72 72 6f 77 | 53 74 6f 72 65 28 67 6d |->Borrow|Store(gm|
|00002260| 78 2c 6d 74 29 3b 0a 58 | 23 65 6c 73 65 0a 58 20 |x,mt);.X|#else.X |
|00002270| 20 20 72 65 74 75 72 6e | 20 67 6d 2d 3e 42 6f 72 | return| gm->Bor|
|00002280| 72 6f 77 53 74 6f 72 65 | 28 67 6d 78 2c 6d 74 29 |rowStore|(gmx,mt)|
|00002290| 3b 0a 58 23 65 6e 64 69 | 66 0a 58 7d 0a 58 0a 58 |;.X#endi|f.X}.X.X|
|000022a0| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 2a 20 4d |GeneralM|atrix* M|
|000022b0| 61 74 65 64 4d 61 74 72 | 69 78 3a 3a 45 76 61 6c |atedMatr|ix::Eval|
|000022c0| 75 61 74 65 28 4d 61 74 | 72 69 78 54 79 70 65 20 |uate(Mat|rixType |
|000022d0| 6d 74 29 0a 58 7b 0a 58 | 20 20 20 54 72 61 63 65 |mt).X{.X| Trace|
|000022e0| 72 20 74 72 28 22 4d 61 | 74 65 64 4d 61 74 72 69 |r tr("Ma|tedMatri|
|000022f0| 78 3a 3a 45 76 61 6c 75 | 61 74 65 22 29 3b 0a 58 |x::Evalu|ate");.X|
|00002300| 20 20 20 67 6d 20 3d 20 | 28 28 42 61 73 65 4d 61 | gm = |((BaseMa|
|00002310| 74 72 69 78 2a 26 29 62 | 6d 29 2d 3e 45 76 61 6c |trix*&)b|m)->Eval|
|00002320| 75 61 74 65 28 29 3b 0a | 58 20 20 20 47 65 6e 65 |uate();.|X Gene|
|00002330| 72 61 6c 4d 61 74 72 69 | 78 2a 20 67 6d 78 20 3d |ralMatri|x* gmx =|
|00002340| 20 6e 65 77 20 4d 61 74 | 72 69 78 3b 20 4d 61 74 | new Mat|rix; Mat|
|00002350| 72 69 78 45 72 72 6f 72 | 4e 6f 53 70 61 63 65 28 |rixError|NoSpace(|
|00002360| 67 6d 78 29 3b 0a 58 20 | 20 20 67 6d 78 2d 3e 6e |gmx);.X | gmx->n|
|00002370| 72 6f 77 73 20 3d 20 6e | 72 3b 20 67 6d 78 2d 3e |rows = n|r; gmx->|
|00002380| 6e 63 6f 6c 73 20 3d 20 | 6e 63 3b 20 67 6d 78 2d |ncols = |nc; gmx-|
|00002390| 3e 73 74 6f 72 61 67 65 | 20 3d 20 67 6d 2d 3e 73 |>storage| = gm->s|
|000023a0| 74 6f 72 61 67 65 3b 0a | 58 20 20 20 69 66 20 28 |torage;.|X if (|
|000023b0| 6e 72 2a 6e 63 20 21 3d | 20 67 6d 78 2d 3e 73 74 |nr*nc !=| gmx->st|
|000023c0| 6f 72 61 67 65 29 0a 58 | 20 20 20 20 20 20 54 68 |orage).X| Th|
|000023d0| 72 6f 77 28 49 6e 63 6f | 6d 70 61 74 69 62 6c 65 |row(Inco|mpatible|
|000023e0| 44 69 6d 65 6e 73 69 6f | 6e 73 45 78 63 65 70 74 |Dimensio|nsExcept|
|000023f0| 69 6f 6e 28 29 29 3b 0a | 58 23 69 66 64 65 66 20 |ion());.|X#ifdef |
|00002400| 54 45 4d 50 53 5f 44 45 | 53 54 52 4f 59 45 44 5f |TEMPS_DE|STROYED_|
|00002410| 51 55 49 43 4b 4c 59 0a | 58 20 20 20 47 65 6e 65 |QUICKLY.|X Gene|
|00002420| 72 61 6c 4d 61 74 72 69 | 78 2a 20 67 6d 79 20 3d |ralMatri|x* gmy =|
|00002430| 20 67 6d 3b 20 64 65 6c | 65 74 65 20 74 68 69 73 | gm; del|ete this|
|00002440| 3b 20 72 65 74 75 72 6e | 20 67 6d 79 2d 3e 42 6f |; return| gmy->Bo|
|00002450| 72 72 6f 77 53 74 6f 72 | 65 28 67 6d 78 2c 6d 74 |rrowStor|e(gmx,mt|
|00002460| 29 3b 0a 58 23 65 6c 73 | 65 0a 58 20 20 20 72 65 |);.X#els|e.X re|
|00002470| 74 75 72 6e 20 67 6d 2d | 3e 42 6f 72 72 6f 77 53 |turn gm-|>BorrowS|
|00002480| 74 6f 72 65 28 67 6d 78 | 2c 6d 74 29 3b 0a 58 23 |tore(gmx|,mt);.X#|
|00002490| 65 6e 64 69 66 0a 58 7d | 0a 58 0a 58 47 65 6e 65 |endif.X}|.X.XGene|
|000024a0| 72 61 6c 4d 61 74 72 69 | 78 2a 20 47 65 74 53 75 |ralMatri|x* GetSu|
|000024b0| 62 4d 61 74 72 69 78 3a | 3a 45 76 61 6c 75 61 74 |bMatrix:|:Evaluat|
|000024c0| 65 28 4d 61 74 72 69 78 | 54 79 70 65 20 6d 74 29 |e(Matrix|Type mt)|
|000024d0| 0a 58 7b 0a 58 20 20 20 | 52 45 50 4f 52 54 0a 58 |.X{.X |REPORT.X|
|000024e0| 20 20 20 54 72 61 63 65 | 72 20 74 72 28 22 53 75 | Trace|r tr("Su|
|000024f0| 62 4d 61 74 72 69 78 28 | 65 76 61 6c 75 61 74 65 |bMatrix(|evaluate|
|00002500| 29 22 29 3b 0a 58 20 20 | 20 67 6d 20 3d 20 28 28 |)");.X | gm = ((|
|00002510| 42 61 73 65 4d 61 74 72 | 69 78 2a 26 29 62 6d 29 |BaseMatr|ix*&)bm)|
|00002520| 2d 3e 45 76 61 6c 75 61 | 74 65 28 29 3b 0a 58 20 |->Evalua|te();.X |
|00002530| 20 20 69 66 20 28 72 6f | 77 5f 6e 75 6d 62 65 72 | if (ro|w_number|
|00002540| 20 3c 20 30 29 20 72 6f | 77 5f 6e 75 6d 62 65 72 | < 0) ro|w_number|
|00002550| 20 3d 20 67 6d 2d 3e 4e | 72 6f 77 73 28 29 3b 0a | = gm->N|rows();.|
|00002560| 58 20 20 20 69 66 20 28 | 63 6f 6c 5f 6e 75 6d 62 |X if (|col_numb|
|00002570| 65 72 20 3c 20 30 29 20 | 63 6f 6c 5f 6e 75 6d 62 |er < 0) |col_numb|
|00002580| 65 72 20 3d 20 67 6d 2d | 3e 4e 63 6f 6c 73 28 29 |er = gm-|>Ncols()|
|00002590| 3b 0a 58 20 20 20 69 66 | 20 28 72 6f 77 5f 73 6b |;.X if| (row_sk|
|000025a0| 69 70 2b 72 6f 77 5f 6e | 75 6d 62 65 72 20 3e 20 |ip+row_n|umber > |
|000025b0| 67 6d 2d 3e 4e 72 6f 77 | 73 28 29 20 7c 7c 20 63 |gm->Nrow|s() || c|
|000025c0| 6f 6c 5f 73 6b 69 70 2b | 63 6f 6c 5f 6e 75 6d 62 |ol_skip+|col_numb|
|000025d0| 65 72 20 3e 20 67 6d 2d | 3e 4e 63 6f 6c 73 28 29 |er > gm-|>Ncols()|
|000025e0| 29 0a 58 20 20 20 20 20 | 20 54 68 72 6f 77 28 53 |).X | Throw(S|
|000025f0| 75 62 4d 61 74 72 69 78 | 44 69 6d 65 6e 73 69 6f |ubMatrix|Dimensio|
|00002600| 6e 45 78 63 65 70 74 69 | 6f 6e 28 29 29 3b 0a 58 |nExcepti|on());.X|
|00002610| 20 20 20 69 66 20 28 49 | 73 53 79 6d 29 20 43 6f | if (I|sSym) Co|
|00002620| 6d 70 61 72 65 28 67 6d | 2d 3e 54 79 70 65 28 29 |mpare(gm|->Type()|
|00002630| 2e 73 73 75 62 28 29 2c | 20 6d 74 29 3b 20 0a 58 |.ssub(),| mt); .X|
|00002640| 20 20 20 65 6c 73 65 20 | 43 6f 6d 70 61 72 65 28 | else |Compare(|
|00002650| 67 6d 2d 3e 54 79 70 65 | 28 29 2e 73 75 62 28 29 |gm->Type|().sub()|
|00002660| 2c 20 6d 74 29 3b 0a 58 | 20 20 20 47 65 6e 65 72 |, mt);.X| Gener|
|00002670| 61 6c 4d 61 74 72 69 78 | 2a 20 67 6d 78 20 3d 20 |alMatrix|* gmx = |
|00002680| 6d 74 2e 4e 65 77 28 72 | 6f 77 5f 6e 75 6d 62 65 |mt.New(r|ow_numbe|
|00002690| 72 2c 20 63 6f 6c 5f 6e | 75 6d 62 65 72 2c 20 74 |r, col_n|umber, t|
|000026a0| 68 69 73 29 3b 0a 58 20 | 20 20 69 6e 74 20 69 20 |his);.X | int i |
|000026b0| 3d 20 72 6f 77 5f 6e 75 | 6d 62 65 72 3b 0a 58 20 |= row_nu|mber;.X |
|000026c0| 20 20 4d 61 74 72 69 78 | 52 6f 77 20 6d 72 28 67 | Matrix|Row mr(g|
|000026d0| 6d 2c 20 4c 6f 61 64 4f | 6e 45 6e 74 72 79 2c 20 |m, LoadO|nEntry, |
|000026e0| 72 6f 77 5f 73 6b 69 70 | 29 3b 20 0a 58 20 20 20 |row_skip|); .X |
|000026f0| 4d 61 74 72 69 78 52 6f | 77 20 6d 72 78 28 67 6d |MatrixRo|w mrx(gm|
|00002700| 78 2c 20 53 74 6f 72 65 | 4f 6e 45 78 69 74 2b 44 |x, Store|OnExit+D|
|00002710| 69 72 65 63 74 50 61 72 | 74 29 3b 0a 58 20 20 20 |irectPar|t);.X |
|00002720| 4d 61 74 72 69 78 52 6f | 77 43 6f 6c 20 73 75 62 |MatrixRo|wCol sub|
|00002730| 3b 0a 58 20 20 20 77 68 | 69 6c 65 20 28 69 2d 2d |;.X wh|ile (i--|
|00002740| 29 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 6d |).X {.|X m|
|00002750| 72 2e 53 75 62 52 6f 77 | 43 6f 6c 28 73 75 62 2c |r.SubRow|Col(sub,|
|00002760| 20 63 6f 6c 5f 73 6b 69 | 70 2c 20 63 6f 6c 5f 6e | col_ski|p, col_n|
|00002770| 75 6d 62 65 72 29 3b 20 | 20 20 2f 2f 20 70 75 74 |umber); | // put|
|00002780| 20 76 61 6c 75 65 73 20 | 69 6e 20 73 75 62 0a 58 | values |in sub.X|
|00002790| 20 20 20 20 20 20 6d 72 | 78 2e 43 6f 70 79 28 73 | mr|x.Copy(s|
|000027a0| 75 62 29 3b 20 6d 72 78 | 2e 4e 65 78 74 28 29 3b |ub); mrx|.Next();|
|000027b0| 20 6d 72 2e 4e 65 78 74 | 28 29 3b 0a 58 20 20 20 | mr.Next|();.X |
|000027c0| 7d 0a 58 20 20 20 67 6d | 78 2d 3e 52 65 6c 65 61 |}.X gm|x->Relea|
|000027d0| 73 65 41 6e 64 44 65 6c | 65 74 65 28 29 3b 20 67 |seAndDel|ete(); g|
|000027e0| 6d 2d 3e 74 44 65 6c 65 | 74 65 28 29 3b 0a 58 23 |m->tDele|te();.X#|
|000027f0| 69 66 64 65 66 20 54 45 | 4d 50 53 5f 44 45 53 54 |ifdef TE|MPS_DEST|
|00002800| 52 4f 59 45 44 5f 51 55 | 49 43 4b 4c 59 0a 58 20 |ROYED_QU|ICKLY.X |
|00002810| 20 20 64 65 6c 65 74 65 | 20 74 68 69 73 3b 0a 58 | delete| this;.X|
|00002820| 23 65 6e 64 69 66 0a 58 | 20 20 20 72 65 74 75 72 |#endif.X| retur|
|00002830| 6e 20 67 6d 78 3b 0a 58 | 7d 20 20 20 0a 58 0a 58 |n gmx;.X|} .X.X|
|00002840| 0a 58 47 65 6e 65 72 61 | 6c 4d 61 74 72 69 78 2a |.XGenera|lMatrix*|
|00002850| 20 52 65 74 75 72 6e 4d | 61 74 72 69 78 58 3a 3a | ReturnM|atrixX::|
|00002860| 45 76 61 6c 75 61 74 65 | 28 4d 61 74 72 69 78 54 |Evaluate|(MatrixT|
|00002870| 79 70 65 20 6d 74 29 0a | 58 7b 0a 58 23 69 66 64 |ype mt).|X{.X#ifd|
|00002880| 65 66 20 54 45 4d 50 53 | 5f 44 45 53 54 52 4f 59 |ef TEMPS|_DESTROY|
|00002890| 45 44 5f 51 55 49 43 4b | 4c 59 0a 58 20 20 20 47 |ED_QUICK|LY.X G|
|000028a0| 65 6e 65 72 61 6c 4d 61 | 74 72 69 78 2a 20 67 6d |eneralMa|trix* gm|
|000028b0| 78 20 3d 20 67 6d 3b 20 | 64 65 6c 65 74 65 20 74 |x = gm; |delete t|
|000028c0| 68 69 73 3b 20 72 65 74 | 75 72 6e 20 67 6d 78 2d |his; ret|urn gmx-|
|000028d0| 3e 45 76 61 6c 75 61 74 | 65 28 6d 74 29 3b 0a 58 |>Evaluat|e(mt);.X|
|000028e0| 23 65 6c 73 65 0a 58 20 | 20 20 72 65 74 75 72 6e |#else.X | return|
|000028f0| 20 67 6d 2d 3e 45 76 61 | 6c 75 61 74 65 28 6d 74 | gm->Eva|luate(mt|
|00002900| 29 3b 0a 58 23 65 6e 64 | 69 66 0a 58 7d 0a 58 0a |);.X#end|if.X}.X.|
|00002910| 58 0a 58 76 6f 69 64 20 | 47 65 6e 65 72 61 6c 4d |X.Xvoid |GeneralM|
|00002920| 61 74 72 69 78 3a 3a 41 | 64 64 28 47 65 6e 65 72 |atrix::A|dd(Gener|
|00002930| 61 6c 4d 61 74 72 69 78 | 2a 20 67 6d 31 2c 20 52 |alMatrix|* gm1, R|
|00002940| 65 61 6c 20 66 29 0a 58 | 7b 0a 58 20 20 20 52 45 |eal f).X|{.X RE|
|00002950| 50 4f 52 54 0a 58 20 20 | 20 52 65 61 6c 2a 20 73 |PORT.X | Real* s|
|00002960| 31 3d 67 6d 31 2d 3e 73 | 74 6f 72 65 3b 20 52 65 |1=gm1->s|tore; Re|
|00002970| 61 6c 2a 20 73 3d 73 74 | 6f 72 65 3b 20 69 6e 74 |al* s=st|ore; int|
|00002980| 20 69 3d 28 73 74 6f 72 | 61 67 65 20 3e 3e 20 32 | i=(stor|age >> 2|
|00002990| 29 3b 0a 58 20 20 20 77 | 68 69 6c 65 20 28 69 2d |);.X w|hile (i-|
|000029a0| 2d 29 0a 58 20 20 20 7b | 20 2a 73 2b 2b 20 3d 20 |-).X {| *s++ = |
|000029b0| 2a 73 31 2b 2b 20 2b 20 | 66 3b 20 2a 73 2b 2b 20 |*s1++ + |f; *s++ |
|000029c0| 3d 20 2a 73 31 2b 2b 20 | 2b 20 66 3b 20 2a 73 2b |= *s1++ |+ f; *s+|
|000029d0| 2b 20 3d 20 2a 73 31 2b | 2b 20 2b 20 66 3b 20 2a |+ = *s1+|+ + f; *|
|000029e0| 73 2b 2b 20 3d 20 2a 73 | 31 2b 2b 20 2b 20 66 3b |s++ = *s|1++ + f;|
|000029f0| 20 7d 0a 58 20 20 20 69 | 20 3d 20 73 74 6f 72 61 | }.X i| = stora|
|00002a00| 67 65 20 26 20 33 3b 20 | 77 68 69 6c 65 20 28 69 |ge & 3; |while (i|
|00002a10| 2d 2d 29 20 2a 73 2b 2b | 20 3d 20 2a 73 31 2b 2b |--) *s++| = *s1++|
|00002a20| 20 2b 20 66 3b 0a 58 7d | 0a 58 20 20 20 0a 58 76 | + f;.X}|.X .Xv|
|00002a30| 6f 69 64 20 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |oid Gene|ralMatri|
|00002a40| 78 3a 3a 41 64 64 28 52 | 65 61 6c 20 66 29 0a 58 |x::Add(R|eal f).X|
|00002a50| 7b 0a 58 20 20 20 52 45 | 50 4f 52 54 0a 58 20 20 |{.X RE|PORT.X |
|00002a60| 20 52 65 61 6c 2a 20 73 | 3d 73 74 6f 72 65 3b 20 | Real* s|=store; |
|00002a70| 69 6e 74 20 69 3d 28 73 | 74 6f 72 61 67 65 20 3e |int i=(s|torage >|
|00002a80| 3e 20 32 29 3b 0a 58 20 | 20 20 77 68 69 6c 65 20 |> 2);.X | while |
|00002a90| 28 69 2d 2d 29 20 7b 20 | 2a 73 2b 2b 20 2b 3d 20 |(i--) { |*s++ += |
|00002aa0| 66 3b 20 2a 73 2b 2b 20 | 2b 3d 20 66 3b 20 2a 73 |f; *s++ |+= f; *s|
|00002ab0| 2b 2b 20 2b 3d 20 66 3b | 20 2a 73 2b 2b 20 2b 3d |++ += f;| *s++ +=|
|00002ac0| 20 66 3b 20 7d 0a 58 20 | 20 20 69 20 3d 20 73 74 | f; }.X | i = st|
|00002ad0| 6f 72 61 67 65 20 26 20 | 33 3b 20 77 68 69 6c 65 |orage & |3; while|
|00002ae0| 20 28 69 2d 2d 29 20 2a | 73 2b 2b 20 2b 3d 20 66 | (i--) *|s++ += f|
|00002af0| 3b 0a 58 7d 0a 58 20 20 | 20 0a 58 76 6f 69 64 20 |;.X}.X | .Xvoid |
|00002b00| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 3a 3a 4e |GeneralM|atrix::N|
|00002b10| 65 67 61 74 65 28 47 65 | 6e 65 72 61 6c 4d 61 74 |egate(Ge|neralMat|
|00002b20| 72 69 78 2a 20 67 6d 31 | 29 0a 58 7b 0a 58 20 20 |rix* gm1|).X{.X |
|00002b30| 20 2f 2f 20 63 68 61 6e | 67 65 20 73 69 67 6e 20 | // chan|ge sign |
|00002b40| 6f 66 20 65 6c 65 6d 65 | 6e 74 73 0a 58 20 20 20 |of eleme|nts.X |
|00002b50| 52 45 50 4f 52 54 0a 58 | 20 20 20 52 65 61 6c 2a |REPORT.X| Real*|
|00002b60| 20 73 31 3d 67 6d 31 2d | 3e 73 74 6f 72 65 3b 20 | s1=gm1-|>store; |
|00002b70| 52 65 61 6c 2a 20 73 3d | 73 74 6f 72 65 3b 20 69 |Real* s=|store; i|
|00002b80| 6e 74 20 69 3d 28 73 74 | 6f 72 61 67 65 20 3e 3e |nt i=(st|orage >>|
|00002b90| 20 32 29 3b 0a 58 20 20 | 20 77 68 69 6c 65 20 28 | 2);.X | while (|
|00002ba0| 69 2d 2d 29 0a 58 20 20 | 20 7b 20 2a 73 2b 2b 20 |i--).X | { *s++ |
|00002bb0| 3d 20 2d 28 2a 73 31 2b | 2b 29 3b 20 2a 73 2b 2b |= -(*s1+|+); *s++|
|00002bc0| 20 3d 20 2d 28 2a 73 31 | 2b 2b 29 3b 20 2a 73 2b | = -(*s1|++); *s+|
|00002bd0| 2b 20 3d 20 2d 28 2a 73 | 31 2b 2b 29 3b 20 2a 73 |+ = -(*s|1++); *s|
|00002be0| 2b 2b 20 3d 20 2d 28 2a | 73 31 2b 2b 29 3b 20 7d |++ = -(*|s1++); }|
|00002bf0| 0a 58 20 20 20 69 20 3d | 20 73 74 6f 72 61 67 65 |.X i =| storage|
|00002c00| 20 26 20 33 3b 20 77 68 | 69 6c 65 28 69 2d 2d 29 | & 3; wh|ile(i--)|
|00002c10| 20 2a 73 2b 2b 20 3d 20 | 2d 28 2a 73 31 2b 2b 29 | *s++ = |-(*s1++)|
|00002c20| 3b 0a 58 7d 0a 58 20 20 | 20 0a 58 76 6f 69 64 20 |;.X}.X | .Xvoid |
|00002c30| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 3a 3a 4e |GeneralM|atrix::N|
|00002c40| 65 67 61 74 65 28 29 0a | 58 7b 0a 58 20 20 20 52 |egate().|X{.X R|
|00002c50| 45 50 4f 52 54 0a 58 20 | 20 20 52 65 61 6c 2a 20 |EPORT.X | Real* |
|00002c60| 73 3d 73 74 6f 72 65 3b | 20 69 6e 74 20 69 3d 28 |s=store;| int i=(|
|00002c70| 73 74 6f 72 61 67 65 20 | 3e 3e 20 32 29 3b 0a 58 |storage |>> 2);.X|
|00002c80| 20 20 20 77 68 69 6c 65 | 20 28 69 2d 2d 29 0a 58 | while| (i--).X|
|00002c90| 20 20 20 7b 20 2a 73 20 | 3d 20 2d 28 2a 73 29 3b | { *s |= -(*s);|
|00002ca0| 20 73 2b 2b 3b 20 2a 73 | 20 3d 20 2d 28 2a 73 29 | s++; *s| = -(*s)|
|00002cb0| 3b 20 73 2b 2b 3b 20 2a | 73 20 3d 20 2d 28 2a 73 |; s++; *|s = -(*s|
|00002cc0| 29 3b 20 73 2b 2b 3b 20 | 2a 73 20 3d 20 2d 28 2a |); s++; |*s = -(*|
|00002cd0| 73 29 3b 20 73 2b 2b 3b | 20 7d 0a 58 20 20 20 69 |s); s++;| }.X i|
|00002ce0| 20 3d 20 73 74 6f 72 61 | 67 65 20 26 20 33 3b 20 | = stora|ge & 3; |
|00002cf0| 77 68 69 6c 65 28 69 2d | 2d 29 20 7b 20 2a 73 20 |while(i-|-) { *s |
|00002d00| 3d 20 2d 28 2a 73 29 3b | 20 73 2b 2b 3b 20 7d 0a |= -(*s);| s++; }.|
|00002d10| 58 7d 0a 58 20 20 20 0a | 58 76 6f 69 64 20 47 65 |X}.X .|Xvoid Ge|
|00002d20| 6e 65 72 61 6c 4d 61 74 | 72 69 78 3a 3a 4d 75 6c |neralMat|rix::Mul|
|00002d30| 74 69 70 6c 79 28 47 65 | 6e 65 72 61 6c 4d 61 74 |tiply(Ge|neralMat|
|00002d40| 72 69 78 2a 20 67 6d 31 | 2c 20 52 65 61 6c 20 66 |rix* gm1|, Real f|
|00002d50| 29 0a 58 7b 0a 58 20 20 | 20 52 45 50 4f 52 54 0a |).X{.X | REPORT.|
|00002d60| 58 20 20 20 52 65 61 6c | 2a 20 73 31 3d 67 6d 31 |X Real|* s1=gm1|
|00002d70| 2d 3e 73 74 6f 72 65 3b | 20 52 65 61 6c 2a 20 73 |->store;| Real* s|
|00002d80| 3d 73 74 6f 72 65 3b 20 | 20 69 6e 74 20 69 3d 28 |=store; | int i=(|
|00002d90| 73 74 6f 72 61 67 65 20 | 3e 3e 20 32 29 3b 0a 58 |storage |>> 2);.X|
|00002da0| 20 20 20 77 68 69 6c 65 | 20 28 69 2d 2d 29 0a 58 | while| (i--).X|
|00002db0| 20 20 20 7b 20 2a 73 2b | 2b 20 3d 20 2a 73 31 2b | { *s+|+ = *s1+|
|00002dc0| 2b 20 2a 20 66 3b 20 2a | 73 2b 2b 20 3d 20 2a 73 |+ * f; *|s++ = *s|
|00002dd0| 31 2b 2b 20 2a 20 66 3b | 20 2a 73 2b 2b 20 3d 20 |1++ * f;| *s++ = |
|00002de0| 2a 73 31 2b 2b 20 2a 20 | 66 3b 20 2a 73 2b 2b 20 |*s1++ * |f; *s++ |
|00002df0| 3d 20 2a 73 31 2b 2b 20 | 2a 20 66 3b 20 7d 0a 58 |= *s1++ |* f; }.X|
|00002e00| 20 20 20 69 20 3d 20 73 | 74 6f 72 61 67 65 20 26 | i = s|torage &|
|00002e10| 20 33 3b 20 77 68 69 6c | 65 20 28 69 2d 2d 29 20 | 3; whil|e (i--) |
|00002e20| 2a 73 2b 2b 20 3d 20 2a | 73 31 2b 2b 20 2a 20 66 |*s++ = *|s1++ * f|
|00002e30| 3b 0a 58 7d 0a 58 20 20 | 20 0a 58 76 6f 69 64 20 |;.X}.X | .Xvoid |
|00002e40| 47 65 6e 65 72 61 6c 4d | 61 74 72 69 78 3a 3a 4d |GeneralM|atrix::M|
|00002e50| 75 6c 74 69 70 6c 79 28 | 52 65 61 6c 20 66 29 0a |ultiply(|Real f).|
|00002e60| 58 7b 0a 58 20 20 20 52 | 45 50 4f 52 54 0a 58 20 |X{.X R|EPORT.X |
|00002e70| 20 20 52 65 61 6c 2a 20 | 73 3d 73 74 6f 72 65 3b | Real* |s=store;|
|00002e80| 20 69 6e 74 20 69 3d 28 | 73 74 6f 72 61 67 65 20 | int i=(|storage |
|00002e90| 3e 3e 20 32 29 3b 0a 58 | 20 20 20 77 68 69 6c 65 |>> 2);.X| while|
|00002ea0| 20 28 69 2d 2d 29 20 7b | 20 2a 73 2b 2b 20 2a 3d | (i--) {| *s++ *=|
|00002eb0| 20 66 3b 20 2a 73 2b 2b | 20 2a 3d 20 66 3b 20 2a | f; *s++| *= f; *|
|00002ec0| 73 2b 2b 20 2a 3d 20 66 | 3b 20 2a 73 2b 2b 20 2a |s++ *= f|; *s++ *|
|00002ed0| 3d 20 66 3b 20 7d 0a 58 | 20 20 20 69 20 3d 20 73 |= f; }.X| i = s|
|00002ee0| 74 6f 72 61 67 65 20 26 | 20 33 3b 20 77 68 69 6c |torage &| 3; whil|
|00002ef0| 65 20 28 69 2d 2d 29 20 | 2a 73 2b 2b 20 2a 3d 20 |e (i--) |*s++ *= |
|00002f00| 66 3b 0a 58 7d 0a 58 20 | 20 20 0a 58 0a 58 2f 2a |f;.X}.X | .X.X/*|
|00002f10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002f20| 2a 2a 2a 2a 2a 2a 2a 20 | 4d 61 74 72 69 78 49 6e |******* |MatrixIn|
|00002f30| 70 75 74 20 72 6f 75 74 | 69 6e 65 73 20 2a 2a 2a |put rout|ines ***|
|00002f40| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002f50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2f 0a 58 0a 58 69 6e |********|*/.X.Xin|
|00002f60| 74 20 4d 61 74 72 69 78 | 49 6e 70 75 74 3a 3a 6e |t Matrix|Input::n|
|00002f70| 20 3d 20 30 3b 20 20 20 | 20 20 20 20 2f 2f 20 6e | = 0; | // n|
|00002f80| 75 6d 62 65 72 20 76 61 | 6c 75 65 73 20 73 74 69 |umber va|lues sti|
|00002f90| 6c 6c 20 74 6f 20 62 65 | 20 72 65 61 64 0a 58 52 |ll to be| read.XR|
|00002fa0| 65 61 6c 2a 20 4d 61 74 | 72 69 78 49 6e 70 75 74 |eal* Mat|rixInput|
|00002fb0| 3a 3a 72 3b 20 20 20 20 | 20 20 20 20 20 2f 2f 20 |::r; | // |
|00002fc0| 70 6f 69 6e 74 65 72 20 | 74 6f 20 6c 61 73 74 20 |pointer |to last |
|00002fd0| 74 68 69 6e 67 20 72 65 | 61 64 0a 58 69 6e 74 20 |thing re|ad.Xint |
|00002fe0| 4d 61 74 72 69 78 49 6e | 70 75 74 3a 3a 64 65 70 |MatrixIn|put::dep|
|00002ff0| 74 68 20 3d 20 30 3b 20 | 20 20 2f 2f 20 6e 75 6d |th = 0; | // num|
|00003000| 62 65 72 20 6f 66 20 6f | 62 6a 65 63 74 73 20 6f |ber of o|bjects o|
|00003010| 66 20 74 68 69 73 20 63 | 6c 61 73 73 20 69 6e 20 |f this c|lass in |
|00003020| 65 78 69 73 74 65 6e 63 | 65 0a 58 0a 58 4d 61 74 |existenc|e.X.XMat|
|00003030| 72 69 78 49 6e 70 75 74 | 20 4d 61 74 72 69 78 49 |rixInput| MatrixI|
|00003040| 6e 70 75 74 3a 3a 6f 70 | 65 72 61 74 6f 72 3c 3c |nput::op|erator<<|
|00003050| 28 52 65 61 6c 20 66 29 | 0a 58 7b 0a 58 20 20 20 |(Real f)|.X{.X |
|00003060| 69 66 20 28 21 28 6e 2d | 2d 29 29 0a 58 20 20 20 |if (!(n-|-)).X |
|00003070| 7b 20 64 65 70 74 68 3d | 30 3b 20 20 6e 3d 30 3b |{ depth=|0; n=0;|
|00003080| 20 54 68 72 6f 77 28 50 | 72 6f 67 72 61 6d 45 78 | Throw(P|rogramEx|
|00003090| 63 65 70 74 69 6f 6e 28 | 22 4c 69 73 74 20 6f 66 |ception(|"List of|
|000030a0| 20 76 61 6c 75 65 73 20 | 74 6f 6f 20 6c 6f 6e 67 | values |too long|
|000030b0| 22 29 29 3b 20 7d 0a 58 | 20 20 20 2a 28 2b 2b 72 |")); }.X| *(++r|
|000030c0| 29 20 3d 20 66 3b 0a 58 | 20 20 20 72 65 74 75 72 |) = f;.X| retur|
|000030d0| 6e 20 4d 61 74 72 69 78 | 49 6e 70 75 74 28 29 3b |n Matrix|Input();|
|000030e0| 0a 58 7d 0a 58 0a 58 0a | 58 4d 61 74 72 69 78 49 |.X}.X.X.|XMatrixI|
|000030f0| 6e 70 75 74 20 42 61 6e | 64 4d 61 74 72 69 78 3a |nput Ban|dMatrix:|
|00003100| 3a 6f 70 65 72 61 74 6f | 72 3c 3c 28 52 65 61 6c |:operato|r<<(Real|
|00003110| 29 0a 58 7b 0a 58 20 20 | 20 54 68 72 6f 77 28 50 |).X{.X | Throw(P|
|00003120| 72 6f 67 72 61 6d 45 78 | 63 65 70 74 69 6f 6e 28 |rogramEx|ception(|
|00003130| 22 43 61 6e 6e 6f 74 20 | 75 73 65 20 6c 69 73 74 |"Cannot |use list|
|00003140| 20 72 65 61 64 20 77 69 | 74 68 20 61 20 42 61 6e | read wi|th a Ban|
|00003150| 64 4d 61 74 72 69 78 22 | 29 29 3b 0a 58 20 20 20 |dMatrix"|));.X |
|00003160| 72 65 74 75 72 6e 20 4d | 61 74 72 69 78 49 6e 70 |return M|atrixInp|
|00003170| 75 74 28 29 3b 20 0a 58 | 7d 0a 58 0a 58 76 6f 69 |ut(); .X|}.X.Xvoi|
|00003180| 64 20 42 61 6e 64 4d 61 | 74 72 69 78 3a 3a 6f 70 |d BandMa|trix::op|
|00003190| 65 72 61 74 6f 72 3c 3c | 28 63 6f 6e 73 74 20 52 |erator<<|(const R|
|000031a0| 65 61 6c 2a 29 0a 58 7b | 20 54 68 72 6f 77 28 50 |eal*).X{| Throw(P|
|000031b0| 72 6f 67 72 61 6d 45 78 | 63 65 70 74 69 6f 6e 28 |rogramEx|ception(|
|000031c0| 22 43 61 6e 6e 6f 74 20 | 75 73 65 20 61 72 72 61 |"Cannot |use arra|
|000031d0| 79 20 72 65 61 64 20 77 | 69 74 68 20 61 20 42 61 |y read w|ith a Ba|
|000031e0| 6e 64 4d 61 74 72 69 78 | 22 29 29 3b 20 7d 0a 58 |ndMatrix|")); }.X|
|000031f0| 0a 58 4d 61 74 72 69 78 | 49 6e 70 75 74 20 47 65 |.XMatrix|Input Ge|
|00003200| 6e 65 72 61 6c 4d 61 74 | 72 69 78 3a 3a 6f 70 65 |neralMat|rix::ope|
|00003210| 72 61 74 6f 72 3c 3c 28 | 52 65 61 6c 20 66 29 0a |rator<<(|Real f).|
|00003220| 58 7b 0a 58 20 20 20 69 | 66 20 28 4d 61 74 72 69 |X{.X i|f (Matri|
|00003230| 78 49 6e 70 75 74 3a 3a | 6e 29 0a 58 20 20 20 7b |xInput::|n).X {|
|00003240| 0a 58 20 20 20 20 20 20 | 4d 61 74 72 69 78 49 6e |.X |MatrixIn|
|00003250| 70 75 74 3a 3a 64 65 70 | 74 68 3d 30 3b 20 20 20 |put::dep|th=0; |
|00003260| 20 20 20 20 20 20 20 20 | 20 2f 2f 20 73 6f 20 77 | | // so w|
|00003270| 65 20 63 61 6e 20 72 65 | 63 6f 76 65 72 0a 58 20 |e can re|cover.X |
|00003280| 20 20 20 20 20 4d 61 74 | 72 69 78 49 6e 70 75 74 | Mat|rixInput|
|00003290| 3a 3a 6e 3d 30 3b 20 20 | 20 20 20 20 20 20 20 20 |::n=0; | |
|000032a0| 20 20 20 20 20 20 2f 2f | 20 73 6f 20 77 65 20 63 | //| so we c|
|000032b0| 61 6e 20 72 65 63 6f 76 | 65 72 0a 58 20 20 20 20 |an recov|er.X |
|000032c0| 20 20 54 68 72 6f 77 28 | 50 72 6f 67 72 61 6d 45 | Throw(|ProgramE|
|000032d0| 78 63 65 70 74 69 6f 6e | 28 22 41 20 6c 69 73 74 |xception|("A list|
|000032e0| 20 6f 66 20 76 61 6c 75 | 65 73 20 77 61 73 20 74 | of valu|es was t|
|000032f0| 6f 6f 20 73 68 6f 72 74 | 22 29 29 3b 0a 58 20 20 |oo short|"));.X |
|00003300| 20 7d 0a 58 20 20 20 4d | 61 74 72 69 78 49 6e 70 | }.X M|atrixInp|
|00003310| 75 74 3a 3a 6e 20 3d 20 | 53 74 6f 72 61 67 65 28 |ut::n = |Storage(|
|00003320| 29 3b 0a 58 20 20 20 69 | 66 20 28 4d 61 74 72 69 |);.X i|f (Matri|
|00003330| 78 49 6e 70 75 74 3a 3a | 6e 3c 3d 30 29 0a 58 20 |xInput::|n<=0).X |
|00003340| 20 20 20 20 20 54 68 72 | 6f 77 28 50 72 6f 67 72 | Thr|ow(Progr|
|00003350| 61 6d 45 78 63 65 70 74 | 69 6f 6e 28 22 4c 6f 61 |amExcept|ion("Loa|
|00003360| 64 69 6e 67 20 64 61 74 | 61 20 74 6f 20 7a 65 72 |ding dat|a to zer|
|00003370| 6f 20 64 69 6d 65 6e 73 | 69 6f 6e 20 6d 61 74 72 |o dimens|ion matr|
|00003380| 69 78 22 29 29 3b 0a 58 | 20 20 20 4d 61 74 72 69 |ix"));.X| Matri|
|00003390| 78 49 6e 70 75 74 3a 3a | 72 20 3d 20 53 74 6f 72 |xInput::|r = Stor|
|000033a0| 65 28 29 3b 20 2a 28 4d | 61 74 72 69 78 49 6e 70 |e(); *(M|atrixInp|
|000033b0| 75 74 3a 3a 72 29 20 3d | 20 66 3b 20 4d 61 74 72 |ut::r) =| f; Matr|
|000033c0| 69 78 49 6e 70 75 74 3a | 3a 6e 2d 2d 3b 0a 58 20 |ixInput:|:n--;.X |
|000033d0| 20 20 72 65 74 75 72 6e | 20 4d 61 74 72 69 78 49 | return| MatrixI|
|000033e0| 6e 70 75 74 28 29 3b 20 | 0a 58 7d 0a 58 0a 58 4d |nput(); |.X}.X.XM|
|000033f0| 61 74 72 69 78 49 6e 70 | 75 74 3a 3a 7e 4d 61 74 |atrixInp|ut::~Mat|
|00003400| 72 69 78 49 6e 70 75 74 | 28 29 0a 58 7b 0a 58 20 |rixInput|().X{.X |
|00003410| 20 20 69 66 20 28 64 65 | 70 74 68 20 26 26 20 2d | if (de|pth && -|
|00003420| 2d 64 65 70 74 68 20 3d | 3d 20 30 20 26 26 20 6e |-depth =|= 0 && n|
|00003430| 20 21 3d 20 30 29 0a 58 | 20 20 20 7b 0a 58 20 20 | != 0).X| {.X |
|00003440| 20 20 20 20 64 65 70 74 | 68 20 3d 20 30 3b 20 6e | dept|h = 0; n|
|00003450| 20 3d 20 30 3b 20 20 20 | 20 20 20 20 20 20 20 20 | = 0; | |
|00003460| 20 20 20 20 20 2f 2f 20 | 73 6f 20 77 65 20 63 61 | // |so we ca|
|00003470| 6e 20 72 65 63 6f 76 65 | 72 0a 58 20 20 20 20 20 |n recove|r.X |
|00003480| 20 54 68 72 6f 77 28 50 | 72 6f 67 72 61 6d 45 78 | Throw(P|rogramEx|
|00003490| 63 65 70 74 69 6f 6e 28 | 22 41 20 6c 69 73 74 20 |ception(|"A list |
|000034a0| 6f 66 20 76 61 6c 75 65 | 73 20 77 61 73 20 74 6f |of value|s was to|
|000034b0| 6f 20 73 68 6f 72 74 22 | 29 29 3b 0a 58 20 20 20 |o short"|));.X |
|000034c0| 7d 0a 58 7d 0a 58 0a 45 | 4e 44 5f 4f 46 5f 46 49 |}.X}.X.E|ND_OF_FI|
|000034d0| 4c 45 0a 69 66 20 74 65 | 73 74 20 31 31 35 36 37 |LE.if te|st 11567|
|000034e0| 20 2d 6e 65 20 60 77 63 | 20 2d 63 20 3c 27 6e 65 | -ne `wc| -c <'ne|
|000034f0| 77 6d 61 74 35 2e 63 78 | 78 27 60 3b 20 74 68 65 |wmat5.cx|x'`; the|
|00003500| 6e 0a 20 20 20 20 65 63 | 68 6f 20 73 68 61 72 3a |n. ec|ho shar:|
|00003510| 20 5c 22 27 6e 65 77 6d | 61 74 35 2e 63 78 78 27 | \"'newm|at5.cxx'|
|00003520| 5c 22 20 75 6e 70 61 63 | 6b 65 64 20 77 69 74 68 |\" unpac|ked with|
|00003530| 20 77 72 6f 6e 67 20 73 | 69 7a 65 21 0a 66 69 0a | wrong s|ize!.fi.|
|00003540| 23 20 65 6e 64 20 6f 66 | 20 27 6e 65 77 6d 61 74 |# end of| 'newmat|
|00003550| 35 2e 63 78 78 27 0a 66 | 69 0a 69 66 20 74 65 73 |5.cxx'.f|i.if tes|
|00003560| 74 20 2d 66 20 27 6e 65 | 77 6d 61 74 36 2e 63 78 |t -f 'ne|wmat6.cx|
|00003570| 78 27 20 2d 61 20 22 24 | 7b 31 7d 22 20 21 3d 20 |x' -a "$|{1}" != |
|00003580| 22 2d 63 22 20 3b 20 74 | 68 65 6e 20 0a 20 20 65 |"-c" ; t|hen . e|
|00003590| 63 68 6f 20 73 68 61 72 | 3a 20 57 69 6c 6c 20 6e |cho shar|: Will n|
|000035a0| 6f 74 20 63 6c 6f 62 62 | 65 72 20 65 78 69 73 74 |ot clobb|er exist|
|000035b0| 69 6e 67 20 66 69 6c 65 | 20 5c 22 27 6e 65 77 6d |ing file| \"'newm|
|000035c0| 61 74 36 2e 63 78 78 27 | 5c 22 0a 65 6c 73 65 0a |at6.cxx'|\".else.|
|000035d0| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|000035e0| 63 74 69 6e 67 20 5c 22 | 27 6e 65 77 6d 61 74 36 |cting \"|'newmat6|
|000035f0| 2e 63 78 78 27 5c 22 20 | 5c 28 31 33 32 36 31 20 |.cxx'\" |\(13261 |
|00003600| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00003610| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 6e 65 77 6d | "s/^X//|" >'newm|
|00003620| 61 74 36 2e 63 78 78 27 | 20 3c 3c 27 45 4e 44 5f |at6.cxx'| <<'END_|
|00003630| 4f 46 5f 46 49 4c 45 27 | 0a 58 2f 2f 24 24 20 6e |OF_FILE'|.X//$$ n|
|00003640| 65 77 6d 61 74 36 2e 63 | 78 78 20 20 20 20 20 20 |ewmat6.c|xx |
|00003650| 20 20 20 20 20 20 4f 70 | 65 72 61 74 6f 72 73 2c | Op|erators,|
|00003660| 20 65 6c 65 6d 65 6e 74 | 20 61 63 63 65 73 73 2c | element| access,|
|00003670| 20 73 75 62 6d 61 74 72 | 69 63 65 73 0a 58 0a 58 | submatr|ices.X.X|
|00003680| 2f 2f 20 43 6f 70 79 72 | 69 67 68 74 20 28 43 29 |// Copyr|ight (C)|
|00003690| 20 31 39 39 31 2c 32 2c | 33 3a 20 52 20 42 20 44 | 1991,2,|3: R B D|
|000036a0| 61 76 69 65 73 0a 58 0a | 58 23 69 6e 63 6c 75 64 |avies.X.|X#includ|
|000036b0| 65 20 22 69 6e 63 6c 75 | 64 65 2e 68 22 0a 58 0a |e "inclu|de.h".X.|
|000036c0| 58 23 69 6e 63 6c 75 64 | 65 20 22 6e 65 77 6d 61 |X#includ|e "newma|
|000036d0| 74 2e 68 22 0a 58 23 69 | 6e 63 6c 75 64 65 20 22 |t.h".X#i|nclude "|
|000036e0| 6e 65 77 6d 61 74 72 63 | 2e 68 22 0a 58 0a 58 0a |newmatrc|.h".X.X.|
|000036f0| 58 2f 2f 23 64 65 66 69 | 6e 65 20 52 45 50 4f 52 |X//#defi|ne REPOR|
|00003700| 54 20 7b 20 73 74 61 74 | 69 63 20 45 78 65 43 6f |T { stat|ic ExeCo|
|00003710| 75 6e 74 65 72 20 45 78 | 65 43 6f 75 6e 74 28 5f |unter Ex|eCount(_|
|00003720| 5f 4c 49 4e 45 5f 5f 2c | 36 29 3b 20 45 78 65 43 |_LINE__,|6); ExeC|
|00003730| 6f 75 6e 74 2b 2b 3b 20 | 7d 0a 58 0a 58 23 64 65 |ount++; |}.X.X#de|
|00003740| 66 69 6e 65 20 52 45 50 | 4f 52 54 20 7b 7d 0a 58 |fine REP|ORT {}.X|
|00003750| 0a 58 2f 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X/*****|********|
|00003760| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 20 67 |********|****** g|
|00003770| 65 6e 65 72 61 6c 20 75 | 74 69 6c 69 74 69 65 73 |eneral u|tilities|
|00003780| 20 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a | *******|********|
|00003790| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2f 0a 58 0a 58 73 |********|**/.X.Xs|
|000037a0| 74 61 74 69 63 20 69 6e | 74 20 74 72 69 73 74 6f |tatic in|t tristo|
|000037b0| 72 65 28 69 6e 74 20 6e | 29 20 20 20 20 20 20 20 |re(int n|) |
|000037c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2f | | /|
|000037d0| 2f 20 65 6c 73 20 69 6e | 20 74 72 69 61 6e 67 75 |/ els in| triangu|
|000037e0| 6c 61 72 20 6d 61 74 72 | 69 78 0a 58 7b 20 72 65 |lar matr|ix.X{ re|
|000037f0| 74 75 72 6e 20 28 6e 2a | 28 6e 2b 31 29 29 2f 32 |turn (n*|(n+1))/2|
|00003800| 3b 20 7d 0a 58 0a 58 0a | 58 2f 2a 2a 2a 2a 2a 2a |; }.X.X.|X/******|
|00003810| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003820| 2a 2a 2a 2a 2a 2a 2a 2a | 20 6f 70 65 72 61 74 6f |********| operato|
|00003830| 72 73 20 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |rs *****|********|
|00003840| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003850| 2a 2a 2f 0a 58 0a 58 52 | 65 61 6c 26 20 4d 61 74 |**/.X.XR|eal& Mat|
|00003860| 72 69 78 3a 3a 6f 70 65 | 72 61 74 6f 72 28 29 28 |rix::ope|rator()(|
|00003870| 69 6e 74 20 6d 2c 20 69 | 6e 74 20 6e 29 0a 58 7b |int m, i|nt n).X{|
|00003880| 0a 58 20 20 20 69 66 20 | 28 6d 3c 3d 30 20 7c 7c |.X if |(m<=0 |||
|00003890| 20 6d 3e 6e 72 6f 77 73 | 20 7c 7c 20 6e 3c 3d 30 | m>nrows| || n<=0|
|000038a0| 20 7c 7c 20 6e 3e 6e 63 | 6f 6c 73 29 0a 58 20 20 | || n>nc|ols).X |
|000038b0| 20 20 20 20 54 68 72 6f | 77 28 49 6e 64 65 78 45 | Thro|w(IndexE|
|000038c0| 78 63 65 70 74 69 6f 6e | 28 6d 2c 6e 2c 2a 74 68 |xception|(m,n,*th|
|000038d0| 69 73 29 29 3b 0a 58 20 | 20 20 72 65 74 75 72 6e |is));.X | return|
|000038e0| 20 73 74 6f 72 65 5b 28 | 6d 2d 31 29 2a 6e 63 6f | store[(|m-1)*nco|
|000038f0| 6c 73 2b 6e 2d 31 5d 3b | 0a 58 7d 0a 58 0a 58 52 |ls+n-1];|.X}.X.XR|
|00003900| 65 61 6c 26 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |eal& Sym|metricMa|
|00003910| 74 72 69 78 3a 3a 6f 70 | 65 72 61 74 6f 72 28 29 |trix::op|erator()|
|00003920| 28 69 6e 74 20 6d 2c 20 | 69 6e 74 20 6e 29 0a 58 |(int m, |int n).X|
|00003930| 7b 0a 58 20 20 20 69 66 | 20 28 6d 3c 3d 30 20 7c |{.X if| (m<=0 ||
|00003940| 7c 20 6e 3c 3d 30 20 7c | 7c 20 6d 3e 6e 72 6f 77 || n<=0 ||| m>nrow|
|00003950| 73 20 7c 7c 20 6e 3e 6e | 63 6f 6c 73 29 0a 58 20 |s || n>n|cols).X |
|00003960| 20 20 20 20 20 54 68 72 | 6f 77 28 49 6e 64 65 78 | Thr|ow(Index|
|00003970| 45 78 63 65 70 74 69 6f | 6e 28 6d 2c 6e 2c 2a 74 |Exceptio|n(m,n,*t|
|00003980| 68 69 73 29 29 3b 0a 58 | 20 20 20 69 66 20 28 6d |his));.X| if (m|
|00003990| 3e 3d 6e 29 20 72 65 74 | 75 72 6e 20 73 74 6f 72 |>=n) ret|urn stor|
|000039a0| 65 5b 74 72 69 73 74 6f | 72 65 28 6d 2d 31 29 2b |e[tristo|re(m-1)+|
|000039b0| 6e 2d 31 5d 3b 0a 58 20 | 20 20 65 6c 73 65 20 72 |n-1];.X | else r|
|000039c0| 65 74 75 72 6e 20 73 74 | 6f 72 65 5b 74 72 69 73 |eturn st|ore[tris|
|000039d0| 74 6f 72 65 28 6e 2d 31 | 29 2b 6d 2d 31 5d 3b 0a |tore(n-1|)+m-1];.|
|000039e0| 58 7d 0a 58 0a 58 52 65 | 61 6c 26 20 55 70 70 65 |X}.X.XRe|al& Uppe|
|000039f0| 72 54 72 69 61 6e 67 75 | 6c 61 72 4d 61 74 72 69 |rTriangu|larMatri|
|00003a00| 78 3a 3a 6f 70 65 72 61 | 74 6f 72 28 29 28 69 6e |x::opera|tor()(in|
|00003a10| 74 20 6d 2c 20 69 6e 74 | 20 6e 29 0a 58 7b 0a 58 |t m, int| n).X{.X|
|00003a20| 20 20 20 69 66 20 28 6d | 3c 3d 30 20 7c 7c 20 6e | if (m|<=0 || n|
|00003a30| 3c 6d 20 7c 7c 20 6e 3e | 6e 63 6f 6c 73 29 0a 58 |<m || n>|ncols).X|
|00003a40| 20 20 20 20 20 20 54 68 | 72 6f 77 28 49 6e 64 65 | Th|row(Inde|
|00003a50| 78 45 78 63 65 70 74 69 | 6f 6e 28 6d 2c 6e 2c 2a |xExcepti|on(m,n,*|
|00003a60| 74 68 69 73 29 29 3b 0a | 58 20 20 20 72 65 74 75 |this));.|X retu|
|00003a70| 72 6e 20 73 74 6f 72 65 | 5b 28 6d 2d 31 29 2a 6e |rn store|[(m-1)*n|
|00003a80| 63 6f 6c 73 2b 6e 2d 31 | 2d 74 72 69 73 74 6f 72 |cols+n-1|-tristor|
|00003a90| 65 28 6d 2d 31 29 5d 3b | 0a 58 7d 0a 58 0a 58 52 |e(m-1)];|.X}.X.XR|
|00003aa0| 65 61 6c 26 20 4c 6f 77 | 65 72 54 72 69 61 6e 67 |eal& Low|erTriang|
|00003ab0| 75 6c 61 72 4d 61 74 72 | 69 78 3a 3a 6f 70 65 72 |ularMatr|ix::oper|
|00003ac0| 61 74 6f 72 28 29 28 69 | 6e 74 20 6d 2c 20 69 6e |ator()(i|nt m, in|
|00003ad0| 74 20 6e 29 0a 58 7b 0a | 58 20 20 20 69 66 20 28 |t n).X{.|X if (|
|00003ae0| 6e 3c 3d 30 20 7c 7c 20 | 6d 3c 6e 20 7c 7c 20 6d |n<=0 || |m<n || m|
|00003af0| 3e 6e 72 6f 77 73 29 0a | 58 20 20 20 20 20 20 54 |>nrows).|X T|
|00003b00| 68 72 6f 77 28 49 6e 64 | 65 78 45 78 63 65 70 74 |hrow(Ind|exExcept|
|00003b10| 69 6f 6e 28 6d 2c 6e 2c | 2a 74 68 69 73 29 29 3b |ion(m,n,|*this));|
|00003b20| 0a 58 20 20 20 72 65 74 | 75 72 6e 20 73 74 6f 72 |.X ret|urn stor|
|00003b30| 65 5b 74 72 69 73 74 6f | 72 65 28 6d 2d 31 29 2b |e[tristo|re(m-1)+|
|00003b40| 6e 2d 31 5d 3b 0a 58 7d | 0a 58 0a 58 52 65 61 6c |n-1];.X}|.X.XReal|
|00003b50| 26 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |& Diagon|alMatrix|
|00003b60| 3a 3a 6f 70 65 72 61 74 | 6f 72 28 29 28 69 6e 74 |::operat|or()(int|
|00003b70| 20 6d 2c 20 69 6e 74 20 | 6e 29 0a 58 7b 0a 58 20 | m, int |n).X{.X |
|00003b80| 20 20 69 66 20 28 6e 3c | 3d 30 20 7c 7c 20 6d 21 | if (n<|=0 || m!|
|00003b90| 3d 6e 20 7c 7c 20 6d 3e | 6e 72 6f 77 73 20 7c 7c |=n || m>|nrows |||
|00003ba0| 20 6e 3e 6e 63 6f 6c 73 | 29 0a 58 20 20 20 20 20 | n>ncols|).X |
|00003bb0| 20 54 68 72 6f 77 28 49 | 6e 64 65 78 45 78 63 65 | Throw(I|ndexExce|
|00003bc0| 70 74 69 6f 6e 28 6d 2c | 6e 2c 2a 74 68 69 73 29 |ption(m,|n,*this)|
|00003bd0| 29 3b 0a 58 20 20 20 72 | 65 74 75 72 6e 20 73 74 |);.X r|eturn st|
|00003be0| 6f 72 65 5b 6e 2d 31 5d | 3b 0a 58 7d 0a 58 0a 58 |ore[n-1]|;.X}.X.X|
|00003bf0| 52 65 61 6c 26 20 44 69 | 61 67 6f 6e 61 6c 4d 61 |Real& Di|agonalMa|
|00003c00| 74 72 69 78 3a 3a 6f 70 | 65 72 61 74 6f 72 28 29 |trix::op|erator()|
|00003c10| 28 69 6e 74 20 6d 29 0a | 58 7b 0a 58 20 20 20 69 |(int m).|X{.X i|
|00003c20| 66 20 28 6d 3c 3d 30 20 | 7c 7c 20 6d 3e 6e 72 6f |f (m<=0 ||| m>nro|
|00003c30| 77 73 29 20 54 68 72 6f | 77 28 49 6e 64 65 78 45 |ws) Thro|w(IndexE|
|00003c40| 78 63 65 70 74 69 6f 6e | 28 6d 2c 2a 74 68 69 73 |xception|(m,*this|
|00003c50| 29 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 73 |));.X |return s|
|00003c60| 74 6f 72 65 5b 6d 2d 31 | 5d 3b 0a 58 7d 0a 58 0a |tore[m-1|];.X}.X.|
|00003c70| 58 52 65 61 6c 26 20 43 | 6f 6c 75 6d 6e 56 65 63 |XReal& C|olumnVec|
|00003c80| 74 6f 72 3a 3a 6f 70 65 | 72 61 74 6f 72 28 29 28 |tor::ope|rator()(|
|00003c90| 69 6e 74 20 6d 29 0a 58 | 7b 0a 58 20 20 20 69 66 |int m).X|{.X if|
|00003ca0| 20 28 6d 3c 3d 30 20 7c | 7c 20 6d 3e 20 6e 72 6f | (m<=0 ||| m> nro|
|00003cb0| 77 73 29 20 54 68 72 6f | 77 28 49 6e 64 65 78 45 |ws) Thro|w(IndexE|
|00003cc0| 78 63 65 70 74 69 6f 6e | 28 6d 2c 2a 74 68 69 73 |xception|(m,*this|
|00003cd0| 29 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 73 |));.X |return s|
|00003ce0| 74 6f 72 65 5b 6d 2d 31 | 5d 3b 0a 58 7d 0a 58 0a |tore[m-1|];.X}.X.|
|00003cf0| 58 52 65 61 6c 26 20 52 | 6f 77 56 65 63 74 6f 72 |XReal& R|owVector|
|00003d00| 3a 3a 6f 70 65 72 61 74 | 6f 72 28 29 28 69 6e 74 |::operat|or()(int|
|00003d10| 20 6e 29 0a 58 7b 0a 58 | 20 20 20 69 66 20 28 6e | n).X{.X| if (n|
|00003d20| 3c 3d 30 20 7c 7c 20 6e | 3e 20 6e 63 6f 6c 73 29 |<=0 || n|> ncols)|
|00003d30| 20 54 68 72 6f 77 28 49 | 6e 64 65 78 45 78 63 65 | Throw(I|ndexExce|
|00003d40| 70 74 69 6f 6e 28 6e 2c | 2a 74 68 69 73 29 29 3b |ption(n,|*this));|
|00003d50| 0a 58 20 20 20 72 65 74 | 75 72 6e 20 73 74 6f 72 |.X ret|urn stor|
|00003d60| 65 5b 6e 2d 31 5d 3b 0a | 58 7d 0a 58 0a 58 52 65 |e[n-1];.|X}.X.XRe|
|00003d70| 61 6c 26 20 42 61 6e 64 | 4d 61 74 72 69 78 3a 3a |al& Band|Matrix::|
|00003d80| 6f 70 65 72 61 74 6f 72 | 28 29 28 69 6e 74 20 6d |operator|()(int m|
|00003d90| 2c 20 69 6e 74 20 6e 29 | 0a 58 7b 0a 58 20 20 20 |, int n)|.X{.X |
|00003da0| 69 6e 74 20 77 20 3d 20 | 75 70 70 65 72 2b 6c 6f |int w = |upper+lo|
|00003db0| 77 65 72 2b 31 3b 20 69 | 6e 74 20 69 20 3d 20 6c |wer+1; i|nt i = l|
|00003dc0| 6f 77 65 72 2b 6e 2d 6d | 3b 0a 58 20 20 20 69 66 |ower+n-m|;.X if|
|00003dd0| 20 28 6d 3c 3d 30 20 7c | 7c 20 6d 3e 6e 72 6f 77 | (m<=0 ||| m>nrow|
|00003de0| 73 20 7c 7c 20 6e 3c 3d | 30 20 7c 7c 20 6e 3e 6e |s || n<=|0 || n>n|
|00003df0| 63 6f 6c 73 20 7c 7c 20 | 69 3c 30 20 7c 7c 20 69 |cols || |i<0 || i|
|00003e00| 3e 3d 77 29 0a 58 20 20 | 20 20 20 20 54 68 72 6f |>=w).X | Thro|
|00003e10| 77 28 49 6e 64 65 78 45 | 78 63 65 70 74 69 6f 6e |w(IndexE|xception|
|00003e20| 28 6d 2c 6e 2c 2a 74 68 | 69 73 29 29 3b 0a 58 20 |(m,n,*th|is));.X |
|00003e30| 20 20 72 65 74 75 72 6e | 20 73 74 6f 72 65 5b 77 | return| store[w|
|00003e40| 2a 28 6d 2d 31 29 2b 69 | 5d 3b 0a 58 7d 0a 58 0a |*(m-1)+i|];.X}.X.|
|00003e50| 58 52 65 61 6c 26 20 55 | 70 70 65 72 42 61 6e 64 |XReal& U|pperBand|
|00003e60| 4d 61 74 72 69 78 3a 3a | 6f 70 65 72 61 74 6f 72 |Matrix::|operator|
|00003e70| 28 29 28 69 6e 74 20 6d | 2c 20 69 6e 74 20 6e 29 |()(int m|, int n)|
|00003e80| 0a 58 7b 0a 58 20 20 20 | 69 6e 74 20 77 20 3d 20 |.X{.X |int w = |
|00003e90| 75 70 70 65 72 2b 31 3b | 20 69 6e 74 20 69 20 3d |upper+1;| int i =|
|00003ea0| 20 6e 2d 6d 3b 0a 58 20 | 20 20 69 66 20 28 6d 3c | n-m;.X | if (m<|
|00003eb0| 3d 30 20 7c 7c 20 6d 3e | 6e 72 6f 77 73 20 7c 7c |=0 || m>|nrows |||
|00003ec0| 20 6e 3c 3d 30 20 7c 7c | 20 6e 3e 6e 63 6f 6c 73 | n<=0 ||| n>ncols|
|00003ed0| 20 7c 7c 20 69 3c 30 20 | 7c 7c 20 69 3e 3d 77 29 | || i<0 ||| i>=w)|
|00003ee0| 0a 58 20 20 20 20 20 20 | 54 68 72 6f 77 28 49 6e |.X |Throw(In|
|00003ef0| 64 65 78 45 78 63 65 70 | 74 69 6f 6e 28 6d 2c 6e |dexExcep|tion(m,n|
|00003f00| 2c 2a 74 68 69 73 29 29 | 3b 0a 58 20 20 20 72 65 |,*this))|;.X re|
|00003f10| 74 75 72 6e 20 73 74 6f | 72 65 5b 77 2a 28 6d 2d |turn sto|re[w*(m-|
|00003f20| 31 29 2b 69 5d 3b 0a 58 | 7d 0a 58 0a 58 52 65 61 |1)+i];.X|}.X.XRea|
|00003f30| 6c 26 20 4c 6f 77 65 72 | 42 61 6e 64 4d 61 74 72 |l& Lower|BandMatr|
|00003f40| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 28 29 28 69 |ix::oper|ator()(i|
|00003f50| 6e 74 20 6d 2c 20 69 6e | 74 20 6e 29 0a 58 7b 0a |nt m, in|t n).X{.|
|00003f60| 58 20 20 20 69 6e 74 20 | 77 20 3d 20 6c 6f 77 65 |X int |w = lowe|
|00003f70| 72 2b 31 3b 20 69 6e 74 | 20 69 20 3d 20 6c 6f 77 |r+1; int| i = low|
|00003f80| 65 72 2b 6e 2d 6d 3b 0a | 58 20 20 20 69 66 20 28 |er+n-m;.|X if (|
|00003f90| 6d 3c 3d 30 20 7c 7c 20 | 6d 3e 6e 72 6f 77 73 20 |m<=0 || |m>nrows |
|00003fa0| 7c 7c 20 6e 3c 3d 30 20 | 7c 7c 20 6e 3e 6e 63 6f ||| n<=0 ||| n>nco|
|00003fb0| 6c 73 20 7c 7c 20 69 3c | 30 20 7c 7c 20 69 3e 3d |ls || i<|0 || i>=|
|00003fc0| 77 29 0a 58 20 20 20 20 | 20 20 54 68 72 6f 77 28 |w).X | Throw(|
|00003fd0| 49 6e 64 65 78 45 78 63 | 65 70 74 69 6f 6e 28 6d |IndexExc|eption(m|
|00003fe0| 2c 6e 2c 2a 74 68 69 73 | 29 29 3b 0a 58 20 20 20 |,n,*this|));.X |
|00003ff0| 72 65 74 75 72 6e 20 73 | 74 6f 72 65 5b 77 2a 28 |return s|tore[w*(|
|00004000| 6d 2d 31 29 2b 69 5d 3b | 0a 58 7d 0a 58 0a 58 52 |m-1)+i];|.X}.X.XR|
|00004010| 65 61 6c 26 20 53 79 6d | 6d 65 74 72 69 63 42 61 |eal& Sym|metricBa|
|00004020| 6e 64 4d 61 74 72 69 78 | 3a 3a 6f 70 65 72 61 74 |ndMatrix|::operat|
|00004030| 6f 72 28 29 28 69 6e 74 | 20 6d 2c 20 69 6e 74 20 |or()(int| m, int |
|00004040| 6e 29 0a 58 7b 0a 58 20 | 20 20 69 6e 74 20 77 20 |n).X{.X | int w |
|00004050| 3d 20 6c 6f 77 65 72 2b | 31 3b 0a 58 20 20 20 69 |= lower+|1;.X i|
|00004060| 66 20 28 6d 3e 3d 6e 29 | 0a 58 20 20 20 7b 0a 58 |f (m>=n)|.X {.X|
|00004070| 20 20 20 20 20 20 69 6e | 74 20 69 20 3d 20 6c 6f | in|t i = lo|
|00004080| 77 65 72 2b 6e 2d 6d 3b | 0a 58 20 20 20 20 20 20 |wer+n-m;|.X |
|00004090| 69 66 20 28 20 6d 3e 6e | 72 6f 77 73 20 7c 7c 20 |if ( m>n|rows || |
|000040a0| 6e 3c 3d 30 20 7c 7c 20 | 69 3c 30 20 29 0a 58 20 |n<=0 || |i<0 ).X |
|000040b0| 20 20 20 20 20 20 20 20 | 54 68 72 6f 77 28 49 6e | |Throw(In|
|000040c0| 64 65 78 45 78 63 65 70 | 74 69 6f 6e 28 6d 2c 6e |dexExcep|tion(m,n|
|000040d0| 2c 2a 74 68 69 73 29 29 | 3b 0a 58 20 20 20 20 20 |,*this))|;.X |
|000040e0| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 77 2a | return |store[w*|
|000040f0| 28 6d 2d 31 29 2b 69 5d | 3b 0a 58 20 20 20 7d 0a |(m-1)+i]|;.X }.|
|00004100| 58 20 20 20 65 6c 73 65 | 0a 58 20 20 20 7b 0a 58 |X else|.X {.X|
|00004110| 20 20 20 20 20 20 69 6e | 74 20 69 20 3d 20 6c 6f | in|t i = lo|
|00004120| 77 65 72 2b 6d 2d 6e 3b | 0a 58 20 20 20 20 20 20 |wer+m-n;|.X |
|00004130| 69 66 20 28 20 6e 3e 6e | 72 6f 77 73 20 7c 7c 20 |if ( n>n|rows || |
|00004140| 6d 3c 3d 30 20 7c 7c 20 | 69 3c 30 20 29 0a 58 20 |m<=0 || |i<0 ).X |
|00004150| 20 20 20 20 20 20 20 20 | 54 68 72 6f 77 28 49 6e | |Throw(In|
|00004160| 64 65 78 45 78 63 65 70 | 74 69 6f 6e 28 6d 2c 6e |dexExcep|tion(m,n|
|00004170| 2c 2a 74 68 69 73 29 29 | 3b 0a 58 20 20 20 20 20 |,*this))|;.X |
|00004180| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 77 2a | return |store[w*|
|00004190| 28 6e 2d 31 29 2b 69 5d | 3b 0a 58 20 20 20 7d 0a |(n-1)+i]|;.X }.|
|000041a0| 58 7d 0a 58 0a 58 23 69 | 66 6e 64 65 66 20 5f 5f |X}.X.X#i|fndef __|
|000041b0| 5a 54 43 5f 5f 0a 58 0a | 58 52 65 61 6c 20 4d 61 |ZTC__.X.|XReal Ma|
|000041c0| 74 72 69 78 3a 3a 6f 70 | 65 72 61 74 6f 72 28 29 |trix::op|erator()|
|000041d0| 28 69 6e 74 20 6d 2c 20 | 69 6e 74 20 6e 29 20 63 |(int m, |int n) c|
|000041e0| 6f 6e 73 74 0a 58 7b 0a | 58 20 20 20 69 66 20 28 |onst.X{.|X if (|
|000041f0| 6d 3c 3d 30 20 7c 7c 20 | 6d 3e 6e 72 6f 77 73 20 |m<=0 || |m>nrows |
|00004200| 7c 7c 20 6e 3c 3d 30 20 | 7c 7c 20 6e 3e 6e 63 6f ||| n<=0 ||| n>nco|
|00004210| 6c 73 29 0a 58 20 20 20 | 20 20 20 54 68 72 6f 77 |ls).X | Throw|
|00004220| 28 49 6e 64 65 78 45 78 | 63 65 70 74 69 6f 6e 28 |(IndexEx|ception(|
|00004230| 6d 2c 6e 2c 2a 74 68 69 | 73 29 29 3b 0a 58 20 20 |m,n,*thi|s));.X |
|00004240| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 28 6d | return |store[(m|
|00004250| 2d 31 29 2a 6e 63 6f 6c | 73 2b 6e 2d 31 5d 3b 0a |-1)*ncol|s+n-1];.|
|00004260| 58 7d 0a 58 0a 58 52 65 | 61 6c 20 53 79 6d 6d 65 |X}.X.XRe|al Symme|
|00004270| 74 72 69 63 4d 61 74 72 | 69 78 3a 3a 6f 70 65 72 |tricMatr|ix::oper|
|00004280| 61 74 6f 72 28 29 28 69 | 6e 74 20 6d 2c 20 69 6e |ator()(i|nt m, in|
|00004290| 74 20 6e 29 20 63 6f 6e | 73 74 0a 58 7b 0a 58 20 |t n) con|st.X{.X |
|000042a0| 20 20 69 66 20 28 6d 3c | 3d 30 20 7c 7c 20 6e 3c | if (m<|=0 || n<|
|000042b0| 3d 30 20 7c 7c 20 6d 3e | 6e 72 6f 77 73 20 7c 7c |=0 || m>|nrows |||
|000042c0| 20 6e 3e 6e 63 6f 6c 73 | 29 0a 58 20 20 20 20 20 | n>ncols|).X |
|000042d0| 20 54 68 72 6f 77 28 49 | 6e 64 65 78 45 78 63 65 | Throw(I|ndexExce|
|000042e0| 70 74 69 6f 6e 28 6d 2c | 6e 2c 2a 74 68 69 73 29 |ption(m,|n,*this)|
|000042f0| 29 3b 0a 58 20 20 20 69 | 66 20 28 6d 3e 3d 6e 29 |);.X i|f (m>=n)|
|00004300| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 74 72 | return |store[tr|
|00004310| 69 73 74 6f 72 65 28 6d | 2d 31 29 2b 6e 2d 31 5d |istore(m|-1)+n-1]|
|00004320| 3b 0a 58 20 20 20 65 6c | 73 65 20 72 65 74 75 72 |;.X el|se retur|
|00004330| 6e 20 73 74 6f 72 65 5b | 74 72 69 73 74 6f 72 65 |n store[|tristore|
|00004340| 28 6e 2d 31 29 2b 6d 2d | 31 5d 3b 0a 58 7d 0a 58 |(n-1)+m-|1];.X}.X|
|00004350| 0a 58 52 65 61 6c 20 55 | 70 70 65 72 54 72 69 61 |.XReal U|pperTria|
|00004360| 6e 67 75 6c 61 72 4d 61 | 74 72 69 78 3a 3a 6f 70 |ngularMa|trix::op|
|00004370| 65 72 61 74 6f 72 28 29 | 28 69 6e 74 20 6d 2c 20 |erator()|(int m, |
|00004380| 69 6e 74 20 6e 29 20 63 | 6f 6e 73 74 0a 58 7b 0a |int n) c|onst.X{.|
|00004390| 58 20 20 20 69 66 20 28 | 6d 3c 3d 30 20 7c 7c 20 |X if (|m<=0 || |
|000043a0| 6e 3c 6d 20 7c 7c 20 6e | 3e 6e 63 6f 6c 73 29 0a |n<m || n|>ncols).|
|000043b0| 58 20 20 20 20 20 20 54 | 68 72 6f 77 28 49 6e 64 |X T|hrow(Ind|
|000043c0| 65 78 45 78 63 65 70 74 | 69 6f 6e 28 6d 2c 6e 2c |exExcept|ion(m,n,|
|000043d0| 2a 74 68 69 73 29 29 3b | 0a 58 20 20 20 72 65 74 |*this));|.X ret|
|000043e0| 75 72 6e 20 73 74 6f 72 | 65 5b 28 6d 2d 31 29 2a |urn stor|e[(m-1)*|
|000043f0| 6e 63 6f 6c 73 2b 6e 2d | 31 2d 74 72 69 73 74 6f |ncols+n-|1-tristo|
|00004400| 72 65 28 6d 2d 31 29 5d | 3b 0a 58 7d 0a 58 0a 58 |re(m-1)]|;.X}.X.X|
|00004410| 52 65 61 6c 20 4c 6f 77 | 65 72 54 72 69 61 6e 67 |Real Low|erTriang|
|00004420| 75 6c 61 72 4d 61 74 72 | 69 78 3a 3a 6f 70 65 72 |ularMatr|ix::oper|
|00004430| 61 74 6f 72 28 29 28 69 | 6e 74 20 6d 2c 20 69 6e |ator()(i|nt m, in|
|00004440| 74 20 6e 29 20 63 6f 6e | 73 74 0a 58 7b 0a 58 20 |t n) con|st.X{.X |
|00004450| 20 20 69 66 20 28 6e 3c | 3d 30 20 7c 7c 20 6d 3c | if (n<|=0 || m<|
|00004460| 6e 20 7c 7c 20 6d 3e 6e | 72 6f 77 73 29 0a 58 20 |n || m>n|rows).X |
|00004470| 20 20 20 20 20 54 68 72 | 6f 77 28 49 6e 64 65 78 | Thr|ow(Index|
|00004480| 45 78 63 65 70 74 69 6f | 6e 28 6d 2c 6e 2c 2a 74 |Exceptio|n(m,n,*t|
|00004490| 68 69 73 29 29 3b 0a 58 | 20 20 20 72 65 74 75 72 |his));.X| retur|
|000044a0| 6e 20 73 74 6f 72 65 5b | 74 72 69 73 74 6f 72 65 |n store[|tristore|
|000044b0| 28 6d 2d 31 29 2b 6e 2d | 31 5d 3b 0a 58 7d 0a 58 |(m-1)+n-|1];.X}.X|
|000044c0| 0a 58 52 65 61 6c 20 44 | 69 61 67 6f 6e 61 6c 4d |.XReal D|iagonalM|
|000044d0| 61 74 72 69 78 3a 3a 6f | 70 65 72 61 74 6f 72 28 |atrix::o|perator(|
|000044e0| 29 28 69 6e 74 20 6d 2c | 20 69 6e 74 20 6e 29 20 |)(int m,| int n) |
|000044f0| 63 6f 6e 73 74 0a 58 7b | 0a 58 20 20 20 69 66 20 |const.X{|.X if |
|00004500| 28 6e 3c 3d 30 20 7c 7c | 20 6d 21 3d 6e 20 7c 7c |(n<=0 ||| m!=n |||
|00004510| 20 6d 3e 6e 72 6f 77 73 | 20 7c 7c 20 6e 3e 6e 63 | m>nrows| || n>nc|
|00004520| 6f 6c 73 29 0a 58 20 20 | 20 20 20 20 54 68 72 6f |ols).X | Thro|
|00004530| 77 28 49 6e 64 65 78 45 | 78 63 65 70 74 69 6f 6e |w(IndexE|xception|
|00004540| 28 6d 2c 6e 2c 2a 74 68 | 69 73 29 29 3b 0a 58 20 |(m,n,*th|is));.X |
|00004550| 20 20 72 65 74 75 72 6e | 20 73 74 6f 72 65 5b 6e | return| store[n|
|00004560| 2d 31 5d 3b 0a 58 7d 0a | 58 0a 58 52 65 61 6c 20 |-1];.X}.|X.XReal |
|00004570| 44 69 61 67 6f 6e 61 6c | 4d 61 74 72 69 78 3a 3a |Diagonal|Matrix::|
|00004580| 6f 70 65 72 61 74 6f 72 | 28 29 28 69 6e 74 20 6d |operator|()(int m|
|00004590| 29 20 63 6f 6e 73 74 0a | 58 7b 0a 58 20 20 20 69 |) const.|X{.X i|
|000045a0| 66 20 28 6d 3c 3d 30 20 | 7c 7c 20 6d 3e 6e 72 6f |f (m<=0 ||| m>nro|
|000045b0| 77 73 29 20 54 68 72 6f | 77 28 49 6e 64 65 78 45 |ws) Thro|w(IndexE|
|000045c0| 78 63 65 70 74 69 6f 6e | 28 6d 2c 2a 74 68 69 73 |xception|(m,*this|
|000045d0| 29 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 73 |));.X |return s|
|000045e0| 74 6f 72 65 5b 6d 2d 31 | 5d 3b 0a 58 7d 0a 58 0a |tore[m-1|];.X}.X.|
|000045f0| 58 52 65 61 6c 20 43 6f | 6c 75 6d 6e 56 65 63 74 |XReal Co|lumnVect|
|00004600| 6f 72 3a 3a 6f 70 65 72 | 61 74 6f 72 28 29 28 69 |or::oper|ator()(i|
|00004610| 6e 74 20 6d 29 20 63 6f | 6e 73 74 0a 58 7b 0a 58 |nt m) co|nst.X{.X|
|00004620| 20 20 20 69 66 20 28 6d | 3c 3d 30 20 7c 7c 20 6d | if (m|<=0 || m|
|00004630| 3e 20 6e 72 6f 77 73 29 | 20 54 68 72 6f 77 28 49 |> nrows)| Throw(I|
|00004640| 6e 64 65 78 45 78 63 65 | 70 74 69 6f 6e 28 6d 2c |ndexExce|ption(m,|
|00004650| 2a 74 68 69 73 29 29 3b | 0a 58 20 20 20 72 65 74 |*this));|.X ret|
|00004660| 75 72 6e 20 73 74 6f 72 | 65 5b 6d 2d 31 5d 3b 0a |urn stor|e[m-1];.|
|00004670| 58 7d 0a 58 0a 58 52 65 | 61 6c 20 52 6f 77 56 65 |X}.X.XRe|al RowVe|
|00004680| 63 74 6f 72 3a 3a 6f 70 | 65 72 61 74 6f 72 28 29 |ctor::op|erator()|
|00004690| 28 69 6e 74 20 6e 29 20 | 63 6f 6e 73 74 0a 58 7b |(int n) |const.X{|
|000046a0| 0a 58 20 20 20 69 66 20 | 28 6e 3c 3d 30 20 7c 7c |.X if |(n<=0 |||
|000046b0| 20 6e 3e 20 6e 63 6f 6c | 73 29 20 54 68 72 6f 77 | n> ncol|s) Throw|
|000046c0| 28 49 6e 64 65 78 45 78 | 63 65 70 74 69 6f 6e 28 |(IndexEx|ception(|
|000046d0| 6e 2c 2a 74 68 69 73 29 | 29 3b 0a 58 20 20 20 72 |n,*this)|);.X r|
|000046e0| 65 74 75 72 6e 20 73 74 | 6f 72 65 5b 6e 2d 31 5d |eturn st|ore[n-1]|
|000046f0| 3b 0a 58 7d 0a 58 0a 58 | 52 65 61 6c 20 42 61 6e |;.X}.X.X|Real Ban|
|00004700| 64 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |dMatrix:|:operato|
|00004710| 72 28 29 28 69 6e 74 20 | 6d 2c 20 69 6e 74 20 6e |r()(int |m, int n|
|00004720| 29 20 63 6f 6e 73 74 0a | 58 7b 0a 58 20 20 20 69 |) const.|X{.X i|
|00004730| 6e 74 20 77 20 3d 20 75 | 70 70 65 72 2b 6c 6f 77 |nt w = u|pper+low|
|00004740| 65 72 2b 31 3b 20 69 6e | 74 20 69 20 3d 20 6c 6f |er+1; in|t i = lo|
|00004750| 77 65 72 2b 6e 2d 6d 3b | 0a 58 20 20 20 69 66 20 |wer+n-m;|.X if |
|00004760| 28 6d 3c 3d 30 20 7c 7c | 20 6d 3e 6e 72 6f 77 73 |(m<=0 ||| m>nrows|
|00004770| 20 7c 7c 20 6e 3c 3d 30 | 20 7c 7c 20 6e 3e 6e 63 | || n<=0| || n>nc|
|00004780| 6f 6c 73 20 7c 7c 20 69 | 3c 30 20 7c 7c 20 69 3e |ols || i|<0 || i>|
|00004790| 3d 77 29 0a 58 20 20 20 | 20 20 20 54 68 72 6f 77 |=w).X | Throw|
|000047a0| 28 49 6e 64 65 78 45 78 | 63 65 70 74 69 6f 6e 28 |(IndexEx|ception(|
|000047b0| 6d 2c 6e 2c 2a 74 68 69 | 73 29 29 3b 0a 58 20 20 |m,n,*thi|s));.X |
|000047c0| 20 72 65 74 75 72 6e 20 | 73 74 6f 72 65 5b 77 2a | return |store[w*|
|000047d0| 28 6d 2d 31 29 2b 69 5d | 3b 0a 58 7d 0a 58 0a 58 |(m-1)+i]|;.X}.X.X|
|000047e0| 52 65 61 6c 20 53 79 6d | 6d 65 74 72 69 63 42 61 |Real Sym|metricBa|
|000047f0| 6e 64 4d 61 74 72 69 78 | 3a 3a 6f 70 65 72 61 74 |ndMatrix|::operat|
|00004800| 6f 72 28 29 28 69 6e 74 | 20 6d 2c 20 69 6e 74 20 |or()(int| m, int |
|00004810| 6e 29 20 63 6f 6e 73 74 | 0a 58 7b 0a 58 20 20 20 |n) const|.X{.X |
|00004820| 69 6e 74 20 77 20 3d 20 | 6c 6f 77 65 72 2b 31 3b |int w = |lower+1;|
|00004830| 0a 58 20 20 20 69 66 20 | 28 6d 3e 3d 6e 29 0a 58 |.X if |(m>=n).X|
|00004840| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 69 6e 74 20 | {.X | int |
|00004850| 69 20 3d 20 6c 6f 77 65 | 72 2b 6e 2d 6d 3b 0a 58 |i = lowe|r+n-m;.X|
|00004860| 20 20 20 20 20 20 69 66 | 20 28 20 6d 3e 6e 72 6f | if| ( m>nro|
|00004870| 77 73 20 7c 7c 20 6e 3c | 3d 30 20 7c 7c 20 69 3c |ws || n<|=0 || i<|
|00004880| 30 20 29 0a 58 20 20 20 | 20 20 20 20 20 20 54 68 |0 ).X | Th|
|00004890| 72 6f 77 28 49 6e 64 65 | 78 45 78 63 65 70 74 69 |row(Inde|xExcepti|
|000048a0| 6f 6e 28 6d 2c 6e 2c 2a | 74 68 69 73 29 29 3b 0a |on(m,n,*|this));.|
|000048b0| 58 20 20 20 20 20 20 72 | 65 74 75 72 6e 20 73 74 |X r|eturn st|
|000048c0| 6f 72 65 5b 77 2a 28 6d | 2d 31 29 2b 69 5d 3b 0a |ore[w*(m|-1)+i];.|
|000048d0| 58 20 20 20 7d 0a 58 20 | 20 20 65 6c 73 65 0a 58 |X }.X | else.X|
|000048e0| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 69 6e 74 20 | {.X | int |
|000048f0| 69 20 3d 20 6c 6f 77 65 | 72 2b 6d 2d 6e 3b 0a 58 |i = lowe|r+m-n;.X|
|00004900| 20 20 20 20 20 20 69 66 | 20 28 20 6e 3e 6e 72 6f | if| ( n>nro|
|00004910| 77 73 20 7c 7c 20 6d 3c | 3d 30 20 7c 7c 20 69 3c |ws || m<|=0 || i<|
|00004920| 30 20 29 0a 58 20 20 20 | 20 20 20 20 20 20 54 68 |0 ).X | Th|
|00004930| 72 6f 77 28 49 6e 64 65 | 78 45 78 63 65 70 74 69 |row(Inde|xExcepti|
|00004940| 6f 6e 28 6d 2c 6e 2c 2a | 74 68 69 73 29 29 3b 0a |on(m,n,*|this));.|
|00004950| 58 20 20 20 20 20 20 72 | 65 74 75 72 6e 20 73 74 |X r|eturn st|
|00004960| 6f 72 65 5b 77 2a 28 6e | 2d 31 29 2b 69 5d 3b 0a |ore[w*(n|-1)+i];.|
|00004970| 58 20 20 20 7d 0a 58 7d | 0a 58 0a 58 23 65 6e 64 |X }.X}|.X.X#end|
|00004980| 69 66 0a 58 0a 58 52 65 | 61 6c 20 42 61 73 65 4d |if.X.XRe|al BaseM|
|00004990| 61 74 72 69 78 3a 3a 41 | 73 53 63 61 6c 61 72 28 |atrix::A|sScalar(|
|000049a0| 29 20 63 6f 6e 73 74 0a | 58 7b 0a 58 20 20 20 52 |) const.|X{.X R|
|000049b0| 45 50 4f 52 54 0a 58 20 | 20 20 47 65 6e 65 72 61 |EPORT.X | Genera|
|000049c0| 6c 4d 61 74 72 69 78 2a | 20 67 6d 20 3d 20 28 28 |lMatrix*| gm = ((|
|000049d0| 42 61 73 65 4d 61 74 72 | 69 78 26 29 2a 74 68 69 |BaseMatr|ix&)*thi|
|000049e0| 73 29 2e 45 76 61 6c 75 | 61 74 65 28 29 3b 0a 58 |s).Evalu|ate();.X|
|000049f0| 20 20 20 69 66 20 28 67 | 6d 2d 3e 6e 72 6f 77 73 | if (g|m->nrows|
|00004a00| 21 3d 31 20 7c 7c 20 67 | 6d 2d 3e 6e 63 6f 6c 73 |!=1 || g|m->ncols|
|00004a10| 21 3d 31 29 0a 58 20 20 | 20 7b 0a 58 20 20 20 20 |!=1).X | {.X |
|00004a20| 20 20 54 72 79 0a 58 20 | 20 20 20 20 20 7b 0a 58 | Try.X | {.X|
|00004a30| 20 20 20 20 20 20 20 20 | 20 54 72 61 63 65 72 20 | | Tracer |
|00004a40| 74 72 28 22 41 73 53 63 | 61 6c 61 72 22 29 3b 0a |tr("AsSc|alar");.|
|00004a50| 58 09 20 54 68 72 6f 77 | 28 50 72 6f 67 72 61 6d |X. Throw|(Program|
|00004a60| 45 78 63 65 70 74 69 6f | 6e 28 22 43 61 6e 6e 6f |Exceptio|n("Canno|
|00004a70| 74 20 63 6f 6e 76 65 72 | 74 20 74 6f 20 73 63 61 |t conver|t to sca|
|00004a80| 6c 61 72 22 2c 20 2a 67 | 6d 29 29 3b 0a 58 20 20 |lar", *g|m));.X |
|00004a90| 20 20 20 20 7d 0a 58 20 | 20 20 20 20 20 43 61 74 | }.X | Cat|
|00004aa0| 63 68 41 6c 6c 20 7b 20 | 67 6d 2d 3e 74 44 65 6c |chAll { |gm->tDel|
|00004ab0| 65 74 65 28 29 3b 20 54 | 68 72 6f 77 28 29 3b 20 |ete(); T|hrow(); |
|00004ac0| 7d 0a 58 20 20 20 7d 0a | 58 20 20 20 52 65 61 6c |}.X }.|X Real|
|00004ad0| 20 78 20 3d 20 2a 28 67 | 6d 2d 3e 73 74 6f 72 65 | x = *(g|m->store|
|00004ae0| 29 3b 20 67 6d 2d 3e 74 | 44 65 6c 65 74 65 28 29 |); gm->t|Delete()|
|00004af0| 3b 20 72 65 74 75 72 6e | 20 78 3b 0a 58 7d 0a 58 |; return| x;.X}.X|
|00004b00| 0a 58 23 69 66 64 65 66 | 20 54 45 4d 50 53 5f 44 |.X#ifdef| TEMPS_D|
|00004b10| 45 53 54 52 4f 59 45 44 | 5f 51 55 49 43 4b 4c 59 |ESTROYED|_QUICKLY|
|00004b20| 0a 58 0a 58 41 64 64 65 | 64 4d 61 74 72 69 78 26 |.X.XAdde|dMatrix&|
|00004b30| 20 42 61 73 65 4d 61 74 | 72 69 78 3a 3a 6f 70 65 | BaseMat|rix::ope|
|00004b40| 72 61 74 6f 72 2b 28 63 | 6f 6e 73 74 20 42 61 73 |rator+(c|onst Bas|
|00004b50| 65 4d 61 74 72 69 78 26 | 20 62 6d 29 20 63 6f 6e |eMatrix&| bm) con|
|00004b60| 73 74 0a 58 7b 0a 58 20 | 20 20 52 45 50 4f 52 54 |st.X{.X | REPORT|
|00004b70| 0a 58 20 20 20 41 64 64 | 65 64 4d 61 74 72 69 78 |.X Add|edMatrix|
|00004b80| 2a 20 78 20 3d 20 6e 65 | 77 20 41 64 64 65 64 4d |* x = ne|w AddedM|
|00004b90| 61 74 72 69 78 28 74 68 | 69 73 2c 20 26 62 6d 29 |atrix(th|is, &bm)|
|00004ba0| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 45 72 72 6f |;.X Ma|trixErro|
|00004bb0| 72 4e 6f 53 70 61 63 65 | 28 78 29 3b 0a 58 20 20 |rNoSpace|(x);.X |
|00004bc0| 20 72 65 74 75 72 6e 20 | 2a 78 3b 0a 58 7d 0a 58 | return |*x;.X}.X|
|00004bd0| 0a 58 4d 75 6c 74 69 70 | 6c 69 65 64 4d 61 74 72 |.XMultip|liedMatr|
|00004be0| 69 78 26 20 42 61 73 65 | 4d 61 74 72 69 78 3a 3a |ix& Base|Matrix::|
|00004bf0| 6f 70 65 72 61 74 6f 72 | 2a 28 63 6f 6e 73 74 20 |operator|*(const |
|00004c00| 42 61 73 65 4d 61 74 72 | 69 78 26 20 62 6d 29 20 |BaseMatr|ix& bm) |
|00004c10| 63 6f 6e 73 74 0a 58 7b | 0a 58 20 20 20 52 45 50 |const.X{|.X REP|
|00004c20| 4f 52 54 0a 58 20 20 20 | 4d 75 6c 74 69 70 6c 69 |ORT.X |Multipli|
|00004c30| 65 64 4d 61 74 72 69 78 | 2a 20 78 20 3d 20 6e 65 |edMatrix|* x = ne|
|00004c40| 77 20 4d 75 6c 74 69 70 | 6c 69 65 64 4d 61 74 72 |w Multip|liedMatr|
|00004c50| 69 78 28 74 68 69 73 2c | 20 26 62 6d 29 3b 0a 58 |ix(this,| &bm);.X|
|00004c60| 20 20 20 4d 61 74 72 69 | 78 45 72 72 6f 72 4e 6f | Matri|xErrorNo|
|00004c70| 53 70 61 63 65 28 78 29 | 3b 0a 58 20 20 20 72 65 |Space(x)|;.X re|
|00004c80| 74 75 72 6e 20 2a 78 3b | 0a 58 7d 0a 58 0a 58 2f |turn *x;|.X}.X.X/|
|00004c90| 2f 53 6f 6c 76 65 64 4d | 61 74 72 69 78 26 20 49 |/SolvedM|atrix& I|
|00004ca0| 6e 76 65 72 74 65 64 4d | 61 74 72 69 78 3a 3a 6f |nvertedM|atrix::o|
|00004cb0| 70 65 72 61 74 6f 72 2a | 28 63 6f 6e 73 74 20 42 |perator*|(const B|
|00004cc0| 61 73 65 4d 61 74 72 69 | 78 26 20 62 6d 78 29 20 |aseMatri|x& bmx) |
|00004cd0| 63 6f 6e 73 74 0a 58 53 | 6f 6c 76 65 64 4d 61 74 |const.XS|olvedMat|
|00004ce0| 72 69 78 26 20 49 6e 76 | 65 72 74 65 64 4d 61 74 |rix& Inv|ertedMat|
|00004cf0| 72 69 78 3a 3a 6f 70 65 | 72 61 74 6f 72 2a 28 63 |rix::ope|rator*(c|
|00004d00| 6f 6e 73 74 20 42 61 73 | 65 4d 61 74 72 69 78 26 |onst Bas|eMatrix&|
|00004d10| 20 62 6d 78 29 0a 58 7b | 0a 58 20 20 20 52 45 50 | bmx).X{|.X REP|
|00004d20| 4f 52 54 0a 58 20 20 20 | 53 6f 6c 76 65 64 4d 61 |ORT.X |SolvedMa|
|00004d30| 74 72 69 78 2a 20 78 3b | 0a 58 20 20 20 54 72 79 |trix* x;|.X Try|
|00004d40| 20 7b 20 78 20 3d 20 6e | 65 77 20 53 6f 6c 76 65 | { x = n|ew Solve|
|00004d50| 64 4d 61 74 72 69 78 28 | 62 6d 2c 20 26 62 6d 78 |dMatrix(|bm, &bmx|
|00004d60| 29 3b 20 4d 61 74 72 69 | 78 45 72 72 6f 72 4e 6f |); Matri|xErrorNo|
|00004d70| 53 70 61 63 65 28 78 29 | 3b 20 7d 0a 58 20 20 20 |Space(x)|; }.X |
|00004d80| 43 61 74 63 68 41 6c 6c | 20 7b 20 64 65 6c 65 74 |CatchAll| { delet|
|00004d90| 65 20 74 68 69 73 3b 20 | 54 68 72 6f 77 28 29 3b |e this; |Throw();|
|00004da0| 20 7d 0a 58 20 20 20 64 | 65 6c 65 74 65 20 74 68 | }.X d|elete th|
|00004db0| 69 73 3b 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |is; | |
|00004dc0| 20 20 20 2f 2f 20 73 69 | 6e 63 65 20 77 65 20 61 | // si|nce we a|
|00004dd0| 72 65 20 75 73 69 6e 67 | 20 62 6d 20 72 61 74 68 |re using| bm rath|
|00004de0| 65 72 20 74 68 61 6e 20 | 74 68 69 73 0a 58 20 20 |er than |this.X |
|00004df0| 20 72 65 74 75 72 6e 20 | 2a 78 3b 0a 58 7d 0a 58 | return |*x;.X}.X|
|00004e00| 0a 58 53 75 62 74 72 61 | 63 74 65 64 4d 61 74 72 |.XSubtra|ctedMatr|
|00004e10| 69 78 26 20 42 61 73 65 | 4d 61 74 72 69 78 3a 3a |ix& Base|Matrix::|
|00004e20| 6f 70 65 72 61 74 6f 72 | 2d 28 63 6f 6e 73 74 20 |operator|-(const |
|00004e30| 42 61 73 65 4d 61 74 72 | 69 78 26 20 62 6d 29 20 |BaseMatr|ix& bm) |
|00004e40| 63 6f 6e 73 74 0a 58 7b | 0a 58 20 20 20 52 45 50 |const.X{|.X REP|
|00004e50| 4f 52 54 0a 58 20 20 20 | 53 75 62 74 72 61 63 74 |ORT.X |Subtract|
|00004e60| 65 64 4d 61 74 72 69 78 | 2a 20 78 20 3d 20 6e 65 |edMatrix|* x = ne|
|00004e70| 77 20 53 75 62 74 72 61 | 63 74 65 64 4d 61 74 72 |w Subtra|ctedMatr|
|00004e80| 69 78 28 74 68 69 73 2c | 20 26 62 6d 29 3b 0a 58 |ix(this,| &bm);.X|
|00004e90| 20 20 20 4d 61 74 72 69 | 78 45 72 72 6f 72 4e 6f | Matri|xErrorNo|
|00004ea0| 53 70 61 63 65 28 78 29 | 3b 0a 58 20 20 20 72 65 |Space(x)|;.X re|
|00004eb0| 74 75 72 6e 20 2a 78 3b | 0a 58 7d 0a 58 0a 58 53 |turn *x;|.X}.X.XS|
|00004ec0| 68 69 66 74 65 64 4d 61 | 74 72 69 78 26 20 42 61 |hiftedMa|trix& Ba|
|00004ed0| 73 65 4d 61 74 72 69 78 | 3a 3a 6f 70 65 72 61 74 |seMatrix|::operat|
|00004ee0| 6f 72 2b 28 52 65 61 6c | 20 66 29 20 63 6f 6e 73 |or+(Real| f) cons|
|00004ef0| 74 0a 58 7b 0a 58 20 20 | 20 52 45 50 4f 52 54 0a |t.X{.X | REPORT.|
|00004f00| 58 20 20 20 53 68 69 66 | 74 65 64 4d 61 74 72 69 |X Shif|tedMatri|
|00004f10| 78 2a 20 78 20 3d 20 6e | 65 77 20 53 68 69 66 74 |x* x = n|ew Shift|
|00004f20| 65 64 4d 61 74 72 69 78 | 28 74 68 69 73 2c 20 66 |edMatrix|(this, f|
|00004f30| 29 3b 0a 58 20 20 20 4d | 61 74 72 69 78 45 72 72 |);.X M|atrixErr|
|00004f40| 6f 72 4e 6f 53 70 61 63 | 65 28 78 29 3b 0a 58 20 |orNoSpac|e(x);.X |
|00004f50| 20 20 72 65 74 75 72 6e | 20 2a 78 3b 0a 58 7d 0a | return| *x;.X}.|
|00004f60| 58 0a 58 53 63 61 6c 65 | 64 4d 61 74 72 69 78 26 |X.XScale|dMatrix&|
|00004f70| 20 42 61 73 65 4d 61 74 | 72 69 78 3a 3a 6f 70 65 | BaseMat|rix::ope|
|00004f80| 72 61 74 6f 72 2a 28 52 | 65 61 6c 20 66 29 20 63 |rator*(R|eal f) c|
|00004f90| 6f 6e 73 74 0a 58 7b 0a | 58 20 20 20 52 45 50 4f |onst.X{.|X REPO|
|00004fa0| 52 54 0a 58 20 20 20 53 | 63 61 6c 65 64 4d 61 74 |RT.X S|caledMat|
|00004fb0| 72 69 78 2a 20 78 20 3d | 20 6e 65 77 20 53 63 61 |rix* x =| new Sca|
|00004fc0| 6c 65 64 4d 61 74 72 69 | 78 28 74 68 69 73 2c 20 |ledMatri|x(this, |
|00004fd0| 66 29 3b 0a 58 20 20 20 | 4d 61 74 72 69 78 45 72 |f);.X |MatrixEr|
|00004fe0| 72 6f 72 4e 6f 53 70 61 | 63 65 28 78 29 3b 0a 58 |rorNoSpa|ce(x);.X|
|00004ff0| 20 20 20 72 65 74 75 72 | 6e 20 2a 78 3b 0a 58 7d | retur|n *x;.X}|
|00005000| 0a 58 0a 58 53 63 61 6c | 65 64 4d 61 74 72 69 78 |.X.XScal|edMatrix|
|00005010| 26 20 42 61 73 65 4d 61 | 74 72 69 78 3a 3a 6f 70 |& BaseMa|trix::op|
|00005020| 65 72 61 74 6f 72 2f 28 | 52 65 61 6c 20 66 29 20 |erator/(|Real f) |
|00005030| 63 6f 6e 73 74 0a 58 7b | 0a 58 20 20 20 52 45 50 |const.X{|.X REP|
|00005040| 4f 52 54 0a 58 20 20 20 | 53 63 61 6c 65 64 4d 61 |ORT.X |ScaledMa|
|00005050| 74 72 69 78 2a 20 78 20 | 3d 20 6e 65 77 20 53 63 |trix* x |= new Sc|
|00005060| 61 6c 65 64 4d 61 74 72 | 69 78 28 74 68 69 73 2c |aledMatr|ix(this,|
|00005070| 20 31 2e 30 2f 66 29 3b | 0a 58 20 20 20 4d 61 74 | 1.0/f);|.X Mat|
|00005080| 72 69 78 45 72 72 6f 72 | 4e 6f 53 70 61 63 65 28 |rixError|NoSpace(|
|00005090| 78 29 3b 0a 58 20 20 20 | 72 65 74 75 72 6e 20 2a |x);.X |return *|
|000050a0| 78 3b 0a 58 7d 0a 58 0a | 58 53 68 69 66 74 65 64 |x;.X}.X.|XShifted|
|000050b0| 4d 61 74 72 69 78 26 20 | 42 61 73 65 4d 61 74 72 |Matrix& |BaseMatr|
|000050c0| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 2d 28 52 65 |ix::oper|ator-(Re|
|000050d0| 61 6c 20 66 29 20 63 6f | 6e 73 74 0a 58 7b 0a 58 |al f) co|nst.X{.X|
|000050e0| 20 20 20 52 45 50 4f 52 | 54 0a 58 20 20 20 53 68 | REPOR|T.X Sh|
|000050f0| 69 66 74 65 64 4d 61 74 | 72 69 78 2a 20 78 20 3d |iftedMat|rix* x =|
|00005100| 20 6e 65 77 20 53 68 69 | 66 74 65 64 4d 61 74 72 | new Shi|ftedMatr|
|00005110| 69 78 28 74 68 69 73 2c | 20 2d 66 29 3b 0a 58 20 |ix(this,| -f);.X |
|00005120| 20 20 4d 61 74 72 69 78 | 45 72 72 6f 72 4e 6f 53 | Matrix|ErrorNoS|
|00005130| 70 61 63 65 28 78 29 3b | 0a 58 20 20 20 72 65 74 |pace(x);|.X ret|
|00005140| 75 72 6e 20 2a 78 3b 0a | 58 7d 0a 58 0a 58 54 72 |urn *x;.|X}.X.XTr|
|00005150| 61 6e 73 70 6f 73 65 64 | 4d 61 74 72 69 78 26 20 |ansposed|Matrix& |
|00005160| 42 61 73 65 4d 61 74 72 | 69 78 3a 3a 74 28 29 20 |BaseMatr|ix::t() |
|00005170| 63 6f 6e 73 74 0a 58 7b | 0a 58 20 20 20 52 45 50 |const.X{|.X REP|
|00005180| 4f 52 54 0a 58 20 20 20 | 54 72 61 6e 73 70 6f 73 |ORT.X |Transpos|
|00005190| 65 64 4d 61 74 72 69 78 | 2a 20 78 20 3d 20 6e 65 |edMatrix|* x = ne|
|000051a0| 77 20 54 72 61 6e 73 70 | 6f 73 65 64 4d 61 74 72 |w Transp|osedMatr|
|000051b0| 69 78 28 74 68 69 73 29 | 3b 0a 58 20 20 20 4d 61 |ix(this)|;.X Ma|
|000051c0| 74 72 69 78 45 72 72 6f | 72 4e 6f 53 70 61 63 65 |trixErro|rNoSpace|
|000051d0| 28 78 29 3b 0a 58 20 20 | 20 72 65 74 75 72 6e 20 |(x);.X | return |
|000051e0| 2a 78 3b 0a 58 7d 0a 58 | 0a 58 4e 65 67 61 74 65 |*x;.X}.X|.XNegate|
|000051f0| 64 4d 61 74 72 69 78 26 | 20 42 61 73 65 4d 61 74 |dMatrix&| BaseMat|
|00005200| 72 69 78 3a 3a 6f 70 65 | 72 61 74 6f 72 2d 28 29 |rix::ope|rator-()|
|00005210| 20 63 6f 6e 73 74 0a 58 | 7b 0a 58 20 20 20 52 45 | const.X|{.X RE|
|00005220| 50 4f 52 54 0a 58 20 20 | 20 4e 65 67 61 74 65 64 |PORT.X | Negated|
|00005230| 4d 61 74 72 69 78 2a 20 | 78 20 3d 20 6e 65 77 20 |Matrix* |x = new |
|00005240| 4e 65 67 61 74 65 64 4d | 61 74 72 69 78 28 74 68 |NegatedM|atrix(th|
|00005250| 69 73 29 3b 0a 58 20 20 | 20 4d 61 74 72 69 78 45 |is);.X | MatrixE|
|00005260| 72 72 6f 72 4e 6f 53 70 | 61 63 65 28 78 29 3b 0a |rrorNoSp|ace(x);.|
|00005270| 58 20 20 20 72 65 74 75 | 72 6e 20 2a 78 3b 0a 58 |X retu|rn *x;.X|
|00005280| 7d 0a 58 0a 58 49 6e 76 | 65 72 74 65 64 4d 61 74 |}.X.XInv|ertedMat|
|00005290| 72 69 78 26 20 42 61 73 | 65 4d 61 74 72 69 78 3a |rix& Bas|eMatrix:|
|000052a0| 3a 69 28 29 20 63 6f 6e | 73 74 0a 58 7b 0a 58 20 |:i() con|st.X{.X |
|000052b0| 20 20 52 45 50 4f 52 54 | 0a 58 20 20 20 49 6e 76 | REPORT|.X Inv|
|000052c0| 65 72 74 65 64 4d 61 74 | 72 69 78 2a 20 78 20 3d |ertedMat|rix* x =|
|000052d0| 20 6e 65 77 20 49 6e 76 | 65 72 74 65 64 4d 61 74 | new Inv|ertedMat|
|000052e0| 72 69 78 28 74 68 69 73 | 29 3b 0a 58 20 20 20 4d |rix(this|);.X M|
|000052f0| 61 74 72 69 78 45 72 72 | 6f 72 4e 6f 53 70 61 63 |atrixErr|orNoSpac|
|00005300| 65 28 78 29 3b 0a 58 20 | 20 20 72 65 74 75 72 6e |e(x);.X | return|
|00005310| 20 2a 78 3b 0a 58 7d 0a | 58 0a 58 43 6f 6e 73 74 | *x;.X}.|X.XConst|
|00005320| 4d 61 74 72 69 78 26 20 | 47 65 6e 65 72 61 6c 4d |Matrix& |GeneralM|
|00005330| 61 74 72 69 78 3a 3a 63 | 28 29 20 63 6f 6e 73 74 |atrix::c|() const|
|00005340| 0a 58 7b 0a 58 20 20 20 | 69 66 20 28 74 61 67 20 |.X{.X |if (tag |
|00005350| 21 3d 20 2d 31 29 0a 58 | 20 20 20 20 20 20 54 68 |!= -1).X| Th|
|00005360| 72 6f 77 28 50 72 6f 67 | 72 61 6d 45 78 63 65 70 |row(Prog|ramExcep|
|00005370| 74 69 6f 6e 28 22 2e 63 | 28 29 20 61 70 70 6c 69 |tion(".c|() appli|
|00005380| 65 64 20 74 6f 20 74 65 | 6d 70 6f 72 61 72 79 20 |ed to te|mporary |
|00005390| 6d 61 74 72 69 78 22 29 | 29 3b 0a 58 20 20 20 52 |matrix")|);.X R|
|000053a0| 45 50 4f 52 54 0a 58 20 | 20 20 43 6f 6e 73 74 4d |EPORT.X | ConstM|
|000053b0| 61 74 72 69 78 2a 20 78 | 20 3d 20 6e 65 77 20 43 |atrix* x| = new C|
|000053c0| 6f 6e 73 74 4d 61 74 72 | 69 78 28 74 68 69 73 29 |onstMatr|ix(this)|
|000053d0| 3b 0a 58 20 20 20 4d 61 | 74 72 69 78 45 72 72 6f |;.X Ma|trixErro|
|000053e0| 72 4e 6f 53 70 61 63 65 | 28 78 29 3b 0a 58 20 20 |rNoSpace|(x);.X |
|000053f0| 20 72 65 74 75 72 6e 20 | 2a 78 3b 0a 58 7d 0a 58 | return |*x;.X}.X|
|00005400| 0a 58 52 6f 77 65 64 4d | 61 74 72 69 78 26 20 42 |.XRowedM|atrix& B|
|00005410| 61 73 65 4d 61 74 72 69 | 78 3a 3a 41 73 52 6f 77 |aseMatri|x::AsRow|
|00005420| 28 29 20 63 6f 6e 73 74 | 0a 58 7b 0a 58 20 20 20 |() const|.X{.X |
|00005430| 52 45 50 4f 52 54 0a 58 | 20 20 20 52 6f 77 65 64 |REPORT.X| Rowed|
|00005440| 4d 61 74 72 69 78 2a 20 | 78 20 3d 20 6e 65 77 20 |Matrix* |x = new |
|00005450| 52 6f 77 65 64 4d 61 74 | 72 69 78 28 74 68 69 73 |RowedMat|rix(this|
|00005460| 29 3b 0a 58 20 20 20 4d | 61 74 72 69 78 45 72 72 |);.X M|atrixErr|
|00005470| 6f 72 4e 6f 53 70 61 63 | 65 28 78 29 3b 0a 58 20 |orNoSpac|e(x);.X |
|00005480| 20 20 72 65 74 75 72 6e | 20 2a 78 3b 0a 58 7d 0a | return| *x;.X}.|
|00005490| 58 0a 58 43 6f 6c 65 64 | 4d 61 74 72 69 78 26 20 |X.XColed|Matrix& |
|000054a0| 42 61 73 65 4d 61 74 72 | 69 78 3a 3a 41 73 43 6f |BaseMatr|ix::AsCo|
|000054b0| 6c 75 6d 6e 28 29 20 63 | 6f 6e 73 74 0a 58 7b 0a |lumn() c|onst.X{.|
|000054c0| 58 20 20 20 52 45 50 4f | 52 54 0a 58 20 20 20 43 |X REPO|RT.X C|
|000054d0| 6f 6c 65 64 4d 61 74 72 | 69 78 2a 20 78 20 3d 20 |oledMatr|ix* x = |
|000054e0| 6e 65 77 20 43 6f 6c 65 | 64 4d 61 74 72 69 78 28 |new Cole|dMatrix(|
|000054f0| 74 68 69 73 29 3b 0a 58 | 20 20 20 4d 61 74 72 69 |this);.X| Matri|
|00005500| 78 45 72 72 6f 72 4e 6f | 53 70 61 63 65 28 78 29 |xErrorNo|Space(x)|
|00005510| 3b 0a 58 20 20 20 72 65 | 74 75 72 6e 20 2a 78 3b |;.X re|turn *x;|
|00005520| 0a 58 7d 0a 58 0a 58 44 | 69 61 67 65 64 4d 61 74 |.X}.X.XD|iagedMat|
|00005530| 72 69 78 26 20 42 61 73 | 65 4d 61 74 72 69 78 3a |rix& Bas|eMatrix:|
|00005540| 3a 41 73 44 69 61 67 6f | 6e 61 6c 28 29 20 63 6f |:AsDiago|nal() co|
|00005550| 6e 73 74 0a 58 7b 0a 58 | 20 20 20 52 45 50 4f 52 |nst.X{.X| REPOR|
|00005560| 54 0a 58 20 20 20 44 69 | 61 67 65 64 4d 61 74 72 |T.X Di|agedMatr|
|00005570| 69 78 2a 20 78 20 3d 20 | 6e 65 77 20 44 69 61 67 |ix* x = |new Diag|
|00005580| 65 64 4d 61 74 72 69 78 | 28 74 68 69 73 29 3b 0a |edMatrix|(this);.|
|00005590| 58 20 20 20 4d 61 74 72 | 69 78 45 72 72 6f 72 4e |X Matr|ixErrorN|
|000055a0| 6f 53 70 61 63 65 28 78 | 29 3b 0a 58 20 20 20 72 |oSpace(x|);.X r|
|000055b0| 65 74 75 72 6e 20 2a 78 | 3b 0a 58 7d 0a 58 0a 58 |eturn *x|;.X}.X.X|
|000055c0| 4d 61 74 65 64 4d 61 74 | 72 69 78 26 20 42 61 73 |MatedMat|rix& Bas|
|000055d0| 65 4d 61 74 72 69 78 3a | 3a 41 73 4d 61 74 72 69 |eMatrix:|:AsMatri|
|000055e0| 78 28 69 6e 74 20 6e 72 | 78 2c 20 69 6e 74 20 6e |x(int nr|x, int n|
|000055f0| 63 78 29 20 63 6f 6e 73 | 74 0a 58 7b 0a 58 20 20 |cx) cons|t.X{.X |
|00005600| 20 52 45 50 4f 52 54 0a | 58 20 20 20 4d 61 74 65 | REPORT.|X Mate|
|00005610| 64 4d 61 74 72 69 78 2a | 20 78 20 3d 20 6e 65 77 |dMatrix*| x = new|
|00005620| 20 4d 61 74 65 64 4d 61 | 74 72 69 78 28 74 68 69 | MatedMa|trix(thi|
|00005630| 73 2c 6e 72 78 2c 6e 63 | 78 29 3b 0a 58 20 20 20 |s,nrx,nc|x);.X |
|00005640| 4d 61 74 72 69 78 45 72 | 72 6f 72 4e 6f 53 70 61 |MatrixEr|rorNoSpa|
|00005650| 63 65 28 78 29 3b 0a 58 | 20 20 20 72 65 74 75 72 |ce(x);.X| retur|
|00005660| 6e 20 2a 78 3b 0a 58 7d | 0a 58 0a 58 23 65 6c 73 |n *x;.X}|.X.X#els|
|00005670| 65 0a 58 0a 58 41 64 64 | 65 64 4d 61 74 72 69 78 |e.X.XAdd|edMatrix|
|00005680| 20 42 61 73 65 4d 61 74 | 72 69 78 3a 3a 6f 70 65 | BaseMat|rix::ope|
|00005690| 72 61 74 6f 72 2b 28 63 | 6f 6e 73 74 20 42 61 73 |rator+(c|onst Bas|
|000056a0| 65 4d 61 74 72 69 78 26 | 20 62 6d 29 20 63 6f 6e |eMatrix&| bm) con|
|000056b0| 73 74 0a 58 7b 20 52 45 | 50 4f 52 54 20 72 65 74 |st.X{ RE|PORT ret|
|000056c0| 75 72 6e 20 41 64 64 65 | 64 4d 61 74 72 69 78 28 |urn Adde|dMatrix(|
|000056d0| 74 68 69 73 2c 20 26 62 | 6d 29 3b 20 7d 0a 58 0a |this, &b|m); }.X.|
|000056e0| 58 4d 75 6c 74 69 70 6c | 69 65 64 4d 61 74 72 69 |XMultipl|iedMatri|
|000056f0| 78 20 42 61 73 65 4d 61 | 74 72 69 78 3a 3a 6f 70 |x BaseMa|trix::op|
|00005700| 65 72 61 74 6f 72 2a 28 | 63 6f 6e 73 74 20 42 61 |erator*(|const Ba|
|00005710| 73 65 4d 61 74 72 69 78 | 26 20 62 6d 29 20 63 6f |seMatrix|& bm) co|
|00005720| 6e 73 74 0a 58 7b 20 52 | 45 50 4f 52 54 20 72 65 |nst.X{ R|EPORT re|
|00005730| 74 75 72 6e 20 4d 75 6c | 74 69 70 6c 69 65 64 4d |turn Mul|tipliedM|
|00005740| 61 74 72 69 78 28 74 68 | 69 73 2c 20 26 62 6d 29 |atrix(th|is, &bm)|
|00005750| 3b 20 7d 0a 58 0a 58 53 | 6f 6c 76 65 64 4d 61 74 |; }.X.XS|olvedMat|
|00005760| 72 69 78 20 49 6e 76 65 | 72 74 65 64 4d 61 74 72 |rix Inve|rtedMatr|
|00005770| 69 78 3a 3a 6f 70 65 72 | 61 74 6f 72 2a 28 63 6f |ix::oper|ator*(co|
|00005780| 6e 73 74 20 42 61 73 65 | 4d 61 74 72 69 78 26 20 |nst Base|Matrix& |
|00005790| 62 6d 78 29 20 63 6f 6e | 73 74 0a 58 7b 20 52 45 |bmx) con|st.X{ RE|
|000057a0| 50 4f 52 54 20 72 65 74 | 75 72 6e 20 53 6f 6c 76 |PORT ret|urn Solv|
|000057b0| 65 64 4d 61 74 72 69 78 | 28 62 6d 2c 20 26 62 6d |edMatrix|(bm, &bm|
|000057c0| 78 29 3b 20 7d 0a 58 0a | 58 53 75 62 74 72 61 63 |x); }.X.|XSubtrac|
|000057d0| 74 65 64 4d 61 74 72 69 | 78 20 42 61 73 65 4d 61 |tedMatri|x BaseMa|
|000057e0| 74 72 69 78 3a 3a 6f 70 | 65 72 61 74 6f 72 2d 28 |trix::op|erator-(|
|000057f0| 63 6f 6e 73 74 20 42 61 | 73 65 4d 61 74 72 69 78 |const Ba|seMatrix|
|00005800| 26 20 62 6d 29 20 63 6f | 6e 73 74 0a 58 7b 20 52 |& bm) co|nst.X{ R|
|00005810| 45 50 4f 52 54 20 72 65 | 74 75 72 6e 20 53 75 62 |EPORT re|turn Sub|
|00005820| 74 72 61 63 74 65 64 4d | 61 74 72 69 78 28 74 68 |tractedM|atrix(th|
|00005830| 69 73 2c 20 26 62 6d 29 | 3b 20 7d 0a 58 0a 58 53 |is, &bm)|; }.X.XS|
|00005840| 68 69 66 74 65 64 4d 61 | 74 72 69 78 20 42 61 73 |hiftedMa|trix Bas|
|00005850| 65 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |eMatrix:|:operato|
|00005860| 72 2b 28 52 65 61 6c 20 | 66 29 20 63 6f 6e 73 74 |r+(Real |f) const|
|00005870| 0a 58 7b 20 52 45 50 4f | 52 54 20 72 65 74 75 72 |.X{ REPO|RT retur|
|00005880| 6e 20 53 68 69 66 74 65 | 64 4d 61 74 72 69 78 28 |n Shifte|dMatrix(|
|00005890| 74 68 69 73 2c 20 66 29 | 3b 20 7d 0a 58 0a 58 53 |this, f)|; }.X.XS|
|000058a0| 63 61 6c 65 64 4d 61 74 | 72 69 78 20 42 61 73 65 |caledMat|rix Base|
|000058b0| 4d 61 74 72 69 78 3a 3a | 6f 70 65 72 61 74 6f 72 |Matrix::|operator|
|000058c0| 2a 28 52 65 61 6c 20 66 | 29 20 63 6f 6e 73 74 0a |*(Real f|) const.|
|000058d0| 58 7b 20 52 45 50 4f 52 | 54 20 72 65 74 75 72 6e |X{ REPOR|T return|
|000058e0| 20 53 63 61 6c 65 64 4d | 61 74 72 69 78 28 74 68 | ScaledM|atrix(th|
|000058f0| 69 73 2c 20 66 29 3b 20 | 7d 0a 58 0a 58 53 63 61 |is, f); |}.X.XSca|
|00005900| 6c 65 64 4d 61 74 72 69 | 78 20 42 61 73 65 4d 61 |ledMatri|x BaseMa|
|00005910| 74 72 69 78 3a 3a 6f 70 | 65 72 61 74 6f 72 2f 28 |trix::op|erator/(|
|00005920| 52 65 61 6c 20 66 29 20 | 63 6f 6e 73 74 0a 58 7b |Real f) |const.X{|
|00005930| 20 52 45 50 4f 52 54 20 | 72 65 74 75 72 6e 20 53 | REPORT |return S|
|00005940| 63 61 6c 65 64 4d 61 74 | 72 69 78 28 74 68 69 73 |caledMat|rix(this|
|00005950| 2c 20 31 2e 30 2f 66 29 | 3b 20 7d 0a 58 0a 58 53 |, 1.0/f)|; }.X.XS|
|00005960| 68 69 66 74 65 64 4d 61 | 74 72 69 78 20 42 61 73 |hiftedMa|trix Bas|
|00005970| 65 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |eMatrix:|:operato|
|00005980| 72 2d 28 52 65 61 6c 20 | 66 29 20 63 6f 6e 73 74 |r-(Real |f) const|
|00005990| 0a 58 7b 20 52 45 50 4f | 52 54 20 72 65 74 75 72 |.X{ REPO|RT retur|
|000059a0| 6e 20 53 68 69 66 74 65 | 64 4d 61 74 72 69 78 28 |n Shifte|dMatrix(|
|000059b0| 74 68 69 73 2c 20 2d 66 | 29 3b 20 7d 0a 58 0a 58 |this, -f|); }.X.X|
|000059c0| 54 72 61 6e 73 70 6f 73 | 65 64 4d 61 74 72 69 78 |Transpos|edMatrix|
|000059d0| 20 42 61 73 65 4d 61 74 | 72 69 78 3a 3a 74 28 29 | BaseMat|rix::t()|
|000059e0| 20 63 6f 6e 73 74 0a 58 | 7b 20 52 45 50 4f 52 54 | const.X|{ REPORT|
|000059f0| 20 72 65 74 75 72 6e 20 | 54 72 61 6e 73 70 6f 73 | return |Transpos|
|00005a00| 65 64 4d 61 74 72 69 78 | 28 74 68 69 73 29 3b 20 |edMatrix|(this); |
|00005a10| 7d 0a 58 0a 58 4e 65 67 | 61 74 65 64 4d 61 74 72 |}.X.XNeg|atedMatr|
|00005a20| 69 78 20 42 61 73 65 4d | 61 74 72 69 78 3a 3a 6f |ix BaseM|atrix::o|
|00005a30| 70 65 72 61 74 6f 72 2d | 28 29 20 63 6f 6e 73 74 |perator-|() const|
|00005a40| 0a 58 7b 20 52 45 50 4f | 52 54 20 72 65 74 75 72 |.X{ REPO|RT retur|
|00005a50| 6e 20 4e 65 67 61 74 65 | 64 4d 61 74 72 69 78 28 |n Negate|dMatrix(|
|00005a60| 74 68 69 73 29 3b 20 7d | 0a 58 0a 58 49 6e 76 65 |this); }|.X.XInve|
|00005a70| 72 74 65 64 4d 61 74 72 | 69 78 20 42 61 73 65 4d |rtedMatr|ix BaseM|
|00005a80| 61 74 72 69 78 3a 3a 69 | 28 29 20 63 6f 6e 73 74 |atrix::i|() const|
|00005a90| 0a 58 7b 20 52 45 50 4f | 52 54 20 72 65 74 75 72 |.X{ REPO|RT retur|
|00005aa0| 6e 20 49 6e 76 65 72 74 | 65 64 4d 61 74 72 69 78 |n Invert|edMatrix|
|00005ab0| 28 74 68 69 73 29 3b 20 | 7d 0a 58 0a 58 43 6f 6e |(this); |}.X.XCon|
|00005ac0| 73 74 4d 61 74 72 69 78 | 20 47 65 6e 65 72 61 6c |stMatrix| General|
|00005ad0| 4d 61 74 72 69 78 3a 3a | 63 28 29 20 63 6f 6e 73 |Matrix::|c() cons|
|00005ae0| 74 0a 58 7b 0a 58 20 20 | 20 69 66 20 28 74 61 67 |t.X{.X | if (tag|
|00005af0| 20 21 3d 20 2d 31 29 0a | 58 20 20 20 20 20 20 54 | != -1).|X T|
|00005b00| 68 72 6f 77 28 50 72 6f | 67 72 61 6d 45 78 63 65 |hrow(Pro|gramExce|
|00005b10| 70 74 69 6f 6e 28 22 2e | 63 28 29 20 61 70 70 6c |ption(".|c() appl|
|00005b20| 69 65 64 20 74 6f 20 74 | 65 6d 70 6f 72 61 72 79 |ied to t|emporary|
|00005b30| 20 6d 61 74 72 69 78 22 | 29 29 3b 0a 58 20 20 20 | matrix"|));.X |
|00005b40| 52 45 50 4f 52 54 20 72 | 65 74 75 72 6e 20 43 6f |REPORT r|eturn Co|
|00005b50| 6e 73 74 4d 61 74 72 69 | 78 28 74 68 69 73 29 3b |nstMatri|x(this);|
|00005b60| 0a 58 7d 0a 58 0a 58 52 | 6f 77 65 64 4d 61 74 72 |.X}.X.XR|owedMatr|
|00005b70| 69 78 20 42 61 73 65 4d | 61 74 72 69 78 3a 3a 41 |ix BaseM|atrix::A|
|00005b80| 73 52 6f 77 28 29 20 63 | 6f 6e 73 74 0a 58 7b 20 |sRow() c|onst.X{ |
|00005b90| 52 45 50 4f 52 54 20 72 | 65 74 75 72 6e 20 52 6f |REPORT r|eturn Ro|
|00005ba0| 77 65 64 4d 61 74 72 69 | 78 28 74 68 69 73 29 3b |wedMatri|x(this);|
|00005bb0| 20 7d 0a 58 0a 58 43 6f | 6c 65 64 4d 61 74 72 69 | }.X.XCo|ledMatri|
|00005bc0| 78 20 42 61 73 65 4d 61 | 74 72 69 78 3a 3a 41 73 |x BaseMa|trix::As|
|00005bd0| 43 6f 6c 75 6d 6e 28 29 | 20 63 6f 6e 73 74 0a 58 |Column()| const.X|
|00005be0| 7b 20 52 45 50 4f 52 54 | 20 72 65 74 75 72 6e 20 |{ REPORT| return |
|00005bf0| 43 6f 6c 65 64 4d 61 74 | 72 69 78 28 74 68 69 73 |ColedMat|rix(this|
|00005c00| 29 3b 20 7d 0a 58 0a 58 | 44 69 61 67 65 64 4d 61 |); }.X.X|DiagedMa|
|00005c10| 74 72 69 78 20 42 61 73 | 65 4d 61 74 72 69 78 3a |trix Bas|eMatrix:|
|00005c20| 3a 41 73 44 69 61 67 6f | 6e 61 6c 28 29 20 63 6f |:AsDiago|nal() co|
|00005c30| 6e 73 74 0a 58 7b 20 52 | 45 50 4f 52 54 20 72 65 |nst.X{ R|EPORT re|
|00005c40| 74 75 72 6e 20 44 69 61 | 67 65 64 4d 61 74 72 69 |turn Dia|gedMatri|
|00005c50| 78 28 74 68 69 73 29 3b | 20 7d 0a 58 0a 58 4d 61 |x(this);| }.X.XMa|
|00005c60| 74 65 64 4d 61 74 72 69 | 78 20 42 61 73 65 4d 61 |tedMatri|x BaseMa|
|00005c70| 74 72 69 78 3a 3a 41 73 | 4d 61 74 72 69 78 28 69 |trix::As|Matrix(i|
|00005c80| 6e 74 20 6e 72 78 2c 20 | 69 6e 74 20 6e 63 78 29 |nt nrx, |int ncx)|
|00005c90| 20 63 6f 6e 73 74 0a 58 | 7b 20 52 45 50 4f 52 54 | const.X|{ REPORT|
|00005ca0| 20 72 65 74 75 72 6e 20 | 4d 61 74 65 64 4d 61 74 | return |MatedMat|
|00005cb0| 72 69 78 28 74 68 69 73 | 2c 6e 72 78 2c 6e 63 78 |rix(this|,nrx,ncx|
|00005cc0| 29 3b 20 7d 0a 58 0a 58 | 23 65 6e 64 69 66 0a 58 |); }.X.X|#endif.X|
|00005cd0| 0a 58 76 6f 69 64 20 47 | 65 6e 65 72 61 6c 4d 61 |.Xvoid G|eneralMa|
|00005ce0| 74 72 69 78 3a 3a 6f 70 | 65 72 61 74 6f 72 3d 28 |trix::op|erator=(|
|00005cf0| 52 65 61 6c 20 66 29 0a | 58 7b 20 52 45 50 4f 52 |Real f).|X{ REPOR|
|00005d00| 54 20 69 6e 74 20 69 3d | 73 74 6f 72 61 67 65 3b |T int i=|storage;|
|00005d10| 20 52 65 61 6c 2a 20 73 | 3d 73 74 6f 72 65 3b 20 | Real* s|=store; |
|00005d20| 77 68 69 6c 65 20 28 69 | 2d 2d 29 20 7b 20 2a 73 |while (i|--) { *s|
|00005d30| 2b 2b 20 3d 20 66 3b 20 | 7d 20 7d 0a 58 0a 58 76 |++ = f; |} }.X.Xv|
|00005d40| 6f 69 64 20 4d 61 74 72 | 69 78 3a 3a 6f 70 65 72 |oid Matr|ix::oper|
|00005d50| 61 74 6f 72 3d 28 63 6f | 6e 73 74 20 42 61 73 65 |ator=(co|nst Base|
|00005d60| 4d 61 74 72 69 78 26 20 | 58 29 0a 58 7b 0a 58 20 |Matrix& |X).X{.X |
|00005d70| 20 20 52 45 50 4f 52 54 | 20 2f 2f 43 68 65 63 6b | REPORT| //Check|
|00005d80| 43 6f 6e 76 65 72 73 69 | 6f 6e 28 58 29 3b 0a 58 |Conversi|on(X);.X|
|00005d90| 20 20 20 4d 61 74 72 69 | 78 43 6f 6e 76 65 72 73 | Matri|xConvers|
|00005da0| 69 6f 6e 43 68 65 63 6b | 20 6d 63 63 3b 0a 58 20 |ionCheck| mcc;.X |
|00005db0| 20 20 45 71 28 58 2c 4d | 61 74 72 69 78 54 79 70 | Eq(X,M|atrixTyp|
|00005dc0| 65 3a 3a 52 74 29 3b 0a | 58 7d 20 0a 58 0a 58 76 |e::Rt);.|X} .X.Xv|
|00005dd0| 6f 69 64 20 52 6f 77 56 | 65 63 74 6f 72 3a 3a 6f |oid RowV|ector::o|
|00005de0| 70 65 72 61 74 6f 72 3d | 28 63 6f 6e 73 74 20 42 |perator=|(const B|
|00005df0| 61 73 65 4d 61 74 72 69 | 78 26 20 58 29 0a 58 7b |aseMatri|x& X).X{|
|00005e00| 0a 58 20 20 20 52 45 50 | 4f 52 54 20 20 2f 2f 20 |.X REP|ORT // |
|00005e10| 43 68 65 63 6b 43 6f 6e | 76 65 72 73 69 6f 6e 28 |CheckCon|version(|
|00005e20| 58 29 3b 0a 58 20 20 20 | 4d 61 74 72 69 78 43 6f |X);.X |MatrixCo|
|00005e30| 6e 76 65 72 73 69 6f 6e | 43 68 65 63 6b 20 6d 63 |nversion|Check mc|
|00005e40| 63 3b 0a 58 20 20 20 45 | 71 28 58 2c 4d 61 74 72 |c;.X E|q(X,Matr|
|00005e50| 69 78 54 79 70 65 3a 3a | 52 56 29 3b 0a 58 20 20 |ixType::|RV);.X |
|00005e60| 20 69 66 20 28 6e 72 6f | 77 73 21 3d 31 29 0a 58 | if (nro|ws!=1).X|
|00005e70| 20 20 20 7b 0a 58 20 20 | 20 20 20 20 54 72 61 63 | {.X | Trac|
|00005e80| 65 72 20 74 72 28 22 52 | 6f 77 56 65 63 74 6f 72 |er tr("R|owVector|
|00005e90| 28 3d 29 22 29 3b 0a 58 | 20 20 20 20 20 20 54 68 |(=)");.X| Th|
|00005ea0| 72 6f 77 28 56 65 63 74 | 6f 72 45 78 63 65 70 74 |row(Vect|orExcept|
|00005eb0| 69 6f 6e 28 2a 74 68 69 | 73 29 29 3b 0a 58 20 20 |ion(*thi|s));.X |
|00005ec0| 20 7d 0a 58 7d 0a 58 0a | 58 76 6f 69 64 20 43 6f | }.X}.X.|Xvoid Co|
|00005ed0| 6c 75 6d 6e 56 65 63 74 | 6f 72 3a 3a 6f 70 65 72 |lumnVect|or::oper|
|00005ee0| 61 74 6f 72 3d 28 63 6f | 6e 73 74 20 42 61 73 65 |ator=(co|nst Base|
|00005ef0| 4d 61 74 72 69 78 26 20 | 58 29 0a 58 7b 0a 58 20 |Matrix& |X).X{.X |
|00005f00| 20 20 52 45 50 4f 52 54 | 20 2f 2f 43 68 65 63 6b | REPORT| //Check|
|00005f10| 43 6f 6e 76 65 72 73 69 | 6f 6e 28 58 29 3b 0a 58 |Conversi|on(X);.X|
|00005f20| 20 20 20 4d 61 74 72 69 | 78 43 6f 6e 76 65 72 73 | Matri|xConvers|
|00005f30| 69 6f 6e 43 68 65 63 6b | 20 6d 63 63 3b 0a 58 20 |ionCheck| mcc;.X |
|00005f40| 20 20 45 71 28 58 2c 4d | 61 74 72 69 78 54 79 70 | Eq(X,M|atrixTyp|
|00005f50| 65 3a 3a 43 56 29 3b 0a | 58 20 20 20 69 66 20 28 |e::CV);.|X if (|
|00005f60| 6e 63 6f 6c 73 21 3d 31 | 29 0a 58 20 20 20 7b 0a |ncols!=1|).X {.|
|00005f70| 58 20 20 20 20 20 20 54 | 72 61 63 65 72 20 74 72 |X T|racer tr|
|00005f80| 28 22 43 6f 6c 75 6d 6e | 56 65 63 74 6f 72 28 3d |("Column|Vector(=|
|00005f90| 29 22 29 3b 0a 58 20 20 | 20 20 20 20 54 68 72 6f |)");.X | Thro|
|00005fa0| 77 28 56 65 63 74 6f 72 | 45 78 63 65 70 74 69 6f |w(Vector|Exceptio|
|00005fb0| 6e 28 2a 74 68 69 73 29 | 29 3b 0a 58 20 20 20 7d |n(*this)|);.X }|
|00005fc0| 0a 58 7d 0a 58 0a 58 76 | 6f 69 64 20 53 79 6d 6d |.X}.X.Xv|oid Symm|
|00005fd0| 65 74 72 69 63 4d 61 74 | 72 69 78 3a 3a 6f 70 65 |etricMat|rix::ope|
|00005fe0| 72 61 74 6f 72 3d 28 63 | 6f 6e 73 74 20 42 61 73 |rator=(c|onst Bas|
|00005ff0| 65 4d 61 74 72 69 78 26 | 20 58 29 0a 58 7b 0a 58 |eMatrix&| X).X{.X|
|00006000| 20 20 20 52 45 50 4f 52 | 54 20 2f 2f 20 43 68 65 | REPOR|T // Che|
|00006010| 63 6b 43 6f 6e 76 65 72 | 73 69 6f 6e 28 58 29 3b |ckConver|sion(X);|
|00006020| 0a 58 20 20 20 4d 61 74 | 72 69 78 43 6f 6e 76 65 |.X Mat|rixConve|
|00006030| 72 73 69 6f 6e 43 68 65 | 63 6b 20 6d 63 63 3b 0a |rsionChe|ck mcc;.|
|00006040| 58 20 20 20 45 71 28 58 | 2c 4d 61 74 72 69 78 54 |X Eq(X|,MatrixT|
|00006050| 79 70 65 3a 3a 53 6d 29 | 3b 0a 58 7d 0a 58 20 0a |ype::Sm)|;.X}.X .|
|00006060| 58 76 6f 69 64 20 55 70 | 70 65 72 54 72 69 61 6e |Xvoid Up|perTrian|
|00006070| 67 75 6c 61 72 4d 61 74 | 72 69 78 3a 3a 6f 70 65 |gularMat|rix::ope|
|00006080| 72 61 74 6f 72 3d 28 63 | 6f 6e 73 74 20 42 61 73 |rator=(c|onst Bas|
|00006090| 65 4d 61 74 72 69 78 26 | 20 58 29 0a 58 7b 0a 58 |eMatrix&| X).X{.X|
|000060a0| 20 20 20 52 45 50 4f 52 | 54 20 2f 2f 43 68 65 63 | REPOR|T //Chec|
|000060b0| 6b 43 6f 6e 76 65 72 73 | 69 6f 6e 28 58 29 3b 0a |kConvers|ion(X);.|
|000060c0| 58 20 20 20 4d 61 74 72 | 69 78 43 6f 6e 76 65 72 |X Matr|ixConver|
|000060d0| 73 69 6f 6e 43 68 65 63 | 6b 20 6d 63 63 3b 0a 58 |sionChec|k mcc;.X|
|000060e0| 20 20 20 45 71 28 58 2c | 4d 61 74 72 69 78 54 79 | Eq(X,|MatrixTy|
|000060f0| 70 65 3a 3a 55 54 29 3b | 0a 58 7d 0a 58 0a 58 76 |pe::UT);|.X}.X.Xv|
|00006100| 6f 69 64 20 4c 6f 77 65 | 72 54 72 69 61 6e 67 75 |oid Lowe|rTriangu|
|00006110| 6c 61 72 4d 61 74 72 69 | 78 3a 3a 6f 70 65 72 61 |larMatri|x::opera|
|00006120| 74 6f 72 3d 28 63 6f 6e | 73 74 20 42 61 73 65 4d |tor=(con|st BaseM|
|00006130| 61 74 72 69 78 26 20 58 | 29 0a 58 7b 0a 58 20 20 |atrix& X|).X{.X |
|00006140| 20 52 45 50 4f 52 54 20 | 2f 2f 43 68 65 63 6b 43 | REPORT |//CheckC|
|00006150| 6f 6e 76 65 72 73 69 6f | 6e 28 58 29 3b 0a 58 20 |onversio|n(X);.X |
|00006160| 20 20 4d 61 74 72 69 78 | 43 6f 6e 76 65 72 73 69 | Matrix|Conversi|
|00006170| 6f 6e 43 68 65 63 6b 20 | 6d 63 63 3b 0a 58 20 20 |onCheck |mcc;.X |
|00006180| 20 45 71 28 58 2c 4d 61 | 74 72 69 78 54 79 70 65 | Eq(X,Ma|trixType|
|00006190| 3a 3a 4c 54 29 3b 0a 58 | 7d 0a 58 0a 58 76 6f 69 |::LT);.X|}.X.Xvoi|
|000061a0| 64 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |d Diagon|alMatrix|
|000061b0| 3a 3a 6f 70 65 72 61 74 | 6f 72 3d 28 63 6f 6e 73 |::operat|or=(cons|
|000061c0| 74 20 42 61 73 65 4d 61 | 74 72 69 78 26 20 58 29 |t BaseMa|trix& X)|
|000061d0| 0a 58 7b 0a 58 20 20 20 | 52 45 50 4f 52 54 20 2f |.X{.X |REPORT /|
|000061e0| 2f 20 43 68 65 63 6b 43 | 6f 6e 76 65 72 73 69 6f |/ CheckC|onversio|
|000061f0| 6e 28 58 29 3b 0a 58 20 | 20 20 4d 61 74 72 69 78 |n(X);.X | Matrix|
|00006200| 43 6f 6e 76 65 72 73 69 | 6f 6e 43 68 65 63 6b 20 |Conversi|onCheck |
|00006210| 6d 63 63 3b 0a 58 20 20 | 20 45 71 28 58 2c 4d 61 |mcc;.X | Eq(X,Ma|
|00006220| 74 72 69 78 54 79 70 65 | 3a 3a 44 67 29 3b 0a 58 |trixType|::Dg);.X|
|00006230| 7d 0a 58 0a 58 76 6f 69 | 64 20 47 65 6e 65 72 61 |}.X.Xvoi|d Genera|
|00006240| 6c 4d 61 74 72 69 78 3a | 3a 6f 70 65 72 61 74 6f |lMatrix:|:operato|
|00006250| 72 3c 3c 28 63 6f 6e 73 | 74 20 52 65 61 6c 2a 20 |r<<(cons|t Real* |
|00006260| 72 29 0a 58 7b 0a 58 20 | 20 20 52 45 50 4f 52 54 |r).X{.X | REPORT|
|00006270| 0a 58 20 20 20 69 6e 74 | 20 69 20 3d 20 73 74 6f |.X int| i = sto|
|00006280| 72 61 67 65 3b 20 52 65 | 61 6c 2a 20 73 3d 73 74 |rage; Re|al* s=st|
|00006290| 6f 72 65 3b 0a 58 20 20 | 20 77 68 69 6c 65 28 69 |ore;.X | while(i|
|000062a0| 2d 2d 29 20 2a 73 2b 2b | 20 3d 20 2a 72 2b 2b 3b |--) *s++| = *r++;|
|000062b0| 0a 58 7d 0a 58 0a 58 0a | 58 2f 2a 2a 2a 2a 2a 2a |.X}.X.X.|X/******|
|000062c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000062d0| 2a 2a 2a 20 65 6c 65 6d | 65 6e 74 20 61 63 63 65 |*** elem|ent acce|
|000062e0| 73 73 20 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |ss *****|********|
|000062f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00006300| 2a 2a 2a 2a 2f 0a 58 0a | 58 52 65 61 6c 26 20 4d |****/.X.|XReal& M|
|00006310| 61 74 72 69 78 3a 3a 65 | 6c 65 6d 65 6e 74 28 69 |atrix::e|lement(i|
|00006320| 6e 74 20 6d 2c 20 69 6e | 74 20 6e 29 0a 58 7b 0a |nt m, in|t n).X{.|
|00006330| 58 20 20 20 69 66 20 28 | 6d 3c 30 20 7c 7c 20 6d |X if (|m<0 || m|
|00006340| 3e 3d 20 6e 72 6f 77 73 | 20 7c 7c 20 6e 3c 30 20 |>= nrows| || n<0 |
|00006350| 7c 7c 20 6e 3e 3d 20 6e | 63 6f 6c 73 29 0a 58 20 ||| n>= n|cols).X |
|00006360| 20 20 20 20 20 54 68 72 | 6f 77 28 49 6e 64 65 78 | Thr|ow(Index|
|00006370| 45 78 63 65 70 74 69 6f | 6e 28 6d 2c 6e 2c 2a 74 |Exceptio|n(m,n,*t|
|00006380| 68 69 73 2c 54 52 55 45 | 29 29 3b 0a 58 20 20 20 |his,TRUE|));.X |
|00006390| 72 65 74 75 72 6e 20 73 | 74 6f 72 65 5b 6d 2a 6e |return s|tore[m*n|
|000063a0| 63 6f 6c 73 2b 6e 5d 3b | 0a 58 7d 0a 58 0a 58 52 |cols+n];|.X}.X.XR|
|000063b0| 65 61 6c 26 20 53 79 6d | 6d 65 74 72 69 63 4d 61 |eal& Sym|metricMa|
|000063c0| 74 72 69 78 3a 3a 65 6c | 65 6d 65 6e 74 28 69 6e |trix::el|ement(in|
|000063d0| 74 20 6d 2c 20 69 6e 74 | 20 6e 29 0a 58 7b 0a 58 |t m, int| n).X{.X|
|000063e0| 20 20 20 69 66 20 28 6d | 3c 30 20 7c 7c 20 6e 3c | if (m|<0 || n<|
|000063f0| 30 20 7c 7c 20 6d 20 3e | 3d 20 6e 72 6f 77 73 20 |0 || m >|= nrows |
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
