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
|
/*
pvread.c:
Copyright (C) 1992, 2000 Richard Karpen, Richard Dobson
This file is part of Csound.
The Csound Library is free software; you can redistribute it
and/or modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
Csound is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with Csound; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA
02110-1301 USA
*/
/**************************************************************/
/***********pvread ********************************************/
/*** By Richard Karpen - July-October 1992*********************/
/**************************************************************/
#include "pvoc.h"
#include <math.h>
/*RWD 10:9:2000 read pvocex file format */
#include "pvfileio.h"
static int32_t pvocex_loadfile(CSOUND *, const char *fname, PVREAD *p);
#define WLN 1 /* time window is WLN*2*ksmps long */
#define OPWLEN (2*WLN*ksmps) /* manifest used for final time wdw */
static void FetchInOne(
float *inp, /* pointer to input data */
MYFLT *buf, /* where to put our nice mag/pha pairs */
int32 fsize, /* frame size we're working with */
MYFLT pos, /* fractional frame we want */
int32 mybin)
{
float *frame0;
float *frame1;
int32 base;
MYFLT frac;
int32 twmybin = mybin+mybin; /* Always used thus */
/***** WITHOUT INFO ON WHERE LAST FRAME IS, MAY 'INTERP' BEYOND IT ****/
base = (int32)pos; /* index of basis frame of interpolation */
frac = ((MYFLT)(pos - (MYFLT)base));
/* & how close to get to next */
frame0 = inp + ((int32) fsize + 2L) * base + twmybin;
frame1 = frame0 + ((int32) fsize + 2L); /* addresses of both frames */
if (frac != 0.0) { /* must have 2 cases to avoid possible segmentation */
/* violations and failed computes, else may interp */
/* beyond valid data */
buf[0] = frame0[0] + frac * (frame1[0] - frame0[0]);
buf[1] = frame0[1] + frac * (frame1[1] - frame0[1]);
}
else {
/* frac is 0.0 i.e. just copy the source frame */
buf[0] = frame0[0];
buf[1] = frame0[1];
}
}
int32_t pvreadset_(CSOUND *csound, PVREAD *p, int32_t stringname)
{
char pvfilnam[256];
if (stringname==0){
if (csound->ISSTRCOD(*p->ifilno))
strNcpy(pvfilnam,get_arg_string(csound, *p->ifilno), MAXNAME);
else csound->strarg2name(csound, pvfilnam, p->ifilno, "pvoc.",0);
}
else strNcpy(pvfilnam, ((STRINGDAT *)p->ifilno)->data, MAXNAME);
if (pvocex_loadfile(csound, pvfilnam, p) == OK) {
p->prFlg = 1;
p->mybin = MYFLT2LRND(*p->ibin);
return OK;
}
return NOTOK;
}
int32_t pvreadset(CSOUND *csound, PVREAD *p){
return pvreadset_(csound,p,0);
}
int32_t pvreadset_S(CSOUND *csound, PVREAD *p){
return pvreadset_(csound,p,1);
}
int32_t pvread(CSOUND *csound, PVREAD *p)
{
MYFLT frIndx;
MYFLT buf[2];
int32_t size = pvfrsiz(p);
if (UNLIKELY((frIndx = *p->ktimpnt * p->frPrtim) < 0)) goto err1;
if (frIndx > p->maxFr) { /* not past last one */
frIndx = (MYFLT)p->maxFr;
if (p->prFlg) {
p->prFlg = 0; /* false */
csound->Warning(csound, Str("PVOC ktimpnt truncated to last frame"));
}
}
FetchInOne(p->frPtr, &(buf[0]), size, frIndx, p->mybin);
*p->kfreq = buf[1];
*p->kamp = buf[0];
return OK;
err1:
return csound->PerfError(csound, &(p->h), Str("PVOC timpnt < 0"));
}
static int32_t pvocex_loadfile(CSOUND *csound, const char *fname, PVREAD *p)
{
PVOCEX_MEMFILE pp;
if (UNLIKELY(csound->PVOCEX_LoadFile(csound, fname, &pp) != 0)) {
return csound->InitError(csound, Str("PVREAD cannot load %s"), fname);
}
/* have to reject m/c files for now, until opcodes upgraded */
if (UNLIKELY(pp.chans > 1)) {
return csound->InitError(csound, Str("pvoc-ex file %s is not mono"), fname);
}
/* ignore the window spec until we can use it! */
p->frSiz = pp.fftsize;
p->frPtr = (float*) pp.data;
p->baseFr = 0; /* point to first data frame */
p->maxFr = pp.nframes - 1;
p->asr = pp.srate;
/* highest possible frame index */
/* factor by which to mult expand phase diffs (ratio of samp spacings) */
p->frPrtim = CS_ESR / ((MYFLT) pp.overlap);
return OK;
}
|
