home *** CD-ROM | disk | FTP | other *** search
view JSON data
|
view as text
|
open on a Mac
|
open on a PC
This file was processed as: SHell self-extracting ARchive
(archive/shar).
You can browse this item here: part05
| 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 46 72 |.sources|.misc.Fr|
|00000020| 6f 6d 3a 20 72 6f 62 65 | 72 74 40 61 6d 2e 64 73 |om: robe|rt@am.ds|
|00000030| 69 72 2e 67 6f 76 74 2e | 6e 7a 20 28 52 6f 62 65 |ir.govt.|nz (Robe|
|00000040| 72 74 20 44 61 76 69 65 | 73 29 0a 53 75 62 6a 65 |rt Davie|s).Subje|
|00000050| 63 74 3a 20 20 76 32 36 | 69 30 39 31 3a 20 20 6e |ct: v26|i091: n|
|00000060| 65 77 6d 61 74 30 34 20 | 2d 20 41 20 6d 61 74 72 |ewmat04 |- A matr|
|00000070| 69 78 20 70 61 63 6b 61 | 67 65 20 69 6e 20 43 2b |ix packa|ge in C+|
|00000080| 2b 2c 20 50 61 72 74 30 | 35 2f 30 35 0a 4d 65 73 |+, Part0|5/05.Mes|
|00000090| 73 61 67 65 2d 49 44 3a | 20 3c 31 39 39 31 4e 6f |sage-ID:| <1991No|
|000000a0| 76 33 30 2e 32 30 35 34 | 32 33 2e 34 30 39 30 40 |v30.2054|23.4090@|
|000000b0| 73 70 61 72 6b 79 2e 69 | 6d 64 2e 73 74 65 72 6c |sparky.i|md.sterl|
|000000c0| 69 6e 67 2e 63 6f 6d 3e | 0a 58 2d 4d 64 34 2d 53 |ing.com>|.X-Md4-S|
|000000d0| 69 67 6e 61 74 75 72 65 | 3a 20 30 36 36 35 32 35 |ignature|: 066525|
|000000e0| 65 36 36 35 64 64 30 33 | 35 66 62 66 39 65 62 33 |e665dd03|5fbf9eb3|
|000000f0| 38 30 65 66 39 33 64 31 | 30 30 0a 44 61 74 65 3a |80ef93d1|00.Date:|
|00000100| 20 53 61 74 2c 20 33 30 | 20 4e 6f 76 20 31 39 39 | Sat, 30| Nov 199|
|00000110| 31 20 32 30 3a 35 34 3a | 32 33 20 47 4d 54 0a 41 |1 20:54:|23 GMT.A|
|00000120| 70 70 72 6f 76 65 64 3a | 20 6b 65 6e 74 40 73 70 |pproved:| kent@sp|
|00000130| 61 72 6b 79 2e 69 6d 64 | 2e 73 74 65 72 6c 69 6e |arky.imd|.sterlin|
|00000140| 67 2e 63 6f 6d 0a 0a 53 | 75 62 6d 69 74 74 65 64 |g.com..S|ubmitted|
|00000150| 2d 62 79 3a 20 72 6f 62 | 65 72 74 40 61 6d 2e 64 |-by: rob|ert@am.d|
|00000160| 73 69 72 2e 67 6f 76 74 | 2e 6e 7a 20 28 52 6f 62 |sir.govt|.nz (Rob|
|00000170| 65 72 74 20 44 61 76 69 | 65 73 29 0a 50 6f 73 74 |ert Davi|es).Post|
|00000180| 69 6e 67 2d 6e 75 6d 62 | 65 72 3a 20 56 6f 6c 75 |ing-numb|er: Volu|
|00000190| 6d 65 20 32 36 2c 20 49 | 73 73 75 65 20 39 31 0a |me 26, I|ssue 91.|
|000001a0| 41 72 63 68 69 76 65 2d | 6e 61 6d 65 3a 20 6e 65 |Archive-|name: ne|
|000001b0| 77 6d 61 74 30 34 2f 70 | 61 72 74 30 35 0a 45 6e |wmat04/p|art05.En|
|000001c0| 76 69 72 6f 6e 6d 65 6e | 74 3a 20 43 2b 2b 0a 53 |vironmen|t: C++.S|
|000001d0| 75 70 65 72 73 65 64 65 | 73 3a 20 6e 65 77 6d 61 |upersede|s: newma|
|000001e0| 74 30 32 3a 20 56 6f 6c | 75 6d 65 20 32 31 2c 20 |t02: Vol|ume 21, |
|000001f0| 49 73 73 75 65 20 34 39 | 2d 35 32 0a 0a 23 21 20 |Issue 49|-52..#! |
|00000200| 2f 62 69 6e 2f 73 68 0a | 23 20 54 68 69 73 20 69 |/bin/sh.|# This i|
|00000210| 73 20 61 20 73 68 65 6c | 6c 20 61 72 63 68 69 76 |s a shel|l archiv|
|00000220| 65 2e 20 20 52 65 6d 6f | 76 65 20 61 6e 79 74 68 |e. Remo|ve anyth|
|00000230| 69 6e 67 20 62 65 66 6f | 72 65 20 74 68 69 73 20 |ing befo|re this |
|00000240| 6c 69 6e 65 2c 20 74 68 | 65 6e 20 75 6e 70 61 63 |line, th|en unpac|
|00000250| 6b 0a 23 20 69 74 20 62 | 79 20 73 61 76 69 6e 67 |k.# it b|y saving|
|00000260| 20 69 74 20 69 6e 74 6f | 20 61 20 66 69 6c 65 20 | it into| a file |
|00000270| 61 6e 64 20 74 79 70 69 | 6e 67 20 22 73 68 20 66 |and typi|ng "sh f|
|00000280| 69 6c 65 22 2e 20 20 54 | 6f 20 6f 76 65 72 77 72 |ile". T|o overwr|
|00000290| 69 74 65 20 65 78 69 73 | 74 69 6e 67 0a 23 20 66 |ite exis|ting.# f|
|000002a0| 69 6c 65 73 2c 20 74 79 | 70 65 20 22 73 68 20 66 |iles, ty|pe "sh f|
|000002b0| 69 6c 65 20 2d 63 22 2e | 20 20 59 6f 75 20 63 61 |ile -c".| You ca|
|000002c0| 6e 20 61 6c 73 6f 20 66 | 65 65 64 20 74 68 69 73 |n also f|eed this|
|000002d0| 20 61 73 20 73 74 61 6e | 64 61 72 64 20 69 6e 70 | as stan|dard inp|
|000002e0| 75 74 20 76 69 61 0a 23 | 20 75 6e 73 68 61 72 2c |ut via.#| unshar,|
|000002f0| 20 6f 72 20 62 79 20 74 | 79 70 69 6e 67 20 22 73 | or by t|yping "s|
|00000300| 68 20 3c 66 69 6c 65 22 | 2c 20 65 2e 67 2e 2e 20 |h <file"|, e.g.. |
|00000310| 20 49 66 20 74 68 69 73 | 20 61 72 63 68 69 76 65 | If this| archive|
|00000320| 20 69 73 20 63 6f 6d 70 | 6c 65 74 65 2c 20 79 6f | is comp|lete, yo|
|00000330| 75 0a 23 20 77 69 6c 6c | 20 73 65 65 20 74 68 65 |u.# will| see the|
|00000340| 20 66 6f 6c 6c 6f 77 69 | 6e 67 20 6d 65 73 73 61 | followi|ng messa|
|00000350| 67 65 20 61 74 20 74 68 | 65 20 65 6e 64 3a 0a 23 |ge at th|e end:.#|
|00000360| 09 09 22 45 6e 64 20 6f | 66 20 73 68 65 6c 6c 20 |.."End o|f shell |
|00000370| 61 72 63 68 69 76 65 2e | 22 0a 23 20 43 6f 6e 74 |archive.|".# Cont|
|00000380| 65 6e 74 73 3a 20 20 63 | 68 6f 6c 65 73 6b 79 2e |ents: c|holesky.|
|00000390| 63 78 78 20 63 6c 61 73 | 73 2e 74 78 74 20 65 76 |cxx clas|s.txt ev|
|000003a0| 61 6c 75 65 2e 63 78 78 | 20 65 78 61 6d 70 6c 65 |alue.cxx| example|
|000003b0| 2e 63 78 78 20 65 78 61 | 6d 70 6c 65 2e 64 65 70 |.cxx exa|mple.dep|
|000003c0| 0a 23 20 20 20 65 78 61 | 6d 70 6c 65 2e 74 78 74 |.# exa|mple.txt|
|000003d0| 20 66 66 74 2e 63 78 78 | 20 68 68 6f 6c 64 65 72 | fft.cxx| hholder|
|000003e0| 2e 63 78 78 20 6a 61 63 | 6f 62 69 2e 63 78 78 20 |.cxx jac|obi.cxx |
|000003f0| 6e 65 77 6d 61 74 2e 6c | 69 73 20 72 65 61 64 6d |newmat.l|is readm|
|00000400| 65 0a 23 20 20 20 73 6f | 72 74 2e 63 78 78 20 73 |e.# so|rt.cxx s|
|00000410| 75 62 6d 61 74 2e 63 78 | 78 20 73 76 64 2e 63 78 |ubmat.cx|x svd.cx|
|00000420| 78 0a 23 20 57 72 61 70 | 70 65 64 20 62 79 20 72 |x.# Wrap|ped by r|
|00000430| 6f 62 65 72 74 40 6b 65 | 61 20 6f 6e 20 53 61 74 |obert@ke|a on Sat|
|00000440| 20 4e 6f 76 20 33 30 20 | 31 38 3a 32 39 3a 30 37 | Nov 30 |18:29:07|
|00000450| 20 31 39 39 31 0a 50 41 | 54 48 3d 2f 62 69 6e 3a | 1991.PA|TH=/bin:|
|00000460| 2f 75 73 72 2f 62 69 6e | 3a 2f 75 73 72 2f 75 63 |/usr/bin|:/usr/uc|
|00000470| 62 20 3b 20 65 78 70 6f | 72 74 20 50 41 54 48 0a |b ; expo|rt PATH.|
|00000480| 69 66 20 74 65 73 74 20 | 2d 66 20 27 63 68 6f 6c |if test |-f 'chol|
|00000490| 65 73 6b 79 2e 63 78 78 | 27 20 2d 61 20 22 24 7b |esky.cxx|' -a "${|
|000004a0| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|000004b0| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en . ec|ho shar:|
|000004c0| 20 57 69 6c 6c 20 6e 6f | 74 20 63 6c 6f 62 62 65 | Will no|t clobbe|
|000004d0| 72 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |r existi|ng file |
|000004e0| 5c 22 27 63 68 6f 6c 65 | 73 6b 79 2e 63 78 78 27 |\"'chole|sky.cxx'|
|000004f0| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|00000500| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|00000510| 27 63 68 6f 6c 65 73 6b | 79 2e 63 78 78 27 5c 22 |'cholesk|y.cxx'\"|
|00000520| 20 5c 28 39 31 36 20 63 | 68 61 72 61 63 74 65 72 | \(916 c|haracter|
|00000530| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|00000540| 20 3e 27 63 68 6f 6c 65 | 73 6b 79 2e 63 78 78 27 | >'chole|sky.cxx'|
|00000550| 20 3c 3c 27 45 4e 44 5f | 4f 46 5f 46 49 4c 45 27 | <<'END_|OF_FILE'|
|00000560| 0a 58 2f 2f 24 24 20 63 | 68 6f 6c 65 73 6b 79 2e |.X//$$ c|holesky.|
|00000570| 63 78 78 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |cxx | |
|00000580| 20 20 20 20 20 20 20 20 | 63 68 6f 6c 65 73 6b 79 | |cholesky|
|00000590| 20 64 65 63 6f 6d 70 6f | 73 69 74 69 6f 6e 0a 58 | decompo|sition.X|
|000005a0| 0a 58 2f 2f 20 43 6f 70 | 79 72 69 67 68 74 20 28 |.X// Cop|yright (|
|000005b0| 43 29 20 31 39 39 31 3a | 20 52 20 42 20 44 61 76 |C) 1991:| R B Dav|
|000005c0| 69 65 73 20 61 6e 64 20 | 44 53 49 52 0a 58 0a 58 |ies and |DSIR.X.X|
|000005d0| 23 64 65 66 69 6e 65 20 | 57 41 4e 54 5f 4d 41 54 |#define |WANT_MAT|
|000005e0| 48 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 20 22 69 |H.X.X#in|clude "i|
|000005f0| 6e 63 6c 75 64 65 2e 68 | 78 78 22 0a 58 0a 58 23 |nclude.h|xx".X.X#|
|00000600| 69 6e 63 6c 75 64 65 20 | 22 6e 65 77 6d 61 74 2e |include |"newmat.|
|00000610| 68 78 78 22 0a 58 0a 58 | 0a 58 2f 2a 2a 2a 2a 2a |hxx".X.X|.X/*****|
|00000620| 2a 2a 2a 2a 20 43 68 6f | 6c 65 73 6b 79 20 64 65 |**** Cho|lesky de|
|00000630| 63 6f 6d 70 6f 73 69 74 | 69 6f 6e 20 6f 66 20 61 |composit|ion of a|
|00000640| 20 70 6f 73 69 74 69 76 | 65 20 64 65 66 69 6e 69 | positiv|e defini|
|00000650| 74 65 20 6d 61 74 72 69 | 78 20 2a 2a 2a 2a 2a 2a |te matri|x ******|
|00000660| 2a 2a 2a 2a 2a 2a 2a 2f | 0a 58 0a 58 2f 2f 20 53 |*******/|.X.X// S|
|00000670| 75 70 70 6f 73 65 20 53 | 20 69 73 20 73 79 6d 6d |uppose S| is symm|
|00000680| 65 74 72 69 78 20 61 6e | 64 20 70 6f 73 69 74 69 |etrix an|d positi|
|00000690| 76 65 20 64 65 66 69 6e | 69 74 65 2e 20 54 68 65 |ve defin|ite. The|
|000006a0| 6e 20 74 68 65 72 65 20 | 65 78 69 73 74 73 20 61 |n there |exists a|
|000006b0| 20 75 6e 69 71 75 65 0a | 58 2f 2f 20 6c 6f 77 65 | unique.|X// lowe|
|000006c0| 72 20 74 72 69 61 6e 67 | 75 6c 61 72 20 6d 61 74 |r triang|ular mat|
|000006d0| 72 69 78 20 4c 20 73 75 | 63 68 20 74 68 61 74 20 |rix L su|ch that |
|000006e0| 4c 20 4c 2e 74 28 29 20 | 3d 20 53 3b 0a 58 0a 58 |L L.t() |= S;.X.X|
|000006f0| 52 65 74 75 72 6e 4d 61 | 74 72 69 78 20 43 68 6f |ReturnMa|trix Cho|
|00000700| 6c 65 73 6b 79 28 63 6f | 6e 73 74 20 53 79 6d 6d |lesky(co|nst Symm|
|00000710| 65 74 72 69 63 4d 61 74 | 72 69 78 26 20 53 29 0a |etricMat|rix& S).|
|00000720| 58 7b 0a 58 20 20 20 69 | 6e 74 20 6e 72 20 3d 20 |X{.X i|nt nr = |
|00000730| 53 2e 4e 72 6f 77 73 28 | 29 3b 0a 58 20 20 20 4c |S.Nrows(|);.X L|
|00000740| 6f 77 65 72 54 72 69 61 | 6e 67 75 6c 61 72 4d 61 |owerTria|ngularMa|
|00000750| 74 72 69 78 20 54 28 6e | 72 29 3b 0a 58 20 20 20 |trix T(n|r);.X |
|00000760| 72 65 61 6c 2a 20 73 20 | 3d 20 53 2e 53 74 6f 72 |real* s |= S.Stor|
|00000770| 65 28 29 3b 20 72 65 61 | 6c 2a 20 74 20 3d 20 54 |e(); rea|l* t = T|
|00000780| 2e 53 74 6f 72 65 28 29 | 3b 20 72 65 61 6c 2a 20 |.Store()|; real* |
|00000790| 74 6b 20 3d 20 74 3b 0a | 58 0a 58 20 20 20 66 6f |tk = t;.|X.X fo|
|000007a0| 72 20 28 69 6e 74 20 69 | 3d 30 3b 20 69 3c 6e 72 |r (int i|=0; i<nr|
|000007b0| 3b 20 69 2b 2b 29 0a 58 | 20 20 20 7b 0a 58 20 20 |; i++).X| {.X |
|000007c0| 20 20 20 20 72 65 61 6c | 2a 20 74 6a 20 3d 20 74 | real|* tj = t|
|000007d0| 3b 20 72 65 61 6c 2a 20 | 74 69 20 3d 20 74 6b 3b |; real* |ti = tk;|
|000007e0| 0a 58 20 20 20 20 20 20 | 66 6f 72 20 28 69 6e 74 |.X |for (int|
|000007f0| 20 6a 3d 30 3b 20 6a 3c | 3d 69 3b 20 6a 2b 2b 29 | j=0; j<|=i; j++)|
|00000800| 0a 58 20 20 20 20 20 20 | 7b 0a 58 09 20 74 6b 20 |.X |{.X. tk |
|00000810| 3d 20 74 69 3b 20 72 65 | 61 6c 20 73 75 6d 20 3d |= ti; re|al sum =|
|00000820| 20 30 2e 30 3b 20 69 6e | 74 20 6b 20 3d 20 6a 3b | 0.0; in|t k = j;|
|00000830| 0a 58 09 20 77 68 69 6c | 65 20 28 6b 2d 2d 29 20 |.X. whil|e (k--) |
|00000840| 7b 20 73 75 6d 20 2b 3d | 20 2a 74 6a 2b 2b 20 2a |{ sum +=| *tj++ *|
|00000850| 20 2a 74 6b 2b 2b 3b 20 | 7d 0a 58 09 20 69 66 20 | *tk++; |}.X. if |
|00000860| 28 69 3d 3d 6a 29 20 2a | 74 6b 2b 2b 20 3d 20 73 |(i==j) *|tk++ = s|
|00000870| 71 72 74 28 2a 73 2b 2b | 20 2d 20 73 75 6d 29 3b |qrt(*s++| - sum);|
|00000880| 0a 58 20 20 20 20 20 20 | 20 20 20 65 6c 73 65 20 |.X | else |
|00000890| 2a 74 6b 2b 2b 20 3d 20 | 28 2a 73 2b 2b 20 2d 20 |*tk++ = |(*s++ - |
|000008a0| 73 75 6d 29 20 2f 20 2a | 74 6a 2b 2b 3b 0a 58 20 |sum) / *|tj++;.X |
|000008b0| 20 20 20 20 20 7d 0a 58 | 20 20 20 7d 0a 58 20 20 | }.X| }.X |
|000008c0| 20 54 2e 52 65 6c 65 61 | 73 65 28 29 3b 20 20 20 | T.Relea|se(); |
|000008d0| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 6e 6f | | // no|
|000008e0| 74 20 77 61 6e 74 65 64 | 20 69 66 20 72 6f 75 74 |t wanted| if rout|
|000008f0| 69 6e 65 20 61 76 6f 69 | 64 73 20 72 65 74 75 72 |ine avoi|ds retur|
|00000900| 6e 20 69 6e 69 74 0a 58 | 20 20 20 72 65 74 75 72 |n init.X| retur|
|00000910| 6e 20 54 3b 0a 58 7d 0a | 58 0a 45 4e 44 5f 4f 46 |n T;.X}.|X.END_OF|
|00000920| 5f 46 49 4c 45 0a 69 66 | 20 74 65 73 74 20 39 31 |_FILE.if| test 91|
|00000930| 36 20 2d 6e 65 20 60 77 | 63 20 2d 63 20 3c 27 63 |6 -ne `w|c -c <'c|
|00000940| 68 6f 6c 65 73 6b 79 2e | 63 78 78 27 60 3b 20 74 |holesky.|cxx'`; t|
|00000950| 68 65 6e 0a 20 20 20 20 | 65 63 68 6f 20 73 68 61 |hen. |echo sha|
|00000960| 72 3a 20 5c 22 27 63 68 | 6f 6c 65 73 6b 79 2e 63 |r: \"'ch|olesky.c|
|00000970| 78 78 27 5c 22 20 75 6e | 70 61 63 6b 65 64 20 77 |xx'\" un|packed w|
|00000980| 69 74 68 20 77 72 6f 6e | 67 20 73 69 7a 65 21 0a |ith wron|g size!.|
|00000990| 66 69 0a 23 20 65 6e 64 | 20 6f 66 20 27 63 68 6f |fi.# end| of 'cho|
|000009a0| 6c 65 73 6b 79 2e 63 78 | 78 27 0a 66 69 0a 69 66 |lesky.cx|x'.fi.if|
|000009b0| 20 74 65 73 74 20 2d 66 | 20 27 63 6c 61 73 73 2e | test -f| 'class.|
|000009c0| 74 78 74 27 20 2d 61 20 | 22 24 7b 31 7d 22 20 21 |txt' -a |"${1}" !|
|000009d0| 3d 20 22 2d 63 22 20 3b | 20 74 68 65 6e 20 0a 20 |= "-c" ;| then . |
|000009e0| 20 65 63 68 6f 20 73 68 | 61 72 3a 20 57 69 6c 6c | echo sh|ar: Will|
|000009f0| 20 6e 6f 74 20 63 6c 6f | 62 62 65 72 20 65 78 69 | not clo|bber exi|
|00000a00| 73 74 69 6e 67 20 66 69 | 6c 65 20 5c 22 27 63 6c |sting fi|le \"'cl|
|00000a10| 61 73 73 2e 74 78 74 27 | 5c 22 0a 65 6c 73 65 0a |ass.txt'|\".else.|
|00000a20| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|00000a30| 63 74 69 6e 67 20 5c 22 | 27 63 6c 61 73 73 2e 74 |cting \"|'class.t|
|00000a40| 78 74 27 5c 22 20 5c 28 | 32 33 39 31 20 63 68 61 |xt'\" \(|2391 cha|
|00000a50| 72 61 63 74 65 72 73 5c | 29 0a 73 65 64 20 22 73 |racters\|).sed "s|
|00000a60| 2f 5e 58 2f 2f 22 20 3e | 27 63 6c 61 73 73 2e 74 |/^X//" >|'class.t|
|00000a70| 78 74 27 20 3c 3c 27 45 | 4e 44 5f 4f 46 5f 46 49 |xt' <<'E|ND_OF_FI|
|00000a80| 4c 45 27 0a 58 2f 2f 24 | 24 20 63 6c 61 73 73 2e |LE'.X//$|$ class.|
|00000a90| 74 78 74 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |txt | |
|00000aa0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000ab0| 20 20 20 6c 69 73 74 20 | 6f 66 20 63 6c 61 73 73 | list |of class|
|00000ac0| 65 73 0a 58 0a 58 0a 58 | 4c 69 73 74 20 6f 66 20 |es.X.X.X|List of |
|00000ad0| 63 6c 61 73 73 65 73 20 | 69 6e 20 6e 65 77 6d 61 |classes |in newma|
|00000ae0| 74 30 34 0a 58 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |t04.X===|========|
|00000af0| 3d 3d 3d 3d 3d 3d 3d 3d | 3d 3d 3d 3d 3d 3d 3d 3d |========|========|
|00000b00| 0a 58 0a 58 0a 58 42 4f | 4f 4c 20 20 20 20 20 20 |.X.X.XBO|OL |
|00000b10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000b20| 20 20 20 20 20 20 42 6f | 6f 6c 65 61 6e 20 63 6c | Bo|olean cl|
|00000b30| 61 73 73 0a 58 0a 58 43 | 6f 6e 74 72 6f 6c 57 6f |ass.X.XC|ontrolWo|
|00000b40| 72 64 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |rd | |
|00000b50| 20 20 20 20 20 20 20 69 | 6e 74 20 69 6e 74 65 72 | i|nt inter|
|00000b60| 70 72 65 74 65 64 20 61 | 73 20 73 65 71 75 65 6e |preted a|s sequen|
|00000b70| 63 65 20 6f 66 20 62 69 | 74 73 0a 58 0a 58 4c 6f |ce of bi|ts.X.XLo|
|00000b80| 67 41 6e 64 53 69 67 6e | 20 20 20 20 20 20 20 20 |gAndSign| |
|00000b90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 72 65 | | re|
|00000ba0| 74 75 72 6e 20 66 72 6f | 6d 20 64 65 74 65 72 6d |turn fro|m determ|
|00000bb0| 69 6e 61 6e 74 20 72 6f | 75 74 69 6e 65 0a 58 0a |inant ro|utine.X.|
|00000bc0| 58 42 61 73 65 4d 61 74 | 72 69 78 20 20 20 20 20 |XBaseMat|rix |
|00000bd0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000be0| 20 62 61 73 65 20 6f 66 | 20 61 6c 6c 20 6d 61 74 | base of| all mat|
|00000bf0| 72 69 78 20 63 6c 61 73 | 73 65 73 0a 58 20 7c 0a |rix clas|ses.X |.|
|00000c00| 58 20 2b 2d 47 65 6e 65 | 72 61 6c 4d 61 74 72 69 |X +-Gene|ralMatri|
|00000c10| 78 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |x | |
|00000c20| 20 62 61 73 65 20 6f 66 | 20 6d 61 74 72 69 78 20 | base of| matrix |
|00000c30| 63 6c 61 73 73 65 73 20 | 61 63 63 65 73 73 69 62 |classes |accessib|
|00000c40| 6c 65 20 74 6f 20 61 20 | 75 73 65 72 0a 58 20 7c |le to a |user.X ||
|00000c50| 20 20 7c 0a 58 20 7c 20 | 20 2b 2d 4d 61 74 72 69 | |.X | | +-Matri|
|00000c60| 78 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |x | |
|00000c70| 20 20 20 20 20 72 65 63 | 74 61 6e 67 75 6c 61 72 | rec|tangular|
|00000c80| 20 6d 61 74 72 69 78 0a | 58 20 7c 20 20 7c 20 20 | matrix.|X | | |
|00000c90| 7c 0a 58 20 7c 20 20 7c | 20 20 2b 2d 52 6f 77 56 ||.X | || +-RowV|
|00000ca0| 65 63 74 6f 72 0a 58 20 | 7c 20 20 7c 20 20 7c 0a |ector.X || | |.|
|00000cb0| 58 20 7c 20 20 7c 20 20 | 2b 2d 43 6f 6c 75 6d 6e |X | | |+-Column|
|00000cc0| 56 65 63 74 6f 72 0a 58 | 20 7c 20 20 7c 20 20 7c |Vector.X| | | ||
|00000cd0| 0a 58 20 7c 20 20 7c 20 | 20 2b 2d 6e 72 69 63 4d |.X | | | +-nricM|
|00000ce0| 61 74 72 69 78 20 20 20 | 20 20 20 20 20 20 20 20 |atrix | |
|00000cf0| 20 20 66 6f 72 20 63 6f | 6e 6e 65 63 74 69 6e 67 | for co|nnecting|
|00000d00| 20 77 69 74 68 20 22 4e | 75 6d 65 72 69 63 61 6c | with "N|umerical|
|00000d10| 20 52 65 63 69 70 65 73 | 20 69 6e 20 43 22 0a 58 | Recipes| in C".X|
|00000d20| 20 7c 20 20 7c 0a 58 20 | 7c 20 20 2b 2d 53 79 6d | | |.X || +-Sym|
|00000d30| 6d 65 74 72 69 63 4d 61 | 74 72 69 78 0a 58 20 7c |metricMa|trix.X ||
|00000d40| 20 20 7c 0a 58 20 7c 20 | 20 2b 2d 55 70 70 65 72 | |.X | | +-Upper|
|00000d50| 54 72 69 61 6e 67 75 6c | 61 72 4d 61 74 72 69 78 |Triangul|arMatrix|
|00000d60| 0a 58 20 7c 20 20 7c 0a | 58 20 7c 20 20 2b 2d 4c |.X | |.|X | +-L|
|00000d70| 6f 77 65 72 54 72 69 61 | 6e 67 75 6c 61 72 4d 61 |owerTria|ngularMa|
|00000d80| 74 72 69 78 0a 58 20 7c | 20 20 7c 0a 58 20 7c 20 |trix.X || |.X | |
|00000d90| 20 2b 2d 44 69 61 67 6f | 6e 61 6c 4d 61 74 72 69 | +-Diago|nalMatri|
|00000da0| 78 0a 58 20 7c 20 20 7c | 0a 58 20 7c 20 20 2b 2d |x.X | ||.X | +-|
|00000db0| 43 72 6f 75 74 4d 61 74 | 72 69 78 20 20 20 20 20 |CroutMat|rix |
|00000dc0| 20 20 20 20 20 20 20 20 | 20 20 43 72 6f 75 74 20 | | Crout |
|00000dd0| 64 65 63 6f 6d 70 6f 73 | 69 74 69 6f 6e 20 6f 66 |decompos|ition of|
|00000de0| 20 61 20 72 65 63 74 61 | 6e 67 75 6c 61 72 20 6d | a recta|ngular m|
|00000df0| 61 74 72 69 78 0a 58 20 | 7c 20 20 0a 58 20 2b 2d |atrix.X || .X +-|
|00000e00| 52 65 74 75 72 6e 4d 61 | 74 72 69 78 20 20 20 20 |ReturnMa|trix |
|00000e10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 6d 61 74 | | mat|
|00000e20| 72 69 78 20 74 6f 20 62 | 65 20 72 65 74 75 72 6e |rix to b|e return|
|00000e30| 65 64 20 66 72 6f 6d 20 | 66 75 6e 63 74 69 6f 6e |ed from |function|
|00000e40| 0a 58 20 7c 0a 58 2e 7c | 2e 2e 2e 2e 2e 2e 2e 2e |.X |.X.||........|
|00000e50| 2e 2e 2e 2e 2e 2e 2e 2e | 2e 2e 2e 2e 2e 2e 2e 2e |........|........|
|00000e60| 2e 2e 2e 2e 2e 2e 43 6c | 61 73 73 65 73 20 62 65 |......Cl|asses be|
|00000e70| 6c 6f 77 20 68 65 72 65 | 20 61 72 65 20 6e 6f 74 |low here| are not|
|00000e80| 20 69 6e 74 65 6e 64 65 | 64 20 66 6f 72 0a 58 20 | intende|d for.X |
|00000e90| 7c 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 || | |
|00000ea0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 75 | | u|
|00000eb0| 73 65 72 20 61 63 63 65 | 73 73 2e 0a 58 20 7c 0a |ser acce|ss..X |.|
|00000ec0| 58 20 2b 2d 41 64 64 65 | 64 4d 61 74 72 69 78 20 |X +-Adde|dMatrix |
|00000ed0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000ee0| 20 63 6c 61 73 73 20 6f | 66 20 41 20 2b 20 42 0a | class o|f A + B.|
|00000ef0| 58 20 7c 20 20 7c 0a 58 | 20 7c 20 20 2b 2d 4d 75 |X | |.X| | +-Mu|
|00000f00| 6c 74 69 70 6c 69 65 64 | 4d 61 74 72 69 78 20 20 |ltiplied|Matrix |
|00000f10| 20 20 20 20 20 20 20 20 | 63 6c 61 73 73 20 6f 66 | |class of|
|00000f20| 20 41 20 2a 20 42 0a 58 | 20 7c 20 20 7c 0a 58 20 | A * B.X| | |.X |
|00000f30| 7c 20 20 2b 2d 53 6f 6c | 76 65 64 4d 61 74 72 69 || +-Sol|vedMatri|
|00000f40| 78 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 63 |x | c|
|00000f50| 6c 61 73 73 20 6f 66 20 | 41 2e 69 28 29 20 2a 20 |lass of |A.i() * |
|00000f60| 42 0a 58 20 7c 20 20 7c | 0a 58 20 7c 20 20 2b 2d |B.X | ||.X | +-|
|00000f70| 53 75 62 74 72 61 63 74 | 65 64 4d 61 74 72 69 78 |Subtract|edMatrix|
|00000f80| 20 20 20 20 20 20 20 20 | 20 20 63 6c 61 73 73 20 | | class |
|00000f90| 6f 66 20 41 20 2d 20 42 | 0a 58 20 7c 0a 58 20 2b |of A - B|.X |.X +|
|00000fa0| 2d 53 68 69 66 74 65 64 | 4d 61 74 72 69 78 20 20 |-Shifted|Matrix |
|00000fb0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 63 6c | | cl|
|00000fc0| 61 73 73 20 6f 66 20 41 | 20 2b 20 66 0a 58 20 7c |ass of A| + f.X ||
|00000fd0| 20 20 7c 0a 58 20 7c 20 | 20 2b 2d 53 63 61 6c 65 | |.X | | +-Scale|
|00000fe0| 64 4d 61 74 72 69 78 20 | 20 20 20 20 20 20 20 20 |dMatrix | |
|00000ff0| 20 20 20 20 20 63 6c 61 | 73 73 20 6f 66 20 41 20 | cla|ss of A |
|00001000| 2a 20 66 0a 58 20 7c 20 | 20 0a 58 20 2b 2d 4e 65 |* f.X | | .X +-Ne|
|00001010| 67 61 74 65 64 4d 61 74 | 72 69 78 20 20 20 20 20 |gatedMat|rix |
|00001020| 20 20 20 20 20 20 20 20 | 20 20 20 63 6c 61 73 73 | | class|
|00001030| 20 6f 66 20 2d 41 0a 58 | 20 7c 20 20 7c 0a 58 20 | of -A.X| | |.X |
|00001040| 7c 20 20 2b 2d 54 72 61 | 6e 73 70 6f 73 65 64 4d || +-Tra|nsposedM|
|00001050| 61 74 72 69 78 20 20 20 | 20 20 20 20 20 20 20 63 |atrix | c|
|00001060| 6c 61 73 73 20 6f 66 20 | 41 2e 74 28 29 0a 58 20 |lass of |A.t().X |
|00001070| 7c 20 20 7c 0a 58 20 7c | 20 20 2b 2d 49 6e 76 65 || |.X || +-Inve|
|00001080| 72 74 65 64 4d 61 74 72 | 69 78 20 20 20 20 20 20 |rtedMatr|ix |
|00001090| 20 20 20 20 20 20 63 6c | 61 73 73 20 6f 66 20 41 | cl|ass of A|
|000010a0| 2e 69 28 29 0a 58 20 7c | 20 20 7c 0a 58 20 7c 20 |.i().X || |.X | |
|000010b0| 20 2b 2d 52 6f 77 65 64 | 4d 61 74 72 69 78 20 20 | +-Rowed|Matrix |
|000010c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 63 6c 61 | | cla|
|000010d0| 73 73 20 6f 66 20 41 2e | 43 6f 70 79 54 6f 52 6f |ss of A.|CopyToRo|
|000010e0| 77 28 29 0a 58 20 7c 20 | 20 7c 0a 58 20 7c 20 20 |w().X | | |.X | |
|000010f0| 2b 2d 43 6f 6c 65 64 4d | 61 74 72 69 78 20 20 20 |+-ColedM|atrix |
|00001100| 20 20 20 20 20 20 20 20 | 20 20 20 20 63 6c 61 73 | | clas|
|00001110| 73 20 6f 66 20 41 2e 43 | 6f 70 79 54 6f 43 6f 6c |s of A.C|opyToCol|
|00001120| 75 6d 6e 28 29 0a 58 20 | 7c 20 20 7c 0a 58 20 7c |umn().X || |.X ||
|00001130| 20 20 2b 2d 44 69 61 67 | 65 64 4d 61 74 72 69 78 | +-Diag|edMatrix|
|00001140| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 63 6c | | cl|
|00001150| 61 73 73 20 6f 66 20 41 | 2e 43 6f 70 79 54 6f 44 |ass of A|.CopyToD|
|00001160| 69 61 67 6f 6e 61 6c 28 | 29 0a 58 20 7c 20 20 7c |iagonal(|).X | ||
|00001170| 0a 58 20 7c 20 20 2b 2d | 4d 61 74 65 64 4d 61 74 |.X | +-|MatedMat|
|00001180| 72 69 78 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |rix | |
|00001190| 20 20 63 6c 61 73 73 20 | 6f 66 20 41 2e 43 6f 70 | class |of A.Cop|
|000011a0| 79 74 6f 4d 61 74 72 69 | 78 28 72 2c 63 29 0a 58 |ytoMatri|x(r,c).X|
|000011b0| 20 7c 20 20 7c 0a 58 20 | 7c 20 20 2b 2d 47 65 74 | | |.X || +-Get|
|000011c0| 53 75 62 4d 61 74 72 69 | 78 20 20 20 20 20 20 20 |SubMatri|x |
|000011d0| 20 20 20 20 20 20 20 63 | 6c 61 73 73 20 6f 66 20 | c|lass of |
|000011e0| 41 2e 53 75 62 4d 61 74 | 72 69 78 28 66 72 2c 6c |A.SubMat|rix(fr,l|
|000011f0| 72 2c 66 63 2c 6c 63 29 | 0a 58 20 7c 20 20 0a 58 |r,fc,lc)|.X | .X|
|00001200| 20 2b 2d 43 6f 6e 73 74 | 4d 61 74 72 69 78 20 20 | +-Const|Matrix |
|00001210| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001220| 63 6c 61 73 73 20 6f 66 | 20 41 2e 63 28 29 0a 58 |class of| A.c().X|
|00001230| 0a 58 46 6c 6f 61 74 69 | 6e 67 50 6f 69 6e 74 50 |.XFloati|ngPointP|
|00001240| 72 65 63 69 73 69 6f 6e | 20 20 20 20 20 20 20 20 |recision| |
|00001250| 20 20 63 6f 6e 74 61 69 | 6e 73 20 66 6c 6f 61 74 | contai|ns float|
|00001260| 69 6e 67 20 70 6f 69 6e | 74 20 70 72 65 63 69 73 |ing poin|t precis|
|00001270| 69 6f 6e 20 63 6f 6e 73 | 74 61 6e 74 73 0a 58 0a |ion cons|tants.X.|
|00001280| 58 53 79 6d 6d 65 74 72 | 69 63 45 69 67 65 6e 41 |XSymmetr|icEigenA|
|00001290| 6e 61 6c 79 73 69 73 20 | 20 20 20 20 20 20 20 20 |nalysis | |
|000012a0| 20 6e 6f 74 20 75 73 65 | 64 20 79 65 74 0a 58 0a | not use|d yet.X.|
|000012b0| 58 4d 61 74 72 69 78 52 | 6f 77 63 6f 6c 20 20 20 |XMatrixR|owcol |
|000012c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000012d0| 20 72 6f 77 20 61 6e 64 | 20 63 6f 6c 75 6d 6e 20 | row and| column |
|000012e0| 63 6c 61 73 73 65 73 20 | 66 6f 72 20 67 65 6e 65 |classes |for gene|
|000012f0| 72 61 6c 20 6d 61 74 72 | 69 78 0a 58 20 7c 0a 58 |ral matr|ix.X |.X|
|00001300| 20 2b 2d 4d 61 74 72 69 | 78 52 6f 77 0a 58 20 7c | +-Matri|xRow.X ||
|00001310| 0a 58 20 2b 2d 4d 61 74 | 72 69 78 43 6f 6c 0a 58 |.X +-Mat|rixCol.X|
|00001320| 0a 58 52 65 63 74 4d 61 | 74 72 69 78 52 6f 77 43 |.XRectMa|trixRowC|
|00001330| 6f 6c 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |ol | |
|00001340| 20 20 72 6f 77 20 61 6e | 64 20 63 6f 6c 75 6d 6e | row an|d column|
|00001350| 20 63 6c 61 73 73 65 73 | 20 66 6f 72 20 72 65 63 | classes| for rec|
|00001360| 74 61 6e 67 75 6c 61 72 | 20 6d 61 74 72 69 78 0a |tangular| matrix.|
|00001370| 58 20 7c 0a 58 20 2b 2d | 52 65 63 74 4d 61 74 72 |X |.X +-|RectMatr|
|00001380| 69 78 52 6f 77 0a 58 20 | 7c 0a 58 20 2b 2d 52 65 |ixRow.X ||.X +-Re|
|00001390| 63 74 4d 61 74 72 69 78 | 43 6f 6c 0a 58 20 7c 0a |ctMatrix|Col.X |.|
|000013a0| 58 20 2b 2d 52 65 63 74 | 4d 61 74 72 69 78 44 69 |X +-Rect|MatrixDi|
|000013b0| 61 67 0a 58 0a 58 45 78 | 65 43 6f 75 6e 74 65 72 |ag.X.XEx|eCounter|
|000013c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000013d0| 20 20 20 20 20 20 63 6f | 75 6e 74 73 20 6e 75 6d | co|unts num|
|000013e0| 62 65 72 20 6f 66 20 65 | 6e 74 72 69 65 73 20 74 |ber of e|ntries t|
|000013f0| 6f 20 61 20 70 69 65 63 | 65 20 6f 66 20 63 6f 64 |o a piec|e of cod|
|00001400| 65 0a 58 0a 58 4d 61 74 | 72 69 78 54 79 70 65 20 |e.X.XMat|rixType |
|00001410| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001420| 20 20 20 20 20 74 79 70 | 65 20 6f 66 20 61 20 6d | typ|e of a m|
|00001430| 61 74 72 69 78 0a 45 4e | 44 5f 4f 46 5f 46 49 4c |atrix.EN|D_OF_FIL|
|00001440| 45 0a 69 66 20 74 65 73 | 74 20 32 33 39 31 20 2d |E.if tes|t 2391 -|
|00001450| 6e 65 20 60 77 63 20 2d | 63 20 3c 27 63 6c 61 73 |ne `wc -|c <'clas|
|00001460| 73 2e 74 78 74 27 60 3b | 20 74 68 65 6e 0a 20 20 |s.txt'`;| then. |
|00001470| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 27 | echo s|har: \"'|
|00001480| 63 6c 61 73 73 2e 74 78 | 74 27 5c 22 20 75 6e 70 |class.tx|t'\" unp|
|00001490| 61 63 6b 65 64 20 77 69 | 74 68 20 77 72 6f 6e 67 |acked wi|th wrong|
|000014a0| 20 73 69 7a 65 21 0a 66 | 69 0a 23 20 65 6e 64 20 | size!.f|i.# end |
|000014b0| 6f 66 20 27 63 6c 61 73 | 73 2e 74 78 74 27 0a 66 |of 'clas|s.txt'.f|
|000014c0| 69 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 65 76 |i.if tes|t -f 'ev|
|000014d0| 61 6c 75 65 2e 63 78 78 | 27 20 2d 61 20 22 24 7b |alue.cxx|' -a "${|
|000014e0| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|000014f0| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en . ec|ho shar:|
|00001500| 20 57 69 6c 6c 20 6e 6f | 74 20 63 6c 6f 62 62 65 | Will no|t clobbe|
|00001510| 72 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |r existi|ng file |
|00001520| 5c 22 27 65 76 61 6c 75 | 65 2e 63 78 78 27 5c 22 |\"'evalu|e.cxx'\"|
|00001530| 0a 65 6c 73 65 0a 65 63 | 68 6f 20 73 68 61 72 3a |.else.ec|ho shar:|
|00001540| 20 45 78 74 72 61 63 74 | 69 6e 67 20 5c 22 27 65 | Extract|ing \"'e|
|00001550| 76 61 6c 75 65 2e 63 78 | 78 27 5c 22 20 5c 28 37 |value.cx|x'\" \(7|
|00001560| 35 37 32 20 63 68 61 72 | 61 63 74 65 72 73 5c 29 |572 char|acters\)|
|00001570| 0a 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 27 |.sed "s/|^X//" >'|
|00001580| 65 76 61 6c 75 65 2e 63 | 78 78 27 20 3c 3c 27 45 |evalue.c|xx' <<'E|
|00001590| 4e 44 5f 4f 46 5f 46 49 | 4c 45 27 0a 58 2f 2f 24 |ND_OF_FI|LE'.X//$|
|000015a0| 24 65 76 61 6c 75 65 2e | 63 78 78 20 20 20 20 20 |$evalue.|cxx |
|000015b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000015c0| 20 20 20 20 20 20 65 69 | 67 65 6e 2d 76 61 6c 75 | ei|gen-valu|
|000015d0| 65 20 64 65 63 6f 6d 70 | 6f 73 69 74 69 6f 6e 0a |e decomp|osition.|
|000015e0| 58 0a 58 2f 2f 20 43 6f | 70 79 72 69 67 68 74 20 |X.X// Co|pyright |
|000015f0| 28 43 29 20 31 39 39 31 | 3a 20 52 20 42 20 44 61 |(C) 1991|: R B Da|
|00001600| 76 69 65 73 20 61 6e 64 | 20 44 53 49 52 0a 58 0a |vies and| DSIR.X.|
|00001610| 58 23 64 65 66 69 6e 65 | 20 57 41 4e 54 5f 4d 41 |X#define| WANT_MA|
|00001620| 54 48 0a 58 0a 58 23 69 | 6e 63 6c 75 64 65 20 22 |TH.X.X#i|nclude "|
|00001630| 69 6e 63 6c 75 64 65 2e | 68 78 78 22 0a 58 23 69 |include.|hxx".X#i|
|00001640| 6e 63 6c 75 64 65 20 22 | 6e 65 77 6d 61 74 2e 68 |nclude "|newmat.h|
|00001650| 78 78 22 0a 58 23 69 6e | 63 6c 75 64 65 20 22 6e |xx".X#in|clude "n|
|00001660| 65 77 6d 61 74 72 6d 2e | 68 78 78 22 0a 58 23 69 |ewmatrm.|hxx".X#i|
|00001670| 6e 63 6c 75 64 65 20 22 | 70 72 65 63 69 73 69 6f |nclude "|precisio|
|00001680| 2e 68 78 78 22 0a 58 0a | 58 0a 58 73 74 61 74 69 |.hxx".X.|X.Xstati|
|00001690| 63 20 76 6f 69 64 20 74 | 72 65 64 32 28 63 6f 6e |c void t|red2(con|
|000016a0| 73 74 20 53 79 6d 6d 65 | 74 72 69 63 4d 61 74 72 |st Symme|tricMatr|
|000016b0| 69 78 26 20 41 2c 20 44 | 69 61 67 6f 6e 61 6c 4d |ix& A, D|iagonalM|
|000016c0| 61 74 72 69 78 26 20 44 | 2c 0a 58 20 20 20 44 69 |atrix& D|,.X Di|
|000016d0| 61 67 6f 6e 61 6c 4d 61 | 74 72 69 78 26 20 45 2c |agonalMa|trix& E,|
|000016e0| 20 4d 61 74 72 69 78 26 | 20 5a 29 0a 58 7b 0a 58 | Matrix&| Z).X{.X|
|000016f0| 20 20 20 72 65 61 6c 20 | 74 6f 6c 20 3d 0a 58 20 | real |tol =.X |
|00001700| 20 20 20 20 20 46 6c 6f | 61 74 69 6e 67 50 6f 69 | Flo|atingPoi|
|00001710| 6e 74 50 72 65 63 69 73 | 69 6f 6e 3a 3a 4d 69 6e |ntPrecis|ion::Min|
|00001720| 69 6d 75 6d 28 29 2f 46 | 6c 6f 61 74 69 6e 67 50 |imum()/F|loatingP|
|00001730| 6f 69 6e 74 50 72 65 63 | 69 73 69 6f 6e 3a 3a 45 |ointPrec|ision::E|
|00001740| 70 73 69 6c 6f 6e 28 29 | 3b 0a 58 20 20 20 69 6e |psilon()|;.X in|
|00001750| 74 20 6e 20 3d 20 41 2e | 4e 72 6f 77 73 28 29 3b |t n = A.|Nrows();|
|00001760| 20 5a 2e 52 65 44 69 6d | 65 6e 73 69 6f 6e 28 6e | Z.ReDim|ension(n|
|00001770| 2c 6e 29 3b 20 5a 2e 49 | 6e 6a 65 63 74 28 41 29 |,n); Z.I|nject(A)|
|00001780| 3b 0a 58 20 20 20 44 2e | 52 65 44 69 6d 65 6e 73 |;.X D.|ReDimens|
|00001790| 69 6f 6e 28 6e 29 3b 20 | 45 2e 52 65 44 69 6d 65 |ion(n); |E.ReDime|
|000017a0| 6e 73 69 6f 6e 28 6e 29 | 3b 20 72 65 61 6c 2a 20 |nsion(n)|; real* |
|000017b0| 7a 20 3d 20 5a 2e 53 74 | 6f 72 65 28 29 3b 0a 58 |z = Z.St|ore();.X|
|000017c0| 0a 58 20 20 20 66 6f 72 | 20 28 69 6e 74 20 69 3d |.X for| (int i=|
|000017d0| 6e 2d 31 3b 20 69 20 3e | 20 30 3b 20 69 2d 2d 29 |n-1; i >| 0; i--)|
|000017e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000017f0| 20 20 20 2f 2f 20 69 3d | 30 20 69 73 20 65 78 63 | // i=|0 is exc|
|00001800| 6c 75 64 65 64 0a 58 20 | 20 20 7b 0a 58 20 20 20 |luded.X | {.X |
|00001810| 20 20 20 72 65 61 6c 20 | 66 20 3d 20 5a 2e 65 6c | real |f = Z.el|
|00001820| 65 6d 65 6e 74 28 69 2c | 69 2d 31 29 3b 20 72 65 |ement(i,|i-1); re|
|00001830| 61 6c 20 67 20 3d 20 30 | 2e 30 3b 0a 58 20 20 20 |al g = 0|.0;.X |
|00001840| 20 20 20 69 6e 74 20 6b | 20 3d 20 69 2d 31 3b 20 | int k| = i-1; |
|00001850| 72 65 61 6c 2a 20 7a 69 | 6b 20 3d 20 7a 20 2b 20 |real* zi|k = z + |
|00001860| 69 2a 6e 3b 0a 58 20 20 | 20 20 20 20 77 68 69 6c |i*n;.X | whil|
|00001870| 65 20 28 6b 2d 2d 29 20 | 67 20 2b 3d 20 73 71 75 |e (k--) |g += squ|
|00001880| 61 72 65 28 2a 7a 69 6b | 2b 2b 29 3b 0a 58 20 20 |are(*zik|++);.X |
|00001890| 20 20 20 20 72 65 61 6c | 20 68 20 3d 20 67 20 2b | real| h = g +|
|000018a0| 20 73 71 75 61 72 65 28 | 66 29 3b 0a 58 20 20 20 | square(|f);.X |
|000018b0| 20 20 20 69 66 20 28 67 | 20 3c 3d 20 74 6f 6c 29 | if (g| <= tol)|
|000018c0| 20 7b 20 45 2e 65 6c 65 | 6d 65 6e 74 28 69 29 20 | { E.ele|ment(i) |
|000018d0| 3d 20 66 3b 20 68 20 3d | 20 30 2e 30 3b 20 7d 0a |= f; h =| 0.0; }.|
|000018e0| 58 20 20 20 20 20 20 65 | 6c 73 65 0a 58 20 20 20 |X e|lse.X |
|000018f0| 20 20 20 7b 0a 58 09 20 | 67 20 3d 20 73 69 67 6e | {.X. |g = sign|
|00001900| 28 2d 73 71 72 74 28 68 | 29 2c 20 66 29 3b 20 45 |(-sqrt(h|), f); E|
|00001910| 2e 65 6c 65 6d 65 6e 74 | 28 69 29 20 3d 20 67 3b |.element|(i) = g;|
|00001920| 20 68 20 2d 3d 20 66 2a | 67 3b 0a 58 09 20 5a 2e | h -= f*|g;.X. Z.|
|00001930| 65 6c 65 6d 65 6e 74 28 | 69 2c 69 2d 31 29 20 3d |element(|i,i-1) =|
|00001940| 20 66 2d 67 3b 20 66 20 | 3d 20 30 2e 30 3b 0a 58 | f-g; f |= 0.0;.X|
|00001950| 20 20 20 20 20 20 20 20 | 20 72 65 61 6c 2a 20 7a | | real* z|
|00001960| 6a 69 20 3d 20 7a 20 2b | 20 69 3b 20 72 65 61 6c |ji = z +| i; real|
|00001970| 2a 20 7a 69 6a 20 3d 20 | 7a 20 2b 20 69 2a 6e 3b |* zij = |z + i*n;|
|00001980| 20 72 65 61 6c 2a 20 65 | 6a 20 3d 20 45 2e 53 74 | real* e|j = E.St|
|00001990| 6f 72 65 28 29 3b 0a 58 | 09 20 66 6f 72 20 28 69 |ore();.X|. for (i|
|000019a0| 6e 74 20 6a 3d 30 3b 20 | 6a 3c 69 3b 20 6a 2b 2b |nt j=0; |j<i; j++|
|000019b0| 29 0a 58 09 20 7b 0a 58 | 09 20 20 20 20 2a 7a 6a |).X. {.X|. *zj|
|000019c0| 69 20 3d 20 28 2a 7a 69 | 6a 2b 2b 29 2f 68 3b 20 |i = (*zi|j++)/h; |
|000019d0| 67 20 3d 20 30 2e 30 3b | 0a 58 20 20 20 20 20 20 |g = 0.0;|.X |
|000019e0| 20 20 20 20 20 20 72 65 | 61 6c 2a 20 7a 6a 6b 20 | re|al* zjk |
|000019f0| 3d 20 7a 20 2b 20 6a 2a | 6e 3b 20 7a 69 6b 20 3d |= z + j*|n; zik =|
|00001a00| 20 7a 20 2b 20 69 2a 6e | 3b 0a 58 20 20 20 20 20 | z + i*n|;.X |
|00001a10| 20 20 20 20 20 20 20 6b | 20 3d 20 6a 3b 20 77 68 | k| = j; wh|
|00001a20| 69 6c 65 20 28 6b 2d 2d | 29 20 67 20 2b 3d 20 2a |ile (k--|) g += *|
|00001a30| 7a 6a 6b 2b 2b 20 2a 20 | 28 2a 7a 69 6b 2b 2b 29 |zjk++ * |(*zik++)|
|00001a40| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 20 20 20 6b |;.X | k|
|00001a50| 20 3d 20 69 2d 6a 3b 20 | 77 68 69 6c 65 20 28 6b | = i-j; |while (k|
|00001a60| 2d 2d 29 20 7b 20 67 20 | 2b 3d 20 2a 7a 6a 6b 20 |--) { g |+= *zjk |
|00001a70| 2a 20 28 2a 7a 69 6b 2b | 2b 29 3b 20 7a 6a 6b 20 |* (*zik+|+); zjk |
|00001a80| 2b 3d 20 6e 3b 20 7d 0a | 58 09 20 20 20 20 2a 65 |+= n; }.|X. *e|
|00001a90| 6a 2b 2b 20 3d 20 67 2f | 68 3b 20 66 20 2b 3d 20 |j++ = g/|h; f += |
|00001aa0| 67 20 2a 20 28 2a 7a 6a | 69 29 3b 20 7a 6a 69 20 |g * (*zj|i); zji |
|00001ab0| 2b 3d 20 6e 3b 0a 58 09 | 20 7d 0a 58 09 20 72 65 |+= n;.X.| }.X. re|
|00001ac0| 61 6c 20 68 68 20 3d 20 | 66 20 2f 20 28 68 20 2b |al hh = |f / (h +|
|00001ad0| 20 68 29 3b 20 7a 69 6a | 20 3d 20 7a 20 2b 20 69 | h); zij| = z + i|
|00001ae0| 2a 6e 3b 20 65 6a 20 3d | 20 45 2e 53 74 6f 72 65 |*n; ej =| E.Store|
|00001af0| 28 29 3b 0a 58 09 20 66 | 6f 72 20 28 6a 3d 30 3b |();.X. f|or (j=0;|
|00001b00| 20 6a 3c 69 3b 20 6a 2b | 2b 29 0a 58 09 20 7b 0a | j<i; j+|+).X. {.|
|00001b10| 58 09 20 20 20 20 66 20 | 3d 20 2a 7a 69 6a 2b 2b |X. f |= *zij++|
|00001b20| 3b 20 67 20 3d 20 2a 65 | 6a 20 2d 20 68 68 20 2a |; g = *e|j - hh *|
|00001b30| 20 66 3b 20 2a 65 6a 2b | 2b 20 3d 20 67 3b 0a 58 | f; *ej+|+ = g;.X|
|00001b40| 20 20 20 20 20 20 20 20 | 20 20 20 20 72 65 61 6c | | real|
|00001b50| 2a 20 7a 6a 6b 20 3d 20 | 7a 20 2b 20 6a 2a 6e 3b |* zjk = |z + j*n;|
|00001b60| 20 72 65 61 6c 2a 20 7a | 69 6b 20 3d 20 7a 20 2b | real* z|ik = z +|
|00001b70| 20 69 2a 6e 3b 0a 58 20 | 20 20 20 20 20 20 20 20 | i*n;.X | |
|00001b80| 20 20 20 72 65 61 6c 2a | 20 65 6b 20 3d 20 45 2e | real*| ek = E.|
|00001b90| 53 74 6f 72 65 28 29 3b | 20 6b 20 3d 20 6a 2b 31 |Store();| k = j+1|
|00001ba0| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 20 20 20 77 |;.X | w|
|00001bb0| 68 69 6c 65 20 28 6b 2d | 2d 29 20 20 2a 7a 6a 6b |hile (k-|-) *zjk|
|00001bc0| 2b 2b 20 2d 3d 20 28 20 | 66 2a 28 2a 65 6b 2b 2b |++ -= ( |f*(*ek++|
|00001bd0| 29 20 2b 20 67 2a 28 2a | 7a 69 6b 2b 2b 29 20 29 |) + g*(*|zik++) )|
|00001be0| 3b 20 0a 58 09 20 7d 0a | 58 20 20 20 20 20 20 7d |; .X. }.|X }|
|00001bf0| 0a 58 20 20 20 20 20 20 | 44 2e 65 6c 65 6d 65 6e |.X |D.elemen|
|00001c00| 74 28 69 29 20 3d 20 68 | 3b 0a 58 20 20 20 7d 0a |t(i) = h|;.X }.|
|00001c10| 58 0a 58 20 20 20 44 2e | 65 6c 65 6d 65 6e 74 28 |X.X D.|element(|
|00001c20| 30 29 20 3d 20 30 2e 30 | 3b 20 45 2e 65 6c 65 6d |0) = 0.0|; E.elem|
|00001c30| 65 6e 74 28 30 29 20 3d | 20 30 2e 30 3b 0a 58 20 |ent(0) =| 0.0;.X |
|00001c40| 20 20 66 6f 72 20 28 69 | 3d 30 3b 20 69 3c 6e 3b | for (i|=0; i<n;|
|00001c50| 20 69 2b 2b 29 0a 58 20 | 20 20 7b 0a 58 20 20 20 | i++).X | {.X |
|00001c60| 20 20 20 69 66 20 28 44 | 2e 65 6c 65 6d 65 6e 74 | if (D|.element|
|00001c70| 28 69 29 20 21 3d 20 30 | 2e 30 29 0a 58 20 20 20 |(i) != 0|.0).X |
|00001c80| 20 20 20 7b 0a 58 09 20 | 66 6f 72 20 28 69 6e 74 | {.X. |for (int|
|00001c90| 20 6a 3d 30 3b 20 6a 3c | 69 3b 20 6a 2b 2b 29 0a | j=0; j<|i; j++).|
|00001ca0| 58 09 20 7b 0a 58 09 20 | 20 20 20 72 65 61 6c 20 |X. {.X. | real |
|00001cb0| 67 20 3d 20 30 2e 30 3b | 0a 58 20 20 20 20 20 20 |g = 0.0;|.X |
|00001cc0| 20 20 20 20 20 20 72 65 | 61 6c 2a 20 7a 69 6b 20 | re|al* zik |
|00001cd0| 3d 20 7a 20 2b 20 69 2a | 6e 3b 20 72 65 61 6c 2a |= z + i*|n; real*|
|00001ce0| 20 7a 6b 6a 20 3d 20 7a | 20 2b 20 6a 3b 0a 58 20 | zkj = z| + j;.X |
|00001cf0| 20 20 20 20 20 20 20 20 | 20 20 20 69 6e 74 20 6b | | int k|
|00001d00| 20 3d 20 69 3b 20 77 68 | 69 6c 65 20 28 6b 2d 2d | = i; wh|ile (k--|
|00001d10| 29 20 7b 20 67 20 2b 3d | 20 2a 7a 69 6b 2b 2b 20 |) { g +=| *zik++ |
|00001d20| 2a 20 28 2a 7a 6b 6a 29 | 3b 20 7a 6b 6a 20 2b 3d |* (*zkj)|; zkj +=|
|00001d30| 20 6e 3b 20 7d 0a 58 20 | 20 20 20 20 20 20 20 20 | n; }.X | |
|00001d40| 20 20 20 72 65 61 6c 2a | 20 7a 6b 69 20 3d 20 7a | real*| zki = z|
|00001d50| 20 2b 20 69 3b 20 7a 6b | 6a 20 3d 20 7a 20 2b 20 | + i; zk|j = z + |
|00001d60| 6a 3b 0a 58 20 20 20 20 | 20 20 20 20 20 20 20 20 |j;.X | |
|00001d70| 6b 20 3d 20 69 3b 20 77 | 68 69 6c 65 20 28 6b 2d |k = i; w|hile (k-|
|00001d80| 2d 29 20 7b 20 2a 7a 6b | 6a 20 2d 3d 20 67 20 2a |-) { *zk|j -= g *|
|00001d90| 20 28 2a 7a 6b 69 29 3b | 20 7a 6b 6a 20 2b 3d 20 | (*zki);| zkj += |
|00001da0| 6e 3b 20 7a 6b 69 20 2b | 3d 20 6e 3b 20 7d 0a 58 |n; zki +|= n; }.X|
|00001db0| 09 20 7d 0a 58 20 20 20 | 20 20 20 7d 0a 58 20 20 |. }.X | }.X |
|00001dc0| 20 20 20 20 72 65 61 6c | 2a 20 7a 69 6a 20 3d 20 | real|* zij = |
|00001dd0| 7a 20 2b 20 69 2a 6e 3b | 20 72 65 61 6c 2a 20 7a |z + i*n;| real* z|
|00001de0| 6a 69 20 3d 20 7a 20 2b | 20 69 3b 0a 58 20 20 20 |ji = z +| i;.X |
|00001df0| 20 20 20 69 6e 74 20 6a | 20 3d 20 69 3b 20 77 68 | int j| = i; wh|
|00001e00| 69 6c 65 20 28 6a 2d 2d | 29 20 20 7b 20 2a 7a 69 |ile (j--|) { *zi|
|00001e10| 6a 2b 2b 20 3d 20 30 2e | 30 3b 20 2a 7a 6a 69 20 |j++ = 0.|0; *zji |
|00001e20| 3d 20 30 2e 30 3b 20 7a | 6a 69 20 2b 3d 20 6e 3b |= 0.0; z|ji += n;|
|00001e30| 20 7d 0a 58 20 20 20 20 | 20 20 44 2e 65 6c 65 6d | }.X | D.elem|
|00001e40| 65 6e 74 28 69 29 20 3d | 20 2a 7a 69 6a 3b 20 2a |ent(i) =| *zij; *|
|00001e50| 7a 69 6a 20 3d 20 31 2e | 30 3b 0a 58 20 20 20 7d |zij = 1.|0;.X }|
|00001e60| 0a 58 7d 0a 58 0a 58 73 | 74 61 74 69 63 20 76 6f |.X}.X.Xs|tatic vo|
|00001e70| 69 64 20 74 71 6c 32 28 | 44 69 61 67 6f 6e 61 6c |id tql2(|Diagonal|
|00001e80| 4d 61 74 72 69 78 26 20 | 44 2c 20 44 69 61 67 6f |Matrix& |D, Diago|
|00001e90| 6e 61 6c 4d 61 74 72 69 | 78 26 20 45 2c 20 4d 61 |nalMatri|x& E, Ma|
|00001ea0| 74 72 69 78 26 20 5a 29 | 0a 58 7b 0a 58 20 20 20 |trix& Z)|.X{.X |
|00001eb0| 72 65 61 6c 20 65 70 73 | 20 3d 20 46 6c 6f 61 74 |real eps| = Float|
|00001ec0| 69 6e 67 50 6f 69 6e 74 | 50 72 65 63 69 73 69 6f |ingPoint|Precisio|
|00001ed0| 6e 3a 3a 45 70 73 69 6c | 6f 6e 28 29 3b 0a 58 20 |n::Epsil|on();.X |
|00001ee0| 20 20 69 6e 74 20 6e 20 | 3d 20 44 2e 4e 72 6f 77 | int n |= D.Nrow|
|00001ef0| 73 28 29 3b 20 72 65 61 | 6c 2a 20 7a 20 3d 20 5a |s(); rea|l* z = Z|
|00001f00| 2e 53 74 6f 72 65 28 29 | 3b 0a 58 20 20 20 66 6f |.Store()|;.X fo|
|00001f10| 72 20 28 69 6e 74 20 6c | 3d 31 3b 20 6c 3c 6e 3b |r (int l|=1; l<n;|
|00001f20| 20 6c 2b 2b 29 20 45 2e | 65 6c 65 6d 65 6e 74 28 | l++) E.|element(|
|00001f30| 6c 2d 31 29 20 3d 20 45 | 2e 65 6c 65 6d 65 6e 74 |l-1) = E|.element|
|00001f40| 28 6c 29 3b 0a 58 20 20 | 20 72 65 61 6c 20 62 20 |(l);.X | real b |
|00001f50| 3d 20 30 2e 30 3b 20 72 | 65 61 6c 20 66 20 3d 20 |= 0.0; r|eal f = |
|00001f60| 30 2e 30 3b 20 45 2e 65 | 6c 65 6d 65 6e 74 28 6e |0.0; E.e|lement(n|
|00001f70| 2d 31 29 20 3d 20 30 2e | 30 3b 0a 58 20 20 20 66 |-1) = 0.|0;.X f|
|00001f80| 6f 72 20 28 6c 3d 30 3b | 20 6c 3c 6e 3b 20 6c 2b |or (l=0;| l<n; l+|
|00001f90| 2b 29 0a 58 20 20 20 7b | 0a 58 20 20 20 20 20 20 |+).X {|.X |
|00001fa0| 69 6e 74 20 69 2c 6a 3b | 0a 58 20 20 20 20 20 20 |int i,j;|.X |
|00001fb0| 72 65 61 6c 26 20 64 6c | 20 3d 20 44 2e 65 6c 65 |real& dl| = D.ele|
|00001fc0| 6d 65 6e 74 28 6c 29 3b | 20 72 65 61 6c 26 20 65 |ment(l);| real& e|
|00001fd0| 6c 20 3d 20 45 2e 65 6c | 65 6d 65 6e 74 28 6c 29 |l = E.el|ement(l)|
|00001fe0| 3b 0a 58 20 20 20 20 20 | 20 72 65 61 6c 20 68 20 |;.X | real h |
|00001ff0| 3d 20 65 70 73 20 2a 20 | 28 20 66 61 62 73 28 64 |= eps * |( fabs(d|
|00002000| 6c 29 20 2b 20 66 61 62 | 73 28 65 6c 29 20 29 3b |l) + fab|s(el) );|
|00002010| 0a 58 20 20 20 20 20 20 | 69 66 20 28 62 20 3c 20 |.X |if (b < |
|00002020| 68 29 20 62 20 3d 20 68 | 3b 0a 58 20 20 20 20 20 |h) b = h|;.X |
|00002030| 20 66 6f 72 20 28 69 6e | 74 20 6d 3d 6c 3b 20 6d | for (in|t m=l; m|
|00002040| 3c 6e 3b 20 6d 2b 2b 29 | 20 69 66 20 28 66 61 62 |<n; m++)| if (fab|
|00002050| 73 28 45 2e 65 6c 65 6d | 65 6e 74 28 6d 29 29 20 |s(E.elem|ent(m)) |
|00002060| 3c 3d 20 62 29 20 62 72 | 65 61 6b 3b 0a 58 20 20 |<= b) br|eak;.X |
|00002070| 20 20 20 20 66 6f 72 20 | 28 6a 3d 30 3b 20 6a 3c | for |(j=0; j<|
|00002080| 33 30 3b 20 6a 2b 2b 29 | 0a 58 20 20 20 20 20 20 |30; j++)|.X |
|00002090| 7b 0a 58 09 20 69 66 20 | 28 6d 3d 3d 6c 29 20 67 |{.X. if |(m==l) g|
|000020a0| 6f 74 6f 20 72 6f 6f 74 | 3b 0a 58 09 20 72 65 61 |oto root|;.X. rea|
|000020b0| 6c 26 20 64 6c 31 20 3d | 20 44 2e 65 6c 65 6d 65 |l& dl1 =| D.eleme|
|000020c0| 6e 74 28 6c 2b 31 29 3b | 0a 58 09 20 72 65 61 6c |nt(l+1);|.X. real|
|000020d0| 20 67 20 3d 20 64 6c 3b | 20 72 65 61 6c 20 70 20 | g = dl;| real p |
|000020e0| 3d 20 28 64 6c 31 2d 67 | 29 20 2f 20 28 32 2e 30 |= (dl1-g|) / (2.0|
|000020f0| 2a 65 6c 29 3b 20 72 65 | 61 6c 20 72 20 3d 20 73 |*el); re|al r = s|
|00002100| 71 72 74 28 70 2a 70 20 | 2b 20 31 2e 30 29 3b 0a |qrt(p*p |+ 1.0);.|
|00002110| 58 09 20 64 6c 20 3d 20 | 65 6c 20 2f 20 28 70 20 |X. dl = |el / (p |
|00002120| 3c 20 30 2e 30 20 3f 20 | 70 2d 72 20 3a 20 70 2b |< 0.0 ? |p-r : p+|
|00002130| 72 29 3b 20 72 65 61 6c | 20 68 20 3d 20 67 20 2d |r); real| h = g -|
|00002140| 20 64 6c 3b 20 66 20 2b | 3d 20 68 3b 0a 58 09 20 | dl; f +|= h;.X. |
|00002150| 72 65 61 6c 2a 20 64 6c | 78 20 3d 20 26 64 6c 31 |real* dl|x = &dl1|
|00002160| 3b 20 69 20 3d 20 6e 2d | 6c 2d 31 3b 20 77 68 69 |; i = n-|l-1; whi|
|00002170| 6c 65 20 28 69 2d 2d 29 | 20 2a 64 6c 78 2b 2b 20 |le (i--)| *dlx++ |
|00002180| 2d 3d 20 68 3b 0a 58 0a | 58 09 20 70 20 3d 20 44 |-= h;.X.|X. p = D|
|00002190| 2e 65 6c 65 6d 65 6e 74 | 28 6d 29 3b 20 72 65 61 |.element|(m); rea|
|000021a0| 6c 20 63 20 3d 20 31 2e | 30 3b 20 72 65 61 6c 20 |l c = 1.|0; real |
|000021b0| 73 20 3d 20 30 2e 30 3b | 0a 58 09 20 66 6f 72 20 |s = 0.0;|.X. for |
|000021c0| 28 69 3d 6d 2d 31 3b 20 | 69 3e 3d 6c 3b 20 69 2d |(i=m-1; |i>=l; i-|
|000021d0| 2d 29 0a 58 09 20 7b 0a | 58 09 20 20 20 20 72 65 |-).X. {.|X. re|
|000021e0| 61 6c 20 65 69 20 3d 20 | 45 2e 65 6c 65 6d 65 6e |al ei = |E.elemen|
|000021f0| 74 28 69 29 3b 20 72 65 | 61 6c 20 64 69 20 3d 20 |t(i); re|al di = |
|00002200| 44 2e 65 6c 65 6d 65 6e | 74 28 69 29 3b 0a 58 09 |D.elemen|t(i);.X.|
|00002210| 20 20 20 20 72 65 61 6c | 26 20 65 69 31 20 3d 20 | real|& ei1 = |
|00002220| 45 2e 65 6c 65 6d 65 6e | 74 28 69 2b 31 29 3b 0a |E.elemen|t(i+1);.|
|00002230| 58 09 20 20 20 20 67 20 | 3d 20 63 20 2a 20 65 69 |X. g |= c * ei|
|00002240| 3b 20 68 20 3d 20 63 20 | 2a 20 70 3b 0a 58 09 20 |; h = c |* p;.X. |
|00002250| 20 20 20 69 66 20 28 20 | 66 61 62 73 28 70 29 20 | if ( |fabs(p) |
|00002260| 3e 3d 20 66 61 62 73 28 | 65 69 29 29 0a 58 09 20 |>= fabs(|ei)).X. |
|00002270| 20 20 20 7b 0a 58 09 20 | 20 20 20 20 20 20 63 20 | {.X. | c |
|00002280| 3d 20 65 69 20 2f 20 70 | 3b 20 72 20 3d 20 73 71 |= ei / p|; r = sq|
|00002290| 72 74 28 63 2a 63 20 2b | 20 31 2e 30 29 3b 0a 58 |rt(c*c +| 1.0);.X|
|000022a0| 09 20 20 20 20 20 20 20 | 65 69 31 20 3d 20 73 2a |. |ei1 = s*|
|000022b0| 70 2a 72 3b 20 73 20 3d | 20 63 2f 72 3b 20 63 20 |p*r; s =| c/r; c |
|000022c0| 3d 20 31 2e 30 2f 72 3b | 0a 58 09 20 20 20 20 7d |= 1.0/r;|.X. }|
|000022d0| 0a 58 09 20 20 20 20 65 | 6c 73 65 0a 58 09 20 20 |.X. e|lse.X. |
|000022e0| 20 20 7b 0a 58 09 20 20 | 20 20 20 20 20 63 20 3d | {.X. | c =|
|000022f0| 20 70 20 2f 20 65 69 3b | 20 72 20 3d 20 73 71 72 | p / ei;| r = sqr|
|00002300| 74 28 63 2a 63 20 2b 20 | 31 2e 30 29 3b 0a 58 09 |t(c*c + |1.0);.X.|
|00002310| 20 20 20 20 20 20 20 65 | 69 31 20 3d 20 73 20 2a | e|i1 = s *|
|00002320| 20 65 69 20 2a 20 72 3b | 20 73 20 3d 20 31 2e 30 | ei * r;| s = 1.0|
|00002330| 2f 72 3b 20 63 20 2f 3d | 20 72 3b 0a 58 09 20 20 |/r; c /=| r;.X. |
|00002340| 20 20 7d 0a 58 09 20 20 | 20 20 70 20 3d 20 63 20 | }.X. | p = c |
|00002350| 2a 20 64 69 20 2d 20 73 | 2a 67 3b 20 44 2e 65 6c |* di - s|*g; D.el|
|00002360| 65 6d 65 6e 74 28 69 2b | 31 29 20 3d 20 68 20 2b |ement(i+|1) = h +|
|00002370| 20 73 20 2a 20 28 63 2a | 67 20 2b 20 73 2a 64 69 | s * (c*|g + s*di|
|00002380| 29 3b 0a 58 0a 58 09 20 | 20 20 20 72 65 61 6c 2a |);.X.X. | real*|
|00002390| 20 7a 6b 69 20 3d 20 7a | 20 2b 20 69 3b 20 72 65 | zki = z| + i; re|
|000023a0| 61 6c 2a 20 7a 6b 69 31 | 20 3d 20 7a 6b 69 20 2b |al* zki1| = zki +|
|000023b0| 20 31 3b 20 69 6e 74 20 | 6b 20 3d 20 6e 3b 0a 58 | 1; int |k = n;.X|
|000023c0| 09 20 20 20 20 77 68 69 | 6c 65 20 28 6b 2d 2d 29 |. whi|le (k--)|
|000023d0| 0a 58 09 20 20 20 20 7b | 0a 58 09 20 20 20 20 20 |.X. {|.X. |
|000023e0| 20 20 68 20 3d 20 2a 7a | 6b 69 31 3b 20 2a 7a 6b | h = *z|ki1; *zk|
|000023f0| 69 31 20 3d 20 73 2a 28 | 2a 7a 6b 69 29 20 2b 20 |i1 = s*(|*zki) + |
|00002400| 63 2a 68 3b 20 2a 7a 6b | 69 20 3d 20 63 2a 28 2a |c*h; *zk|i = c*(*|
|00002410| 7a 6b 69 29 20 2d 20 73 | 2a 68 3b 0a 58 09 20 20 |zki) - s|*h;.X. |
|00002420| 20 20 20 20 20 7a 6b 69 | 20 2b 3d 20 6e 3b 20 7a | zki| += n; z|
|00002430| 6b 69 31 20 2b 3d 20 6e | 3b 0a 58 09 20 20 20 20 |ki1 += n|;.X. |
|00002440| 7d 0a 58 09 20 7d 0a 58 | 09 20 65 6c 20 3d 20 73 |}.X. }.X|. el = s|
|00002450| 2a 70 3b 20 64 6c 20 3d | 20 63 2a 70 3b 0a 58 09 |*p; dl =| c*p;.X.|
|00002460| 20 69 66 20 28 66 61 62 | 73 28 65 6c 29 20 3c 3d | if (fab|s(el) <=|
|00002470| 20 62 29 20 67 6f 74 6f | 20 72 6f 6f 74 3b 0a 58 | b) goto| root;.X|
|00002480| 20 20 20 20 20 20 7d 0a | 58 20 20 20 20 20 20 4d | }.|X M|
|00002490| 61 74 72 69 78 45 72 72 | 6f 72 28 22 45 69 67 65 |atrixErr|or("Eige|
|000024a0| 6e 76 61 6c 75 65 20 72 | 6f 75 74 69 6e 65 20 66 |nvalue r|outine f|
|000024b0| 61 69 6c 73 22 29 3b 0a | 58 20 20 20 72 6f 6f 74 |ails");.|X root|
|000024c0| 3a 0a 58 20 20 20 20 20 | 20 64 6c 20 2b 3d 20 66 |:.X | dl += f|
|000024d0| 3b 0a 58 20 20 20 7d 0a | 58 0a 58 20 20 20 66 6f |;.X }.|X.X fo|
|000024e0| 72 20 28 69 6e 74 20 69 | 3d 30 3b 20 69 3c 6e 3b |r (int i|=0; i<n;|
|000024f0| 20 69 2b 2b 29 0a 58 20 | 20 20 7b 0a 58 20 20 20 | i++).X | {.X |
|00002500| 20 20 20 69 6e 74 20 6b | 20 3d 20 69 3b 20 72 65 | int k| = i; re|
|00002510| 61 6c 20 70 20 3d 20 44 | 2e 65 6c 65 6d 65 6e 74 |al p = D|.element|
|00002520| 28 69 29 3b 0a 58 20 20 | 20 20 20 20 66 6f 72 20 |(i);.X | for |
|00002530| 28 69 6e 74 20 6a 3d 69 | 2b 31 3b 20 6a 3c 6e 3b |(int j=i|+1; j<n;|
|00002540| 20 6a 2b 2b 29 0a 58 20 | 20 20 20 20 20 20 20 20 | j++).X | |
|00002550| 7b 20 69 66 20 28 44 2e | 65 6c 65 6d 65 6e 74 28 |{ if (D.|element(|
|00002560| 6a 29 20 3c 20 70 29 20 | 7b 20 6b 20 3d 20 6a 3b |j) < p) |{ k = j;|
|00002570| 20 70 20 3d 20 44 2e 65 | 6c 65 6d 65 6e 74 28 6a | p = D.e|lement(j|
|00002580| 29 3b 20 7d 20 7d 0a 58 | 20 20 20 20 20 20 69 66 |); } }.X| if|
|00002590| 20 28 6b 20 21 3d 20 69 | 29 0a 58 20 20 20 20 20 | (k != i|).X |
|000025a0| 20 7b 0a 58 20 20 20 20 | 20 20 20 20 20 44 2e 65 | {.X | D.e|
|000025b0| 6c 65 6d 65 6e 74 28 6b | 29 20 3d 20 44 2e 65 6c |lement(k|) = D.el|
|000025c0| 65 6d 65 6e 74 28 69 29 | 3b 20 44 2e 65 6c 65 6d |ement(i)|; D.elem|
|000025d0| 65 6e 74 28 69 29 20 3d | 20 70 3b 20 69 6e 74 20 |ent(i) =| p; int |
|000025e0| 6a 20 3d 20 6e 3b 0a 58 | 09 20 72 65 61 6c 2a 20 |j = n;.X|. real* |
|000025f0| 7a 6a 69 20 3d 20 7a 20 | 2b 20 69 3b 20 72 65 61 |zji = z |+ i; rea|
|00002600| 6c 2a 20 7a 6a 6b 20 3d | 20 7a 20 2b 20 6b 3b 0a |l* zjk =| z + k;.|
|00002610| 58 20 20 20 20 20 20 20 | 20 20 77 68 69 6c 65 20 |X | while |
|00002620| 28 6a 2d 2d 29 20 7b 20 | 70 20 3d 20 2a 7a 6a 69 |(j--) { |p = *zji|
|00002630| 3b 20 2a 7a 6a 69 20 3d | 20 2a 7a 6a 6b 3b 20 2a |; *zji =| *zjk; *|
|00002640| 7a 6a 6b 20 3d 20 70 3b | 20 7a 6a 69 20 2b 3d 20 |zjk = p;| zji += |
|00002650| 6e 3b 20 7a 6a 6b 20 2b | 3d 20 6e 3b 20 7d 0a 58 |n; zjk +|= n; }.X|
|00002660| 20 20 20 20 20 20 7d 0a | 58 20 20 20 7d 0a 58 0a | }.|X }.X.|
|00002670| 58 7d 0a 58 0a 58 73 74 | 61 74 69 63 20 76 6f 69 |X}.X.Xst|atic voi|
|00002680| 64 20 74 72 65 64 33 28 | 63 6f 6e 73 74 20 53 79 |d tred3(|const Sy|
|00002690| 6d 6d 65 74 72 69 63 4d | 61 74 72 69 78 26 20 58 |mmetricM|atrix& X|
|000026a0| 2c 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |, Diagon|alMatrix|
|000026b0| 26 20 44 2c 0a 58 20 20 | 20 44 69 61 67 6f 6e 61 |& D,.X | Diagona|
|000026c0| 6c 4d 61 74 72 69 78 26 | 20 45 2c 20 53 79 6d 6d |lMatrix&| E, Symm|
|000026d0| 65 74 72 69 63 4d 61 74 | 72 69 78 26 20 41 29 0a |etricMat|rix& A).|
|000026e0| 58 7b 0a 58 20 20 20 72 | 65 61 6c 20 74 6f 6c 20 |X{.X r|eal tol |
|000026f0| 3d 0a 58 20 20 20 20 20 | 20 46 6c 6f 61 74 69 6e |=.X | Floatin|
|00002700| 67 50 6f 69 6e 74 50 72 | 65 63 69 73 69 6f 6e 3a |gPointPr|ecision:|
|00002710| 3a 4d 69 6e 69 6d 75 6d | 28 29 2f 46 6c 6f 61 74 |:Minimum|()/Float|
|00002720| 69 6e 67 50 6f 69 6e 74 | 50 72 65 63 69 73 69 6f |ingPoint|Precisio|
|00002730| 6e 3a 3a 45 70 73 69 6c | 6f 6e 28 29 3b 0a 58 20 |n::Epsil|on();.X |
|00002740| 20 20 69 6e 74 20 6e 20 | 3d 20 58 2e 4e 72 6f 77 | int n |= X.Nrow|
|00002750| 73 28 29 3b 20 41 20 3d | 20 58 2e 63 28 29 3b 20 |s(); A =| X.c(); |
|00002760| 44 2e 52 65 44 69 6d 65 | 6e 73 69 6f 6e 28 6e 29 |D.ReDime|nsion(n)|
|00002770| 3b 20 45 2e 52 65 44 69 | 6d 65 6e 73 69 6f 6e 28 |; E.ReDi|mension(|
|00002780| 6e 29 3b 0a 58 20 20 20 | 72 65 61 6c 2a 20 65 69 |n);.X |real* ei|
|00002790| 20 3d 20 45 2e 53 74 6f | 72 65 28 29 20 2b 20 6e | = E.Sto|re() + n|
|000027a0| 3b 0a 58 20 20 20 66 6f | 72 20 28 69 6e 74 20 69 |;.X fo|r (int i|
|000027b0| 20 3d 20 6e 2d 31 3b 20 | 69 20 3e 3d 20 30 3b 20 | = n-1; |i >= 0; |
|000027c0| 69 2d 2d 29 0a 58 20 20 | 20 7b 0a 58 20 20 20 20 |i--).X | {.X |
|000027d0| 20 20 72 65 61 6c 20 68 | 20 3d 20 30 2e 30 3b 20 | real h| = 0.0; |
|000027e0| 72 65 61 6c 20 66 3b 0a | 58 20 20 20 20 20 20 72 |real f;.|X r|
|000027f0| 65 61 6c 2a 20 64 20 3d | 20 44 2e 53 74 6f 72 65 |eal* d =| D.Store|
|00002800| 28 29 3b 20 72 65 61 6c | 2a 20 61 20 3d 20 41 2e |(); real|* a = A.|
|00002810| 53 74 6f 72 65 28 29 20 | 2b 20 28 69 2a 28 69 2b |Store() |+ (i*(i+|
|00002820| 31 29 29 2f 32 3b 20 69 | 6e 74 20 6b 20 3d 20 69 |1))/2; i|nt k = i|
|00002830| 3b 0a 58 20 20 20 20 20 | 20 77 68 69 6c 65 20 28 |;.X | while (|
|00002840| 6b 2d 2d 29 20 7b 20 66 | 20 3d 20 2a 61 2b 2b 3b |k--) { f| = *a++;|
|00002850| 20 2a 64 2b 2b 20 3d 20 | 66 3b 20 68 20 2b 3d 20 | *d++ = |f; h += |
|00002860| 73 71 75 61 72 65 28 66 | 29 3b 20 7d 0a 58 20 20 |square(f|); }.X |
|00002870| 20 20 20 20 69 66 20 28 | 68 20 3c 3d 20 74 6f 6c | if (|h <= tol|
|00002880| 29 20 7b 20 2a 28 2d 2d | 65 69 29 20 3d 20 30 2e |) { *(--|ei) = 0.|
|00002890| 30 3b 20 68 20 3d 20 30 | 2e 30 3b 20 7d 0a 58 20 |0; h = 0|.0; }.X |
|000028a0| 20 20 20 20 20 65 6c 73 | 65 0a 58 20 20 20 20 20 | els|e.X |
|000028b0| 20 7b 0a 58 09 20 72 65 | 61 6c 20 67 20 3d 20 73 | {.X. re|al g = s|
|000028c0| 69 67 6e 28 2d 73 71 72 | 74 28 68 29 2c 20 66 29 |ign(-sqr|t(h), f)|
|000028d0| 3b 20 2a 28 2d 2d 65 69 | 29 20 3d 20 67 3b 20 68 |; *(--ei|) = g; h|
|000028e0| 20 2d 3d 20 66 2a 67 3b | 0a 58 20 20 20 20 20 20 | -= f*g;|.X |
|000028f0| 20 20 20 66 20 2d 3d 20 | 67 3b 20 2a 28 64 2d 31 | f -= |g; *(d-1|
|00002900| 29 20 3d 20 66 3b 20 2a | 28 61 2d 31 29 20 3d 20 |) = f; *|(a-1) = |
|00002910| 66 3b 20 66 20 3d 20 30 | 2e 30 3b 0a 58 20 20 20 |f; f = 0|.0;.X |
|00002920| 20 20 20 20 20 20 72 65 | 61 6c 2a 20 64 6a 20 3d | re|al* dj =|
|00002930| 20 44 2e 53 74 6f 72 65 | 28 29 3b 20 72 65 61 6c | D.Store|(); real|
|00002940| 2a 20 65 6a 20 3d 20 45 | 2e 53 74 6f 72 65 28 29 |* ej = E|.Store()|
|00002950| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 66 6f 72 20 |;.X | for |
|00002960| 28 69 6e 74 20 6a 20 3d | 20 30 3b 20 6a 20 3c 20 |(int j =| 0; j < |
|00002970| 69 3b 20 6a 2b 2b 29 0a | 58 20 20 20 20 20 20 20 |i; j++).|X |
|00002980| 20 20 7b 0a 58 20 20 20 | 20 20 20 20 20 20 20 20 | {.X | |
|00002990| 20 72 65 61 6c 2a 20 64 | 6b 20 3d 20 44 2e 53 74 | real* d|k = D.St|
|000029a0| 6f 72 65 28 29 3b 20 72 | 65 61 6c 2a 20 61 6b 20 |ore(); r|eal* ak |
|000029b0| 3d 20 41 2e 53 74 6f 72 | 65 28 29 2b 28 6a 2a 28 |= A.Stor|e()+(j*(|
|000029c0| 6a 2b 31 29 29 2f 32 3b | 0a 58 20 20 20 20 20 20 |j+1))/2;|.X |
|000029d0| 20 20 20 20 20 20 72 65 | 61 6c 20 67 20 3d 20 30 | re|al g = 0|
|000029e0| 2e 30 3b 20 6b 20 3d 20 | 6a 3b 0a 58 20 20 20 20 |.0; k = |j;.X |
|000029f0| 20 20 20 20 20 20 20 20 | 77 68 69 6c 65 20 28 6b | |while (k|
|00002a00| 2d 2d 29 20 20 67 20 2b | 3d 20 2a 61 6b 2b 2b 20 |--) g +|= *ak++ |
|00002a10| 2a 20 2a 64 6b 2b 2b 3b | 0a 58 20 20 20 20 20 20 |* *dk++;|.X |
|00002a20| 20 20 20 20 20 20 6b 20 | 3d 20 69 2d 6a 3b 20 69 | k |= i-j; i|
|00002a30| 6e 74 20 6c 20 3d 20 6a | 3b 20 0a 58 20 20 20 20 |nt l = j|; .X |
|00002a40| 20 20 20 20 20 20 20 20 | 77 68 69 6c 65 20 28 6b | |while (k|
|00002a50| 2d 2d 29 20 7b 20 67 20 | 2b 3d 20 2a 61 6b 20 2a |--) { g |+= *ak *|
|00002a60| 20 2a 64 6b 2b 2b 3b 20 | 61 6b 20 2b 3d 20 2b 2b | *dk++; |ak += ++|
|00002a70| 6c 3b 20 7d 0a 58 09 20 | 20 20 20 67 20 2f 3d 20 |l; }.X. | g /= |
|00002a80| 68 3b 20 2a 65 6a 2b 2b | 20 3d 20 67 3b 20 66 20 |h; *ej++| = g; f |
|00002a90| 2b 3d 20 67 20 2a 20 2a | 64 6a 2b 2b 3b 0a 58 20 |+= g * *|dj++;.X |
|00002aa0| 20 20 20 20 20 20 20 20 | 7d 20 20 0a 58 09 20 72 | |} .X. r|
|00002ab0| 65 61 6c 20 68 68 20 3d | 20 66 20 2f 20 28 32 20 |eal hh =| f / (2 |
|00002ac0| 2a 20 68 29 3b 20 72 65 | 61 6c 2a 20 61 6b 20 3d |* h); re|al* ak =|
|00002ad0| 20 41 2e 53 74 6f 72 65 | 28 29 3b 0a 58 20 20 20 | A.Store|();.X |
|00002ae0| 20 20 20 20 20 20 64 6a | 20 3d 20 44 2e 53 74 6f | dj| = D.Sto|
|00002af0| 72 65 28 29 3b 20 65 6a | 20 3d 20 45 2e 53 74 6f |re(); ej| = E.Sto|
|00002b00| 72 65 28 29 3b 0a 58 20 | 20 20 20 20 20 20 20 20 |re();.X | |
|00002b10| 66 6f 72 20 28 6a 20 3d | 20 30 3b 20 6a 20 3c 20 |for (j =| 0; j < |
|00002b20| 69 3b 20 6a 2b 2b 29 0a | 58 20 20 20 20 20 20 20 |i; j++).|X |
|00002b30| 20 20 7b 0a 58 09 20 20 | 20 20 66 20 3d 20 2a 64 | {.X. | f = *d|
|00002b40| 6a 2b 2b 3b 20 67 20 3d | 20 2a 65 6a 20 2d 20 68 |j++; g =| *ej - h|
|00002b50| 68 20 2a 20 66 3b 20 2a | 65 6a 2b 2b 20 3d 20 67 |h * f; *|ej++ = g|
|00002b60| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 20 20 20 72 |;.X | r|
|00002b70| 65 61 6c 2a 20 64 6b 20 | 3d 20 44 2e 53 74 6f 72 |eal* dk |= D.Stor|
|00002b80| 65 28 29 3b 20 72 65 61 | 6c 2a 20 65 6b 20 3d 20 |e(); rea|l* ek = |
|00002b90| 45 2e 53 74 6f 72 65 28 | 29 3b 20 6b 20 3d 20 6a |E.Store(|); k = j|
|00002ba0| 2b 31 3b 0a 58 09 20 20 | 20 20 77 68 69 6c 65 20 |+1;.X. | while |
|00002bb0| 28 6b 2d 2d 29 20 7b 20 | 2a 61 6b 2b 2b 20 2d 3d |(k--) { |*ak++ -=|
|00002bc0| 20 28 66 20 2a 20 2a 65 | 6b 2b 2b 20 2b 20 67 20 | (f * *e|k++ + g |
|00002bd0| 2a 20 2a 64 6b 2b 2b 29 | 3b 20 7d 0a 58 09 20 7d |* *dk++)|; }.X. }|
|00002be0| 0a 58 20 20 20 20 20 20 | 7d 0a 58 20 20 20 20 20 |.X |}.X |
|00002bf0| 20 2a 64 20 3d 20 2a 61 | 3b 20 2a 61 20 3d 20 68 | *d = *a|; *a = h|
|00002c00| 3b 0a 58 20 20 20 7d 0a | 58 7d 0a 58 0a 58 73 74 |;.X }.|X}.X.Xst|
|00002c10| 61 74 69 63 20 76 6f 69 | 64 20 74 71 6c 31 28 44 |atic voi|d tql1(D|
|00002c20| 69 61 67 6f 6e 61 6c 4d | 61 74 72 69 78 26 20 44 |iagonalM|atrix& D|
|00002c30| 2c 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |, Diagon|alMatrix|
|00002c40| 26 20 45 29 0a 58 7b 0a | 58 20 20 20 72 65 61 6c |& E).X{.|X real|
|00002c50| 20 65 70 73 20 3d 20 46 | 6c 6f 61 74 69 6e 67 50 | eps = F|loatingP|
|00002c60| 6f 69 6e 74 50 72 65 63 | 69 73 69 6f 6e 3a 3a 45 |ointPrec|ision::E|
|00002c70| 70 73 69 6c 6f 6e 28 29 | 3b 0a 58 20 20 20 69 6e |psilon()|;.X in|
|00002c80| 74 20 6e 20 3d 20 44 2e | 4e 72 6f 77 73 28 29 3b |t n = D.|Nrows();|
|00002c90| 0a 58 20 20 20 66 6f 72 | 20 28 69 6e 74 20 6c 3d |.X for| (int l=|
|00002ca0| 31 3b 20 6c 3c 6e 3b 20 | 6c 2b 2b 29 20 45 2e 65 |1; l<n; |l++) E.e|
|00002cb0| 6c 65 6d 65 6e 74 28 6c | 2d 31 29 20 3d 20 45 2e |lement(l|-1) = E.|
|00002cc0| 65 6c 65 6d 65 6e 74 28 | 6c 29 3b 0a 58 20 20 20 |element(|l);.X |
|00002cd0| 72 65 61 6c 20 62 20 3d | 20 30 2e 30 3b 20 72 65 |real b =| 0.0; re|
|00002ce0| 61 6c 20 66 20 3d 20 30 | 2e 30 3b 20 45 2e 65 6c |al f = 0|.0; E.el|
|00002cf0| 65 6d 65 6e 74 28 6e 2d | 31 29 20 3d 20 30 2e 30 |ement(n-|1) = 0.0|
|00002d00| 3b 0a 58 20 20 20 66 6f | 72 20 28 6c 3d 30 3b 20 |;.X fo|r (l=0; |
|00002d10| 6c 3c 6e 3b 20 6c 2b 2b | 29 0a 58 20 20 20 7b 0a |l<n; l++|).X {.|
|00002d20| 58 20 20 20 20 20 20 69 | 6e 74 20 69 2c 6a 3b 0a |X i|nt i,j;.|
|00002d30| 58 20 20 20 20 20 20 72 | 65 61 6c 26 20 64 6c 20 |X r|eal& dl |
|00002d40| 3d 20 44 2e 65 6c 65 6d | 65 6e 74 28 6c 29 3b 20 |= D.elem|ent(l); |
|00002d50| 72 65 61 6c 26 20 65 6c | 20 3d 20 45 2e 65 6c 65 |real& el| = E.ele|
|00002d60| 6d 65 6e 74 28 6c 29 3b | 0a 58 20 20 20 20 20 20 |ment(l);|.X |
|00002d70| 72 65 61 6c 20 68 20 3d | 20 65 70 73 20 2a 20 28 |real h =| eps * (|
|00002d80| 20 66 61 62 73 28 64 6c | 29 20 2b 20 66 61 62 73 | fabs(dl|) + fabs|
|00002d90| 28 65 6c 29 20 29 3b 0a | 58 20 20 20 20 20 20 69 |(el) );.|X i|
|00002da0| 66 20 28 62 20 3c 20 68 | 29 20 62 20 3d 20 68 3b |f (b < h|) b = h;|
|00002db0| 0a 58 20 20 20 20 20 20 | 66 6f 72 20 28 69 6e 74 |.X |for (int|
|00002dc0| 20 6d 3d 6c 3b 20 6d 3c | 6e 3b 20 6d 2b 2b 29 20 | m=l; m<|n; m++) |
|00002dd0| 69 66 20 28 66 61 62 73 | 28 45 2e 65 6c 65 6d 65 |if (fabs|(E.eleme|
|00002de0| 6e 74 28 6d 29 29 20 3c | 3d 20 62 29 20 62 72 65 |nt(m)) <|= b) bre|
|00002df0| 61 6b 3b 0a 58 20 20 20 | 20 20 20 66 6f 72 20 28 |ak;.X | for (|
|00002e00| 6a 3d 30 3b 20 6a 3c 33 | 30 3b 20 6a 2b 2b 29 0a |j=0; j<3|0; j++).|
|00002e10| 58 20 20 20 20 20 20 7b | 0a 58 20 20 20 20 20 20 |X {|.X |
|00002e20| 20 20 20 69 66 20 28 6d | 3d 3d 6c 29 20 67 6f 74 | if (m|==l) got|
|00002e30| 6f 20 72 6f 6f 74 3b 0a | 58 20 20 20 20 20 20 20 |o root;.|X |
|00002e40| 20 20 72 65 61 6c 26 20 | 64 6c 31 20 3d 20 44 2e | real& |dl1 = D.|
|00002e50| 65 6c 65 6d 65 6e 74 28 | 6c 2b 31 29 3b 0a 58 09 |element(|l+1);.X.|
|00002e60| 20 72 65 61 6c 20 67 20 | 3d 20 64 6c 3b 20 72 65 | real g |= dl; re|
|00002e70| 61 6c 20 70 20 3d 20 28 | 64 6c 31 2d 67 29 20 2f |al p = (|dl1-g) /|
|00002e80| 20 28 32 2e 30 2a 65 6c | 29 3b 20 72 65 61 6c 20 | (2.0*el|); real |
|00002e90| 72 20 3d 20 73 71 72 74 | 28 70 2a 70 20 2b 20 31 |r = sqrt|(p*p + 1|
|00002ea0| 2e 30 29 3b 0a 58 09 20 | 64 6c 20 3d 20 65 6c 20 |.0);.X. |dl = el |
|00002eb0| 2f 20 28 70 20 3c 20 30 | 2e 30 20 3f 20 70 2d 72 |/ (p < 0|.0 ? p-r|
|00002ec0| 20 3a 20 70 2b 72 29 3b | 20 72 65 61 6c 20 68 20 | : p+r);| real h |
|00002ed0| 3d 20 67 20 2d 20 64 6c | 3b 20 66 20 2b 3d 20 68 |= g - dl|; f += h|
|00002ee0| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 72 65 61 6c |;.X | real|
|00002ef0| 2a 20 64 6c 78 20 3d 20 | 26 64 6c 31 3b 20 69 20 |* dlx = |&dl1; i |
|00002f00| 3d 20 6e 2d 6c 2d 31 3b | 20 77 68 69 6c 65 20 28 |= n-l-1;| while (|
|00002f10| 69 2d 2d 29 20 2a 64 6c | 78 2b 2b 20 2d 3d 20 68 |i--) *dl|x++ -= h|
|00002f20| 3b 0a 58 0a 58 09 20 70 | 20 3d 20 44 2e 65 6c 65 |;.X.X. p| = D.ele|
|00002f30| 6d 65 6e 74 28 6d 29 3b | 20 72 65 61 6c 20 63 20 |ment(m);| real c |
|00002f40| 3d 20 31 2e 30 3b 20 72 | 65 61 6c 20 73 20 3d 20 |= 1.0; r|eal s = |
|00002f50| 30 2e 30 3b 0a 58 09 20 | 66 6f 72 20 28 69 3d 6d |0.0;.X. |for (i=m|
|00002f60| 2d 31 3b 20 69 3e 3d 6c | 3b 20 69 2d 2d 29 0a 58 |-1; i>=l|; i--).X|
|00002f70| 09 20 7b 0a 58 20 20 20 | 20 20 20 20 20 20 20 20 |. {.X | |
|00002f80| 20 72 65 61 6c 20 65 69 | 20 3d 20 45 2e 65 6c 65 | real ei| = E.ele|
|00002f90| 6d 65 6e 74 28 69 29 3b | 20 72 65 61 6c 20 64 69 |ment(i);| real di|
|00002fa0| 20 3d 20 44 2e 65 6c 65 | 6d 65 6e 74 28 69 29 3b | = D.ele|ment(i);|
|00002fb0| 0a 58 20 20 20 20 20 20 | 20 20 20 20 20 20 72 65 |.X | re|
|00002fc0| 61 6c 26 20 65 69 31 20 | 3d 20 45 2e 65 6c 65 6d |al& ei1 |= E.elem|
|00002fd0| 65 6e 74 28 69 2b 31 29 | 3b 0a 58 09 20 20 20 20 |ent(i+1)|;.X. |
|00002fe0| 67 20 3d 20 63 20 2a 20 | 65 69 3b 20 68 20 3d 20 |g = c * |ei; h = |
|00002ff0| 63 20 2a 20 70 3b 0a 58 | 09 20 20 20 20 69 66 20 |c * p;.X|. if |
|00003000| 28 20 66 61 62 73 28 70 | 29 20 3e 3d 20 66 61 62 |( fabs(p|) >= fab|
|00003010| 73 28 65 69 29 29 0a 58 | 09 20 20 20 20 7b 0a 58 |s(ei)).X|. {.X|
|00003020| 09 20 20 20 20 20 20 20 | 63 20 3d 20 65 69 20 2f |. |c = ei /|
|00003030| 20 70 3b 20 72 20 3d 20 | 73 71 72 74 28 63 2a 63 | p; r = |sqrt(c*c|
|00003040| 20 2b 20 31 2e 30 29 3b | 20 0a 58 20 20 20 20 20 | + 1.0);| .X |
|00003050| 20 20 20 20 20 20 20 20 | 20 20 65 69 31 20 3d 20 | | ei1 = |
|00003060| 73 2a 70 2a 72 3b 20 73 | 20 3d 20 63 2f 72 3b 20 |s*p*r; s| = c/r; |
|00003070| 63 20 3d 20 31 2e 30 2f | 72 3b 0a 58 09 20 20 20 |c = 1.0/|r;.X. |
|00003080| 20 7d 0a 58 09 20 20 20 | 20 65 6c 73 65 0a 58 09 | }.X. | else.X.|
|00003090| 20 20 20 20 7b 0a 58 09 | 20 20 20 20 20 20 20 63 | {.X.| c|
|000030a0| 20 3d 20 70 20 2f 20 65 | 69 3b 20 72 20 3d 20 73 | = p / e|i; r = s|
|000030b0| 71 72 74 28 63 2a 63 20 | 2b 20 31 2e 30 29 3b 0a |qrt(c*c |+ 1.0);.|
|000030c0| 58 09 20 20 20 20 20 20 | 20 65 69 31 20 3d 20 73 |X. | ei1 = s|
|000030d0| 20 2a 20 65 69 20 2a 20 | 72 3b 20 73 20 3d 20 31 | * ei * |r; s = 1|
|000030e0| 2e 30 2f 72 3b 20 63 20 | 2f 3d 20 72 3b 0a 58 09 |.0/r; c |/= r;.X.|
|000030f0| 20 20 20 20 7d 0a 58 09 | 20 20 20 20 70 20 3d 20 | }.X.| p = |
|00003100| 63 20 2a 20 64 69 20 2d | 20 73 2a 67 3b 20 44 2e |c * di -| s*g; D.|
|00003110| 65 6c 65 6d 65 6e 74 28 | 69 2b 31 29 20 3d 20 68 |element(|i+1) = h|
|00003120| 20 2b 20 73 20 2a 20 28 | 63 2a 67 20 2b 20 73 2a | + s * (|c*g + s*|
|00003130| 64 69 29 3b 0a 58 09 20 | 7d 0a 58 09 20 65 6c 20 |di);.X. |}.X. el |
|00003140| 3d 20 73 2a 70 3b 20 64 | 6c 20 3d 20 63 2a 70 3b |= s*p; d|l = c*p;|
|00003150| 0a 58 09 20 69 66 20 28 | 66 61 62 73 28 65 6c 29 |.X. if (|fabs(el)|
|00003160| 20 3c 3d 20 62 29 20 67 | 6f 74 6f 20 72 6f 6f 74 | <= b) g|oto root|
|00003170| 3b 0a 58 20 20 20 20 20 | 20 7d 0a 58 20 20 20 20 |;.X | }.X |
|00003180| 20 20 4d 61 74 72 69 78 | 45 72 72 6f 72 28 22 45 | Matrix|Error("E|
|00003190| 69 67 65 6e 76 61 6c 75 | 65 20 72 6f 75 74 69 6e |igenvalu|e routin|
|000031a0| 65 20 66 61 69 6c 73 22 | 29 3b 0a 58 20 20 20 72 |e fails"|);.X r|
|000031b0| 6f 6f 74 3a 0a 58 20 20 | 20 20 20 20 72 65 61 6c |oot:.X | real|
|000031c0| 20 70 20 3d 20 64 6c 20 | 2b 20 66 3b 0a 58 20 20 | p = dl |+ f;.X |
|000031d0| 20 20 20 20 66 6f 72 20 | 28 69 3d 6c 3b 20 69 3e | for |(i=l; i>|
|000031e0| 30 3b 20 69 2d 2d 29 0a | 58 20 20 20 20 20 20 7b |0; i--).|X {|
|000031f0| 0a 58 20 20 20 20 20 20 | 20 20 20 69 66 20 28 70 |.X | if (p|
|00003200| 20 3c 20 44 2e 65 6c 65 | 6d 65 6e 74 28 69 2d 31 | < D.ele|ment(i-1|
|00003210| 29 29 20 44 2e 65 6c 65 | 6d 65 6e 74 28 69 29 20 |)) D.ele|ment(i) |
|00003220| 3d 20 44 2e 65 6c 65 6d | 65 6e 74 28 69 2d 31 29 |= D.elem|ent(i-1)|
|00003230| 3b 0a 58 20 20 20 20 20 | 20 20 20 20 65 6c 73 65 |;.X | else|
|00003240| 20 67 6f 74 6f 20 63 6f | 6e 74 32 3b 0a 58 20 20 | goto co|nt2;.X |
|00003250| 20 20 20 20 7d 0a 58 20 | 20 20 20 20 20 69 3d 30 | }.X | i=0|
|00003260| 3b 0a 58 20 20 20 63 6f | 6e 74 32 3a 0a 58 20 20 |;.X co|nt2:.X |
|00003270| 20 20 20 20 44 2e 65 6c | 65 6d 65 6e 74 28 69 29 | D.el|ement(i)|
|00003280| 20 3d 20 70 3b 0a 58 20 | 20 20 7d 0a 58 7d 0a 58 | = p;.X | }.X}.X|
|00003290| 0a 58 76 6f 69 64 20 45 | 69 67 65 6e 56 61 6c 75 |.Xvoid E|igenValu|
|000032a0| 65 73 28 63 6f 6e 73 74 | 20 53 79 6d 6d 65 74 72 |es(const| Symmetr|
|000032b0| 69 63 4d 61 74 72 69 78 | 26 20 41 2c 20 44 69 61 |icMatrix|& A, Dia|
|000032c0| 67 6f 6e 61 6c 4d 61 74 | 72 69 78 26 20 44 2c 20 |gonalMat|rix& D, |
|000032d0| 4d 61 74 72 69 78 26 20 | 5a 29 0a 58 7b 20 44 69 |Matrix& |Z).X{ Di|
|000032e0| 61 67 6f 6e 61 6c 4d 61 | 74 72 69 78 20 45 3b 20 |agonalMa|trix E; |
|000032f0| 74 72 65 64 32 28 41 2c | 20 44 2c 20 45 2c 20 5a |tred2(A,| D, E, Z|
|00003300| 29 3b 20 74 71 6c 32 28 | 44 2c 20 45 2c 20 5a 29 |); tql2(|D, E, Z)|
|00003310| 3b 20 7d 0a 58 0a 58 76 | 6f 69 64 20 45 69 67 65 |; }.X.Xv|oid Eige|
|00003320| 6e 56 61 6c 75 65 73 28 | 63 6f 6e 73 74 20 53 79 |nValues(|const Sy|
|00003330| 6d 6d 65 74 72 69 63 4d | 61 74 72 69 78 26 20 58 |mmetricM|atrix& X|
|00003340| 2c 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 |, Diagon|alMatrix|
|00003350| 26 20 44 29 0a 58 7b 20 | 44 69 61 67 6f 6e 61 6c |& D).X{ |Diagonal|
|00003360| 4d 61 74 72 69 78 20 45 | 3b 20 53 79 6d 6d 65 74 |Matrix E|; Symmet|
|00003370| 72 69 63 4d 61 74 72 69 | 78 20 41 3b 20 74 72 65 |ricMatri|x A; tre|
|00003380| 64 33 28 58 2c 44 2c 45 | 2c 41 29 3b 20 74 71 6c |d3(X,D,E|,A); tql|
|00003390| 31 28 44 2c 45 29 3b 20 | 7d 0a 58 0a 58 76 6f 69 |1(D,E); |}.X.Xvoi|
|000033a0| 64 20 45 69 67 65 6e 56 | 61 6c 75 65 73 28 63 6f |d EigenV|alues(co|
|000033b0| 6e 73 74 20 53 79 6d 6d | 65 74 72 69 63 4d 61 74 |nst Symm|etricMat|
|000033c0| 72 69 78 26 20 58 2c 20 | 44 69 61 67 6f 6e 61 6c |rix& X, |Diagonal|
|000033d0| 4d 61 74 72 69 78 26 20 | 44 2c 0a 58 20 20 20 53 |Matrix& |D,.X S|
|000033e0| 79 6d 6d 65 74 72 69 63 | 4d 61 74 72 69 78 26 20 |ymmetric|Matrix& |
|000033f0| 41 29 0a 58 7b 20 44 69 | 61 67 6f 6e 61 6c 4d 61 |A).X{ Di|agonalMa|
|00003400| 74 72 69 78 20 45 3b 20 | 74 72 65 64 33 28 58 2c |trix E; |tred3(X,|
|00003410| 44 2c 45 2c 41 29 3b 20 | 74 71 6c 31 28 44 2c 45 |D,E,A); |tql1(D,E|
|00003420| 29 3b 20 7d 0a 58 0a 45 | 4e 44 5f 4f 46 5f 46 49 |); }.X.E|ND_OF_FI|
|00003430| 4c 45 0a 69 66 20 74 65 | 73 74 20 37 35 37 32 20 |LE.if te|st 7572 |
|00003440| 2d 6e 65 20 60 77 63 20 | 2d 63 20 3c 27 65 76 61 |-ne `wc |-c <'eva|
|00003450| 6c 75 65 2e 63 78 78 27 | 60 3b 20 74 68 65 6e 0a |lue.cxx'|`; then.|
|00003460| 20 20 20 20 65 63 68 6f | 20 73 68 61 72 3a 20 5c | echo| shar: \|
|00003470| 22 27 65 76 61 6c 75 65 | 2e 63 78 78 27 5c 22 20 |"'evalue|.cxx'\" |
|00003480| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|00003490| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|000034a0| 6e 64 20 6f 66 20 27 65 | 76 61 6c 75 65 2e 63 78 |nd of 'e|value.cx|
|000034b0| 78 27 0a 66 69 0a 69 66 | 20 74 65 73 74 20 2d 66 |x'.fi.if| test -f|
|000034c0| 20 27 65 78 61 6d 70 6c | 65 2e 63 78 78 27 20 2d | 'exampl|e.cxx' -|
|000034d0| 61 20 22 24 7b 31 7d 22 | 20 21 3d 20 22 2d 63 22 |a "${1}"| != "-c"|
|000034e0| 20 3b 20 74 68 65 6e 20 | 0a 20 20 65 63 68 6f 20 | ; then |. echo |
|000034f0| 73 68 61 72 3a 20 57 69 | 6c 6c 20 6e 6f 74 20 63 |shar: Wi|ll not c|
|00003500| 6c 6f 62 62 65 72 20 65 | 78 69 73 74 69 6e 67 20 |lobber e|xisting |
|00003510| 66 69 6c 65 20 5c 22 27 | 65 78 61 6d 70 6c 65 2e |file \"'|example.|
|00003520| 63 78 78 27 5c 22 0a 65 | 6c 73 65 0a 65 63 68 6f |cxx'\".e|lse.echo|
|00003530| 20 73 68 61 72 3a 20 45 | 78 74 72 61 63 74 69 6e | shar: E|xtractin|
|00003540| 67 20 5c 22 27 65 78 61 | 6d 70 6c 65 2e 63 78 78 |g \"'exa|mple.cxx|
|00003550| 27 5c 22 20 5c 28 39 32 | 39 33 20 63 68 61 72 61 |'\" \(92|93 chara|
|00003560| 63 74 65 72 73 5c 29 0a | 73 65 64 20 22 73 2f 5e |cters\).|sed "s/^|
|00003570| 58 2f 2f 22 20 3e 27 65 | 78 61 6d 70 6c 65 2e 63 |X//" >'e|xample.c|
|00003580| 78 78 27 20 3c 3c 27 45 | 4e 44 5f 4f 46 5f 46 49 |xx' <<'E|ND_OF_FI|
|00003590| 4c 45 27 0a 58 2f 2f 24 | 24 20 65 78 61 6d 70 6c |LE'.X//$|$ exampl|
|000035a0| 65 2e 63 78 78 20 20 20 | 20 20 20 20 20 20 20 20 |e.cxx | |
|000035b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000035c0| 20 20 45 78 61 6d 70 6c | 65 20 6f 66 20 75 73 65 | Exampl|e of use|
|000035d0| 20 6f 66 20 6d 61 74 72 | 69 78 20 70 61 63 6b 61 | of matr|ix packa|
|000035e0| 67 65 0a 58 0a 58 23 64 | 65 66 69 6e 65 20 57 41 |ge.X.X#d|efine WA|
|000035f0| 4e 54 5f 53 54 52 45 41 | 4d 20 20 20 20 20 20 20 |NT_STREA|M |
|00003600| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 69 6e | | // in|
|00003610| 63 6c 75 64 65 2e 68 78 | 78 20 77 69 6c 6c 20 67 |clude.hx|x will g|
|00003620| 65 74 20 73 74 72 65 61 | 6d 20 66 6e 73 0a 58 23 |et strea|m fns.X#|
|00003630| 64 65 66 69 6e 65 20 57 | 41 4e 54 5f 4d 41 54 48 |define W|ANT_MATH|
|00003640| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003650| 20 20 20 20 2f 2f 20 69 | 6e 63 6c 75 64 65 2e 68 | // i|nclude.h|
|00003660| 78 78 20 77 69 6c 6c 20 | 67 65 74 20 6d 61 74 68 |xx will |get math|
|00003670| 20 66 6e 73 0a 58 0a 58 | 23 69 6e 63 6c 75 64 65 | fns.X.X|#include|
|00003680| 20 22 69 6e 63 6c 75 64 | 65 2e 68 78 78 22 20 20 | "includ|e.hxx" |
|00003690| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2f 2f 20 | | // |
|000036a0| 69 6e 63 6c 75 64 65 20 | 73 74 61 6e 64 61 72 64 |include |standard|
|000036b0| 20 66 69 6c 65 73 0a 58 | 0a 58 23 69 6e 63 6c 75 | files.X|.X#inclu|
|000036c0| 64 65 20 22 6e 65 77 6d | 61 74 61 70 2e 68 78 78 |de "newm|atap.hxx|
|000036d0| 22 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2f |" | /|
|000036e0| 2f 20 6e 65 65 64 20 6d | 61 74 72 69 78 20 61 70 |/ need m|atrix ap|
|000036f0| 70 6c 69 63 61 74 69 6f | 6e 73 0a 58 0a 58 72 65 |plicatio|ns.X.Xre|
|00003700| 61 6c 20 74 33 28 72 65 | 61 6c 29 3b 20 20 20 20 |al t3(re|al); |
|00003710| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003720| 20 20 20 2f 2f 20 72 6f | 75 6e 64 20 74 6f 20 33 | // ro|und to 3|
|00003730| 20 64 65 63 69 6d 61 6c | 20 70 6c 61 63 65 73 0a | decimal| places.|
|00003740| 58 0a 58 0a 58 2f 2f 20 | 64 65 6d 6f 6e 73 74 72 |X.X.X// |demonstr|
|00003750| 61 74 69 6f 6e 20 6f 66 | 20 6d 61 74 72 69 78 20 |ation of| matrix |
|00003760| 70 61 63 6b 61 67 65 20 | 6f 6e 20 6c 69 6e 65 61 |package |on linea|
|00003770| 72 20 72 65 67 72 65 73 | 73 69 6f 6e 20 70 72 6f |r regres|sion pro|
|00003780| 62 6c 65 6d 0a 58 0a 58 | 0a 58 76 6f 69 64 20 74 |blem.X.X|.Xvoid t|
|00003790| 65 73 74 31 28 72 65 61 | 6c 2a 20 79 2c 20 72 65 |est1(rea|l* y, re|
|000037a0| 61 6c 2a 20 78 31 2c 20 | 72 65 61 6c 2a 20 78 32 |al* x1, |real* x2|
|000037b0| 2c 20 69 6e 74 20 6e 6f | 62 73 2c 20 69 6e 74 20 |, int no|bs, int |
|000037c0| 6e 70 72 65 64 29 0a 58 | 7b 0a 58 20 20 20 63 6f |npred).X|{.X co|
|000037d0| 75 74 20 3c 3c 20 22 5c | 6e 5c 6e 54 65 73 74 20 |ut << "\|n\nTest |
|000037e0| 31 20 2d 20 74 72 61 64 | 69 74 69 6f 6e 61 6c 5c |1 - trad|itional\|
|000037f0| 6e 22 3b 0a 58 0a 58 20 | 20 20 2f 2f 20 74 72 61 |n";.X.X | // tra|
|00003800| 64 69 74 69 6f 6e 61 6c | 20 73 75 6d 20 6f 66 20 |ditional| sum of |
|00003810| 73 71 75 61 72 65 73 20 | 61 6e 64 20 70 72 6f 64 |squares |and prod|
|00003820| 75 63 74 73 20 6d 65 74 | 68 6f 64 20 6f 66 20 63 |ucts met|hod of c|
|00003830| 61 6c 63 75 6c 61 74 69 | 6f 6e 0a 58 20 20 20 2f |alculati|on.X /|
|00003840| 2f 20 77 69 74 68 20 73 | 75 62 74 72 61 63 74 69 |/ with s|ubtracti|
|00003850| 6f 6e 20 6f 66 20 6d 65 | 61 6e 73 0a 58 0a 58 20 |on of me|ans.X.X |
|00003860| 20 20 2f 2f 20 6d 61 6b | 65 20 6d 61 74 72 69 78 | // mak|e matrix|
|00003870| 20 6f 66 20 70 72 65 64 | 69 63 74 6f 72 20 76 61 | of pred|ictor va|
|00003880| 6c 75 65 73 0a 58 20 20 | 20 4d 61 74 72 69 78 20 |lues.X | Matrix |
|00003890| 58 28 6e 6f 62 73 2c 6e | 70 72 65 64 29 3b 0a 58 |X(nobs,n|pred);.X|
|000038a0| 0a 58 20 20 20 2f 2f 20 | 6c 6f 61 64 20 78 31 20 |.X // |load x1 |
|000038b0| 61 6e 64 20 78 32 20 69 | 6e 74 6f 20 58 0a 58 20 |and x2 i|nto X.X |
|000038c0| 20 20 2f 2f 20 20 20 20 | 5b 75 73 65 20 3c 3c 20 | // |[use << |
|000038d0| 72 61 74 68 65 72 20 74 | 68 61 6e 20 3d 20 77 69 |rather t|han = wi|
|000038e0| 74 68 20 73 75 62 6d 61 | 74 72 69 63 65 73 20 61 |th subma|trices a|
|000038f0| 6e 64 2f 6f 72 20 6c 6f | 61 64 69 6e 67 20 61 72 |nd/or lo|ading ar|
|00003900| 72 61 79 73 5d 0a 58 20 | 20 20 58 2e 43 6f 6c 75 |rays].X | X.Colu|
|00003910| 6d 6e 28 31 29 20 3c 3c | 20 78 31 3b 20 20 58 2e |mn(1) <<| x1; X.|
|00003920| 43 6f 6c 75 6d 6e 28 32 | 29 20 3c 3c 20 78 32 3b |Column(2|) << x2;|
|00003930| 0a 58 0a 58 20 20 20 2f | 2f 20 76 65 63 74 6f 72 |.X.X /|/ vector|
|00003940| 20 6f 66 20 59 20 76 61 | 6c 75 65 73 0a 58 20 20 | of Y va|lues.X |
|00003950| 20 43 6f 6c 75 6d 6e 56 | 65 63 74 6f 72 20 59 28 | ColumnV|ector Y(|
|00003960| 6e 6f 62 73 29 3b 20 59 | 20 3c 3c 20 79 3b 0a 58 |nobs); Y| << y;.X|
|00003970| 0a 58 20 20 20 2f 2f 20 | 6d 61 6b 65 20 76 65 63 |.X // |make vec|
|00003980| 74 6f 72 20 6f 66 20 31 | 73 0a 58 20 20 20 43 6f |tor of 1|s.X Co|
|00003990| 6c 75 6d 6e 56 65 63 74 | 6f 72 20 4f 6e 65 73 28 |lumnVect|or Ones(|
|000039a0| 6e 6f 62 73 29 3b 20 4f | 6e 65 73 20 3d 20 31 2e |nobs); O|nes = 1.|
|000039b0| 30 3b 0a 58 0a 58 20 20 | 20 2f 2f 20 63 61 6c 63 |0;.X.X | // calc|
|000039c0| 75 6c 61 74 65 20 6d 65 | 61 6e 73 20 28 61 76 65 |ulate me|ans (ave|
|000039d0| 72 61 67 65 73 29 20 6f | 66 20 78 31 20 61 6e 64 |rages) o|f x1 and|
|000039e0| 20 78 32 20 5b 20 2e 74 | 28 29 20 74 61 6b 65 73 | x2 [ .t|() takes|
|000039f0| 20 74 72 61 6e 73 70 6f | 73 65 5d 0a 58 20 20 20 | transpo|se].X |
|00003a00| 52 6f 77 56 65 63 74 6f | 72 20 4d 20 3d 20 4f 6e |RowVecto|r M = On|
|00003a10| 65 73 2e 74 28 29 20 2a | 20 58 20 2f 20 6e 6f 62 |es.t() *| X / nob|
|00003a20| 73 3b 0a 58 0a 58 20 20 | 20 2f 2f 20 61 6e 64 20 |s;.X.X | // and |
|00003a30| 73 75 62 74 72 61 63 74 | 20 6d 65 61 6e 73 20 66 |subtract| means f|
|00003a40| 72 6f 6d 20 78 31 20 61 | 6e 64 20 78 31 0a 58 20 |rom x1 a|nd x1.X |
|00003a50| 20 20 4d 61 74 72 69 78 | 20 58 43 28 6e 6f 62 73 | Matrix| XC(nobs|
|00003a60| 2c 6e 70 72 65 64 29 3b | 0a 58 20 20 20 58 43 20 |,npred);|.X XC |
|00003a70| 3d 20 58 20 2d 20 4f 6e | 65 73 20 2a 20 4d 3b 0a |= X - On|es * M;.|
|00003a80| 58 0a 58 20 20 20 2f 2f | 20 64 6f 20 74 68 65 20 |X.X //| do the |
|00003a90| 73 61 6d 65 20 74 6f 20 | 59 20 5b 6e 65 65 64 20 |same to |Y [need |
|00003aa0| 22 72 65 61 6c 22 20 74 | 6f 20 63 6f 6e 76 65 72 |"real" t|o conver|
|00003ab0| 74 20 31 78 31 20 6d 61 | 74 72 69 78 20 74 6f 20 |t 1x1 ma|trix to |
|00003ac0| 73 63 61 6c 61 72 5d 0a | 58 20 20 20 43 6f 6c 75 |scalar].|X Colu|
|00003ad0| 6d 6e 56 65 63 74 6f 72 | 20 59 43 28 6e 6f 62 73 |mnVector| YC(nobs|
|00003ae0| 29 3b 0a 58 20 20 20 72 | 65 61 6c 20 6d 20 3d 20 |);.X r|eal m = |
|00003af0| 72 65 61 6c 28 4f 6e 65 | 73 2e 74 28 29 20 2a 20 |real(One|s.t() * |
|00003b00| 59 29 20 2f 20 6e 6f 62 | 73 3b 20 20 59 43 20 3d |Y) / nob|s; YC =|
|00003b10| 20 59 20 2d 20 4f 6e 65 | 73 20 2a 20 6d 3b 0a 58 | Y - One|s * m;.X|
|00003b20| 0a 58 20 20 20 2f 2f 20 | 66 6f 72 6d 20 73 75 6d |.X // |form sum|
|00003b30| 20 6f 66 20 73 71 75 61 | 72 65 73 20 61 6e 64 20 | of squa|res and |
|00003b40| 70 72 6f 64 75 63 74 20 | 6d 61 74 72 69 78 0a 58 |product |matrix.X|
|00003b50| 20 20 20 2f 2f 20 20 20 | 20 5b 75 73 65 20 3c 3c | // | [use <<|
|00003b60| 20 72 61 74 68 65 72 20 | 74 68 61 6e 20 3d 20 66 | rather |than = f|
|00003b70| 6f 72 20 63 6f 70 79 69 | 6e 67 20 4d 61 74 72 69 |or copyi|ng Matri|
|00003b80| 78 20 69 6e 74 6f 20 53 | 79 6d 6d 65 74 72 69 63 |x into S|ymmetric|
|00003b90| 4d 61 74 72 69 78 5d 0a | 58 20 20 20 53 79 6d 6d |Matrix].|X Symm|
|00003ba0| 65 74 72 69 63 4d 61 74 | 72 69 78 20 53 53 51 3b |etricMat|rix SSQ;|
|00003bb0| 20 53 53 51 20 3c 3c 20 | 58 43 2e 74 28 29 20 2a | SSQ << |XC.t() *|
|00003bc0| 20 58 43 3b 0a 58 0a 58 | 20 20 20 2f 2f 20 63 61 | XC;.X.X| // ca|
|00003bd0| 6c 63 75 6c 61 74 65 20 | 65 73 74 69 6d 61 74 65 |lculate |estimate|
|00003be0| 0a 58 20 20 20 2f 2f 20 | 20 20 20 5b 62 72 61 63 |.X // | [brac|
|00003bf0| 6b 65 74 20 6c 61 73 74 | 20 74 77 6f 20 74 65 72 |ket last| two ter|
|00003c00| 6d 73 20 74 6f 20 66 6f | 72 63 65 20 74 68 69 73 |ms to fo|rce this|
|00003c10| 20 6d 75 6c 74 69 70 6c | 69 63 61 74 69 6f 6e 20 | multipl|ication |
|00003c20| 66 69 72 73 74 5d 0a 58 | 20 20 20 2f 2f 20 20 20 |first].X| // |
|00003c30| 20 5b 20 2e 69 28 29 20 | 6d 65 61 6e 73 20 69 6e | [ .i() |means in|
|00003c40| 76 65 72 73 65 2c 20 62 | 75 74 20 69 6e 76 65 72 |verse, b|ut inver|
|00003c50| 73 65 20 69 73 20 6e 6f | 74 20 65 78 70 6c 69 63 |se is no|t explic|
|00003c60| 69 74 79 20 63 61 6c 63 | 75 6c 61 74 65 64 5d 0a |ity calc|ulated].|
|00003c70| 58 20 20 20 43 6f 6c 75 | 6d 6e 56 65 63 74 6f 72 |X Colu|mnVector|
|00003c80| 20 41 20 3d 20 53 53 51 | 2e 69 28 29 20 2a 20 28 | A = SSQ|.i() * (|
|00003c90| 58 43 2e 74 28 29 20 2a | 20 59 43 29 3b 0a 58 0a |XC.t() *| YC);.X.|
|00003ca0| 58 20 20 20 2f 2f 20 63 | 61 6c 63 75 6c 61 74 65 |X // c|alculate|
|00003cb0| 20 65 73 74 69 6d 61 74 | 65 20 6f 66 20 63 6f 6e | estimat|e of con|
|00003cc0| 73 74 61 6e 74 20 74 65 | 72 6d 0a 58 20 20 20 72 |stant te|rm.X r|
|00003cd0| 65 61 6c 20 61 20 3d 20 | 6d 20 2d 20 72 65 61 6c |eal a = |m - real|
|00003ce0| 28 4d 20 2a 20 41 29 3b | 0a 58 0a 58 20 20 20 2f |(M * A);|.X.X /|
|00003cf0| 2f 20 47 65 74 20 76 61 | 72 69 61 6e 63 65 73 20 |/ Get va|riances |
|00003d00| 6f 66 20 65 73 74 69 6d | 61 74 65 73 20 66 72 6f |of estim|ates fro|
|00003d10| 6d 20 64 69 61 67 6f 6e | 61 6c 20 65 6c 65 6d 65 |m diagon|al eleme|
|00003d20| 6e 74 73 20 6f 66 20 69 | 6e 76 6f 69 63 65 20 6f |nts of i|nvoice o|
|00003d30| 66 20 53 53 51 0a 58 20 | 20 20 2f 2f 20 20 20 20 |f SSQ.X | // |
|00003d40| 5b 20 77 65 20 61 72 65 | 20 74 61 6b 69 6e 67 20 |[ we are| taking |
|00003d50| 69 6e 76 65 72 73 65 20 | 6f 66 20 53 53 51 3b 20 |inverse |of SSQ; |
|00003d60| 77 6f 75 6c 64 20 68 61 | 76 65 20 62 65 65 6e 20 |would ha|ve been |
|00003d70| 62 65 74 74 65 72 20 74 | 6f 20 75 73 65 0a 58 20 |better t|o use.X |
|00003d80| 20 20 2f 2f 20 20 20 20 | 20 20 20 20 43 72 6f 75 | // | Crou|
|00003d90| 74 4d 61 74 72 69 78 20 | 6d 65 74 68 6f 64 20 2d |tMatrix |method -|
|00003da0| 20 73 65 65 20 64 6f 63 | 75 6d 65 6e 74 61 74 69 | see doc|umentati|
|00003db0| 6f 6e 20 5d 0a 58 20 20 | 20 4d 61 74 72 69 78 20 |on ].X | Matrix |
|00003dc0| 49 53 53 51 20 3d 20 53 | 53 51 2e 69 28 29 3b 20 |ISSQ = S|SQ.i(); |
|00003dd0| 44 69 61 67 6f 6e 61 6c | 4d 61 74 72 69 78 20 44 |Diagonal|Matrix D|
|00003de0| 3b 20 44 20 3c 3c 20 49 | 53 53 51 3b 0a 58 20 20 |; D << I|SSQ;.X |
|00003df0| 20 43 6f 6c 75 6d 6e 56 | 65 63 74 6f 72 20 56 20 | ColumnV|ector V |
|00003e00| 3d 20 44 2e 43 6f 70 79 | 54 6f 43 6f 6c 75 6d 6e |= D.Copy|ToColumn|
|00003e10| 28 29 3b 0a 58 20 20 20 | 72 65 61 6c 20 76 20 3d |();.X |real v =|
|00003e20| 20 31 2e 30 2f 6e 6f 62 | 73 20 2b 20 72 65 61 6c | 1.0/nob|s + real|
|00003e30| 28 4d 20 2a 20 49 53 53 | 51 20 2a 20 4d 2e 74 28 |(M * ISS|Q * M.t(|
|00003e40| 29 29 3b 0a 58 09 09 09 | 09 09 20 20 20 20 2f 2f |));.X...|.. //|
|00003e50| 20 66 6f 72 20 63 61 6c | 63 20 76 61 72 69 61 6e | for cal|c varian|
|00003e60| 63 65 20 63 6f 6e 73 74 | 0a 58 0a 58 20 20 20 2f |ce const|.X.X /|
|00003e70| 2f 20 43 61 6c 63 75 6c | 61 74 65 20 66 69 74 74 |/ Calcul|ate fitt|
|00003e80| 65 64 20 76 61 6c 75 65 | 73 20 61 6e 64 20 72 65 |ed value|s and re|
|00003e90| 73 69 64 75 61 6c 73 0a | 58 20 20 20 69 6e 74 20 |siduals.|X int |
|00003ea0| 6e 70 72 65 64 31 20 3d | 20 6e 70 72 65 64 2b 31 |npred1 =| npred+1|
|00003eb0| 3b 0a 58 20 20 20 43 6f | 6c 75 6d 6e 56 65 63 74 |;.X Co|lumnVect|
|00003ec0| 6f 72 20 46 69 74 74 65 | 64 20 3d 20 58 20 2a 20 |or Fitte|d = X * |
|00003ed0| 41 20 2b 20 61 3b 0a 58 | 20 20 20 43 6f 6c 75 6d |A + a;.X| Colum|
|00003ee0| 6e 56 65 63 74 6f 72 20 | 52 65 73 69 64 75 61 6c |nVector |Residual|
|00003ef0| 20 3d 20 59 20 2d 20 46 | 69 74 74 65 64 3b 0a 58 | = Y - F|itted;.X|
|00003f00| 20 20 20 72 65 61 6c 20 | 52 65 73 56 61 72 20 3d | real |ResVar =|
|00003f10| 20 52 65 73 69 64 75 61 | 6c 2e 53 75 6d 53 71 75 | Residua|l.SumSqu|
|00003f20| 61 72 65 28 29 20 2f 20 | 28 6e 6f 62 73 2d 6e 70 |are() / |(nobs-np|
|00003f30| 72 65 64 31 29 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 |red1);.X|.X // |
|00003f40| 47 65 74 20 64 69 61 67 | 6f 6e 61 6c 73 20 6f 66 |Get diag|onals of|
|00003f50| 20 48 61 74 20 6d 61 74 | 72 69 78 20 28 61 6e 20 | Hat mat|rix (an |
|00003f60| 65 78 70 65 6e 73 69 76 | 65 20 77 61 79 20 6f 66 |expensiv|e way of|
|00003f70| 20 64 6f 69 6e 67 20 74 | 68 69 73 29 0a 58 20 20 | doing t|his).X |
|00003f80| 20 4d 61 74 72 69 78 20 | 58 31 28 6e 6f 62 73 2c | Matrix |X1(nobs,|
|00003f90| 6e 70 72 65 64 31 29 3b | 20 58 31 2e 43 6f 6c 75 |npred1);| X1.Colu|
|00003fa0| 6d 6e 28 31 29 3c 3c 4f | 6e 65 73 3b 20 58 31 2e |mn(1)<<O|nes; X1.|
|00003fb0| 43 6f 6c 75 6d 6e 73 28 | 32 2c 6e 70 72 65 64 31 |Columns(|2,npred1|
|00003fc0| 29 3c 3c 58 3b 0a 58 20 | 20 20 44 69 61 67 6f 6e |)<<X;.X | Diagon|
|00003fd0| 61 6c 4d 61 74 72 69 78 | 20 48 61 74 3b 20 20 48 |alMatrix| Hat; H|
|00003fe0| 61 74 20 3c 3c 20 58 31 | 20 2a 20 28 58 31 2e 74 |at << X1| * (X1.t|
|00003ff0| 28 29 20 2a 20 58 31 29 | 2e 69 28 29 20 2a 20 58 |() * X1)|.i() * X|
|00004000| 31 2e 74 28 29 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 |1.t();.X|.X // |
|00004010| 70 72 69 6e 74 20 6f 75 | 74 20 61 6e 73 77 65 72 |print ou|t answer|
|00004020| 73 0a 58 20 20 20 63 6f | 75 74 20 3c 3c 20 22 5c |s.X co|ut << "\|
|00004030| 6e 45 73 74 69 6d 61 74 | 65 73 20 61 6e 64 20 74 |nEstimat|es and t|
|00004040| 68 65 69 72 20 73 74 61 | 6e 64 61 72 64 20 65 72 |heir sta|ndard er|
|00004050| 72 6f 72 73 5c 6e 5c 6e | 22 3b 0a 58 20 20 20 63 |rors\n\n|";.X c|
|00004060| 6f 75 74 20 3c 3c 20 61 | 20 3c 3c 22 5c 74 22 3c |out << a| <<"\t"<|
|00004070| 3c 20 73 71 72 74 28 76 | 2a 52 65 73 56 61 72 29 |< sqrt(v|*ResVar)|
|00004080| 20 3c 3c 20 22 5c 6e 22 | 3b 0a 58 20 20 20 66 6f | << "\n"|;.X fo|
|00004090| 72 20 28 69 6e 74 20 69 | 3d 31 3b 20 69 3c 3d 6e |r (int i|=1; i<=n|
|000040a0| 70 72 65 64 3b 20 69 2b | 2b 29 0a 58 20 20 20 63 |pred; i+|+).X c|
|000040b0| 6f 75 74 20 3c 3c 20 41 | 28 69 29 20 3c 3c 22 5c |out << A|(i) <<"\|
|000040c0| 74 22 3c 3c 20 73 71 72 | 74 28 56 28 69 29 2a 52 |t"<< sqr|t(V(i)*R|
|000040d0| 65 73 56 61 72 29 20 3c | 3c 20 22 5c 6e 22 3b 0a |esVar) <|< "\n";.|
|000040e0| 58 20 20 20 63 6f 75 74 | 20 3c 3c 20 22 5c 6e 4f |X cout| << "\nO|
|000040f0| 62 73 65 72 76 61 74 69 | 6f 6e 73 2c 20 66 69 74 |bservati|ons, fit|
|00004100| 74 65 64 20 76 61 6c 75 | 65 2c 20 72 65 73 69 64 |ted valu|e, resid|
|00004110| 75 61 6c 20 76 61 6c 75 | 65 2c 20 68 61 74 20 76 |ual valu|e, hat v|
|00004120| 61 6c 75 65 5c 6e 22 3b | 0a 58 20 20 20 66 6f 72 |alue\n";|.X for|
|00004130| 20 28 69 3d 31 3b 20 69 | 3c 3d 6e 6f 62 73 3b 20 | (i=1; i|<=nobs; |
|00004140| 69 2b 2b 29 0a 58 20 20 | 20 20 20 20 63 6f 75 74 |i++).X | cout|
|00004150| 20 3c 3c 20 58 28 69 2c | 31 29 20 3c 3c 22 5c 74 | << X(i,|1) <<"\t|
|00004160| 22 3c 3c 20 58 28 69 2c | 32 29 20 3c 3c 22 5c 74 |"<< X(i,|2) <<"\t|
|00004170| 22 3c 3c 20 59 28 69 29 | 20 3c 3c 22 5c 74 22 3c |"<< Y(i)| <<"\t"<|
|00004180| 3c 0a 58 20 20 20 20 20 | 20 74 33 28 46 69 74 74 |<.X | t3(Fitt|
|00004190| 65 64 28 69 29 29 20 3c | 3c 22 5c 74 22 3c 3c 20 |ed(i)) <|<"\t"<< |
|000041a0| 74 33 28 52 65 73 69 64 | 75 61 6c 28 69 29 29 20 |t3(Resid|ual(i)) |
|000041b0| 3c 3c 22 5c 74 22 3c 3c | 20 74 33 28 48 61 74 28 |<<"\t"<<| t3(Hat(|
|000041c0| 69 29 29 20 3c 3c 22 5c | 6e 22 3b 0a 58 20 20 20 |i)) <<"\|n";.X |
|000041d0| 63 6f 75 74 20 3c 3c 20 | 22 5c 6e 5c 6e 22 3b 0a |cout << |"\n\n";.|
|000041e0| 58 7d 0a 58 0a 58 76 6f | 69 64 20 74 65 73 74 32 |X}.X.Xvo|id test2|
|000041f0| 28 72 65 61 6c 2a 20 79 | 2c 20 72 65 61 6c 2a 20 |(real* y|, real* |
|00004200| 78 31 2c 20 72 65 61 6c | 2a 20 78 32 2c 20 69 6e |x1, real|* x2, in|
|00004210| 74 20 6e 6f 62 73 2c 20 | 69 6e 74 20 6e 70 72 65 |t nobs, |int npre|
|00004220| 64 29 0a 58 7b 0a 58 20 | 20 20 63 6f 75 74 20 3c |d).X{.X | cout <|
|00004230| 3c 20 22 5c 6e 5c 6e 54 | 65 73 74 20 32 20 2d 20 |< "\n\nT|est 2 - |
|00004240| 43 68 6f 6c 65 73 6b 79 | 5c 6e 22 3b 0a 58 0a 58 |Cholesky|\n";.X.X|
|00004250| 20 20 20 2f 2f 20 74 72 | 61 64 69 74 69 6f 6e 61 | // tr|aditiona|
|00004260| 6c 20 73 75 6d 20 6f 66 | 20 73 71 75 61 72 65 73 |l sum of| squares|
|00004270| 20 61 6e 64 20 70 72 6f | 64 75 63 74 73 20 6d 65 | and pro|ducts me|
|00004280| 74 68 6f 64 20 6f 66 20 | 63 61 6c 63 75 6c 61 74 |thod of |calculat|
|00004290| 69 6f 6e 0a 58 20 20 20 | 2f 2f 20 77 69 74 68 20 |ion.X |// with |
|000042a0| 73 75 62 74 72 61 63 74 | 69 6f 6e 20 6f 66 20 6d |subtract|ion of m|
|000042b0| 65 61 6e 73 20 2d 20 75 | 73 69 6e 67 20 43 68 6f |eans - u|sing Cho|
|000042c0| 6c 65 73 6b 79 20 64 65 | 63 6f 6d 70 6f 73 69 74 |lesky de|composit|
|000042d0| 69 6f 6e 0a 58 0a 58 20 | 20 20 4d 61 74 72 69 78 |ion.X.X | Matrix|
|000042e0| 20 58 28 6e 6f 62 73 2c | 6e 70 72 65 64 29 3b 0a | X(nobs,|npred);.|
|000042f0| 58 20 20 20 58 2e 43 6f | 6c 75 6d 6e 28 31 29 20 |X X.Co|lumn(1) |
|00004300| 3c 3c 20 78 31 3b 20 20 | 58 2e 43 6f 6c 75 6d 6e |<< x1; |X.Column|
|00004310| 28 32 29 20 3c 3c 20 78 | 32 3b 0a 58 20 20 20 43 |(2) << x|2;.X C|
|00004320| 6f 6c 75 6d 6e 56 65 63 | 74 6f 72 20 59 28 6e 6f |olumnVec|tor Y(no|
|00004330| 62 73 29 3b 20 59 20 3c | 3c 20 79 3b 0a 58 20 20 |bs); Y <|< y;.X |
|00004340| 20 43 6f 6c 75 6d 6e 56 | 65 63 74 6f 72 20 4f 6e | ColumnV|ector On|
|00004350| 65 73 28 6e 6f 62 73 29 | 3b 20 4f 6e 65 73 20 3d |es(nobs)|; Ones =|
|00004360| 20 31 2e 30 3b 0a 58 20 | 20 20 52 6f 77 56 65 63 | 1.0;.X | RowVec|
|00004370| 74 6f 72 20 4d 20 3d 20 | 4f 6e 65 73 2e 74 28 29 |tor M = |Ones.t()|
|00004380| 20 2a 20 58 20 2f 20 6e | 6f 62 73 3b 0a 58 20 20 | * X / n|obs;.X |
|00004390| 20 4d 61 74 72 69 78 20 | 58 43 28 6e 6f 62 73 2c | Matrix |XC(nobs,|
|000043a0| 6e 70 72 65 64 29 3b 0a | 58 20 20 20 58 43 20 3d |npred);.|X XC =|
|000043b0| 20 58 20 2d 20 4f 6e 65 | 73 20 2a 20 4d 3b 0a 58 | X - One|s * M;.X|
|000043c0| 20 20 20 43 6f 6c 75 6d | 6e 56 65 63 74 6f 72 20 | Colum|nVector |
|000043d0| 59 43 28 6e 6f 62 73 29 | 3b 0a 58 20 20 20 72 65 |YC(nobs)|;.X re|
|000043e0| 61 6c 20 6d 20 3d 20 72 | 65 61 6c 28 4f 6e 65 73 |al m = r|eal(Ones|
|000043f0| 2e 74 28 29 20 2a 20 59 | 29 20 2f 20 6e 6f 62 73 |.t() * Y|) / nobs|
|00004400| 3b 20 20 59 43 20 3d 20 | 59 20 2d 20 4f 6e 65 73 |; YC = |Y - Ones|
|00004410| 20 2a 20 6d 3b 0a 58 20 | 20 20 53 79 6d 6d 65 74 | * m;.X | Symmet|
|00004420| 72 69 63 4d 61 74 72 69 | 78 20 53 53 51 3b 20 53 |ricMatri|x SSQ; S|
|00004430| 53 51 20 3c 3c 20 58 43 | 2e 74 28 29 20 2a 20 58 |SQ << XC|.t() * X|
|00004440| 43 3b 0a 58 0a 58 20 20 | 20 2f 2f 20 43 68 6f 6c |C;.X.X | // Chol|
|00004450| 65 73 6b 79 20 64 65 63 | 6f 6d 70 6f 73 69 74 69 |esky dec|ompositi|
|00004460| 6f 6e 20 6f 66 20 53 53 | 51 0a 58 20 20 20 4c 6f |on of SS|Q.X Lo|
|00004470| 77 65 72 54 72 69 61 6e | 67 75 6c 61 72 4d 61 74 |werTrian|gularMat|
|00004480| 72 69 78 20 4c 20 3d 20 | 43 68 6f 6c 65 73 6b 79 |rix L = |Cholesky|
|00004490| 28 53 53 51 29 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 |(SSQ);.X|.X // |
|000044a0| 63 61 6c 63 75 6c 61 74 | 65 20 65 73 74 69 6d 61 |calculat|e estima|
|000044b0| 74 65 0a 58 20 20 20 43 | 6f 6c 75 6d 6e 56 65 63 |te.X C|olumnVec|
|000044c0| 74 6f 72 20 41 20 3d 20 | 4c 2e 74 28 29 2e 69 28 |tor A = |L.t().i(|
|000044d0| 29 20 2a 20 28 4c 2e 69 | 28 29 20 2a 20 28 58 43 |) * (L.i|() * (XC|
|000044e0| 2e 74 28 29 20 2a 20 59 | 43 29 29 3b 0a 58 0a 58 |.t() * Y|C));.X.X|
|000044f0| 20 20 20 2f 2f 20 63 61 | 6c 63 75 6c 61 74 65 20 | // ca|lculate |
|00004500| 65 73 74 69 6d 61 74 65 | 20 6f 66 20 63 6f 6e 73 |estimate| of cons|
|00004510| 74 61 6e 74 20 74 65 72 | 6d 0a 58 20 20 20 72 65 |tant ter|m.X re|
|00004520| 61 6c 20 61 20 3d 20 6d | 20 2d 20 72 65 61 6c 28 |al a = m| - real(|
|00004530| 4d 20 2a 20 41 29 3b 0a | 58 0a 58 20 20 20 2f 2f |M * A);.|X.X //|
|00004540| 20 47 65 74 20 76 61 72 | 69 61 6e 63 65 73 20 6f | Get var|iances o|
|00004550| 66 20 65 73 74 69 6d 61 | 74 65 73 20 66 72 6f 6d |f estima|tes from|
|00004560| 20 64 69 61 67 6f 6e 61 | 6c 20 65 6c 65 6d 65 6e | diagona|l elemen|
|00004570| 74 73 20 6f 66 20 69 6e | 76 6f 69 63 65 20 6f 66 |ts of in|voice of|
|00004580| 20 53 53 51 0a 58 20 20 | 20 44 69 61 67 6f 6e 61 | SSQ.X | Diagona|
|00004590| 6c 4d 61 74 72 69 78 20 | 44 3b 20 44 20 3c 3c 20 |lMatrix |D; D << |
|000045a0| 4c 2e 74 28 29 2e 69 28 | 29 20 2a 20 4c 2e 69 28 |L.t().i(|) * L.i(|
|000045b0| 29 3b 0a 58 20 20 20 43 | 6f 6c 75 6d 6e 56 65 63 |);.X C|olumnVec|
|000045c0| 74 6f 72 20 56 20 3d 20 | 44 2e 43 6f 70 79 54 6f |tor V = |D.CopyTo|
|000045d0| 43 6f 6c 75 6d 6e 28 29 | 3b 0a 58 20 20 20 72 65 |Column()|;.X re|
|000045e0| 61 6c 20 76 20 3d 20 31 | 2e 30 2f 6e 6f 62 73 20 |al v = 1|.0/nobs |
|000045f0| 2b 20 28 4c 2e 69 28 29 | 20 2a 20 4d 2e 74 28 29 |+ (L.i()| * M.t()|
|00004600| 29 2e 53 75 6d 53 71 75 | 61 72 65 28 29 3b 0a 58 |).SumSqu|are();.X|
|00004610| 0a 58 20 20 20 2f 2f 20 | 43 61 6c 63 75 6c 61 74 |.X // |Calculat|
|00004620| 65 20 66 69 74 74 65 64 | 20 76 61 6c 75 65 73 20 |e fitted| values |
|00004630| 61 6e 64 20 72 65 73 69 | 64 75 61 6c 73 0a 58 20 |and resi|duals.X |
|00004640| 20 20 69 6e 74 20 6e 70 | 72 65 64 31 20 3d 20 6e | int np|red1 = n|
|00004650| 70 72 65 64 2b 31 3b 0a | 58 20 20 20 43 6f 6c 75 |pred+1;.|X Colu|
|00004660| 6d 6e 56 65 63 74 6f 72 | 20 46 69 74 74 65 64 20 |mnVector| Fitted |
|00004670| 3d 20 58 20 2a 20 41 20 | 2b 20 61 3b 0a 58 20 20 |= X * A |+ a;.X |
|00004680| 20 43 6f 6c 75 6d 6e 56 | 65 63 74 6f 72 20 52 65 | ColumnV|ector Re|
|00004690| 73 69 64 75 61 6c 20 3d | 20 59 20 2d 20 46 69 74 |sidual =| Y - Fit|
|000046a0| 74 65 64 3b 0a 58 20 20 | 20 72 65 61 6c 20 52 65 |ted;.X | real Re|
|000046b0| 73 56 61 72 20 3d 20 52 | 65 73 69 64 75 61 6c 2e |sVar = R|esidual.|
|000046c0| 53 75 6d 53 71 75 61 72 | 65 28 29 20 2f 20 28 6e |SumSquar|e() / (n|
|000046d0| 6f 62 73 2d 6e 70 72 65 | 64 31 29 3b 0a 58 0a 58 |obs-npre|d1);.X.X|
|000046e0| 20 20 20 2f 2f 20 47 65 | 74 20 64 69 61 67 6f 6e | // Ge|t diagon|
|000046f0| 61 6c 73 20 6f 66 20 48 | 61 74 20 6d 61 74 72 69 |als of H|at matri|
|00004700| 78 20 28 61 6e 20 65 78 | 70 65 6e 73 69 76 65 20 |x (an ex|pensive |
|00004710| 77 61 79 20 6f 66 20 64 | 6f 69 6e 67 20 74 68 69 |way of d|oing thi|
|00004720| 73 29 0a 58 20 20 20 4d | 61 74 72 69 78 20 58 31 |s).X M|atrix X1|
|00004730| 28 6e 6f 62 73 2c 6e 70 | 72 65 64 31 29 3b 20 58 |(nobs,np|red1); X|
|00004740| 31 2e 43 6f 6c 75 6d 6e | 28 31 29 3c 3c 4f 6e 65 |1.Column|(1)<<One|
|00004750| 73 3b 20 58 31 2e 43 6f | 6c 75 6d 6e 73 28 32 2c |s; X1.Co|lumns(2,|
|00004760| 6e 70 72 65 64 31 29 3c | 3c 58 3b 0a 58 20 20 20 |npred1)<|<X;.X |
|00004770| 44 69 61 67 6f 6e 61 6c | 4d 61 74 72 69 78 20 48 |Diagonal|Matrix H|
|00004780| 61 74 3b 20 20 48 61 74 | 20 3c 3c 20 58 31 20 2a |at; Hat| << X1 *|
|00004790| 20 28 58 31 2e 74 28 29 | 20 2a 20 58 31 29 2e 69 | (X1.t()| * X1).i|
|000047a0| 28 29 20 2a 20 58 31 2e | 74 28 29 3b 0a 58 0a 58 |() * X1.|t();.X.X|
|000047b0| 20 20 20 2f 2f 20 70 72 | 69 6e 74 20 6f 75 74 20 | // pr|int out |
|000047c0| 61 6e 73 77 65 72 73 0a | 58 20 20 20 63 6f 75 74 |answers.|X cout|
|000047d0| 20 3c 3c 20 22 5c 6e 45 | 73 74 69 6d 61 74 65 73 | << "\nE|stimates|
|000047e0| 20 61 6e 64 20 74 68 65 | 69 72 20 73 74 61 6e 64 | and the|ir stand|
|000047f0| 61 72 64 20 65 72 72 6f | 72 73 5c 6e 5c 6e 22 3b |ard erro|rs\n\n";|
|00004800| 0a 58 20 20 20 63 6f 75 | 74 20 3c 3c 20 61 20 3c |.X cou|t << a <|
|00004810| 3c 22 5c 74 22 3c 3c 20 | 73 71 72 74 28 76 2a 52 |<"\t"<< |sqrt(v*R|
|00004820| 65 73 56 61 72 29 20 3c | 3c 20 22 5c 6e 22 3b 0a |esVar) <|< "\n";.|
|00004830| 58 20 20 20 66 6f 72 20 | 28 69 6e 74 20 69 3d 31 |X for |(int i=1|
|00004840| 3b 20 69 3c 3d 6e 70 72 | 65 64 3b 20 69 2b 2b 29 |; i<=npr|ed; i++)|
|00004850| 0a 58 20 20 20 20 20 20 | 63 6f 75 74 20 3c 3c 20 |.X |cout << |
|00004860| 41 28 69 29 20 3c 3c 22 | 5c 74 22 3c 3c 20 73 71 |A(i) <<"|\t"<< sq|
|00004870| 72 74 28 56 28 69 29 2a | 52 65 73 56 61 72 29 20 |rt(V(i)*|ResVar) |
|00004880| 3c 3c 20 22 5c 6e 22 3b | 0a 58 20 20 20 63 6f 75 |<< "\n";|.X cou|
|00004890| 74 20 3c 3c 20 22 5c 6e | 4f 62 73 65 72 76 61 74 |t << "\n|Observat|
|000048a0| 69 6f 6e 73 2c 20 66 69 | 74 74 65 64 20 76 61 6c |ions, fi|tted val|
|000048b0| 75 65 2c 20 72 65 73 69 | 64 75 61 6c 20 76 61 6c |ue, resi|dual val|
|000048c0| 75 65 2c 20 68 61 74 20 | 76 61 6c 75 65 5c 6e 22 |ue, hat |value\n"|
|000048d0| 3b 0a 58 20 20 20 66 6f | 72 20 28 69 3d 31 3b 20 |;.X fo|r (i=1; |
|000048e0| 69 3c 3d 6e 6f 62 73 3b | 20 69 2b 2b 29 0a 58 20 |i<=nobs;| i++).X |
|000048f0| 20 20 20 20 20 63 6f 75 | 74 20 3c 3c 20 58 28 69 | cou|t << X(i|
|00004900| 2c 31 29 20 3c 3c 22 5c | 74 22 3c 3c 20 58 28 69 |,1) <<"\|t"<< X(i|
|00004910| 2c 32 29 20 3c 3c 22 5c | 74 22 3c 3c 20 59 28 69 |,2) <<"\|t"<< Y(i|
|00004920| 29 20 3c 3c 22 5c 74 22 | 3c 3c 0a 58 20 20 20 20 |) <<"\t"|<<.X |
|00004930| 20 20 74 33 28 46 69 74 | 74 65 64 28 69 29 29 20 | t3(Fit|ted(i)) |
|00004940| 3c 3c 22 5c 74 22 3c 3c | 20 74 33 28 52 65 73 69 |<<"\t"<<| t3(Resi|
|00004950| 64 75 61 6c 28 69 29 29 | 20 3c 3c 22 5c 74 22 3c |dual(i))| <<"\t"<|
|00004960| 3c 20 74 33 28 48 61 74 | 28 69 29 29 20 3c 3c 22 |< t3(Hat|(i)) <<"|
|00004970| 5c 6e 22 3b 0a 58 20 20 | 20 63 6f 75 74 20 3c 3c |\n";.X | cout <<|
|00004980| 20 22 5c 6e 5c 6e 22 3b | 0a 58 7d 0a 58 0a 58 76 | "\n\n";|.X}.X.Xv|
|00004990| 6f 69 64 20 74 65 73 74 | 33 28 72 65 61 6c 2a 20 |oid test|3(real* |
|000049a0| 79 2c 20 72 65 61 6c 2a | 20 78 31 2c 20 72 65 61 |y, real*| x1, rea|
|000049b0| 6c 2a 20 78 32 2c 20 69 | 6e 74 20 6e 6f 62 73 2c |l* x2, i|nt nobs,|
|000049c0| 20 69 6e 74 20 6e 70 72 | 65 64 29 0a 58 7b 0a 58 | int npr|ed).X{.X|
|000049d0| 20 20 20 63 6f 75 74 20 | 3c 3c 20 22 5c 6e 5c 6e | cout |<< "\n\n|
|000049e0| 54 65 73 74 20 33 20 2d | 20 48 6f 75 73 65 68 6f |Test 3 -| Househo|
|000049f0| 6c 64 65 72 20 74 72 69 | 61 6e 67 75 6c 61 72 69 |lder tri|angulari|
|00004a00| 73 61 74 69 6f 6e 5c 6e | 22 3b 0a 58 0a 58 20 20 |sation\n|";.X.X |
|00004a10| 20 2f 2f 20 48 6f 75 73 | 65 68 6f 6c 64 65 72 20 | // Hous|eholder |
|00004a20| 74 72 69 61 6e 67 75 6c | 61 72 69 73 61 74 69 6f |triangul|arisatio|
|00004a30| 6e 20 6d 65 74 68 6f 64 | 0a 58 20 0a 58 20 20 20 |n method|.X .X |
|00004a40| 2f 2f 20 6c 6f 61 64 20 | 64 61 74 61 20 2d 20 31 |// load |data - 1|
|00004a50| 73 20 69 6e 74 6f 20 63 | 6f 6c 20 31 20 6f 66 20 |s into c|ol 1 of |
|00004a60| 6d 61 74 72 69 78 0a 58 | 20 20 20 69 6e 74 20 6e |matrix.X| int n|
|00004a70| 70 72 65 64 31 20 3d 20 | 6e 70 72 65 64 2b 31 3b |pred1 = |npred+1;|
|00004a80| 0a 58 20 20 20 4d 61 74 | 72 69 78 20 58 28 6e 6f |.X Mat|rix X(no|
|00004a90| 62 73 2c 6e 70 72 65 64 | 31 29 3b 20 43 6f 6c 75 |bs,npred|1); Colu|
|00004aa0| 6d 6e 56 65 63 74 6f 72 | 20 59 28 6e 6f 62 73 29 |mnVector| Y(nobs)|
|00004ab0| 3b 0a 58 20 20 20 58 2e | 43 6f 6c 75 6d 6e 28 31 |;.X X.|Column(1|
|00004ac0| 29 20 3c 3c 20 31 2e 30 | 3b 20 20 58 2e 43 6f 6c |) << 1.0|; X.Col|
|00004ad0| 75 6d 6e 28 32 29 20 3c | 3c 20 78 31 3b 20 20 58 |umn(2) <|< x1; X|
|00004ae0| 2e 43 6f 6c 75 6d 6e 28 | 33 29 20 3c 3c 20 78 32 |.Column(|3) << x2|
|00004af0| 3b 20 20 59 20 3c 3c 20 | 79 3b 0a 58 0a 58 20 20 |; Y << |y;.X.X |
|00004b00| 20 2f 2f 20 64 6f 20 48 | 6f 75 73 65 68 6f 6c 64 | // do H|ousehold|
|00004b10| 65 72 20 74 72 69 61 6e | 67 75 6c 61 72 69 73 61 |er trian|gularisa|
|00004b20| 74 69 6f 6e 0a 58 20 20 | 20 2f 2f 20 6e 6f 20 6e |tion.X | // no n|
|00004b30| 65 65 64 20 74 6f 20 64 | 65 61 6c 20 77 69 74 68 |eed to d|eal with|
|00004b40| 20 63 6f 6e 73 74 61 6e | 74 20 74 65 72 6d 20 73 | constan|t term s|
|00004b50| 65 70 61 72 61 74 65 6c | 79 0a 58 20 20 20 4d 61 |eparatel|y.X Ma|
|00004b60| 74 72 69 78 20 58 54 20 | 3d 20 58 2e 74 28 29 3b |trix XT |= X.t();|
|00004b70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2f 2f 20 | | // |
|00004b80| 57 61 6e 74 20 64 61 74 | 61 20 74 6f 20 62 65 20 |Want dat|a to be |
|00004b90| 61 6c 6f 6e 67 20 72 6f | 77 73 0a 58 20 20 20 52 |along ro|ws.X R|
|00004ba0| 6f 77 56 65 63 74 6f 72 | 20 59 54 20 3d 20 59 2e |owVector| YT = Y.|
|00004bb0| 74 28 29 3b 0a 58 20 20 | 20 4c 6f 77 65 72 54 72 |t();.X | LowerTr|
|00004bc0| 69 61 6e 67 75 6c 61 72 | 4d 61 74 72 69 78 20 4c |iangular|Matrix L|
|00004bd0| 3b 20 52 6f 77 56 65 63 | 74 6f 72 20 4d 3b 0a 58 |; RowVec|tor M;.X|
|00004be0| 20 20 20 48 48 44 65 63 | 6f 6d 70 6f 73 65 28 58 | HHDec|ompose(X|
|00004bf0| 54 2c 20 4c 29 3b 20 48 | 48 44 65 63 6f 6d 70 6f |T, L); H|HDecompo|
|00004c00| 73 65 28 58 54 2c 20 59 | 54 2c 20 4d 29 3b 20 2f |se(XT, Y|T, M); /|
|00004c10| 2f 20 59 54 20 6e 6f 77 | 20 63 6f 6e 74 61 69 6e |/ YT now| contain|
|00004c20| 73 20 72 65 73 69 64 73 | 0a 58 20 20 20 43 6f 6c |s resids|.X Col|
|00004c30| 75 6d 6e 56 65 63 74 6f | 72 20 41 20 3d 20 4c 2e |umnVecto|r A = L.|
|00004c40| 74 28 29 2e 69 28 29 20 | 2a 20 4d 2e 74 28 29 3b |t().i() |* M.t();|
|00004c50| 0a 58 20 20 20 43 6f 6c | 75 6d 6e 56 65 63 74 6f |.X Col|umnVecto|
|00004c60| 72 20 46 69 74 74 65 64 | 20 3d 20 58 20 2a 20 41 |r Fitted| = X * A|
|00004c70| 3b 0a 58 20 20 20 72 65 | 61 6c 20 52 65 73 56 61 |;.X re|al ResVa|
|00004c80| 72 20 3d 20 59 54 2e 53 | 75 6d 53 71 75 61 72 65 |r = YT.S|umSquare|
|00004c90| 28 29 20 2f 20 28 6e 6f | 62 73 2d 6e 70 72 65 64 |() / (no|bs-npred|
|00004ca0| 31 29 3b 0a 58 0a 58 20 | 20 20 2f 2f 20 67 65 74 |1);.X.X | // get|
|00004cb0| 20 76 61 72 69 61 6e 63 | 65 73 20 6f 66 20 65 73 | varianc|es of es|
|00004cc0| 74 69 6d 61 74 65 73 0a | 58 20 20 20 4c 20 3d 20 |timates.|X L = |
|00004cd0| 4c 2e 69 28 29 3b 20 44 | 69 61 67 6f 6e 61 6c 4d |L.i(); D|iagonalM|
|00004ce0| 61 74 72 69 78 20 44 3b | 20 44 20 3c 3c 20 4c 2e |atrix D;| D << L.|
|00004cf0| 74 28 29 20 2a 20 4c 3b | 0a 58 0a 58 20 20 20 2f |t() * L;|.X.X /|
|00004d00| 2f 20 47 65 74 20 64 69 | 61 67 6f 6e 61 6c 73 20 |/ Get di|agonals |
|00004d10| 6f 66 20 48 61 74 20 6d | 61 74 72 69 78 0a 58 20 |of Hat m|atrix.X |
|00004d20| 20 20 44 69 61 67 6f 6e | 61 6c 4d 61 74 72 69 78 | Diagon|alMatrix|
|00004d30| 20 48 61 74 3b 20 20 48 | 61 74 20 3c 3c 20 58 54 | Hat; H|at << XT|
|00004d40| 2e 74 28 29 20 2a 20 58 | 54 3b 0a 58 0a 58 20 20 |.t() * X|T;.X.X |
|00004d50| 20 2f 2f 20 70 72 69 6e | 74 20 6f 75 74 20 61 6e | // prin|t out an|
|00004d60| 73 77 65 72 73 0a 58 20 | 20 20 63 6f 75 74 20 3c |swers.X | cout <|
|00004d70| 3c 20 22 5c 6e 45 73 74 | 69 6d 61 74 65 73 20 61 |< "\nEst|imates a|
|00004d80| 6e 64 20 74 68 65 69 72 | 20 73 74 61 6e 64 61 72 |nd their| standar|
|00004d90| 64 20 65 72 72 6f 72 73 | 5c 6e 5c 6e 22 3b 0a 58 |d errors|\n\n";.X|
|00004da0| 20 20 20 66 6f 72 20 28 | 69 6e 74 20 69 3d 31 3b | for (|int i=1;|
|00004db0| 20 69 3c 3d 6e 70 72 65 | 64 31 3b 20 69 2b 2b 29 | i<=npre|d1; i++)|
|00004dc0| 0a 58 20 20 20 20 20 20 | 63 6f 75 74 20 3c 3c 20 |.X |cout << |
|00004dd0| 41 28 69 29 20 3c 3c 22 | 5c 74 22 3c 3c 20 73 71 |A(i) <<"|\t"<< sq|
|00004de0| 72 74 28 44 28 69 29 2a | 52 65 73 56 61 72 29 20 |rt(D(i)*|ResVar) |
|00004df0| 3c 3c 20 22 5c 6e 22 3b | 0a 58 20 20 20 63 6f 75 |<< "\n";|.X cou|
|00004e00| 74 20 3c 3c 20 22 5c 6e | 4f 62 73 65 72 76 61 74 |t << "\n|Observat|
|00004e10| 69 6f 6e 73 2c 20 66 69 | 74 74 65 64 20 76 61 6c |ions, fi|tted val|
|00004e20| 75 65 2c 20 72 65 73 69 | 64 75 61 6c 20 76 61 6c |ue, resi|dual val|
|00004e30| 75 65 2c 20 68 61 74 20 | 76 61 6c 75 65 5c 6e 22 |ue, hat |value\n"|
|00004e40| 3b 0a 58 20 20 20 66 6f | 72 20 28 69 3d 31 3b 20 |;.X fo|r (i=1; |
|00004e50| 69 3c 3d 6e 6f 62 73 3b | 20 69 2b 2b 29 0a 58 20 |i<=nobs;| i++).X |
|00004e60| 20 20 20 20 20 63 6f 75 | 74 20 3c 3c 20 58 28 69 | cou|t << X(i|
|00004e70| 2c 32 29 20 3c 3c 22 5c | 74 22 3c 3c 20 58 28 69 |,2) <<"\|t"<< X(i|
|00004e80| 2c 33 29 20 3c 3c 22 5c | 74 22 3c 3c 20 59 28 69 |,3) <<"\|t"<< Y(i|
|00004e90| 29 20 3c 3c 22 5c 74 22 | 3c 3c 0a 58 20 20 20 20 |) <<"\t"|<<.X |
|00004ea0| 20 20 74 33 28 46 69 74 | 74 65 64 28 69 29 29 20 | t3(Fit|ted(i)) |
|00004eb0| 3c 3c 22 5c 74 22 3c 3c | 20 74 33 28 59 54 28 69 |<<"\t"<<| t3(YT(i|
|00004ec0| 29 29 20 3c 3c 22 5c 74 | 22 3c 3c 20 74 33 28 48 |)) <<"\t|"<< t3(H|
|00004ed0| 61 74 28 69 29 29 20 3c | 3c 22 5c 6e 22 3b 0a 58 |at(i)) <|<"\n";.X|
|00004ee0| 20 20 20 63 6f 75 74 20 | 3c 3c 20 22 5c 6e 5c 6e | cout |<< "\n\n|
|00004ef0| 22 3b 0a 58 7d 0a 58 0a | 58 76 6f 69 64 20 74 65 |";.X}.X.|Xvoid te|
|00004f00| 73 74 34 28 72 65 61 6c | 2a 20 79 2c 20 72 65 61 |st4(real|* y, rea|
|00004f10| 6c 2a 20 78 31 2c 20 72 | 65 61 6c 2a 20 78 32 2c |l* x1, r|eal* x2,|
|00004f20| 20 69 6e 74 20 6e 6f 62 | 73 2c 20 69 6e 74 20 6e | int nob|s, int n|
|00004f30| 70 72 65 64 29 0a 58 7b | 0a 58 20 20 20 63 6f 75 |pred).X{|.X cou|
|00004f40| 74 20 3c 3c 20 22 5c 6e | 5c 6e 54 65 73 74 20 34 |t << "\n|\nTest 4|
|00004f50| 20 2d 20 73 69 6e 67 75 | 6c 61 72 20 76 61 6c 75 | - singu|lar valu|
|00004f60| 65 5c 6e 22 3b 0a 58 0a | 58 20 20 20 2f 2f 20 53 |e\n";.X.|X // S|
|00004f70| 69 6e 67 75 6c 61 72 20 | 76 61 6c 75 65 20 64 65 |ingular |value de|
|00004f80| 63 6f 6d 70 6f 73 69 74 | 69 6f 6e 20 6d 65 74 68 |composit|ion meth|
|00004f90| 6f 64 0a 58 20 0a 58 20 | 20 20 2f 2f 20 6c 6f 61 |od.X .X | // loa|
|00004fa0| 64 20 64 61 74 61 20 2d | 20 31 73 20 69 6e 74 6f |d data -| 1s into|
|00004fb0| 20 63 6f 6c 20 31 20 6f | 66 20 6d 61 74 72 69 78 | col 1 o|f matrix|
|00004fc0| 0a 58 20 20 20 69 6e 74 | 20 6e 70 72 65 64 31 20 |.X int| npred1 |
|00004fd0| 3d 20 6e 70 72 65 64 2b | 31 3b 0a 58 20 20 20 4d |= npred+|1;.X M|
|00004fe0| 61 74 72 69 78 20 58 28 | 6e 6f 62 73 2c 6e 70 72 |atrix X(|nobs,npr|
|00004ff0| 65 64 31 29 3b 20 43 6f | 6c 75 6d 6e 56 65 63 74 |ed1); Co|lumnVect|
|00005000| 6f 72 20 59 28 6e 6f 62 | 73 29 3b 0a 58 20 20 20 |or Y(nob|s);.X |
|00005010| 58 2e 43 6f 6c 75 6d 6e | 28 31 29 20 3c 3c 20 31 |X.Column|(1) << 1|
|00005020| 2e 30 3b 20 20 58 2e 43 | 6f 6c 75 6d 6e 28 32 29 |.0; X.C|olumn(2)|
|00005030| 20 3c 3c 20 78 31 3b 20 | 20 58 2e 43 6f 6c 75 6d | << x1; | X.Colum|
|00005040| 6e 28 33 29 20 3c 3c 20 | 78 32 3b 20 20 59 20 3c |n(3) << |x2; Y <|
|00005050| 3c 20 79 3b 0a 58 0a 58 | 20 20 20 2f 2f 20 64 6f |< y;.X.X| // do|
|00005060| 20 53 56 44 0a 58 20 20 | 20 4d 61 74 72 69 78 20 | SVD.X | Matrix |
|00005070| 55 2c 20 56 3b 20 44 69 | 61 67 6f 6e 61 6c 4d 61 |U, V; Di|agonalMa|
|00005080| 74 72 69 78 20 44 3b 0a | 58 20 20 20 53 56 44 28 |trix D;.|X SVD(|
|00005090| 58 2c 44 2c 55 2c 56 29 | 3b 20 20 20 20 20 20 20 |X,D,U,V)|; |
|000050a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000050b0| 20 20 20 20 20 20 20 2f | 2f 20 58 20 3d 20 55 20 | /|/ X = U |
|000050c0| 2a 20 44 20 2a 20 56 2e | 74 28 29 0a 58 20 20 20 |* D * V.|t().X |
|000050d0| 43 6f 6c 75 6d 6e 56 65 | 63 74 6f 72 20 46 69 74 |ColumnVe|ctor Fit|
|000050e0| 74 65 64 20 3d 20 55 2e | 74 28 29 20 2a 20 59 3b |ted = U.|t() * Y;|
|000050f0| 0a 58 20 20 20 43 6f 6c | 75 6d 6e 56 65 63 74 6f |.X Col|umnVecto|
|00005100| 72 20 41 20 3d 20 56 20 | 2a 20 28 20 44 2e 69 28 |r A = V |* ( D.i(|
|00005110| 29 20 2a 20 46 69 74 74 | 65 64 20 29 3b 0a 58 20 |) * Fitt|ed );.X |
|00005120| 20 20 46 69 74 74 65 64 | 20 3d 20 55 20 2a 20 46 | Fitted| = U * F|
|00005130| 69 74 74 65 64 3b 0a 58 | 20 20 20 43 6f 6c 75 6d |itted;.X| Colum|
|00005140| 6e 56 65 63 74 6f 72 20 | 52 65 73 69 64 75 61 6c |nVector |Residual|
|00005150| 20 3d 20 59 20 2d 20 46 | 69 74 74 65 64 3b 0a 58 | = Y - F|itted;.X|
|00005160| 20 20 20 72 65 61 6c 20 | 52 65 73 56 61 72 20 3d | real |ResVar =|
|00005170| 20 52 65 73 69 64 75 61 | 6c 2e 53 75 6d 53 71 75 | Residua|l.SumSqu|
|00005180| 61 72 65 28 29 20 2f 20 | 28 6e 6f 62 73 2d 6e 70 |are() / |(nobs-np|
|00005190| 72 65 64 31 29 3b 0a 58 | 0a 58 20 20 20 2f 2f 20 |red1);.X|.X // |
|000051a0| 67 65 74 20 76 61 72 69 | 61 6e 63 65 73 20 6f 66 |get vari|ances of|
|000051b0| 20 65 73 74 69 6d 61 74 | 65 73 0a 58 20 20 20 44 | estimat|es.X D|
|000051c0| 20 3c 3c 20 56 20 2a 20 | 28 44 20 2a 20 44 29 2e | << V * |(D * D).|
|000051d0| 69 28 29 20 2a 20 56 2e | 74 28 29 3b 0a 58 0a 58 |i() * V.|t();.X.X|
|000051e0| 20 20 20 2f 2f 20 47 65 | 74 20 64 69 61 67 6f 6e | // Ge|t diagon|
|000051f0| 61 6c 73 20 6f 66 20 48 | 61 74 20 6d 61 74 72 69 |als of H|at matri|
|00005200| 78 0a 58 20 20 20 44 69 | 61 67 6f 6e 61 6c 4d 61 |x.X Di|agonalMa|
|00005210| 74 72 69 78 20 48 61 74 | 3b 20 20 48 61 74 20 3c |trix Hat|; Hat <|
|00005220| 3c 20 55 20 2a 20 55 2e | 74 28 29 3b 0a 58 0a 58 |< U * U.|t();.X.X|
|00005230| 20 20 20 2f 2f 20 70 72 | 69 6e 74 20 6f 75 74 20 | // pr|int out |
|00005240| 61 6e 73 77 65 72 73 0a | 58 20 20 20 63 6f 75 74 |answers.|X cout|
|00005250| 20 3c 3c 20 22 5c 6e 45 | 73 74 69 6d 61 74 65 73 | << "\nE|stimates|
|00005260| 20 61 6e 64 20 74 68 65 | 69 72 20 73 74 61 6e 64 | and the|ir stand|
|00005270| 61 72 64 20 65 72 72 6f | 72 73 5c 6e 5c 6e 22 3b |ard erro|rs\n\n";|
|00005280| 0a 58 20 20 20 66 6f 72 | 20 28 69 6e 74 20 69 3d |.X for| (int i=|
|00005290| 31 3b 20 69 3c 3d 6e 70 | 72 65 64 31 3b 20 69 2b |1; i<=np|red1; i+|
|000052a0| 2b 29 0a 58 20 20 20 20 | 20 20 63 6f 75 74 20 3c |+).X | cout <|
|000052b0| 3c 20 41 28 69 29 20 3c | 3c 22 5c 74 22 3c 3c 20 |< A(i) <|<"\t"<< |
|000052c0| 73 71 72 74 28 44 28 69 | 29 2a 52 65 73 56 61 72 |sqrt(D(i|)*ResVar|
|000052d0| 29 20 3c 3c 20 22 5c 6e | 22 3b 0a 58 20 20 20 63 |) << "\n|";.X c|
|000052e0| 6f 75 74 20 3c 3c 20 22 | 5c 6e 4f 62 73 65 72 76 |out << "|\nObserv|
|000052f0| 61 74 69 6f 6e 73 2c 20 | 66 69 74 74 65 64 20 76 |ations, |fitted v|
|00005300| 61 6c 75 65 2c 20 72 65 | 73 69 64 75 61 6c 20 76 |alue, re|sidual v|
|00005310| 61 6c 75 65 2c 20 68 61 | 74 20 76 61 6c 75 65 5c |alue, ha|t value\|
|00005320| 6e 22 3b 0a 58 20 20 20 | 66 6f 72 20 28 69 3d 31 |n";.X |for (i=1|
|00005330| 3b 20 69 3c 3d 6e 6f 62 | 73 3b 20 69 2b 2b 29 0a |; i<=nob|s; i++).|
|00005340| 58 20 20 20 20 20 20 63 | 6f 75 74 20 3c 3c 20 58 |X c|out << X|
|00005350| 28 69 2c 32 29 20 3c 3c | 22 5c 74 22 3c 3c 20 58 |(i,2) <<|"\t"<< X|
|00005360| 28 69 2c 33 29 20 3c 3c | 22 5c 74 22 3c 3c 20 59 |(i,3) <<|"\t"<< Y|
|00005370| 28 69 29 20 3c 3c 22 5c | 74 22 3c 3c 0a 58 20 20 |(i) <<"\|t"<<.X |
|00005380| 20 20 20 20 74 33 28 46 | 69 74 74 65 64 28 69 29 | t3(F|itted(i)|
|00005390| 29 20 3c 3c 22 5c 74 22 | 3c 3c 20 74 33 28 52 65 |) <<"\t"|<< t3(Re|
|000053a0| 73 69 64 75 61 6c 28 69 | 29 29 20 3c 3c 22 5c 74 |sidual(i|)) <<"\t|
|000053b0| 22 3c 3c 20 74 33 28 48 | 61 74 28 69 29 29 20 3c |"<< t3(H|at(i)) <|
|000053c0| 3c 22 5c 6e 22 3b 0a 58 | 20 20 20 63 6f 75 74 20 |<"\n";.X| cout |
|000053d0| 3c 3c 20 22 5c 6e 5c 6e | 22 3b 0a 58 7d 0a 58 0a |<< "\n\n|";.X}.X.|
|000053e0| 58 6d 61 69 6e 28 29 0a | 58 7b 0a 58 20 20 20 63 |Xmain().|X{.X c|
|000053f0| 6f 75 74 20 3c 3c 20 22 | 5c 6e 44 65 6d 6f 6e 73 |out << "|\nDemons|
|00005400| 74 72 61 74 69 6f 6e 20 | 6f 66 20 4d 61 74 72 69 |tration |of Matri|
|00005410| 78 20 70 61 63 6b 61 67 | 65 5c 6e 5c 6e 22 3b 0a |x packag|e\n\n";.|
|00005420| 58 0a 58 20 20 20 2f 2f | 20 54 65 73 74 20 66 6f |X.X //| Test fo|
|00005430| 72 20 61 6e 79 20 6d 65 | 6d 6f 72 79 20 6e 6f 74 |r any me|mory not|
|00005440| 20 64 65 61 6c 6c 6f 63 | 61 74 65 64 20 61 66 74 | dealloc|ated aft|
|00005450| 65 72 20 72 75 6e 6e 69 | 6e 67 20 74 68 69 73 20 |er runni|ng this |
|00005460| 70 72 6f 67 72 61 6d 0a | 58 20 20 20 72 65 61 6c |program.|X real|
|00005470| 2a 20 73 31 3b 20 7b 20 | 43 6f 6c 75 6d 6e 56 65 |* s1; { |ColumnVe|
|00005480| 63 74 6f 72 20 41 28 38 | 30 30 30 29 3b 20 73 31 |ctor A(8|000); s1|
|00005490| 20 3d 20 41 2e 53 74 6f | 72 65 28 29 3b 20 7d 0a | = A.Sto|re(); }.|
|000054a0| 58 0a 58 20 20 20 7b 0a | 58 20 20 20 20 20 20 2f |X.X {.|X /|
|000054b0| 2f 20 74 68 65 20 64 61 | 74 61 0a 58 20 20 20 20 |/ the da|ta.X |
|000054c0| 20 20 2f 2f 20 79 6f 75 | 20 6d 61 79 20 6e 65 65 | // you| may nee|
|000054d0| 64 20 74 6f 20 72 65 61 | 64 20 74 68 69 73 20 64 |d to rea|d this d|
|000054e0| 61 74 61 20 75 73 69 6e | 67 20 63 69 6e 20 69 66 |ata usin|g cin if|
|000054f0| 20 79 6f 75 20 61 72 65 | 20 75 73 69 6e 67 20 61 | you are| using a|
|00005500| 0a 58 20 20 20 20 20 20 | 2f 2f 20 63 6f 6d 70 69 |.X |// compi|
|00005510| 6c 65 72 20 74 68 61 74 | 20 64 6f 65 73 6e 27 74 |ler that| doesn't|
|00005520| 20 75 6e 64 65 72 73 74 | 61 6e 64 20 61 67 67 72 | underst|and aggr|
|00005530| 65 67 61 74 65 73 0a 58 | 23 69 66 6e 64 65 66 20 |egates.X|#ifndef |
|00005540| 41 54 61 6e 64 54 0a 58 | 20 20 20 20 20 20 72 65 |ATandT.X| re|
|00005550| 61 6c 20 79 5b 5d 20 20 | 3d 20 7b 20 38 2e 33 2c |al y[] |= { 8.3,|
|00005560| 20 35 2e 35 2c 20 38 2e | 30 2c 20 38 2e 35 2c 20 | 5.5, 8.|0, 8.5, |
|00005570| 35 2e 37 2c 20 34 2e 34 | 2c 20 36 2e 33 2c 20 37 |5.7, 4.4|, 6.3, 7|
|00005580| 2e 39 2c 20 39 2e 31 20 | 7d 3b 0a 58 20 20 20 20 |.9, 9.1 |};.X |
|00005590| 20 20 72 65 61 6c 20 78 | 31 5b 5d 20 3d 20 7b 20 | real x|1[] = { |
|000055a0| 32 2e 34 2c 20 31 2e 38 | 2c 20 32 2e 34 2c 20 33 |2.4, 1.8|, 2.4, 3|
|000055b0| 2e 30 2c 20 32 2e 30 2c | 20 31 2e 32 2c 20 32 2e |.0, 2.0,| 1.2, 2.|
|000055c0| 30 2c 20 32 2e 37 2c 20 | 33 2e 36 20 7d 3b 0a 58 |0, 2.7, |3.6 };.X|
|000055d0| 20 20 20 20 20 20 72 65 | 61 6c 20 78 32 5b 5d 20 | re|al x2[] |
|000055e0| 3d 20 7b 20 31 2e 37 2c | 20 30 2e 39 2c 20 31 2e |= { 1.7,| 0.9, 1.|
|000055f0| 36 2c 20 31 2e 39 2c 20 | 30 2e 35 2c 20 30 2e 36 |6, 1.9, |0.5, 0.6|
|00005600| 2c 20 31 2e 31 2c 20 31 | 2e 30 2c 20 30 2e 35 20 |, 1.1, 1|.0, 0.5 |
|00005610| 7d 3b 0a 58 23 65 6c 73 | 65 0a 58 20 20 20 20 20 |};.X#els|e.X |
|00005620| 20 72 65 61 6c 20 79 5b | 39 5d 2c 20 78 31 5b 39 | real y[|9], x1[9|
|00005630| 5d 2c 20 78 32 5b 39 5d | 3b 0a 58 20 20 20 20 20 |], x2[9]|;.X |
|00005640| 20 79 5b 30 5d 3d 38 2e | 33 3b 20 79 5b 31 5d 3d | y[0]=8.|3; y[1]=|
|00005650| 35 2e 35 3b 20 79 5b 32 | 5d 3d 38 2e 30 3b 20 79 |5.5; y[2|]=8.0; y|
|00005660| 5b 33 5d 3d 38 2e 35 3b | 20 79 5b 34 5d 3d 35 2e |[3]=8.5;| y[4]=5.|
|00005670| 37 3b 0a 58 20 20 20 20 | 20 20 79 5b 35 5d 3d 34 |7;.X | y[5]=4|
|00005680| 2e 34 3b 20 79 5b 36 5d | 3d 36 2e 33 3b 20 79 5b |.4; y[6]|=6.3; y[|
|00005690| 37 5d 3d 37 2e 39 3b 20 | 79 5b 38 5d 3d 39 2e 31 |7]=7.9; |y[8]=9.1|
|000056a0| 3b 0a 58 20 20 20 20 20 | 20 78 31 5b 30 5d 3d 32 |;.X | x1[0]=2|
|000056b0| 2e 34 3b 20 78 31 5b 31 | 5d 3d 31 2e 38 3b 20 78 |.4; x1[1|]=1.8; x|
|000056c0| 31 5b 32 5d 3d 32 2e 34 | 3b 20 78 31 5b 33 5d 3d |1[2]=2.4|; x1[3]=|
|000056d0| 33 2e 30 3b 20 78 31 5b | 34 5d 3d 32 2e 30 3b 0a |3.0; x1[|4]=2.0;.|
|000056e0| 58 20 20 20 20 20 20 78 | 31 5b 35 5d 3d 31 2e 32 |X x|1[5]=1.2|
|000056f0| 3b 20 78 31 5b 36 5d 3d | 32 2e 30 3b 20 78 31 5b |; x1[6]=|2.0; x1[|
|00005700| 37 5d 3d 32 2e 37 3b 20 | 78 31 5b 38 5d 3d 33 2e |7]=2.7; |x1[8]=3.|
|00005710| 36 3b 0a 58 20 20 20 20 | 20 20 78 32 5b 30 5d 3d |6;.X | x2[0]=|
|00005720| 31 2e 37 3b 20 78 32 5b | 31 5d 3d 30 2e 39 3b 20 |1.7; x2[|1]=0.9; |
|00005730| 78 32 5b 32 5d 3d 31 2e | 36 3b 20 78 32 5b 33 5d |x2[2]=1.|6; x2[3]|
|00005740| 3d 31 2e 39 3b 20 78 32 | 5b 34 5d 3d 30 2e 35 3b |=1.9; x2|[4]=0.5;|
|00005750| 0a 58 20 20 20 20 20 20 | 78 32 5b 35 5d 3d 30 2e |.X |x2[5]=0.|
|00005760| 36 3b 20 78 32 5b 36 5d | 3d 31 2e 31 3b 20 78 32 |6; x2[6]|=1.1; x2|
|00005770| 5b 37 5d 3d 31 2e 30 3b | 20 78 32 5b 38 5d 3d 30 |[7]=1.0;| x2[8]=0|
|00005780| 2e 35 3b 0a 58 23 65 6e | 64 69 66 0a 58 20 20 20 |.5;.X#en|dif.X |
|00005790| 20 20 20 69 6e 74 20 6e | 6f 62 73 20 3d 20 39 3b | int n|obs = 9;|
|000057a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000057b0| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 6e 75 | | // nu|
|000057c0| 6d 62 65 72 20 6f 66 20 | 6f 62 73 65 72 76 61 74 |mber of |observat|
|000057d0| 69 6f 6e 73 0a 58 20 20 | 20 20 20 20 69 6e 74 20 |ions.X | int |
|000057e0| 6e 70 72 65 64 20 3d 20 | 32 3b 20 20 20 20 20 20 |npred = |2; |
|000057f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005800| 20 20 20 20 2f 2f 20 6e | 75 6d 62 65 72 20 6f 66 | // n|umber of|
|00005810| 20 70 72 65 64 69 63 74 | 6f 72 20 76 61 6c 75 65 | predict|or value|
|00005820| 73 0a 58 0a 58 20 20 20 | 20 20 20 2f 2f 20 77 65 |s.X.X | // we|
|00005830| 20 77 61 6e 74 20 74 6f | 20 66 69 6e 64 20 74 68 | want to| find th|
|00005840| 65 20 76 61 6c 75 65 73 | 20 6f 66 20 61 2c 61 31 |e values| of a,a1|
|00005850| 2c 61 32 20 74 6f 20 67 | 69 76 65 20 74 68 65 20 |,a2 to g|ive the |
|00005860| 62 65 73 74 0a 58 20 20 | 20 20 20 20 2f 2f 20 66 |best.X | // f|
|00005870| 69 74 20 6f 66 20 79 5b | 69 5d 20 77 69 74 68 20 |it of y[|i] with |
|00005880| 61 30 20 2b 20 61 31 2a | 78 31 5b 69 5d 20 2b 20 |a0 + a1*|x1[i] + |
|00005890| 61 32 2a 78 32 5b 69 5d | 0a 58 20 20 20 20 20 20 |a2*x2[i]|.X |
|000058a0| 2f 2f 20 41 6c 73 6f 20 | 70 72 69 6e 74 20 64 69 |// Also |print di|
|000058b0| 61 67 6f 6e 61 6c 20 65 | 6c 65 6d 65 6e 74 73 20 |agonal e|lements |
|000058c0| 6f 66 20 68 61 74 20 6d | 61 74 72 69 78 2c 20 58 |of hat m|atrix, X|
|000058d0| 2a 28 58 2e 74 28 29 2a | 58 29 2e 69 28 29 2a 58 |*(X.t()*|X).i()*X|
|000058e0| 2e 74 28 29 0a 58 0a 58 | 20 20 20 20 20 20 2f 2f |.t().X.X| //|
|000058f0| 20 74 68 69 73 20 65 78 | 61 6d 70 6c 65 20 64 65 | this ex|ample de|
|00005900| 6d 6f 6e 73 74 72 61 74 | 65 73 20 66 6f 75 72 20 |monstrat|es four |
|00005910| 6d 65 74 68 6f 64 73 20 | 6f 66 20 63 61 6c 63 75 |methods |of calcu|
|00005920| 6c 61 74 69 6f 6e 0a 58 | 0a 58 0a 58 20 20 20 20 |lation.X|.X.X |
|00005930| 20 20 74 65 73 74 31 28 | 79 2c 20 78 31 2c 20 78 | test1(|y, x1, x|
|00005940| 32 2c 20 6e 6f 62 73 2c | 20 6e 70 72 65 64 29 3b |2, nobs,| npred);|
|00005950| 0a 58 20 20 20 20 20 20 | 74 65 73 74 32 28 79 2c |.X |test2(y,|
|00005960| 20 78 31 2c 20 78 32 2c | 20 6e 6f 62 73 2c 20 6e | x1, x2,| nobs, n|
|00005970| 70 72 65 64 29 3b 0a 58 | 20 20 20 20 20 20 74 65 |pred);.X| te|
|00005980| 73 74 33 28 79 2c 20 78 | 31 2c 20 78 32 2c 20 6e |st3(y, x|1, x2, n|
|00005990| 6f 62 73 2c 20 6e 70 72 | 65 64 29 3b 0a 58 20 20 |obs, npr|ed);.X |
|000059a0| 20 20 20 20 74 65 73 74 | 34 28 79 2c 20 78 31 2c | test|4(y, x1,|
|000059b0| 20 78 32 2c 20 6e 6f 62 | 73 2c 20 6e 70 72 65 64 | x2, nob|s, npred|
|000059c0| 29 3b 0a 58 20 20 20 7d | 0a 58 0a 58 20 20 20 72 |);.X }|.X.X r|
|000059d0| 65 61 6c 2a 20 73 32 3b | 20 7b 20 43 6f 6c 75 6d |eal* s2;| { Colum|
|000059e0| 6e 56 65 63 74 6f 72 20 | 41 28 38 30 30 30 29 3b |nVector |A(8000);|
|000059f0| 20 73 32 20 3d 20 41 2e | 53 74 6f 72 65 28 29 3b | s2 = A.|Store();|
|00005a00| 20 7d 0a 58 20 20 20 63 | 6f 75 74 20 3c 3c 20 22 | }.X c|out << "|
|00005a10| 5c 6e 5c 6e 43 68 65 63 | 6b 69 6e 67 20 66 6f 72 |\n\nChec|king for|
|00005a20| 20 6c 6f 73 74 20 6d 65 | 6d 6f 72 79 3a 20 22 0a | lost me|mory: ".|
|00005a30| 58 20 20 20 20 20 20 3c | 3c 20 28 75 6e 73 69 67 |X <|< (unsig|
|00005a40| 6e 65 64 20 6c 6f 6e 67 | 29 73 31 20 3c 3c 20 22 |ned long|)s1 << "|
|00005a50| 20 22 20 3c 3c 20 28 75 | 6e 73 69 67 6e 65 64 20 | " << (u|nsigned |
|00005a60| 6c 6f 6e 67 29 73 32 20 | 3c 3c 20 22 20 22 3b 0a |long)s2 |<< " ";.|
|00005a70| 58 20 20 20 69 66 20 28 | 73 31 20 21 3d 20 73 32 |X if (|s1 != s2|
|00005a80| 29 20 63 6f 75 74 20 3c | 3c 20 22 20 2d 20 65 72 |) cout <|< " - er|
|00005a90| 72 6f 72 5c 6e 22 3b 20 | 65 6c 73 65 20 63 6f 75 |ror\n"; |else cou|
|00005aa0| 74 20 3c 3c 20 22 20 2d | 20 6f 6b 5c 6e 22 3b 0a |t << " -| ok\n";.|
|00005ab0| 58 0a 58 0a 58 7d 0a 58 | 0a 58 72 65 61 6c 20 74 |X.X.X}.X|.Xreal t|
|00005ac0| 33 28 72 65 61 6c 20 72 | 29 20 7b 20 72 65 74 75 |3(real r|) { retu|
|00005ad0| 72 6e 20 69 6e 74 28 72 | 2a 31 30 30 30 29 20 2f |rn int(r|*1000) /|
|00005ae0| 20 31 30 30 30 2e 30 3b | 20 7d 0a 45 4e 44 5f 4f | 1000.0;| }.END_O|
|00005af0| 46 5f 46 49 4c 45 0a 69 | 66 20 74 65 73 74 20 39 |F_FILE.i|f test 9|
|00005b00| 32 39 33 20 2d 6e 65 20 | 60 77 63 20 2d 63 20 3c |293 -ne |`wc -c <|
|00005b10| 27 65 78 61 6d 70 6c 65 | 2e 63 78 78 27 60 3b 20 |'example|.cxx'`; |
|00005b20| 74 68 65 6e 0a 20 20 20 | 20 65 63 68 6f 20 73 68 |then. | echo sh|
|00005b30| 61 72 3a 20 5c 22 27 65 | 78 61 6d 70 6c 65 2e 63 |ar: \"'e|xample.c|
|00005b40| 78 78 27 5c 22 20 75 6e | 70 61 63 6b 65 64 20 77 |xx'\" un|packed w|
|00005b50| 69 74 68 20 77 72 6f 6e | 67 20 73 69 7a 65 21 0a |ith wron|g size!.|
|00005b60| 66 69 0a 23 20 65 6e 64 | 20 6f 66 20 27 65 78 61 |fi.# end| of 'exa|
|00005b70| 6d 70 6c 65 2e 63 78 78 | 27 0a 66 69 0a 69 66 20 |mple.cxx|'.fi.if |
|00005b80| 74 65 73 74 20 2d 66 20 | 27 65 78 61 6d 70 6c 65 |test -f |'example|
|00005b90| 2e 64 65 70 27 20 2d 61 | 20 22 24 7b 31 7d 22 20 |.dep' -a| "${1}" |
|00005ba0| 21 3d 20 22 2d 63 22 20 | 3b 20 74 68 65 6e 20 0a |!= "-c" |; then .|
|00005bb0| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 57 69 6c | echo s|har: Wil|
|00005bc0| 6c 20 6e 6f 74 20 63 6c | 6f 62 62 65 72 20 65 78 |l not cl|obber ex|
|00005bd0| 69 73 74 69 6e 67 20 66 | 69 6c 65 20 5c 22 27 65 |isting f|ile \"'e|
|00005be0| 78 61 6d 70 6c 65 2e 64 | 65 70 27 5c 22 0a 65 6c |xample.d|ep'\".el|
|00005bf0| 73 65 0a 65 63 68 6f 20 | 73 68 61 72 3a 20 45 78 |se.echo |shar: Ex|
|00005c00| 74 72 61 63 74 69 6e 67 | 20 5c 22 27 65 78 61 6d |tracting| \"'exam|
|00005c10| 70 6c 65 2e 64 65 70 27 | 5c 22 20 5c 28 31 37 33 |ple.dep'|\" \(173|
|00005c20| 34 20 63 68 61 72 61 63 | 74 65 72 73 5c 29 0a 73 |4 charac|ters\).s|
|00005c30| 65 64 20 22 73 2f 5e 58 | 2f 2f 22 20 3e 27 65 78 |ed "s/^X|//" >'ex|
|00005c40| 61 6d 70 6c 65 2e 64 65 | 70 27 20 3c 3c 27 45 4e |ample.de|p' <<'EN|
|00005c50| 44 5f 4f 46 5f 46 49 4c | 45 27 0a 58 0a 58 0a 58 |D_OF_FIL|E'.X.X.X|
|00005c60| 4f 42 4a 20 20 3d 20 20 | 65 78 61 6d 70 6c 65 2e |OBJ = |example.|
|00005c70| 6f 20 5c 0a 58 20 20 20 | 20 20 20 20 20 63 68 6f |o \.X | cho|
|00005c80| 6c 65 73 6b 79 2e 6f 20 | 65 76 61 6c 75 65 2e 6f |lesky.o |evalue.o|
|00005c90| 20 66 66 74 2e 6f 20 68 | 68 6f 6c 64 65 72 2e 6f | fft.o h|holder.o|
|00005ca0| 20 6a 61 63 6f 62 69 2e | 6f 20 5c 0a 58 20 20 20 | jacobi.|o \.X |
|00005cb0| 20 20 20 20 20 6e 65 77 | 6d 61 74 31 2e 6f 20 6e | new|mat1.o n|
|00005cc0| 65 77 6d 61 74 32 2e 6f | 20 6e 65 77 6d 61 74 33 |ewmat2.o| newmat3|
|00005cd0| 2e 6f 20 6e 65 77 6d 61 | 74 34 2e 6f 20 6e 65 77 |.o newma|t4.o new|
|00005ce0| 6d 61 74 35 2e 6f 20 5c | 0a 58 20 20 20 20 20 20 |mat5.o \|.X |
|00005cf0| 20 20 6e 65 77 6d 61 74 | 36 2e 6f 20 6e 65 77 6d | newmat|6.o newm|
|00005d00| 61 74 37 2e 6f 20 6e 65 | 77 6d 61 74 38 2e 6f 20 |at7.o ne|wmat8.o |
|00005d10| 6e 65 77 6d 61 74 39 2e | 6f 20 6e 65 77 6d 61 74 |newmat9.|o newmat|
|00005d20| 72 6d 2e 6f 20 5c 0a 58 | 20 20 20 20 20 20 20 20 |rm.o \.X| |
|00005d30| 73 6f 72 74 2e 6f 20 73 | 75 62 6d 61 74 2e 6f 20 |sort.o s|ubmat.o |
|00005d40| 73 76 64 2e 6f 0a 58 0a | 58 0a 58 65 78 61 6d 70 |svd.o.X.|X.Xexamp|
|00005d50| 6c 65 2e 6f 3a 20 20 69 | 6e 63 6c 75 64 65 2e 68 |le.o: i|nclude.h|
|00005d60| 78 78 20 6e 65 77 6d 61 | 74 61 70 2e 68 78 78 20 |xx newma|tap.hxx |
|00005d70| 6e 65 77 6d 61 74 2e 68 | 78 78 20 62 6f 6f 6c 65 |newmat.h|xx boole|
|00005d80| 61 6e 2e 68 78 78 20 65 | 78 61 6d 70 6c 65 2e 63 |an.hxx e|xample.c|
|00005d90| 78 78 0a 58 0a 58 63 68 | 6f 6c 65 73 6b 79 2e 6f |xx.X.Xch|olesky.o|
|00005da0| 3a 20 69 6e 63 6c 75 64 | 65 2e 68 78 78 20 6e 65 |: includ|e.hxx ne|
|00005db0| 77 6d 61 74 2e 68 78 78 | 20 62 6f 6f 6c 65 61 6e |wmat.hxx| boolean|
|00005dc0| 2e 68 78 78 20 63 68 6f | 6c 65 73 6b 79 2e 63 78 |.hxx cho|lesky.cx|
|00005dd0| 78 0a 58 0a 58 65 76 61 | 6c 75 65 2e 6f 3a 20 20 |x.X.Xeva|lue.o: |
|00005de0| 20 69 6e 63 6c 75 64 65 | 2e 68 78 78 20 6e 65 77 | include|.hxx new|
|00005df0| 6d 61 74 2e 68 78 78 20 | 62 6f 6f 6c 65 61 6e 2e |mat.hxx |boolean.|
|00005e00| 68 78 78 20 6e 65 77 6d | 61 74 72 6d 2e 68 78 78 |hxx newm|atrm.hxx|
|00005e10| 20 70 72 65 63 69 73 69 | 6f 2e 68 78 78 20 65 76 | precisi|o.hxx ev|
|00005e20| 61 6c 75 65 2e 63 78 78 | 0a 58 0a 58 66 66 74 2e |alue.cxx|.X.Xfft.|
|00005e30| 6f 3a 20 20 20 20 20 20 | 69 6e 63 6c 75 64 65 2e |o: |include.|
|00005e40| 68 78 78 20 6e 65 77 6d | 61 74 61 70 2e 68 78 78 |hxx newm|atap.hxx|
|00005e50| 20 6e 65 77 6d 61 74 2e | 68 78 78 20 62 6f 6f 6c | newmat.|hxx bool|
|00005e60| 65 61 6e 2e 68 78 78 20 | 66 66 74 2e 63 78 78 0a |ean.hxx |fft.cxx.|
|00005e70| 58 0a 58 68 68 6f 6c 64 | 65 72 2e 6f 3a 20 20 69 |X.Xhhold|er.o: i|
|00005e80| 6e 63 6c 75 64 65 2e 68 | 78 78 20 6e 65 77 6d 61 |nclude.h|xx newma|
|00005e90| 74 61 70 2e 68 78 78 20 | 6e 65 77 6d 61 74 2e 68 |tap.hxx |newmat.h|
|00005ea0| 78 78 20 62 6f 6f 6c 65 | 61 6e 2e 68 78 78 20 68 |xx boole|an.hxx h|
|00005eb0| 68 6f 6c 64 65 72 2e 63 | 78 78 0a 58 0a 58 6a 61 |holder.c|xx.X.Xja|
|00005ec0| 63 6f 62 69 2e 6f 3a 20 | 20 20 69 6e 63 6c 75 64 |cobi.o: | includ|
|00005ed0| 65 2e 68 78 78 20 6e 65 | 77 6d 61 74 2e 68 78 78 |e.hxx ne|wmat.hxx|
|00005ee0| 20 62 6f 6f 6c 65 61 6e | 2e 68 78 78 20 70 72 65 | boolean|.hxx pre|
|00005ef0| 63 69 73 69 6f 2e 68 78 | 78 20 6e 65 77 6d 61 74 |cisio.hx|x newmat|
|00005f00| 72 6d 2e 68 78 78 20 6a | 61 63 6f 62 69 2e 63 78 |rm.hxx j|acobi.cx|
|00005f10| 78 0a 58 0a 58 6e 65 77 | 6d 61 74 31 2e 6f 3a 20 |x.X.Xnew|mat1.o: |
|00005f20| 20 69 6e 63 6c 75 64 65 | 2e 68 78 78 20 6e 65 77 | include|.hxx new|
|00005f30| 6d 61 74 2e 68 78 78 20 | 62 6f 6f 6c 65 61 6e 2e |mat.hxx |boolean.|
|00005f40| 68 78 78 20 6e 65 77 6d | 61 74 31 2e 63 78 78 0a |hxx newm|at1.cxx.|
|00005f50| 58 0a 58 6e 65 77 6d 61 | 74 32 2e 6f 3a 20 20 69 |X.Xnewma|t2.o: i|
|00005f60| 6e 63 6c 75 64 65 2e 68 | 78 78 20 6e 65 77 6d 61 |nclude.h|xx newma|
|00005f70| 74 2e 68 78 78 20 62 6f | 6f 6c 65 61 6e 2e 68 78 |t.hxx bo|olean.hx|
|00005f80| 78 20 6e 65 77 6d 61 74 | 72 63 2e 68 78 78 20 63 |x newmat|rc.hxx c|
|00005f90| 6f 6e 74 72 6f 6c 77 2e | 68 78 78 20 6e 65 77 6d |ontrolw.|hxx newm|
|00005fa0| 61 74 32 2e 63 78 78 0a | 58 0a 58 6e 65 77 6d 61 |at2.cxx.|X.Xnewma|
|00005fb0| 74 33 2e 6f 3a 20 20 69 | 6e 63 6c 75 64 65 2e 68 |t3.o: i|nclude.h|
|00005fc0| 78 78 20 6e 65 77 6d 61 | 74 2e 68 78 78 20 62 6f |xx newma|t.hxx bo|
|00005fd0| 6f 6c 65 61 6e 2e 68 78 | 78 20 6e 65 77 6d 61 74 |olean.hx|x newmat|
|00005fe0| 72 63 2e 68 78 78 20 63 | 6f 6e 74 72 6f 6c 77 2e |rc.hxx c|ontrolw.|
|00005ff0| 68 78 78 20 6e 65 77 6d | 61 74 33 2e 63 78 78 0a |hxx newm|at3.cxx.|
|00006000| 58 0a 58 6e 65 77 6d 61 | 74 34 2e 6f 3a 20 20 69 |X.Xnewma|t4.o: i|
|00006010| 6e 63 6c 75 64 65 2e 68 | 78 78 20 6e 65 77 6d 61 |nclude.h|xx newma|
|00006020| 74 2e 68 78 78 20 62 6f | 6f 6c 65 61 6e 2e 68 78 |t.hxx bo|olean.hx|
|00006030| 78 20 6e 65 77 6d 61 74 | 72 63 2e 68 78 78 20 63 |x newmat|rc.hxx c|
|00006040| 6f 6e 74 72 6f 6c 77 2e | 68 78 78 20 6e 65 77 6d |ontrolw.|hxx newm|
|00006050| 61 74 34 2e 63 78 78 0a | 58 0a 58 6e 65 77 6d 61 |at4.cxx.|X.Xnewma|
|00006060| 74 35 2e 6f 3a 20 20 69 | 6e 63 6c 75 64 65 2e 68 |t5.o: i|nclude.h|
|00006070| 78 78 20 6e 65 77 6d 61 | 74 2e 68 78 78 20 62 6f |xx newma|t.hxx bo|
|00006080| 6f 6c 65 61 6e 2e 68 78 | 78 20 6e 65 77 6d 61 74 |olean.hx|x newmat|
|00006090| 72 63 2e 68 78 78 20 63 | 6f 6e 74 72 6f 6c 77 2e |rc.hxx c|ontrolw.|
|000060a0| 68 78 78 20 6e 65 77 6d | 61 74 35 2e 63 78 78 0a |hxx newm|at5.cxx.|
|000060b0| 58 0a 58 6e 65 77 6d 61 | 74 36 2e 6f 3a 20 20 69 |X.Xnewma|t6.o: i|
|000060c0| 6e 63 6c 75 64 65 2e 68 | 78 78 20 6e 65 77 6d 61 |nclude.h|xx newma|
|000060d0| 74 2e 68 78 78 20 62 6f | 6f 6c 65 61 6e 2e 68 78 |t.hxx bo|olean.hx|
|000060e0| 78 20 6e 65 77 6d 61 74 | 72 63 2e 68 78 78 20 63 |x newmat|rc.hxx c|
|000060f0| 6f 6e 74 72 6f 6c 77 2e | 68 78 78 20 6e 65 77 6d |ontrolw.|hxx newm|
|00006100| 61 74 36 2e 63 78 78 0a | 58 0a 58 6e 65 77 6d 61 |at6.cxx.|X.Xnewma|
|00006110| 74 37 2e 6f 3a 20 20 69 | 6e 63 6c 75 64 65 2e 68 |t7.o: i|nclude.h|
|00006120| 78 78 20 6e 65 77 6d 61 | 74 2e 68 78 78 20 62 6f |xx newma|t.hxx bo|
|00006130| 6f 6c 65 61 6e 2e 68 78 | 78 20 6e 65 77 6d 61 74 |olean.hx|x newmat|
|00006140| 72 63 2e 68 78 78 20 63 | 6f 6e 74 72 6f 6c 77 2e |rc.hxx c|ontrolw.|
|00006150| 68 78 78 20 6e 65 77 6d | 61 74 37 2e 63 78 78 0a |hxx newm|at7.cxx.|
|00006160| 58 0a 58 6e 65 77 6d 61 | 74 38 2e 6f 3a 20 20 69 |X.Xnewma|t8.o: i|
|00006170| 6e 63 6c 75 64 65 2e 68 | 78 78 20 6e 65 77 6d 61 |nclude.h|xx newma|
|00006180| 74 61 70 2e 68 78 78 20 | 6e 65 77 6d 61 74 2e 68 |tap.hxx |newmat.h|
|00006190| 78 78 20 62 6f 6f 6c 65 | 61 6e 2e 68 78 78 20 6e |xx boole|an.hxx n|
|000061a0| 65 77 6d 61 74 38 2e 63 | 78 78 0a 58 0a 58 6e 65 |ewmat8.c|xx.X.Xne|
|000061b0| 77 6d 61 74 39 2e 6f 3a | 20 20 69 6e 63 6c 75 64 |wmat9.o:| includ|
|000061c0| 65 2e 68 78 78 20 6e 65 | 77 6d 61 74 2e 68 78 78 |e.hxx ne|wmat.hxx|
|000061d0| 20 62 6f 6f 6c 65 61 6e | 2e 68 78 78 20 6e 65 77 | boolean|.hxx new|
|000061e0| 6d 61 74 72 63 2e 68 78 | 78 20 63 6f 6e 74 72 6f |matrc.hx|x contro|
|000061f0| 6c 77 2e 68 78 78 20 6e | 65 77 6d 61 74 69 6f 2e |lw.hxx n|ewmatio.|
|00006200| 68 78 78 20 6e 65 77 6d | 61 74 39 2e 63 78 78 0a |hxx newm|at9.cxx.|
|00006210| 58 0a 58 6e 65 77 6d 61 | 74 72 6d 2e 6f 3a 20 69 |X.Xnewma|trm.o: i|
|00006220| 6e 63 6c 75 64 65 2e 68 | 78 78 20 6e 65 77 6d 61 |nclude.h|xx newma|
|00006230| 74 2e 68 78 78 20 62 6f | 6f 6c 65 61 6e 2e 68 78 |t.hxx bo|olean.hx|
|00006240| 78 20 6e 65 77 6d 61 74 | 72 6d 2e 68 78 78 20 6e |x newmat|rm.hxx n|
|00006250| 65 77 6d 61 74 72 6d 2e | 63 78 78 0a 58 0a 58 73 |ewmatrm.|cxx.X.Xs|
|00006260| 6f 72 74 2e 6f 3a 20 20 | 20 20 20 69 6e 63 6c 75 |ort.o: | inclu|
|00006270| 64 65 2e 68 78 78 20 6e | 65 77 6d 61 74 61 70 2e |de.hxx n|ewmatap.|
|00006280| 68 78 78 20 6e 65 77 6d | 61 74 2e 68 78 78 20 62 |hxx newm|at.hxx b|
|00006290| 6f 6f 6c 65 61 6e 2e 68 | 78 78 20 73 6f 72 74 2e |oolean.h|xx sort.|
|000062a0| 63 78 78 0a 58 0a 58 73 | 75 62 6d 61 74 2e 6f 3a |cxx.X.Xs|ubmat.o:|
|000062b0| 20 20 20 69 6e 63 6c 75 | 64 65 2e 68 78 78 20 6e | inclu|de.hxx n|
|000062c0| 65 77 6d 61 74 2e 68 78 | 78 20 62 6f 6f 6c 65 61 |ewmat.hx|x boolea|
|000062d0| 6e 2e 68 78 78 20 6e 65 | 77 6d 61 74 72 63 2e 68 |n.hxx ne|wmatrc.h|
|000062e0| 78 78 20 63 6f 6e 74 72 | 6f 6c 77 2e 68 78 78 20 |xx contr|olw.hxx |
|000062f0| 73 75 62 6d 61 74 2e 63 | 78 78 0a 58 0a 58 73 76 |submat.c|xx.X.Xsv|
|00006300| 64 2e 6f 3a 20 20 20 20 | 20 20 69 6e 63 6c 75 64 |d.o: | includ|
|00006310| 65 2e 68 78 78 20 6e 65 | 77 6d 61 74 2e 68 78 78 |e.hxx ne|wmat.hxx|
|00006320| 20 62 6f 6f 6c 65 61 6e | 2e 68 78 78 20 6e 65 77 | boolean|.hxx new|
|00006330| 6d 61 74 72 6d 2e 68 78 | 78 20 70 72 65 63 69 73 |matrm.hx|x precis|
|00006340| 69 6f 2e 68 78 78 20 73 | 76 64 2e 63 78 78 0a 45 |io.hxx s|vd.cxx.E|
|00006350| 4e 44 5f 4f 46 5f 46 49 | 4c 45 0a 69 66 20 74 65 |ND_OF_FI|LE.if te|
|00006360| 73 74 20 31 37 33 34 20 | 2d 6e 65 20 60 77 63 20 |st 1734 |-ne `wc |
|00006370| 2d 63 20 3c 27 65 78 61 | 6d 70 6c 65 2e 64 65 70 |-c <'exa|mple.dep|
|00006380| 27 60 3b 20 74 68 65 6e | 0a 20 20 20 20 65 63 68 |'`; then|. ech|
|00006390| 6f 20 73 68 61 72 3a 20 | 5c 22 27 65 78 61 6d 70 |o shar: |\"'examp|
|000063a0| 6c 65 2e 64 65 70 27 5c | 22 20 75 6e 70 61 63 6b |le.dep'\|" unpack|
|000063b0| 65 64 20 77 69 74 68 20 | 77 72 6f 6e 67 20 73 69 |ed with |wrong si|
|000063c0| 7a 65 21 0a 66 69 0a 23 | 20 65 6e 64 20 6f 66 20 |ze!.fi.#| end of |
|000063d0| 27 65 78 61 6d 70 6c 65 | 2e 64 65 70 27 0a 66 69 |'example|.dep'.fi|
|000063e0| 0a 69 66 20 74 65 73 74 | 20 2d 66 20 27 65 78 61 |.if test| -f 'exa|
|000063f0| 6d 70 6c 65 2e 74 78 74 | 27 20 2d 61 20 22 24 7b |mple.txt|' -a "${|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
