static void pluck(t_bowedbar *x, float amplitude)
{
long i,j;
int pluckLen;
float temp;
pluckLen = (int)(x->length/x->modes[x->NR_MODES-1]);
for (j=1;j<pluckLen/2;j++) {
temp = amplitude*2.0*Noise_tick();
for(i=0; i<x->NR_MODES; i++)
DLineN_tick(&x->delay[i], temp*j/pluckLen*x->gains[i]);
}
for (j=pluckLen/2;j>0;j--) {
temp = amplitude*2.0*Noise_tick();
for(i=0; i<x->NR_MODES; i++)
DLineN_tick(&x->delay[i], temp*j/pluckLen*x->gains[i]);;
}
}
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 SubNoise_tick(CSOUND *csound, SubNoise *p)
{
MYFLT last;
if (p->counter==0) {
last = p->lastOutput = Noise_tick(csound, &p->lastOutput);
p->counter = p->howOften;
}
else {
(p->counter)--;
last = p->lastOutput;
}
return last;
}
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;
}
int clarin(CSOUND *csound, CLARIN *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 nGain = *p->noiseGain;
int v_len = (int)p->vibr->flen;
MYFLT *v_data = p->vibr->ftable;
MYFLT vibGain = *p->vibAmt;
MYFLT vTime = p->v_time;
if (p->envelope.rate==FL(0.0)) {
p->envelope.rate = amp /(*p->attack*CS_ESR);
p->envelope.value = p->envelope.target = FL(0.55) + amp*FL(0.30);
}
p->outputGain = amp + FL(0.001);
DLineL_setDelay(&p->delayLine, /* length - approx filter delay */
(CS_ESR/ *p->frequency) * FL(0.5) - FL(1.5));
p->v_rate = *p->vibFreq * p->vibr->flen * csound->onedsr;
/* Check to see if into decay yet */
if (p->kloop>0 && p->h.insdshead->relesing) p->kloop=1;
if ((--p->kloop) == 0) {
p->envelope.state = 1; /* Start change */
p->envelope.rate = p->envelope.value / (*p->dettack * CS_ESR);
p->envelope.target = FL(0.0);
#ifdef BETA
csound->Message(csound, "Set off phase time = %f Breath v,r = %f, %f\n",
(MYFLT) CS_KCNT * CS_ONEDKR,
p->envelope.value, p->envelope.rate);
#endif
}
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 pressureDiff;
MYFLT breathPressure;
int32 temp;
MYFLT temp_time, alpha;
MYFLT nextsamp;
MYFLT v_lastOutput;
MYFLT lastOutput;
breathPressure = Envelope_tick(&p->envelope);
breathPressure += breathPressure * nGain * Noise_tick(csound,&p->noise);
/* Tick on vibrato table */
vTime += p->v_rate; /* Update current time */
while (vTime >= v_len) /* Check for end of sound */
vTime -= v_len; /* loop back to beginning */
while (vTime < FL(0.0)) /* Check for end of sound */
vTime += v_len; /* loop back to beginning */
temp_time = vTime;
#ifdef have_phase
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 */
breathPressure += breathPressure * vibGain * v_lastOutput;
pressureDiff = OneZero_tick(&p->filter, /* differential pressure */
DLineL_lastOut(&p->delayLine));
pressureDiff = (-FL(0.95)*pressureDiff) - breathPressure;
/* of reflected and mouth */
nextsamp = pressureDiff * ReedTabl_LookUp(&p->reedTable,pressureDiff);
nextsamp = breathPressure + nextsamp;
/* perform scattering in economical way */
lastOutput = DLineL_tick(&p->delayLine, nextsamp);
lastOutput *= p->outputGain;
ar[n] = lastOutput*AMP_SCALE;
}
p->v_time = vTime;
return OK;
}
