home *** CD-ROM | disk | FTP | other *** search
open in:
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: LaTeX Document
(document/latex).
Confidence | Program | Detection | Match Type | Support
|
---|
100%
| dexvert
| LaTeX Document (document/latex)
| magic
| Supported |
1%
| dexvert
| Corel 10 Texture (image/corel10Texture)
| ext
| Unsupported |
1%
| dexvert
| Croteam texture file (image/croteamTextureFile)
| ext
| Unsupported |
1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
100%
| file
| LaTeX document text
| default
| |
99%
| file
| C source text
| default
| |
98%
| file
| C source, ASCII text, with CRLF line terminators
| default
| |
100%
| checkBytes
| Printable ASCII
| default
| |
100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
100%
| siegfried
| fmt/281 LaTeX (Subdocument)
| default
| |
100%
| detectItEasy
| Format: plain text[CRLF]
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 5c 61 70 70 65 6e 64 69 | 78 0d 0a 5c 63 68 61 70 |\appendi|x..\chap|
|00000010| 74 65 72 7b 74 65 73 74 | 2e 63 70 70 7d 0d 0a 54 |ter{test|.cpp}..T|
|00000020| 68 69 73 20 69 73 20 74 | 68 65 20 74 65 73 74 20 |his is t|he test |
|00000030| 70 72 6f 67 72 61 6d 20 | 66 6f 72 20 74 68 65 20 |program |for the |
|00000040| 56 4d 61 74 72 69 78 20 | 73 79 73 74 65 6d 2e 20 |VMatrix |system. |
|00000050| 49 74 0d 0a 72 65 71 75 | 69 72 65 73 20 74 68 65 |It..requ|ires the|
|00000060| 20 64 61 74 61 20 73 65 | 74 20 5c 76 65 72 62 2b | data se|t \verb+|
|00000070| 20 63 61 74 63 68 76 2e | 64 61 74 2b 2e 0d 0a 5c | catchv.|dat+...\|
|00000080| 62 65 67 69 6e 7b 76 65 | 72 62 61 74 69 6d 7d 0d |begin{ve|rbatim}.|
|00000090| 0a 0d 0a 2f 2f 20 6d 61 | 6b 65 20 73 75 72 65 20 |...// ma|ke sure |
|000000a0| 74 6f 20 64 65 66 69 6e | 65 20 74 68 65 20 67 6c |to defin|e the gl|
|000000b0| 6f 62 61 6c 20 64 65 66 | 69 6e 65 20 76 61 72 69 |obal def|ine vari|
|000000c0| 61 62 6c 65 20 20 49 4e | 5f 52 41 4d 20 69 66 20 |able IN|_RAM if |
|000000d0| 79 6f 75 20 77 61 6e 74 | 0d 0a 2f 2f 20 74 6f 20 |you want|..// to |
|000000e0| 75 73 65 20 74 68 65 20 | 6d 65 64 69 75 6d 20 6d |use the |medium m|
|000000f0| 6f 64 65 6c 2e 20 54 68 | 69 73 20 63 61 6e 20 62 |odel. Th|is can b|
|00000100| 65 20 64 6f 6e 65 20 69 | 6e 20 74 68 65 20 63 6f |e done i|n the co|
|00000110| 6d 70 69 6c 65 72 20 63 | 6f 64 65 0d 0a 2f 2f 20 |mpiler c|ode..// |
|00000120| 67 65 6e 65 72 61 74 69 | 6f 6e 20 77 69 6e 64 6f |generati|on windo|
|00000130| 77 20 6f 66 20 63 6f 6d | 70 69 6c 65 72 20 6f 70 |w of com|piler op|
|00000140| 74 69 6f 6e 73 2e 20 49 | 74 20 6d 61 79 20 61 6c |tions. I|t may al|
|00000150| 73 6f 20 62 65 20 64 65 | 66 69 6e 65 64 20 66 6f |so be de|fined fo|
|00000160| 72 0d 0a 2f 2f 20 74 68 | 65 20 63 6f 6d 6d 61 6e |r..// th|e comman|
|00000170| 64 20 6c 69 6e 65 20 63 | 6f 6d 70 69 6c 65 72 0d |d line c|ompiler.|
|00000180| 0a 0d 0a 23 69 6e 63 6c | 75 64 65 20 22 76 69 72 |...#incl|ude "vir|
|00000190| 74 2e 68 22 0d 0a 0d 0a | 2f 2a 20 72 65 71 75 69 |t.h"....|/* requi|
|000001a0| 72 65 64 20 67 6c 6f 62 | 61 6c 20 64 65 63 6c 61 |red glob|al decla|
|000001b0| 72 61 74 69 6f 6e 20 66 | 6f 72 20 74 68 65 20 6d |ration f|or the m|
|000001c0| 61 74 72 69 78 20 73 74 | 61 63 6b 20 6f 62 6a 65 |atrix st|ack obje|
|000001d0| 63 74 20 2a 2f 0d 0a 0d | 0a 2f 2f 75 6e 73 69 67 |ct */...|.//unsig|
|000001e0| 6e 65 64 20 69 6e 74 20 | 5f 73 74 6b 6c 65 6e 20 |ned int |_stklen |
|000001f0| 3d 20 53 54 41 43 4b 4c | 45 4e 47 54 48 3b 0d 0a |= STACKL|ENGTH;..|
|00000200| 0d 0a 4d 53 74 61 63 6b | 20 2a 44 69 73 70 61 74 |..MStack| *Dispat|
|00000210| 63 68 20 3d 20 6e 65 77 | 20 4d 53 74 61 63 6b 3b |ch = new| MStack;|
|00000220| 0d 0a 0d 0a 0d 0a 56 4d | 61 74 72 69 78 26 20 66 |......VM|atrix& f|
|00000230| 75 6e 63 74 69 6f 6e 31 | 28 20 56 4d 61 74 72 69 |unction1|( VMatri|
|00000240| 78 20 26 61 2c 20 56 4d | 61 74 72 69 78 20 26 62 |x &a, VM|atrix &b|
|00000250| 20 29 0d 0a 7b 0d 0a 20 | 20 20 2f 2f 20 54 68 69 | )..{.. | // Thi|
|00000260| 73 20 66 75 6e 63 74 69 | 6f 6e 20 74 65 73 74 73 |s functi|on tests|
|00000270| 20 74 68 65 20 66 72 65 | 69 6e 64 20 66 75 6e 63 | the fre|ind func|
|00000280| 74 69 6f 6e 73 20 61 6e | 64 20 72 65 74 75 72 6e |tions an|d return|
|00000290| 73 20 61 20 76 61 6c 75 | 65 0d 0a 0d 0a 20 20 20 |s a valu|e.... |
|000002a0| 61 2e 47 61 72 62 61 67 | 65 28 22 66 75 6e 63 74 |a.Garbag|e("funct|
|000002b0| 69 6f 6e 31 22 29 3b 20 | 20 20 20 20 20 20 20 2f |ion1"); | /|
|000002c0| 2f 20 63 68 65 63 6b 20 | 61 20 61 6e 64 20 62 0d |/ check |a and b.|
|000002d0| 0a 20 20 20 62 2e 47 61 | 72 62 61 67 65 28 22 66 |. b.Ga|rbage("f|
|000002e0| 75 6e 63 74 69 6f 6e 31 | 22 29 3b 0d 0a 0d 0a 20 |unction1|");.... |
|000002f0| 20 20 44 69 73 70 61 74 | 63 68 2d 3e 49 6e 63 6c | Dispat|ch->Incl|
|00000300| 65 76 65 6c 28 29 3b 20 | 20 20 20 20 20 20 20 20 |evel(); | |
|00000310| 20 2f 2f 20 69 6e 63 72 | 65 6d 65 6e 74 20 70 75 | // incr|ement pu|
|00000320| 73 68 2d 70 6f 70 20 6c | 65 76 65 6c 0d 0a 20 20 |sh-pop l|evel.. |
|00000330| 20 56 4d 61 74 72 69 78 | 20 63 28 22 63 22 2c 31 | VMatrix| c("c",1|
|00000340| 2c 31 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |,1); | |
|00000350| 2f 2f 20 63 72 65 61 74 | 65 20 61 20 6c 6f 63 61 |// creat|e a loca|
|00000360| 6c 20 6d 61 74 72 69 78 | 0d 0a 0d 0a 20 20 20 63 |l matrix|.... c|
|00000370| 20 3d 20 61 20 2b 20 62 | 20 2b 20 61 20 2b 20 62 | = a + b| + a + b|
|00000380| 3b 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2f 2f |; | //|
|00000390| 20 63 68 65 63 6b 20 61 | 20 72 65 70 65 61 74 65 | check a| repeate|
|000003a0| 64 20 6d 61 74 72 69 78 | 20 61 64 64 69 74 69 6f |d matrix| additio|
|000003b0| 6e 0d 0a 20 20 20 63 2e | 44 69 73 70 6c 61 79 4d |n.. c.|DisplayM|
|000003c0| 61 74 28 29 3b 20 20 20 | 20 20 20 20 20 20 20 20 |at(); | |
|000003d0| 20 20 20 20 20 2f 2f 20 | 70 72 69 6e 74 20 63 0d | // |print c.|
|000003e0| 0a 20 20 20 63 20 3d 20 | 34 33 31 2e 32 20 2b 20 |. c = |431.2 + |
|000003f0| 54 72 61 6e 28 61 29 2a | 62 20 2b 20 32 2e 31 33 |Tran(a)*|b + 2.13|
|00000400| 34 3b 20 2f 2f 20 63 68 | 65 63 6b 20 63 6f 6d 6d |4; // ch|eck comm|
|00000410| 75 74 69 76 69 74 79 20 | 6f 66 20 73 63 61 6c 61 |utivity |of scala|
|00000420| 72 20 61 64 64 69 74 69 | 6f 6e 0d 0a 20 20 20 63 |r additi|on.. c|
|00000430| 2e 44 69 73 70 6c 61 79 | 4d 61 74 28 29 3b 20 20 |.Display|Mat(); |
|00000440| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2f 2f | | //|
|00000450| 20 70 72 69 6e 74 20 63 | 0d 0a 20 20 20 63 20 3d | print c|.. c =|
|00000460| 20 2d 63 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 | -c; | |
|00000470| 20 20 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 63 | | // c|
|00000480| 68 65 63 6b 20 75 6e 69 | 61 72 79 20 6d 69 6e 75 |heck uni|ary minu|
|00000490| 73 0d 0a 20 20 20 63 2e | 44 69 73 70 6c 61 79 4d |s.. c.|DisplayM|
|000004a0| 61 74 28 29 3b 20 20 20 | 20 20 20 20 20 20 20 20 |at(); | |
|000004b0| 20 20 20 20 20 2f 2f 20 | 70 72 69 6e 74 20 63 0d | // |print c.|
|000004c0| 0a 20 20 20 63 20 3d 20 | 61 20 2d 20 62 3b 20 20 |. c = |a - b; |
|000004d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000004e0| 20 20 20 2f 2f 20 6d 61 | 74 72 69 78 20 73 75 62 | // ma|trix sub|
|000004f0| 72 61 63 74 69 6f 6e 0d | 0a 20 20 20 63 2e 44 69 |raction.|. c.Di|
|00000500| 73 70 6c 61 79 4d 61 74 | 28 29 3b 20 20 20 20 20 |splayMat|(); |
|00000510| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 70 72 | | // pr|
|00000520| 69 6e 74 20 63 0d 0a 20 | 20 20 63 20 3d 20 35 20 |int c.. | c = 5 |
|00000530| 2d 20 61 20 2d 20 35 3b | 20 20 20 20 20 20 20 20 |- a - 5;| |
|00000540| 20 20 20 20 20 20 20 20 | 20 2f 2f 20 63 68 65 63 | | // chec|
|00000550| 6b 20 63 6f 6d 6d 75 74 | 69 76 69 74 79 20 6f 66 |k commut|ivity of|
|00000560| 20 73 63 61 6c 61 72 20 | 73 75 62 74 72 61 63 74 | scalar |subtract|
|00000570| 69 6f 6e 0d 0a 20 20 20 | 63 2e 44 69 73 70 6c 61 |ion.. |c.Displa|
|00000580| 79 4d 61 74 28 29 3b 20 | 20 20 20 20 20 20 20 20 |yMat(); | |
|00000590| 20 20 20 20 20 20 20 2f | 2f 20 70 72 69 6e 74 20 | /|/ print |
|000005a0| 63 0d 0a 20 20 20 63 20 | 3d 20 35 2a 61 2a 35 3b |c.. c |= 5*a*5;|
|000005b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000005c0| 20 20 20 20 20 2f 2f 20 | 63 68 65 63 6b 20 63 6f | // |check co|
|000005d0| 6d 6d 75 74 69 76 69 74 | 79 20 6f 66 20 73 63 61 |mmutivit|y of sca|
|000005e0| 6c 61 72 20 6d 75 6c 74 | 69 70 6c 69 63 61 74 69 |lar mult|iplicati|
|000005f0| 6f 6e 0d 0a 20 20 20 63 | 2e 44 69 73 70 6c 61 79 |on.. c|.Display|
|00000600| 4d 61 74 28 29 3b 20 20 | 20 20 20 20 20 20 20 20 |Mat(); | |
|00000610| 20 20 20 20 20 20 2f 2f | 20 70 72 69 6e 74 20 63 | //| print c|
|00000620| 0d 0a 20 20 20 63 20 3d | 20 61 25 61 3b 20 20 20 |.. c =| a%a; |
|00000630| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000640| 20 20 20 20 2f 2f 20 63 | 68 65 63 6b 20 65 6c 65 | // c|heck ele|
|00000650| 6d 65 6e 74 77 69 73 65 | 20 6d 75 6c 74 69 70 6c |mentwise| multipl|
|00000660| 69 63 61 74 69 6f 6e 0d | 0a 20 20 20 63 2e 44 69 |ication.|. c.Di|
|00000670| 73 70 6c 61 79 4d 61 74 | 28 29 3b 20 20 20 20 20 |splayMat|(); |
|00000680| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 70 72 | | // pr|
|00000690| 69 6e 74 20 63 0d 0a 20 | 20 20 63 20 3d 20 61 2f |int c.. | c = a/|
|000006a0| 31 32 33 34 3b 20 20 20 | 20 20 20 20 20 20 20 20 |1234; | |
|000006b0| 20 20 20 20 20 20 20 20 | 20 2f 2f 20 63 68 65 63 | | // chec|
|000006c0| 6b 20 73 63 61 6c 61 72 | 20 64 69 76 69 73 69 6f |k scalar| divisio|
|000006d0| 6e 0d 0a 20 20 20 63 2e | 44 69 73 70 6c 61 79 4d |n.. c.|DisplayM|
|000006e0| 61 74 28 29 3b 20 20 20 | 20 20 20 20 20 20 20 20 |at(); | |
|000006f0| 20 20 20 20 20 2f 2f 20 | 70 72 69 6e 74 20 63 0d | // |print c.|
|00000700| 0a 20 20 20 63 20 3d 20 | 61 2f 62 3b 20 20 20 20 |. c = |a/b; |
|00000710| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000720| 20 20 20 2f 2f 20 63 68 | 65 63 6b 20 65 6c 65 6d | // ch|eck elem|
|00000730| 65 6e 74 77 69 73 65 20 | 64 69 76 69 73 69 6f 6e |entwise |division|
|00000740| 0d 0a 20 20 20 63 2e 44 | 69 73 70 6c 61 79 4d 61 |.. c.D|isplayMa|
|00000750| 74 28 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |t(); | |
|00000760| 20 20 20 20 2f 2f 20 70 | 72 69 6e 74 20 63 0d 0a | // p|rint c..|
|00000770| 0d 0a 20 20 20 44 69 73 | 70 61 74 63 68 2d 3e 50 |.. Dis|patch->P|
|00000780| 72 69 6e 74 53 74 61 63 | 6b 28 29 3b 20 20 20 20 |rintStac|k(); |
|00000790| 20 20 20 20 2f 2f 20 73 | 68 6f 77 20 73 74 61 63 | // s|how stac|
|000007a0| 6b 20 62 65 66 6f 72 65 | 20 70 75 73 68 0d 0a 20 |k before| push.. |
|000007b0| 20 20 44 69 73 70 61 74 | 63 68 2d 3e 50 75 73 68 | Dispat|ch->Push|
|000007c0| 28 63 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |(c); | |
|000007d0| 20 2f 2f 20 70 75 73 68 | 20 63 20 6f 6e 74 6f 20 | // push| c onto |
|000007e0| 73 74 61 63 6b 0d 0a 20 | 20 20 44 69 73 70 61 74 |stack.. | Dispat|
|000007f0| 63 68 2d 3e 50 72 69 6e | 74 53 74 61 63 6b 28 29 |ch->Prin|tStack()|
|00000800| 3b 20 20 20 20 20 20 20 | 20 2f 2f 20 65 78 61 6d |; | // exam|
|00000810| 69 6e 65 20 73 74 61 63 | 6b 20 61 66 74 65 72 20 |ine stac|k after |
|00000820| 70 75 73 68 0d 0a 20 20 | 20 72 65 74 75 72 6e 20 |push.. | return |
|00000830| 44 69 73 70 61 74 63 68 | 2d 3e 44 65 63 52 65 74 |Dispatch|->DecRet|
|00000840| 75 72 6e 28 29 3b 20 20 | 2f 2f 20 64 65 63 72 65 |urn(); |// decre|
|00000850| 6d 65 6e 74 20 70 75 73 | 68 2d 70 6f 70 20 6c 65 |ment pus|h-pop le|
|00000860| 76 65 6c 2c 20 61 6e 64 | 20 72 65 74 75 72 6e 0d |vel, and| return.|
|00000870| 0a 09 09 09 09 20 20 2f | 2f 20 73 74 61 63 6b 20 |..... /|/ stack |
|00000880| 74 6f 70 0d 0a 7d 0d 0a | 0d 0a 56 4d 61 74 72 69 |top..}..|..VMatri|
|00000890| 78 20 26 66 75 6e 63 74 | 69 6f 6e 30 20 28 20 76 |x &funct|ion0 ( v|
|000008a0| 6f 69 64 20 29 0d 0a 7b | 0d 0a 20 20 20 20 2f 2f |oid )..{|.. //|
|000008b0| 20 74 65 73 74 20 73 6f | 6d 65 20 6f 66 20 74 68 | test so|me of th|
|000008c0| 65 20 6f 75 74 70 75 74 | 20 66 75 6e 63 74 69 6f |e output| functio|
|000008d0| 6e 73 20 61 6e 64 20 72 | 61 77 20 6d 61 74 72 69 |ns and r|aw matri|
|000008e0| 78 20 66 75 6e 63 74 69 | 6f 6e 73 0d 0a 0d 0a 20 |x functi|ons.... |
|000008f0| 20 20 20 44 69 73 70 61 | 74 63 68 2d 3e 49 6e 63 | Dispa|tch->Inc|
|00000900| 6c 65 76 65 6c 28 29 3b | 0d 0a 20 20 20 20 56 4d |level();|.. VM|
|00000910| 61 74 72 69 78 20 64 20 | 3d 20 56 4d 61 74 72 69 |atrix d |= VMatri|
|00000920| 78 28 22 64 22 2c 31 2c | 31 29 3b 0d 0a 20 20 20 |x("d",1,|1);.. |
|00000930| 20 56 4d 61 74 72 69 78 | 20 61 2c 63 2c 48 3b 0d | VMatrix| a,c,H;.|
|00000940| 0a 0d 0a 20 20 20 20 61 | 20 3d 20 52 65 61 64 61 |... a| = Reada|
|00000950| 28 20 22 63 61 74 63 68 | 76 2e 64 61 74 22 20 29 |( "catch|v.dat" )|
|00000960| 3b 20 20 20 20 20 20 2f | 2f 20 72 65 61 64 20 61 |; /|/ read a|
|00000970| 6e 20 61 73 63 69 69 20 | 6d 61 74 72 69 78 0d 0a |n ascii |matrix..|
|00000980| 20 20 20 20 61 2e 44 69 | 73 70 6c 61 79 4d 61 74 | a.Di|splayMat|
|00000990| 28 29 3b 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |(); | |
|000009a0| 20 20 20 20 2f 2f 20 64 | 69 73 70 6c 61 79 20 61 | // d|isplay a|
|000009b0| 0d 0a 20 20 20 20 57 72 | 69 74 65 61 28 20 22 6a |.. Wr|itea( "j|
|000009c0| 75 6e 6b 2e 64 61 74 22 | 2c 20 61 29 3b 20 20 20 |unk.dat"|, a); |
|000009d0| 20 20 20 20 20 20 2f 2f | 20 77 72 69 74 65 20 61 | //| write a|
|000009e0| 73 63 69 69 20 6d 61 74 | 72 69 78 0d 0a 20 20 20 |scii mat|rix.. |
|000009f0| 20 61 2e 57 72 69 74 65 | 62 28 20 22 6a 75 6e 6b | a.Write|b( "junk|
|00000a00| 2e 62 69 6e 22 2c 20 61 | 29 3b 20 20 20 20 20 20 |.bin", a|); |
|00000a10| 20 2f 2f 20 77 72 69 74 | 65 20 62 69 6e 61 72 79 | // writ|e binary|
|00000a20| 20 6d 61 74 72 69 78 0d | 0a 20 20 20 20 61 20 3d | matrix.|. a =|
|00000a30| 20 52 65 61 64 62 28 20 | 22 6a 75 6e 6b 2e 62 69 | Readb( |"junk.bi|
|00000a40| 6e 22 20 29 3b 20 20 20 | 20 20 20 20 20 2f 2f 20 |n" ); | // |
|00000a50| 72 65 61 64 20 62 69 6e | 61 72 79 20 6d 61 74 72 |read bin|ary matr|
|00000a60| 69 78 0d 0a 20 20 20 20 | 61 2e 49 6e 66 6f 4d 61 |ix.. |a.InfoMa|
|00000a70| 74 28 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |t(); | |
|00000a80| 20 20 20 20 20 20 20 20 | 2f 2f 20 64 69 73 70 6c | |// displ|
|00000a90| 61 79 20 6d 61 74 72 69 | 78 20 69 6e 66 6f 0d 0a |ay matri|x info..|
|00000aa0| 20 20 20 20 61 2e 44 69 | 73 70 6c 61 79 4d 61 74 | a.Di|splayMat|
|00000ab0| 28 29 3b 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |(); | |
|00000ac0| 20 20 20 20 2f 2f 20 64 | 69 73 70 6c 61 79 20 61 | // d|isplay a|
|00000ad0| 0d 0a 20 20 20 20 64 20 | 3d 20 53 75 62 6d 61 74 |.. d |= Submat|
|00000ae0| 28 20 61 2c 20 61 2e 72 | 2c 20 34 2c 20 31 2c 20 |( a, a.r|, 4, 1, |
|00000af0| 32 20 29 3b 20 20 2f 2f | 20 74 61 6b 65 20 61 20 |2 ); //| take a |
|00000b00| 73 75 62 6d 61 74 72 69 | 78 20 6f 66 20 61 0d 0a |submatri|x of a..|
|00000b10| 20 20 20 20 64 2e 44 69 | 73 70 6c 61 79 4d 61 74 | d.Di|splayMat|
|00000b20| 28 29 3b 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |(); | |
|00000b30| 20 20 20 20 2f 2f 20 64 | 69 73 70 6c 61 79 20 69 | // d|isplay i|
|00000b40| 74 0d 0a 20 20 20 20 63 | 20 3d 20 43 68 28 20 64 |t.. c| = Ch( d|
|00000b50| 2c 64 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |,d); | |
|00000b60| 20 20 20 20 20 20 20 2f | 2f 20 68 6f 72 69 7a 6f | /|/ horizo|
|00000b70| 6e 74 61 6c 20 63 6f 6e | 63 61 74 65 6e 61 74 69 |ntal con|catenati|
|00000b80| 6f 6e 0d 0a 20 20 20 20 | 63 2e 44 69 73 70 6c 61 |on.. |c.Displa|
|00000b90| 79 4d 61 74 28 29 3b 0d | 0a 20 20 20 20 63 20 3d |yMat();.|. c =|
|00000ba0| 20 43 76 28 20 64 2c 64 | 29 3b 20 20 20 20 20 20 | Cv( d,d|); |
|00000bb0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2f 2f 20 | | // |
|00000bc0| 76 65 72 74 69 63 61 6c | 20 63 6f 6e 74 61 74 65 |vertical| contate|
|00000bd0| 6e 61 74 69 6f 6e 0d 0a | 20 20 20 20 63 2e 44 69 |nation..| c.Di|
|00000be0| 73 70 6c 61 79 4d 61 74 | 28 29 3b 0d 0a 20 20 20 |splayMat|();.. |
|00000bf0| 20 63 20 3d 20 4b 72 6f | 6e 28 49 64 65 6e 74 28 | c = Kro|n(Ident(|
|00000c00| 20 33 20 29 2c 64 29 3b | 20 20 20 20 20 20 20 20 | 3 ),d);| |
|00000c10| 20 2f 2f 20 4b 72 6f 6e | 69 6b 65 72 27 73 20 70 | // Kron|iker's p|
|00000c20| 72 6f 64 75 63 74 0d 0a | 20 20 20 20 53 65 74 64 |roduct..| Setd|
|00000c30| 65 63 28 20 31 20 29 3b | 20 20 20 20 20 20 20 20 |ec( 1 );| |
|00000c40| 20 20 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 73 | | // s|
|00000c50| 65 74 20 6e 75 6d 62 65 | 72 20 6f 66 20 64 65 63 |et numbe|r of dec|
|00000c60| 69 6d 61 6c 73 20 74 6f | 20 70 72 69 6e 74 0d 0a |imals to| print..|
|00000c70| 20 20 20 20 53 65 74 77 | 69 64 28 20 35 20 29 3b | Setw|id( 5 );|
|00000c80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000c90| 20 20 20 20 2f 2f 20 73 | 65 74 20 70 72 69 6e 74 | // s|et print|
|00000ca0| 20 77 69 64 74 68 0d 0a | 20 20 20 20 63 2e 44 69 | width..| c.Di|
|00000cb0| 73 70 6c 61 79 4d 61 74 | 28 29 3b 0d 0a 20 20 20 |splayMat|();.. |
|00000cc0| 20 48 20 3d 20 48 65 6c | 6d 28 20 34 20 29 3b 20 | H = Hel|m( 4 ); |
|00000cd0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000ce0| 20 2f 2f 20 6d 61 6b 65 | 20 61 20 68 65 6c 6d 65 | // make| a helme|
|00000cf0| 72 74 20 6d 61 74 72 69 | 78 0d 0a 20 20 20 20 48 |rt matri|x.. H|
|00000d00| 2e 44 69 73 70 6c 61 79 | 4d 61 74 28 29 3b 20 20 |.Display|Mat(); |
|00000d10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2f | | /|
|00000d20| 2f 20 73 68 6f 77 20 74 | 68 61 74 20 69 74 20 69 |/ show t|hat it i|
|00000d30| 73 20 6f 72 74 68 6f 6e | 6f 72 6d 61 6c 0d 0a 20 |s orthon|ormal.. |
|00000d40| 20 20 20 64 20 3d 20 54 | 72 61 6e 28 20 48 20 29 | d = T|ran( H )|
|00000d50| 2a 48 3b 0d 0a 20 20 20 | 20 64 2e 44 69 73 70 6c |*H;.. | d.Displ|
|00000d60| 61 79 4d 61 74 28 29 3b | 0d 0a 20 20 20 20 64 20 |ayMat();|.. d |
|00000d70| 3d 20 48 2a 54 72 61 6e | 28 48 29 3b 0d 0a 20 20 |= H*Tran|(H);.. |
|00000d80| 20 20 64 2e 44 69 73 70 | 6c 61 79 4d 61 74 28 29 | d.Disp|layMat()|
|00000d90| 3b 0d 0a 0d 0a 20 20 20 | 20 61 20 3d 20 49 64 65 |;.... | a = Ide|
|00000da0| 6e 74 28 20 34 20 29 20 | 2b 20 46 69 6c 6c 28 34 |nt( 4 ) |+ Fill(4|
|00000db0| 2c 34 2c 30 2e 35 29 3b | 20 2f 2f 20 72 65 64 65 |,4,0.5);| // rede|
|00000dc0| 66 69 6e 65 20 61 0d 0a | 20 20 20 20 61 2e 44 69 |fine a..| a.Di|
|00000dd0| 73 70 6c 61 79 4d 61 74 | 28 29 3b 20 20 20 20 20 |splayMat|(); |
|00000de0| 20 20 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 70 | | // p|
|00000df0| 72 69 6e 74 20 61 0d 0a | 20 20 20 20 63 20 3d 20 |rint a..| c = |
|00000e00| 54 72 61 6e 28 48 65 6c | 6d 28 34 29 29 2a 61 2a |Tran(Hel|m(4))*a*|
|00000e10| 48 3b 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 64 |H; | // d|
|00000e20| 69 61 67 6f 6e 61 6c 69 | 7a 65 20 61 0d 0a 20 20 |iagonali|ze a.. |
|00000e30| 20 20 63 2e 44 69 73 70 | 6c 61 79 4d 61 74 28 29 | c.Disp|layMat()|
|00000e40| 3b 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |; | |
|00000e50| 20 20 2f 2f 20 64 69 73 | 70 6c 61 79 20 63 0d 0a | // dis|play c..|
|00000e60| 20 20 20 20 63 20 3d 20 | 49 6e 76 28 20 61 20 29 | c = |Inv( a )|
|00000e70| 3b 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |; | |
|00000e80| 20 20 20 20 2f 2f 20 69 | 6e 76 65 72 74 20 61 0d | // i|nvert a.|
|00000e90| 0a 20 20 20 20 63 2e 44 | 69 73 70 6c 61 79 4d 61 |. c.D|isplayMa|
|00000ea0| 74 28 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |t(); | |
|00000eb0| 20 20 20 20 20 2f 2f 20 | 64 69 73 70 6c 61 79 20 | // |display |
|00000ec0| 63 0d 0a 20 20 20 20 56 | 4d 61 74 72 69 78 20 62 |c.. V|Matrix b|
|00000ed0| 3b 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |; | |
|00000ee0| 20 20 20 20 20 20 20 2f | 2f 20 63 72 65 61 74 65 | /|/ create|
|00000ef0| 20 62 0d 0a 20 20 20 20 | 62 20 3d 20 63 2a 61 3b | b.. |b = c*a;|
|00000f00| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000f10| 20 20 20 20 20 20 20 20 | 2f 2f 20 72 65 64 65 66 | |// redef|
|00000f20| 69 6e 65 20 62 0d 0a 20 | 20 20 20 62 2e 44 69 73 |ine b.. | b.Dis|
|00000f30| 70 6c 61 79 4d 61 74 28 | 29 3b 20 20 20 20 20 20 |playMat(|); |
|00000f40| 20 20 20 20 20 20 20 20 | 20 20 20 2f 2f 20 64 69 | | // di|
|00000f50| 73 70 6c 61 79 20 62 0d | 0a 20 20 20 20 63 20 3d |splay b.|. c =|
|00000f60| 20 66 75 6e 63 74 69 6f | 6e 31 28 61 2c 61 29 3b | functio|n1(a,a);|
|00000f70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2f 2f 20 | | // |
|00000f80| 63 61 6c 6c 20 61 20 66 | 75 6e 63 74 69 6f 6e 0d |call a f|unction.|
|00000f90| 0a 20 20 20 20 63 2e 44 | 69 73 70 6c 61 79 4d 61 |. c.D|isplayMa|
|00000fa0| 74 28 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 20 |t(); | |
|00000fb0| 20 20 20 20 20 2f 2f 20 | 64 69 73 70 6c 61 79 20 | // |display |
|00000fc0| 63 0d 0a 20 20 20 20 44 | 69 73 70 61 74 63 68 2d |c.. D|ispatch-|
|00000fd0| 3e 50 75 73 68 28 63 29 | 3b 0d 0a 20 20 20 20 72 |>Push(c)|;.. r|
|00000fe0| 65 74 75 72 6e 20 44 69 | 73 70 61 74 63 68 2d 3e |eturn Di|spatch->|
|00000ff0| 44 65 63 52 65 74 75 72 | 6e 28 29 3b 0d 0a 7d 0d |DecRetur|n();..}.|
|00001000| 0a 0d 0a 56 4d 61 74 72 | 69 78 20 26 72 65 67 72 |...VMatr|ix ®r|
|00001010| 65 73 73 69 6f 6e 28 20 | 76 6f 69 64 20 29 20 20 |ession( |void ) |
|00001020| 2f 2f 20 64 6f 20 61 20 | 6d 75 6c 74 69 70 6c 65 |// do a |multiple|
|00001030| 20 6c 69 6e 65 61 72 20 | 72 65 67 72 65 73 73 69 | linear |regressi|
|00001040| 6f 6e 0d 0a 20 20 7b 0d | 0a 20 20 20 20 20 20 44 |on.. {.|. D|
|00001050| 69 73 70 61 74 63 68 2d | 3e 49 6e 63 6c 65 76 65 |ispatch-|>Incleve|
|00001060| 6c 28 29 3b 0d 0a 20 20 | 20 20 20 20 56 4d 61 74 |l();.. | VMat|
|00001070| 72 69 78 20 61 2c 78 79 | 2c 72 65 67 3b 0d 0a 20 |rix a,xy|,reg;.. |
|00001080| 20 20 20 20 20 61 20 3d | 20 52 65 61 64 61 28 20 | a =| Reada( |
|00001090| 22 63 61 74 63 68 76 2e | 64 61 74 22 29 3b 20 20 |"catchv.|dat"); |
|000010a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2f | | /|
|000010b0| 2f 20 72 65 61 64 20 64 | 61 74 61 0d 0a 20 20 20 |/ read d|ata.. |
|000010c0| 20 20 20 56 4d 61 74 72 | 69 78 20 62 20 3d 20 56 | VMatr|ix b = V|
|000010d0| 4d 61 74 72 69 78 28 22 | 62 22 2c 20 61 2e 72 2c |Matrix("|b", a.r,|
|000010e0| 20 61 2e 63 29 3b 20 20 | 20 20 20 20 20 2f 2f 73 | a.c); | //s|
|000010f0| 69 6d 70 6c 69 66 79 20 | 69 6e 64 65 78 65 73 0d |implify |indexes.|
|00001100| 0a 20 20 20 20 20 20 69 | 6e 74 20 4e 3d 61 2e 72 |. i|nt N=a.r|
|00001110| 3b 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |; | |
|00001120| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001130| 20 2f 2f 20 4e 0d 0a 20 | 20 20 20 20 20 69 6e 74 | // N.. | int|
|00001140| 20 70 3d 33 3b 20 20 20 | 20 20 20 20 20 20 20 20 | p=3; | |
|00001150| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001160| 20 20 20 20 20 20 20 2f | 2f 20 70 61 72 61 6d 73 | /|/ params|
|00001170| 0d 0a 20 20 20 20 20 20 | 78 79 20 3d 20 43 68 28 |.. |xy = Ch(|
|00001180| 20 46 69 6c 6c 28 4e 2c | 31 2c 31 29 2c 20 53 75 | Fill(N,|1,1), Su|
|00001190| 62 6d 61 74 28 61 2c 4e | 2c 34 2c 31 2c 32 29 20 |bmat(a,N|,4,1,2) |
|000011a0| 29 3b 2f 2f 20 6d 61 6b | 65 20 78 20 79 0d 0a 20 |);// mak|e x y.. |
|000011b0| 20 20 20 20 20 72 65 67 | 20 3d 20 53 77 65 65 70 | reg| = Sweep|
|000011c0| 28 20 31 2c 70 2c 20 54 | 72 61 6e 28 78 79 29 2a |( 1,p, T|ran(xy)*|
|000011d0| 78 79 29 3b 20 20 20 20 | 20 20 20 20 20 20 20 2f |xy); | /|
|000011e0| 2f 20 73 6f 6c 76 65 20 | 72 65 67 72 65 73 73 69 |/ solve |regressi|
|000011f0| 6f 6e 0d 0a 09 09 09 09 | 09 09 2f 2f 20 75 73 69 |on......|..// usi|
|00001200| 6e 67 20 73 77 65 65 70 | 0d 0a 20 20 20 20 20 20 |ng sweep|.. |
|00001210| 72 65 67 2e 4d 28 70 2b | 31 2c 70 2b 31 29 20 3d |reg.M(p+|1,p+1) =|
|00001220| 20 72 65 67 2e 6d 28 70 | 2b 31 2c 70 2b 31 29 2f | reg.m(p|+1,p+1)/|
|00001230| 28 28 64 6f 75 62 6c 65 | 20 29 28 20 4e 2d 70 29 |((double| )( N-p)|
|00001240| 29 3b 0d 0a 09 09 09 09 | 09 09 2f 2f 20 64 69 76 |);......|..// div|
|00001250| 69 64 65 20 74 6f 20 67 | 65 74 20 6d 73 65 0d 0a |ide to g|et mse..|
|00001260| 0d 0a 20 20 20 20 20 20 | 56 4d 61 74 72 69 78 20 |.. |VMatrix |
|00001270| 78 3d 20 53 75 62 6d 61 | 74 28 78 79 2c 4e 2c 33 |x= Subma|t(xy,N,3|
|00001280| 2c 31 2c 31 29 3b 0d 0a | 20 20 20 20 20 20 56 4d |,1,1);..| VM|
|00001290| 61 74 72 69 78 20 79 3d | 20 53 75 62 6d 61 74 28 |atrix y=| Submat(|
|000012a0| 78 79 2c 4e 2c 34 2c 31 | 2c 34 29 3b 0d 0a 20 20 |xy,N,4,1|,4);.. |
|000012b0| 20 20 20 20 56 4d 61 74 | 72 69 78 20 62 65 74 61 | VMat|rix beta|
|000012c0| 68 61 74 20 3d 20 49 6e | 76 28 54 72 61 6e 28 78 |hat = In|v(Tran(x|
|000012d0| 29 2a 78 29 2a 54 72 61 | 6e 28 78 29 2a 79 3b 20 |)*x)*Tra|n(x)*y; |
|000012e0| 2f 2f 73 6f 6c 76 65 20 | 72 65 67 72 65 73 73 69 |//solve |regressi|
|000012f0| 6f 6e 20 75 73 69 6e 67 | 0d 0a 09 09 09 09 09 09 |on using|........|
|00001300| 20 20 2f 2f 6e 6f 72 6d | 61 6c 20 65 71 75 61 74 | //norm|al equat|
|00001310| 69 6f 6e 73 0d 0a 20 20 | 20 20 20 20 62 65 74 61 |ions.. | beta|
|00001320| 68 61 74 2e 44 69 73 70 | 6c 61 79 4d 61 74 28 29 |hat.Disp|layMat()|
|00001330| 3b 0d 0a 20 20 20 20 20 | 20 72 65 67 2e 44 69 73 |;.. | reg.Dis|
|00001340| 70 6c 61 79 4d 61 74 28 | 29 3b 0d 0a 20 20 20 20 |playMat(|);.. |
|00001350| 20 20 44 69 73 70 61 74 | 63 68 2d 3e 50 75 73 68 | Dispat|ch->Push|
|00001360| 28 20 72 65 67 20 29 3b | 20 20 20 20 20 20 20 20 |( reg );| |
|00001370| 20 20 20 20 20 20 20 20 | 20 20 20 20 2f 2f 20 70 | | // p|
|00001380| 75 74 20 72 65 74 75 72 | 6e 20 6d 61 74 20 6f 6e |ut retur|n mat on|
|00001390| 20 73 74 61 63 6b 0d 0a | 20 20 20 20 20 20 72 65 | stack..| re|
|000013a0| 74 75 72 6e 20 44 69 73 | 70 61 74 63 68 2d 3e 44 |turn Dis|patch->D|
|000013b0| 65 63 52 65 74 75 72 6e | 28 29 3b 20 20 20 20 20 |ecReturn|(); |
|000013c0| 20 20 20 20 20 20 20 20 | 2f 2f 20 64 65 63 20 6c | |// dec l|
|000013d0| 65 76 65 6c 20 26 20 72 | 65 74 75 72 6e 0d 0a 7d |evel & r|eturn..}|
|000013e0| 0d 0a 0d 0a 76 6f 69 64 | 20 61 6c 74 65 72 6d 61 |....void| alterma|
|000013f0| 74 72 69 78 28 20 56 4d | 61 74 72 69 78 20 26 74 |trix( VM|atrix &t|
|00001400| 20 29 20 2f 2f 20 61 6c | 74 65 72 20 61 20 6d 61 | ) // al|ter a ma|
|00001410| 74 72 69 78 20 65 6c 65 | 6d 65 6e 74 0d 0a 20 20 |trix ele|ment.. |
|00001420| 7b 0d 0a 20 20 20 20 20 | 74 2e 4d 28 31 2c 31 29 |{.. |t.M(1,1)|
|00001430| 20 3d 20 32 35 32 35 2e | 30 3b 0d 0a 20 20 20 20 | = 2525.|0;.. |
|00001440| 20 56 4d 61 74 72 69 78 | 20 74 65 6d 70 20 3d 20 | VMatrix| temp = |
|00001450| 4d 53 6f 72 74 28 20 74 | 2c 20 31 20 29 3b 0d 0a |MSort( t|, 1 );..|
|00001460| 20 20 20 20 20 74 65 6d | 70 2e 44 69 73 70 6c 61 | tem|p.Displa|
|00001470| 79 4d 61 74 28 29 3b 0d | 0a 20 20 7d 0d 0a 0d 0a |yMat();.|. }....|
|00001480| 76 6f 69 64 20 74 65 73 | 74 68 75 67 65 28 20 76 |void tes|thuge( v|
|00001490| 6f 69 64 20 29 0d 0a 7b | 0d 0a 20 20 56 4d 61 74 |oid )..{|.. VMat|
|000014a0| 72 69 78 20 56 65 72 79 | 42 69 67 28 22 76 62 22 |rix Very|Big("vb"|
|000014b0| 2c 33 30 30 2c 33 30 30 | 29 3b 0d 0a 20 20 56 65 |,300,300|);.. Ve|
|000014c0| 72 79 42 69 67 2e 49 6e | 66 6f 4d 61 74 28 29 3b |ryBig.In|foMat();|
|000014d0| 0d 0a 0d 0a 7d 0d 0a 0d | 0a 76 6f 69 64 20 76 65 |....}...|.void ve|
|000014e0| 72 73 69 6f 6e 31 70 31 | 20 28 20 76 6f 69 64 20 |rsion1p1| ( void |
|000014f0| 29 0d 0a 7b 0d 0a 20 20 | 20 2f 2f 0d 0a 20 20 20 |)..{.. | //.. |
|00001500| 2f 2f 20 74 65 73 74 72 | 65 67 2e 63 70 70 20 66 |// testr|eg.cpp f|
|00001510| 6f 72 20 74 65 73 74 20 | 6f 66 20 72 65 67 72 65 |or test |of regre|
|00001520| 73 73 69 6f 6e 20 66 75 | 6e 63 74 69 6f 6e 73 2e |ssion fu|nctions.|
|00001530| 0d 0a 20 20 20 2f 2f 0d | 0a 20 20 20 44 69 73 70 |.. //.|. Disp|
|00001540| 61 74 63 68 2d 3e 49 6e | 63 6c 65 76 65 6c 28 29 |atch->In|clevel()|
|00001550| 3b 0d 0a 20 20 20 56 4d | 61 74 72 69 78 20 41 3d |;.. VM|atrix A=|
|00001560| 46 69 6c 6c 28 36 2c 36 | 2c 31 2e 30 29 3b 0d 0a |Fill(6,6|,1.0);..|
|00001570| 20 20 20 56 4d 61 74 72 | 69 78 20 42 3d 46 69 6c | VMatr|ix B=Fil|
|00001580| 6c 28 36 2c 31 2c 32 2e | 30 29 3b 0d 0a 0d 0a 20 |l(6,1,2.|0);.... |
|00001590| 20 20 28 56 65 63 28 41 | 29 29 2e 44 69 73 70 6c | (Vec(A|)).Displ|
|000015a0| 61 79 4d 61 74 28 29 3b | 0d 0a 20 20 20 28 56 65 |ayMat();|.. (Ve|
|000015b0| 63 64 69 61 67 28 41 29 | 29 2e 44 69 73 70 6c 61 |cdiag(A)|).Displa|
|000015c0| 79 4d 61 74 28 29 3b 0d | 0a 20 20 20 28 44 69 61 |yMat();.|. (Dia|
|000015d0| 67 28 42 29 29 2e 44 69 | 73 70 6c 61 79 4d 61 74 |g(B)).Di|splayMat|
|000015e0| 28 29 3b 0d 0a 20 20 20 | 28 53 68 61 70 65 28 41 |();.. |(Shape(A|
|000015f0| 2c 33 29 29 2e 44 69 73 | 70 6c 61 79 4d 61 74 28 |,3)).Dis|playMat(|
|00001600| 29 3b 0d 0a 20 20 20 28 | 53 75 6d 28 41 2c 52 4f |);.. (|Sum(A,RO|
|00001610| 57 53 29 29 2e 44 69 73 | 70 6c 61 79 4d 61 74 28 |WS)).Dis|playMat(|
|00001620| 29 3b 0d 0a 20 20 20 28 | 53 75 6d 73 71 28 41 2c |);.. (|Sumsq(A,|
|00001630| 43 4f 4c 55 4d 4e 53 29 | 29 2e 44 69 73 70 6c 61 |COLUMNS)|).Displa|
|00001640| 79 4d 61 74 28 29 3b 0d | 0a 20 20 20 28 43 75 73 |yMat();.|. (Cus|
|00001650| 75 6d 28 41 29 29 2e 44 | 69 73 70 6c 61 79 4d 61 |um(A)).D|isplayMa|
|00001660| 74 28 29 3b 0d 0a 20 20 | 20 28 4d 6d 69 6e 28 42 |t();.. | (Mmin(B|
|00001670| 29 29 2e 44 69 73 70 6c | 61 79 4d 61 74 28 29 3b |)).Displ|ayMat();|
|00001680| 0d 0a 20 20 20 28 4d 6d | 61 78 28 42 2c 52 4f 57 |.. (Mm|ax(B,ROW|
|00001690| 53 29 29 2e 44 69 73 70 | 6c 61 79 4d 61 74 28 29 |S)).Disp|layMat()|
|000016a0| 3b 0d 0a 0d 0a 20 20 20 | 56 4d 61 74 72 69 78 20 |;.... |VMatrix |
|000016b0| 43 20 3d 20 43 68 28 41 | 2c 44 69 61 67 28 42 29 |C = Ch(A|,Diag(B)|
|000016c0| 29 3b 0d 0a 20 20 20 53 | 65 74 77 69 64 28 35 29 |);.. S|etwid(5)|
|000016d0| 3b 0d 0a 20 20 20 53 65 | 74 64 65 63 28 32 29 3b |;.. Se|tdec(2);|
|000016e0| 0d 0a 20 20 20 43 2e 44 | 69 73 70 6c 61 79 4d 61 |.. C.D|isplayMa|
|000016f0| 74 28 29 3b 0d 0a 20 20 | 20 43 72 6f 77 28 20 43 |t();.. | Crow( C|
|00001700| 2c 20 31 2c 20 30 2e 36 | 20 29 3b 0d 0a 20 20 20 |, 1, 0.6| );.. |
|00001710| 43 2e 44 69 73 70 6c 61 | 79 4d 61 74 28 29 3b 0d |C.Displa|yMat();.|
|00001720| 0a 20 20 20 53 72 6f 77 | 28 20 43 2c 20 31 2c 20 |. Srow|( C, 1, |
|00001730| 36 20 29 3b 0d 0a 20 20 | 20 43 2e 44 69 73 70 6c |6 );.. | C.Displ|
|00001740| 61 79 4d 61 74 28 29 3b | 0d 0a 20 20 20 4c 72 6f |ayMat();|.. Lro|
|00001750| 77 28 20 43 2c 20 32 2c | 33 2c 20 2e 35 29 3b 0d |w( C, 2,|3, .5);.|
|00001760| 0a 20 20 20 43 2e 44 69 | 73 70 6c 61 79 4d 61 74 |. C.Di|splayMat|
|00001770| 28 29 3b 0d 0a 20 20 20 | 43 63 6f 6c 28 20 43 2c |();.. |Ccol( C,|
|00001780| 20 31 2c 20 30 2e 36 20 | 29 3b 0d 0a 20 20 20 43 | 1, 0.6 |);.. C|
|00001790| 2e 44 69 73 70 6c 61 79 | 4d 61 74 28 29 3b 0d 0a |.Display|Mat();..|
|000017a0| 20 20 20 53 63 6f 6c 28 | 20 43 2c 20 31 2c 36 20 | Scol(| C, 1,6 |
|000017b0| 29 3b 0d 0a 20 20 20 43 | 2e 44 69 73 70 6c 61 79 |);.. C|.Display|
|000017c0| 4d 61 74 28 29 3b 0d 0a | 20 20 20 4c 63 6f 6c 28 |Mat();..| Lcol(|
|000017d0| 20 43 2c 20 32 2c 33 2c | 20 2e 35 29 3b 0d 0a 20 | C, 2,3,| .5);.. |
|000017e0| 20 20 43 2e 44 69 73 70 | 6c 61 79 4d 61 74 28 29 | C.Disp|layMat()|
|000017f0| 3b 0d 0a 20 20 20 44 69 | 73 70 61 74 63 68 2d 3e |;.. Di|spatch->|
|00001800| 44 65 63 6c 65 76 65 6c | 28 29 3b 0d 0a 7d 0d 0a |Declevel|();..}..|
|00001810| 0d 0a 6d 61 69 6e 28 29 | 0d 0a 7b 0d 0a 20 20 56 |..main()|..{.. V|
|00001820| 4d 61 74 72 69 78 20 74 | 65 73 74 69 74 3b 0d 0a |Matrix t|estit;..|
|00001830| 0d 0a 20 20 70 72 69 6e | 74 66 28 22 25 64 5c 6e |.. prin|tf("%d\n|
|00001840| 22 2c 20 73 69 7a 65 6f | 66 28 74 65 73 74 69 74 |", sizeo|f(testit|
|00001850| 29 20 29 3b 0d 0a 0d 0a | 20 20 74 65 73 74 69 74 |) );....| testit|
|00001860| 20 3d 20 72 65 67 72 65 | 73 73 69 6f 6e 28 29 3b | = regre|ssion();|
|00001870| 0d 0a 20 20 74 65 73 74 | 69 74 2e 44 69 73 70 6c |.. test|it.Displ|
|00001880| 61 79 4d 61 74 28 29 3b | 0d 0a 0d 0a 20 20 56 4d |ayMat();|.... VM|
|00001890| 61 74 72 69 78 20 62 20 | 3d 20 74 65 73 74 69 74 |atrix b |= testit|
|000018a0| 3b 0d 0a 20 20 56 4d 61 | 74 72 69 78 20 61 20 3d |;.. VMa|trix a =|
|000018b0| 20 32 2a 74 65 73 74 69 | 74 3b 0d 0a 0d 0a 20 20 | 2*testi|t;.... |
|000018c0| 56 4d 61 74 72 69 78 20 | 63 20 3d 20 61 20 2b 20 |VMatrix |c = a + |
|000018d0| 62 20 2b 20 61 20 2b 20 | 62 3b 0d 0a 20 20 63 2e |b + a + |b;.. c.|
|000018e0| 44 69 73 70 6c 61 79 4d | 61 74 28 29 3b 0d 0a 0d |DisplayM|at();...|
|000018f0| 0a 20 20 56 4d 61 74 72 | 69 78 20 74 20 3d 20 66 |. VMatr|ix t = f|
|00001900| 75 6e 63 74 69 6f 6e 30 | 28 29 20 2b 20 66 75 6e |unction0|() + fun|
|00001910| 63 74 69 6f 6e 30 28 29 | 3b 0d 0a 20 20 74 2e 44 |ction0()|;.. t.D|
|00001920| 69 73 70 6c 61 79 4d 61 | 74 28 29 3b 0d 0a 0d 0a |isplayMa|t();....|
|00001930| 20 20 61 6c 74 65 72 6d | 61 74 72 69 78 28 20 74 | alterm|atrix( t|
|00001940| 65 73 74 69 74 20 29 3b | 0d 0a 20 20 74 65 73 74 |estit );|.. test|
|00001950| 69 74 2e 44 69 73 70 6c | 61 79 4d 61 74 28 29 3b |it.Displ|ayMat();|
|00001960| 0d 0a 0d 0a 20 20 76 65 | 72 73 69 6f 6e 31 70 31 |.... ve|rsion1p1|
|00001970| 28 29 3b 0d 0a 0d 0a 23 | 69 66 6e 64 65 66 20 49 |();....#|ifndef I|
|00001980| 4e 5f 52 41 4d 0d 0a 20 | 20 20 20 20 20 20 74 65 |N_RAM.. | te|
|00001990| 73 74 68 75 67 65 28 29 | 3b 0d 0a 23 65 6e 64 69 |sthuge()|;..#endi|
|000019a0| 66 0d 0a 0d 0a 20 20 76 | 63 6c 6f 73 65 28 29 3b |f.... v|close();|
|000019b0| 0d 0a 7d 0d 0a 0d 0a 5c | 65 6e 64 7b 76 65 72 62 |..}....\|end{verb|
|000019c0| 61 74 69 6d 7d | |atim} | |
+--------+-------------------------+-------------------------+--------+--------+