1 |
#define _FILE_OFFSET_BITS 64 |
2 |
|
3 |
#include <stdio.h> |
4 |
#include <stdlib.h> |
5 |
#include <sys/types.h> |
6 |
#include <sys/stat.h> |
7 |
#include <string.h> |
8 |
|
9 |
#include "params.h" |
10 |
#include "tcProps.h" |
11 |
#include "readWrite.h" |
12 |
|
13 |
|
14 |
#define BUFFER_SIZE 2000 |
15 |
int isScanned; |
16 |
FILE* inFile; |
17 |
|
18 |
|
19 |
struct linkedPos{ |
20 |
|
21 |
fpos_t *myPos; |
22 |
double timeStamp; |
23 |
double Hmat[3][3]; |
24 |
struct linkedPos* next; |
25 |
}; |
26 |
|
27 |
struct staticPos{ |
28 |
|
29 |
fpos_t *myPos; |
30 |
double timeStamp; |
31 |
double Hmat[3][3]; |
32 |
}; |
33 |
struct staticPos* posArray; |
34 |
|
35 |
|
36 |
void closeFile( void ){ |
37 |
|
38 |
fclose( inFile ); |
39 |
} |
40 |
|
41 |
int setFrames( char* inFile ){ |
42 |
|
43 |
int nFrames = 0; |
44 |
int i,j,k; |
45 |
struct linkedPos* headPos; |
46 |
struct linkedPos* currPos; |
47 |
fpos_t *currPT; |
48 |
char readBuffer[BUFFER_SIZE]; |
49 |
char trash[BUFFER_SIZE]; |
50 |
char* foo; |
51 |
int lineNum = 0; |
52 |
|
53 |
|
54 |
|
55 |
inFile = fopen(inName); |
56 |
if(inFile ==NULL){ |
57 |
fprintf(stderr, |
58 |
"Error opening file \"%s\"\n", |
59 |
inName); |
60 |
exit(0); |
61 |
} |
62 |
|
63 |
|
64 |
headPos = (struct linkedPos*)malloc(sizeof(struct linkedPos)); |
65 |
currPos = headPos; |
66 |
while( !feof( inFile ) ){ |
67 |
|
68 |
currPT = (fpos_t *)malloc(sizeof(fpos_t)); |
69 |
fgetpos(inFile, currPT); |
70 |
|
71 |
fgets( readBuffer, sizeof( readBuffer ), inFile ); |
72 |
lineNum++; |
73 |
if( feof( inFile ) ){ |
74 |
fprintf( stderr, |
75 |
"File \"%s\" ended unexpectedly at line %d\n", |
76 |
inName, |
77 |
lineNum ); |
78 |
exit(0); |
79 |
} |
80 |
|
81 |
currPos->next = (struct linkedPos*)malloc(sizeof(struct linkedPos)); |
82 |
currPos = currPos->next; |
83 |
currPos->myPos = currPT; |
84 |
nFrames++; |
85 |
|
86 |
i = atoi(readBuffer); |
87 |
|
88 |
fgets( readBuffer, sizeof( readBuffer ), inFile ); |
89 |
lineNum++; |
90 |
if( feof( inFile ) ){ |
91 |
fprintf( stderr, |
92 |
"File \"%s\" ended unexpectedly at line %d\n", |
93 |
inName, |
94 |
lineNum ); |
95 |
exit(0); |
96 |
} |
97 |
|
98 |
|
99 |
// parse the comment line |
100 |
|
101 |
foo = strtok(readBuffer, " ,;\t"); |
102 |
|
103 |
if(foo == NULL){ |
104 |
fprintf(stderr, |
105 |
"error in reading time from %s at line %d\n", |
106 |
inName, lineNum ); |
107 |
exit(0); |
108 |
} |
109 |
currPos->timeStamp = atof( foo ); |
110 |
|
111 |
// get the Hx vector |
112 |
|
113 |
foo = strtok(NULL, " ,;\t"); |
114 |
if(foo == NULL){ |
115 |
fprintf(stderr, |
116 |
"error in reading Hx[0] from %s at line %d\n", |
117 |
inName, lineNum ); |
118 |
exit(0); |
119 |
} |
120 |
currPos->Hmat[0][0] = atof( foo ); |
121 |
|
122 |
foo = strtok(NULL, " ,;\t"); |
123 |
if(foo == NULL){ |
124 |
fprintf(stderr, |
125 |
"error in reading Hx[1] from %s at line %d\n", |
126 |
inName, lineNum ); |
127 |
exit(0); |
128 |
} |
129 |
currPos->Hmat[1][0] = atof( foo ); |
130 |
|
131 |
foo = strtok(NULL, " ,;\t"); |
132 |
if(foo == NULL){ |
133 |
fprintf(stderr, |
134 |
"error in reading Hx[2] from %s at line %d\n", |
135 |
inName, lineNum ); |
136 |
exit(0); |
137 |
} |
138 |
currPos->Hmat[2][0] = atof( foo ); |
139 |
|
140 |
// get the Hy vector |
141 |
|
142 |
foo = strtok(NULL, " ,;\t"); |
143 |
if(foo == NULL){ |
144 |
fprintf(stderr, |
145 |
"error in reading Hy[0] from %s at line %d\n", |
146 |
inName, lineNum ); |
147 |
exit(0); |
148 |
} |
149 |
currPos->Hmat[0][1] = atof( foo ); |
150 |
|
151 |
foo = strtok(NULL, " ,;\t"); |
152 |
if(foo == NULL){ |
153 |
fprintf(stderr, |
154 |
"error in reading Hy[1] from %s at line %d\n", |
155 |
inName, lineNum ); |
156 |
exit(0); |
157 |
} |
158 |
currPos->Hmat[1][1] = atof( foo ); |
159 |
|
160 |
foo = strtok(NULL, " ,;\t"); |
161 |
if(foo == NULL){ |
162 |
fprintf(stderr, |
163 |
"error in reading Hy[2] from %s at line %d\n", |
164 |
inName, lineNum ); |
165 |
exit(0); |
166 |
} |
167 |
currPos->Hmat[2][1] = atof( foo ); |
168 |
|
169 |
// get the Hz vector |
170 |
|
171 |
foo = strtok(NULL, " ,;\t"); |
172 |
if(foo == NULL){ |
173 |
fprintf(stderr, |
174 |
"error in reading Hz[0] from %s at line %d\n", |
175 |
inName, lineNum ); |
176 |
exit(0); |
177 |
} |
178 |
currPos->Hmat[0][2] = atof( foo ); |
179 |
|
180 |
foo = strtok(NULL, " ,;\t"); |
181 |
if(foo == NULL){ |
182 |
fprintf(stderr, |
183 |
"error in reading Hz[1] from %s at line %d\n", |
184 |
inName, lineNum ); |
185 |
exit(0); |
186 |
} |
187 |
currPos->Hmat[1][2] = atof( foo ); |
188 |
|
189 |
foo = strtok(NULL, " ,;\t"); |
190 |
if(foo == NULL){ |
191 |
fprintf(stderr, |
192 |
"error in reading Hz[2] from %s at line %d\n", |
193 |
inName, lineNum ); |
194 |
exit(0); |
195 |
} |
196 |
currPos->Hmat[2][2] = atof( foo ); |
197 |
|
198 |
// skip the atoms |
199 |
|
200 |
for(j=0; j<i; j++){ |
201 |
|
202 |
fgets( readBuffer, sizeof( readBuffer ), inFile ); |
203 |
lineNum++; |
204 |
if( feof( inFile ) ){ |
205 |
fprintf( stderr, |
206 |
"File \"%s\" ended unexpectedly at line %d," |
207 |
" with atom %d\n", |
208 |
inName, |
209 |
lineNum, |
210 |
j ); |
211 |
exit(0); |
212 |
} |
213 |
} |
214 |
} |
215 |
|
216 |
isScanned = 1; |
217 |
return nFrames; |
218 |
} |