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
|
#include <stdio.h>
#include <string.h>
#include <math.h>
#include "../machine.h"
extern void sciprint __PARAMS((char *fmt,...));
void mput2 __PARAMS((FILE *fa, integer swap, double *res, integer n, char *type, integer *ierr));
void
writeau(flag,nevprt,t,xd,x,nx,z,nz,tvec,ntvec,rpar,nrpar,
ipar,nipar,inptr,insz,nin,outptr,outsz,nout)
integer *flag,*nevprt,*nx,*nz,*ntvec,*nrpar,ipar[],*nipar,insz[],*nin,outsz[],*nout;
double x[],xd[],z[],tvec[],rpar[];
double *inptr[],*outptr[],*t;
/*
ipar[1] = lfil : file name length
ipar[2:4] = fmt : numbers type ascii code
ipar[5] = n : buffer length in number of records
ipar[6] = swap
ipar[7:6+lfil] = character codes for file name
*/
{
FILE *fd;
int n, k,/* m,*/ i, ierr;
double *buffer,*record;
/* long offset;*/
int SCALE = 32768;
int BIAS = 132;
int CLIP = 32635;
int OFFSET = 335;
double y;
int sig;
int e;
double f;
--ipar;
--z;
fd=(FILE *)(long)z[2];
buffer = (z+3);
ierr=0;
/*
k : record counter within the buffer
*/
if (*flag==2&&*nevprt>0)
{ /* add a new record to the buffer */
n = ipar[5];
k = (int)z[1];
/* copy current record to output
printf("%i\n",k);*/
record=buffer+(k-1)*(*nin);
for (i=0;i<*nin;i++)
{
y= *inptr[i];
y=SCALE*y;
if (y<0.0)
{
y=-y;
sig=-1;
}
else
sig=1;
if(y>CLIP)
y=CLIP;
y=y+BIAS;
f=frexp(y,&e);
y=64*sig-16*e- (int) (32*f)+OFFSET;
record[i] = y;
}
if (k<n)
z[1] = z[1]+1.0;
else {
mput2(fd,ipar[6],buffer,ipar[5]*(*nin),"uc",&ierr);
if(ierr!=0) {
*flag = -3;
return;
}
z[1] = 1.0;
}
}
else if (*flag==4) {
fd = fopen("/dev/audio","wb");
if (!fd ) {
sciprint("Could not open /dev/audio!\n");
*flag = -3;
return;
}
z[2]=(double)(long)fd;
z[1] = 1.0;
}
else if (*flag==5) {
if(z[2]==0) return;
k =(int) z[1];
if (k>1) {/* flush rest of buffer */
mput2(fd,ipar[6],buffer,(k-1)*(*nin),"uc",&ierr);
if(ierr!=0) {
*flag = -3;
return;
}
}
fclose(fd);
z[2] = 0.0;
}
return;
}
|
