1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
|
#include <stdio.h>
#include <string.h>
#include "../machine.h"
extern int C2F(cvstr) __PARAMS((integer *,integer *,char *,integer *,unsigned long int));
extern void sciprint __PARAMS((char *fmt,...));
extern void mget2 __PARAMS((FILE *fa, integer swap, double *res, integer n, char *type, integer *ierr));
int worldsize(type)
char type[4];
{
char c;
int wsz;
c=type[0];
if (c=='u') c=type[1];
switch ( c )
{
case 'l' : wsz=sizeof(long);
break;
case 's' : wsz=sizeof(short);
break;
case 'c' : wsz=sizeof(char);
break;
case 'd' : wsz=sizeof(double);
break;
case 'f' : wsz=sizeof(float);
break;
}
return wsz;
}
void
readc(flag,nevprt,t,xd,x,nx,z,nz,tvec,ntvec,rpar,nrpar,
ipar,nipar,inptr,insz,nin,outptr,outsz,nout)
/*
ipar[1] = lfil : file name length
ipar[2:4] = fmt : numbers type ascii code
ipar[5] = is there a time record
ipar[6] = n : buffer length in number of records
ipar[7] = maxvoie : record size
ipar[8] = swap
ipar[9] = first : first record to read
ipar[10:9+lfil] = character codes for file name
ipar[10+lfil:9+lfil++ny+ievt] = reading mask
*/
integer *flag,*nevprt,*nx,*nz,*ntvec,*nrpar,ipar[],*nipar;
integer insz[],*nin,outsz[],*nout;
double x[],xd[],z[],tvec[],rpar[];
double *inptr[],*outptr[],*t;
{
char str[100],type[4];
int job = 1,three=3;
FILE *fd;
int n, k, ievt, kmax,/* no,*//* lfil,*/ m, i, irep,/* nm,*/ ierr;
double *buffer,*record;
int *mask;
long offset;
--ipar;
--z;
fd=(FILE *)(long)z[3];
buffer = (z+4);
mask = ipar+11+ipar[1]-ipar[5];
/*
k : record counter within the buffer
kmax : number of records in the buffer
*/
if (*flag==1) {
n = ipar[6];
ievt = ipar[5];
k = (int)z[1];
/* copy current record to output */
record=buffer+(k-1)*ipar[7]-1;
for (i=0;i<outsz[0];i++)
*(outptr[0]+i)=record[mask[ievt+i]];
if (*nevprt>0) {
/* discrete state */
kmax = (int)z[2];
if (k>=kmax&&kmax==n) {
/* read a new buffer */
m=ipar[6]*ipar[7];
F2C(cvstr)(&three,&(ipar[2]),type,&job, strlen(type));
for (i=2;i>=0;i--)
if (type[i]!=' ') { type[i+1]='\0';break;}
ierr=0;
mget2(fd,ipar[8],buffer,m,type,&ierr);
if (ierr>0) {
sciprint("Read error!\n");
fclose(fd);
z[3] = 0.0;
*flag = -1;
return;
}
else if (ierr<0) { /* EOF reached */
kmax=-(ierr+1)/ipar[7];
}
else
kmax=ipar[6];
z[1] = 1.0;
z[2] = kmax;
}
else if (k<kmax)
z[1] = z[1]+1.0;
}
}
else if (*flag==3) {
ievt = ipar[5];
n = ipar[6];
k = (int)z[1];
kmax = (int) z[2];
if (k > kmax && kmax < n) {
if(ievt)
tvec[0] = *t-1.0;
else
tvec[0] = *t*(1.0+0.0000000001);
}
else {
record=buffer+(k-1)*ipar[7]-1;
if(ievt) tvec[0] = record[mask[0]];
}
}
else if (*flag==4) {
F2C(cvstr)(&(ipar[1]),&(ipar[10]),str,&job,strlen(str));
str[ipar[1]] = '\0';
fd = fopen(str,"rb");
if (!fd ) {
sciprint("Could not open the file!\n");
*flag = -1;
return;
}
z[3]=(long)fd;
/* skip first records */
if (ipar[9]>1) {
F2C(cvstr)(&three,&(ipar[2]),type,&job,strlen(type));
for (i=2;i>=0;i--)
if (type[i]!=' ') { type[i+1]='\0';break;}
offset=(ipar[9]-1)*ipar[7]*worldsize(type);
irep = fseek(fd,offset,0) ;
if ( irep != 0 )
{
sciprint("Read error\r\n");
*flag = -1;
fclose(fd);
z[3] = 0.0;
return;
}
}
/* read first buffer */
m=ipar[6]*ipar[7];
F2C(cvstr)(&three,&(ipar[2]),type,&job,strlen(type));
for (i=2;i>=0;i--)
if (type[i]!=' ') { type[i+1]='\0';break;}
mget2(fd,ipar[8],buffer,m,type,&ierr);
if (ierr>0) {
sciprint("Read error!\n");
*flag = -1;
fclose(fd);
z[3] = 0.0;
return;
}
else if (ierr<0) { /* EOF reached */
kmax=-(ierr+1)/ipar[7];
}
else
kmax=ipar[6];
z[1] = 1.0;
z[2] = kmax;
}
else if (*flag==5) {
if(z[3]==0) return;
fclose(fd);
z[3] = 0.0;
}
return;
}
|
