int32_t voicform(CSOUND *csound, VOICF *p)
{
MYFLT *ar = p->ar;
uint32_t offset = p->h.insdshead->ksmps_offset;
uint32_t early = p->h.insdshead->ksmps_no_end;
uint32_t n, nsmps = CS_KSMPS;
if (p->basef != *p->frequency) {
p->basef = *p->frequency;
SingWave_setFreq(csound, &p->voiced, p->basef);
}
/* OnePole_setPole(&p->onepole, 0.95 - (amp * 0.1)); */
/* Envelope_setTarget(&(p->voiced.envelope), amp); */
/* Envelope_setTarget(&p->noiseEnv, 0.95 - (amp * 0.1)); */
SingWave_setVibFreq(p->voiced, *p->vibf);
Modulatr_setVibAmt(p->voiced.modulator, *p->vibAmt);
/* Set phoneme */
if (p->oldform != *p->formant || p->ph != (int32_t)(0.5+*p->phoneme)) {
p->oldform = *p->formant;
p->ph = (int32_t)(0.5 + *p->phoneme);
csound->Warning(csound, Str("Setting Phoneme: %d %f\n"),
p->ph, p->oldform);
VoicForm_setPhoneme(csound, p, (int32_t) *p->phoneme, p->oldform);
}
/* voicprint(csound, p); */
if (UNLIKELY(offset)) memset(ar, '\0', offset*sizeof(MYFLT));
if (UNLIKELY(early)) {
nsmps -= early;
memset(&ar[nsmps], '\0', early*sizeof(MYFLT));
}
for (n=offset; n<nsmps; n++) {
MYFLT temp;
MYFLT lastOutput;
temp = OnePole_tick(&p->onepole,
OneZero_tick(&p->onezero,
SingWave_tick(csound, &p->voiced)));
// printf("%d: temp=%f ", n, temp);
temp += Envelope_tick(&p->noiseEnv) * Noise_tick(csound, &p->noise);
// printf("%f\n", temp);
lastOutput = FormSwep_tick(csound, &p->filters[0], temp);
// printf("%d: output=%f ", lastOutput);
lastOutput = FormSwep_tick(csound, &p->filters[1], lastOutput);
// printf("%f ", lastOutput);
lastOutput = FormSwep_tick(csound, &p->filters[2], lastOutput);
// printf("%f ", lastOutput);
lastOutput = FormSwep_tick(csound, &p->filters[3], lastOutput);
// printf("%f ", lastOutput);
lastOutput *= p->lastGain;
// printf("%f ", lastOutput);
// printf("->%f\n", lastOutput* AMP_SCALE);
ar[n] = lastOutput * FL(0.22) * AMP_SCALE * *p->amp;
}
return OK;
}
static MYFLT Modulatr_tick(CSOUND *csound, Modulatr *p)
{
MYFLT lastOutput;
lastOutput = Wave_tick(&p->v_time, p->wave->flen, p->wave->ftable,
p->v_rate, FL(0.0));
lastOutput *= p->vibAmt; /* Compute periodic and */
/* random modulations */
lastOutput += OnePole_tick(&p->onepole, SubNoise_tick(csound, &p->noise));
return lastOutput;
}
int bowed(CSOUND *csound, BOWED *p)
{
MYFLT *ar = p->ar;
uint32_t offset = p->h.insdshead->ksmps_offset;
uint32_t early = p->h.insdshead->ksmps_no_end;
uint32_t n, nsmps = CS_KSMPS;
MYFLT amp = (*p->amp)*AMP_RSCALE; /* Normalise */
MYFLT maxVel;
int freq_changed = 0;
if (amp != p->lastamp) {
p->maxVelocity = FL(0.03) + (FL(0.2) * amp);
p->lastamp = amp;
}
maxVel = p->maxVelocity;
if (p->lastpress != *p->bowPress)
p->bowTabl.slope = p->lastpress = *p->bowPress;
/* Set Frequency if changed */
if (p->lastfreq != *p->frequency) {
/* delay - approx. filter delay */
if (p->limit<=*p->frequency)
p->lastfreq = *p->frequency;
else {
p->lastfreq = p->limit;
csound->Warning(csound, Str("frequency too low, set to minimum"));
}
p->baseDelay = CS_ESR / p->lastfreq - FL(4.0);
freq_changed = 1;
}
if (p->lastbeta != *p->betaRatio ||
freq_changed) { /* Reset delays if changed */
p->lastbeta = *p->betaRatio;
DLineL_setDelay(&p->bridgeDelay, /* bow to bridge length */
p->baseDelay * p->lastbeta);
DLineL_setDelay(&p->neckDelay, /* bow to nut (finger) length */
p->baseDelay *(FL(1.0) - p->lastbeta));
}
p->v_rate = *p->vibFreq * p->vibr->flen * csound->onedsr;
if (p->kloop>0 && p->h.insdshead->relesing) p->kloop=1;
if ((--p->kloop) == 0) {
ADSR_setDecayRate(csound, &p->adsr, (FL(1.0) - p->adsr.value) * FL(0.005));
p->adsr.target = FL(0.0);
p->adsr.rate = p->adsr.releaseRate;
p->adsr.state = RELEASE;
}
if (UNLIKELY(offset)) memset(ar, '\0', offset*sizeof(MYFLT));
if (UNLIKELY(early)) {
nsmps -= early;
memset(&ar[nsmps], '\0', early*sizeof(MYFLT));
}
for (n=offset;n<nsmps;n++) {
MYFLT bowVelocity;
MYFLT bridgeRefl=FL(0.0), nutRefl=FL(0.0);
MYFLT newVel=FL(0.0), velDiff=FL(0.0), stringVel=FL(0.0);
MYFLT lastOutput;
bowVelocity = maxVel * ADSR_tick(&p->adsr);
/* Bridge Reflection */
bridgeRefl = - OnePole_tick(&p->reflFilt, p->bridgeDelay.lastOutput);
nutRefl = - p->neckDelay.lastOutput; /* Nut Reflection */
stringVel = bridgeRefl + nutRefl; /* Sum is String Velocity */
velDiff = bowVelocity - stringVel; /* Differential Velocity */
/* Non-Lin Bow Function */
newVel = velDiff * BowTabl_lookup(csound, &p->bowTabl, velDiff);
DLineL_tick(&p->neckDelay, bridgeRefl + newVel); /* Do string */
DLineL_tick(&p->bridgeDelay, nutRefl + newVel); /* propagations */
if (*p->vibAmt > FL(0.0)) {
int32 temp;
MYFLT temp_time, alpha;
/* Tick on vibrato table */
p->v_time += p->v_rate; /* Update current time */
while (p->v_time >= p->vibr->flen) /* Check for end of sound */
p->v_time -= p->vibr->flen; /* loop back to beginning */
while (p->v_time < FL(0.0)) /* Check for end of sound */
p->v_time += p->vibr->flen; /* loop back to beginning */
temp_time = p->v_time;
#ifdef phase_offset
if (p->v_phaseOffset != FL(0.0)) {
temp_time += p->v_phaseOffset; /* Add phase offset */
while (temp_time >= p->vibr->flen) /* Check for end of sound */
temp_time -= p->vibr->flen; /* loop back to beginning */
while (temp_time < FL(0.0)) /* Check for end of sound */
temp_time += p->vibr->flen; /* loop back to beginning */
}
#endif
temp = (int32) temp_time; /* Integer part of time address */
/* fractional part of time address */
alpha = temp_time - (MYFLT)temp;
p->v_lastOutput = p->vibr->ftable[temp]; /* Do linear interpolation */
/* same as alpha*data[temp+1] + (1-alpha)data[temp] */
p->v_lastOutput = p->v_lastOutput +
(alpha * (p->vibr->ftable[temp+1] - p->v_lastOutput));
/* End of vibrato tick */
DLineL_setDelay(&p->neckDelay,
(p->baseDelay * (FL(1.0) - p->lastbeta)) +
(p->baseDelay * *p->vibAmt * p->v_lastOutput));
}
else
DLineL_setDelay(&p->neckDelay,
(p->baseDelay * (FL(1.0) - p->lastbeta)));
lastOutput = BiQuad_tick(&p->bodyFilt, p->bridgeDelay.lastOutput);
ar[n] = lastOutput*AMP_SCALE * amp *FL(1.8);
}
return OK;
}
static MYFLT Modal4_tick(Modal4 *m)
{
MYFLT temp,temp2;
int32 itemp;
MYFLT temp_time, alpha, lastOutput;
int length = (int)m->wave->flen;
m->w_time += m->w_rate; /* Update current time */
if (m->w_time >= length) { /* Check for end of sound */
m->w_time = (MYFLT)(length-1); /* stick at end */
m->w_allDone = 1; /* Information for one-shot use */
}
else if (m->w_time < FL(0.0)) /* Check for end of sound */
m->w_time = FL(0.0); /* stick at beg */
temp_time = m->w_time;
#ifdef phase_offset
if (m->w_phaseOffset != FL(0.0)) {
temp_time += m->w_phaseOffset; /* Add phase offset */
if (temp_time >= length) /* Check for end of sound */
temp_time = length-1; /* stick at end */
else if (temp_time < FL(0.0)) /* check for end of sound */
temp_time = FL(0.0); /* stick at beg */
}
#endif
itemp = (int32) temp_time; /* Integer part of time address */
alpha = temp_time - (MYFLT)itemp; /* fractional part of time address */
lastOutput = m->wave->ftable[itemp]; /* Do linear interpolation */
lastOutput = lastOutput + /* same as alpha*data[temp+1] */
(alpha * (m->wave->ftable[itemp+1] -
lastOutput)); /* + (1-alpha)data[temp] */
temp = m->masterGain *
OnePole_tick(&m->onepole, lastOutput * Envelope_tick(&m->envelope));
temp2 = BiQuad_tick(&m->filters[0], temp);
temp2 += BiQuad_tick(&m->filters[1], temp);
temp2 += BiQuad_tick(&m->filters[2], temp);
temp2 += BiQuad_tick(&m->filters[3], temp);
temp2 = temp2 - (temp2 * m->directGain);
temp2 += m->directGain * temp;
if (m->vibrGain != 0.0) {
/* Tick on vibrato table */
m->v_time += m->v_rate; /* Update current time */
while (m->v_time >= m->vibr->flen) /* Check for end of sound */
m->v_time -= m->vibr->flen; /* loop back to beginning */
while (m->v_time < FL(0.0)) /* Check for end of sound */
m->v_time += m->vibr->flen; /* loop back to beginning */
temp_time = m->v_time;
#ifdef phase_offset
if (m->v_phaseOffset != FL(0.0)) {
temp_time += m->v_phaseOffset; /* Add phase offset */
while (temp_time >= m->vibr->flen) /* Check for end of sound */
temp_time -= m->vibr->flen; /* loop back to beginning */
while (temp_time < FL(0.0)) /* Check for end of sound */
temp_time += m->vibr->flen; /* loop back to beginning */
}
#endif
itemp = (int32) temp_time; /* Integer part of time address */
/* fractional part of time address */
alpha = temp_time - (MYFLT)itemp;
lastOutput = m->vibr->ftable[itemp]; /* Do linear interpolation */
/* same as alpha*data[itemp+1] + (1-alpha)data[temp] */
lastOutput = /*m->v)*/lastOutput +
(alpha * (m->vibr->ftable[itemp+1] - lastOutput));
/* End of vibrato tick */
temp = FL(1.0) + (lastOutput * m->vibrGain); /* Calculate AM */
temp2 = temp * temp2; /* and apply to master out */
}
return (temp2 + temp2);
}
int flute(CSOUND *csound, FLUTE *p)
{
MYFLT *ar = p->ar;
uint32_t offset = p->h.insdshead->ksmps_offset;
uint32_t early = p->h.insdshead->ksmps_no_end;
uint32_t n, nsmps = CS_KSMPS;
MYFLT amp = (*p->amp)*AMP_RSCALE; /* Normalise */
MYFLT temp;
int v_len = (int)p->vibr->flen;
MYFLT *v_data = p->vibr->ftable;
MYFLT v_time = p->v_time;
MYFLT vibGain = *p->vibAmt;
MYFLT jetRefl, endRefl, noisegain;
if (amp!=p->lastamp) { /* If amplitude has changed */
/* This should be controlled by attack */
ADSR_setAttackRate(csound, &p->adsr, amp * FL(0.02));
p->maxPress = (FL(1.1) + (amp * FL(0.20))) / FL(0.8);
p->outputGain = amp + FL(0.001);
p->lastamp = amp;
}
p->v_rate = *p->vibFreq * v_len * csound->onedsr;
/* Start SetFreq */
if (p->lastFreq != *p->frequency) { /* It changed */
p->lastFreq = *p->frequency;
if (p->limit>p->lastFreq) {
p->lastFreq = p->limit;
csound->Warning(csound, Str("frequency too low, set to minimum"));
}
p->lastJet = *p->jetRatio;
/* freq = (2/3)*p->frequency as we're overblowing here */
/* but 1/(2/3) is 1.5 so multiply for speed */
/* Length - approx. filter delay */
temp = FL(1.5)* CS_ESR / p->lastFreq - FL(2.0);
DLineL_setDelay(&p->boreDelay, temp); /* Length of bore tube */
DLineL_setDelay(&p->jetDelay, temp * p->lastJet); /* jet delay shorter */
}
else if (*p->jetRatio != p->lastJet) { /* Freq same but jet changed */
p->lastJet = *p->jetRatio;
/* Length - approx. filter delay */
temp = FL(1.5)* CS_ESR / p->lastFreq - FL(2.0);
DLineL_setDelay(&p->jetDelay, temp * p->lastJet); /* jet delay shorter */
}
/* End SetFreq */
if (p->kloop>FL(0.0) && p->h.insdshead->relesing) p->kloop=FL(1.0);
if ((--p->kloop) == 0) {
p->adsr.releaseRate = p->adsr.value / (*p->dettack * CS_ESR);
p->adsr.target = FL(0.0);
p->adsr.rate = p->adsr.releaseRate;
p->adsr.state = RELEASE;
}
noisegain = *p->noiseGain; jetRefl = *p->jetRefl; endRefl = *p->endRefl;
if (UNLIKELY(offset)) memset(ar, '\0', offset*sizeof(MYFLT));
if (UNLIKELY(early)) {
nsmps -= early;
memset(&ar[nsmps], '\0', early*sizeof(MYFLT));
}
for (n=offset;n<nsmps;n++) {
int32 temp;
MYFLT temf;
MYFLT temp_time, alpha;
MYFLT pressDiff;
MYFLT randPress;
MYFLT breathPress;
MYFLT lastOutput;
MYFLT v_lastOutput;
breathPress = p->maxPress * ADSR_tick(&p->adsr); /* Breath Pressure */
randPress = noisegain*Noise_tick(csound,&p->noise); /* Random Deviation */
/* Tick on vibrato table */
v_time += p->v_rate; /* Update current time */
while (v_time >= v_len) /* Check for end of sound */
v_time -= v_len; /* loop back to beginning */
while (v_time < FL(0.0)) /* Check for end of sound */
v_time += v_len; /* loop back to beginning */
temp_time = v_time;
#ifdef phase_offset
if (p->v_phaseOffset != FL(0.0)) {
temp_time += p->v_phaseOffset;/* Add phase offset */
while (temp_time >= v_len) /* Check for end of sound */
temp_time -= v_len; /* loop back to beginning */
while (temp_time < FL(0.0)) /* Check for end of sound */
temp_time += v_len; /* loop back to beginning */
}
#endif
temp = (int32) temp_time; /* Integer part of time address */
/* fractional part of time address */
alpha = temp_time - (MYFLT)temp;
v_lastOutput = v_data[temp]; /* Do linear interpolation */
/* same as alpha*data[temp+1] + (1-alpha)data[temp] */
v_lastOutput += (alpha * (v_data[temp+1] - v_lastOutput));
/* End of vibrato tick */
randPress += vibGain * v_lastOutput; /* + breath vibrato */
randPress *= breathPress; /* All scaled by Breath Pressure */
temf = OnePole_tick(&p->filter, DLineL_lastOut(&p->boreDelay));
temf = DCBlock_tick(&p->dcBlock, temf); /* Block DC on reflection */
pressDiff = breathPress + randPress /* Breath Pressure */
- (jetRefl * temf); /* - reflected */
pressDiff = DLineL_tick(&p->jetDelay, pressDiff); /* Jet Delay Line */
pressDiff = JetTabl_lookup(pressDiff) /* Non-Lin Jet + reflected */
+ (endRefl * temf);
/* Bore Delay and "bell" filter */
lastOutput = FL(0.3) * DLineL_tick(&p->boreDelay, pressDiff);
lastOutput *= p->outputGain;
ar[n] = lastOutput*AMP_SCALE*FL(1.4);
}
p->v_time = v_time;
return OK;
}