/ Mac Power 1997 January / macpower199701.bin / AMUG / Programming_10 / Videotoolbox 96-08-22-c.sit / VideoToolbox / VideoToolboxSources / ReadMATLABFile.c

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 10 52 65 61 64 4d 41 | 54 4c 41 42 46 69 6c 65 |..ReadMA|TLABFile|
|00000010| 2e 63 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |.c......|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 54 45 58 54 43 57 49 | 45 01 00 00 00 00 00 00 |.TEXTCWI|E.......|
|00000050| 00 00 00 00 00 30 91 00 | 00 01 c8 ac 34 24 64 ac |.....0..|....4$d.|
|00000060| 34 24 64 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |4$d.....|........|
|00000070| 00 00 00 00 00 00 00 00 | 00 00 81 81 4a 9a 00 00 |........|....J...|
|00000080| 2f 2a 0d 52 65 61 64 4d | 41 54 4c 41 42 46 69 6c |/*.ReadM|ATLABFil|
|00000090| 65 2e 63 20 2d 20 6c 6f | 61 64 2f 73 61 76 65 20 |e.c - lo|ad/save |
|000000a0| 61 20 6d 61 74 72 69 78 | 20 66 72 6f 6d 2f 74 6f |a matrix| from/to|
|000000b0| 20 61 20 4d 41 54 4c 41 | 42 20 64 61 74 61 20 66 | a MATLA|B data f|
|000000c0| 69 6c 65 2e 0d 54 68 65 | 73 65 20 66 75 6e 63 74 |ile..The|se funct|
|000000d0| 69 6f 6e 73 20 61 6c 6c | 6f 77 20 61 6e 79 20 43 |ions all|ow any C|
|000000e0| 20 70 72 6f 67 72 61 6d | 20 74 6f 20 72 65 61 64 | program| to read|
|000000f0| 20 61 6e 64 20 77 72 69 | 74 65 20 4d 41 54 4c 41 | and wri|te MATLA|
|00000100| 42 20 64 61 74 61 20 66 | 69 6c 65 73 2e 0d 0d 69 |B data f|iles...i|
|00000110| 6e 74 20 53 61 76 65 4d | 61 74 44 6f 75 62 6c 65 |nt SaveM|atDouble|
|00000120| 73 28 46 49 4c 45 20 2a | 66 2c 63 68 61 72 20 2a |s(FILE *|f,char *|
|00000130| 6e 61 6d 65 2c 6c 6f 6e | 67 20 72 6f 77 73 2c 6c |name,lon|g rows,l|
|00000140| 6f 6e 67 20 63 6f 6c 73 | 0d 09 2c 64 6f 75 62 6c |ong cols|..,doubl|
|00000150| 65 20 2a 72 65 61 6c 2c | 64 6f 75 62 6c 65 20 2a |e *real,|double *|
|00000160| 69 6d 61 67 29 3b 09 09 | 09 2f 2f 20 77 72 69 74 |imag);..|.// writ|
|00000170| 65 20 64 6f 75 62 6c 65 | 73 20 61 73 20 38 2d 62 |e double|s as 8-b|
|00000180| 79 74 65 20 44 4f 55 42 | 4c 45 53 0d 69 6e 74 20 |yte DOUB|LES.int |
|00000190| 53 61 76 65 4d 61 74 53 | 68 6f 72 74 73 28 46 49 |SaveMatS|horts(FI|
|000001a0| 4c 45 20 2a 66 2c 63 68 | 61 72 20 2a 6e 61 6d 65 |LE *f,ch|ar *name|
|000001b0| 2c 6c 6f 6e 67 20 72 6f | 77 73 2c 6c 6f 6e 67 20 |,long ro|ws,long |
|000001c0| 63 6f 6c 73 0d 09 2c 73 | 68 6f 72 74 20 2a 72 65 |cols..,s|hort *re|
|000001d0| 61 6c 2c 73 68 6f 72 74 | 20 2a 69 6d 61 67 29 3b |al,short| *imag);|
|000001e0| 09 09 09 09 2f 2f 20 77 | 72 69 74 65 20 73 68 6f |....// w|rite sho|
|000001f0| 72 74 73 20 61 73 20 73 | 68 6f 72 74 73 0d 69 6e |rts as s|horts.in|
|00000200| 74 20 4c 6f 61 64 4e 65 | 78 74 4d 61 74 49 6e 74 |t LoadNe|xtMatInt|
|00000210| 6f 44 6f 75 62 6c 65 73 | 28 46 49 4c 45 20 2a 66 |oDoubles|(FILE *f|
|00000220| 2c 63 68 61 72 20 2a 6e | 61 6d 65 2c 6c 6f 6e 67 |,char *n|ame,long|
|00000230| 20 2a 72 6f 77 73 2c 6c | 6f 6e 67 20 2a 63 6f 6c | *rows,l|ong *col|
|00000240| 73 0d 09 2c 64 6f 75 62 | 6c 65 20 2a 2a 72 65 61 |s..,doub|le **rea|
|00000250| 6c 2c 64 6f 75 62 6c 65 | 20 2a 2a 69 6d 61 67 29 |l,double| **imag)|
|00000260| 3b 09 09 09 2f 2f 20 72 | 65 61 64 20 38 2d 62 79 |;...// r|ead 8-by|
|00000270| 74 65 20 44 4f 55 42 4c | 45 53 20 69 6e 74 6f 20 |te DOUBL|ES into |
|00000280| 64 6f 75 62 6c 65 73 0d | 69 6e 74 20 4c 6f 61 64 |doubles.|int Load|
|00000290| 4e 65 78 74 4d 61 74 49 | 6e 74 6f 53 68 6f 72 74 |NextMatI|ntoShort|
|000002a0| 73 28 46 49 4c 45 20 2a | 66 2c 63 68 61 72 20 2a |s(FILE *|f,char *|
|000002b0| 6e 61 6d 65 2c 6c 6f 6e | 67 20 2a 72 6f 77 73 2c |name,lon|g *rows,|
|000002c0| 6c 6f 6e 67 20 2a 63 6f | 6c 73 0d 09 2c 73 68 6f |long *co|ls..,sho|
|000002d0| 72 74 20 2a 2a 72 65 61 | 6c 2c 73 68 6f 72 74 20 |rt **rea|l,short |
|000002e0| 2a 2a 69 6d 61 67 29 3b | 09 09 09 2f 2f 20 72 65 |**imag);|...// re|
|000002f0| 61 64 20 38 2d 62 79 74 | 65 20 44 4f 55 42 4c 45 |ad 8-byt|e DOUBLE|
|00000300| 53 20 69 6e 74 6f 20 73 | 68 6f 72 74 73 0d 69 6e |S into s|horts.in|
|00000310| 74 20 4c 6f 61 64 4e 61 | 6d 65 64 4d 61 74 49 6e |t LoadNa|medMatIn|
|00000320| 74 6f 44 6f 75 62 6c 65 | 73 28 46 49 4c 45 20 2a |toDouble|s(FILE *|
|00000330| 66 2c 63 68 61 72 20 2a | 6e 61 6d 65 2c 6c 6f 6e |f,char *|name,lon|
|00000340| 67 20 2a 72 6f 77 73 2c | 6c 6f 6e 67 20 2a 63 6f |g *rows,|long *co|
|00000350| 6c 73 0d 09 2c 64 6f 75 | 62 6c 65 20 2a 2a 72 65 |ls..,dou|ble **re|
|00000360| 61 6c 2c 64 6f 75 62 6c | 65 20 2a 2a 69 6d 61 67 |al,doubl|e **imag|
|00000370| 29 3b 09 09 09 2f 2f 20 | 72 65 61 64 20 38 2d 62 |);...// |read 8-b|
|00000380| 79 74 65 20 44 4f 55 42 | 4c 45 53 20 69 6e 74 6f |yte DOUB|LES into|
|00000390| 20 64 6f 75 62 6c 65 73 | 0d 69 6e 74 20 4c 6f 61 | doubles|.int Loa|
|000003a0| 64 4e 61 6d 65 64 4d 61 | 74 49 6e 74 6f 53 68 6f |dNamedMa|tIntoSho|
|000003b0| 72 74 73 28 46 49 4c 45 | 20 2a 66 2c 63 68 61 72 |rts(FILE| *f,char|
|000003c0| 20 2a 6e 61 6d 65 2c 6c | 6f 6e 67 20 2a 72 6f 77 | *name,l|ong *row|
|000003d0| 73 2c 6c 6f 6e 67 20 2a | 63 6f 6c 73 0d 09 2c 73 |s,long *|cols..,s|
|000003e0| 68 6f 72 74 20 2a 2a 72 | 65 61 6c 2c 73 68 6f 72 |hort **r|eal,shor|
|000003f0| 74 20 2a 2a 69 6d 61 67 | 29 3b 09 09 09 2f 2f 20 |t **imag|);...// |
|00000400| 72 65 61 64 20 38 2d 62 | 79 74 65 20 44 4f 55 42 |read 8-b|yte DOUB|
|00000410| 4c 45 53 20 69 6e 74 6f | 20 73 68 6f 72 74 73 0d |LES into| shorts.|
|00000420| 0d 41 6c 6c 20 74 68 65 | 20 72 6f 75 74 69 6e 65 |.All the| routine|
|00000430| 73 20 61 73 73 75 6d 65 | 20 61 6e 20 6f 70 65 6e |s assume| an open|
|00000440| 20 66 69 6c 65 2c 20 61 | 6e 64 20 72 65 74 75 72 | file, a|nd retur|
|00000450| 6e 20 61 6e 20 65 72 72 | 6f 72 20 69 66 20 74 68 |n an err|or if th|
|00000460| 65 79 20 68 69 74 20 65 | 6e 64 20 6f 66 0d 66 69 |ey hit e|nd of.fi|
|00000470| 6c 65 2e 20 54 68 65 20 | 53 61 76 65 58 58 58 20 |le. The |SaveXXX |
|00000480| 72 6f 75 74 69 6e 65 73 | 20 61 70 70 65 6e 64 20 |routines| append |
|00000490| 61 20 6d 61 74 72 69 78 | 20 74 6f 20 74 68 65 20 |a matrix| to the |
|000004a0| 65 78 69 73 74 69 6e 67 | 20 66 69 6c 65 2e 20 54 |existing| file. T|
|000004b0| 68 65 20 4c 6f 61 64 4e | 65 78 74 58 58 58 0d 72 |he LoadN|extXXX.r|
|000004c0| 6f 75 74 69 6e 65 73 20 | 72 65 61 64 20 74 68 65 |outines |read the|
|000004d0| 20 6e 65 78 74 20 6d 61 | 74 72 69 78 20 66 72 6f | next ma|trix fro|
|000004e0| 6d 20 74 68 65 20 66 69 | 6c 65 2c 20 61 64 76 61 |m the fi|le, adva|
|000004f0| 6e 63 69 6e 67 20 74 68 | 65 20 66 69 6c 65 20 70 |ncing th|e file p|
|00000500| 6f 69 6e 74 65 72 2e 0d | 53 61 76 65 58 58 58 20 |ointer..|SaveXXX |
|00000510| 61 6e 64 20 4c 6f 61 64 | 4e 65 78 74 58 58 58 20 |and Load|NextXXX |
|00000520| 63 61 6e 20 62 65 20 63 | 61 6c 6c 65 64 20 72 65 |can be c|alled re|
|00000530| 70 65 61 74 65 64 6c 79 | 20 75 6e 74 69 6c 20 74 |peatedly| until t|
|00000540| 68 65 79 20 66 61 69 6c | 20 28 65 2e 67 2e 20 65 |hey fail| (e.g. e|
|00000550| 6e 64 20 6f 66 0d 66 69 | 6c 65 20 69 73 20 72 65 |nd of.fi|le is re|
|00000560| 61 63 68 65 64 29 2e 20 | 54 68 65 20 4c 6f 61 64 |ached). |The Load|
|00000570| 4e 61 6d 65 64 58 58 58 | 20 72 6f 75 74 69 6e 65 |NamedXXX| routine|
|00000580| 73 20 66 69 72 73 74 20 | 72 65 77 69 6e 64 20 74 |s first |rewind t|
|00000590| 68 65 20 66 69 6c 65 20 | 61 6e 64 20 74 68 65 6e |he file |and then|
|000005a0| 20 73 63 61 6e 0d 69 74 | 20 75 6e 74 69 6c 20 74 | scan.it| until t|
|000005b0| 68 65 79 20 66 69 6e 64 | 20 61 20 6d 61 74 72 69 |hey find| a matri|
|000005c0| 78 20 77 68 6f 73 65 20 | 6e 61 6d 65 20 6d 61 74 |x whose |name mat|
|000005d0| 63 68 65 73 20 74 68 61 | 74 20 73 75 70 70 6c 69 |ches tha|t suppli|
|000005e0| 65 64 2c 20 6f 72 20 65 | 6e 64 20 6f 66 20 66 69 |ed, or e|nd of fi|
|000005f0| 6c 65 20 69 73 0d 72 65 | 61 63 68 65 64 2e 20 54 |le is.re|ached. T|
|00000600| 68 65 20 66 69 6c 65 20 | 70 6f 69 6e 74 65 72 20 |he file |pointer |
|00000610| 69 73 20 6c 65 66 74 20 | 61 74 20 74 68 65 20 65 |is left |at the e|
|00000620| 6e 64 20 6f 66 20 74 68 | 61 74 20 6d 61 74 72 69 |nd of th|at matri|
|00000630| 78 2e 20 43 61 6c 6c 69 | 6e 67 20 4c 6f 61 64 4e |x. Calli|ng LoadN|
|00000640| 61 6d 65 58 58 58 0d 72 | 65 70 65 61 74 65 64 6c |ameXXX.r|epeatedl|
|00000650| 79 20 77 69 6c 6c 20 72 | 65 61 64 20 74 68 65 20 |y will r|ead the |
|00000660| 73 61 6d 65 20 6d 61 74 | 72 69 78 20 65 76 65 72 |same mat|rix ever|
|00000670| 79 20 74 69 6d 65 3a 20 | 74 68 65 20 66 69 72 73 |y time: |the firs|
|00000680| 74 20 69 6e 73 74 61 6e | 63 65 20 6f 66 20 74 68 |t instan|ce of th|
|00000690| 65 20 6e 61 6d 65 64 0d | 6d 61 74 72 69 78 2e 0d |e named.|matrix..|
|000006a0| 0d 54 68 65 20 4c 6f 61 | 64 58 58 58 20 61 6e 64 |.The Loa|dXXX and|
|000006b0| 20 53 61 76 65 58 58 58 | 20 72 6f 75 74 69 6e 65 | SaveXXX| routine|
|000006c0| 73 20 72 65 74 75 72 6e | 20 30 20 69 66 20 73 75 |s return| 0 if su|
|000006d0| 63 63 65 73 73 66 75 6c | 20 61 6e 64 20 61 20 70 |ccessful| and a p|
|000006e0| 6f 73 69 74 69 76 65 20 | 6e 6f 6e 7a 65 72 6f 0d |ositive |nonzero.|
|000006f0| 69 6e 74 65 67 65 72 20 | 69 66 20 74 68 65 79 20 |integer |if they |
|00000700| 66 61 69 6c 2e 20 20 54 | 68 65 20 74 77 6f 20 61 |fail.  T|he two a|
|00000710| 72 72 61 79 20 70 6f 69 | 6e 74 65 72 73 20 77 69 |rray poi|nters wi|
|00000720| 6c 6c 20 62 65 20 4e 55 | 4c 4c 20 61 66 74 65 72 |ll be NU|LL after|
|00000730| 20 63 61 6c 6c 69 6e 67 | 20 4c 6f 61 64 58 58 58 | calling| LoadXXX|
|00000740| 0d 75 6e 6c 65 73 73 20 | 74 68 65 79 20 61 72 65 |.unless |they are|
|00000750| 20 6c 6f 61 64 65 64 2e | 20 54 68 65 20 73 65 63 | loaded.| The sec|
|00000760| 6f 6e 64 20 28 69 6d 61 | 67 69 6e 61 72 79 29 20 |ond (ima|ginary) |
|00000770| 70 6f 69 6e 74 65 72 20 | 77 69 6c 6c 20 62 65 20 |pointer |will be |
|00000780| 4e 55 4c 4c 20 75 6e 6c | 65 73 73 20 74 68 65 0d |NULL unl|ess the.|
|00000790| 6d 61 74 72 69 78 20 77 | 61 73 20 63 6f 6d 70 6c |matrix w|as compl|
|000007a0| 65 78 2e 20 53 69 6d 69 | 6c 61 72 6c 79 2c 20 77 |ex. Simi|larly, w|
|000007b0| 68 65 6e 20 79 6f 75 20 | 63 61 6c 6c 20 53 61 76 |hen you |call Sav|
|000007c0| 65 58 58 58 20 79 6f 75 | 20 69 6e 64 69 63 61 74 |eXXX you| indicat|
|000007d0| 65 20 61 20 72 65 61 6c | 20 6d 61 74 72 69 78 0d |e a real| matrix.|
|000007e0| 62 79 20 73 75 70 70 6c | 79 69 6e 67 20 61 20 4e |by suppl|ying a N|
|000007f0| 55 4c 4c 20 73 65 63 6f | 6e 64 20 70 6f 69 6e 74 |ULL seco|nd point|
|00000800| 65 72 2e 0d 0d 54 68 65 | 72 65 20 69 73 20 61 20 |er...The|re is a |
|00000810| 73 6c 69 67 68 74 20 61 | 73 79 6d 6d 65 74 72 79 |slight a|symmetry|
|00000820| 20 68 65 72 65 20 74 68 | 61 74 20 79 6f 75 20 73 | here th|at you s|
|00000830| 68 6f 75 6c 64 20 62 65 | 20 61 77 61 72 65 20 6f |hould be| aware o|
|00000840| 66 2e 20 53 61 76 65 4d | 61 74 53 68 6f 72 74 73 |f. SaveM|atShorts|
|00000850| 28 29 0d 70 72 6f 64 75 | 63 65 73 20 61 20 28 63 |().produ|ces a (c|
|00000860| 6f 6d 70 61 63 74 29 20 | 4d 41 54 4c 41 42 20 66 |ompact) |MATLAB f|
|00000870| 69 6c 65 20 6d 61 64 65 | 20 75 70 20 6f 66 20 73 |ile made| up of s|
|00000880| 68 6f 72 74 73 2c 20 62 | 75 74 20 74 68 65 0d 4c |horts, b|ut the.L|
|00000890| 6f 61 64 58 58 58 4d 61 | 74 49 6e 74 6f 53 68 6f |oadXXXMa|tIntoSho|
|000008a0| 72 74 73 28 29 20 77 69 | 6c 6c 20 6f 6e 6c 79 20 |rts() wi|ll only |
|000008b0| 72 65 61 64 20 61 20 28 | 62 69 67 29 20 4d 41 54 |read a (|big) MAT|
|000008c0| 4c 41 42 20 66 69 6c 65 | 20 6d 61 64 65 20 75 70 |LAB file| made up|
|000008d0| 20 6f 66 20 38 2d 62 79 | 74 65 0d 64 6f 75 62 6c | of 8-by|te.doubl|
|000008e0| 65 73 2e 20 53 6f 20 79 | 6f 75 20 77 6f 6e 27 74 |es. So y|ou won't|
|000008f0| 20 61 6c 77 61 79 73 20 | 62 65 20 61 62 6c 65 20 | always |be able |
|00000900| 74 6f 20 72 65 61 64 20 | 66 72 6f 6d 20 43 20 74 |to read |from C t|
|00000910| 68 65 20 66 69 6c 65 20 | 79 6f 75 20 77 72 6f 74 |he file |you wrot|
|00000920| 65 20 66 72 6f 6d 20 43 | 2e 0d 54 68 65 20 72 65 |e from C|..The re|
|00000930| 61 73 6f 6e 20 66 6f 72 | 20 74 68 69 73 20 69 73 |ason for| this is|
|00000940| 20 74 68 61 74 20 77 68 | 69 6c 65 20 74 68 65 20 | that wh|ile the |
|00000950| 4d 41 54 4c 41 42 20 22 | 4c 6f 61 64 22 20 63 6f |MATLAB "|Load" co|
|00000960| 6d 6d 61 6e 64 20 73 75 | 70 70 6f 72 74 73 20 61 |mmand su|pports a|
|00000970| 6c 6c 20 74 68 65 0d 4d | 41 54 4c 41 42 20 66 69 |ll the.M|ATLAB fi|
|00000980| 6c 65 20 66 6f 72 6d 61 | 74 73 2c 20 74 68 65 20 |le forma|ts, the |
|00000990| 4d 41 54 4c 41 42 20 22 | 53 61 76 65 22 20 63 6f |MATLAB "|Save" co|
|000009a0| 6d 6d 61 6e 64 20 68 61 | 73 20 6e 6f 20 70 72 6f |mmand ha|s no pro|
|000009b0| 76 69 73 69 6f 6e 20 66 | 6f 72 20 73 70 65 63 69 |vision f|or speci|
|000009c0| 66 79 69 6e 67 0d 28 62 | 69 6e 61 72 79 29 20 66 |fying.(b|inary) f|
|000009d0| 69 6c 65 20 66 6f 72 6d | 61 74 73 20 61 6e 64 20 |ile form|ats and |
|000009e0| 74 68 75 73 20 61 6c 77 | 61 79 73 20 75 73 65 73 |thus alw|ays uses|
|000009f0| 20 66 6f 72 6d 61 74 20 | 74 79 70 65 20 30 2c 20 | format |type 0, |
|00000a00| 38 2d 62 79 74 65 20 73 | 68 6f 72 74 20 64 6f 75 |8-byte s|hort dou|
|00000a10| 62 6c 65 73 2c 0d 73 6f | 20 49 20 64 69 64 6e 27 |bles,.so| I didn'|
|00000a20| 74 20 62 6f 74 68 65 72 | 20 74 6f 20 69 6d 70 6c |t bother| to impl|
|00000a30| 65 6d 65 6e 74 20 61 6e | 79 20 72 6f 75 74 69 6e |ement an|y routin|
|00000a40| 65 73 20 74 6f 20 72 65 | 61 64 20 61 6e 79 20 6f |es to re|ad any o|
|00000a50| 74 68 65 72 20 66 6f 72 | 6d 61 74 2e 20 28 49 74 |ther for|mat. (It|
|00000a60| 20 77 6f 75 6c 64 0d 62 | 65 20 65 61 73 79 20 74 | would.b|e easy t|
|00000a70| 6f 20 61 64 64 20 74 68 | 69 73 20 63 61 70 61 62 |o add th|is capab|
|00000a80| 69 6c 69 74 79 2c 20 69 | 66 20 69 74 20 77 65 72 |ility, i|f it wer|
|00000a90| 65 20 6e 65 65 64 65 64 | 2e 29 0d 0d 45 58 41 4d |e needed|.)..EXAM|
|00000aa0| 50 4c 45 53 3a 0d 54 68 | 69 73 20 65 78 61 6d 70 |PLES:.Th|is examp|
|00000ab0| 6c 65 20 6c 6f 61 64 73 | 20 74 68 65 20 66 69 72 |le loads| the fir|
|00000ac0| 73 74 20 74 77 6f 20 6d | 61 74 72 69 63 65 73 20 |st two m|atrices |
|00000ad0| 66 72 6f 6d 20 61 20 4d | 41 54 4c 41 42 20 66 69 |from a M|ATLAB fi|
|00000ae0| 6c 65 3a 0d 0d 09 46 49 | 4c 45 20 2a 66 3b 0d 09 |le:...FI|LE *f;..|
|00000af0| 63 68 61 72 20 6e 61 6d | 65 5b 36 34 5d 3b 0d 09 |char nam|e[64];..|
|00000b00| 6c 6f 6e 67 20 72 6f 77 | 73 2c 63 6f 6c 73 3b 0d |long row|s,cols;.|
|00000b10| 09 64 6f 75 62 6c 65 20 | 2a 78 72 2c 2a 78 69 3b |.double |*xr,*xi;|
|00000b20| 0d 09 73 68 6f 72 74 20 | 2a 6d 72 2c 2a 6d 69 3b |..short |*mr,*mi;|
|00000b30| 0d 09 69 6e 74 20 65 72 | 72 6f 72 3b 0d 09 66 20 |..int er|ror;..f |
|00000b40| 3d 20 66 6f 70 65 6e 28 | 22 66 6f 6f 2e 6d 61 74 |= fopen(|"foo.mat|
|00000b50| 22 2c 22 72 62 22 29 3b | 0d 09 65 72 72 6f 72 3d |","rb");|..error=|
|00000b60| 4c 6f 61 64 4e 65 78 74 | 4d 61 74 49 6e 74 6f 44 |LoadNext|MatIntoD|
|00000b70| 6f 75 62 6c 65 73 28 66 | 2c 6e 61 6d 65 2c 26 72 |oubles(f|,name,&r|
|00000b80| 6f 77 73 2c 26 63 6f 6c | 73 2c 26 78 72 2c 26 78 |ows,&col|s,&xr,&x|
|00000b90| 69 29 3b 0d 20 09 65 72 | 72 6f 72 3d 4c 6f 61 64 |i);. .er|ror=Load|
|00000ba0| 4e 65 78 74 4d 61 74 49 | 6e 74 6f 53 68 6f 72 74 |NextMatI|ntoShort|
|00000bb0| 73 28 66 2c 6e 61 6d 65 | 2c 26 72 6f 77 73 2c 26 |s(f,name|,&rows,&|
|00000bc0| 63 6f 6c 73 2c 26 6d 72 | 2c 26 6d 69 29 3b 0d 09 |cols,&mr|,&mi);..|
|00000bd0| 66 63 6c 6f 73 65 28 66 | 29 3b 0d 09 66 72 65 65 |fclose(f|);..free|
|00000be0| 28 78 72 29 3b 0d 09 66 | 72 65 65 28 78 69 29 3b |(xr);..f|ree(xi);|
|00000bf0| 0d 09 66 72 65 65 28 6d | 72 29 3b 0d 09 66 72 65 |..free(m|r);..fre|
|00000c00| 65 28 6d 69 29 3b 0d 0d | 41 6c 74 65 72 6e 61 74 |e(mi);..|Alternat|
|00000c10| 69 76 65 6c 79 2c 20 79 | 6f 75 20 63 61 6e 20 72 |ively, y|ou can r|
|00000c20| 65 61 64 20 61 20 73 70 | 65 63 69 66 69 63 20 6d |ead a sp|ecific m|
|00000c30| 61 74 72 69 78 2c 20 62 | 79 20 6e 61 6d 65 2c 20 |atrix, b|y name, |
|00000c40| 66 72 6f 6d 20 74 68 65 | 20 66 69 6c 65 3a 0d 0d |from the| file:..|
|00000c50| 09 65 72 72 6f 72 3d 4c | 6f 61 64 4e 61 6d 65 64 |.error=L|oadNamed|
|00000c60| 4d 61 74 49 6e 74 6f 44 | 6f 75 62 6c 65 73 28 66 |MatIntoD|oubles(f|
|00000c70| 2c 22 50 75 74 20 6e 61 | 6d 65 20 68 65 72 65 22 |,"Put na|me here"|
|00000c80| 2c 26 72 6f 77 73 2c 26 | 63 6f 6c 73 2c 26 78 72 |,&rows,&|cols,&xr|
|00000c90| 2c 26 78 69 29 3b 0d 20 | 09 65 72 72 6f 72 3d 4c |,&xi);. |.error=L|
|00000ca0| 6f 61 64 4e 61 6d 65 64 | 4d 61 74 49 6e 74 6f 53 |oadNamed|MatIntoS|
|00000cb0| 68 6f 72 74 73 28 66 2c | 22 50 75 74 20 6e 61 6d |horts(f,|"Put nam|
|00000cc0| 65 20 68 65 72 65 22 2c | 26 72 6f 77 73 2c 26 63 |e here",|&rows,&c|
|00000cd0| 6f 6c 73 2c 26 6d 72 2c | 26 6d 69 29 3b 0d 0d 54 |ols,&mr,|&mi);..T|
|00000ce0| 68 69 73 20 65 78 61 6d | 70 6c 65 20 63 72 65 61 |his exam|ple crea|
|00000cf0| 74 65 73 20 61 20 4d 41 | 54 4c 41 42 20 66 69 6c |tes a MA|TLAB fil|
|00000d00| 65 20 63 6f 6e 74 61 69 | 6e 69 6e 67 20 66 6f 75 |e contai|ning fou|
|00000d10| 72 20 6d 61 74 72 69 63 | 65 73 3a 0d 0d 09 46 49 |r matric|es:...FI|
|00000d20| 4c 45 20 2a 66 70 3b 0d | 09 64 6f 75 62 6c 65 20 |LE *fp;.|.double |
|00000d30| 78 79 7a 5b 31 30 30 30 | 5d 2c 61 72 5b 31 30 30 |xyz[1000|],ar[100|
|00000d40| 30 5d 2c 61 69 5b 31 30 | 30 30 5d 3b 0d 09 73 68 |0],ai[10|00];..sh|
|00000d50| 6f 72 74 20 6d 5b 31 30 | 30 5d 2c 6e 5b 31 30 30 |ort m[10|0],n[100|
|00000d60| 5d 3b 0d 09 66 70 20 3d | 20 66 6f 70 65 6e 28 22 |];..fp =| fopen("|
|00000d70| 62 61 72 2e 6d 61 74 22 | 2c 22 77 62 22 29 3b 0d |bar.mat"|,"wb");.|
|00000d80| 09 53 65 74 46 69 6c 65 | 49 6e 66 6f 28 22 62 61 |.SetFile|Info("ba|
|00000d90| 72 2e 6d 61 74 22 2c 27 | 4d 41 54 57 27 2c 27 4d |r.mat",'|MATW','M|
|00000da0| 41 54 4c 27 29 3b 0d 09 | 65 72 72 6f 72 3d 53 61 |ATL');..|error=Sa|
|00000db0| 76 65 4d 61 74 44 6f 75 | 62 6c 65 73 28 66 70 2c |veMatDou|bles(fp,|
|00000dc0| 22 78 79 7a 22 2c 32 2c | 33 2c 78 79 7a 2c 4e 55 |"xyz",2,|3,xyz,NU|
|00000dd0| 4c 4c 29 3b 09 2f 2f 20 | 72 65 61 6c 0d 09 65 72 |LL);.// |real..er|
|00000de0| 72 6f 72 3d 53 61 76 65 | 4d 61 74 53 68 6f 72 74 |ror=Save|MatShort|
|00000df0| 73 28 66 70 2c 22 6d 22 | 2c 32 2c 33 2c 6d 2c 4e |s(fp,"m"|,2,3,m,N|
|00000e00| 55 4c 4c 29 3b 09 09 09 | 2f 2f 20 72 65 61 6c 0d |ULL);...|// real.|
|00000e10| 09 65 72 72 6f 72 3d 53 | 61 76 65 4d 61 74 44 6f |.error=S|aveMatDo|
|00000e20| 75 62 6c 65 73 28 66 70 | 2c 22 61 22 2c 35 2c 35 |ubles(fp|,"a",5,5|
|00000e30| 2c 61 72 2c 61 69 29 3b | 09 09 09 2f 2f 20 63 6f |,ar,ai);|...// co|
|00000e40| 6d 70 6c 65 78 0d 09 65 | 72 72 6f 72 3d 53 61 76 |mplex..e|rror=Sav|
|00000e50| 65 4d 61 74 53 68 6f 72 | 74 73 28 66 70 2c 22 6d |eMatShor|ts(fp,"m|
|00000e60| 63 22 2c 32 2c 33 2c 6d | 2c 6e 29 3b 09 09 09 2f |c",2,3,m|,n);.../|
|00000e70| 2f 20 63 6f 6d 70 6c 65 | 78 0d 09 66 63 6c 6f 73 |/ comple|x..fclos|
|00000e80| 65 28 66 70 29 3b 0d 0d | 0d 48 49 53 54 4f 52 59 |e(fp);..|.HISTORY|
|00000e90| 3a 0d 31 31 2d 33 2d 38 | 36 20 4a 2e 4e 2e 20 4c |:.11-3-8|6 J.N. L|
|00000ea0| 69 74 74 6c 65 20 77 72 | 6f 74 65 20 6c 6f 61 64 |ittle wr|ote load|
|00000eb0| 6d 61 74 2e 63 20 61 6e | 64 20 73 61 76 65 6d 61 |mat.c an|d savema|
|00000ec0| 74 2e 63 0d 0d 5b 31 34 | 2d 46 65 62 2d 39 31 5d |t.c..[14|-Feb-91]|
|00000ed0| 20 6a 6d 62 20 2d 2d 20 | 41 64 64 65 64 20 73 75 | jmb -- |Added su|
|00000ee0| 70 70 6f 72 74 20 66 6f | 72 20 4d 50 57 20 43 20 |pport fo|r MPW C |
|00000ef0| 33 2e 78 20 61 6e 64 20 | 54 48 49 4e 4b 20 43 2e |3.x and |THINK C.|
|00000f00| 0d 0d 31 2f 34 2f 39 33 | 09 64 67 70 20 52 65 6e |..1/4/93|.dgp Ren|
|00000f10| 61 6d 65 64 20 74 6f 20 | 4d 61 74 4c 61 62 2e 63 |amed to |MatLab.c|
|00000f20| 2e 20 44 4f 55 42 4c 45 | 20 69 73 20 6e 6f 77 20 |. DOUBLE| is now |
|00000f30| 75 73 65 64 20 73 6f 6c | 65 79 20 77 69 74 68 69 |used sol|ey withi|
|00000f40| 6e 20 74 68 65 20 6e 65 | 77 0d 57 72 69 74 65 58 |n the ne|w.WriteX|
|00000f50| 58 58 20 61 6e 64 20 52 | 65 61 64 58 58 58 20 73 |XX and R|eadXXX s|
|00000f60| 75 62 72 6f 75 74 69 6e | 65 73 2e 20 43 20 75 73 |ubroutin|es. C us|
|00000f70| 65 72 73 20 6e 6f 77 20 | 6f 6e 6c 79 20 64 65 61 |ers now |only dea|
|00000f80| 6c 20 77 69 74 68 20 64 | 6f 75 62 6c 65 20 61 6e |l with d|ouble an|
|00000f90| 64 20 73 68 6f 72 74 0d | 61 72 72 61 79 73 3b 20 |d short.|arrays; |
|00000fa0| 6e 65 63 65 73 73 61 72 | 79 20 63 6f 6e 76 65 72 |necessar|y conver|
|00000fb0| 73 69 6f 6e 73 20 61 72 | 65 20 64 6f 6e 65 20 62 |sions ar|e done b|
|00000fc0| 65 68 69 6e 64 20 74 68 | 65 20 73 63 65 6e 65 73 |ehind th|e scenes|
|00000fd0| 2e 20 64 6f 75 62 6c 65 | 20 69 73 20 66 61 73 74 |. double| is fast|
|00000fe0| 65 72 20 61 6e 64 0d 6d | 6f 72 65 20 63 6f 6e 76 |er and.m|ore conv|
|00000ff0| 65 6e 69 65 6e 74 20 66 | 6f 72 20 4d 61 63 69 6e |enient f|or Macin|
|00001000| 74 6f 73 68 20 70 72 6f | 67 72 61 6d 6d 69 6e 67 |tosh pro|gramming|
|00001010| 20 74 68 61 6e 20 74 68 | 65 20 38 2d 62 79 74 65 | than th|e 8-byte|
|00001020| 20 73 68 6f 72 74 20 64 | 6f 75 62 6c 65 20 75 73 | short d|ouble us|
|00001030| 65 64 20 69 6e 0d 74 68 | 65 20 4d 41 54 4c 41 42 |ed in.th|e MATLAB|
|00001040| 20 66 69 6c 65 73 2e 20 | 41 64 64 65 64 20 73 75 | files. |Added su|
|00001050| 70 70 6f 72 74 20 66 6f | 72 20 66 69 6c 65 20 66 |pport fo|r file f|
|00001060| 6f 72 6d 61 74 20 33 3a | 20 73 69 67 6e 65 64 20 |ormat 3:| signed |
|00001070| 73 68 6f 72 74 20 69 6e | 74 73 2e 20 4f 6d 69 74 |short in|ts. Omit|
|00001080| 74 65 64 0d 74 68 65 20 | 69 6d 61 67 66 20 66 6c |ted.the |imagf fl|
|00001090| 61 67 20 73 69 6e 63 65 | 20 69 74 27 73 20 72 65 |ag since| it's re|
|000010a0| 64 75 6e 64 61 6e 74 2e | 20 4f 6d 69 74 74 65 64 |dundant.| Omitted|
|000010b0| 20 74 68 65 20 74 79 70 | 65 20 61 72 67 75 6d 65 | the typ|e argume|
|000010c0| 6e 74 20 62 79 20 73 75 | 70 70 6c 79 69 6e 67 0d |nt by su|pplying.|
|000010d0| 73 65 76 65 72 61 6c 20 | 75 73 65 72 2d 63 61 6c |several |user-cal|
|000010e0| 6c 61 62 6c 65 20 66 72 | 6f 6e 74 20 65 6e 64 73 |lable fr|ont ends|
|000010f0| 2c 20 65 61 63 68 20 74 | 61 69 6c 6f 72 65 64 20 |, each t|ailored |
|00001100| 74 6f 20 61 20 70 61 72 | 74 69 63 75 6c 61 72 20 |to a par|ticular |
|00001110| 6e 75 6d 62 65 72 20 74 | 79 70 65 2e 0d 0d 32 2f |number t|ype...2/|
|00001120| 39 33 20 64 67 70 20 52 | 65 6e 61 6d 65 64 20 74 |93 dgp R|enamed t|
|00001130| 6f 20 52 65 61 64 4d 61 | 74 4c 61 62 46 69 6c 65 |o ReadMa|tLabFile|
|00001140| 2e 63 0d 0d 36 2f 32 32 | 2d 32 35 2f 39 33 20 6a |.c..6/22|-25/93 j|
|00001150| 61 73 20 26 20 64 67 70 | 20 77 72 6f 74 65 20 4c |as & dgp| wrote L|
|00001160| 6f 61 64 4e 61 6d 65 64 | 4d 61 74 2e 20 52 65 6e |oadNamed|Mat. Ren|
|00001170| 61 6d 65 64 20 4c 6f 61 | 64 4d 61 74 53 68 6f 72 |amed Loa|dMatShor|
|00001180| 74 73 20 61 6e 64 0d 4c | 6f 61 64 4d 61 74 44 6f |ts and.L|oadMatDo|
|00001190| 75 62 6c 65 73 20 74 6f | 20 4c 6f 61 64 4e 65 78 |ubles to| LoadNex|
|000011a0| 74 4d 61 74 49 6e 74 6f | 53 68 6f 72 74 73 20 61 |tMatInto|Shorts a|
|000011b0| 6e 64 20 4c 6f 61 64 4e | 65 78 74 4d 61 74 49 6e |nd LoadN|extMatIn|
|000011c0| 74 6f 44 6f 75 62 6c 65 | 73 2e 20 41 64 64 65 64 |toDouble|s. Added|
|000011d0| 0d 4c 6f 61 64 4e 61 6d | 65 64 4d 61 74 49 6e 74 |.LoadNam|edMatInt|
|000011e0| 6f 53 68 6f 72 74 73 20 | 61 6e 64 20 4c 6f 61 64 |oShorts |and Load|
|000011f0| 4e 61 6d 65 64 4d 61 74 | 49 6e 74 6f 44 6f 75 62 |NamedMat|IntoDoub|
|00001200| 6c 65 73 2e 0d 0d 31 32 | 2f 39 33 20 64 67 70 20 |les...12|/93 dgp |
|00001210| 52 65 6e 61 6d 65 64 20 | 74 6f 20 52 65 61 64 4d |Renamed |to ReadM|
|00001220| 41 54 4c 41 42 46 69 6c | 65 2e 63 0d 0d 37 2f 32 |ATLABFil|e.c..7/2|
|00001230| 38 2f 39 34 20 64 67 70 | 20 61 64 64 65 64 20 73 |8/94 dgp| added s|
|00001240| 75 70 70 6f 72 74 20 66 | 6f 72 20 4d 65 74 72 6f |upport f|or Metro|
|00001250| 77 65 72 6b 73 20 43 6f | 64 65 57 61 72 72 69 6f |werks Co|deWarrio|
|00001260| 72 20 43 20 63 6f 6d 70 | 69 6c 65 72 2e 0d 0d 36 |r C comp|iler...6|
|00001270| 2f 38 2f 39 35 20 64 67 | 70 20 4d 61 64 65 20 73 |/8/95 dg|p Made s|
|00001280| 75 72 65 20 74 68 61 74 | 20 74 68 65 20 4d 41 54 |ure that| the MAT|
|00001290| 4c 41 42 20 73 74 72 75 | 63 74 20 22 46 6d 61 74 |LAB stru|ct "Fmat|
|000012a0| 72 69 78 22 20 69 73 20 | 61 6c 77 61 79 73 20 6e |rix" is |always n|
|000012b0| 61 74 69 76 65 20 61 6c | 69 67 6e 65 64 2e 0d 2a |ative al|igned..*|
|000012c0| 2f 0d 23 69 6e 63 6c 75 | 64 65 20 22 56 69 64 65 |/.#inclu|de "Vide|
|000012d0| 6f 54 6f 6f 6c 62 6f 78 | 2e 68 22 0d 0d 64 6f 75 |oToolbox|.h"..dou|
|000012e0| 62 6c 65 20 2a 52 65 61 | 64 44 6f 75 62 6c 65 73 |ble *Rea|dDoubles|
|000012f0| 28 46 49 4c 45 20 2a 66 | 2c 6c 6f 6e 67 20 65 6c |(FILE *f|,long el|
|00001300| 65 6d 65 6e 74 73 2c 63 | 68 61 72 20 2a 6e 61 6d |ements,c|har *nam|
|00001310| 65 29 3b 0d 73 68 6f 72 | 74 20 2a 52 65 61 64 44 |e);.shor|t *ReadD|
|00001320| 6f 75 62 6c 65 73 49 6e | 74 6f 53 68 6f 72 74 73 |oublesIn|toShorts|
|00001330| 28 46 49 4c 45 20 2a 66 | 2c 6c 6f 6e 67 20 65 6c |(FILE *f|,long el|
|00001340| 65 6d 65 6e 74 73 2c 63 | 68 61 72 20 2a 6e 61 6d |ements,c|har *nam|
|00001350| 65 29 3b 0d 69 6e 74 20 | 57 72 69 74 65 44 6f 75 |e);.int |WriteDou|
|00001360| 62 6c 65 73 28 46 49 4c | 45 20 2a 66 2c 6c 6f 6e |bles(FIL|E *f,lon|
|00001370| 67 20 65 6c 65 6d 65 6e | 74 73 2c 63 68 61 72 20 |g elemen|ts,char |
|00001380| 2a 6e 61 6d 65 2c 64 6f | 75 62 6c 65 20 2a 64 29 |*name,do|uble *d)|
|00001390| 3b 0d 69 6e 74 20 57 72 | 69 74 65 53 68 6f 72 74 |;.int Wr|iteShort|
|000013a0| 73 28 46 49 4c 45 20 2a | 66 2c 6c 6f 6e 67 20 65 |s(FILE *|f,long e|
|000013b0| 6c 65 6d 65 6e 74 73 2c | 63 68 61 72 20 2a 6e 61 |lements,|char *na|
|000013c0| 6d 65 2c 73 68 6f 72 74 | 20 2a 64 29 3b 0d 69 6e |me,short| *d);.in|
|000013d0| 74 20 53 61 76 65 4d 61 | 74 28 46 49 4c 45 20 2a |t SaveMa|t(FILE *|
|000013e0| 66 2c 63 68 61 72 20 2a | 6e 61 6d 65 2c 6c 6f 6e |f,char *|name,lon|
|000013f0| 67 20 72 6f 77 73 2c 6c | 6f 6e 67 20 63 6f 6c 73 |g rows,l|ong cols|
|00001400| 2c 76 6f 69 64 20 2a 72 | 65 61 6c 2c 76 6f 69 64 |,void *r|eal,void|
|00001410| 20 2a 69 6d 61 67 2c 6c | 6f 6e 67 20 74 79 70 65 | *imag,l|ong type|
|00001420| 29 3b 0d 69 6e 74 20 4c | 6f 61 64 4e 61 6d 65 64 |);.int L|oadNamed|
|00001430| 4d 61 74 28 46 49 4c 45 | 20 2a 66 2c 63 6f 6e 73 |Mat(FILE| *f,cons|
|00001440| 74 20 63 68 61 72 20 2a | 6e 61 6d 65 2c 6c 6f 6e |t char *|name,lon|
|00001450| 67 20 2a 72 6f 77 73 2c | 6c 6f 6e 67 20 2a 63 6f |g *rows,|long *co|
|00001460| 6c 73 0d 09 2c 76 6f 69 | 64 20 2a 2a 72 65 61 6c |ls..,voi|d **real|
|00001470| 2c 76 6f 69 64 20 2a 2a | 69 6d 61 67 2c 6c 6f 6e |,void **|imag,lon|
|00001480| 67 20 2a 66 69 6c 65 54 | 79 70 65 2c 6c 6f 6e 67 |g *fileT|ype,long|
|00001490| 20 64 65 73 69 72 65 64 | 4e 75 6d 62 65 72 54 79 | desired|NumberTy|
|000014a0| 70 65 29 3b 0d 69 6e 74 | 20 4c 6f 61 64 4e 65 78 |pe);.int| LoadNex|
|000014b0| 74 4d 61 74 28 46 49 4c | 45 20 2a 66 2c 63 68 61 |tMat(FIL|E *f,cha|
|000014c0| 72 20 2a 6e 61 6d 65 2c | 6c 6f 6e 67 20 2a 72 6f |r *name,|long *ro|
|000014d0| 77 73 2c 6c 6f 6e 67 20 | 2a 63 6f 6c 73 0d 09 2c |ws,long |*cols..,|
|000014e0| 76 6f 69 64 20 2a 2a 72 | 65 61 6c 2c 76 6f 69 64 |void **r|eal,void|
|000014f0| 20 2a 2a 69 6d 61 67 2c | 6c 6f 6e 67 20 2a 66 69 | **imag,|long *fi|
|00001500| 6c 65 54 79 70 65 2c 6c | 6f 6e 67 20 64 65 73 69 |leType,l|ong desi|
|00001510| 72 65 64 4e 75 6d 62 65 | 72 54 79 70 65 29 3b 0d |redNumbe|rType);.|
|00001520| 0d 23 69 66 20 50 52 41 | 47 4d 41 5f 41 4c 49 47 |.#if PRA|GMA_ALIG|
|00001530| 4e 5f 53 55 50 50 4f 52 | 54 45 44 20 7c 7c 20 5f |N_SUPPOR|TED || _|
|00001540| 5f 4d 57 45 52 4b 53 5f | 5f 0d 09 23 70 72 61 67 |_MWERKS_|_..#prag|
|00001550| 6d 61 20 6f 70 74 69 6f | 6e 73 20 61 6c 69 67 6e |ma optio|ns align|
|00001560| 3d 6e 61 74 69 76 65 0d | 23 65 6e 64 69 66 0d 0d |=native.|#endif..|
|00001570| 74 79 70 65 64 65 66 20 | 73 74 72 75 63 74 20 7b |typedef |struct {|
|00001580| 0d 20 20 20 20 20 6c 6f | 6e 67 20 74 79 70 65 3b |.     lo|ng type;|
|00001590| 09 09 2f 2a 20 74 79 70 | 65 20 2a 2f 0d 20 20 20 |../* typ|e */.   |
|000015a0| 20 20 6c 6f 6e 67 20 72 | 6f 77 73 3b 09 09 2f 2a |  long r|ows;../*|
|000015b0| 20 72 6f 77 20 64 69 6d | 65 6e 73 69 6f 6e 20 2a | row dim|ension *|
|000015c0| 2f 0d 20 20 20 20 20 6c | 6f 6e 67 20 63 6f 6c 73 |/.     l|ong cols|
|000015d0| 3b 09 09 2f 2a 20 63 6f | 6c 75 6d 6e 20 64 69 6d |;../* co|lumn dim|
|000015e0| 65 6e 73 69 6f 6e 20 2a | 2f 0d 20 20 20 20 20 6c |ension *|/.     l|
|000015f0| 6f 6e 67 20 69 6d 61 67 | 66 3b 09 2f 2a 20 66 6c |ong imag|f;./* fl|
|00001600| 61 67 20 69 6e 64 69 63 | 61 74 69 6e 67 20 69 6d |ag indic|ating im|
|00001610| 61 67 20 70 61 72 74 20 | 2a 2f 0d 20 20 20 20 20 |ag part |*/.     |
|00001620| 6c 6f 6e 67 20 6e 61 6d | 6c 65 6e 3b 09 2f 2a 20 |long nam|len;./* |
|00001630| 6e 61 6d 65 20 6c 65 6e | 67 74 68 20 28 69 6e 63 |name len|gth (inc|
|00001640| 6c 75 64 69 6e 67 20 4e | 55 4c 4c 29 20 2a 2f 0d |luding N|ULL) */.|
|00001650| 7d 20 46 6d 61 74 72 69 | 78 3b 0d 0d 23 69 66 20 |} Fmatri|x;..#if |
|00001660| 50 52 41 47 4d 41 5f 41 | 4c 49 47 4e 5f 53 55 50 |PRAGMA_A|LIGN_SUP|
|00001670| 50 4f 52 54 45 44 20 7c | 7c 20 5f 5f 4d 57 45 52 |PORTED ||| __MWER|
|00001680| 4b 53 5f 5f 0d 09 23 70 | 72 61 67 6d 61 20 6f 70 |KS__..#p|ragma op|
|00001690| 74 69 6f 6e 73 20 61 6c | 69 67 6e 3d 72 65 73 65 |tions al|ign=rese|
|000016a0| 74 0d 23 65 6e 64 69 66 | 0d 0d 2f 2a 20 46 72 6f |t.#endif|../* Fro|
|000016b0| 6d 20 63 6d 65 78 2e 68 | 20 2a 2f 0d 2f 2a 0d 4d |m cmex.h| */./*.M|
|000016c0| 41 54 4c 41 42 20 66 69 | 6c 65 73 20 77 69 74 68 |ATLAB fi|les with|
|000016d0| 20 6e 75 6d 62 65 72 20 | 66 6f 72 6d 61 74 20 30 | number |format 0|
|000016e0| 20 75 73 65 20 38 2d 62 | 79 74 65 20 66 6c 6f 61 | use 8-b|yte floa|
|000016f0| 74 69 6e 67 20 70 6f 69 | 6e 74 20 6e 75 6d 62 65 |ting poi|nt numbe|
|00001700| 72 73 2c 20 77 68 69 63 | 68 0d 77 65 27 6c 6c 20 |rs, whic|h.we'll |
|00001710| 63 61 6c 6c 20 44 4f 55 | 42 4c 45 2e 20 49 74 27 |call DOU|BLE. It'|
|00001720| 73 20 64 65 66 69 6e 65 | 64 20 69 6e 20 56 69 64 |s define|d in Vid|
|00001730| 65 6f 54 6f 6f 6c 62 6f | 78 2e 68 0d 2a 2f 0d 0d |eoToolbo|x.h.*/..|
|00001740| 69 6e 74 20 4c 6f 61 64 | 4e 65 78 74 4d 61 74 49 |int Load|NextMatI|
|00001750| 6e 74 6f 44 6f 75 62 6c | 65 73 28 46 49 4c 45 20 |ntoDoubl|es(FILE |
|00001760| 2a 66 2c 63 68 61 72 20 | 2a 6e 61 6d 65 2c 6c 6f |*f,char |*name,lo|
|00001770| 6e 67 20 2a 72 6f 77 73 | 2c 6c 6f 6e 67 20 2a 63 |ng *rows|,long *c|
|00001780| 6f 6c 73 0d 09 2c 64 6f | 75 62 6c 65 20 2a 2a 72 |ols..,do|uble **r|
|00001790| 65 61 6c 2c 64 6f 75 62 | 6c 65 20 2a 2a 69 6d 61 |eal,doub|le **ima|
|000017a0| 67 29 0d 7b 0d 09 6c 6f | 6e 67 20 66 69 6c 65 54 |g).{..lo|ng fileT|
|000017b0| 79 70 65 2c 64 65 73 69 | 72 65 64 4e 75 6d 62 65 |ype,desi|redNumbe|
|000017c0| 72 54 79 70 65 3d 30 3b | 0d 09 0d 09 72 65 74 75 |rType=0;|....retu|
|000017d0| 72 6e 20 4c 6f 61 64 4e | 65 78 74 4d 61 74 28 66 |rn LoadN|extMat(f|
|000017e0| 2c 6e 61 6d 65 2c 72 6f | 77 73 2c 63 6f 6c 73 2c |,name,ro|ws,cols,|
|000017f0| 28 76 6f 69 64 20 2a 2a | 29 72 65 61 6c 2c 28 76 |(void **|)real,(v|
|00001800| 6f 69 64 20 2a 2a 29 69 | 6d 61 67 0d 09 09 2c 26 |oid **)i|mag...,&|
|00001810| 66 69 6c 65 54 79 70 65 | 2c 64 65 73 69 72 65 64 |fileType|,desired|
|00001820| 4e 75 6d 62 65 72 54 79 | 70 65 29 3b 0d 7d 0d 09 |NumberTy|pe);.}..|
|00001830| 0d 69 6e 74 20 4c 6f 61 | 64 4e 65 78 74 4d 61 74 |.int Loa|dNextMat|
|00001840| 49 6e 74 6f 53 68 6f 72 | 74 73 28 46 49 4c 45 20 |IntoShor|ts(FILE |
|00001850| 2a 66 2c 63 68 61 72 20 | 2a 6e 61 6d 65 2c 6c 6f |*f,char |*name,lo|
|00001860| 6e 67 20 2a 72 6f 77 73 | 2c 6c 6f 6e 67 20 2a 63 |ng *rows|,long *c|
|00001870| 6f 6c 73 0d 09 2c 73 68 | 6f 72 74 20 2a 2a 72 65 |ols..,sh|ort **re|
|00001880| 61 6c 2c 73 68 6f 72 74 | 20 2a 2a 69 6d 61 67 29 |al,short| **imag)|
|00001890| 0d 7b 0d 09 6c 6f 6e 67 | 20 66 69 6c 65 54 79 70 |.{..long| fileTyp|
|000018a0| 65 2c 64 65 73 69 72 65 | 64 4e 75 6d 62 65 72 54 |e,desire|dNumberT|
|000018b0| 79 70 65 3d 33 3b 0d 09 | 0d 09 72 65 74 75 72 6e |ype=3;..|..return|
|000018c0| 20 4c 6f 61 64 4e 65 78 | 74 4d 61 74 28 66 2c 6e | LoadNex|tMat(f,n|
|000018d0| 61 6d 65 2c 72 6f 77 73 | 2c 63 6f 6c 73 2c 28 76 |ame,rows|,cols,(v|
|000018e0| 6f 69 64 20 2a 2a 29 72 | 65 61 6c 2c 28 76 6f 69 |oid **)r|eal,(voi|
|000018f0| 64 20 2a 2a 29 69 6d 61 | 67 0d 09 09 2c 26 66 69 |d **)ima|g...,&fi|
|00001900| 6c 65 54 79 70 65 2c 64 | 65 73 69 72 65 64 4e 75 |leType,d|esiredNu|
|00001910| 6d 62 65 72 54 79 70 65 | 29 3b 0d 7d 0d 0d 69 6e |mberType|);.}..in|
|00001920| 74 20 4c 6f 61 64 4e 61 | 6d 65 64 4d 61 74 49 6e |t LoadNa|medMatIn|
|00001930| 74 6f 44 6f 75 62 6c 65 | 73 28 46 49 4c 45 20 2a |toDouble|s(FILE *|
|00001940| 66 2c 63 6f 6e 73 74 20 | 63 68 61 72 20 2a 6e 61 |f,const |char *na|
|00001950| 6d 65 2c 6c 6f 6e 67 20 | 2a 72 6f 77 73 2c 6c 6f |me,long |*rows,lo|
|00001960| 6e 67 20 2a 63 6f 6c 73 | 0d 09 2c 64 6f 75 62 6c |ng *cols|..,doubl|
|00001970| 65 20 2a 2a 72 65 61 6c | 2c 64 6f 75 62 6c 65 20 |e **real|,double |
|00001980| 2a 2a 69 6d 61 67 29 0d | 7b 0d 09 6c 6f 6e 67 20 |**imag).|{..long |
|00001990| 66 69 6c 65 54 79 70 65 | 2c 64 65 73 69 72 65 64 |fileType|,desired|
|000019a0| 4e 75 6d 62 65 72 54 79 | 70 65 3d 30 3b 0d 09 0d |NumberTy|pe=0;...|
|000019b0| 09 72 65 74 75 72 6e 20 | 4c 6f 61 64 4e 61 6d 65 |.return |LoadName|
|000019c0| 64 4d 61 74 28 66 2c 6e | 61 6d 65 2c 72 6f 77 73 |dMat(f,n|ame,rows|
|000019d0| 2c 63 6f 6c 73 2c 28 76 | 6f 69 64 20 2a 2a 29 72 |,cols,(v|oid **)r|
|000019e0| 65 61 6c 2c 28 76 6f 69 | 64 20 2a 2a 29 69 6d 61 |eal,(voi|d **)ima|
|000019f0| 67 0d 09 09 2c 26 66 69 | 6c 65 54 79 70 65 2c 64 |g...,&fi|leType,d|
|00001a00| 65 73 69 72 65 64 4e 75 | 6d 62 65 72 54 79 70 65 |esiredNu|mberType|
|00001a10| 29 3b 0d 7d 0d 09 0d 69 | 6e 74 20 4c 6f 61 64 4e |);.}...i|nt LoadN|
|00001a20| 61 6d 65 64 4d 61 74 49 | 6e 74 6f 53 68 6f 72 74 |amedMatI|ntoShort|
|00001a30| 73 28 46 49 4c 45 20 2a | 66 2c 63 6f 6e 73 74 20 |s(FILE *|f,const |
|00001a40| 63 68 61 72 20 2a 6e 61 | 6d 65 2c 6c 6f 6e 67 20 |char *na|me,long |
|00001a50| 2a 72 6f 77 73 2c 6c 6f | 6e 67 20 2a 63 6f 6c 73 |*rows,lo|ng *cols|
|00001a60| 0d 09 2c 73 68 6f 72 74 | 20 2a 2a 72 65 61 6c 2c |..,short| **real,|
|00001a70| 73 68 6f 72 74 20 2a 2a | 69 6d 61 67 29 0d 7b 0d |short **|imag).{.|
|00001a80| 09 6c 6f 6e 67 20 66 69 | 6c 65 54 79 70 65 2c 64 |.long fi|leType,d|
|00001a90| 65 73 69 72 65 64 4e 75 | 6d 62 65 72 54 79 70 65 |esiredNu|mberType|
|00001aa0| 3d 33 3b 0d 09 0d 09 72 | 65 74 75 72 6e 20 4c 6f |=3;....r|eturn Lo|
|00001ab0| 61 64 4e 61 6d 65 64 4d | 61 74 28 66 2c 6e 61 6d |adNamedM|at(f,nam|
|00001ac0| 65 2c 72 6f 77 73 2c 63 | 6f 6c 73 2c 28 76 6f 69 |e,rows,c|ols,(voi|
|00001ad0| 64 20 2a 2a 29 72 65 61 | 6c 2c 28 76 6f 69 64 20 |d **)rea|l,(void |
|00001ae0| 2a 2a 29 69 6d 61 67 0d | 09 09 2c 26 66 69 6c 65 |**)imag.|..,&file|
|00001af0| 54 79 70 65 2c 64 65 73 | 69 72 65 64 4e 75 6d 62 |Type,des|iredNumb|
|00001b00| 65 72 54 79 70 65 29 3b | 0d 7d 0d 0d 69 6e 74 20 |erType);|.}..int |
|00001b10| 53 61 76 65 4d 61 74 44 | 6f 75 62 6c 65 73 28 46 |SaveMatD|oubles(F|
|00001b20| 49 4c 45 20 2a 66 2c 63 | 68 61 72 20 2a 6e 61 6d |ILE *f,c|har *nam|
|00001b30| 65 2c 6c 6f 6e 67 20 72 | 6f 77 73 2c 6c 6f 6e 67 |e,long r|ows,long|
|00001b40| 20 63 6f 6c 73 2c 64 6f | 75 62 6c 65 20 2a 72 65 | cols,do|uble *re|
|00001b50| 61 6c 2c 64 6f 75 62 6c | 65 20 2a 69 6d 61 67 29 |al,doubl|e *imag)|
|00001b60| 0d 7b 0d 09 72 65 74 75 | 72 6e 20 53 61 76 65 4d |.{..retu|rn SaveM|
|00001b70| 61 74 28 66 2c 6e 61 6d | 65 2c 72 6f 77 73 2c 63 |at(f,nam|e,rows,c|
|00001b80| 6f 6c 73 2c 72 65 61 6c | 2c 69 6d 61 67 2c 30 29 |ols,real|,imag,0)|
|00001b90| 3b 0d 7d 0d 0d 69 6e 74 | 20 53 61 76 65 4d 61 74 |;.}..int| SaveMat|
|00001ba0| 53 68 6f 72 74 73 28 46 | 49 4c 45 20 2a 66 2c 63 |Shorts(F|ILE *f,c|
|00001bb0| 68 61 72 20 2a 6e 61 6d | 65 2c 6c 6f 6e 67 20 72 |har *nam|e,long r|
|00001bc0| 6f 77 73 2c 6c 6f 6e 67 | 20 63 6f 6c 73 2c 73 68 |ows,long| cols,sh|
|00001bd0| 6f 72 74 20 2a 72 65 61 | 6c 2c 73 68 6f 72 74 20 |ort *rea|l,short |
|00001be0| 2a 69 6d 61 67 29 0d 7b | 0d 09 72 65 74 75 72 6e |*imag).{|..return|
|00001bf0| 20 53 61 76 65 4d 61 74 | 28 66 2c 6e 61 6d 65 2c | SaveMat|(f,name,|
|00001c00| 72 6f 77 73 2c 63 6f 6c | 73 2c 72 65 61 6c 2c 69 |rows,col|s,real,i|
|00001c10| 6d 61 67 2c 33 30 29 3b | 0d 7d 0d 09 0d 69 6e 74 |mag,30);|.}...int|
|00001c20| 20 4c 6f 61 64 4e 61 6d | 65 64 4d 61 74 28 46 49 | LoadNam|edMat(FI|
|00001c30| 4c 45 20 2a 66 2c 63 6f | 6e 73 74 20 63 68 61 72 |LE *f,co|nst char|
|00001c40| 20 2a 6e 61 6d 65 2c 6c | 6f 6e 67 20 2a 72 6f 77 | *name,l|ong *row|
|00001c50| 73 2c 6c 6f 6e 67 20 2a | 63 6f 6c 73 0d 09 2c 76 |s,long *|cols..,v|
|00001c60| 6f 69 64 20 2a 2a 72 65 | 61 6c 2c 76 6f 69 64 20 |oid **re|al,void |
|00001c70| 2a 2a 69 6d 61 67 2c 6c | 6f 6e 67 20 2a 66 69 6c |**imag,l|ong *fil|
|00001c80| 65 54 79 70 65 2c 6c 6f | 6e 67 20 64 65 73 69 72 |eType,lo|ng desir|
|00001c90| 65 64 4e 75 6d 62 65 72 | 54 79 70 65 29 0d 7b 0d |edNumber|Type).{.|
|00001ca0| 09 63 68 61 72 20 63 75 | 72 6e 61 6d 65 5b 31 30 |.char cu|rname[10|
|00001cb0| 30 5d 3b 0d 09 69 6e 74 | 20 65 72 72 6f 72 3b 0d |0];..int| error;.|
|00001cc0| 09 0d 09 72 65 77 69 6e | 64 28 66 29 3b 0d 09 77 |...rewin|d(f);..w|
|00001cd0| 68 69 6c 65 28 31 29 7b | 0d 09 09 65 72 72 6f 72 |hile(1){|...error|
|00001ce0| 3d 4c 6f 61 64 4e 65 78 | 74 4d 61 74 28 66 2c 63 |=LoadNex|tMat(f,c|
|00001cf0| 75 72 6e 61 6d 65 2c 72 | 6f 77 73 2c 63 6f 6c 73 |urname,r|ows,cols|
|00001d00| 2c 72 65 61 6c 2c 69 6d | 61 67 2c 66 69 6c 65 54 |,real,im|ag,fileT|
|00001d10| 79 70 65 2c 64 65 73 69 | 72 65 64 4e 75 6d 62 65 |ype,desi|redNumbe|
|00001d20| 72 54 79 70 65 29 3b 0d | 09 09 69 66 28 65 72 72 |rType);.|..if(err|
|00001d30| 6f 72 29 72 65 74 75 72 | 6e 20 65 72 72 6f 72 3b |or)retur|n error;|
|00001d40| 0d 09 09 69 66 28 73 74 | 72 63 6d 70 28 63 75 72 |...if(st|rcmp(cur|
|00001d50| 6e 61 6d 65 2c 6e 61 6d | 65 29 3d 3d 30 29 62 72 |name,nam|e)==0)br|
|00001d60| 65 61 6b 3b 0d 09 09 69 | 66 28 2a 72 65 61 6c 21 |eak;...i|f(*real!|
|00001d70| 3d 4e 55 4c 4c 29 7b 0d | 09 09 09 66 72 65 65 28 |=NULL){.|...free(|
|00001d80| 2a 72 65 61 6c 29 3b 0d | 09 09 09 2a 72 65 61 6c |*real);.|...*real|
|00001d90| 3d 4e 55 4c 4c 3b 0d 09 | 09 7d 0d 09 09 69 66 28 |=NULL;..|.}...if(|
|00001da0| 2a 69 6d 61 67 21 3d 4e | 55 4c 4c 29 7b 0d 09 09 |*imag!=N|ULL){...|
|00001db0| 09 66 72 65 65 28 2a 69 | 6d 61 67 29 3b 0d 09 09 |.free(*i|mag);...|
|00001dc0| 09 2a 69 6d 61 67 3d 4e | 55 4c 4c 3b 0d 09 09 7d |.*imag=N|ULL;...}|
|00001dd0| 0d 09 7d 0d 09 72 65 74 | 75 72 6e 20 30 3b 0d 7d |..}..ret|urn 0;.}|
|00001de0| 0d 0d 69 6e 74 20 4c 6f | 61 64 4e 65 78 74 4d 61 |..int Lo|adNextMa|
|00001df0| 74 28 46 49 4c 45 20 2a | 66 2c 63 68 61 72 20 2a |t(FILE *|f,char *|
|00001e00| 6e 61 6d 65 2c 6c 6f 6e | 67 20 2a 72 6f 77 73 2c |name,lon|g *rows,|
|00001e10| 6c 6f 6e 67 20 2a 63 6f | 6c 73 0d 09 2c 76 6f 69 |long *co|ls..,voi|
|00001e20| 64 20 2a 2a 72 65 61 6c | 2c 76 6f 69 64 20 2a 2a |d **real|,void **|
|00001e30| 69 6d 61 67 2c 6c 6f 6e | 67 20 2a 66 69 6c 65 54 |imag,lon|g *fileT|
|00001e40| 79 70 65 2c 6c 6f 6e 67 | 20 64 65 73 69 72 65 64 |ype,long| desired|
|00001e50| 4e 75 6d 62 65 72 54 79 | 70 65 29 0d 7b 0d 09 46 |NumberTy|pe).{..F|
|00001e60| 6d 61 74 72 69 78 09 78 | 3b 0d 09 6c 6f 6e 67 20 |matrix.x|;..long |
|00001e70| 65 6c 65 6d 65 6e 74 73 | 2c 6e 61 6d 6c 65 6e 2c |elements|,namlen,|
|00001e80| 6e 75 6d 62 65 72 46 6f | 72 6d 61 74 3b 0d 09 0d |numberFo|rmat;...|
|00001e90| 09 2a 72 65 61 6c 3d 2a | 69 6d 61 67 3d 4e 55 4c |.*real=*|imag=NUL|
|00001ea0| 4c 3b 0d 09 6e 61 6d 65 | 5b 30 5d 3d 30 3b 0d 09 |L;..name|[0]=0;..|
|00001eb0| 69 66 28 66 65 6f 66 28 | 66 29 29 72 65 74 75 72 |if(feof(|f))retur|
|00001ec0| 6e 20 31 3b 0d 0d 09 2f | 2a 20 47 65 74 20 46 6d |n 1;.../|* Get Fm|
|00001ed0| 61 74 72 69 78 20 73 74 | 72 75 63 74 75 72 65 20 |atrix st|ructure |
|00001ee0| 66 72 6f 6d 20 66 69 6c | 65 20 2a 2f 0d 09 69 66 |from fil|e */..if|
|00001ef0| 20 28 66 72 65 61 64 28 | 28 63 68 61 72 20 2a 29 | (fread(|(char *)|
|00001f00| 26 78 2c 73 69 7a 65 6f | 66 28 46 6d 61 74 72 69 |&x,sizeo|f(Fmatri|
|00001f10| 78 29 2c 31 2c 66 29 20 | 21 3d 20 31 29 20 72 65 |x),1,f) |!= 1) re|
|00001f20| 74 75 72 6e 20 32 3b 0d | 09 2a 66 69 6c 65 54 79 |turn 2;.|.*fileTy|
|00001f30| 70 65 20 3d 20 78 2e 74 | 79 70 65 3b 0d 09 2a 72 |pe = x.t|ype;..*r|
|00001f40| 6f 77 73 20 3d 20 78 2e | 72 6f 77 73 3b 0d 09 2a |ows = x.|rows;..*|
|00001f50| 63 6f 6c 73 20 3d 20 78 | 2e 63 6f 6c 73 3b 0d 09 |cols = x|.cols;..|
|00001f60| 6e 61 6d 6c 65 6e 20 3d | 20 78 2e 6e 61 6d 6c 65 |namlen =| x.namle|
|00001f70| 6e 3b 0d 09 65 6c 65 6d | 65 6e 74 73 20 3d 20 78 |n;..elem|ents = x|
|00001f80| 2e 72 6f 77 73 20 2a 20 | 78 2e 63 6f 6c 73 3b 0d |.rows * |x.cols;.|
|00001f90| 0d 09 69 66 28 78 2e 74 | 79 70 65 2f 31 30 30 30 |..if(x.t|ype/1000|
|00001fa0| 21 3d 31 29 72 65 74 75 | 72 6e 20 33 3b 09 09 2f |!=1)retu|rn 3;../|
|00001fb0| 2a 20 6e 6f 74 20 4d 61 | 63 69 6e 74 6f 73 68 2d |* not Ma|cintosh-|
|00001fc0| 63 6f 6d 70 61 74 69 62 | 6c 65 20 62 69 6e 61 72 |compatib|le binar|
|00001fd0| 79 20 66 6f 72 6d 61 74 | 20 2a 2f 0d 09 69 66 28 |y format| */..if(|
|00001fe0| 78 2e 74 79 70 65 2f 31 | 30 30 25 31 30 21 3d 30 |x.type/1|00%10!=0|
|00001ff0| 29 72 65 74 75 72 6e 20 | 34 3b 09 2f 2a 20 74 72 |)return |4;./* tr|
|00002000| 61 6e 73 70 6f 73 65 64 | 20 2a 2f 0d 09 6e 75 6d |ansposed| */..num|
|00002010| 62 65 72 46 6f 72 6d 61 | 74 3d 78 2e 74 79 70 65 |berForma|t=x.type|
|00002020| 2f 31 30 25 31 30 3b 0d | 09 69 66 28 6e 75 6d 62 |/10%10;.|.if(numb|
|00002030| 65 72 46 6f 72 6d 61 74 | 21 3d 30 29 72 65 74 75 |erFormat|!=0)retu|
|00002040| 72 6e 20 35 3b 09 2f 2a | 20 6e 6f 74 20 38 2d 62 |rn 5;./*| not 8-b|
|00002050| 79 74 65 20 64 6f 75 62 | 6c 65 73 20 2a 2f 0d 0d |yte doub|les */..|
|00002060| 09 2f 2a 20 47 65 74 20 | 6d 61 74 72 69 78 20 6e |./* Get |matrix n|
|00002070| 61 6d 65 20 66 72 6f 6d | 20 66 69 6c 65 20 2a 2f |ame from| file */|
|00002080| 0d 09 69 66 20 28 66 72 | 65 61 64 28 6e 61 6d 65 |..if (fr|ead(name|
|00002090| 2c 73 69 7a 65 6f 66 28 | 63 68 61 72 29 2c 6e 61 |,sizeof(|char),na|
|000020a0| 6d 6c 65 6e 2c 66 29 20 | 21 3d 20 6e 61 6d 6c 65 |mlen,f) |!= namle|
|000020b0| 6e 29 20 72 65 74 75 72 | 6e 20 36 3b 0d 09 0d 09 |n) retur|n 6;....|
|000020c0| 2f 2a 20 47 65 74 20 52 | 65 61 6c 20 70 61 72 74 |/* Get R|eal part|
|000020d0| 20 6f 66 20 6d 61 74 72 | 69 78 20 66 72 6f 6d 20 | of matr|ix from |
|000020e0| 66 69 6c 65 20 2a 2f 0d | 09 73 77 69 74 63 68 28 |file */.|.switch(|
|000020f0| 64 65 73 69 72 65 64 4e | 75 6d 62 65 72 54 79 70 |desiredN|umberTyp|
|00002100| 65 29 7b 0d 09 63 61 73 | 65 20 30 3a 0d 09 09 2a |e){..cas|e 0:...*|
|00002110| 72 65 61 6c 3d 52 65 61 | 64 44 6f 75 62 6c 65 73 |real=Rea|dDoubles|
|00002120| 28 66 2c 65 6c 65 6d 65 | 6e 74 73 2c 6e 61 6d 65 |(f,eleme|nts,name|
|00002130| 29 3b 0d 09 09 62 72 65 | 61 6b 3b 0d 09 63 61 73 |);...bre|ak;..cas|
|00002140| 65 20 33 3a 0d 09 09 2a | 72 65 61 6c 3d 52 65 61 |e 3:...*|real=Rea|
|00002150| 64 44 6f 75 62 6c 65 73 | 49 6e 74 6f 53 68 6f 72 |dDoubles|IntoShor|
|00002160| 74 73 28 66 2c 65 6c 65 | 6d 65 6e 74 73 2c 6e 61 |ts(f,ele|ments,na|
|00002170| 6d 65 29 3b 0d 09 09 62 | 72 65 61 6b 3b 0d 09 64 |me);...b|reak;..d|
|00002180| 65 66 61 75 6c 74 3a 0d | 09 09 2a 72 65 61 6c 3d |efault:.|..*real=|
|00002190| 4e 55 4c 4c 3b 0d 09 7d | 0d 09 69 66 28 2a 72 65 |NULL;..}|..if(*re|
|000021a0| 61 6c 3d 3d 4e 55 4c 4c | 29 72 65 74 75 72 6e 20 |al==NULL|)return |
|000021b0| 37 3b 0d 0d 09 2f 2a 20 | 47 65 74 20 49 6d 61 67 |7;.../* |Get Imag|
|000021c0| 20 70 61 72 74 20 6f 66 | 20 6d 61 74 72 69 78 20 | part of| matrix |
|000021d0| 66 72 6f 6d 20 66 69 6c | 65 2c 20 69 66 20 69 74 |from fil|e, if it|
|000021e0| 20 65 78 69 73 74 73 20 | 2a 2f 0d 09 69 66 20 28 | exists |*/..if (|
|000021f0| 78 2e 69 6d 61 67 66 29 | 20 7b 0d 09 09 73 77 69 |x.imagf)| {...swi|
|00002200| 74 63 68 28 64 65 73 69 | 72 65 64 4e 75 6d 62 65 |tch(desi|redNumbe|
|00002210| 72 54 79 70 65 29 7b 0d | 09 09 63 61 73 65 20 30 |rType){.|..case 0|
|00002220| 3a 0d 09 09 09 2a 69 6d | 61 67 3d 52 65 61 64 44 |:....*im|ag=ReadD|
|00002230| 6f 75 62 6c 65 73 28 66 | 2c 65 6c 65 6d 65 6e 74 |oubles(f|,element|
|00002240| 73 2c 6e 61 6d 65 29 3b | 0d 09 09 09 62 72 65 61 |s,name);|....brea|
|00002250| 6b 3b 0d 09 09 63 61 73 | 65 20 33 3a 0d 09 09 09 |k;...cas|e 3:....|
|00002260| 2a 69 6d 61 67 3d 52 65 | 61 64 44 6f 75 62 6c 65 |*imag=Re|adDouble|
|00002270| 73 49 6e 74 6f 53 68 6f | 72 74 73 28 66 2c 65 6c |sIntoSho|rts(f,el|
|00002280| 65 6d 65 6e 74 73 2c 6e | 61 6d 65 29 3b 0d 09 09 |ements,n|ame);...|
|00002290| 09 62 72 65 61 6b 3b 0d | 09 09 7d 0d 09 09 69 66 |.break;.|..}...if|
|000022a0| 28 2a 69 6d 61 67 3d 3d | 4e 55 4c 4c 29 7b 0d 09 |(*imag==|NULL){..|
|000022b0| 09 09 66 72 65 65 28 2a | 72 65 61 6c 29 3b 0d 09 |..free(*|real);..|
|000022c0| 09 09 2a 72 65 61 6c 3d | 4e 55 4c 4c 3b 0d 09 09 |..*real=|NULL;...|
|000022d0| 09 72 65 74 75 72 6e 20 | 38 3b 0d 09 09 7d 0d 09 |.return |8;...}..|
|000022e0| 7d 0d 09 72 65 74 75 72 | 6e 20 30 3b 0d 7d 0d 0d |}..retur|n 0;.}..|
|000022f0| 69 6e 74 20 53 61 76 65 | 4d 61 74 28 46 49 4c 45 |int Save|Mat(FILE|
|00002300| 20 2a 66 2c 63 68 61 72 | 20 2a 6e 61 6d 65 2c 6c | *f,char| *name,l|
|00002310| 6f 6e 67 20 72 6f 77 73 | 2c 6c 6f 6e 67 20 63 6f |ong rows|,long co|
|00002320| 6c 73 2c 76 6f 69 64 20 | 2a 72 65 61 6c 2c 76 6f |ls,void |*real,vo|
|00002330| 69 64 20 2a 69 6d 61 67 | 2c 6c 6f 6e 67 20 74 79 |id *imag|,long ty|
|00002340| 70 65 29 0d 7b 0d 09 46 | 6d 61 74 72 69 78 09 78 |pe).{..F|matrix.x|
|00002350| 3b 0d 09 6c 6f 6e 67 20 | 65 6c 65 6d 65 6e 74 73 |;..long |elements|
|00002360| 3b 0d 09 69 6e 74 20 65 | 72 72 6f 72 3b 0d 09 69 |;..int e|rror;..i|
|00002370| 6e 74 20 66 69 6c 65 46 | 6f 72 6d 61 74 3b 0d 09 |nt fileF|ormat;..|
|00002380| 0d 09 74 79 70 65 25 3d | 31 30 30 3b 09 2f 2a 20 |..type%=|100;./* |
|00002390| 6e 6f 74 20 74 72 61 6e | 73 70 6f 73 65 64 20 2a |not tran|sposed *|
|000023a0| 2f 0d 09 74 79 70 65 2b | 3d 31 30 30 30 3b 09 2f |/..type+|=1000;./|
|000023b0| 2a 20 69 6e 64 69 63 61 | 74 65 20 4d 61 63 69 6e |* indica|te Macin|
|000023c0| 74 6f 73 68 2d 63 6f 6d | 70 61 74 69 62 6c 65 20 |tosh-com|patible |
|000023d0| 62 69 6e 61 72 79 20 66 | 6f 72 6d 61 74 20 2a 2f |binary f|ormat */|
|000023e0| 0d 09 78 2e 74 79 70 65 | 20 3d 20 74 79 70 65 3b |..x.type| = type;|
|000023f0| 0d 09 78 2e 72 6f 77 73 | 20 3d 20 72 6f 77 73 3b |..x.rows| = rows;|
|00002400| 0d 09 78 2e 63 6f 6c 73 | 20 3d 20 63 6f 6c 73 3b |..x.cols| = cols;|
|00002410| 0d 09 69 66 28 69 6d 61 | 67 3d 3d 4e 55 4c 4c 29 |..if(ima|g==NULL)|
|00002420| 78 2e 69 6d 61 67 66 3d | 30 3b 0d 09 65 6c 73 65 |x.imagf=|0;..else|
|00002430| 20 78 2e 69 6d 61 67 66 | 3d 31 3b 0d 09 78 2e 6e | x.imagf|=1;..x.n|
|00002440| 61 6d 6c 65 6e 20 3d 20 | 73 74 72 6c 65 6e 28 6e |amlen = |strlen(n|
|00002450| 61 6d 65 29 20 2b 20 31 | 3b 0d 09 65 6c 65 6d 65 |ame) + 1|;..eleme|
|00002460| 6e 74 73 20 3d 20 78 2e | 72 6f 77 73 20 2a 20 78 |nts = x.|rows * x|
|00002470| 2e 63 6f 6c 73 3b 0d 0d | 09 69 66 28 78 2e 74 79 |.cols;..|.if(x.ty|
|00002480| 70 65 2f 31 30 30 30 21 | 3d 31 29 72 65 74 75 72 |pe/1000!|=1)retur|
|00002490| 6e 20 39 3b 09 09 2f 2a | 20 6e 6f 74 20 4d 61 63 |n 9;../*| not Mac|
|000024a0| 69 6e 74 6f 73 68 2d 63 | 6f 6d 70 61 74 69 62 6c |intosh-c|ompatibl|
|000024b0| 65 20 62 69 6e 61 72 79 | 20 66 6f 72 6d 61 74 20 |e binary| format |
|000024c0| 2a 2f 0d 09 69 66 28 78 | 2e 74 79 70 65 2f 31 30 |*/..if(x|.type/10|
|000024d0| 30 25 31 30 21 3d 30 29 | 72 65 74 75 72 6e 20 31 |0%10!=0)|return 1|
|000024e0| 30 3b 09 2f 2a 20 74 72 | 61 6e 73 70 6f 73 65 64 |0;./* tr|ansposed|
|000024f0| 20 2a 2f 0d 09 66 69 6c | 65 46 6f 72 6d 61 74 3d | */..fil|eFormat=|
|00002500| 78 2e 74 79 70 65 2f 31 | 30 25 31 30 3b 0d 09 69 |x.type/1|0%10;..i|
|00002510| 66 28 66 69 6c 65 46 6f | 72 6d 61 74 21 3d 30 20 |f(fileFo|rmat!=0 |
|00002520| 26 26 20 66 69 6c 65 46 | 6f 72 6d 61 74 21 3d 33 |&& fileF|ormat!=3|
|00002530| 29 72 65 74 75 72 6e 20 | 31 31 3b 20 2f 2a 20 6e |)return |11; /* n|
|00002540| 65 69 74 68 65 72 20 38 | 2d 62 79 74 65 20 64 6f |either 8|-byte do|
|00002550| 75 62 6c 65 2c 20 6e 6f | 72 20 73 68 6f 72 74 20 |uble, no|r short |
|00002560| 2a 2f 0d 0d 09 66 77 72 | 69 74 65 28 26 78 2c 73 |*/...fwr|ite(&x,s|
|00002570| 69 7a 65 6f 66 28 46 6d | 61 74 72 69 78 29 2c 31 |izeof(Fm|atrix),1|
|00002580| 2c 66 29 3b 0d 09 66 77 | 72 69 74 65 28 6e 61 6d |,f);..fw|rite(nam|
|00002590| 65 2c 73 69 7a 65 6f 66 | 28 63 68 61 72 29 2c 28 |e,sizeof|(char),(|
|000025a0| 6c 6f 6e 67 29 78 2e 6e | 61 6d 6c 65 6e 2c 66 29 |long)x.n|amlen,f)|
|000025b0| 3b 0d 09 73 77 69 74 63 | 68 28 66 69 6c 65 46 6f |;..switc|h(fileFo|
|000025c0| 72 6d 61 74 29 7b 0d 09 | 63 61 73 65 20 30 3a 0d |rmat){..|case 0:.|
|000025d0| 09 09 65 72 72 6f 72 3d | 57 72 69 74 65 44 6f 75 |..error=|WriteDou|
|000025e0| 62 6c 65 73 28 66 2c 65 | 6c 65 6d 65 6e 74 73 2c |bles(f,e|lements,|
|000025f0| 6e 61 6d 65 2c 72 65 61 | 6c 29 3b 0d 09 09 62 72 |name,rea|l);...br|
|00002600| 65 61 6b 3b 0d 09 63 61 | 73 65 20 33 3a 0d 09 09 |eak;..ca|se 3:...|
|00002610| 65 72 72 6f 72 3d 57 72 | 69 74 65 53 68 6f 72 74 |error=Wr|iteShort|
|00002620| 73 28 66 2c 65 6c 65 6d | 65 6e 74 73 2c 6e 61 6d |s(f,elem|ents,nam|
|00002630| 65 2c 72 65 61 6c 29 3b | 0d 09 09 62 72 65 61 6b |e,real);|...break|
|00002640| 3b 0d 09 64 65 66 61 75 | 6c 74 3a 0d 09 09 65 72 |;..defau|lt:...er|
|00002650| 72 6f 72 3d 31 3b 0d 09 | 7d 0d 09 69 66 28 65 72 |ror=1;..|}..if(er|
|00002660| 72 6f 72 29 72 65 74 75 | 72 6e 20 31 32 3b 0d 09 |ror)retu|rn 12;..|
|00002670| 69 66 20 28 69 6d 61 67 | 21 3d 4e 55 4c 4c 29 20 |if (imag|!=NULL) |
|00002680| 7b 0d 09 09 73 77 69 74 | 63 68 28 66 69 6c 65 46 |{...swit|ch(fileF|
|00002690| 6f 72 6d 61 74 29 7b 0d | 09 09 63 61 73 65 20 30 |ormat){.|..case 0|
|000026a0| 3a 0d 09 09 09 65 72 72 | 6f 72 3d 57 72 69 74 65 |:....err|or=Write|
|000026b0| 44 6f 75 62 6c 65 73 28 | 66 2c 65 6c 65 6d 65 6e |Doubles(|f,elemen|
|000026c0| 74 73 2c 6e 61 6d 65 2c | 69 6d 61 67 29 3b 0d 09 |ts,name,|imag);..|
|000026d0| 09 09 62 72 65 61 6b 3b | 0d 09 09 63 61 73 65 20 |..break;|...case |
|000026e0| 33 3a 0d 09 09 09 65 72 | 72 6f 72 3d 57 72 69 74 |3:....er|ror=Writ|
|000026f0| 65 53 68 6f 72 74 73 28 | 66 2c 65 6c 65 6d 65 6e |eShorts(|f,elemen|
|00002700| 74 73 2c 6e 61 6d 65 2c | 69 6d 61 67 29 3b 0d 09 |ts,name,|imag);..|
|00002710| 09 09 62 72 65 61 6b 3b | 0d 09 09 64 65 66 61 75 |..break;|...defau|
|00002720| 6c 74 3a 0d 09 09 09 65 | 72 72 6f 72 3d 31 3b 0d |lt:....e|rror=1;.|
|00002730| 09 09 7d 0d 09 09 69 66 | 28 65 72 72 6f 72 29 72 |..}...if|(error)r|
|00002740| 65 74 75 72 6e 20 31 33 | 3b 0d 09 7d 0d 09 72 65 |eturn 13|;..}..re|
|00002750| 74 75 72 6e 20 30 3b 0d | 7d 0d 0d 2f 2a 20 52 65 |turn 0;.|}../* Re|
|00002760| 61 64 73 20 44 4f 55 42 | 4c 45 53 20 28 38 2d 62 |ads DOUB|LES (8-b|
|00002770| 79 74 65 20 64 6f 75 62 | 6c 65 73 29 20 69 6e 74 |yte doub|les) int|
|00002780| 6f 20 64 6f 75 62 6c 65 | 73 20 2a 2f 0d 64 6f 75 |o double|s */.dou|
|00002790| 62 6c 65 20 2a 52 65 61 | 64 44 6f 75 62 6c 65 73 |ble *Rea|dDoubles|
|000027a0| 28 46 49 4c 45 20 2a 66 | 2c 6c 6f 6e 67 20 65 6c |(FILE *f|,long el|
|000027b0| 65 6d 65 6e 74 73 2c 63 | 68 61 72 20 2a 6e 61 6d |ements,c|har *nam|
|000027c0| 65 29 0d 7b 0d 09 44 4f | 55 42 4c 45 20 2a 44 3b |e).{..DO|UBLE *D;|
|000027d0| 0d 09 64 6f 75 62 6c 65 | 20 2a 64 3b 0d 09 6c 6f |..double| *d;..lo|
|000027e0| 6e 67 20 69 2c 6e 2c 6a | 3b 0d 09 63 6f 6e 73 74 |ng i,n,j|;..const|
|000027f0| 20 64 53 69 7a 65 3d 32 | 30 34 38 2f 73 69 7a 65 | dSize=2|048/size|
|00002800| 6f 66 28 44 4f 55 42 4c | 45 29 3b 09 2f 2a 20 43 |of(DOUBL|E);./* C|
|00002810| 61 6e 20 62 65 20 77 68 | 61 74 65 76 65 72 20 79 |an be wh|atever y|
|00002820| 6f 75 20 77 61 6e 74 2e | 20 42 69 67 67 65 72 20 |ou want.| Bigger |
|00002830| 69 73 20 66 61 73 74 65 | 72 2e 20 2a 2f 0d 09 0d |is faste|r. */...|
|00002840| 09 61 73 73 65 72 74 28 | 73 69 7a 65 6f 66 28 44 |.assert(|sizeof(D|
|00002850| 4f 55 42 4c 45 29 3d 3d | 38 29 3b 0d 09 44 3d 28 |OUBLE)==|8);..D=(|
|00002860| 44 4f 55 42 4c 45 20 2a | 29 6d 61 6c 6c 6f 63 28 |DOUBLE *|)malloc(|
|00002870| 64 53 69 7a 65 2a 73 69 | 7a 65 6f 66 28 44 4f 55 |dSize*si|zeof(DOU|
|00002880| 42 4c 45 29 29 3b 0d 09 | 69 66 28 44 3d 3d 4e 55 |BLE));..|if(D==NU|
|00002890| 4c 4c 29 7b 0d 09 09 70 | 72 69 6e 74 66 28 22 5c |LL){...p|rintf("\|
|000028a0| 6e 45 72 72 6f 72 3a 20 | 6e 6f 74 20 65 6e 6f 75 |nError: |not enou|
|000028b0| 67 68 20 72 6f 6f 6d 20 | 66 6f 72 20 62 75 66 66 |gh room |for buff|
|000028c0| 65 72 5c 6e 22 29 3b 0d | 09 09 72 65 74 75 72 6e |er\n");.|..return|
|000028d0| 20 4e 55 4c 4c 3b 0d 09 | 7d 0d 09 69 66 20 28 21 | NULL;..|}..if (!|
|000028e0| 28 64 20 3d 20 28 64 6f | 75 62 6c 65 20 2a 29 6d |(d = (do|uble *)m|
|000028f0| 61 6c 6c 6f 63 28 65 6c | 65 6d 65 6e 74 73 2a 73 |alloc(el|ements*s|
|00002900| 69 7a 65 6f 66 28 2a 64 | 29 29 29 29 20 7b 0d 09 |izeof(*d|)))) {..|
|00002910| 09 70 72 69 6e 74 66 28 | 22 5c 6e 45 72 72 6f 72 |.printf(|"\nError|
|00002920| 3a 20 56 61 72 69 61 62 | 6c 65 20 25 73 20 74 6f |: Variab|le %s to|
|00002930| 6f 20 62 69 67 20 74 6f | 20 6c 6f 61 64 5c 6e 22 |o big to| load\n"|
|00002940| 2c 6e 61 6d 65 29 3b 0d | 09 09 72 65 74 75 72 6e |,name);.|..return|
|00002950| 20 4e 55 4c 4c 3b 0d 09 | 7d 0d 09 66 6f 72 28 69 | NULL;..|}..for(i|
|00002960| 3d 30 3b 69 3c 65 6c 65 | 6d 65 6e 74 73 3b 29 7b |=0;i<ele|ments;){|
|00002970| 0d 09 09 69 66 28 64 53 | 69 7a 65 3c 65 6c 65 6d |...if(dS|ize<elem|
|00002980| 65 6e 74 73 2d 69 29 6e | 3d 64 53 69 7a 65 3b 0d |ents-i)n|=dSize;.|
|00002990| 09 09 65 6c 73 65 20 6e | 3d 65 6c 65 6d 65 6e 74 |..else n|=element|
|000029a0| 73 2d 69 3b 0d 09 09 69 | 66 28 66 72 65 61 64 28 |s-i;...i|f(fread(|
|000029b0| 44 2c 73 69 7a 65 6f 66 | 28 44 4f 55 42 4c 45 29 |D,sizeof|(DOUBLE)|
|000029c0| 2c 6e 2c 66 29 20 21 3d | 20 6e 29 20 7b 0d 09 09 |,n,f) !=| n) {...|
|000029d0| 09 66 72 65 65 28 64 29 | 3b 0d 09 09 09 66 72 65 |.free(d)|;....fre|
|000029e0| 65 28 44 29 3b 0d 09 09 | 09 72 65 74 75 72 6e 20 |e(D);...|.return |
|000029f0| 4e 55 4c 4c 3b 0d 09 09 | 7d 0d 09 09 66 6f 72 28 |NULL;...|}...for(|
|00002a00| 6a 3d 30 3b 6a 3c 6e 3b | 6a 2b 2b 29 64 5b 69 2b |j=0;j<n;|j++)d[i+|
|00002a10| 2b 5d 3d 44 5b 6a 5d 3b | 09 2f 2a 20 63 6f 6e 76 |+]=D[j];|./* conv|
|00002a20| 65 72 74 20 44 4f 55 42 | 4c 45 20 74 6f 20 64 6f |ert DOUB|LE to do|
|00002a30| 75 62 6c 65 20 2a 2f 0d | 09 7d 0d 09 66 72 65 65 |uble */.|.}..free|
|00002a40| 28 44 29 3b 0d 09 72 65 | 74 75 72 6e 20 64 3b 0d |(D);..re|turn d;.|
|00002a50| 7d 0d 0d 2f 2a 20 52 65 | 61 64 73 20 44 4f 55 42 |}../* Re|ads DOUB|
|00002a60| 4c 45 53 20 28 38 2d 62 | 79 74 65 20 73 68 6f 72 |LES (8-b|yte shor|
|00002a70| 74 20 64 6f 75 62 6c 65 | 73 29 20 69 6e 74 6f 20 |t double|s) into |
|00002a80| 73 68 6f 72 74 73 20 2a | 2f 0d 73 68 6f 72 74 20 |shorts *|/.short |
|00002a90| 2a 52 65 61 64 44 6f 75 | 62 6c 65 73 49 6e 74 6f |*ReadDou|blesInto|
|00002aa0| 53 68 6f 72 74 73 28 46 | 49 4c 45 20 2a 66 2c 6c |Shorts(F|ILE *f,l|
|00002ab0| 6f 6e 67 20 65 6c 65 6d | 65 6e 74 73 2c 63 68 61 |ong elem|ents,cha|
|00002ac0| 72 20 2a 6e 61 6d 65 29 | 0d 7b 0d 09 44 4f 55 42 |r *name)|.{..DOUB|
|00002ad0| 4c 45 20 2a 44 3b 0d 09 | 73 68 6f 72 74 20 2a 64 |LE *D;..|short *d|
|00002ae0| 3b 0d 09 6c 6f 6e 67 20 | 69 2c 6e 2c 6a 3b 0d 09 |;..long |i,n,j;..|
|00002af0| 63 6f 6e 73 74 20 64 53 | 69 7a 65 3d 32 30 34 38 |const dS|ize=2048|
|00002b00| 2f 73 69 7a 65 6f 66 28 | 44 4f 55 42 4c 45 29 3b |/sizeof(|DOUBLE);|
|00002b10| 09 2f 2a 20 43 61 6e 20 | 62 65 20 77 68 61 74 65 |./* Can |be whate|
|00002b20| 76 65 72 20 79 6f 75 20 | 77 61 6e 74 2e 20 42 69 |ver you |want. Bi|
|00002b30| 67 67 65 72 20 69 73 20 | 66 61 73 74 65 72 2e 20 |gger is |faster. |
|00002b40| 2a 2f 0d 09 0d 09 61 73 | 73 65 72 74 28 73 69 7a |*/....as|sert(siz|
|00002b50| 65 6f 66 28 44 4f 55 42 | 4c 45 29 3d 3d 38 29 3b |eof(DOUB|LE)==8);|
|00002b60| 0d 09 44 3d 28 44 4f 55 | 42 4c 45 20 2a 29 6d 61 |..D=(DOU|BLE *)ma|
|00002b70| 6c 6c 6f 63 28 64 53 69 | 7a 65 2a 73 69 7a 65 6f |lloc(dSi|ze*sizeo|
|00002b80| 66 28 44 4f 55 42 4c 45 | 29 29 3b 0d 09 69 66 28 |f(DOUBLE|));..if(|
|00002b90| 44 3d 3d 4e 55 4c 4c 29 | 7b 0d 09 09 70 72 69 6e |D==NULL)|{...prin|
|00002ba0| 74 66 28 22 5c 6e 45 72 | 72 6f 72 3a 20 6e 6f 74 |tf("\nEr|ror: not|
|00002bb0| 20 65 6e 6f 75 67 68 20 | 72 6f 6f 6d 20 66 6f 72 | enough |room for|
|00002bc0| 20 62 75 66 66 65 72 5c | 6e 22 29 3b 0d 09 09 72 | buffer\|n");...r|
|00002bd0| 65 74 75 72 6e 20 4e 55 | 4c 4c 3b 0d 09 7d 0d 09 |eturn NU|LL;..}..|
|00002be0| 69 66 20 28 21 28 64 20 | 3d 20 28 73 68 6f 72 74 |if (!(d |= (short|
|00002bf0| 20 2a 29 6d 61 6c 6c 6f | 63 28 65 6c 65 6d 65 6e | *)mallo|c(elemen|
|00002c00| 74 73 2a 73 69 7a 65 6f | 66 28 2a 64 29 29 29 29 |ts*sizeo|f(*d))))|
|00002c10| 20 7b 0d 09 09 70 72 69 | 6e 74 66 28 22 5c 6e 45 | {...pri|ntf("\nE|
|00002c20| 72 72 6f 72 3a 20 56 61 | 72 69 61 62 6c 65 20 25 |rror: Va|riable %|
|00002c30| 73 20 74 6f 6f 20 62 69 | 67 20 74 6f 20 6c 6f 61 |s too bi|g to loa|
|00002c40| 64 5c 6e 22 2c 6e 61 6d | 65 29 3b 0d 09 09 72 65 |d\n",nam|e);...re|
|00002c50| 74 75 72 6e 20 4e 55 4c | 4c 3b 0d 09 7d 0d 09 66 |turn NUL|L;..}..f|
|00002c60| 6f 72 28 69 3d 30 3b 69 | 3c 65 6c 65 6d 65 6e 74 |or(i=0;i|<element|
|00002c70| 73 3b 29 7b 0d 09 09 69 | 66 28 64 53 69 7a 65 3c |s;){...i|f(dSize<|
|00002c80| 65 6c 65 6d 65 6e 74 73 | 2d 69 29 6e 3d 64 53 69 |elements|-i)n=dSi|
|00002c90| 7a 65 3b 0d 09 09 65 6c | 73 65 20 6e 3d 65 6c 65 |ze;...el|se n=ele|
|00002ca0| 6d 65 6e 74 73 2d 69 3b | 0d 09 09 69 66 28 66 72 |ments-i;|...if(fr|
|00002cb0| 65 61 64 28 44 2c 73 69 | 7a 65 6f 66 28 44 4f 55 |ead(D,si|zeof(DOU|
|00002cc0| 42 4c 45 29 2c 6e 2c 66 | 29 20 21 3d 20 6e 29 20 |BLE),n,f|) != n) |
|00002cd0| 7b 0d 09 09 09 66 72 65 | 65 28 64 29 3b 0d 09 09 |{....fre|e(d);...|
|00002ce0| 09 66 72 65 65 28 44 29 | 3b 0d 09 09 09 72 65 74 |.free(D)|;....ret|
|00002cf0| 75 72 6e 20 4e 55 4c 4c | 3b 0d 09 09 7d 0d 09 09 |urn NULL|;...}...|
|00002d00| 66 6f 72 28 6a 3d 30 3b | 6a 3c 6e 3b 6a 2b 2b 29 |for(j=0;|j<n;j++)|
|00002d10| 64 5b 69 2b 2b 5d 3d 44 | 5b 6a 5d 3b 09 2f 2a 20 |d[i++]=D|[j];./* |
|00002d20| 63 6f 6e 76 65 72 74 20 | 44 4f 55 42 4c 45 20 74 |convert |DOUBLE t|
|00002d30| 6f 20 73 68 6f 72 74 20 | 2a 2f 0d 09 7d 0d 09 66 |o short |*/..}..f|
|00002d40| 72 65 65 28 44 29 3b 0d | 09 72 65 74 75 72 6e 20 |ree(D);.|.return |
|00002d50| 64 3b 0d 7d 0d 0d 2f 2a | 20 57 72 69 74 65 73 20 |d;.}../*| Writes |
|00002d60| 64 6f 75 62 6c 65 73 20 | 61 73 20 44 4f 55 42 4c |doubles |as DOUBL|
|00002d70| 45 53 20 28 38 2d 62 79 | 74 65 20 73 68 6f 72 74 |ES (8-by|te short|
|00002d80| 20 64 6f 75 62 6c 65 73 | 29 20 2a 2f 0d 69 6e 74 | doubles|) */.int|
|00002d90| 20 57 72 69 74 65 44 6f | 75 62 6c 65 73 28 46 49 | WriteDo|ubles(FI|
|00002da0| 4c 45 20 2a 66 2c 6c 6f | 6e 67 20 65 6c 65 6d 65 |LE *f,lo|ng eleme|
|00002db0| 6e 74 73 2c 63 68 61 72 | 20 2a 6e 61 6d 65 2c 64 |nts,char| *name,d|
|00002dc0| 6f 75 62 6c 65 20 2a 64 | 29 0d 7b 0d 09 44 4f 55 |ouble *d|).{..DOU|
|00002dd0| 42 4c 45 20 2a 44 3b 0d | 09 6c 6f 6e 67 20 69 2c |BLE *D;.|.long i,|
|00002de0| 6e 2c 6a 3b 0d 09 63 6f | 6e 73 74 20 64 53 69 7a |n,j;..co|nst dSiz|
|00002df0| 65 3d 32 30 34 38 2f 73 | 69 7a 65 6f 66 28 44 4f |e=2048/s|izeof(DO|
|00002e00| 55 42 4c 45 29 3b 09 2f | 2a 20 43 61 6e 20 62 65 |UBLE);./|* Can be|
|00002e10| 20 77 68 61 74 65 76 65 | 72 20 79 6f 75 20 77 61 | whateve|r you wa|
|00002e20| 6e 74 2e 20 42 69 67 67 | 65 72 20 69 73 20 66 61 |nt. Bigg|er is fa|
|00002e30| 73 74 65 72 2e 20 2a 2f | 0d 09 0d 09 6e 61 6d 65 |ster. */|....name|
|00002e40| 3b 09 2f 2a 20 70 72 65 | 76 65 6e 74 20 22 75 6e |;./* pre|vent "un|
|00002e50| 75 73 65 64 20 61 72 67 | 75 6d 65 6e 74 22 20 77 |used arg|ument" w|
|00002e60| 61 72 6e 69 6e 67 20 2a | 2f 0d 09 61 73 73 65 72 |arning *|/..asser|
|00002e70| 74 28 73 69 7a 65 6f 66 | 28 44 4f 55 42 4c 45 29 |t(sizeof|(DOUBLE)|
|00002e80| 3d 3d 38 29 3b 0d 09 69 | 66 20 28 64 3d 3d 4e 55 |==8);..i|f (d==NU|
|00002e90| 4c 4c 29 20 72 65 74 75 | 72 6e 20 31 34 3b 0d 09 |LL) retu|rn 14;..|
|00002ea0| 44 3d 28 44 4f 55 42 4c | 45 20 2a 29 6d 61 6c 6c |D=(DOUBL|E *)mall|
|00002eb0| 6f 63 28 64 53 69 7a 65 | 2a 73 69 7a 65 6f 66 28 |oc(dSize|*sizeof(|
|00002ec0| 44 4f 55 42 4c 45 29 29 | 3b 0d 09 69 66 28 44 3d |DOUBLE))|;..if(D=|
|00002ed0| 3d 4e 55 4c 4c 29 7b 0d | 09 09 70 72 69 6e 74 66 |=NULL){.|..printf|
|00002ee0| 28 22 5c 6e 45 72 72 6f | 72 3a 20 6e 6f 74 20 65 |("\nErro|r: not e|
|00002ef0| 6e 6f 75 67 68 20 72 6f | 6f 6d 20 66 6f 72 20 62 |nough ro|om for b|
|00002f00| 75 66 66 65 72 5c 6e 22 | 29 3b 0d 09 09 72 65 74 |uffer\n"|);...ret|
|00002f10| 75 72 6e 20 31 35 3b 0d | 09 7d 0d 09 66 6f 72 28 |urn 15;.|.}..for(|
|00002f20| 69 3d 30 3b 69 3c 65 6c | 65 6d 65 6e 74 73 3b 29 |i=0;i<el|ements;)|
|00002f30| 7b 0d 09 09 69 66 28 64 | 53 69 7a 65 3c 65 6c 65 |{...if(d|Size<ele|
|00002f40| 6d 65 6e 74 73 2d 69 29 | 6e 3d 64 53 69 7a 65 3b |ments-i)|n=dSize;|
|00002f50| 0d 09 09 65 6c 73 65 20 | 6e 3d 65 6c 65 6d 65 6e |...else |n=elemen|
|00002f60| 74 73 2d 69 3b 0d 09 09 | 66 6f 72 28 6a 3d 30 3b |ts-i;...|for(j=0;|
|00002f70| 6a 3c 6e 3b 6a 2b 2b 29 | 44 5b 6a 5d 3d 64 5b 69 |j<n;j++)|D[j]=d[i|
|00002f80| 2b 2b 5d 3b 09 2f 2a 20 | 63 6f 6e 76 65 72 74 20 |++];./* |convert |
|00002f90| 64 6f 75 62 6c 65 20 74 | 6f 20 44 4f 55 42 4c 45 |double t|o DOUBLE|
|00002fa0| 20 2a 2f 0d 09 09 69 66 | 28 66 77 72 69 74 65 28 | */...if|(fwrite(|
|00002fb0| 44 2c 73 69 7a 65 6f 66 | 28 44 4f 55 42 4c 45 29 |D,sizeof|(DOUBLE)|
|00002fc0| 2c 6e 2c 66 29 20 21 3d | 20 6e 29 20 7b 0d 09 09 |,n,f) !=| n) {...|
|00002fd0| 09 66 72 65 65 28 44 29 | 3b 0d 09 09 09 72 65 74 |.free(D)|;....ret|
|00002fe0| 75 72 6e 20 31 36 3b 0d | 09 09 7d 0d 09 7d 0d 09 |urn 16;.|..}..}..|
|00002ff0| 66 72 65 65 28 44 29 3b | 0d 09 72 65 74 75 72 6e |free(D);|..return|
|00003000| 20 30 3b 0d 7d 0d 0d 2f | 2a 20 57 72 69 74 65 73 | 0;.}../|* Writes|
|00003010| 20 73 68 6f 72 74 73 20 | 61 73 20 73 68 6f 72 74 | shorts |as short|
|00003020| 73 20 2a 2f 0d 69 6e 74 | 20 57 72 69 74 65 53 68 |s */.int| WriteSh|
|00003030| 6f 72 74 73 28 46 49 4c | 45 20 2a 66 2c 6c 6f 6e |orts(FIL|E *f,lon|
|00003040| 67 20 65 6c 65 6d 65 6e | 74 73 2c 63 68 61 72 20 |g elemen|ts,char |
|00003050| 2a 6e 61 6d 65 2c 73 68 | 6f 72 74 20 2a 64 29 0d |*name,sh|ort *d).|
|00003060| 7b 0d 09 6e 61 6d 65 3b | 09 2f 2a 20 70 72 65 76 |{..name;|./* prev|
|00003070| 65 6e 74 20 22 75 6e 75 | 73 65 64 20 61 72 67 75 |ent "unu|sed argu|
|00003080| 6d 65 6e 74 22 20 77 61 | 72 6e 69 6e 67 20 2a 2f |ment" wa|rning */|
|00003090| 0d 09 61 73 73 65 72 74 | 28 73 69 7a 65 6f 66 28 |..assert|(sizeof(|
|000030a0| 73 68 6f 72 74 29 3d 3d | 32 29 3b 0d 09 69 66 20 |short)==|2);..if |
|000030b0| 28 64 3d 3d 4e 55 4c 4c | 29 20 72 65 74 75 72 6e |(d==NULL|) return|
|000030c0| 20 31 37 3b 0d 09 69 66 | 28 66 77 72 69 74 65 28 | 17;..if|(fwrite(|
|000030d0| 64 2c 73 69 7a 65 6f 66 | 28 73 68 6f 72 74 29 2c |d,sizeof|(short),|
|000030e0| 65 6c 65 6d 65 6e 74 73 | 2c 66 29 20 21 3d 20 65 |elements|,f) != e|
|000030f0| 6c 65 6d 65 6e 74 73 29 | 20 72 65 74 75 72 6e 20 |lements)| return |
|00003100| 31 38 3b 0d 09 72 65 74 | 75 72 6e 20 30 3b 0d 7d |18;..ret|urn 0;.}|
|00003110| 0d 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003120| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003130| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003140| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003150| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003160| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003170| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003180| 00 00 01 00 00 00 01 76 | 00 00 00 76 00 00 00 52 |.......v|...v...R|
|00003190| 22 24 00 42 40 48 40 48 | 42 42 42 26 01 72 00 78 |"$.B@H@H|BBB&.r.x|
|000031a0| 0f d4 82 d1 80 d2 81 b0 | 83 6d 04 90 83 52 01 51 |........|.m...R.Q|
|000031b0| 10 52 65 61 64 4d 41 54 | 4c 41 42 46 69 6c 65 2e |.ReadMAT|LABFile.|
|000031c0| 63 46 02 00 00 00 54 45 | 58 54 43 57 49 45 01 00 |cF....TE|XTCWIE..|
|000031d0| 00 00 54 45 58 54 43 57 | 49 45 01 00 00 78 04 80 |..TEXTCW|IE...x..|
|000031e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000031f0| 00 00 a7 af 02 a3 00 00 | 30 91 00 00 01 c8 fc 51 |........|0......Q|
|00003200| c8 ff de 20 61 10 81 30 | 21 20 49 d0 40 d0 40 d0 |... a..0|! I.@.@.|
|00003210| c0 54 48 20 a1 22 88 4c | df 07 07 2e 5f 4e 75 48 |.TH .".L|...._NuH|
|00003220| e7 c0 c0 72 01 60 06 48 | e7 c0 c0 42 41 20 6f 00 |...r.`.H|...BA o.|
|00003230| 18 22 6f 00 14 42 40 10 | 10 b3 08 56 c8 ff fc 66 |."o..B@.|...V...f|
|00003240| 04 0a 41 00 01 1f 41 00 | 1a 2f 6f 00 10 00 16 4c |..A...A.|./o....L|
|00003250| df 03 03 5c 4f 4e 75 00 | 00 07 78 60 00 02 66 4e |...\ONu.|..x`..fN|
|00003260| 56 fd d2 48 e7 0f 08 28 | 6e 00 08 2f 0c 20 5f a0 |V..H...(|n../. _.|
|00003270| 29 42 67 2f 0c a9 50 3d | 5f ff fe 3e 2e ff fe 70 |)Bg/..P=|_..>...p|
|00003280| 00 00 00 1e 00 29 01 42 | 01 6a 02 79 00 29 01 42 |.....).B|.j.y.).B|
|00003290| 01 6a 02 79 00 00 00 00 | 00 00 00 00 00 00 00 00 |.j.y....|........|
|000032a0| 00 00 00 00 00 48 00 09 | 4d 6f 6e 61 63 6f 00 00 |.....H..|Monaco..|
|000032b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000032c0| 00 00 00 00 00 00 00 00 | 00 03 00 04 00 3d 00 03 |........|.....=..|
|000032d0| 02 6c 03 3d 00 3d 00 03 | 02 6c 03 3d ac 33 ec 32 |.l.=.=..|.l.=.3.2|
|000032e0| 00 00 30 10 00 00 30 10 | 00 00 2c 29 00 00 00 00 |..0...0.|..,)....|
|000032f0| 00 04 00 01 00 01 00 00 | 01 00 00 00 01 76 00 00 |........|.....v..|
|00003300| 00 76 00 00 00 52 01 a1 | 5f e8 1d 62 00 00 00 1c |.v...R..|_..b....|
|00003310| 00 52 00 01 4d 50 53 52 | 00 01 00 12 4d 57 42 42 |.R..MPSR|....MWBB|
|00003320| 00 00 00 2a 03 f0 ff ff | 00 00 00 00 00 00 00 00 |...*....|........|
|00003330| 03 ed ff ff 00 00 00 22 | 00 00 00 00 03 f0 ff ff |......."|........|
|00003340| 00 00 00 6e 00 00 00 00 | 00 00 00 00 00 00 00 00 |...n....|........|
|00003350| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003360| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00003370| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
+--------+-------------------------+-------------------------+--------+--------+