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 30 37 3a 20 |bject: v|34i107: |
|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 31 2f 30 38 0a 4d |C++, Par|t01/08.M|
|000000a0| 65 73 73 61 67 65 2d 49 | 44 3a 20 3c 63 73 6d 2d |essage-I|D: <csm-|
|000000b0| 76 33 34 69 31 30 37 3d | 6e 65 77 6d 61 74 30 37 |v34i107=|newmat07|
|000000c0| 2e 30 39 32 37 30 36 40 | 73 70 61 72 6b 79 2e 49 |.092706@|sparky.I|
|000000d0| 4d 44 2e 53 74 65 72 6c | 69 6e 67 2e 43 4f 4d 3e |MD.Sterl|ing.COM>|
|000000e0| 0a 46 6f 6c 6c 6f 77 75 | 70 2d 54 6f 3a 20 63 6f |.Followu|p-To: co|
|000000f0| 6d 70 2e 73 6f 75 72 63 | 65 73 2e 64 0a 58 2d 4d |mp.sourc|es.d.X-M|
|00000100| 64 34 2d 53 69 67 6e 61 | 74 75 72 65 3a 20 31 37 |d4-Signa|ture: 17|
|00000110| 62 65 38 36 31 64 37 38 | 32 31 31 38 37 37 35 33 |be861d78|21187753|
|00000120| 36 64 64 30 34 35 65 34 | 35 36 31 61 32 66 0a 53 |6dd045e4|561a2f.S|
|00000130| 65 6e 64 65 72 3a 20 6b | 65 6e 74 40 73 70 61 72 |ender: k|ent@spar|
|00000140| 6b 79 2e 69 6d 64 2e 73 | 74 65 72 6c 69 6e 67 2e |ky.imd.s|terling.|
|00000150| 63 6f 6d 20 28 4b 65 6e | 74 20 4c 61 6e 64 66 69 |com (Ken|t Landfi|
|00000160| 65 6c 64 29 0a 4f 72 67 | 61 6e 69 7a 61 74 69 6f |eld).Org|anizatio|
|00000170| 6e 3a 20 53 74 65 72 6c | 69 6e 67 20 53 6f 66 74 |n: Sterl|ing Soft|
|00000180| 77 61 72 65 0a 44 61 74 | 65 3a 20 4d 6f 6e 2c 20 |ware.Dat|e: Mon, |
|00000190| 31 31 20 4a 61 6e 20 31 | 39 39 33 20 31 35 3a 32 |11 Jan 1|993 15:2|
|000001a0| 38 3a 32 33 20 47 4d 54 | 0a 41 70 70 72 6f 76 65 |8:23 GMT|.Approve|
|000001b0| 64 3a 20 6b 65 6e 74 40 | 73 70 61 72 6b 79 2e 69 |d: kent@|sparky.i|
|000001c0| 6d 64 2e 73 74 65 72 6c | 69 6e 67 2e 63 6f 6d 0a |md.sterl|ing.com.|
|000001d0| 4c 69 6e 65 73 3a 20 31 | 37 37 37 0a 0a 53 75 62 |Lines: 1|777..Sub|
|000001e0| 6d 69 74 74 65 64 2d 62 | 79 3a 20 72 6f 62 65 72 |mitted-b|y: rober|
|000001f0| 74 64 40 6b 61 75 72 69 | 2e 76 75 77 2e 61 63 2e |td@kauri|.vuw.ac.|
|00000200| 6e 7a 20 28 52 6f 62 65 | 72 74 20 44 61 76 69 65 |nz (Robe|rt Davie|
|00000210| 73 29 0a 50 6f 73 74 69 | 6e 67 2d 6e 75 6d 62 65 |s).Posti|ng-numbe|
|00000220| 72 3a 20 56 6f 6c 75 6d | 65 20 33 34 2c 20 49 73 |r: Volum|e 34, Is|
|00000230| 73 75 65 20 31 30 37 0a | 41 72 63 68 69 76 65 2d |sue 107.|Archive-|
|00000240| 6e 61 6d 65 3a 20 6e 65 | 77 6d 61 74 30 37 2f 70 |name: ne|wmat07/p|
|00000250| 61 72 74 30 31 0a 45 6e | 76 69 72 6f 6e 6d 65 6e |art01.En|vironmen|
|00000260| 74 3a 20 43 2b 2b 0a 53 | 75 70 65 72 73 65 64 65 |t: C++.S|upersede|
|00000270| 73 3a 20 6e 65 77 6d 61 | 74 30 36 3a 20 56 6f 6c |s: newma|t06: Vol|
|00000280| 75 6d 65 20 33 34 2c 20 | 49 73 73 75 65 20 37 2d |ume 34, |Issue 7-|
|00000290| 31 33 0a 0a 20 20 20 4e | 65 77 6d 61 74 20 69 73 |13.. N|ewmat is|
|000002a0| 20 61 6e 20 65 78 70 65 | 72 69 6d 65 6e 74 61 6c | an expe|rimental|
|000002b0| 20 20 6d 61 74 72 69 78 | 20 20 70 61 63 6b 61 67 | matrix| packag|
|000002c0| 65 20 69 6e 20 20 43 2b | 2b 2e 20 20 49 74 20 73 |e in C+|+. It s|
|000002d0| 75 70 70 6f 72 74 73 20 | 20 0a 6d 61 74 72 69 78 |upports | .matrix|
|000002e0| 20 74 79 70 65 73 3a 20 | 4d 61 74 72 69 78 2c 20 | types: |Matrix, |
|000002f0| 55 70 70 65 72 54 72 69 | 61 6e 67 75 6c 61 72 4d |UpperTri|angularM|
|00000300| 61 74 72 69 78 2c 20 4c | 6f 77 65 72 54 72 69 61 |atrix, L|owerTria|
|00000310| 6e 67 75 6c 61 72 4d 61 | 74 72 69 78 2c 20 0a 44 |ngularMa|trix, .D|
|00000320| 69 61 67 6f 6e 61 6c 4d | 61 74 72 69 78 2c 20 53 |iagonalM|atrix, S|
|00000330| 79 6d 6d 65 74 72 69 63 | 4d 61 74 72 69 78 2c 20 |ymmetric|Matrix, |
|00000340| 52 6f 77 56 65 63 74 6f | 72 2c 20 43 6f 6c 75 6d |RowVecto|r, Colum|
|00000350| 6e 56 65 63 74 6f 72 2c | 20 42 61 6e 64 4d 61 74 |nVector,| BandMat|
|00000360| 72 69 78 2c 0a 55 70 70 | 65 72 42 61 6e 64 4d 61 |rix,.Upp|erBandMa|
|00000370| 74 72 69 78 2c 20 4c 6f | 77 65 72 42 61 6e 64 4d |trix, Lo|werBandM|
|00000380| 61 74 72 69 78 20 61 6e | 64 20 53 79 6d 6d 65 74 |atrix an|d Symmet|
|00000390| 72 69 63 42 61 6e 64 4d | 61 74 72 69 78 2e 20 20 |ricBandM|atrix. |
|000003a0| 20 4f 6e 6c 79 20 6f 6e | 65 20 0a 65 6c 65 6d 65 | Only on|e .eleme|
|000003b0| 6e 74 20 74 79 70 65 20 | 28 66 6c 6f 61 74 20 6f |nt type |(float o|
|000003c0| 72 20 64 6f 75 62 6c 65 | 29 20 69 73 20 73 75 70 |r double|) is sup|
|000003d0| 70 6f 72 74 65 64 2e 0a | 0a 20 20 20 54 68 65 20 |ported..|. The |
|000003e0| 70 61 63 6b 61 67 65 20 | 69 6e 63 6c 75 64 65 73 |package |includes|
|000003f0| 20 74 68 65 20 6f 70 65 | 72 61 74 69 6f 6e 73 20 | the ope|rations |
|00000400| 2a 2c 20 2b 2c 20 2d 2c | 20 28 64 65 66 69 6e 65 |*, +, -,| (define|
|00000410| 64 20 61 73 20 6f 70 65 | 72 61 74 6f 72 73 29 0a |d as ope|rators).|
|00000420| 69 6e 76 65 72 73 65 2c | 20 74 72 61 6e 73 70 6f |inverse,| transpo|
|00000430| 73 65 2c 20 63 6f 6e 76 | 65 72 73 69 6f 6e 20 62 |se, conv|ersion b|
|00000440| 65 74 77 65 65 6e 20 74 | 79 70 65 73 2c 20 73 75 |etween t|ypes, su|
|00000450| 62 6d 61 74 72 69 78 2c | 20 64 65 74 65 72 6d 69 |bmatrix,| determi|
|00000460| 6e 61 6e 74 2c 0a 43 68 | 6f 6c 65 73 6b 79 20 64 |nant,.Ch|olesky d|
|00000470| 65 63 6f 6d 70 6f 73 69 | 74 69 6f 6e 2c 20 48 6f |ecomposi|tion, Ho|
|00000480| 75 73 65 68 6f 6c 64 65 | 72 20 74 72 69 61 6e 67 |useholde|r triang|
|00000490| 75 6c 61 72 69 73 61 74 | 69 6f 6e 2c 20 73 69 6e |ularisat|ion, sin|
|000004a0| 67 75 6c 61 72 20 76 61 | 6c 75 65 0a 64 65 63 6f |gular va|lue.deco|
|000004b0| 6d 70 6f 73 69 74 69 6f | 6e 2c 20 65 69 67 65 6e |mpositio|n, eigen|
|000004c0| 76 61 6c 75 65 73 20 6f | 66 20 61 20 73 79 6d 6d |values o|f a symm|
|000004d0| 65 74 72 69 63 20 6d 61 | 74 72 69 78 2c 20 73 6f |etric ma|trix, so|
|000004e0| 72 74 69 6e 67 2c 20 66 | 61 73 74 20 46 6f 75 72 |rting, f|ast Four|
|000004f0| 69 65 72 0a 74 72 61 6e | 73 66 6f 72 6d 2c 20 70 |ier.tran|sform, p|
|00000500| 72 69 6e 74 69 6e 67 2c | 20 61 6e 20 69 6e 74 65 |rinting,| an inte|
|00000510| 72 66 61 63 65 20 77 69 | 74 68 20 22 4e 75 6d 65 |rface wi|th "Nume|
|00000520| 72 69 63 61 6c 20 52 65 | 63 69 70 65 73 20 69 6e |rical Re|cipes in|
|00000530| 20 43 22 20 61 6e 64 20 | 61 6e 0a 65 6d 75 6c 61 | C" and |an.emula|
|00000540| 74 69 6f 6e 20 6f 66 20 | 65 78 63 65 70 74 69 6f |tion of |exceptio|
|00000550| 6e 73 2e 0a 0a 20 20 20 | 49 74 20 69 73 20 69 6e |ns... |It is in|
|00000560| 74 65 6e 64 65 64 20 66 | 6f 72 20 6d 61 74 72 69 |tended f|or matri|
|00000570| 63 65 73 20 69 6e 20 74 | 68 65 20 72 61 6e 67 65 |ces in t|he range|
|00000580| 20 31 35 20 78 20 31 35 | 20 75 70 20 74 6f 20 74 | 15 x 15| up to t|
|00000590| 68 65 20 73 69 7a 65 20 | 79 6f 75 72 0a 63 6f 6d |he size |your.com|
|000005a0| 70 75 74 65 72 20 63 61 | 6e 20 73 74 6f 72 65 20 |puter ca|n store |
|000005b0| 69 6e 20 6f 6e 65 20 62 | 6c 6f 63 6b 2e 0a 2d 2d |in one b|lock..--|
|000005c0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000005d0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000005e0| 2d 2d 2d 2d 2d 2d 2d 2d | 0a 23 21 20 2f 62 69 6e |--------|.#! /bin|
|000005f0| 2f 73 68 0a 23 20 54 68 | 69 73 20 69 73 20 61 20 |/sh.# Th|is is a |
|00000600| 73 68 65 6c 6c 20 61 72 | 63 68 69 76 65 2e 20 20 |shell ar|chive. |
|00000610| 52 65 6d 6f 76 65 20 61 | 6e 79 74 68 69 6e 67 20 |Remove a|nything |
|00000620| 62 65 66 6f 72 65 20 74 | 68 69 73 20 6c 69 6e 65 |before t|his line|
|00000630| 2c 20 74 68 65 6e 20 75 | 6e 70 61 63 6b 0a 23 20 |, then u|npack.# |
|00000640| 69 74 20 62 79 20 73 61 | 76 69 6e 67 20 69 74 20 |it by sa|ving it |
|00000650| 69 6e 74 6f 20 61 20 66 | 69 6c 65 20 61 6e 64 20 |into a f|ile and |
|00000660| 74 79 70 69 6e 67 20 22 | 73 68 20 66 69 6c 65 22 |typing "|sh file"|
|00000670| 2e 20 20 54 6f 20 6f 76 | 65 72 77 72 69 74 65 20 |. To ov|erwrite |
|00000680| 65 78 69 73 74 69 6e 67 | 0a 23 20 66 69 6c 65 73 |existing|.# files|
|00000690| 2c 20 74 79 70 65 20 22 | 73 68 20 66 69 6c 65 20 |, type "|sh file |
|000006a0| 2d 63 22 2e 20 20 59 6f | 75 20 63 61 6e 20 61 6c |-c". Yo|u can al|
|000006b0| 73 6f 20 66 65 65 64 20 | 74 68 69 73 20 61 73 20 |so feed |this as |
|000006c0| 73 74 61 6e 64 61 72 64 | 20 69 6e 70 75 74 20 76 |standard| input v|
|000006d0| 69 61 0a 23 20 75 6e 73 | 68 61 72 2c 20 6f 72 20 |ia.# uns|har, or |
|000006e0| 62 79 20 74 79 70 69 6e | 67 20 22 73 68 20 3c 66 |by typin|g "sh <f|
|000006f0| 69 6c 65 22 2c 20 65 2e | 67 2e 2e 20 20 49 66 20 |ile", e.|g.. If |
|00000700| 74 68 69 73 20 61 72 63 | 68 69 76 65 20 69 73 20 |this arc|hive is |
|00000710| 63 6f 6d 70 6c 65 74 65 | 2c 20 79 6f 75 0a 23 20 |complete|, you.# |
|00000720| 77 69 6c 6c 20 73 65 65 | 20 74 68 65 20 66 6f 6c |will see| the fol|
|00000730| 6c 6f 77 69 6e 67 20 6d | 65 73 73 61 67 65 20 61 |lowing m|essage a|
|00000740| 74 20 74 68 65 20 65 6e | 64 3a 0a 23 09 09 22 45 |t the en|d:.#.."E|
|00000750| 6e 64 20 6f 66 20 61 72 | 63 68 69 76 65 20 31 20 |nd of ar|chive 1 |
|00000760| 28 6f 66 20 38 29 2e 22 | 0a 23 20 43 6f 6e 74 65 |(of 8)."|.# Conte|
|00000770| 6e 74 73 3a 20 20 65 78 | 61 6d 70 6c 65 2e 63 78 |nts: ex|ample.cx|
|00000780| 78 20 6e 65 77 6d 61 74 | 61 2e 74 78 74 20 72 65 |x newmat|a.txt re|
|00000790| 61 64 6d 65 0a 23 20 57 | 72 61 70 70 65 64 20 62 |adme.# W|rapped b|
|000007a0| 79 20 72 6f 62 65 72 74 | 40 6b 65 61 20 6f 6e 20 |y robert|@kea on |
|000007b0| 53 75 6e 20 4a 61 6e 20 | 31 30 20 32 33 3a 35 37 |Sun Jan |10 23:57|
|000007c0| 3a 30 39 20 31 39 39 33 | 0a 50 41 54 48 3d 2f 62 |:09 1993|.PATH=/b|
|000007d0| 69 6e 3a 2f 75 73 72 2f | 62 69 6e 3a 2f 75 73 72 |in:/usr/|bin:/usr|
|000007e0| 2f 75 63 62 20 3b 20 65 | 78 70 6f 72 74 20 50 41 |/ucb ; e|xport PA|
|000007f0| 54 48 0a 69 66 20 74 65 | 73 74 20 2d 66 20 27 65 |TH.if te|st -f 'e|
|00000800| 78 61 6d 70 6c 65 2e 63 | 78 78 27 20 2d 61 20 22 |xample.c|xx' -a "|
|00000810| 24 7b 31 7d 22 20 21 3d | 20 22 2d 63 22 20 3b 20 |${1}" !=| "-c" ; |
|00000820| 74 68 65 6e 20 0a 20 20 | 65 63 68 6f 20 73 68 61 |then . |echo sha|
|00000830| 72 3a 20 57 69 6c 6c 20 | 6e 6f 74 20 63 6c 6f 62 |r: Will |not clob|
|00000840| 62 65 72 20 65 78 69 73 | 74 69 6e 67 20 66 69 6c |ber exis|ting fil|
|00000850| 65 20 5c 22 27 65 78 61 | 6d 70 6c 65 2e 63 78 78 |e \"'exa|mple.cxx|
|00000860| 27 5c 22 0a 65 6c 73 65 | 0a 65 63 68 6f 20 73 68 |'\".else|.echo sh|
|00000870| 61 72 3a 20 45 78 74 72 | 61 63 74 69 6e 67 20 5c |ar: Extr|acting \|
|00000880| 22 27 65 78 61 6d 70 6c | 65 2e 63 78 78 27 5c 22 |"'exampl|e.cxx'\"|
|00000890| 20 5c 28 39 35 34 33 20 | 63 68 61 72 61 63 74 65 | \(9543 |characte|
|000008a0| 72 73 5c 29 0a 73 65 64 | 20 22 73 2f 5e 58 2f 2f |rs\).sed| "s/^X//|
|000008b0| 22 20 3e 27 65 78 61 6d | 70 6c 65 2e 63 78 78 27 |" >'exam|ple.cxx'|
|000008c0| 20 3c 3c 27 45 4e 44 5f | 4f 46 5f 46 49 4c 45 27 | <<'END_|OF_FILE'|
|000008d0| 0a 58 2f 2f 24 24 20 65 | 78 61 6d 70 6c 65 2e 63 |.X//$$ e|xample.c|
|000008e0| 78 78 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |xx | |
|000008f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 45 | | E|
|00000900| 78 61 6d 70 6c 65 20 6f | 66 20 75 73 65 20 6f 66 |xample o|f use of|
|00000910| 20 6d 61 74 72 69 78 20 | 70 61 63 6b 61 67 65 0a | matrix |package.|
|00000920| 58 0a 58 23 64 65 66 69 | 6e 65 20 57 41 4e 54 5f |X.X#defi|ne WANT_|
|00000930| 53 54 52 45 41 4d 20 20 | 20 20 20 20 20 20 20 20 |STREAM | |
|00000940| 20 20 20 20 20 20 20 20 | 2f 2f 20 69 6e 63 6c 75 | |// inclu|
|00000950| 64 65 2e 68 20 77 69 6c | 6c 20 67 65 74 20 73 74 |de.h wil|l get st|
|00000960| 72 65 61 6d 20 66 6e 73 | 0a 58 23 64 65 66 69 6e |ream fns|.X#defin|
|00000970| 65 20 57 41 4e 54 5f 4d | 41 54 48 20 20 20 20 20 |e WANT_M|ATH |
|00000980| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2f | | /|
|00000990| 2f 20 69 6e 63 6c 75 64 | 65 2e 68 20 77 69 6c 6c |/ includ|e.h will|
|000009a0| 20 67 65 74 20 6d 61 74 | 68 20 66 6e 73 0a 58 0a | get mat|h fns.X.|
|000009b0| 58 23 69 6e 63 6c 75 64 | 65 20 22 69 6e 63 6c 75 |X#includ|e "inclu|
|000009c0| 64 65 2e 68 22 20 20 20 | 20 20 20 20 20 20 20 20 |de.h" | |
|000009d0| 20 20 20 20 20 20 2f 2f | 20 69 6e 63 6c 75 64 65 | //| include|
|000009e0| 20 73 74 61 6e 64 61 72 | 64 20 66 69 6c 65 73 0a | standar|d files.|
|000009f0| 58 0a 58 23 69 6e 63 6c | 75 64 65 20 22 6e 65 77 |X.X#incl|ude "new|
|00000a00| 6d 61 74 61 70 2e 68 22 | 20 20 20 20 20 20 20 20 |matap.h"| |
|00000a10| 20 20 20 20 20 20 20 20 | 2f 2f 20 6e 65 65 64 20 | |// need |
|00000a20| 6d 61 74 72 69 78 20 61 | 70 70 6c 69 63 61 74 69 |matrix a|pplicati|
|00000a30| 6f 6e 73 0a 58 0a 58 52 | 65 61 6c 20 74 33 28 52 |ons.X.XR|eal t3(R|
|00000a40| 65 61 6c 29 3b 20 20 20 | 20 20 20 20 20 20 20 20 |eal); | |
|00000a50| 20 20 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 72 | | // r|
|00000a60| 6f 75 6e 64 20 74 6f 20 | 33 20 64 65 63 69 6d 61 |ound to |3 decima|
|00000a70| 6c 20 70 6c 61 63 65 73 | 0a 58 0a 58 0a 58 2f 2f |l places|.X.X.X//|
|00000a80| 20 64 65 6d 6f 6e 73 74 | 72 61 74 69 6f 6e 20 6f | demonst|ration o|
|00000a90| 66 20 6d 61 74 72 69 78 | 20 70 61 63 6b 61 67 65 |f matrix| package|
|00000aa0| 20 6f 6e 20 6c 69 6e 65 | 61 72 20 72 65 67 72 65 | on line|ar regre|
|00000ab0| 73 73 69 6f 6e 20 70 72 | 6f 62 6c 65 6d 0a 58 0a |ssion pr|oblem.X.|
|00000ac0| 58 0a 58 76 6f 69 64 20 | 74 65 73 74 31 28 52 65 |X.Xvoid |test1(Re|
|00000ad0| 61 6c 2a 20 79 2c 20 52 | 65 61 6c 2a 20 78 31 2c |al* y, R|eal* x1,|
|00000ae0| 20 52 65 61 6c 2a 20 78 | 32 2c 20 69 6e 74 20 6e | Real* x|2, int n|
|00000af0| 6f 62 73 2c 20 69 6e 74 | 20 6e 70 72 65 64 29 0a |obs, int| npred).|
|00000b00| 58 7b 0a 58 20 20 20 63 | 6f 75 74 20 3c 3c 20 22 |X{.X c|out << "|
|00000b10| 5c 6e 5c 6e 54 65 73 74 | 20 31 20 2d 20 74 72 61 |\n\nTest| 1 - tra|
|00000b20| 64 69 74 69 6f 6e 61 6c | 5c 6e 22 3b 0a 58 0a 58 |ditional|\n";.X.X|
|00000b30| 20 20 20 2f 2f 20 74 72 | 61 64 69 74 69 6f 6e 61 | // tr|aditiona|
|00000b40| 6c 20 73 75 6d 20 6f 66 | 20 73 71 75 61 72 65 73 |l sum of| squares|
|00000b50| 20 61 6e 64 20 70 72 6f | 64 75 63 74 73 20 6d 65 | and pro|ducts me|
|00000b60| 74 68 6f 64 20 6f 66 20 | 63 61 6c 63 75 6c 61 74 |thod of |calculat|
|00000b70| 69 6f 6e 0a 58 20 20 20 | 2f 2f 20 77 69 74 68 20 |ion.X |// with |
|00000b80| 73 75 62 74 72 61 63 74 | 69 6f 6e 20 6f 66 20 6d |subtract|ion of m|
|00000b90| 65 61 6e 73 0a 58 0a 58 | 20 20 20 2f 2f 20 6d 61 |eans.X.X| // ma|
|00000ba0| 6b 65 20 6d 61 74 72 69 | 78 20 6f 66 20 70 72 65 |ke matri|x of pre|
|00000bb0| 64 69 63 74 6f 72 20 76 | 61 6c 75 65 73 0a 58 20 |dictor v|alues.X |
|00000bc0| 20 20 4d 61 74 72 69 78 | 20 58 28 6e 6f 62 73 2c | Matrix| X(nobs,|
|00000bd0| 6e 70 72 65 64 29 3b 0a | 58 0a 58 20 20 20 2f 2f |npred);.|X.X //|
|00000be0| 20 6c 6f 61 64 20 78 31 | 20 61 6e 64 20 78 32 20 | load x1| and x2 |
|00000bf0| 69 6e 74 6f 20 58 0a 58 | 20 20 20 2f 2f 20 20 20 |into X.X| // |
|00000c00| 20 5b 75 73 65 20 3c 3c | 20 72 61 74 68 65 72 20 | [use <<| rather |
|00000c10| 74 68 61 6e 20 3d 20 77 | 69 74 68 20 73 75 62 6d |than = w|ith subm|
|00000c20| 61 74 72 69 63 65 73 20 | 61 6e 64 2f 6f 72 20 6c |atrices |and/or l|
|00000c30| 6f 61 64 69 6e 67 20 61 | 72 72 61 79 73 5d 0a 58 |oading a|rrays].X|
|00000c40| 20 20 20 58 2e 43 6f 6c | 75 6d 6e 28 31 29 20 3c | X.Col|umn(1) <|
|00000c50| 3c 20 78 31 3b 20 20 58 | 2e 43 6f 6c 75 6d 6e 28 |< x1; X|.Column(|
|00000c60| 32 29 20 3c 3c 20 78 32 | 3b 0a 58 0a 58 20 20 20 |2) << x2|;.X.X |
|00000c70| 2f 2f 20 76 65 63 74 6f | 72 20 6f 66 20 59 20 76 |// vecto|r of Y v|
|00000c80| 61 6c 75 65 73 0a 58 20 | 20 20 43 6f 6c 75 6d 6e |alues.X | Column|
|00000c90| 56 65 63 74 6f 72 20 59 | 28 6e 6f 62 73 29 3b 20 |Vector Y|(nobs); |
|00000ca0| 59 20 3c 3c 20 79 3b 0a | 58 0a 58 20 20 20 2f 2f |Y << y;.|X.X //|
|00000cb0| 20 6d 61 6b 65 20 76 65 | 63 74 6f 72 20 6f 66 20 | make ve|ctor of |
|00000cc0| 31 73 0a 58 20 20 20 43 | 6f 6c 75 6d 6e 56 65 63 |1s.X C|olumnVec|
|00000cd0| 74 6f 72 20 4f 6e 65 73 | 28 6e 6f 62 73 29 3b 20 |tor Ones|(nobs); |
|00000ce0| 4f 6e 65 73 20 3d 20 31 | 2e 30 3b 0a 58 0a 58 20 |Ones = 1|.0;.X.X |
|00000cf0| 20 20 2f 2f 20 63 61 6c | 63 75 6c 61 74 65 20 6d | // cal|culate m|
|00000d00| 65 61 6e 73 20 28 61 76 | 65 72 61 67 65 73 29 20 |eans (av|erages) |
|00000d10| 6f 66 20 78 31 20 61 6e | 64 20 78 32 20 5b 20 2e |of x1 an|d x2 [ .|
|00000d20| 74 28 29 20 74 61 6b 65 | 73 20 74 72 61 6e 73 70 |t() take|s transp|
|00000d30| 6f 73 65 5d 0a 58 20 20 | 20 52 6f 77 56 65 63 74 |ose].X | RowVect|
|00000d40| 6f 72 20 4d 20 3d 20 4f | 6e 65 73 2e 74 28 29 20 |or M = O|nes.t() |
|00000d50| 2a 20 58 20 2f 20 6e 6f | 62 73 3b 0a 58 0a 58 20 |* X / no|bs;.X.X |
|00000d60| 20 20 2f 2f 20 61 6e 64 | 20 73 75 62 74 72 61 63 | // and| subtrac|
|00000d70| 74 20 6d 65 61 6e 73 20 | 66 72 6f 6d 20 78 31 20 |t means |from x1 |
|00000d80| 61 6e 64 20 78 31 0a 58 | 20 20 20 4d 61 74 72 69 |and x1.X| Matri|
|00000d90| 78 20 58 43 28 6e 6f 62 | 73 2c 6e 70 72 65 64 29 |x XC(nob|s,npred)|
|00000da0| 3b 0a 58 20 20 20 58 43 | 20 3d 20 58 20 2d 20 4f |;.X XC| = X - O|
|00000db0| 6e 65 73 20 2a 20 4d 3b | 0a 58 0a 58 20 20 20 2f |nes * M;|.X.X /|
|00000dc0| 2f 20 64 6f 20 74 68 65 | 20 73 61 6d 65 20 74 6f |/ do the| same to|
|00000dd0| 20 59 20 5b 6e 65 65 64 | 20 22 52 65 61 6c 22 20 | Y [need| "Real" |
|00000de0| 74 6f 20 63 6f 6e 76 65 | 72 74 20 31 78 31 20 6d |to conve|rt 1x1 m|
|00000df0| 61 74 72 69 78 20 74 6f | 20 73 63 61 6c 61 72 5d |atrix to| scalar]|
|00000e00| 0a 58 20 20 20 43 6f 6c | 75 6d 6e 56 65 63 74 6f |.X Col|umnVecto|
|00000e10| 72 20 59 43 28 6e 6f 62 | 73 29 3b 0a 58 20 20 20 |r YC(nob|s);.X |
|00000e20| 52 65 61 6c 20 6d 20 3d | 20 28 4f 6e 65 73 2e 74 |Real m =| (Ones.t|
|00000e30| 28 29 20 2a 20 59 29 2e | 41 73 53 63 61 6c 61 72 |() * Y).|AsScalar|
|00000e40| 28 29 20 2f 20 6e 6f 62 | 73 3b 20 20 59 43 20 3d |() / nob|s; YC =|
|00000e50| 20 59 20 2d 20 4f 6e 65 | 73 20 2a 20 6d 3b 0a 58 | Y - One|s * m;.X|
|00000e60| 0a 58 20 20 20 2f 2f 20 | 66 6f 72 6d 20 73 75 6d |.X // |form sum|
|00000e70| 20 6f 66 20 73 71 75 61 | 72 65 73 20 61 6e 64 20 | of squa|res and |
|00000e80| 70 72 6f 64 75 63 74 20 | 6d 61 74 72 69 78 0a 58 |product |matrix.X|
|00000e90| 20 20 20 2f 2f 20 20 20 | 20 5b 75 73 65 20 3c 3c | // | [use <<|
|00000ea0| 20 72 61 74 68 65 72 20 | 74 68 61 6e 20 3d 20 66 | rather |than = f|
|00000eb0| 6f 72 20 63 6f 70 79 69 | 6e 67 20 4d 61 74 72 69 |or copyi|ng Matri|
|00000ec0| 78 20 69 6e 74 6f 20 53 | 79 6d 6d 65 74 72 69 63 |x into S|ymmetric|
|00000ed0| 4d 61 74 72 69 78 5d 0a | 58 20 20 20 53 79 6d 6d |Matrix].|X Symm|
|00000ee0| 65 74 72 69 63 4d 61 74 | 72 69 78 20 53 53 51 3b |etricMat|rix SSQ;|
|00000ef0| 20 53 53 51 20 3c 3c 20 | 58 43 2e 74 28 29 20 2a | SSQ << |XC.t() *|
|00000f00| 20 58 43 3b 0a 58 0a 58 | 20 20 20 2f 2f 20 63 61 | XC;.X.X| // ca|
|00000f10| 6c 63 75 6c 61 74 65 20 | 65 73 74 69 6d 61 74 65 |lculate |estimate|
|00000f20| 0a 58 20 20 20 2f 2f 20 | 20 20 20 5b 62 72 61 63 |.X // | [brac|
|00000f30| 6b 65 74 20 6c 61 73 74 | 20 74 77 6f 20 74 65 72 |ket last| two ter|
|00000f40| 6d 73 20 74 6f 20 66 6f | 72 63 65 20 74 68 69 73 |ms to fo|rce this|
|00000f50| 20 6d 75 6c 74 69 70 6c | 69 63 61 74 69 6f 6e 20 | multipl|ication |
|00000f60| 66 69 72 73 74 5d 0a 58 | 20 20 20 2f 2f 20 20 20 |first].X| // |
|00000f70| 20 5b 20 2e 69 28 29 20 | 6d 65 61 6e 73 20 69 6e | [ .i() |means in|
|00000f80| 76 65 72 73 65 2c 20 62 | 75 74 20 69 6e 76 65 72 |verse, b|ut inver|
|00000f90| 73 65 20 69 73 20 6e 6f | 74 20 65 78 70 6c 69 63 |se is no|t explic|
|00000fa0| 69 74 79 20 63 61 6c 63 | 75 6c 61 74 65 64 5d 0a |ity calc|ulated].|
|00000fb0| 58 20 20 20 43 6f 6c 75 | 6d 6e 56 65 63 74 6f 72 |X Colu|mnVector|
|00000fc0| 20 41 20 3d 20 53 53 51 | 2e 69 28 29 20 2a 20 28 | A = SSQ|.i() * (|
|00000fd0| 58 43 2e 74 28 29 20 2a | 20 59 43 29 3b 0a 58 0a |XC.t() *| YC);.X.|
|00000fe0| 58 20 20 20 2f 2f 20 63 | 61 6c 63 75 6c 61 74 65 |X // c|alculate|
|00000ff0| 20 65 73 74 69 6d 61 74 | 65 20 6f 66 20 63 6f 6e | estimat|e of con|
|00001000| 73 74 61 6e 74 20 74 65 | 72 6d 0a 58 20 20 20 52 |stant te|rm.X R|
|00001010| 65 61 6c 20 61 20 3d 20 | 6d 20 2d 20 28 4d 20 2a |eal a = |m - (M *|
|00001020| 20 41 29 2e 41 73 53 63 | 61 6c 61 72 28 29 3b 0a | A).AsSc|alar();.|
|00001030| 58 0a 58 20 20 20 2f 2f | 20 47 65 74 20 76 61 72 |X.X //| Get var|
|00001040| 69 61 6e 63 65 73 20 6f | 66 20 65 73 74 69 6d 61 |iances o|f estima|
|00001050| 74 65 73 20 66 72 6f 6d | 20 64 69 61 67 6f 6e 61 |tes from| diagona|
|00001060| 6c 20 65 6c 65 6d 65 6e | 74 73 20 6f 66 20 69 6e |l elemen|ts of in|
|00001070| 76 6f 69 63 65 20 6f 66 | 20 53 53 51 0a 58 20 20 |voice of| SSQ.X |
|00001080| 20 2f 2f 20 20 20 20 5b | 20 77 65 20 61 72 65 20 | // [| we are |
|00001090| 74 61 6b 69 6e 67 20 69 | 6e 76 65 72 73 65 20 6f |taking i|nverse o|
|000010a0| 66 20 53 53 51 3b 20 77 | 6f 75 6c 64 20 68 61 76 |f SSQ; w|ould hav|
|000010b0| 65 20 62 65 65 6e 20 62 | 65 74 74 65 72 20 74 6f |e been b|etter to|
|000010c0| 20 75 73 65 0a 58 20 20 | 20 2f 2f 20 20 20 20 20 | use.X | // |
|000010d0| 20 20 20 43 72 6f 75 74 | 4d 61 74 72 69 78 20 6d | Crout|Matrix m|
|000010e0| 65 74 68 6f 64 20 2d 20 | 73 65 65 20 64 6f 63 75 |ethod - |see docu|
|000010f0| 6d 65 6e 74 61 74 69 6f | 6e 20 5d 0a 58 20 20 20 |mentatio|n ].X |
|00001100| 4d 61 74 72 69 78 20 49 | 53 53 51 20 3d 20 53 53 |Matrix I|SSQ = SS|
|00001110| 51 2e 69 28 29 3b 20 44 | 69 61 67 6f 6e 61 6c 4d |Q.i(); D|iagonalM|
|00001120| 61 74 72 69 78 20 44 3b | 20 44 20 3c 3c 20 49 53 |atrix D;| D << IS|
|00001130| 53 51 3b 0a 58 20 20 20 | 43 6f 6c 75 6d 6e 56 65 |SQ;.X |ColumnVe|
|00001140| 63 74 6f 72 20 56 20 3d | 20 44 2e 41 73 43 6f 6c |ctor V =| D.AsCol|
|00001150| 75 6d 6e 28 29 3b 0a 58 | 20 20 20 52 65 61 6c 20 |umn();.X| Real |
|00001160| 76 20 3d 20 31 2e 30 2f | 6e 6f 62 73 20 2b 20 28 |v = 1.0/|nobs + (|
|00001170| 4d 20 2a 20 49 53 53 51 | 20 2a 20 4d 2e 74 28 29 |M * ISSQ| * M.t()|
|00001180| 29 2e 41 73 53 63 61 6c | 61 72 28 29 3b 0a 58 09 |).AsScal|ar();.X.|
|00001190| 09 09 09 09 20 20 20 20 | 2f 2f 20 66 6f 72 20 63 |.... |// for c|
|000011a0| 61 6c 63 20 76 61 72 69 | 61 6e 63 65 20 63 6f 6e |alc vari|ance con|
|000011b0| 73 74 0a 58 0a 58 20 20 | 20 2f 2f 20 43 61 6c 63 |st.X.X | // Calc|
|000011c0| 75 6c 61 74 65 20 66 69 | 74 74 65 64 20 76 61 6c |ulate fi|tted val|
|000011d0| 75 65 73 20 61 6e 64 20 | 72 65 73 69 64 75 61 6c |ues and |residual|
|000011e0| 73 0a 58 20 20 20 69 6e | 74 20 6e 70 72 65 64 31 |s.X in|t npred1|
|000011f0| 20 3d 20 6e 70 72 65 64 | 2b 31 3b 0a 58 20 20 20 | = npred|+1;.X |
|00001200| 43 6f 6c 75 6d 6e 56 65 | 63 74 6f 72 20 46 69 74 |ColumnVe|ctor Fit|
|00001210| 74 65 64 20 3d 20 58 20 | 2a 20 41 20 2b 20 61 3b |ted = X |* A + a;|
|00001220| 0a 58 20 20 20 43 6f 6c | 75 6d 6e 56 65 63 74 6f |.X Col|umnVecto|
|00001230| 72 20 52 65 73 69 64 75 | 61 6c 20 3d 20 59 20 2d |r Residu|al = Y -|
|00001240| 20 46 69 74 74 65 64 3b | 0a 58 20 20 20 52 65 61 | Fitted;|.X Rea|
|00001250| 6c 20 52 65 73 56 61 72 | 20 3d 20 52 65 73 69 64 |l ResVar| = Resid|
|00001260| 75 61 6c 2e 53 75 6d 53 | 71 75 61 72 65 28 29 20 |ual.SumS|quare() |
|00001270| 2f 20 28 6e 6f 62 73 2d | 6e 70 72 65 64 31 29 3b |/ (nobs-|npred1);|
|00001280| 0a 58 0a 58 20 20 20 2f | 2f 20 47 65 74 20 64 69 |.X.X /|/ Get di|
|00001290| 61 67 6f 6e 61 6c 73 20 | 6f 66 20 48 61 74 20 6d |agonals |of Hat m|
|000012a0| 61 74 72 69 78 20 28 61 | 6e 20 65 78 70 65 6e 73 |atrix (a|n expens|
|000012b0| 69 76 65 20 77 61 79 20 | 6f 66 20 64 6f 69 6e 67 |ive way |of doing|
|000012c0| 20 74 68 69 73 29 0a 58 | 20 20 20 4d 61 74 72 69 | this).X| Matri|
|000012d0| 78 20 58 31 28 6e 6f 62 | 73 2c 6e 70 72 65 64 31 |x X1(nob|s,npred1|
|000012e0| 29 3b 20 58 31 2e 43 6f | 6c 75 6d 6e 28 31 29 3c |); X1.Co|lumn(1)<|
|000012f0| 3c 4f 6e 65 73 3b 20 58 | 31 2e 43 6f 6c 75 6d 6e |<Ones; X|1.Column|
|00001300| 73 28 32 2c 6e 70 72 65 | 64 31 29 3c 3c 58 3b 0a |s(2,npre|d1)<<X;.|
|00001310| 58 20 20 20 44 69 61 67 | 6f 6e 61 6c 4d 61 74 72 |X Diag|onalMatr|
|00001320| 69 78 20 48 61 74 3b 20 | 20 48 61 74 20 3c 3c 20 |ix Hat; | Hat << |
|00001330| 58 31 20 2a 20 28 58 31 | 2e 74 28 29 20 2a 20 58 |X1 * (X1|.t() * X|
|00001340| 31 29 2e 69 28 29 20 2a | 20 58 31 2e 74 28 29 3b |1).i() *| X1.t();|
|00001350| 0a 58 0a 58 20 20 20 2f | 2f 20 70 72 69 6e 74 20 |.X.X /|/ print |
|00001360| 6f 75 74 20 61 6e 73 77 | 65 72 73 0a 58 20 20 20 |out answ|ers.X |
|00001370| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 45 73 74 69 6d |cout << |"\nEstim|
|00001380| 61 74 65 73 20 61 6e 64 | 20 74 68 65 69 72 20 73 |ates and| their s|
|00001390| 74 61 6e 64 61 72 64 20 | 65 72 72 6f 72 73 5c 6e |tandard |errors\n|
|000013a0| 5c 6e 22 3b 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c |\n";.X | cout <<|
|000013b0| 20 61 20 3c 3c 22 5c 74 | 22 3c 3c 20 73 71 72 74 | a <<"\t|"<< sqrt|
|000013c0| 28 76 2a 52 65 73 56 61 | 72 29 20 3c 3c 20 22 5c |(v*ResVa|r) << "\|
|000013d0| 6e 22 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 6e 74 |n";.X |for (int|
|000013e0| 20 69 3d 31 3b 20 69 3c | 3d 6e 70 72 65 64 3b 20 | i=1; i<|=npred; |
|000013f0| 69 2b 2b 29 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c |i++).X | cout <<|
|00001400| 20 41 28 69 29 20 3c 3c | 22 5c 74 22 3c 3c 20 73 | A(i) <<|"\t"<< s|
|00001410| 71 72 74 28 56 28 69 29 | 2a 52 65 73 56 61 72 29 |qrt(V(i)|*ResVar)|
|00001420| 20 3c 3c 20 22 5c 6e 22 | 3b 0a 58 20 20 20 63 6f | << "\n"|;.X co|
|00001430| 75 74 20 3c 3c 20 22 5c | 6e 4f 62 73 65 72 76 61 |ut << "\|nObserva|
|00001440| 74 69 6f 6e 73 2c 20 66 | 69 74 74 65 64 20 76 61 |tions, f|itted va|
|00001450| 6c 75 65 2c 20 72 65 73 | 69 64 75 61 6c 20 76 61 |lue, res|idual va|
|00001460| 6c 75 65 2c 20 68 61 74 | 20 76 61 6c 75 65 5c 6e |lue, hat| value\n|
|00001470| 22 3b 0a 58 20 20 20 66 | 6f 72 20 28 69 3d 31 3b |";.X f|or (i=1;|
|00001480| 20 69 3c 3d 6e 6f 62 73 | 3b 20 69 2b 2b 29 0a 58 | i<=nobs|; i++).X|
|00001490| 20 20 20 20 20 20 63 6f | 75 74 20 3c 3c 20 58 28 | co|ut << X(|
|000014a0| 69 2c 31 29 20 3c 3c 22 | 5c 74 22 3c 3c 20 58 28 |i,1) <<"|\t"<< X(|
|000014b0| 69 2c 32 29 20 3c 3c 22 | 5c 74 22 3c 3c 20 59 28 |i,2) <<"|\t"<< Y(|
|000014c0| 69 29 20 3c 3c 22 5c 74 | 22 3c 3c 0a 58 20 20 20 |i) <<"\t|"<<.X |
|000014d0| 20 20 20 74 33 28 46 69 | 74 74 65 64 28 69 29 29 | t3(Fi|tted(i))|
|000014e0| 20 3c 3c 22 5c 74 22 3c | 3c 20 74 33 28 52 65 73 | <<"\t"<|< t3(Res|
|000014f0| 69 64 75 61 6c 28 69 29 | 29 20 3c 3c 22 5c 74 22 |idual(i)|) <<"\t"|
|00001500| 3c 3c 20 74 33 28 48 61 | 74 28 69 29 29 20 3c 3c |<< t3(Ha|t(i)) <<|
|00001510| 22 5c 6e 22 3b 0a 58 20 | 20 20 63 6f 75 74 20 3c |"\n";.X | cout <|
|00001520| 3c 20 22 5c 6e 5c 6e 22 | 3b 0a 58 7d 0a 58 0a 58 |< "\n\n"|;.X}.X.X|
|00001530| 76 6f 69 64 20 74 65 73 | 74 32 28 52 65 61 6c 2a |void tes|t2(Real*|
|00001540| 20 79 2c 20 52 65 61 6c | 2a 20 78 31 2c 20 52 65 | y, Real|* x1, Re|
|00001550| 61 6c 2a 20 78 32 2c 20 | 69 6e 74 20 6e 6f 62 73 |al* x2, |int nobs|
|00001560| 2c 20 69 6e 74 20 6e 70 | 72 65 64 29 0a 58 7b 0a |, int np|red).X{.|
|00001570| 58 20 20 20 63 6f 75 74 | 20 3c 3c 20 22 5c 6e 5c |X cout| << "\n\|
|00001580| 6e 54 65 73 74 20 32 20 | 2d 20 43 68 6f 6c 65 73 |nTest 2 |- Choles|
|00001590| 6b 79 5c 6e 22 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 |ky\n";.X|.X // |
|000015a0| 74 72 61 64 69 74 69 6f | 6e 61 6c 20 73 75 6d 20 |traditio|nal sum |
|000015b0| 6f 66 20 73 71 75 61 72 | 65 73 20 61 6e 64 20 70 |of squar|es and p|
|000015c0| 72 6f 64 75 63 74 73 20 | 6d 65 74 68 6f 64 20 6f |roducts |method o|
|000015d0| 66 20 63 61 6c 63 75 6c | 61 74 69 6f 6e 0a 58 20 |f calcul|ation.X |
|000015e0| 20 20 2f 2f 20 77 69 74 | 68 20 73 75 62 74 72 61 | // wit|h subtra|
|000015f0| 63 74 69 6f 6e 20 6f 66 | 20 6d 65 61 6e 73 20 2d |ction of| means -|
|00001600| 20 75 73 69 6e 67 20 43 | 68 6f 6c 65 73 6b 79 20 | using C|holesky |
|00001610| 64 65 63 6f 6d 70 6f 73 | 69 74 69 6f 6e 0a 58 0a |decompos|ition.X.|
|00001620| 58 20 20 20 4d 61 74 72 | 69 78 20 58 28 6e 6f 62 |X Matr|ix X(nob|
|00001630| 73 2c 6e 70 72 65 64 29 | 3b 0a 58 20 20 20 58 2e |s,npred)|;.X X.|
|00001640| 43 6f 6c 75 6d 6e 28 31 | 29 20 3c 3c 20 78 31 3b |Column(1|) << x1;|
|00001650| 20 20 58 2e 43 6f 6c 75 | 6d 6e 28 32 29 20 3c 3c | X.Colu|mn(2) <<|
|00001660| 20 78 32 3b 0a 58 20 20 | 20 43 6f 6c 75 6d 6e 56 | x2;.X | ColumnV|
|00001670| 65 63 74 6f 72 20 59 28 | 6e 6f 62 73 29 3b 20 59 |ector Y(|nobs); Y|
|00001680| 20 3c 3c 20 79 3b 0a 58 | 20 20 20 43 6f 6c 75 6d | << y;.X| Colum|
|00001690| 6e 56 65 63 74 6f 72 20 | 4f 6e 65 73 28 6e 6f 62 |nVector |Ones(nob|
|000016a0| 73 29 3b 20 4f 6e 65 73 | 20 3d 20 31 2e 30 3b 0a |s); Ones| = 1.0;.|
|000016b0| 58 20 20 20 52 6f 77 56 | 65 63 74 6f 72 20 4d 20 |X RowV|ector M |
|000016c0| 3d 20 4f 6e 65 73 2e 74 | 28 29 20 2a 20 58 20 2f |= Ones.t|() * X /|
|000016d0| 20 6e 6f 62 73 3b 0a 58 | 20 20 20 4d 61 74 72 69 | nobs;.X| Matri|
|000016e0| 78 20 58 43 28 6e 6f 62 | 73 2c 6e 70 72 65 64 29 |x XC(nob|s,npred)|
|000016f0| 3b 0a 58 20 20 20 58 43 | 20 3d 20 58 20 2d 20 4f |;.X XC| = X - O|
|00001700| 6e 65 73 20 2a 20 4d 3b | 0a 58 20 20 20 43 6f 6c |nes * M;|.X Col|
|00001710| 75 6d 6e 56 65 63 74 6f | 72 20 59 43 28 6e 6f 62 |umnVecto|r YC(nob|
|00001720| 73 29 3b 0a 58 20 20 20 | 52 65 61 6c 20 6d 20 3d |s);.X |Real m =|
|00001730| 20 28 4f 6e 65 73 2e 74 | 28 29 20 2a 20 59 29 2e | (Ones.t|() * Y).|
|00001740| 41 73 53 63 61 6c 61 72 | 28 29 20 2f 20 6e 6f 62 |AsScalar|() / nob|
|00001750| 73 3b 20 20 59 43 20 3d | 20 59 20 2d 20 4f 6e 65 |s; YC =| Y - One|
|00001760| 73 20 2a 20 6d 3b 0a 58 | 20 20 20 53 79 6d 6d 65 |s * m;.X| Symme|
|00001770| 74 72 69 63 4d 61 74 72 | 69 78 20 53 53 51 3b 20 |tricMatr|ix SSQ; |
|00001780| 53 53 51 20 3c 3c 20 58 | 43 2e 74 28 29 20 2a 20 |SSQ << X|C.t() * |
|00001790| 58 43 3b 0a 58 0a 58 20 | 20 20 2f 2f 20 43 68 6f |XC;.X.X | // Cho|
|000017a0| 6c 65 73 6b 79 20 64 65 | 63 6f 6d 70 6f 73 69 74 |lesky de|composit|
|000017b0| 69 6f 6e 20 6f 66 20 53 | 53 51 0a 58 20 20 20 4c |ion of S|SQ.X L|
|000017c0| 6f 77 65 72 54 72 69 61 | 6e 67 75 6c 61 72 4d 61 |owerTria|ngularMa|
|000017d0| 74 72 69 78 20 4c 20 3d | 20 43 68 6f 6c 65 73 6b |trix L =| Cholesk|
|000017e0| 79 28 53 53 51 29 3b 0a | 58 0a 58 20 20 20 2f 2f |y(SSQ);.|X.X //|
|000017f0| 20 63 61 6c 63 75 6c 61 | 74 65 20 65 73 74 69 6d | calcula|te estim|
|00001800| 61 74 65 0a 58 20 20 20 | 43 6f 6c 75 6d 6e 56 65 |ate.X |ColumnVe|
|00001810| 63 74 6f 72 20 41 20 3d | 20 4c 2e 74 28 29 2e 69 |ctor A =| L.t().i|
|00001820| 28 29 20 2a 20 28 4c 2e | 69 28 29 20 2a 20 28 58 |() * (L.|i() * (X|
|00001830| 43 2e 74 28 29 20 2a 20 | 59 43 29 29 3b 0a 58 0a |C.t() * |YC));.X.|
|00001840| 58 20 20 20 2f 2f 20 63 | 61 6c 63 75 6c 61 74 65 |X // c|alculate|
|00001850| 20 65 73 74 69 6d 61 74 | 65 20 6f 66 20 63 6f 6e | estimat|e of con|
|00001860| 73 74 61 6e 74 20 74 65 | 72 6d 0a 58 20 20 20 52 |stant te|rm.X R|
|00001870| 65 61 6c 20 61 20 3d 20 | 6d 20 2d 20 28 4d 20 2a |eal a = |m - (M *|
|00001880| 20 41 29 2e 41 73 53 63 | 61 6c 61 72 28 29 3b 0a | A).AsSc|alar();.|
|00001890| 58 0a 58 20 20 20 2f 2f | 20 47 65 74 20 76 61 72 |X.X //| Get var|
|000018a0| 69 61 6e 63 65 73 20 6f | 66 20 65 73 74 69 6d 61 |iances o|f estima|
|000018b0| 74 65 73 20 66 72 6f 6d | 20 64 69 61 67 6f 6e 61 |tes from| diagona|
|000018c0| 6c 20 65 6c 65 6d 65 6e | 74 73 20 6f 66 20 69 6e |l elemen|ts of in|
|000018d0| 76 6f 69 63 65 20 6f 66 | 20 53 53 51 0a 58 20 20 |voice of| SSQ.X |
|000018e0| 20 44 69 61 67 6f 6e 61 | 6c 4d 61 74 72 69 78 20 | Diagona|lMatrix |
|000018f0| 44 3b 20 44 20 3c 3c 20 | 4c 2e 74 28 29 2e 69 28 |D; D << |L.t().i(|
|00001900| 29 20 2a 20 4c 2e 69 28 | 29 3b 0a 58 20 20 20 43 |) * L.i(|);.X C|
|00001910| 6f 6c 75 6d 6e 56 65 63 | 74 6f 72 20 56 20 3d 20 |olumnVec|tor V = |
|00001920| 44 2e 41 73 43 6f 6c 75 | 6d 6e 28 29 3b 0a 58 20 |D.AsColu|mn();.X |
|00001930| 20 20 52 65 61 6c 20 76 | 20 3d 20 31 2e 30 2f 6e | Real v| = 1.0/n|
|00001940| 6f 62 73 20 2b 20 28 4c | 2e 69 28 29 20 2a 20 4d |obs + (L|.i() * M|
|00001950| 2e 74 28 29 29 2e 53 75 | 6d 53 71 75 61 72 65 28 |.t()).Su|mSquare(|
|00001960| 29 3b 0a 58 0a 58 20 20 | 20 2f 2f 20 43 61 6c 63 |);.X.X | // Calc|
|00001970| 75 6c 61 74 65 20 66 69 | 74 74 65 64 20 76 61 6c |ulate fi|tted val|
|00001980| 75 65 73 20 61 6e 64 20 | 72 65 73 69 64 75 61 6c |ues and |residual|
|00001990| 73 0a 58 20 20 20 69 6e | 74 20 6e 70 72 65 64 31 |s.X in|t npred1|
|000019a0| 20 3d 20 6e 70 72 65 64 | 2b 31 3b 0a 58 20 20 20 | = npred|+1;.X |
|000019b0| 43 6f 6c 75 6d 6e 56 65 | 63 74 6f 72 20 46 69 74 |ColumnVe|ctor Fit|
|000019c0| 74 65 64 20 3d 20 58 20 | 2a 20 41 20 2b 20 61 3b |ted = X |* A + a;|
|000019d0| 0a 58 20 20 20 43 6f 6c | 75 6d 6e 56 65 63 74 6f |.X Col|umnVecto|
|000019e0| 72 20 52 65 73 69 64 75 | 61 6c 20 3d 20 59 20 2d |r Residu|al = Y -|
|000019f0| 20 46 69 74 74 65 64 3b | 0a 58 20 20 20 52 65 61 | Fitted;|.X Rea|
|00001a00| 6c 20 52 65 73 56 61 72 | 20 3d 20 52 65 73 69 64 |l ResVar| = Resid|
|00001a10| 75 61 6c 2e 53 75 6d 53 | 71 75 61 72 65 28 29 20 |ual.SumS|quare() |
|00001a20| 2f 20 28 6e 6f 62 73 2d | 6e 70 72 65 64 31 29 3b |/ (nobs-|npred1);|
|00001a30| 0a 58 0a 58 20 20 20 2f | 2f 20 47 65 74 20 64 69 |.X.X /|/ Get di|
|00001a40| 61 67 6f 6e 61 6c 73 20 | 6f 66 20 48 61 74 20 6d |agonals |of Hat m|
|00001a50| 61 74 72 69 78 20 28 61 | 6e 20 65 78 70 65 6e 73 |atrix (a|n expens|
|00001a60| 69 76 65 20 77 61 79 20 | 6f 66 20 64 6f 69 6e 67 |ive way |of doing|
|00001a70| 20 74 68 69 73 29 0a 58 | 20 20 20 4d 61 74 72 69 | this).X| Matri|
|00001a80| 78 20 58 31 28 6e 6f 62 | 73 2c 6e 70 72 65 64 31 |x X1(nob|s,npred1|
|00001a90| 29 3b 20 58 31 2e 43 6f | 6c 75 6d 6e 28 31 29 3c |); X1.Co|lumn(1)<|
|00001aa0| 3c 4f 6e 65 73 3b 20 58 | 31 2e 43 6f 6c 75 6d 6e |<Ones; X|1.Column|
|00001ab0| 73 28 32 2c 6e 70 72 65 | 64 31 29 3c 3c 58 3b 0a |s(2,npre|d1)<<X;.|
|00001ac0| 58 20 20 20 44 69 61 67 | 6f 6e 61 6c 4d 61 74 72 |X Diag|onalMatr|
|00001ad0| 69 78 20 48 61 74 3b 20 | 20 48 61 74 20 3c 3c 20 |ix Hat; | Hat << |
|00001ae0| 58 31 20 2a 20 28 58 31 | 2e 74 28 29 20 2a 20 58 |X1 * (X1|.t() * X|
|00001af0| 31 29 2e 69 28 29 20 2a | 20 58 31 2e 74 28 29 3b |1).i() *| X1.t();|
|00001b00| 0a 58 0a 58 20 20 20 2f | 2f 20 70 72 69 6e 74 20 |.X.X /|/ print |
|00001b10| 6f 75 74 20 61 6e 73 77 | 65 72 73 0a 58 20 20 20 |out answ|ers.X |
|00001b20| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 45 73 74 69 6d |cout << |"\nEstim|
|00001b30| 61 74 65 73 20 61 6e 64 | 20 74 68 65 69 72 20 73 |ates and| their s|
|00001b40| 74 61 6e 64 61 72 64 20 | 65 72 72 6f 72 73 5c 6e |tandard |errors\n|
|00001b50| 5c 6e 22 3b 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c |\n";.X | cout <<|
|00001b60| 20 61 20 3c 3c 22 5c 74 | 22 3c 3c 20 73 71 72 74 | a <<"\t|"<< sqrt|
|00001b70| 28 76 2a 52 65 73 56 61 | 72 29 20 3c 3c 20 22 5c |(v*ResVa|r) << "\|
|00001b80| 6e 22 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 6e 74 |n";.X |for (int|
|00001b90| 20 69 3d 31 3b 20 69 3c | 3d 6e 70 72 65 64 3b 20 | i=1; i<|=npred; |
|00001ba0| 69 2b 2b 29 0a 58 20 20 | 20 20 20 20 63 6f 75 74 |i++).X | cout|
|00001bb0| 20 3c 3c 20 41 28 69 29 | 20 3c 3c 22 5c 74 22 3c | << A(i)| <<"\t"<|
|00001bc0| 3c 20 73 71 72 74 28 56 | 28 69 29 2a 52 65 73 56 |< sqrt(V|(i)*ResV|
|00001bd0| 61 72 29 20 3c 3c 20 22 | 5c 6e 22 3b 0a 58 20 20 |ar) << "|\n";.X |
|00001be0| 20 63 6f 75 74 20 3c 3c | 20 22 5c 6e 4f 62 73 65 | cout <<| "\nObse|
|00001bf0| 72 76 61 74 69 6f 6e 73 | 2c 20 66 69 74 74 65 64 |rvations|, fitted|
|00001c00| 20 76 61 6c 75 65 2c 20 | 72 65 73 69 64 75 61 6c | value, |residual|
|00001c10| 20 76 61 6c 75 65 2c 20 | 68 61 74 20 76 61 6c 75 | value, |hat valu|
|00001c20| 65 5c 6e 22 3b 0a 58 20 | 20 20 66 6f 72 20 28 69 |e\n";.X | for (i|
|00001c30| 3d 31 3b 20 69 3c 3d 6e | 6f 62 73 3b 20 69 2b 2b |=1; i<=n|obs; i++|
|00001c40| 29 0a 58 20 20 20 20 20 | 20 63 6f 75 74 20 3c 3c |).X | cout <<|
|00001c50| 20 58 28 69 2c 31 29 20 | 3c 3c 22 5c 74 22 3c 3c | X(i,1) |<<"\t"<<|
|00001c60| 20 58 28 69 2c 32 29 20 | 3c 3c 22 5c 74 22 3c 3c | X(i,2) |<<"\t"<<|
|00001c70| 20 59 28 69 29 20 3c 3c | 22 5c 74 22 3c 3c 0a 58 | Y(i) <<|"\t"<<.X|
|00001c80| 20 20 20 20 20 20 74 33 | 28 46 69 74 74 65 64 28 | t3|(Fitted(|
|00001c90| 69 29 29 20 3c 3c 22 5c | 74 22 3c 3c 20 74 33 28 |i)) <<"\|t"<< t3(|
|00001ca0| 52 65 73 69 64 75 61 6c | 28 69 29 29 20 3c 3c 22 |Residual|(i)) <<"|
|00001cb0| 5c 74 22 3c 3c 20 74 33 | 28 48 61 74 28 69 29 29 |\t"<< t3|(Hat(i))|
|00001cc0| 20 3c 3c 22 5c 6e 22 3b | 0a 58 20 20 20 63 6f 75 | <<"\n";|.X cou|
|00001cd0| 74 20 3c 3c 20 22 5c 6e | 5c 6e 22 3b 0a 58 7d 0a |t << "\n|\n";.X}.|
|00001ce0| 58 0a 58 76 6f 69 64 20 | 74 65 73 74 33 28 52 65 |X.Xvoid |test3(Re|
|00001cf0| 61 6c 2a 20 79 2c 20 52 | 65 61 6c 2a 20 78 31 2c |al* y, R|eal* x1,|
|00001d00| 20 52 65 61 6c 2a 20 78 | 32 2c 20 69 6e 74 20 6e | Real* x|2, int n|
|00001d10| 6f 62 73 2c 20 69 6e 74 | 20 6e 70 72 65 64 29 0a |obs, int| npred).|
|00001d20| 58 7b 0a 58 20 20 20 63 | 6f 75 74 20 3c 3c 20 22 |X{.X c|out << "|
|00001d30| 5c 6e 5c 6e 54 65 73 74 | 20 33 20 2d 20 48 6f 75 |\n\nTest| 3 - Hou|
|00001d40| 73 65 68 6f 6c 64 65 72 | 20 74 72 69 61 6e 67 75 |seholder| triangu|
|00001d50| 6c 61 72 69 73 61 74 69 | 6f 6e 5c 6e 22 3b 0a 58 |larisati|on\n";.X|
|00001d60| 0a 58 20 20 20 2f 2f 20 | 48 6f 75 73 65 68 6f 6c |.X // |Househol|
|00001d70| 64 65 72 20 74 72 69 61 | 6e 67 75 6c 61 72 69 73 |der tria|ngularis|
|00001d80| 61 74 69 6f 6e 20 6d 65 | 74 68 6f 64 0a 58 20 0a |ation me|thod.X .|
|00001d90| 58 20 20 20 2f 2f 20 6c | 6f 61 64 20 64 61 74 61 |X // l|oad data|
|00001da0| 20 2d 20 31 73 20 69 6e | 74 6f 20 63 6f 6c 20 31 | - 1s in|to col 1|
|00001db0| 20 6f 66 20 6d 61 74 72 | 69 78 0a 58 20 20 20 69 | of matr|ix.X i|
|00001dc0| 6e 74 20 6e 70 72 65 64 | 31 20 3d 20 6e 70 72 65 |nt npred|1 = npre|
|00001dd0| 64 2b 31 3b 0a 58 20 20 | 20 4d 61 74 72 69 78 20 |d+1;.X | Matrix |
|00001de0| 58 28 6e 6f 62 73 2c 6e | 70 72 65 64 31 29 3b 20 |X(nobs,n|pred1); |
|00001df0| 43 6f 6c 75 6d 6e 56 65 | 63 74 6f 72 20 59 28 6e |ColumnVe|ctor Y(n|
|00001e00| 6f 62 73 29 3b 0a 58 20 | 20 20 58 2e 43 6f 6c 75 |obs);.X | X.Colu|
|00001e10| 6d 6e 28 31 29 20 3d 20 | 31 2e 30 3b 20 20 58 2e |mn(1) = |1.0; X.|
|00001e20| 43 6f 6c 75 6d 6e 28 32 | 29 20 3c 3c 20 78 31 3b |Column(2|) << x1;|
|00001e30| 20 20 58 2e 43 6f 6c 75 | 6d 6e 28 33 29 20 3c 3c | X.Colu|mn(3) <<|
|00001e40| 20 78 32 3b 20 20 59 20 | 3c 3c 20 79 3b 0a 58 0a | x2; Y |<< y;.X.|
|00001e50| 58 20 20 20 2f 2f 20 64 | 6f 20 48 6f 75 73 65 68 |X // d|o Househ|
|00001e60| 6f 6c 64 65 72 20 74 72 | 69 61 6e 67 75 6c 61 72 |older tr|iangular|
|00001e70| 69 73 61 74 69 6f 6e 0a | 58 20 20 20 2f 2f 20 6e |isation.|X // n|
|00001e80| 6f 20 6e 65 65 64 20 74 | 6f 20 64 65 61 6c 20 77 |o need t|o deal w|
|00001e90| 69 74 68 20 63 6f 6e 73 | 74 61 6e 74 20 74 65 72 |ith cons|tant ter|
|00001ea0| 6d 20 73 65 70 61 72 61 | 74 65 6c 79 0a 58 20 20 |m separa|tely.X |
|00001eb0| 20 4d 61 74 72 69 78 20 | 58 54 20 3d 20 58 2e 74 | Matrix |XT = X.t|
|00001ec0| 28 29 3b 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |(); | |
|00001ed0| 2f 2f 20 57 61 6e 74 20 | 64 61 74 61 20 74 6f 20 |// Want |data to |
|00001ee0| 62 65 20 61 6c 6f 6e 67 | 20 72 6f 77 73 0a 58 20 |be along| rows.X |
|00001ef0| 20 20 52 6f 77 56 65 63 | 74 6f 72 20 59 54 20 3d | RowVec|tor YT =|
|00001f00| 20 59 2e 74 28 29 3b 0a | 58 20 20 20 4c 6f 77 65 | Y.t();.|X Lowe|
|00001f10| 72 54 72 69 61 6e 67 75 | 6c 61 72 4d 61 74 72 69 |rTriangu|larMatri|
|00001f20| 78 20 4c 3b 20 52 6f 77 | 56 65 63 74 6f 72 20 4d |x L; Row|Vector M|
|00001f30| 3b 0a 58 20 20 20 48 48 | 44 65 63 6f 6d 70 6f 73 |;.X HH|Decompos|
|00001f40| 65 28 58 54 2c 20 4c 29 | 3b 20 48 48 44 65 63 6f |e(XT, L)|; HHDeco|
|00001f50| 6d 70 6f 73 65 28 58 54 | 2c 20 59 54 2c 20 4d 29 |mpose(XT|, YT, M)|
|00001f60| 3b 20 2f 2f 20 59 54 20 | 6e 6f 77 20 63 6f 6e 74 |; // YT |now cont|
|00001f70| 61 69 6e 73 20 72 65 73 | 69 64 73 0a 58 20 20 20 |ains res|ids.X |
|00001f80| 43 6f 6c 75 6d 6e 56 65 | 63 74 6f 72 20 41 20 3d |ColumnVe|ctor A =|
|00001f90| 20 4c 2e 74 28 29 2e 69 | 28 29 20 2a 20 4d 2e 74 | L.t().i|() * M.t|
|00001fa0| 28 29 3b 0a 58 20 20 20 | 43 6f 6c 75 6d 6e 56 65 |();.X |ColumnVe|
|00001fb0| 63 74 6f 72 20 46 69 74 | 74 65 64 20 3d 20 58 20 |ctor Fit|ted = X |
|00001fc0| 2a 20 41 3b 0a 58 20 20 | 20 52 65 61 6c 20 52 65 |* A;.X | Real Re|
|00001fd0| 73 56 61 72 20 3d 20 59 | 54 2e 53 75 6d 53 71 75 |sVar = Y|T.SumSqu|
|00001fe0| 61 72 65 28 29 20 2f 20 | 28 6e 6f 62 73 2d 6e 70 |are() / |(nobs-np|
|00001ff0| 72 65 64 31 29 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 |red1);.X|.X // |
|00002000| 67 65 74 20 76 61 72 69 | 61 6e 63 65 73 20 6f 66 |get vari|ances of|
|00002010| 20 65 73 74 69 6d 61 74 | 65 73 0a 58 20 20 20 4c | estimat|es.X L|
|00002020| 20 3d 20 4c 2e 69 28 29 | 3b 20 44 69 61 67 6f 6e | = L.i()|; Diagon|
|00002030| 61 6c 4d 61 74 72 69 78 | 20 44 3b 20 44 20 3c 3c |alMatrix| D; D <<|
|00002040| 20 4c 2e 74 28 29 20 2a | 20 4c 3b 0a 58 0a 58 20 | L.t() *| L;.X.X |
|00002050| 20 20 2f 2f 20 47 65 74 | 20 64 69 61 67 6f 6e 61 | // Get| diagona|
|00002060| 6c 73 20 6f 66 20 48 61 | 74 20 6d 61 74 72 69 78 |ls of Ha|t matrix|
|00002070| 0a 58 20 20 20 44 69 61 | 67 6f 6e 61 6c 4d 61 74 |.X Dia|gonalMat|
|00002080| 72 69 78 20 48 61 74 3b | 20 20 48 61 74 20 3c 3c |rix Hat;| Hat <<|
|00002090| 20 58 54 2e 74 28 29 20 | 2a 20 58 54 3b 0a 58 0a | XT.t() |* XT;.X.|
|000020a0| 58 20 20 20 2f 2f 20 70 | 72 69 6e 74 20 6f 75 74 |X // p|rint out|
|000020b0| 20 61 6e 73 77 65 72 73 | 0a 58 20 20 20 63 6f 75 | answers|.X cou|
|000020c0| 74 20 3c 3c 20 22 5c 6e | 45 73 74 69 6d 61 74 65 |t << "\n|Estimate|
|000020d0| 73 20 61 6e 64 20 74 68 | 65 69 72 20 73 74 61 6e |s and th|eir stan|
|000020e0| 64 61 72 64 20 65 72 72 | 6f 72 73 5c 6e 5c 6e 22 |dard err|ors\n\n"|
|000020f0| 3b 0a 58 20 20 20 66 6f | 72 20 28 69 6e 74 20 69 |;.X fo|r (int i|
|00002100| 3d 31 3b 20 69 3c 3d 6e | 70 72 65 64 31 3b 20 69 |=1; i<=n|pred1; i|
|00002110| 2b 2b 29 0a 58 20 20 20 | 20 20 20 63 6f 75 74 20 |++).X | cout |
|00002120| 3c 3c 20 41 28 69 29 20 | 3c 3c 22 5c 74 22 3c 3c |<< A(i) |<<"\t"<<|
|00002130| 20 73 71 72 74 28 44 28 | 69 29 2a 52 65 73 56 61 | sqrt(D(|i)*ResVa|
|00002140| 72 29 20 3c 3c 20 22 5c | 6e 22 3b 0a 58 20 20 20 |r) << "\|n";.X |
|00002150| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 4f 62 73 65 72 |cout << |"\nObser|
|00002160| 76 61 74 69 6f 6e 73 2c | 20 66 69 74 74 65 64 20 |vations,| fitted |
|00002170| 76 61 6c 75 65 2c 20 72 | 65 73 69 64 75 61 6c 20 |value, r|esidual |
|00002180| 76 61 6c 75 65 2c 20 68 | 61 74 20 76 61 6c 75 65 |value, h|at value|
|00002190| 5c 6e 22 3b 0a 58 20 20 | 20 66 6f 72 20 28 69 3d |\n";.X | for (i=|
|000021a0| 31 3b 20 69 3c 3d 6e 6f | 62 73 3b 20 69 2b 2b 29 |1; i<=no|bs; i++)|
|000021b0| 0a 58 20 20 20 20 20 20 | 63 6f 75 74 20 3c 3c 20 |.X |cout << |
|000021c0| 58 28 69 2c 32 29 20 3c | 3c 22 5c 74 22 3c 3c 20 |X(i,2) <|<"\t"<< |
|000021d0| 58 28 69 2c 33 29 20 3c | 3c 22 5c 74 22 3c 3c 20 |X(i,3) <|<"\t"<< |
|000021e0| 59 28 69 29 20 3c 3c 22 | 5c 74 22 3c 3c 0a 58 20 |Y(i) <<"|\t"<<.X |
|000021f0| 20 20 20 20 20 74 33 28 | 46 69 74 74 65 64 28 69 | t3(|Fitted(i|
|00002200| 29 29 20 3c 3c 22 5c 74 | 22 3c 3c 20 74 33 28 59 |)) <<"\t|"<< t3(Y|
|00002210| 54 28 69 29 29 20 3c 3c | 22 5c 74 22 3c 3c 20 74 |T(i)) <<|"\t"<< t|
|00002220| 33 28 48 61 74 28 69 29 | 29 20 3c 3c 22 5c 6e 22 |3(Hat(i)|) <<"\n"|
|00002230| 3b 0a 58 20 20 20 63 6f | 75 74 20 3c 3c 20 22 5c |;.X co|ut << "\|
|00002240| 6e 5c 6e 22 3b 0a 58 7d | 0a 58 0a 58 76 6f 69 64 |n\n";.X}|.X.Xvoid|
|00002250| 20 74 65 73 74 34 28 52 | 65 61 6c 2a 20 79 2c 20 | test4(R|eal* y, |
|00002260| 52 65 61 6c 2a 20 78 31 | 2c 20 52 65 61 6c 2a 20 |Real* x1|, Real* |
|00002270| 78 32 2c 20 69 6e 74 20 | 6e 6f 62 73 2c 20 69 6e |x2, int |nobs, in|
|00002280| 74 20 6e 70 72 65 64 29 | 0a 58 7b 0a 58 20 20 20 |t npred)|.X{.X |
|00002290| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 5c 6e 54 65 73 |cout << |"\n\nTes|
|000022a0| 74 20 34 20 2d 20 73 69 | 6e 67 75 6c 61 72 20 76 |t 4 - si|ngular v|
|000022b0| 61 6c 75 65 5c 6e 22 3b | 0a 58 0a 58 20 20 20 2f |alue\n";|.X.X /|
|000022c0| 2f 20 53 69 6e 67 75 6c | 61 72 20 76 61 6c 75 65 |/ Singul|ar value|
|000022d0| 20 64 65 63 6f 6d 70 6f | 73 69 74 69 6f 6e 20 6d | decompo|sition m|
|000022e0| 65 74 68 6f 64 0a 58 20 | 0a 58 20 20 20 2f 2f 20 |ethod.X |.X // |
|000022f0| 6c 6f 61 64 20 64 61 74 | 61 20 2d 20 31 73 20 69 |load dat|a - 1s i|
|00002300| 6e 74 6f 20 63 6f 6c 20 | 31 20 6f 66 20 6d 61 74 |nto col |1 of mat|
|00002310| 72 69 78 0a 58 20 20 20 | 69 6e 74 20 6e 70 72 65 |rix.X |int npre|
|00002320| 64 31 20 3d 20 6e 70 72 | 65 64 2b 31 3b 0a 58 20 |d1 = npr|ed+1;.X |
|00002330| 20 20 4d 61 74 72 69 78 | 20 58 28 6e 6f 62 73 2c | Matrix| X(nobs,|
|00002340| 6e 70 72 65 64 31 29 3b | 20 43 6f 6c 75 6d 6e 56 |npred1);| ColumnV|
|00002350| 65 63 74 6f 72 20 59 28 | 6e 6f 62 73 29 3b 0a 58 |ector Y(|nobs);.X|
|00002360| 20 20 20 58 2e 43 6f 6c | 75 6d 6e 28 31 29 20 3d | X.Col|umn(1) =|
|00002370| 20 31 2e 30 3b 20 20 58 | 2e 43 6f 6c 75 6d 6e 28 | 1.0; X|.Column(|
|00002380| 32 29 20 3c 3c 20 78 31 | 3b 20 20 58 2e 43 6f 6c |2) << x1|; X.Col|
|00002390| 75 6d 6e 28 33 29 20 3c | 3c 20 78 32 3b 20 20 59 |umn(3) <|< x2; Y|
|000023a0| 20 3c 3c 20 79 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 | << y;.X|.X // |
|000023b0| 64 6f 20 53 56 44 0a 58 | 20 20 20 4d 61 74 72 69 |do SVD.X| Matri|
|000023c0| 78 20 55 2c 20 56 3b 20 | 44 69 61 67 6f 6e 61 6c |x U, V; |Diagonal|
|000023d0| 4d 61 74 72 69 78 20 44 | 3b 0a 58 20 20 20 53 56 |Matrix D|;.X SV|
|000023e0| 44 28 58 2c 44 2c 55 2c | 56 29 3b 20 20 20 20 20 |D(X,D,U,|V); |
|000023f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002400| 20 20 20 20 20 20 20 20 | 20 2f 2f 20 58 20 3d 20 | | // X = |
|00002410| 55 20 2a 20 44 20 2a 20 | 56 2e 74 28 29 0a 58 20 |U * D * |V.t().X |
|00002420| 20 20 43 6f 6c 75 6d 6e | 56 65 63 74 6f 72 20 46 | Column|Vector F|
|00002430| 69 74 74 65 64 20 3d 20 | 55 2e 74 28 29 20 2a 20 |itted = |U.t() * |
|00002440| 59 3b 0a 58 20 20 20 43 | 6f 6c 75 6d 6e 56 65 63 |Y;.X C|olumnVec|
|00002450| 74 6f 72 20 41 20 3d 20 | 56 20 2a 20 28 20 44 2e |tor A = |V * ( D.|
|00002460| 69 28 29 20 2a 20 46 69 | 74 74 65 64 20 29 3b 0a |i() * Fi|tted );.|
|00002470| 58 20 20 20 46 69 74 74 | 65 64 20 3d 20 55 20 2a |X Fitt|ed = U *|
|00002480| 20 46 69 74 74 65 64 3b | 0a 58 20 20 20 43 6f 6c | Fitted;|.X Col|
|00002490| 75 6d 6e 56 65 63 74 6f | 72 20 52 65 73 69 64 75 |umnVecto|r Residu|
|000024a0| 61 6c 20 3d 20 59 20 2d | 20 46 69 74 74 65 64 3b |al = Y -| Fitted;|
|000024b0| 0a 58 20 20 20 52 65 61 | 6c 20 52 65 73 56 61 72 |.X Rea|l ResVar|
|000024c0| 20 3d 20 52 65 73 69 64 | 75 61 6c 2e 53 75 6d 53 | = Resid|ual.SumS|
|000024d0| 71 75 61 72 65 28 29 20 | 2f 20 28 6e 6f 62 73 2d |quare() |/ (nobs-|
|000024e0| 6e 70 72 65 64 31 29 3b | 0a 58 0a 58 20 20 20 2f |npred1);|.X.X /|
|000024f0| 2f 20 67 65 74 20 76 61 | 72 69 61 6e 63 65 73 20 |/ get va|riances |
|00002500| 6f 66 20 65 73 74 69 6d | 61 74 65 73 0a 58 20 20 |of estim|ates.X |
|00002510| 20 44 20 3c 3c 20 56 20 | 2a 20 28 44 20 2a 20 44 | D << V |* (D * D|
|00002520| 29 2e 69 28 29 20 2a 20 | 56 2e 74 28 29 3b 0a 58 |).i() * |V.t();.X|
|00002530| 0a 58 20 20 20 2f 2f 20 | 47 65 74 20 64 69 61 67 |.X // |Get diag|
|00002540| 6f 6e 61 6c 73 20 6f 66 | 20 48 61 74 20 6d 61 74 |onals of| Hat mat|
|00002550| 72 69 78 0a 58 20 20 20 | 44 69 61 67 6f 6e 61 6c |rix.X |Diagonal|
|00002560| 4d 61 74 72 69 78 20 48 | 61 74 3b 20 20 48 61 74 |Matrix H|at; Hat|
|00002570| 20 3c 3c 20 55 20 2a 20 | 55 2e 74 28 29 3b 0a 58 | << U * |U.t();.X|
|00002580| 0a 58 20 20 20 2f 2f 20 | 70 72 69 6e 74 20 6f 75 |.X // |print ou|
|00002590| 74 20 61 6e 73 77 65 72 | 73 0a 58 20 20 20 63 6f |t answer|s.X co|
|000025a0| 75 74 20 3c 3c 20 22 5c | 6e 45 73 74 69 6d 61 74 |ut << "\|nEstimat|
|000025b0| 65 73 20 61 6e 64 20 74 | 68 65 69 72 20 73 74 61 |es and t|heir sta|
|000025c0| 6e 64 61 72 64 20 65 72 | 72 6f 72 73 5c 6e 5c 6e |ndard er|rors\n\n|
|000025d0| 22 3b 0a 58 20 20 20 66 | 6f 72 20 28 69 6e 74 20 |";.X f|or (int |
|000025e0| 69 3d 31 3b 20 69 3c 3d | 6e 70 72 65 64 31 3b 20 |i=1; i<=|npred1; |
|000025f0| 69 2b 2b 29 0a 58 20 20 | 20 20 20 20 63 6f 75 74 |i++).X | cout|
|00002600| 20 3c 3c 20 41 28 69 29 | 20 3c 3c 22 5c 74 22 3c | << A(i)| <<"\t"<|
|00002610| 3c 20 73 71 72 74 28 44 | 28 69 29 2a 52 65 73 56 |< sqrt(D|(i)*ResV|
|00002620| 61 72 29 20 3c 3c 20 22 | 5c 6e 22 3b 0a 58 20 20 |ar) << "|\n";.X |
|00002630| 20 63 6f 75 74 20 3c 3c | 20 22 5c 6e 4f 62 73 65 | cout <<| "\nObse|
|00002640| 72 76 61 74 69 6f 6e 73 | 2c 20 66 69 74 74 65 64 |rvations|, fitted|
|00002650| 20 76 61 6c 75 65 2c 20 | 72 65 73 69 64 75 61 6c | value, |residual|
|00002660| 20 76 61 6c 75 65 2c 20 | 68 61 74 20 76 61 6c 75 | value, |hat valu|
|00002670| 65 5c 6e 22 3b 0a 58 20 | 20 20 66 6f 72 20 28 69 |e\n";.X | for (i|
|00002680| 3d 31 3b 20 69 3c 3d 6e | 6f 62 73 3b 20 69 2b 2b |=1; i<=n|obs; i++|
|00002690| 29 0a 58 20 20 20 20 20 | 20 63 6f 75 74 20 3c 3c |).X | cout <<|
|000026a0| 20 58 28 69 2c 32 29 20 | 3c 3c 22 5c 74 22 3c 3c | X(i,2) |<<"\t"<<|
|000026b0| 20 58 28 69 2c 33 29 20 | 3c 3c 22 5c 74 22 3c 3c | X(i,3) |<<"\t"<<|
|000026c0| 20 59 28 69 29 20 3c 3c | 22 5c 74 22 3c 3c 0a 58 | Y(i) <<|"\t"<<.X|
|000026d0| 20 20 20 20 20 20 74 33 | 28 46 69 74 74 65 64 28 | t3|(Fitted(|
|000026e0| 69 29 29 20 3c 3c 22 5c | 74 22 3c 3c 20 74 33 28 |i)) <<"\|t"<< t3(|
|000026f0| 52 65 73 69 64 75 61 6c | 28 69 29 29 20 3c 3c 22 |Residual|(i)) <<"|
|00002700| 5c 74 22 3c 3c 20 74 33 | 28 48 61 74 28 69 29 29 |\t"<< t3|(Hat(i))|
|00002710| 20 3c 3c 22 5c 6e 22 3b | 0a 58 20 20 20 63 6f 75 | <<"\n";|.X cou|
|00002720| 74 20 3c 3c 20 22 5c 6e | 5c 6e 22 3b 0a 58 7d 0a |t << "\n|\n";.X}.|
|00002730| 58 0a 58 6d 61 69 6e 28 | 29 0a 58 7b 0a 58 20 20 |X.Xmain(|).X{.X |
|00002740| 20 63 6f 75 74 20 3c 3c | 20 22 5c 6e 44 65 6d 6f | cout <<| "\nDemo|
|00002750| 6e 73 74 72 61 74 69 6f | 6e 20 6f 66 20 4d 61 74 |nstratio|n of Mat|
|00002760| 72 69 78 20 70 61 63 6b | 61 67 65 5c 6e 5c 6e 22 |rix pack|age\n\n"|
|00002770| 3b 0a 58 0a 58 20 20 20 | 2f 2f 20 54 65 73 74 20 |;.X.X |// Test |
|00002780| 66 6f 72 20 61 6e 79 20 | 6d 65 6d 6f 72 79 20 6e |for any |memory n|
|00002790| 6f 74 20 64 65 61 6c 6c | 6f 63 61 74 65 64 20 61 |ot deall|ocated a|
|000027a0| 66 74 65 72 20 72 75 6e | 6e 69 6e 67 20 74 68 69 |fter run|ning thi|
|000027b0| 73 20 70 72 6f 67 72 61 | 6d 0a 58 20 20 20 52 65 |s progra|m.X Re|
|000027c0| 61 6c 2a 20 73 31 3b 20 | 7b 20 43 6f 6c 75 6d 6e |al* s1; |{ Column|
|000027d0| 56 65 63 74 6f 72 20 41 | 28 38 30 30 30 29 3b 20 |Vector A|(8000); |
|000027e0| 73 31 20 3d 20 41 2e 53 | 74 6f 72 65 28 29 3b 20 |s1 = A.S|tore(); |
|000027f0| 7d 0a 58 0a 58 20 20 20 | 7b 0a 58 20 20 20 20 20 |}.X.X |{.X |
|00002800| 20 2f 2f 20 74 68 65 20 | 64 61 74 61 0a 58 0a 58 | // the |data.X.X|
|00002810| 23 69 66 6e 64 65 66 20 | 41 54 61 6e 64 54 0a 58 |#ifndef |ATandT.X|
|00002820| 20 20 20 20 20 20 52 65 | 61 6c 20 79 5b 5d 20 20 | Re|al y[] |
|00002830| 3d 20 7b 20 38 2e 33 2c | 20 35 2e 35 2c 20 38 2e |= { 8.3,| 5.5, 8.|
|00002840| 30 2c 20 38 2e 35 2c 20 | 35 2e 37 2c 20 34 2e 34 |0, 8.5, |5.7, 4.4|
|00002850| 2c 20 36 2e 33 2c 20 37 | 2e 39 2c 20 39 2e 31 20 |, 6.3, 7|.9, 9.1 |
|00002860| 7d 3b 0a 58 20 20 20 20 | 20 20 52 65 61 6c 20 78 |};.X | Real x|
|00002870| 31 5b 5d 20 3d 20 7b 20 | 32 2e 34 2c 20 31 2e 38 |1[] = { |2.4, 1.8|
|00002880| 2c 20 32 2e 34 2c 20 33 | 2e 30 2c 20 32 2e 30 2c |, 2.4, 3|.0, 2.0,|
|00002890| 20 31 2e 32 2c 20 32 2e | 30 2c 20 32 2e 37 2c 20 | 1.2, 2.|0, 2.7, |
|000028a0| 33 2e 36 20 7d 3b 0a 58 | 20 20 20 20 20 20 52 65 |3.6 };.X| Re|
|000028b0| 61 6c 20 78 32 5b 5d 20 | 3d 20 7b 20 31 2e 37 2c |al x2[] |= { 1.7,|
|000028c0| 20 30 2e 39 2c 20 31 2e | 36 2c 20 31 2e 39 2c 20 | 0.9, 1.|6, 1.9, |
|000028d0| 30 2e 35 2c 20 30 2e 36 | 2c 20 31 2e 31 2c 20 31 |0.5, 0.6|, 1.1, 1|
|000028e0| 2e 30 2c 20 30 2e 35 20 | 7d 3b 0a 58 23 65 6c 73 |.0, 0.5 |};.X#els|
|000028f0| 65 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2f 2f |e | //|
|00002900| 20 66 6f 72 20 63 6f 6d | 70 69 6c 65 72 73 20 74 | for com|pilers t|
|00002910| 68 61 74 20 64 6f 6e 27 | 74 20 75 6e 64 65 72 73 |hat don'|t unders|
|00002920| 74 61 6e 64 20 61 67 67 | 72 65 67 72 61 74 65 73 |tand agg|regrates|
|00002930| 0a 58 20 20 20 20 20 20 | 52 65 61 6c 20 79 5b 39 |.X |Real y[9|
|00002940| 5d 2c 20 78 31 5b 39 5d | 2c 20 78 32 5b 39 5d 3b |], x1[9]|, x2[9];|
|00002950| 0a 58 20 20 20 20 20 20 | 79 5b 30 5d 3d 38 2e 33 |.X |y[0]=8.3|
|00002960| 3b 20 79 5b 31 5d 3d 35 | 2e 35 3b 20 79 5b 32 5d |; y[1]=5|.5; y[2]|
|00002970| 3d 38 2e 30 3b 20 79 5b | 33 5d 3d 38 2e 35 3b 20 |=8.0; y[|3]=8.5; |
|00002980| 79 5b 34 5d 3d 35 2e 37 | 3b 0a 58 20 20 20 20 20 |y[4]=5.7|;.X |
|00002990| 20 79 5b 35 5d 3d 34 2e | 34 3b 20 79 5b 36 5d 3d | y[5]=4.|4; y[6]=|
|000029a0| 36 2e 33 3b 20 79 5b 37 | 5d 3d 37 2e 39 3b 20 79 |6.3; y[7|]=7.9; y|
|000029b0| 5b 38 5d 3d 39 2e 31 3b | 0a 58 20 20 20 20 20 20 |[8]=9.1;|.X |
|000029c0| 78 31 5b 30 5d 3d 32 2e | 34 3b 20 78 31 5b 31 5d |x1[0]=2.|4; x1[1]|
|000029d0| 3d 31 2e 38 3b 20 78 31 | 5b 32 5d 3d 32 2e 34 3b |=1.8; x1|[2]=2.4;|
|000029e0| 20 78 31 5b 33 5d 3d 33 | 2e 30 3b 20 78 31 5b 34 | x1[3]=3|.0; x1[4|
|000029f0| 5d 3d 32 2e 30 3b 0a 58 | 20 20 20 20 20 20 78 31 |]=2.0;.X| x1|
|00002a00| 5b 35 5d 3d 31 2e 32 3b | 20 78 31 5b 36 5d 3d 32 |[5]=1.2;| x1[6]=2|
|00002a10| 2e 30 3b 20 78 31 5b 37 | 5d 3d 32 2e 37 3b 20 78 |.0; x1[7|]=2.7; x|
|00002a20| 31 5b 38 5d 3d 33 2e 36 | 3b 0a 58 20 20 20 20 20 |1[8]=3.6|;.X |
|00002a30| 20 78 32 5b 30 5d 3d 31 | 2e 37 3b 20 78 32 5b 31 | x2[0]=1|.7; x2[1|
|00002a40| 5d 3d 30 2e 39 3b 20 78 | 32 5b 32 5d 3d 31 2e 36 |]=0.9; x|2[2]=1.6|
|00002a50| 3b 20 78 32 5b 33 5d 3d | 31 2e 39 3b 20 78 32 5b |; x2[3]=|1.9; x2[|
|00002a60| 34 5d 3d 30 2e 35 3b 0a | 58 20 20 20 20 20 20 78 |4]=0.5;.|X x|
|00002a70| 32 5b 35 5d 3d 30 2e 36 | 3b 20 78 32 5b 36 5d 3d |2[5]=0.6|; x2[6]=|
|00002a80| 31 2e 31 3b 20 78 32 5b | 37 5d 3d 31 2e 30 3b 20 |1.1; x2[|7]=1.0; |
|00002a90| 78 32 5b 38 5d 3d 30 2e | 35 3b 0a 58 23 65 6e 64 |x2[8]=0.|5;.X#end|
|00002aa0| 69 66 0a 58 0a 58 20 20 | 20 20 20 20 69 6e 74 20 |if.X.X | int |
|00002ab0| 6e 6f 62 73 20 3d 20 39 | 3b 20 20 20 20 20 20 20 |nobs = 9|; |
|00002ac0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002ad0| 20 20 20 20 2f 2f 20 6e | 75 6d 62 65 72 20 6f 66 | // n|umber of|
|00002ae0| 20 6f 62 73 65 72 76 61 | 74 69 6f 6e 73 0a 58 20 | observa|tions.X |
|00002af0| 20 20 20 20 20 69 6e 74 | 20 6e 70 72 65 64 20 3d | int| npred =|
|00002b00| 20 32 3b 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | 2; | |
|00002b10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2f 2f 20 | | // |
|00002b20| 6e 75 6d 62 65 72 20 6f | 66 20 70 72 65 64 69 63 |number o|f predic|
|00002b30| 74 6f 72 20 76 61 6c 75 | 65 73 0a 58 0a 58 20 20 |tor valu|es.X.X |
|00002b40| 20 20 20 20 2f 2f 20 77 | 65 20 77 61 6e 74 20 74 | // w|e want t|
|00002b50| 6f 20 66 69 6e 64 20 74 | 68 65 20 76 61 6c 75 65 |o find t|he value|
|00002b60| 73 20 6f 66 20 61 2c 61 | 31 2c 61 32 20 74 6f 20 |s of a,a|1,a2 to |
|00002b70| 67 69 76 65 20 74 68 65 | 20 62 65 73 74 0a 58 20 |give the| best.X |
|00002b80| 20 20 20 20 20 2f 2f 20 | 66 69 74 20 6f 66 20 79 | // |fit of y|
|00002b90| 5b 69 5d 20 77 69 74 68 | 20 61 30 20 2b 20 61 31 |[i] with| a0 + a1|
|00002ba0| 2a 78 31 5b 69 5d 20 2b | 20 61 32 2a 78 32 5b 69 |*x1[i] +| a2*x2[i|
|00002bb0| 5d 0a 58 20 20 20 20 20 | 20 2f 2f 20 41 6c 73 6f |].X | // Also|
|00002bc0| 20 70 72 69 6e 74 20 64 | 69 61 67 6f 6e 61 6c 20 | print d|iagonal |
|00002bd0| 65 6c 65 6d 65 6e 74 73 | 20 6f 66 20 68 61 74 20 |elements| of hat |
|00002be0| 6d 61 74 72 69 78 2c 20 | 58 2a 28 58 2e 74 28 29 |matrix, |X*(X.t()|
|00002bf0| 2a 58 29 2e 69 28 29 2a | 58 2e 74 28 29 0a 58 0a |*X).i()*|X.t().X.|
|00002c00| 58 20 20 20 20 20 20 2f | 2f 20 74 68 69 73 20 65 |X /|/ this e|
|00002c10| 78 61 6d 70 6c 65 20 64 | 65 6d 6f 6e 73 74 72 61 |xample d|emonstra|
|00002c20| 74 65 73 20 66 6f 75 72 | 20 6d 65 74 68 6f 64 73 |tes four| methods|
|00002c30| 20 6f 66 20 63 61 6c 63 | 75 6c 61 74 69 6f 6e 0a | of calc|ulation.|
|00002c40| 58 0a 58 20 20 20 20 20 | 20 54 72 79 0a 58 20 20 |X.X | Try.X |
|00002c50| 20 20 20 20 7b 0a 58 20 | 20 20 20 20 20 20 20 20 | {.X | |
|00002c60| 74 65 73 74 31 28 79 2c | 20 78 31 2c 20 78 32 2c |test1(y,| x1, x2,|
|00002c70| 20 6e 6f 62 73 2c 20 6e | 70 72 65 64 29 3b 0a 58 | nobs, n|pred);.X|
|00002c80| 20 20 20 20 20 20 20 20 | 20 74 65 73 74 32 28 79 | | test2(y|
|00002c90| 2c 20 78 31 2c 20 78 32 | 2c 20 6e 6f 62 73 2c 20 |, x1, x2|, nobs, |
|00002ca0| 6e 70 72 65 64 29 3b 0a | 58 20 20 20 20 20 20 20 |npred);.|X |
|00002cb0| 20 20 74 65 73 74 33 28 | 79 2c 20 78 31 2c 20 78 | test3(|y, x1, x|
|00002cc0| 32 2c 20 6e 6f 62 73 2c | 20 6e 70 72 65 64 29 3b |2, nobs,| npred);|
|00002cd0| 0a 58 20 20 20 20 20 20 | 20 20 20 74 65 73 74 34 |.X | test4|
|00002ce0| 28 79 2c 20 78 31 2c 20 | 78 32 2c 20 6e 6f 62 73 |(y, x1, |x2, nobs|
|00002cf0| 2c 20 6e 70 72 65 64 29 | 3b 0a 58 20 20 20 20 20 |, npred)|;.X |
|00002d00| 20 7d 0a 58 20 20 20 20 | 20 20 43 61 74 63 68 28 | }.X | Catch(|
|00002d10| 44 61 74 61 45 78 63 65 | 70 74 69 6f 6e 29 20 7b |DataExce|ption) {|
|00002d20| 20 63 6f 75 74 20 3c 3c | 20 22 5c 6e 49 6e 76 61 | cout <<| "\nInva|
|00002d30| 6c 69 64 20 64 61 74 61 | 5c 6e 22 3b 20 7d 0a 58 |lid data|\n"; }.X|
|00002d40| 20 20 20 20 20 20 43 61 | 74 63 68 28 53 70 61 63 | Ca|tch(Spac|
|00002d50| 65 45 78 63 65 70 74 69 | 6f 6e 29 20 7b 20 63 6f |eExcepti|on) { co|
|00002d60| 75 74 20 3c 3c 20 22 5c | 6e 4d 65 6d 6f 72 79 20 |ut << "\|nMemory |
|00002d70| 65 78 68 61 75 73 74 65 | 64 5c 6e 22 3b 20 7d 0a |exhauste|d\n"; }.|
|00002d80| 58 20 20 20 20 20 20 43 | 61 74 63 68 41 6c 6c 20 |X C|atchAll |
|00002d90| 7b 20 63 6f 75 74 20 3c | 3c 20 22 5c 6e 55 6e 65 |{ cout <|< "\nUne|
|00002da0| 78 70 65 63 74 65 64 20 | 70 72 6f 67 72 61 6d 20 |xpected |program |
|00002db0| 66 61 69 6c 75 72 65 5c | 6e 22 3b 20 7d 0a 58 20 |failure\|n"; }.X |
|00002dc0| 20 20 7d 0a 58 0a 58 23 | 69 66 64 65 66 20 44 4f | }.X.X#|ifdef DO|
|00002dd0| 5f 46 52 45 45 5f 43 48 | 45 43 4b 0a 58 20 20 20 |_FREE_CH|ECK.X |
|00002de0| 46 72 65 65 43 68 65 63 | 6b 3a 3a 53 74 61 74 75 |FreeChec|k::Statu|
|00002df0| 73 28 29 3b 0a 58 23 65 | 6e 64 69 66 0a 58 20 20 |s();.X#e|ndif.X |
|00002e00| 20 52 65 61 6c 2a 20 73 | 32 3b 20 7b 20 43 6f 6c | Real* s|2; { Col|
|00002e10| 75 6d 6e 56 65 63 74 6f | 72 20 41 28 38 30 30 30 |umnVecto|r A(8000|
|00002e20| 29 3b 20 73 32 20 3d 20 | 41 2e 53 74 6f 72 65 28 |); s2 = |A.Store(|
|00002e30| 29 3b 20 7d 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c |); }.X | cout <<|
|00002e40| 20 22 5c 6e 5c 6e 43 68 | 65 63 6b 69 6e 67 20 66 | "\n\nCh|ecking f|
|00002e50| 6f 72 20 6c 6f 73 74 20 | 6d 65 6d 6f 72 79 3a 20 |or lost |memory: |
|00002e60| 22 0a 58 20 20 20 20 20 | 20 3c 3c 20 28 75 6e 73 |".X | << (uns|
|00002e70| 69 67 6e 65 64 20 6c 6f | 6e 67 29 73 31 20 3c 3c |igned lo|ng)s1 <<|
|00002e80| 20 22 20 22 20 3c 3c 20 | 28 75 6e 73 69 67 6e 65 | " " << |(unsigne|
|00002e90| 64 20 6c 6f 6e 67 29 73 | 32 20 3c 3c 20 22 20 22 |d long)s|2 << " "|
|00002ea0| 3b 0a 58 20 20 20 69 66 | 20 28 73 31 20 21 3d 20 |;.X if| (s1 != |
|00002eb0| 73 32 29 20 63 6f 75 74 | 20 3c 3c 20 22 20 2d 20 |s2) cout| << " - |
|00002ec0| 65 72 72 6f 72 5c 6e 22 | 3b 20 65 6c 73 65 20 63 |error\n"|; else c|
|00002ed0| 6f 75 74 20 3c 3c 20 22 | 20 2d 20 6f 6b 5c 6e 22 |out << "| - ok\n"|
|00002ee0| 3b 0a 58 0a 58 20 20 20 | 72 65 74 75 72 6e 20 30 |;.X.X |return 0|
|00002ef0| 3b 0a 58 0a 58 7d 0a 58 | 0a 58 52 65 61 6c 20 74 |;.X.X}.X|.XReal t|
|00002f00| 33 28 52 65 61 6c 20 72 | 29 20 7b 20 72 65 74 75 |3(Real r|) { retu|
|00002f10| 72 6e 20 69 6e 74 28 72 | 2a 31 30 30 30 29 20 2f |rn int(r|*1000) /|
|00002f20| 20 31 30 30 30 2e 30 3b | 20 7d 0a 45 4e 44 5f 4f | 1000.0;| }.END_O|
|00002f30| 46 5f 46 49 4c 45 0a 69 | 66 20 74 65 73 74 20 39 |F_FILE.i|f test 9|
|00002f40| 35 34 33 20 2d 6e 65 20 | 60 77 63 20 2d 63 20 3c |543 -ne |`wc -c <|
|00002f50| 27 65 78 61 6d 70 6c 65 | 2e 63 78 78 27 60 3b 20 |'example|.cxx'`; |
|00002f60| 74 68 65 6e 0a 20 20 20 | 20 65 63 68 6f 20 73 68 |then. | echo sh|
|00002f70| 61 72 3a 20 5c 22 27 65 | 78 61 6d 70 6c 65 2e 63 |ar: \"'e|xample.c|
|00002f80| 78 78 27 5c 22 20 75 6e | 70 61 63 6b 65 64 20 77 |xx'\" un|packed w|
|00002f90| 69 74 68 20 77 72 6f 6e | 67 20 73 69 7a 65 21 0a |ith wron|g size!.|
|00002fa0| 66 69 0a 23 20 65 6e 64 | 20 6f 66 20 27 65 78 61 |fi.# end| of 'exa|
|00002fb0| 6d 70 6c 65 2e 63 78 78 | 27 0a 66 69 0a 69 66 20 |mple.cxx|'.fi.if |
|00002fc0| 74 65 73 74 20 2d 66 20 | 27 6e 65 77 6d 61 74 61 |test -f |'newmata|
|00002fd0| 2e 74 78 74 27 20 2d 61 | 20 22 24 7b 31 7d 22 20 |.txt' -a| "${1}" |
|00002fe0| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|00002ff0| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c | echo s|har: Wil|
|00003000| 6c 20 6e 6f 74 20 63 6c | 6f 62 62 65 72 20 65 78 |l not cl|obber ex|
|00003010| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 27 6e |isting f|ile \"'n|
|00003020| 65 77 6d 61 74 61 2e 74 | 78 74 27 5c 22 0a 65 6c |ewmata.t|xt'\".el|
|00003030| 73 65 0a 65 63 68 6f 20 | 73 68 61 72 3a 20 45 78 |se.echo |shar: Ex|
|00003040| 74 72 61 63 74 69 6e 67 | 20 5c 22 27 6e 65 77 6d |tracting| \"'newm|
|00003050| 61 74 61 2e 74 78 74 27 | 5c 22 20 5c 28 34 38 31 |ata.txt'|\" \(481|
|00003060| 30 33 20 63 68 61 72 61 | 63 74 65 72 73 5c 29 0a |03 chara|cters\).|
|00003070| 73 65 64 20 22 73 2f 5e | 58 2f 2f 22 20 3e 27 6e |sed "s/^|X//" >'n|
|00003080| 65 77 6d 61 74 61 2e 74 | 78 74 27 20 3c 3c 27 45 |ewmata.t|xt' <<'E|
|00003090| 4e 44 5f 4f 46 5f 46 49 | 4c 45 27 0a 58 2f 2f 24 |ND_OF_FI|LE'.X//$|
|000030a0| 24 20 6e 65 77 6d 61 74 | 61 2e 74 78 74 20 20 20 |$ newmat|a.txt |
|000030b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000030c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000030d0| 44 6f 63 75 6d 65 6e 74 | 61 74 69 6f 6e 20 66 69 |Document|ation fi|
|000030e0| 6c 65 0a 58 0a 58 0a 58 | 20 20 20 44 6f 63 75 6d |le.X.X.X| Docum|
|000030f0| 65 6e 74 61 74 69 6f 6e | 20 66 6f 72 20 6e 65 77 |entation| for new|
|00003100| 6d 61 74 30 37 2c 20 61 | 6e 20 65 78 70 65 72 69 |mat07, a|n experi|
|00003110| 6d 65 6e 74 61 6c 20 6d | 61 74 72 69 78 20 70 61 |mental m|atrix pa|
|00003120| 63 6b 61 67 65 20 69 6e | 20 43 2b 2b 2e 0a 58 20 |ckage in| C++..X |
|00003130| 20 20 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d | ======|========|
|00003140| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00003150| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00003160| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00003170| 3d 3d 3d 3d 0a 58 0a 58 | 0a 58 4d 41 54 52 49 58 |====.X.X|.XMATRIX|
|00003180| 20 50 41 43 4b 41 47 45 | 20 20 20 20 20 20 20 20 | PACKAGE| |
|00003190| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000031a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000031b0| 20 20 31 20 4a 61 6e 75 | 61 72 79 2c 20 31 39 39 | 1 Janu|ary, 199|
|000031c0| 33 0a 58 0a 58 43 6f 70 | 79 72 69 67 68 74 20 28 |3.X.XCop|yright (|
|000031d0| 43 29 20 31 39 39 31 2c | 32 2c 33 3a 20 52 20 42 |C) 1991,|2,3: R B|
|000031e0| 20 44 61 76 69 65 73 0a | 58 0a 58 50 65 72 6d 69 | Davies.|X.XPermi|
|000031f0| 73 73 69 6f 6e 20 69 73 | 20 67 72 61 6e 74 65 64 |ssion is| granted|
|00003200| 20 74 6f 20 75 73 65 20 | 62 75 74 20 6e 6f 74 20 | to use |but not |
|00003210| 74 6f 20 73 65 6c 6c 2e | 0a 58 0a 58 0a 58 43 6f |to sell.|.X.X.XCo|
|00003220| 6e 74 65 6e 74 73 0a 58 | 3d 3d 3d 3d 3d 3d 3d 3d |ntents.X|========|
|00003230| 0a 58 0a 58 47 65 6e 65 | 72 61 6c 20 64 65 73 63 |.X.XGene|ral desc|
|00003240| 72 69 70 74 69 6f 6e 0a | 58 49 73 20 74 68 69 73 |ription.|XIs this|
|00003250| 20 74 68 65 20 70 61 63 | 6b 61 67 65 20 79 6f 75 | the pac|kage you|
|00003260| 20 6e 65 65 64 3f 0a 58 | 43 68 61 6e 67 65 73 0a | need?.X|Changes.|
|00003270| 58 57 68 65 72 65 20 79 | 6f 75 20 63 61 6e 20 67 |XWhere y|ou can g|
|00003280| 65 74 20 61 20 63 6f 70 | 79 20 6f 66 20 74 68 69 |et a cop|y of thi|
|00003290| 73 20 70 61 63 6b 61 67 | 65 0a 58 43 6f 6d 70 69 |s packag|e.XCompi|
|000032a0| 6c 65 72 20 70 65 72 66 | 6f 72 6d 61 6e 63 65 0a |ler perf|ormance.|
|000032b0| 58 45 78 61 6d 70 6c 65 | 0a 58 44 65 74 61 69 6c |XExample|.XDetail|
|000032c0| 65 64 20 64 6f 63 75 6d | 65 6e 74 61 74 69 6f 6e |ed docum|entation|
|000032d0| 0a 58 20 20 20 43 75 73 | 74 6f 6d 69 73 69 6e 67 |.X Cus|tomising|
|000032e0| 0a 58 20 20 20 43 6f 6e | 73 74 72 75 63 74 6f 72 |.X Con|structor|
|000032f0| 73 0a 58 20 20 20 45 6c | 65 6d 65 6e 74 73 20 6f |s.X El|ements o|
|00003300| 66 20 6d 61 74 72 69 63 | 65 73 0a 58 20 20 20 4d |f matric|es.X M|
|00003310| 61 74 72 69 78 20 63 6f | 70 79 0a 58 20 20 20 45 |atrix co|py.X E|
|00003320| 6e 74 65 72 69 6e 67 20 | 76 61 6c 75 65 73 0a 58 |ntering |values.X|
|00003330| 20 20 20 55 6e 61 72 79 | 20 6f 70 65 72 61 74 6f | Unary| operato|
|00003340| 72 73 0a 58 20 20 20 42 | 69 6e 61 72 79 20 6f 70 |rs.X B|inary op|
|00003350| 65 72 61 74 6f 72 73 0a | 58 20 20 20 43 6f 6d 62 |erators.|X Comb|
|00003360| 69 6e 61 74 69 6f 6e 20 | 6f 66 20 61 20 6d 61 74 |ination |of a mat|
|00003370| 72 69 78 20 61 6e 64 20 | 73 63 61 6c 61 72 0a 58 |rix and |scalar.X|
|00003380| 20 20 20 53 63 61 6c 61 | 72 20 66 75 6e 63 74 69 | Scala|r functi|
|00003390| 6f 6e 73 20 6f 66 20 6d | 61 74 72 69 63 65 73 0a |ons of m|atrices.|
|000033a0| 58 20 20 20 53 75 62 6d | 61 74 72 69 78 20 6f 70 |X Subm|atrix op|
|000033b0| 65 72 61 74 69 6f 6e 73 | 0a 58 20 20 20 43 68 61 |erations|.X Cha|
|000033c0| 6e 67 65 20 64 69 6d 65 | 6e 73 69 6f 6e 73 0a 58 |nge dime|nsions.X|
|000033d0| 20 20 20 43 68 61 6e 67 | 65 20 74 79 70 65 0a 58 | Chang|e type.X|
|000033e0| 20 20 20 4d 75 6c 74 69 | 70 6c 65 20 6d 61 74 72 | Multi|ple matr|
|000033f0| 69 78 20 73 6f 6c 76 65 | 0a 58 20 20 20 4d 65 6d |ix solve|.X Mem|
|00003400| 6f 72 79 20 6d 61 6e 61 | 67 65 6d 65 6e 74 0a 58 |ory mana|gement.X|
|00003410| 20 20 20 45 66 66 69 63 | 69 65 6e 63 79 0a 58 20 | Effic|iency.X |
|00003420| 20 20 4f 75 74 70 75 74 | 0a 58 20 20 20 41 63 63 | Output|.X Acc|
|00003430| 65 73 73 69 6e 67 20 6d | 61 74 72 69 63 65 73 20 |essing m|atrices |
|00003440| 6f 66 20 75 6e 73 70 65 | 63 69 66 69 65 64 20 74 |of unspe|cified t|
|00003450| 79 70 65 0a 58 20 20 20 | 43 68 6f 6c 65 73 6b 79 |ype.X |Cholesky|
|00003460| 20 64 65 63 6f 6d 70 6f | 73 69 74 69 6f 6e 0a 58 | decompo|sition.X|
|00003470| 20 20 20 48 6f 75 73 65 | 68 6f 6c 64 65 72 20 74 | House|holder t|
|00003480| 72 69 61 6e 67 75 6c 61 | 72 69 73 61 74 69 6f 6e |riangula|risation|
|00003490| 0a 58 20 20 20 53 69 6e | 67 75 6c 61 72 20 56 61 |.X Sin|gular Va|
|000034a0| 6c 75 65 20 44 65 63 6f | 6d 70 6f 73 69 74 69 6f |lue Deco|mpositio|
|000034b0| 6e 0a 58 20 20 20 45 69 | 67 65 6e 76 61 6c 75 65 |n.X Ei|genvalue|
|000034c0| 73 0a 58 20 20 20 53 6f | 72 74 69 6e 67 0a 58 20 |s.X So|rting.X |
|000034d0| 20 20 46 61 73 74 20 46 | 6f 75 72 69 65 72 20 54 | Fast F|ourier T|
|000034e0| 72 61 6e 73 66 6f 72 6d | 0a 58 20 20 20 49 6e 74 |ransform|.X Int|
|000034f0| 65 72 66 61 63 65 20 74 | 6f 20 4e 75 6d 65 72 69 |erface t|o Numeri|
|00003500| 63 61 6c 20 52 65 63 69 | 70 65 73 20 69 6e 20 43 |cal Reci|pes in C|
|00003510| 0a 58 20 20 20 45 78 63 | 65 70 74 69 6f 6e 73 0a |.X Exc|eptions.|
|00003520| 58 20 20 20 43 6c 65 61 | 6e 20 75 70 20 66 6f 6c |X Clea|n up fol|
|00003530| 6c 6f 77 69 6e 67 20 61 | 6e 20 65 78 63 65 70 74 |lowing a|n except|
|00003540| 69 6f 6e 0a 58 4c 69 73 | 74 20 6f 66 20 66 69 6c |ion.XLis|t of fil|
|00003550| 65 73 0a 58 50 72 6f 62 | 6c 65 6d 20 72 65 70 6f |es.XProb|lem repo|
|00003560| 72 74 20 66 6f 72 6d 0a | 58 0a 58 0a 58 2d 2d 2d |rt form.|X.X.X---|
|00003570| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003580| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003590| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000035a0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000035b0| 2d 2d 2d 2d 2d 2d 2d 2d | 0a 58 0a 58 0a 58 47 65 |--------|.X.X.XGe|
|000035c0| 6e 65 72 61 6c 20 64 65 | 73 63 72 69 70 74 69 6f |neral de|scriptio|
|000035d0| 6e 0a 58 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |n.X=====|========|
|000035e0| 3d 3d 3d 3d 3d 3d 0a 58 | 0a 58 54 68 65 20 70 61 |======.X|.XThe pa|
|000035f0| 63 6b 61 67 65 20 69 73 | 20 69 6e 74 65 6e 74 65 |ckage is| intente|
|00003600| 64 20 66 6f 72 20 73 63 | 69 65 6e 74 69 73 74 73 |d for sc|ientists|
|00003610| 20 61 6e 64 20 65 6e 67 | 69 6e 65 65 72 73 20 77 | and eng|ineers w|
|00003620| 68 6f 20 6e 65 65 64 20 | 74 6f 0a 58 6d 61 6e 69 |ho need |to.Xmani|
|00003630| 70 75 6c 61 74 65 20 61 | 20 76 61 72 69 65 74 79 |pulate a| variety|
|00003640| 20 6f 66 20 74 79 70 65 | 73 20 6f 66 20 6d 61 74 | of type|s of mat|
|00003650| 72 69 63 65 73 20 75 73 | 69 6e 67 20 73 74 61 6e |rices us|ing stan|
|00003660| 64 61 72 64 20 6d 61 74 | 72 69 78 0a 58 6f 70 65 |dard mat|rix.Xope|
|00003670| 72 61 74 69 6f 6e 73 2e | 20 45 6d 70 68 61 73 69 |rations.| Emphasi|
|00003680| 73 20 69 73 20 6f 6e 20 | 74 68 65 20 6b 69 6e 64 |s is on |the kind|
|00003690| 20 6f 66 20 6f 70 65 72 | 61 74 69 6f 6e 73 20 6e | of oper|ations n|
|000036a0| 65 65 64 65 64 20 69 6e | 20 73 74 61 74 69 73 74 |eeded in| statist|
|000036b0| 69 63 61 6c 0a 58 63 61 | 6c 63 75 6c 61 74 69 6f |ical.Xca|lculatio|
|000036c0| 6e 73 20 73 75 63 68 20 | 61 73 20 6c 65 61 73 74 |ns such |as least|
|000036d0| 20 73 71 75 61 72 65 73 | 2c 20 6c 69 6e 65 61 72 | squares|, linear|
|000036e0| 20 65 71 75 61 74 69 6f | 6e 20 73 6f 6c 76 65 20 | equatio|n solve |
|000036f0| 61 6e 64 0a 58 65 69 67 | 65 6e 76 61 6c 75 65 73 |and.Xeig|envalues|
|00003700| 2e 0a 58 0a 58 49 74 20 | 73 75 70 70 6f 72 74 73 |..X.XIt |supports|
|00003710| 20 6d 61 74 72 69 78 20 | 74 79 70 65 73 0a 58 0a | matrix |types.X.|
|00003720| 58 20 20 20 20 4d 61 74 | 72 69 78 20 20 20 20 20 |X Mat|rix |
|00003730| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003740| 20 20 28 72 65 63 74 61 | 6e 67 75 6c 61 72 20 6d | (recta|ngular m|
|00003750| 61 74 72 69 78 29 0a 58 | 20 20 20 20 6e 72 69 63 |atrix).X| nric|
|00003760| 4d 61 74 72 69 78 20 20 | 20 20 20 20 20 20 20 20 |Matrix | |
|00003770| 20 20 20 20 20 20 20 20 | 20 28 76 61 72 69 61 6e | | (varian|
|00003780| 74 20 6f 66 20 72 65 63 | 74 61 6e 67 75 6c 61 72 |t of rec|tangular|
|00003790| 20 6d 61 74 72 69 78 29 | 0a 58 20 20 20 20 55 70 | matrix)|.X Up|
|000037a0| 70 65 72 54 72 69 61 6e | 67 75 6c 61 72 4d 61 74 |perTrian|gularMat|
|000037b0| 72 69 78 0a 58 20 20 20 | 20 4c 6f 77 65 72 54 72 |rix.X | LowerTr|
|000037c0| 69 61 6e 67 75 6c 61 72 | 4d 61 74 72 69 78 0a 58 |iangular|Matrix.X|
|000037d0| 20 20 20 20 44 69 61 67 | 6f 6e 61 6c 4d 61 74 72 | Diag|onalMatr|
|000037e0| 69 78 0a 58 20 20 20 20 | 53 79 6d 6d 65 74 72 69 |ix.X |Symmetri|
|000037f0| 63 4d 61 74 72 69 78 0a | 58 20 20 20 20 42 61 6e |cMatrix.|X Ban|
|00003800| 64 4d 61 74 72 69 78 0a | 58 20 20 20 20 55 70 70 |dMatrix.|X Upp|
|00003810| 65 72 42 61 6e 64 4d 61 | 74 72 69 78 20 20 20 20 |erBandMa|trix |
|00003820| 20 20 20 20 20 20 20 20 | 20 20 28 75 70 70 65 72 | | (upper|
|00003830| 20 74 72 69 61 6e 67 75 | 6c 61 72 20 62 61 6e 64 | triangu|lar band|
|00003840| 20 6d 61 74 72 69 78 29 | 0a 58 20 20 20 20 4c 6f | matrix)|.X Lo|
|00003850| 77 65 72 42 61 6e 64 4d | 61 74 72 69 78 20 20 20 |werBandM|atrix |
|00003860| 20 20 20 20 20 20 20 20 | 20 20 20 28 6c 6f 77 65 | | (lowe|
|00003870| 72 20 74 72 69 61 6e 67 | 75 6c 61 72 20 62 61 6e |r triang|ular ban|
|00003880| 64 20 6d 61 74 72 69 78 | 29 0a 58 20 20 20 20 53 |d matrix|).X S|
|00003890| 79 6d 6d 65 74 72 69 63 | 42 61 6e 64 4d 61 74 72 |ymmetric|BandMatr|
|000038a0| 69 78 0a 58 20 20 20 20 | 52 6f 77 56 65 63 74 6f |ix.X |RowVecto|
|000038b0| 72 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |r | |
|000038c0| 20 20 20 20 20 28 64 65 | 72 69 76 65 64 20 66 72 | (de|rived fr|
|000038d0| 6f 6d 20 4d 61 74 72 69 | 78 29 0a 58 20 20 20 20 |om Matri|x).X |
|000038e0| 43 6f 6c 75 6d 6e 56 65 | 63 74 6f 72 20 20 20 20 |ColumnVe|ctor |
|000038f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 28 64 65 | | (de|
|00003900| 72 69 76 65 64 20 66 72 | 6f 6d 20 4d 61 74 72 69 |rived fr|om Matri|
|00003910| 78 29 2e 0a 58 0a 58 4f | 6e 6c 79 20 6f 6e 65 20 |x)..X.XO|nly one |
|00003920| 65 6c 65 6d 65 6e 74 20 | 74 79 70 65 20 28 66 6c |element |type (fl|
|00003930| 6f 61 74 20 6f 72 20 64 | 6f 75 62 6c 65 29 20 69 |oat or d|ouble) i|
|00003940| 73 20 73 75 70 70 6f 72 | 74 65 64 2e 0a 58 0a 58 |s suppor|ted..X.X|
|00003950| 54 68 65 20 70 61 63 6b | 61 67 65 20 69 6e 63 6c |The pack|age incl|
|00003960| 75 64 65 73 20 74 68 65 | 20 6f 70 65 72 61 74 69 |udes the| operati|
|00003970| 6f 6e 73 20 2a 2c 20 2b | 2c 20 2d 2c 20 69 6e 76 |ons *, +|, -, inv|
|00003980| 65 72 73 65 2c 20 74 72 | 61 6e 73 70 6f 73 65 2c |erse, tr|anspose,|
|00003990| 0a 58 63 6f 6e 76 65 72 | 73 69 6f 6e 20 62 65 74 |.Xconver|sion bet|
|000039a0| 77 65 65 6e 20 74 79 70 | 65 73 2c 20 73 75 62 6d |ween typ|es, subm|
|000039b0| 61 74 72 69 78 2c 20 64 | 65 74 65 72 6d 69 6e 61 |atrix, d|etermina|
|000039c0| 6e 74 2c 20 43 68 6f 6c | 65 73 6b 79 0a 58 64 65 |nt, Chol|esky.Xde|
|000039d0| 63 6f 6d 70 6f 73 69 74 | 69 6f 6e 2c 20 48 6f 75 |composit|ion, Hou|
|000039e0| 73 65 68 6f 6c 64 65 72 | 20 74 72 69 61 6e 67 75 |seholder| triangu|
|000039f0| 6c 61 72 69 73 61 74 69 | 6f 6e 2c 20 73 69 6e 67 |larisati|on, sing|
|00003a00| 75 6c 61 72 20 76 61 6c | 75 65 0a 58 64 65 63 6f |ular val|ue.Xdeco|
|00003a10| 6d 70 6f 73 69 74 69 6f | 6e 2c 20 65 69 67 65 6e |mpositio|n, eigen|
|00003a20| 76 61 6c 75 65 73 20 6f | 66 20 61 20 73 79 6d 6d |values o|f a symm|
|00003a30| 65 74 72 69 63 20 6d 61 | 74 72 69 78 2c 20 73 6f |etric ma|trix, so|
|00003a40| 72 74 69 6e 67 2c 20 66 | 61 73 74 20 66 6f 75 72 |rting, f|ast four|
|00003a50| 69 65 72 0a 58 74 72 61 | 6e 73 66 6f 72 6d 2c 20 |ier.Xtra|nsform, |
|00003a60| 70 72 69 6e 74 69 6e 67 | 20 61 6e 64 20 61 6e 20 |printing| and an |
|00003a70| 69 6e 74 65 72 66 61 63 | 65 20 77 69 74 68 20 22 |interfac|e with "|
|00003a80| 4e 75 6d 65 72 69 63 61 | 6c 20 52 65 63 69 70 65 |Numerica|l Recipe|
|00003a90| 73 20 69 6e 20 43 22 2e | 0a 58 0a 58 49 74 20 69 |s in C".|.X.XIt i|
|00003aa0| 73 20 69 6e 74 65 6e 64 | 65 64 20 66 6f 72 20 6d |s intend|ed for m|
|00003ab0| 61 74 72 69 63 65 73 20 | 69 6e 20 74 68 65 20 72 |atrices |in the r|
|00003ac0| 61 6e 67 65 20 31 35 20 | 78 20 31 35 20 74 6f 20 |ange 15 |x 15 to |
|00003ad0| 74 68 65 20 6d 61 78 69 | 6d 75 6d 20 73 69 7a 65 |the maxi|mum size|
|00003ae0| 0a 58 79 6f 75 72 20 6d | 61 63 68 69 6e 65 20 77 |.Xyour m|achine w|
|00003af0| 69 6c 6c 20 61 63 63 6f | 6d 6f 64 61 74 65 20 69 |ill acco|modate i|
|00003b00| 6e 20 61 20 73 69 6e 67 | 6c 65 20 61 72 72 61 79 |n a sing|le array|
|00003b10| 2e 20 46 6f 72 20 65 78 | 61 6d 70 6c 65 20 39 30 |. For ex|ample 90|
|00003b20| 20 78 20 39 30 20 28 31 | 32 35 0a 58 78 20 31 32 | x 90 (1|25.Xx 12|
|00003b30| 35 20 66 6f 72 20 74 72 | 69 61 6e 67 75 6c 61 72 |5 for tr|iangular|
|00003b40| 20 6d 61 74 72 69 63 65 | 73 29 20 69 6e 20 6d 61 | matrice|s) in ma|
|00003b50| 63 68 69 6e 65 73 20 74 | 68 61 74 20 68 61 76 65 |chines t|hat have|
|00003b60| 20 38 31 39 32 20 64 6f | 75 62 6c 65 73 20 61 73 | 8192 do|ubles as|
|00003b70| 20 74 68 65 0a 58 6d 61 | 78 69 6d 75 6d 20 73 69 | the.Xma|ximum si|
|00003b80| 7a 65 20 6f 66 20 61 6e | 20 61 72 72 61 79 2e 20 |ze of an| array. |
|00003b90| 54 68 65 20 6e 75 6d 62 | 65 72 20 6f 66 20 65 6c |The numb|er of el|
|00003ba0| 65 6d 65 6e 74 73 20 69 | 6e 20 61 6e 20 61 72 72 |ements i|n an arr|
|00003bb0| 61 79 20 63 61 6e 6e 6f | 74 0a 58 65 78 63 65 65 |ay canno|t.Xexcee|
|00003bc0| 64 20 74 68 65 20 6d 61 | 78 69 6d 75 6d 20 73 69 |d the ma|ximum si|
|00003bd0| 7a 65 20 6f 66 20 61 6e | 20 22 69 6e 74 22 2e 20 |ze of an| "int". |
|00003be0| 54 68 65 20 70 61 63 6b | 61 67 65 20 77 69 6c 6c |The pack|age will|
|00003bf0| 20 77 6f 72 6b 20 66 6f | 72 20 76 65 72 79 0a 58 | work fo|r very.X|
|00003c00| 73 6d 61 6c 6c 20 6d 61 | 74 72 69 63 65 73 20 62 |small ma|trices b|
|00003c10| 75 74 20 62 65 63 6f 6d | 65 73 20 72 61 74 68 65 |ut becom|es rathe|
|00003c20| 72 20 69 6e 65 66 66 69 | 63 69 65 6e 74 2e 0a 58 |r ineffi|cient..X|
|00003c30| 0a 58 41 20 74 77 6f 2d | 73 74 61 67 65 20 61 70 |.XA two-|stage ap|
|00003c40| 70 72 6f 61 63 68 20 74 | 6f 20 65 76 61 6c 75 61 |proach t|o evalua|
|00003c50| 74 69 6e 67 20 6d 61 74 | 72 69 78 20 65 78 70 72 |ting mat|rix expr|
|00003c60| 65 73 73 69 6f 6e 73 20 | 69 73 20 75 73 65 64 20 |essions |is used |
|00003c70| 74 6f 20 69 6d 70 72 6f | 76 65 0a 58 65 66 66 69 |to impro|ve.Xeffi|
|00003c80| 63 69 65 6e 63 79 20 61 | 6e 64 20 72 65 64 75 63 |ciency a|nd reduc|
|00003c90| 65 20 75 73 65 20 6f 66 | 20 74 65 6d 70 6f 72 61 |e use of| tempora|
|00003ca0| 72 79 20 73 74 6f 72 61 | 67 65 2e 0a 58 0a 58 54 |ry stora|ge..X.XT|
|00003cb0| 68 65 20 70 61 63 6b 61 | 67 65 20 69 73 20 64 65 |he packa|ge is de|
|00003cc0| 73 69 67 6e 65 64 20 66 | 6f 72 20 76 65 72 73 69 |signed f|or versi|
|00003cd0| 6f 6e 20 32 20 6f 72 20 | 33 20 6f 66 20 43 2b 2b |on 2 or |3 of C++|
|00003ce0| 2e 20 49 74 20 77 6f 72 | 6b 73 20 77 69 74 68 20 |. It wor|ks with |
|00003cf0| 42 6f 72 6c 61 6e 64 0a | 58 28 33 2e 31 29 20 61 |Borland.|X(3.1) a|
|00003d00| 6e 64 20 4d 69 63 72 6f | 73 6f 66 74 20 28 37 2e |nd Micro|soft (7.|
|00003d10| 30 29 20 43 2b 2b 20 6f | 6e 20 61 20 50 43 20 61 |0) C++ o|n a PC a|
|00003d20| 6e 64 20 41 54 26 54 20 | 43 2b 2b 20 28 32 2e 31 |nd AT&T |C++ (2.1|
|00003d30| 20 26 20 33 29 20 61 6e | 64 20 47 6e 75 20 43 2b | & 3) an|d Gnu C+|
|00003d40| 2b 0a 58 28 32 2e 32 29 | 2e 20 49 74 20 77 6f 72 |+.X(2.2)|. It wor|
|00003d50| 6b 73 20 77 69 74 68 20 | 73 6f 6d 65 20 70 72 6f |ks with |some pro|
|00003d60| 62 6c 65 6d 73 20 77 69 | 74 68 20 5a 6f 72 74 65 |blems wi|th Zorte|
|00003d70| 63 68 20 43 2b 2b 20 28 | 76 65 72 73 69 6f 6e 20 |ch C++ (|version |
|00003d80| 33 2e 30 34 29 2e 0a 58 | 0a 58 0a 58 2d 2d 2d 2d |3.04)..X|.X.X----|
|00003d90| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003da0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003db0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003dc0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00003dd0| 2d 2d 2d 2d 2d 2d 2d 0a | 58 0a 58 0a 58 49 73 20 |-------.|X.X.XIs |
|00003de0| 74 68 69 73 20 74 68 65 | 20 70 61 63 6b 61 67 65 |this the| package|
|00003df0| 20 79 6f 75 20 6e 65 65 | 64 3f 0a 58 3d 3d 3d 3d | you nee|d?.X====|
|00003e00| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00003e10| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 0a 58 0a 58 44 6f 20 |========|=.X.XDo |
|00003e20| 79 6f 75 0a 58 0a 58 31 | 2e 20 20 20 6e 65 65 64 |you.X.X1|. need|
|00003e30| 20 6d 61 74 72 69 78 20 | 6f 70 65 72 61 74 6f 72 | matrix |operator|
|00003e40| 73 20 73 75 63 68 20 61 | 73 20 2a 20 61 6e 64 20 |s such a|s * and |
|00003e50| 2b 20 64 65 66 69 6e 65 | 64 20 61 73 20 6f 70 65 |+ define|d as ope|
|00003e60| 72 61 74 6f 72 73 20 73 | 6f 20 79 6f 75 0a 58 20 |rators s|o you.X |
|00003e70| 20 20 20 20 63 61 6e 20 | 77 72 69 74 65 20 74 68 | can |write th|
|00003e80| 69 6e 67 73 20 6c 69 6b | 65 0a 58 0a 58 20 20 20 |ings lik|e.X.X |
|00003e90| 20 20 20 20 20 58 20 20 | 3d 20 41 20 2a 20 28 42 | X |= A * (B|
|00003ea0| 20 2b 20 43 29 3b 0a 58 | 0a 58 32 2e 20 20 20 6e | + C);.X|.X2. n|
|00003eb0| 65 65 64 20 61 20 76 61 | 72 69 65 74 79 20 6f 66 |eed a va|riety of|
|00003ec0| 20 74 79 70 65 73 20 6f | 66 20 6d 61 74 72 69 63 | types o|f matric|
|00003ed0| 65 73 0a 58 0a 58 33 2e | 20 20 20 6e 65 65 64 20 |es.X.X3.| need |
|00003ee0| 6f 6e 6c 79 20 6f 6e 65 | 20 65 6c 65 6d 65 6e 74 |only one| element|
|00003ef0| 20 74 79 70 65 20 28 66 | 6c 6f 61 74 20 6f 72 20 | type (f|loat or |
|00003f00| 64 6f 75 62 6c 65 29 0a | 58 0a 58 34 2e 20 20 20 |double).|X.X4. |
|00003f10| 77 6f 72 6b 20 77 69 74 | 68 20 6d 61 74 72 69 63 |work wit|h matric|
|00003f20| 65 73 20 69 6e 20 74 68 | 65 20 72 61 6e 67 65 20 |es in th|e range |
|00003f30| 31 30 20 78 20 31 30 20 | 75 70 20 74 6f 20 77 68 |10 x 10 |up to wh|
|00003f40| 61 74 20 63 61 6e 20 62 | 65 20 73 74 6f 72 65 64 |at can b|e stored|
|00003f50| 20 69 6e 0a 58 20 20 20 | 20 20 6f 6e 65 20 6d 65 | in.X | one me|
|00003f60| 6d 6f 72 79 20 62 6c 6f | 63 6b 0a 58 0a 58 35 2e |mory blo|ck.X.X5.|
|00003f70| 20 20 20 74 6f 6c 65 72 | 61 74 65 20 61 20 6c 61 | toler|ate a la|
|00003f80| 72 67 65 20 70 61 63 6b | 61 67 65 0a 58 0a 58 0a |rge pack|age.X.X.|
|00003f90| 58 54 68 65 6e 20 6d 61 | 79 62 65 20 74 68 69 73 |XThen ma|ybe this|
|00003fa0| 20 69 73 20 74 68 65 20 | 72 69 67 68 74 20 70 61 | is the |right pa|
|00003fb0| 63 6b 61 67 65 20 66 6f | 72 20 79 6f 75 2e 20 0a |ckage fo|r you. .|
|00003fc0| 58 0a 58 49 66 20 79 6f | 75 20 64 6f 6e 27 74 20 |X.XIf yo|u don't |
|00003fd0| 6e 65 65 64 20 28 31 29 | 20 74 68 65 6e 20 74 68 |need (1)| then th|
|00003fe0| 65 72 65 20 6d 61 79 20 | 62 65 20 62 65 74 74 65 |ere may |be bette|
|00003ff0| 72 20 6f 70 74 69 6f 6e | 73 2e 20 4c 69 6b 65 77 |r option|s. Likew|
|00004000| 69 73 65 20 69 66 20 79 | 6f 75 0a 58 64 6f 6e 27 |ise if y|ou.Xdon'|
|00004010| 74 20 6e 65 65 64 20 28 | 32 29 20 74 68 65 72 65 |t need (|2) there|
|00004020| 20 6d 61 79 20 62 65 20 | 62 65 74 74 65 72 20 6f | may be |better o|
|00004030| 70 74 69 6f 6e 73 2e 20 | 49 66 20 79 6f 75 20 72 |ptions. |If you r|
|00004040| 65 71 75 69 72 65 20 22 | 6e 6f 74 20 28 35 29 22 |equire "|not (5)"|
|00004050| 0a 58 74 68 65 6e 20 74 | 68 69 73 20 69 73 20 6e |.Xthen t|his is n|
|00004060| 6f 74 20 74 68 65 20 70 | 61 63 6b 61 67 65 20 66 |ot the p|ackage f|
|00004070| 6f 72 20 79 6f 75 2e 0a | 58 0a 58 0a 58 49 66 20 |or you..|X.X.XIf |
|00004080| 79 6f 75 20 6e 65 65 64 | 20 28 32 29 20 61 6e 64 |you need| (2) and|
|00004090| 20 22 6e 6f 74 20 28 33 | 29 22 20 61 6e 64 20 68 | "not (3|)" and h|
|000040a0| 61 76 65 20 73 6f 6d 65 | 20 73 70 61 72 65 20 6d |ave some| spare m|
|000040b0| 6f 6e 65 79 2c 20 74 68 | 65 6e 20 6d 61 79 62 65 |oney, th|en maybe|
|000040c0| 20 79 6f 75 0a 58 73 68 | 6f 75 6c 64 20 6c 6f 6f | you.Xsh|ould loo|
|000040d0| 6b 20 61 74 20 4d 2b 2b | 20 66 72 6f 6d 20 44 79 |k at M++| from Dy|
|000040e0| 61 64 20 5b 70 68 6f 6e | 65 20 69 6e 20 74 68 65 |ad [phon|e in the|
|000040f0| 20 55 53 41 20 28 38 30 | 30 29 33 36 36 2d 31 35 | USA (80|0)366-15|
|00004100| 37 33 2c 0a 58 28 32 30 | 36 29 36 33 37 2d 39 34 |73,.X(20|6)637-94|
|00004110| 32 37 2c 20 66 61 78 20 | 28 32 30 36 29 36 33 37 |27, fax |(206)637|
|00004120| 2d 39 34 32 38 5d 2c 20 | 6f 72 20 74 68 65 20 52 |-9428], |or the R|
|00004130| 6f 67 75 65 20 57 61 76 | 65 20 6d 61 74 72 69 78 |ogue Wav|e matrix|
|00004140| 20 70 61 63 6b 61 67 65 | 0a 58 5b 28 38 30 30 29 | package|.X[(800)|
|00004150| 34 38 37 2d 33 32 31 37 | 2c 20 28 35 30 33 29 37 |487-3217|, (503)7|
|00004160| 35 34 2d 33 30 31 30 2c | 20 66 61 78 20 28 35 30 |54-3010,| fax (50|
|00004170| 33 29 37 35 37 2d 36 36 | 35 30 5d 2e 0a 58 0a 58 |3)757-66|50]..X.X|
|00004180| 0a 58 49 66 20 79 6f 75 | 20 6e 65 65 64 20 6e 6f |.XIf you| need no|
|00004190| 74 20 28 34 29 3b 20 74 | 68 61 74 20 69 73 20 76 |t (4); t|hat is v|
|000041a0| 65 72 79 20 6c 61 72 67 | 65 20 6d 61 74 72 69 63 |ery larg|e matric|
|000041b0| 65 73 20 74 68 61 74 20 | 77 69 6c 6c 20 6e 65 65 |es that |will nee|
|000041c0| 64 20 74 6f 20 62 65 0a | 58 73 74 6f 72 65 64 20 |d to be.|Xstored |
|000041d0| 6f 6e 20 64 69 73 6b 2c | 20 74 68 65 72 65 20 69 |on disk,| there i|
|000041e0| 73 20 61 20 70 61 63 6b | 61 67 65 20 59 41 4d 50 |s a pack|age YAMP|
|000041f0| 20 6f 6e 20 43 6f 6d 70 | 75 53 65 72 76 65 20 75 | on Comp|uServe u|
|00004200| 6e 64 65 72 20 74 68 65 | 20 42 6f 72 6c 61 6e 64 |nder the| Borland|
|00004210| 0a 58 43 2b 2b 20 6c 69 | 62 72 61 72 79 20 74 68 |.XC++ li|brary th|
|00004220| 61 74 20 6d 69 67 68 74 | 20 62 65 20 6f 66 20 69 |at might| be of i|
|00004230| 6e 74 65 72 65 73 74 2e | 0a 58 0a 58 0a 58 44 65 |nterest.|.X.X.XDe|
|00004240| 74 61 69 6c 73 20 6f 66 | 20 73 6f 6d 65 20 6f 74 |tails of| some ot|
|00004250| 68 65 72 20 66 72 65 65 | 20 43 20 6f 72 20 43 2b |her free| C or C+|
|00004260| 2b 20 6d 61 74 72 69 78 | 20 70 61 63 6b 61 67 65 |+ matrix| package|
|00004270| 73 20 66 6f 6c 6c 6f 77 | 20 2d 20 65 78 74 72 61 |s follow| - extra|
|00004280| 63 74 65 64 0a 58 66 72 | 6f 6d 20 74 68 65 20 6c |cted.Xfr|om the l|
|00004290| 69 73 74 20 61 73 73 65 | 6d 62 6c 65 64 20 62 79 |ist asse|mbled by|
|000042a0| 20 61 6a 61 79 73 68 61 | 68 40 75 73 63 2e 65 64 | ajaysha|h@usc.ed|
|000042b0| 75 2e 0a 58 0a 58 4e 61 | 6d 65 3a 20 53 50 41 52 |u..X.XNa|me: SPAR|
|000042c0| 53 45 0a 58 57 68 65 72 | 65 3a 20 69 6e 20 73 70 |SE.XWher|e: in sp|
|000042d0| 61 72 73 65 20 6f 6e 20 | 4e 65 74 6c 69 62 0a 58 |arse on |Netlib.X|
|000042e0| 44 65 73 63 72 69 70 74 | 69 6f 6e 3a 20 6c 69 62 |Descript|ion: lib|
|000042f0| 72 61 72 79 20 66 6f 72 | 20 4c 55 20 66 61 63 74 |rary for| LU fact|
|00004300| 6f 72 69 73 61 74 69 6f | 6e 20 66 6f 72 20 6c 61 |orisatio|n for la|
|00004310| 72 67 65 20 73 70 61 72 | 73 65 20 6d 61 74 72 69 |rge spar|se matri|
|00004320| 63 65 73 20 0a 58 41 75 | 74 68 6f 72 3a 20 4b 65 |ces .XAu|thor: Ke|
|00004330| 6e 20 4b 75 6e 64 65 72 | 74 2c 20 41 6c 62 65 72 |n Kunder|t, Alber|
|00004340| 74 6f 20 53 61 6e 67 69 | 6f 76 61 6e 6e 69 2d 56 |to Sangi|ovanni-V|
|00004350| 69 6e 63 65 6e 74 65 6c | 6c 69 2c 0a 58 20 20 20 |incentel|li,.X |
|00004360| 20 20 73 70 61 72 73 65 | 40 69 63 2e 62 65 72 6b | sparse|@ic.berk|
|00004370| 65 6c 65 79 2e 65 64 75 | 0a 58 0a 58 4e 61 6d 65 |eley.edu|.X.XName|
|00004380| 3a 20 6d 61 74 72 69 78 | 2e 74 61 72 2e 5a 0a 58 |: matrix|.tar.Z.X|
|00004390| 57 68 65 72 65 3a 20 69 | 6e 20 66 74 70 2d 72 61 |Where: i|n ftp-ra|
|000043a0| 69 6d 75 6e 64 2f 70 75 | 62 2f 73 72 63 2f 4d 61 |imund/pu|b/src/Ma|
|000043b0| 74 68 20 6f 6e 20 6e 65 | 73 74 72 6f 79 2e 77 75 |th on ne|stroy.wu|
|000043c0| 2d 77 69 65 6e 2e 61 63 | 2e 61 74 0a 58 20 20 20 |-wien.ac|.at.X |
|000043d0| 20 20 28 31 33 37 2e 32 | 30 38 2e 33 2e 34 29 0a | (137.2|08.3.4).|
|000043e0| 58 41 75 74 68 6f 72 3a | 20 50 61 75 6c 20 53 63 |XAuthor:| Paul Sc|
|000043f0| 68 6d 69 64 74 2c 20 54 | 49 0a 58 44 65 73 63 72 |hmidt, T|I.XDescr|
|00004400| 69 70 74 69 6f 6e 3a 20 | 53 6d 61 6c 6c 20 6d 61 |iption: |Small ma|
|00004410| 74 72 69 78 20 6c 69 62 | 72 61 72 79 2c 20 69 6e |trix lib|rary, in|
|00004420| 63 6c 75 64 69 6e 67 20 | 53 4f 52 2c 20 57 4c 53 |cluding |SOR, WLS|
|00004430| 0a 58 0a 58 4e 61 6d 65 | 3a 20 6d 61 74 72 69 78 |.X.XName|: matrix|
|00004440| 30 34 2e 7a 69 70 0a 58 | 57 68 65 72 65 3a 20 69 |04.zip.X|Where: i|
|00004450| 6e 20 6d 69 72 72 6f 72 | 73 2f 6d 73 64 6f 73 2f |n mirror|s/msdos/|
|00004460| 63 20 6f 6e 20 77 75 61 | 72 63 68 69 76 65 2e 77 |c on wua|rchive.w|
|00004470| 75 73 74 6c 2e 65 64 75 | 0a 58 44 65 73 63 72 69 |ustl.edu|.XDescri|
|00004480| 70 74 69 6f 6e 3a 20 53 | 6d 61 6c 6c 20 6d 61 74 |ption: S|mall mat|
|00004490| 72 69 78 20 74 6f 6f 6c | 62 6f 78 0a 58 0a 58 4e |rix tool|box.X.XN|
|000044a0| 61 6d 65 3a 20 4d 61 74 | 72 69 78 2e 74 61 72 2e |ame: Mat|rix.tar.|
|000044b0| 5a 0a 58 57 68 65 72 65 | 3a 20 69 6e 20 70 75 62 |Z.XWhere|: in pub|
|000044c0| 20 66 74 70 2e 63 73 2e | 75 63 6c 61 2e 65 64 75 | ftp.cs.|ucla.edu|
|000044d0| 0a 58 44 65 73 63 72 69 | 70 74 69 6f 6e 3a 20 54 |.XDescri|ption: T|
|000044e0| 68 65 20 43 2b 2b 20 4d | 61 74 72 69 78 20 63 6c |he C++ M|atrix cl|
|000044f0| 61 73 73 2c 20 69 6e 63 | 6c 75 64 69 6e 67 20 61 |ass, inc|luding a|
|00004500| 20 6d 61 74 72 69 78 20 | 69 6d 70 6c 65 6d 65 6e | matrix |implemen|
|00004510| 74 61 74 69 6f 6e 20 6f | 66 0a 58 20 20 20 20 20 |tation o|f.X |
|00004520| 74 68 65 20 62 61 63 6b | 77 61 72 64 20 65 72 72 |the back|ward err|
|00004530| 6f 72 20 70 72 6f 70 61 | 67 61 74 69 6f 6e 20 28 |or propa|gation (|
|00004540| 62 61 63 6b 70 72 6f 70 | 29 20 61 6c 67 6f 72 69 |backprop|) algori|
|00004550| 74 68 6d 20 66 6f 72 20 | 74 72 61 69 6e 69 6e 67 |thm for |training|
|00004560| 0a 58 20 20 20 20 20 6d | 75 6c 74 69 2d 6c 61 79 |.X m|ulti-lay|
|00004570| 65 72 2c 20 66 65 65 64 | 2d 66 6f 72 77 61 72 64 |er, feed|-forward|
|00004580| 20 61 72 74 69 66 69 63 | 69 61 6c 20 6e 65 75 72 | artific|ial neur|
|00004590| 61 6c 20 6e 65 74 77 6f | 72 6b 73 0a 58 41 75 74 |al netwo|rks.XAut|
|000045a0| 68 6f 72 3a 20 45 2e 20 | 52 6f 62 65 72 74 20 28 |hor: E. |Robert (|
|000045b0| 42 6f 62 29 20 54 69 73 | 64 61 6c 65 2c 20 65 64 |Bob) Tis|dale, ed|
|000045c0| 77 69 6e 40 63 73 2e 75 | 63 6c 61 2e 65 64 75 0a |win@cs.u|cla.edu.|
|000045d0| 58 0a 58 4e 61 6d 65 3a | 20 6d 65 73 63 68 61 63 |X.XName:| meschac|
|000045e0| 68 0a 58 57 68 65 72 65 | 3a 20 69 6e 20 63 2f 6d |h.XWhere|: in c/m|
|000045f0| 65 73 63 68 61 63 68 20 | 6f 6e 20 6e 65 74 6c 69 |eschach |on netli|
|00004600| 62 0a 58 53 79 73 74 65 | 6d 73 3a 20 55 6e 69 78 |b.XSyste|ms: Unix|
|00004610| 2c 20 50 43 0a 58 44 65 | 73 63 72 69 70 74 69 6f |, PC.XDe|scriptio|
|00004620| 6e 3a 20 61 20 6c 69 62 | 72 61 72 79 20 66 6f 72 |n: a lib|rary for|
|00004630| 20 6d 61 74 72 69 78 20 | 63 6f 6d 70 75 74 61 74 | matrix |computat|
|00004640| 69 6f 6e 3b 20 6d 6f 72 | 65 20 66 75 6e 63 74 69 |ion; mor|e functi|
|00004650| 6f 6e 61 6c 69 74 79 20 | 74 68 61 6e 0a 58 20 20 |onality |than.X |
|00004660| 20 20 20 4c 69 6e 70 61 | 63 6b 3b 20 6e 6f 6e 73 | Linpa|ck; nons|
|00004670| 74 61 6e 64 61 72 64 20 | 6d 61 74 72 69 63 65 73 |tandard |matrices|
|00004680| 0a 58 41 75 74 68 6f 72 | 3a 20 44 61 76 69 64 20 |.XAuthor|: David |
|00004690| 45 2e 20 53 74 65 77 61 | 72 74 2c 20 64 65 73 40 |E. Stewa|rt, des@|
|000046a0| 74 68 72 61 69 6e 2e 61 | 6e 75 2e 65 64 75 2e 61 |thrain.a|nu.edu.a|
|000046b0| 75 0a 58 56 65 72 73 69 | 6f 6e 3a 20 31 2e 30 2c |u.XVersi|on: 1.0,|
|000046c0| 20 46 65 62 20 31 39 39 | 32 0a 58 0a 58 4e 61 6d | Feb 199|2.X.XNam|
|000046d0| 65 3a 20 6e 6c 6d 64 6c | 0a 58 57 68 65 72 65 3a |e: nlmdl|.XWhere:|
|000046e0| 20 69 6e 20 70 75 62 2f | 61 72 67 2f 6e 6c 6d 64 | in pub/|arg/nlmd|
|000046f0| 6c 20 61 74 20 63 63 76 | 72 31 2e 63 63 2e 6e 63 |l at ccv|r1.cc.nc|
|00004700| 73 75 2e 65 64 75 20 28 | 31 32 38 2e 31 30 39 2e |su.edu (|128.109.|
|00004710| 32 31 32 2e 32 30 29 0a | 58 4c 61 6e 67 75 61 67 |212.20).|XLanguag|
|00004720| 65 3a 20 43 2b 2b 0a 58 | 53 79 73 74 65 6d 73 3a |e: C++.X|Systems:|
|00004730| 20 55 6e 69 78 2c 20 4d | 53 2d 44 4f 53 20 28 54 | Unix, M|S-DOS (T|
|00004740| 75 72 62 6f 20 43 2b 2b | 29 0a 58 44 65 73 63 72 |urbo C++|).XDescr|
|00004750| 69 70 74 69 6f 6e 3a 20 | 61 20 6c 69 62 72 61 72 |iption: |a librar|
|00004760| 79 20 66 6f 72 20 65 73 | 74 69 6d 61 74 69 6f 6e |y for es|timation|
|00004770| 20 6f 66 20 6e 6f 6e 6c | 69 6e 65 61 72 20 6d 6f | of nonl|inear mo|
|00004780| 64 65 6c 73 0a 58 41 75 | 74 68 6f 72 3a 20 41 2e |dels.XAu|thor: A.|
|00004790| 20 52 6f 6e 61 6c 64 20 | 47 61 6c 6c 61 6e 74 2c | Ronald |Gallant,|
|000047a0| 20 61 72 67 40 63 63 76 | 72 31 2e 63 63 2e 6e 63 | arg@ccv|r1.cc.nc|
|000047b0| 73 75 2e 65 64 75 0a 58 | 43 6f 6d 6d 65 6e 74 73 |su.edu.X|Comments|
|000047c0| 3a 20 6e 6f 6e 6c 69 6e | 65 61 72 20 6d 61 78 69 |: nonlin|ear maxi|
|000047d0| 6d 69 73 61 74 69 6f 6e | 2c 20 65 73 74 69 6d 61 |misation|, estima|
|000047e0| 74 69 6f 6e 2c 20 69 6e | 63 6c 75 64 65 73 20 61 |tion, in|cludes a|
|000047f0| 20 72 65 61 6c 20 6d 61 | 74 72 69 78 0a 58 20 20 | real ma|trix.X |
|00004800| 20 20 20 63 6c 61 73 73 | 0a 58 56 65 72 73 69 6f | class|.XVersio|
|00004810| 6e 3a 20 4a 61 6e 75 61 | 72 79 20 31 39 39 31 0a |n: Janua|ry 1991.|
|00004820| 58 0a 58 0a 58 0a 58 2d | 2d 2d 2d 2d 2d 2d 2d 2d |X.X.X.X-|--------|
|00004830| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004840| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004850| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004860| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004870| 2d 2d 0a 58 0a 58 0a 58 | 43 68 61 6e 67 65 73 0a |--.X.X.X|Changes.|
|00004880| 58 3d 3d 3d 3d 3d 3d 3d | 0a 58 0a 58 4e 65 77 6d |X=======|.X.XNewm|
|00004890| 61 74 30 37 20 2d 20 4a | 61 6e 75 61 72 79 2c 20 |at07 - J|anuary, |
|000048a0| 31 39 39 33 0a 58 0a 58 | 4d 69 6e 6f 72 20 63 6f |1993.X.X|Minor co|
|000048b0| 72 72 65 63 74 69 6f 6e | 73 20 74 6f 20 69 6d 70 |rrection|s to imp|
|000048c0| 72 6f 76 65 20 63 6f 6d | 70 61 74 69 62 69 6c 69 |rove com|patibili|
|000048d0| 74 79 20 77 69 74 68 20 | 5a 6f 72 74 65 63 68 2c |ty with |Zortech,|
|000048e0| 20 4d 69 63 72 6f 73 6f | 66 74 20 61 6e 64 0a 58 | Microso|ft and.X|
|000048f0| 47 6e 75 2e 20 43 6f 72 | 72 65 63 74 69 6f 6e 20 |Gnu. Cor|rection |
|00004900| 74 6f 20 65 78 63 65 70 | 74 69 6f 6e 20 6d 6f 64 |to excep|tion mod|
|00004910| 75 6c 65 2e 20 41 64 64 | 69 74 69 6f 6e 61 6c 20 |ule. Add|itional |
|00004920| 46 46 54 20 66 75 6e 63 | 74 69 6f 6e 73 2e 20 53 |FFT func|tions. S|
|00004930| 6f 6d 65 0a 58 6d 69 6e | 6f 72 20 69 6e 63 72 65 |ome.Xmin|or incre|
|00004940| 61 73 65 73 20 69 6e 20 | 65 66 66 69 63 69 65 6e |ases in |efficien|
|00004950| 63 79 2e 20 53 75 62 6d | 61 74 72 69 63 65 73 20 |cy. Subm|atrices |
|00004960| 63 61 6e 20 6e 6f 77 20 | 62 65 20 75 73 65 64 20 |can now |be used |
|00004970| 6f 6e 20 52 48 53 20 6f | 66 20 3d 2e 0a 58 4f 70 |on RHS o|f =..XOp|
|00004980| 74 69 6f 6e 20 66 6f 72 | 20 61 6c 6c 6f 77 69 6e |tion for| allowin|
|00004990| 67 20 43 20 74 79 70 65 | 20 73 75 62 73 63 72 69 |g C type| subscri|
|000049a0| 70 74 73 2e 20 4d 65 74 | 68 6f 64 20 66 6f 72 20 |pts. Met|hod for |
|000049b0| 6c 6f 61 64 69 6e 67 20 | 73 68 6f 72 74 20 6c 69 |loading |short li|
|000049c0| 73 74 73 20 6f 66 0a 58 | 6e 75 6d 62 65 72 73 2e |sts of.X|numbers.|
|000049d0| 0a 58 0a 58 0a 58 4e 65 | 77 6d 61 74 30 36 20 2d |.X.X.XNe|wmat06 -|
|000049e0| 20 44 65 63 65 6d 62 65 | 72 20 31 39 39 32 3a 0a | Decembe|r 1992:.|
|000049f0| 58 0a 58 41 64 64 65 64 | 20 62 61 6e 64 20 6d 61 |X.XAdded| band ma|
|00004a00| 74 72 69 63 65 73 3b 20 | 27 72 65 61 6c 27 20 63 |trices; |'real' c|
|00004a10| 68 61 6e 67 65 64 20 74 | 6f 20 27 52 65 61 6c 27 |hanged t|o 'Real'|
|00004a20| 20 28 74 6f 20 61 76 6f | 69 64 20 70 6f 74 65 6e | (to avo|id poten|
|00004a30| 74 69 61 6c 0a 58 63 6f | 6e 66 6c 69 63 74 20 69 |tial.Xco|nflict i|
|00004a40| 6e 20 63 6f 6d 70 6c 65 | 78 20 63 6c 61 73 73 29 |n comple|x class)|
|00004a50| 3b 20 49 6e 6a 65 63 74 | 20 64 6f 65 73 6e 27 74 |; Inject| doesn't|
|00004a60| 20 63 68 65 63 6b 20 66 | 6f 72 20 6e 6f 20 6c 6f | check f|or no lo|
|00004a70| 73 73 20 6f 66 0a 58 69 | 6e 66 6f 72 6d 61 74 69 |ss of.Xi|nformati|
|00004a80| 6f 6e 3b 20 20 66 69 78 | 65 73 20 66 6f 72 20 41 |on; fix|es for A|
|00004a90| 54 26 54 20 43 2b 2b 20 | 76 65 72 73 69 6f 6e 20 |T&T C++ |version |
|00004aa0| 33 2e 30 3b 20 72 65 61 | 6c 28 41 29 20 62 65 63 |3.0; rea|l(A) bec|
|00004ab0| 6f 6d 65 73 0a 58 41 2e | 41 73 53 63 61 6c 61 72 |omes.XA.|AsScalar|
|00004ac0| 28 29 3b 20 43 6f 70 79 | 54 6f 4d 61 74 72 69 78 |(); Copy|ToMatrix|
|00004ad0| 20 62 65 63 6f 6d 65 73 | 20 41 73 4d 61 74 72 69 | becomes| AsMatri|
|00004ae0| 78 2c 20 65 74 63 3b 20 | 2e 63 28 29 20 69 73 20 |x, etc; |.c() is |
|00004af0| 6e 6f 20 6c 6f 6e 67 65 | 72 0a 58 72 65 71 75 69 |no longe|r.Xrequi|
|00004b00| 72 65 64 20 28 74 6f 20 | 62 65 20 64 65 6c 65 74 |red (to |be delet|
|00004b10| 65 64 20 69 6e 20 6e 65 | 78 74 20 76 65 72 73 69 |ed in ne|xt versi|
|00004b20| 6f 6e 29 3b 20 6f 70 74 | 69 6f 6e 20 66 6f 72 20 |on); opt|ion for |
|00004b30| 76 65 72 73 69 6f 6e 20 | 32 2e 31 20 6f 72 0a 58 |version |2.1 or.X|
|00004b40| 6c 61 74 65 72 2e 20 53 | 75 66 66 69 78 20 66 6f |later. S|uffix fo|
|00004b50| 72 20 69 6e 63 6c 75 64 | 65 20 66 69 6c 65 73 20 |r includ|e files |
|00004b60| 63 68 61 6e 67 65 64 20 | 74 6f 20 2e 68 3b 20 42 |changed |to .h; B|
|00004b70| 4f 4f 4c 20 63 68 61 6e | 67 65 64 20 74 6f 20 42 |OOL chan|ged to B|
|00004b80| 6f 6f 6c 65 61 6e 0a 58 | 28 42 4f 4f 4c 20 64 6f |oolean.X|(BOOL do|
|00004b90| 65 73 6e 27 74 20 77 6f | 72 6b 20 69 6e 20 67 2b |esn't wo|rk in g+|
|00004ba0| 2b 20 76 20 32 2e 30 29 | 3b 20 6d 6f 64 66 69 63 |+ v 2.0)|; modfic|
|00004bb0| 61 74 69 6f 6e 73 20 74 | 6f 20 61 6c 6c 6f 77 20 |ations t|o allow |
|00004bc0| 66 6f 72 20 63 6f 6d 70 | 69 6c 65 72 73 0a 58 74 |for comp|ilers.Xt|
|00004bd0| 68 61 74 20 64 65 73 74 | 72 6f 79 20 74 65 6d 70 |hat dest|roy temp|
|00004be0| 6f 72 61 72 69 65 73 20 | 76 65 72 79 20 71 75 69 |oraries |very qui|
|00004bf0| 63 6b 6c 79 3b 20 28 47 | 6e 75 20 75 73 65 72 73 |ckly; (G|nu users|
|00004c00| 20 2d 20 73 65 65 20 74 | 68 65 20 73 65 63 74 69 | - see t|he secti|
|00004c10| 6f 6e 20 6f 66 0a 58 63 | 6f 6d 70 69 6c 65 72 20 |on of.Xc|ompiler |
|00004c20| 70 65 72 66 6f 72 6d 61 | 6e 63 65 29 2e 20 41 64 |performa|nce). Ad|
|00004c30| 64 65 64 20 43 6c 65 61 | 6e 55 70 2c 20 4c 69 6e |ded Clea|nUp, Lin|
|00004c40| 65 61 72 45 71 75 61 74 | 69 6f 6e 53 6f 6c 76 65 |earEquat|ionSolve|
|00004c50| 72 2c 20 70 72 69 6d 69 | 74 69 76 65 0a 58 76 65 |r, primi|tive.Xve|
|00004c60| 72 73 69 6f 6e 20 6f 66 | 20 65 78 63 65 70 74 69 |rsion of| excepti|
|00004c70| 6f 6e 73 2e 0a 58 0a 58 | 0a 58 4e 65 77 6d 61 74 |ons..X.X|.XNewmat|
|00004c80| 30 35 20 2d 20 4a 75 6e | 65 20 31 39 39 32 3a 0a |05 - Jun|e 1992:.|
|00004c90| 58 0a 58 46 6f 72 20 70 | 72 69 76 61 74 65 20 72 |X.XFor p|rivate r|
|00004ca0| 65 6c 65 61 73 65 20 6f | 6e 6c 79 20 0a 58 0a 58 |elease o|nly .X.X|
|00004cb0| 0a 58 4e 65 77 6d 61 74 | 30 34 20 2d 20 44 65 63 |.XNewmat|04 - Dec|
|00004cc0| 65 6d 62 65 72 20 31 39 | 39 31 3a 0a 58 0a 58 46 |ember 19|91:.X.XF|
|00004cd0| 69 78 20 70 72 6f 62 6c | 65 6d 20 77 69 74 68 20 |ix probl|em with |
|00004ce0| 47 2b 2b 31 2e 34 30 2c | 20 73 6f 6d 65 20 65 78 |G++1.40,| some ex|
|00004cf0| 74 72 61 20 64 6f 63 75 | 6d 65 6e 74 61 74 69 6f |tra docu|mentatio|
|00004d00| 6e 0a 58 0a 58 0a 58 4e | 65 77 6d 61 74 30 33 20 |n.X.X.XN|ewmat03 |
|00004d10| 2d 20 4e 6f 76 65 6d 62 | 65 72 20 31 39 39 31 3a |- Novemb|er 1991:|
|00004d20| 0a 58 0a 58 43 6f 6c 20 | 61 6e 64 20 43 6f 6c 73 |.X.XCol |and Cols|
|00004d30| 20 62 65 63 6f 6d 65 20 | 43 6f 6c 75 6d 6e 20 61 | become |Column a|
|00004d40| 6e 64 20 43 6f 6c 75 6d | 6e 73 2e 20 41 64 64 65 |nd Colum|ns. Adde|
|00004d50| 64 20 53 6f 72 74 2c 20 | 53 56 44 2c 20 4a 61 63 |d Sort, |SVD, Jac|
|00004d60| 6f 62 69 2c 0a 58 45 69 | 67 65 6e 76 61 6c 75 65 |obi,.XEi|genvalue|
|00004d70| 73 2c 20 46 46 54 2c 20 | 72 65 61 6c 20 63 6f 6e |s, FFT, |real con|
|00004d80| 76 65 72 73 69 6f 6e 20 | 6f 66 20 31 78 31 20 6d |version |of 1x1 m|
|00004d90| 61 74 72 69 78 2c 20 22 | 4e 75 6d 65 72 69 63 61 |atrix, "|Numerica|
|00004da0| 6c 20 52 65 63 69 70 65 | 73 20 69 6e 0a 58 43 22 |l Recipe|s in.XC"|
|00004db0| 20 69 6e 74 65 72 66 61 | 63 65 2c 20 6f 75 74 70 | interfa|ce, outp|
|00004dc0| 75 74 20 6f 70 65 72 61 | 74 69 6f 6e 73 2c 20 76 |ut opera|tions, v|
|00004dd0| 61 72 69 6f 75 73 20 73 | 63 61 6c 61 72 20 66 75 |arious s|calar fu|
|00004de0| 6e 63 74 69 6f 6e 73 2e | 20 49 6d 70 72 6f 76 65 |nctions.| Improve|
|00004df0| 64 0a 58 72 65 74 75 72 | 6e 20 66 72 6f 6d 20 66 |d.Xretur|n from f|
|00004e00| 75 6e 63 74 69 6f 6e 73 | 2e 20 52 65 6f 72 67 61 |unctions|. Reorga|
|00004e10| 6e 69 73 65 64 20 73 65 | 74 74 69 6e 67 20 6f 70 |nised se|tting op|
|00004e20| 74 69 6f 6e 73 20 69 6e | 20 22 69 6e 63 6c 75 64 |tions in| "includ|
|00004e30| 65 2e 68 78 78 22 2e 0a | 58 0a 58 0a 58 4e 65 77 |e.hxx"..|X.X.XNew|
|00004e40| 6d 61 74 30 32 20 2d 20 | 4a 75 6c 79 20 31 39 39 |mat02 - |July 199|
|00004e50| 31 3a 0a 58 0a 58 56 65 | 72 73 69 6f 6e 20 77 69 |1:.X.XVe|rsion wi|
|00004e60| 74 68 20 6d 61 74 72 69 | 78 20 72 6f 77 2f 63 6f |th matri|x row/co|
|00004e70| 6c 75 6d 6e 20 6f 70 65 | 72 61 74 69 6f 6e 73 20 |lumn ope|rations |
|00004e80| 61 6e 64 20 6e 75 6d 65 | 72 6f 75 73 20 61 64 64 |and nume|rous add|
|00004e90| 69 74 69 6f 6e 61 6c 0a | 58 66 75 6e 63 74 69 6f |itional.|Xfunctio|
|00004ea0| 6e 73 2e 0a 58 0a 58 0a | 58 4d 61 74 72 69 78 20 |ns..X.X.|XMatrix |
|00004eb0| 2d 20 4f 63 74 6f 62 65 | 72 20 31 39 39 30 3a 0a |- Octobe|r 1990:.|
|00004ec0| 58 0a 58 45 61 72 6c 79 | 20 76 65 72 73 69 6f 6e |X.XEarly| version|
|00004ed0| 20 6f 66 20 70 61 63 6b | 61 67 65 2e 0a 58 0a 58 | of pack|age..X.X|
|00004ee0| 0a 58 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |.X------|--------|
|00004ef0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004f00| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004f10| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00004f20| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 0a 58 0a |--------|-----.X.|
|00004f30| 58 0a 58 48 6f 77 20 74 | 6f 20 67 65 74 20 61 20 |X.XHow t|o get a |
|00004f40| 63 6f 70 79 20 6f 66 20 | 74 68 69 73 20 70 61 63 |copy of |this pac|
|00004f50| 6b 61 67 65 0a 58 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |kage.X==|========|
|00004f60| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00004f70| 3d 3d 3d 3d 3d 3d 3d 0a | 58 0a 58 49 20 61 6d 20 |=======.|X.XI am |
|00004f80| 70 75 74 74 69 6e 67 20 | 63 6f 70 69 65 73 20 6f |putting |copies o|
|00004f90| 6e 20 43 6f 6d 70 75 73 | 65 72 76 65 20 28 42 6f |n Compus|erve (Bo|
|00004fa0| 72 6c 61 6e 64 20 6c 69 | 62 72 61 72 79 2c 20 7a |rland li|brary, z|
|00004fb0| 69 70 20 66 6f 72 6d 61 | 74 29 2c 0a 58 53 49 4d |ip forma|t),.XSIM|
|00004fc0| 54 45 4c 32 30 20 28 4d | 73 44 6f 73 20 6c 69 62 |TEL20 (M|sDos lib|
|00004fd0| 72 61 72 79 2c 20 7a 69 | 70 20 66 6f 72 6d 61 74 |rary, zi|p format|
|00004fe0| 29 2c 20 63 6f 6d 70 2e | 73 6f 75 72 63 65 73 2e |), comp.|sources.|
|00004ff0| 6d 69 73 63 20 6f 6e 20 | 49 6e 74 65 72 6e 65 74 |misc on |Internet|
|00005000| 0a 58 28 73 68 61 72 20 | 66 6f 72 6d 61 74 29 2e |.X(shar |format).|
|00005010| 0a 58 0a 58 0a 58 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |.X.X.X--|--------|
|00005020| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005030| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005040| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005050| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00005060| 2d 0a 58 0a 58 0a 58 43 | 6f 6d 70 69 6c 65 72 20 |-.X.X.XC|ompiler |
|00005070| 70 65 72 66 6f 72 6d 61 | 6e 63 65 0a 58 3d 3d 3d |performa|nce.X===|
|00005080| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00005090| 3d 0a 58 0a 58 49 20 68 | 61 76 65 20 74 65 73 74 |=.X.XI h|ave test|
|000050a0| 65 64 20 74 68 69 73 20 | 70 61 63 6b 61 67 65 20 |ed this |package |
|000050b0| 6f 6e 20 61 20 6e 75 6d | 62 65 72 20 6f 66 20 63 |on a num|ber of c|
|000050c0| 6f 6d 70 69 6c 65 72 73 | 2e 20 48 65 72 65 20 61 |ompilers|. Here a|
|000050d0| 72 65 20 74 68 65 20 6c | 65 76 65 6c 73 0a 58 6f |re the l|evels.Xo|
|000050e0| 66 20 73 75 63 63 65 73 | 73 20 77 69 74 68 20 74 |f succes|s with t|
|000050f0| 68 69 73 20 70 61 63 6b | 61 67 65 2e 20 49 6e 20 |his pack|age. In |
|00005100| 6d 6f 73 74 20 63 61 73 | 65 73 20 49 20 68 61 76 |most cas|es I hav|
|00005110| 65 20 63 68 6f 73 65 6e | 20 63 6f 64 65 20 74 68 |e chosen| code th|
|00005120| 61 74 0a 58 77 6f 72 6b | 73 20 75 6e 64 65 72 20 |at.Xwork|s under |
|00005130| 61 6c 6c 20 74 68 65 20 | 63 6f 6d 70 69 6c 65 72 |all the |compiler|
|00005140| 73 20 49 20 68 61 76 65 | 20 61 63 63 65 73 73 20 |s I have| access |
|00005150| 74 6f 2c 20 62 75 74 20 | 49 20 68 61 76 65 20 68 |to, but |I have h|
|00005160| 61 64 20 74 6f 0a 58 69 | 6e 63 6c 75 64 65 20 73 |ad to.Xi|nclude s|
|00005170| 6f 6d 65 20 73 70 65 63 | 69 66 69 63 20 77 6f 72 |ome spec|ific wor|
|00005180| 6b 2d 61 72 6f 75 6e 64 | 73 20 66 6f 72 20 73 6f |k-around|s for so|
|00005190| 6d 65 20 63 6f 6d 70 69 | 6c 65 72 73 2e 20 46 6f |me compi|lers. Fo|
|000051a0| 72 20 74 68 65 20 4d 73 | 44 6f 73 0a 58 76 65 72 |r the Ms|Dos.Xver|
|000051b0| 73 69 6f 6e 73 2c 20 49 | 20 75 73 65 20 61 20 34 |sions, I| use a 4|
|000051c0| 38 36 64 78 20 63 6f 6d | 70 75 74 65 72 20 72 75 |86dx com|puter ru|
|000051d0| 6e 6e 69 6e 67 20 4d 73 | 44 6f 73 20 35 2e 20 54 |nning Ms|Dos 5. T|
|000051e0| 68 65 20 75 6e 69 78 20 | 76 65 72 73 69 6f 6e 73 |he unix |versions|
|000051f0| 20 61 72 65 0a 58 6f 6e | 20 61 20 53 75 6e 20 53 | are.Xon| a Sun S|
|00005200| 70 61 72 63 20 73 74 61 | 74 69 6f 6e 20 6f 72 20 |parc sta|tion or |
|00005210| 61 20 53 69 6c 69 63 6f | 6e 20 47 72 61 70 68 69 |a Silico|n Graphi|
|00005220| 63 73 20 6f 72 20 61 20 | 48 50 20 75 6e 69 78 20 |cs or a |HP unix |
|00005230| 77 6f 72 6b 73 74 61 74 | 69 6f 6e 2e 0a 58 54 68 |workstat|ion..XTh|
|00005240| 61 6e 6b 73 20 74 6f 20 | 56 69 63 74 6f 72 69 61 |anks to |Victoria|
|00005250| 20 55 6e 69 76 65 72 73 | 69 74 79 20 61 6e 64 20 | Univers|ity and |
|00005260| 49 6e 64 75 73 74 72 69 | 61 6c 20 52 65 73 65 61 |Industri|al Resea|
|00005270| 72 63 68 20 4c 74 64 20 | 66 6f 72 20 61 63 63 65 |rch Ltd |for acce|
|00005280| 73 73 20 74 6f 0a 58 74 | 68 65 20 55 6e 69 78 20 |ss to.Xt|he Unix |
|00005290| 6d 61 63 68 69 6e 65 73 | 2e 0a 58 0a 58 41 20 73 |machines|..X.XA s|
|000052a0| 65 72 69 65 73 20 6f 66 | 20 23 64 65 66 69 6e 65 |eries of| #define|
|000052b0| 73 20 61 74 20 74 68 65 | 20 62 65 67 69 6e 6e 69 |s at the| beginni|
|000052c0| 6e 67 20 6f 66 20 22 69 | 6e 63 6c 75 64 65 2e 68 |ng of "i|nclude.h|
|000052d0| 22 20 63 75 73 74 6f 6d | 69 73 65 73 20 74 68 65 |" custom|ises the|
|000052e0| 0a 58 70 61 63 6b 61 67 | 65 20 66 6f 72 20 74 68 |.Xpackag|e for th|
|000052f0| 65 20 63 6f 6d 70 69 6c | 65 72 20 79 6f 75 20 61 |e compil|er you a|
|00005300| 72 65 20 75 73 69 6e 67 | 2e 20 54 75 72 62 6f 2c |re using|. Turbo,|
|00005310| 20 42 6f 72 6c 61 6e 64 | 2c 20 47 6e 75 20 61 6e | Borland|, Gnu an|
|00005320| 64 20 5a 6f 72 74 65 63 | 68 0a 58 61 72 65 20 72 |d Zortec|h.Xare r|
|00005330| 65 63 6f 67 6e 69 73 65 | 64 20 61 75 74 6f 6d 61 |ecognise|d automa|
|00005340| 74 69 63 61 6c 6c 79 2c | 20 6f 74 68 65 72 77 69 |tically,| otherwi|
|00005350| 73 65 20 79 6f 75 20 68 | 61 76 65 20 74 6f 20 73 |se you h|ave to s|
|00005360| 65 74 20 74 68 65 20 61 | 70 70 72 6f 70 72 69 61 |et the a|ppropria|
|00005370| 74 65 0a 58 23 64 65 66 | 69 6e 65 20 73 74 61 74 |te.X#def|ine stat|
|00005380| 65 6d 65 6e 74 2e 20 41 | 63 74 69 76 61 74 65 20 |ement. A|ctivate |
|00005390| 74 68 65 20 6f 70 74 69 | 6f 6e 20 66 6f 72 20 76 |the opti|on for v|
|000053a0| 65 72 73 69 6f 6e 20 32 | 2e 31 20 69 66 20 79 6f |ersion 2|.1 if yo|
|000053b0| 75 20 61 72 65 20 75 73 | 69 6e 67 0a 58 76 65 72 |u are us|ing.Xver|
|000053c0| 73 69 6f 6e 20 32 2e 31 | 20 6f 66 20 43 2b 2b 20 |sion 2.1| of C++ |
|000053d0| 6f 72 20 6c 61 74 65 72 | 2e 0a 58 0a 58 42 6f 72 |or later|..X.XBor|
|000053e0| 6c 61 6e 64 20 43 2b 2b | 20 33 2e 31 3a 20 52 65 |land C++| 3.1: Re|
|000053f0| 63 65 6e 74 6c 79 20 74 | 68 69 73 20 68 61 73 20 |cently t|his has |
|00005400| 62 65 65 6e 20 6d 79 20 | 6d 61 69 6e 20 64 65 76 |been my |main dev|
|00005410| 65 6c 6f 70 6d 65 6e 74 | 20 70 6c 61 74 66 6f 72 |elopment| platfor|
|00005420| 6d 2c 0a 58 73 6f 20 6e | 61 74 75 72 61 6c 6c 79 |m,.Xso n|aturally|
|00005430| 20 61 6c 6d 6f 73 74 20 | 65 76 65 72 79 74 68 69 | almost |everythi|
|00005440| 6e 67 20 77 6f 72 6b 73 | 20 77 69 74 68 20 74 68 |ng works| with th|
|00005450| 69 73 20 63 6f 6d 70 69 | 6c 65 72 2e 20 4d 61 6b |is compi|ler. Mak|
|00005460| 65 20 73 75 72 65 20 79 | 6f 75 0a 58 68 61 76 65 |e sure y|ou.Xhave|
|00005470| 20 74 68 65 20 63 6f 6d | 70 69 6c 65 72 20 6f 70 | the com|piler op|
|00005480| 74 69 6f 6e 20 22 74 72 | 65 61 74 20 65 6e 75 6d |tion "tr|eat enum|
|00005490| 73 20 61 73 20 69 6e 74 | 73 22 20 73 65 74 2e 20 |s as int|s" set. |
|000054a0| 54 68 65 72 65 20 77 61 | 73 20 61 20 70 72 6f 62 |There wa|s a prob|
|000054b0| 6c 65 6d 0a 58 77 69 74 | 68 20 74 68 65 20 6c 69 |lem.Xwit|h the li|
|000054c0| 62 72 61 72 79 20 75 74 | 69 6c 69 74 79 20 69 6e |brary ut|ility in|
|000054d0| 20 76 65 72 73 69 6f 6e | 20 32 2e 30 20 77 68 69 | version| 2.0 whi|
|000054e0| 63 68 20 69 73 20 6e 6f | 77 20 66 69 78 65 64 2e |ch is no|w fixed.|
|000054f0| 20 59 6f 75 20 77 69 6c | 6c 0a 58 6e 65 65 64 20 | You wil|l.Xneed |
|00005500| 74 6f 20 75 73 65 20 74 | 68 65 20 6c 61 72 67 65 |to use t|he large|
|00005510| 20 6d 6f 64 65 6c 2e 20 | 49 66 20 79 6f 75 20 61 | model. |If you a|
|00005520| 72 65 20 6e 6f 74 20 64 | 65 62 75 67 67 69 6e 67 |re not d|ebugging|
|00005530| 2c 20 74 75 72 6e 20 6f | 66 66 20 74 68 65 0a 58 |, turn o|ff the.X|
|00005540| 6f 70 74 69 6f 6e 73 20 | 74 68 61 74 20 63 6f 6c |options |that col|
|00005550| 6c 65 63 74 20 64 65 62 | 75 67 67 69 6e 67 20 69 |lect deb|ugging i|
|00005560| 6e 66 6f 72 6d 61 74 69 | 6f 6e 2e 0a 58 0a 58 4d |nformati|on..X.XM|
|00005570| 69 63 72 6f 73 6f 66 74 | 20 43 2b 2b 20 28 37 2e |icrosoft| C++ (7.|
|00005580| 30 29 3a 20 53 65 65 6d | 73 20 74 6f 20 77 6f 72 |0): Seem|s to wor|
|00005590| 6b 20 4f 4b 2e 20 59 6f | 75 20 6d 75 73 74 20 23 |k OK. Yo|u must #|
|000055a0| 64 65 66 69 6e 65 0a 58 | 54 45 4d 50 53 5f 44 45 |define.X|TEMPS_DE|
|000055b0| 53 54 52 4f 59 45 44 5f | 51 55 49 43 4b 4c 59 20 |STROYED_|QUICKLY |
|000055c0| 6f 77 69 6e 67 20 74 6f | 20 61 20 62 75 67 20 69 |owing to| a bug i|
|000055d0| 6e 20 74 68 65 20 63 75 | 72 72 65 6e 74 20 76 65 |n the cu|rrent ve|
|000055e0| 72 73 69 6f 6e 20 6f 66 | 20 4d 53 43 2e 0a 58 0a |rsion of| MSC..X.|
|000055f0| 58 5a 6f 72 74 65 63 68 | 20 43 2b 2b 20 33 2e 30 |XZortech| C++ 3.0|
|00005600| 3a 20 22 63 6f 6e 73 74 | 22 20 64 6f 65 73 6e 27 |: "const|" doesn'|
|00005610| 74 20 77 6f 72 6b 20 63 | 6f 72 72 65 63 74 6c 79 |t work c|orrectly|
|00005620| 20 77 69 74 68 20 74 68 | 69 73 20 63 6f 6d 70 69 | with th|is compi|
|00005630| 6c 65 72 2c 20 73 6f 0a | 58 74 68 65 20 70 61 63 |ler, so.|Xthe pac|
|00005640| 6b 61 67 65 20 73 6b 69 | 70 73 20 61 6c 6c 20 6f |kage ski|ps all o|
|00005650| 66 20 74 68 65 20 73 74 | 61 74 65 6d 65 6e 74 73 |f the st|atements|
|00005660| 20 5a 6f 72 74 65 63 68 | 20 63 61 6e 27 74 20 68 | Zortech| can't h|
|00005670| 61 6e 64 6c 65 2e 20 5a | 6f 72 74 65 63 68 0a 58 |andle. Z|ortech.X|
|00005680| 6c 65 61 76 65 73 20 72 | 75 62 62 69 73 68 20 6f |leaves r|ubbish o|
|00005690| 6e 20 74 68 65 20 68 65 | 61 70 2e 20 49 20 64 6f |n the he|ap. I do|
|000056a0| 6e 27 74 20 6b 6e 6f 77 | 20 77 68 65 74 68 65 72 |n't know| whether|
|000056b0| 20 74 68 69 73 20 69 73 | 20 6d 79 20 70 72 6f 67 | this is| my prog|
|000056c0| 72 61 6d 6d 69 6e 67 0a | 58 65 72 72 6f 72 20 6f |ramming.|Xerror o|
|000056d0| 72 20 61 20 5a 6f 72 74 | 65 63 68 20 65 72 72 6f |r a Zort|ech erro|
|000056e0| 72 20 6f 72 20 61 64 64 | 69 74 69 6f 6e 61 6c 20 |r or add|itional |
|000056f0| 70 72 69 6e 74 65 72 20 | 62 75 66 66 65 72 73 2e |printer |buffers.|
|00005700| 20 44 65 61 63 74 69 76 | 61 74 65 20 74 68 65 0a | Deactiv|ate the.|
|00005710| 58 6f 70 74 69 6f 6e 20 | 66 6f 72 20 76 65 72 73 |Xoption |for vers|
|00005720| 69 6f 6e 20 32 2e 31 20 | 69 6e 20 69 6e 63 6c 75 |ion 2.1 |in inclu|
|00005730| 64 65 2e 68 2e 20 44 6f | 65 73 20 6e 6f 74 20 73 |de.h. Do|es not s|
|00005740| 75 70 70 6f 72 74 20 49 | 4f 20 6d 61 6e 69 70 75 |upport I|O manipu|
|00005750| 6c 61 74 6f 72 73 2e 0a | 58 4f 74 68 65 72 77 69 |lators..|XOtherwi|
|00005760| 73 65 20 74 68 65 20 70 | 61 63 6b 61 67 65 20 6d |se the p|ackage m|
|00005770| 6f 73 74 6c 79 20 77 6f | 72 6b 73 2c 20 62 75 74 |ostly wo|rks, but|
|00005780| 20 6e 6f 74 20 63 6f 6d | 70 6c 65 74 65 6c 79 2e | not com|pletely.|
|00005790| 20 42 65 73 74 20 64 6f | 6e 27 74 0a 58 23 64 65 | Best do|n't.X#de|
|000057a0| 66 69 6e 65 20 54 45 4d | 50 53 5f 44 45 53 54 52 |fine TEM|PS_DESTR|
|000057b0| 4f 59 45 44 5f 51 55 49 | 43 4b 4c 59 2e 20 45 78 |OYED_QUI|CKLY. Ex|
|000057c0| 63 65 70 74 69 6f 6e 73 | 20 61 6e 64 20 74 68 65 |ceptions| and the|
|000057d0| 20 6e 72 69 63 20 69 6e | 74 65 72 66 61 63 65 20 | nric in|terface |
|000057e0| 64 6f 6e 27 74 0a 58 77 | 6f 72 6b 2e 20 49 20 74 |don't.Xw|ork. I t|
|000057f0| 68 69 6e 6b 20 74 68 65 | 20 70 72 6f 62 6c 65 6d |hink the| problem|
|00005800| 73 20 61 72 65 20 62 65 | 63 61 75 73 65 20 5a 6f |s are be|cause Zo|
|00005810| 72 74 65 63 68 20 64 6f | 65 73 6e 27 74 20 68 61 |rtech do|esn't ha|
|00005820| 6e 64 6c 65 0a 58 63 6f | 6e 76 65 72 73 69 6f 6e |ndle.Xco|nversion|
|00005830| 73 20 63 6f 72 72 65 63 | 74 6c 79 2c 20 70 61 72 |s correc|tly, par|
|00005840| 74 69 63 75 6c 61 72 6c | 79 20 61 75 74 6f 6d 61 |ticularl|y automa|
|00005850| 74 69 63 20 63 6f 6e 76 | 65 72 73 69 6f 6e 73 2e |tic conv|ersions.|
|00005860| 20 5a 6f 72 74 65 63 68 | 20 72 75 6e 73 0a 58 6d | Zortech| runs.Xm|
|00005870| 75 63 68 20 6d 6f 72 65 | 20 73 6c 6f 77 6c 79 20 |uch more| slowly |
|00005880| 74 68 61 6e 20 42 6f 72 | 6c 61 6e 64 20 61 6e 64 |than Bor|land and|
|00005890| 20 4d 69 63 72 6f 73 6f | 66 74 2e 20 55 73 65 20 | Microso|ft. Use |
|000058a0| 74 68 65 20 6c 61 72 67 | 65 20 6d 6f 64 65 6c 20 |the larg|e model |
|000058b0| 61 6e 64 0a 58 6f 70 74 | 69 6d 69 73 61 74 69 6f |and.Xopt|imisatio|
|000058c0| 6e 2e 0a 58 0a 58 47 6c | 6f 63 6b 65 6e 73 70 69 |n..X.XGl|ockenspi|
|000058d0| 65 6c 20 43 2b 2b 20 28 | 32 2e 30 30 61 20 66 6f |el C++ (|2.00a fo|
|000058e0| 72 20 4d 73 44 6f 73 20 | 6c 6f 61 64 69 6e 67 20 |r MsDos |loading |
|000058f0| 69 6e 74 6f 20 4d 69 63 | 72 6f 73 6f 66 74 20 43 |into Mic|rosoft C|
|00005900| 20 35 2e 31 29 3a 20 49 | 0a 58 68 61 76 65 6e 27 | 5.1): I|.Xhaven'|
|00005910| 74 20 74 65 73 74 65 64 | 20 74 68 65 20 6c 61 74 |t tested| the lat|
|00005920| 65 73 74 20 76 65 72 73 | 69 6f 6e 20 6f 66 20 6d |est vers|ion of m|
|00005930| 79 20 70 61 63 6b 61 67 | 65 20 77 69 74 68 20 47 |y packag|e with G|
|00005940| 6c 6f 63 6b 65 6e 73 70 | 69 65 6c 2e 20 49 20 68 |lockensp|iel. I h|
|00005950| 61 64 0a 58 74 6f 20 23 | 64 65 66 69 6e 65 20 74 |ad.Xto #|define t|
|00005960| 68 65 20 6d 61 74 72 69 | 78 20 6e 61 6d 65 73 20 |he matri|x names |
|00005970| 74 6f 20 73 68 6f 72 74 | 65 72 20 6e 61 6d 65 73 |to short|er names|
|00005980| 20 74 6f 20 61 76 6f 69 | 64 20 61 6d 62 69 67 75 | to avoi|d ambigu|
|00005990| 69 74 69 65 73 20 61 6e | 64 0a 58 68 61 64 20 71 |ities an|d.Xhad q|
|000059a0| 75 69 74 65 20 61 20 62 | 69 74 20 6f 66 20 64 69 |uite a b|it of di|
|000059b0| 66 66 69 63 75 6c 74 79 | 20 73 74 6f 70 70 69 6e |fficulty| stoppin|
|000059c0| 67 20 74 68 65 20 63 6f | 6d 70 69 6c 65 73 20 66 |g the co|mpiles f|
|000059d0| 72 6f 6d 20 72 75 6e 6e | 69 6e 67 20 6f 75 74 20 |rom runn|ing out |
|000059e0| 6f 66 0a 58 73 70 61 63 | 65 20 61 6e 64 20 6e 6f |of.Xspac|e and no|
|000059f0| 74 20 65 78 63 65 65 64 | 69 6e 67 20 4d 69 63 72 |t exceed|ing Micr|
|00005a00| 6f 73 6f 66 74 27 73 20 | 62 6c 6f 63 6b 20 6e 65 |osoft's |block ne|
|00005a10| 73 74 69 6e 67 20 6c 69 | 6d 69 74 2e 20 41 20 63 |sting li|mit. A c|
|00005a20| 6f 75 70 6c 65 20 6f 66 | 20 6d 79 0a 58 74 65 73 |ouple of| my.Xtes|
|00005a30| 74 20 73 74 61 74 65 6d | 65 6e 74 73 20 70 72 6f |t statem|ents pro|
|00005a40| 64 75 63 65 64 20 73 74 | 61 74 65 6d 65 6e 74 73 |duced st|atements|
|00005a50| 20 74 6f 6f 20 63 6f 6d | 70 6c 65 78 20 66 6f 72 | too com|plex for|
|00005a60| 20 4d 69 63 72 6f 73 6f | 66 74 2c 20 62 75 74 0a | Microso|ft, but.|
|00005a70| 58 62 61 73 69 63 61 6c | 6c 79 20 74 68 65 20 70 |Xbasical|ly the p|
|00005a80| 61 63 6b 61 67 65 20 77 | 6f 72 6b 65 64 2e 20 54 |ackage w|orked. T|
|00005a90| 68 69 73 20 77 61 73 20 | 6d 79 20 6f 72 69 67 69 |his was |my origi|
|00005aa0| 6e 61 6c 20 64 65 76 65 | 6c 6f 70 6d 65 6e 74 20 |nal deve|lopment |
|00005ab0| 70 6c 61 74 66 6f 72 6d | 0a 58 61 6e 64 20 49 20 |platform|.Xand I |
|00005ac0| 73 74 69 6c 6c 20 75 73 | 65 20 2e 63 78 78 20 61 |still us|e .cxx a|
|00005ad0| 73 20 6d 79 20 66 69 6c | 65 20 6e 61 6d 65 20 65 |s my fil|e name e|
|00005ae0| 78 74 65 6e 73 69 6f 6e | 73 20 66 6f 72 20 74 68 |xtension|s for th|
|00005af0| 65 20 43 2b 2b 20 66 69 | 6c 65 73 2e 0a 58 0a 58 |e C++ fi|les..X.X|
|00005b00| 53 75 6e 20 41 54 26 54 | 20 43 2b 2b 20 32 2e 31 |Sun AT&T| C++ 2.1|
|00005b10| 3b 33 2e 30 3a 20 54 68 | 69 73 20 77 6f 72 6b 73 |;3.0: Th|is works|
|00005b20| 20 66 69 6e 65 2e 20 45 | 78 63 65 70 74 20 61 67 | fine. E|xcept ag|
|00005b30| 67 72 65 67 61 74 65 73 | 20 61 72 65 20 6e 6f 74 |gregates| are not|
|00005b40| 0a 58 73 75 70 70 6f 72 | 74 65 64 20 69 6e 20 32 |.Xsuppor|ted in 2|
|00005b50| 2e 31 20 61 6e 64 20 73 | 65 74 6a 6d 70 2e 68 20 |.1 and s|etjmp.h |
|00005b60| 67 65 6e 65 72 61 74 65 | 64 20 61 20 77 61 72 6e |generate|d a warn|
|00005b70| 69 6e 67 20 6d 65 73 73 | 61 67 65 2e 20 4e 65 69 |ing mess|age. Nei|
|00005b80| 74 68 65 72 0a 58 63 6f | 6d 70 69 6c 65 72 20 77 |ther.Xco|mpiler w|
|00005b90| 6f 75 6c 64 20 63 6f 6d | 70 69 6c 65 20 77 68 65 |ould com|pile whe|
|00005ba0| 6e 20 49 20 73 65 74 20 | 44 4f 5f 46 52 45 45 5f |n I set |DO_FREE_|
|00005bb0| 43 48 45 43 4b 20 28 73 | 65 65 20 6d 79 20 66 69 |CHECK (s|ee my fi|
|00005bc0| 6c 65 0a 58 6e 65 77 6d | 61 74 63 2e 74 78 74 29 |le.Xnewm|atc.txt)|
|00005bd0| 2e 20 49 66 20 79 6f 75 | 20 61 72 65 20 75 73 69 |. If you| are usi|
|00005be0| 6e 67 20 22 69 6e 74 65 | 72 76 69 65 77 73 22 20 |ng "inte|rviews" |
|00005bf0| 79 6f 75 20 6d 61 79 20 | 67 65 74 20 61 20 63 6f |you may |get a co|
|00005c00| 6e 66 6c 69 63 74 20 77 | 69 74 68 0a 58 43 61 74 |nflict w|ith.XCat|
|00005c10| 63 68 2e 20 45 69 74 68 | 65 72 20 23 75 6e 64 65 |ch. Eith|er #unde|
|00005c20| 66 69 6e 65 20 43 61 74 | 63 68 20 6f 72 20 72 65 |fine Cat|ch or re|
|00005c30| 70 6c 61 63 65 20 43 61 | 74 63 68 20 77 69 74 68 |place Ca|tch with|
|00005c40| 20 43 41 54 43 48 20 74 | 68 72 6f 75 67 68 6f 75 | CATCH t|hroughou|
|00005c50| 74 20 6d 79 0a 58 70 61 | 63 6b 61 67 65 2e 20 49 |t my.Xpa|ckage. I|
|00005c60| 6e 20 41 54 26 54 20 32 | 2e 31 20 79 6f 75 20 6d |n AT&T 2|.1 you m|
|00005c70| 61 79 20 67 65 74 20 61 | 6e 20 65 72 72 6f 72 20 |ay get a|n error |
|00005c80| 77 68 65 6e 20 79 6f 75 | 20 75 73 65 20 61 6e 20 |when you| use an |
|00005c90| 65 78 70 72 65 73 73 69 | 6f 6e 20 66 6f 72 0a 58 |expressi|on for.X|
|00005ca0| 74 68 65 20 73 69 6e 67 | 6c 65 20 61 72 67 75 6d |the sing|le argum|
|00005cb0| 65 6e 74 20 77 68 65 6e | 20 63 6f 6e 73 74 72 75 |ent when| constru|
|00005cc0| 63 74 69 6e 67 20 61 20 | 56 65 63 74 6f 72 20 6f |cting a |Vector o|
|00005cd0| 72 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |r Diagon|alMatrix|
|00005ce0| 20 6f 72 20 6f 6e 65 0a | 58 6f 66 20 74 68 65 20 | or one.|Xof the |
|00005cf0| 54 72 69 61 6e 67 75 6c | 61 72 20 4d 61 74 72 69 |Triangul|ar Matri|
|00005d00| 63 65 73 2e 20 59 6f 75 | 20 6e 65 65 64 20 74 6f |ces. You| need to|
|00005d10| 20 65 76 61 6c 75 61 74 | 65 20 74 68 65 20 65 78 | evaluat|e the ex|
|00005d20| 70 72 65 73 73 69 6f 6e | 0a 58 73 65 70 61 72 61 |pression|.Xsepara|
|00005d30| 74 65 6c 79 2e 0a 58 0a | 58 47 6e 75 20 47 2b 2b |tely..X.|XGnu G++|
|00005d40| 20 32 2e 32 3a 20 20 54 | 68 69 73 20 6d 6f 73 74 | 2.2: T|his most|
|00005d50| 6c 79 20 77 6f 72 6b 73 | 2e 20 59 6f 75 20 63 61 |ly works|. You ca|
|00005d60| 6e 27 74 20 75 73 65 20 | 65 78 70 72 65 73 73 69 |n't use |expressi|
|00005d70| 6f 6e 73 20 6c 69 6b 65 | 0a 58 4d 61 74 72 69 78 |ons like|.XMatrix|
|00005d80| 28 58 2a 59 29 20 69 6e | 20 74 68 65 20 6d 69 64 |(X*Y) in| the mid|
|00005d90| 64 6c 65 20 6f 66 20 61 | 6e 20 65 78 70 72 65 73 |dle of a|n expres|
|00005da0| 73 69 6f 6e 20 61 6e 64 | 20 28 4d 61 74 72 69 78 |sion and| (Matrix|
|00005db0| 29 28 58 2a 59 29 20 69 | 73 0a 58 75 6e 72 65 6c |)(X*Y) i|s.Xunrel|
|00005dc0| 69 61 62 6c 65 2e 20 49 | 66 20 79 6f 75 20 77 72 |iable. I|f you wr|
|00005dd0| 69 74 65 20 61 20 66 75 | 6e 63 74 69 6f 6e 20 72 |ite a fu|nction r|
|00005de0| 65 74 75 72 6e 69 6e 67 | 20 61 20 6d 61 74 72 69 |eturning| a matri|
|00005df0| 78 2c 20 79 6f 75 20 4d | 55 53 54 20 75 73 65 20 |x, you M|UST use |
|00005e00| 74 68 65 0a 58 52 65 74 | 75 72 6e 4d 61 74 72 69 |the.XRet|urnMatri|
|00005e10| 78 20 6d 65 74 68 6f 64 | 20 64 65 73 63 72 69 62 |x method| describ|
|00005e20| 65 64 20 69 6e 20 74 68 | 69 73 20 64 6f 63 75 6d |ed in th|is docum|
|00005e30| 65 6e 74 61 74 69 6f 6e | 2e 20 54 68 69 73 20 69 |entation|. This i|
|00005e40| 73 20 62 65 63 61 75 73 | 65 20 67 2b 2b 0a 58 64 |s becaus|e g++.Xd|
|00005e50| 65 73 74 72 6f 79 73 20 | 74 65 6d 70 6f 72 61 72 |estroys |temporar|
|00005e60| 69 65 73 20 6f 63 63 75 | 72 69 6e 67 20 69 6e 20 |ies occu|ring in |
|00005e70| 61 6e 20 65 78 70 72 65 | 73 73 69 6f 6e 20 74 6f |an expre|ssion to|
|00005e80| 6f 20 73 6f 6f 6e 20 66 | 6f 72 20 74 68 65 20 74 |o soon f|or the t|
|00005e90| 77 6f 0a 58 73 74 61 67 | 65 20 77 61 79 20 6f 66 |wo.Xstag|e way of|
|00005ea0| 20 65 76 61 6c 75 61 74 | 69 6e 67 20 65 78 70 72 | evaluat|ing expr|
|00005eb0| 65 73 73 69 6f 6e 73 20 | 74 68 61 74 20 6e 65 77 |essions |that new|
|00005ec0| 6d 61 74 20 75 73 65 73 | 2e 20 47 6e 75 20 32 2e |mat uses|. Gnu 2.|
|00005ed0| 32 20 64 6f 65 73 20 73 | 65 65 6d 0a 58 74 6f 20 |2 does s|eem.Xto |
|00005ee0| 6c 65 61 76 65 20 73 6f | 6d 65 20 72 75 62 62 69 |leave so|me rubbi|
|00005ef0| 73 68 20 6f 6e 20 74 68 | 65 20 73 74 61 63 6b 2e |sh on th|e stack.|
|00005f00| 20 49 20 73 75 73 70 65 | 63 74 20 74 68 69 73 20 | I suspe|ct this |
|00005f10| 69 73 20 61 20 70 72 69 | 6e 74 65 72 20 62 75 66 |is a pri|nter buf|
|00005f20| 66 65 72 0a 58 73 6f 20 | 69 74 20 6d 61 79 20 6e |fer.Xso |it may n|
|00005f30| 6f 74 20 62 65 20 61 20 | 62 75 67 2e 20 54 68 65 |ot be a |bug. The|
|00005f40| 72 65 20 77 65 72 65 20 | 61 20 6e 75 6d 62 65 72 |re were |a number|
|00005f50| 20 6f 66 20 77 61 72 6e | 69 6e 67 20 6d 65 73 73 | of warn|ing mess|
|00005f60| 61 67 65 73 20 66 72 6f | 6d 20 74 68 65 0a 58 63 |ages fro|m the.Xc|
|00005f70| 6f 6d 70 69 6c 65 72 20 | 61 62 6f 75 74 20 6d 79 |ompiler |about my|
|00005f80| 20 22 75 6e 63 6f 6e 73 | 74 61 6e 74 69 6e 67 22 | "uncons|tanting"|
|00005f90| 20 63 6f 6e 73 74 61 6e | 74 73 3b 20 62 75 74 20 | constan|ts; but |
|00005fa0| 49 20 74 68 69 6e 6b 20 | 74 68 69 73 20 77 61 73 |I think |this was|
|00005fb0| 20 6a 75 73 74 0a 58 67 | 6e 75 20 62 65 69 6e 67 | just.Xg|nu being|
|00005fc0| 20 6f 76 65 72 2d 73 65 | 6e 73 69 74 69 76 65 2e | over-se|nsitive.|
|00005fd0| 20 47 6e 75 20 32 2e 33 | 2e 32 20 73 65 65 6d 73 | Gnu 2.3|.2 seems|
|00005fe0| 20 74 6f 20 72 65 70 6f | 72 74 20 69 6e 74 65 72 | to repo|rt inter|
|00005ff0| 6e 61 6c 20 65 72 72 6f | 72 73 20 2d 0a 58 74 68 |nal erro|rs -.Xth|
|00006000| 65 73 65 20 64 6f 6e 27 | 74 20 73 65 65 6d 20 74 |ese don'|t seem t|
|00006010| 6f 20 62 65 20 63 6f 6e | 73 69 73 74 65 6e 74 2c |o be con|sistent,|
|00006020| 20 64 69 66 66 65 72 65 | 6e 74 20 75 73 65 72 73 | differe|nt users|
|00006030| 20 72 65 70 6f 72 74 20 | 64 69 66 66 65 72 65 6e | report |differen|
|00006040| 74 0a 58 65 78 70 65 72 | 69 65 6e 63 65 73 3b 20 |t.Xexper|iences; |
|00006050| 49 20 73 75 67 67 65 73 | 74 2c 20 69 66 20 70 6f |I sugges|t, if po|
|00006060| 73 73 69 62 6c 65 2c 20 | 79 6f 75 20 73 74 69 63 |ssible, |you stic|
|00006070| 6b 20 74 6f 20 76 65 72 | 73 69 6f 6e 20 32 2e 32 |k to ver|sion 2.2|
|00006080| 20 75 6e 74 69 6c 20 32 | 2e 33 0a 58 73 6f 72 74 | until 2|.3.Xsort|
|00006090| 73 20 69 74 73 65 6c 66 | 20 6f 75 74 2e 0a 58 0a |s itself| out..X.|
|000060a0| 58 4a 50 49 3a 20 54 68 | 65 69 72 20 73 65 63 6f |XJPI: Th|eir seco|
|000060b0| 6e 64 20 72 65 6c 65 61 | 73 65 20 77 6f 72 6b 65 |nd relea|se worke|
|000060c0| 64 20 6f 6e 20 61 20 70 | 72 65 76 69 6f 75 73 20 |d on a p|revious |
|000060d0| 76 65 72 73 69 6f 6e 20 | 6f 66 20 74 68 69 73 20 |version |of this |
|000060e0| 70 61 63 6b 61 67 65 0a | 58 70 72 6f 76 69 64 65 |package.|Xprovide|
|000060f0| 64 20 79 6f 75 20 64 69 | 73 61 62 6c 65 64 20 74 |d you di|sabled t|
|00006100| 68 65 20 73 6d 61 72 74 | 20 6c 69 6e 6b 20 6f 70 |he smart| link op|
|00006110| 74 69 6f 6e 20 2d 20 69 | 74 20 69 73 6e 27 74 20 |tion - i|t isn't |
|00006120| 73 6d 61 72 74 20 65 6e | 6f 75 67 68 2e 20 49 0a |smart en|ough. I.|
|00006130| 58 68 61 76 65 6e 27 74 | 20 74 65 73 74 65 64 20 |Xhaven't| tested |
|00006140| 74 68 65 20 6c 61 74 65 | 73 74 20 76 65 72 73 69 |the late|st versi|
|00006150| 6f 6e 20 6f 66 20 74 68 | 69 73 20 70 61 63 6b 61 |on of th|is packa|
|00006160| 67 65 2e 0a 58 0a 58 0a | 58 2d 2d 2d 2d 2d 2d 2d |ge..X.X.|X-------|
|00006170| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00006180| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00006190| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000061a0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000061b0| 2d 2d 2d 2d 0a 58 0a 58 | 45 78 61 6d 70 6c 65 0a |----.X.X|Example.|
|000061c0| 58 3d 3d 3d 3d 3d 3d 3d | 0a 58 0a 58 41 6e 20 65 |X=======|.X.XAn e|
|000061d0| 78 61 6d 70 6c 65 20 69 | 73 20 67 69 76 65 6e 20 |xample i|s given |
|000061e0| 69 6e 20 20 65 78 61 6d | 70 6c 65 2e 63 78 78 20 |in exam|ple.cxx |
|000061f0| 2e 20 20 54 68 69 73 20 | 67 69 76 65 73 20 61 20 |. This |gives a |
|00006200| 73 69 6d 70 6c 65 20 6c | 69 6e 65 61 72 0a 58 72 |simple l|inear.Xr|
|00006210| 65 67 72 65 73 73 69 6f | 6e 20 65 78 61 6d 70 6c |egressio|n exampl|
|00006220| 65 20 75 73 69 6e 67 20 | 66 6f 75 72 20 64 69 66 |e using |four dif|
|00006230| 66 65 72 65 6e 74 20 61 | 6c 67 6f 72 69 74 68 6d |ferent a|lgorithm|
|00006240| 73 2e 20 54 68 65 20 63 | 6f 72 72 65 63 74 20 6f |s. The c|orrect o|
|00006250| 75 74 70 75 74 0a 58 69 | 73 20 67 69 76 65 6e 20 |utput.Xi|s given |
|00006260| 69 6e 20 65 78 61 6d 70 | 6c 65 2e 74 78 74 2e 20 |in examp|le.txt. |
|00006270| 54 68 65 20 70 72 6f 67 | 72 61 6d 20 63 61 72 72 |The prog|ram carr|
|00006280| 69 65 73 20 6f 75 74 20 | 61 20 63 68 65 63 6b 20 |ies out |a check |
|00006290| 74 68 61 74 20 6e 6f 20 | 6d 65 6d 6f 72 79 0a 58 |that no |memory.X|
|000062a0| 69 73 20 6c 65 66 74 20 | 61 6c 6c 6f 63 61 74 65 |is left |allocate|
|000062b0| 64 20 6f 6e 20 74 68 65 | 20 68 65 61 70 20 77 68 |d on the| heap wh|
|000062c0| 65 6e 20 69 74 20 74 65 | 72 6d 69 6e 61 74 65 73 |en it te|rminates|
|000062d0| 2e 0a 58 0a 58 54 68 65 | 20 66 69 6c 65 20 20 65 |..X.XThe| file e|
|000062e0| 78 61 6d 70 6c 65 2e 6d | 61 6b 20 69 73 20 61 20 |xample.m|ak is a |
|000062f0| 6d 61 6b 65 20 66 69 6c | 65 20 66 6f 72 20 63 6f |make fil|e for co|
|00006300| 6d 70 69 6c 69 6e 67 20 | 65 78 61 6d 70 6c 65 2e |mpiling |example.|
|00006310| 63 78 78 20 75 6e 64 65 | 72 20 67 6e 75 20 67 2b |cxx unde|r gnu g+|
|00006320| 2b 2e 0a 58 55 73 65 20 | 74 68 65 20 67 6e 75 20 |+..XUse |the gnu |
|00006330| 6d 61 6b 65 20 66 61 63 | 69 6c 69 74 79 2e 20 59 |make fac|ility. Y|
|00006340| 6f 75 20 63 61 6e 20 70 | 72 6f 62 61 62 6c 79 20 |ou can p|robably |
|00006350| 61 64 61 70 74 20 69 74 | 20 66 6f 72 20 74 68 65 |adapt it| for the|
|00006360| 20 63 6f 6d 70 69 6c 65 | 72 20 79 6f 75 20 61 72 | compile|r you ar|
|00006370| 65 0a 58 75 73 69 6e 67 | 2e 20 49 20 61 6c 73 6f |e.Xusing|. I also|
|00006380| 20 69 6e 63 6c 75 64 65 | 20 74 68 65 20 6d 61 6b | include| the mak|
|00006390| 65 20 66 69 6c 65 73 20 | 66 6f 72 20 5a 6f 72 74 |e files |for Zort|
|000063a0| 65 63 68 2c 20 42 6f 72 | 6c 61 6e 64 20 61 6e 64 |ech, Bor|land and|
|000063b0| 20 4d 69 63 72 6f 73 6f | 66 74 20 43 20 2d 0a 58 | Microso|ft C -.X|
|000063c0| 73 65 65 20 74 68 65 20 | 6c 69 73 74 20 6f 66 20 |see the |list of |
|000063d0| 66 69 6c 65 73 2e 20 55 | 73 65 20 61 20 63 6f 6d |files. U|se a com|
|000063e0| 6d 61 6e 64 20 6c 69 6b | 65 0a 58 0a 58 67 6d 61 |mand lik|e.X.Xgma|
|000063f0| 6b 65 20 2d 66 20 65 78 | 61 6d 70 6c 65 2e 6d 61 |ke -f ex|ample.ma|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
