float complex
ctanf(float complex z)
{
z = ctanhf(cpackf(-cimagf(z), crealf(z)));
return (cpackf(cimagf(z), -crealf(z)));
}
float complex
ctanf(float complex z)
{
z = ctanhf(CMPLXF(cimagf(z), crealf(z)));
return (CMPLXF(cimagf(z), crealf(z)));
}
fcomplex
ctanf(fcomplex z) {
float x, y;
fcomplex ans, ct;
x = F_RE(z);
y = F_IM(z);
F_RE(z) = y;
F_IM(z) = -x;
ct = ctanhf(z);
F_RE(ans) = -F_IM(ct);
F_IM(ans) = F_RE(ct);
return (ans);
}
ld = casinhl(ld); TEST_TRACE(C99 7.3.6.3) d = catanh(d); f = catanhf(f); ld = catanhl(ld); TEST_TRACE(C99 7.3.6.4) d = ccosh(d); f = ccoshf(f); ld = ccoshl(ld); TEST_TRACE(C99 7.3.6.5) d = csinh(d); f = csinhf(f); ld = csinhl(ld); TEST_TRACE(C99 7.3.6.6) d = ctanh(d); f = ctanhf(f); ld = ctanhl(ld); TEST_TRACE(C99 7.3.7.1) d = cexp(d); f = cexpf(f); ld = cexpl(ld); TEST_TRACE(C99 7.3.7.2) d = clog(d); f = clogf(f); ld = clogl(ld); TEST_TRACE(C99 7.3.8.1) d = cabs(d); f = cabsf(f); ld = cabsl(ld); TEST_TRACE(C99 7.3.8.2) d = cpow(d, d);
void synth::generate_audio(int16_t* audio_buffer, unsigned buffer_length)
{
int16_t* out = audio_buffer;
uint32_t i,j,k;
float out_l=0.0f;
float out_r=0.0f;
float temp_l=0.0f;
float temp_r=0.0f;
PolyHandler g_channel[10];
//initializing channels with instrument parameter sets
g_channel[0].init(pad_channel_params,k_num_synth_params);
g_channel[1].init(fx_channel_params,k_num_synth_params);
g_channel[2].init(hh_channel_params,k_num_synth_params);
g_channel[3].init(sd_channel_params,k_num_synth_params);
g_channel[4].init(bd2_channel_params,k_num_synth_params);
g_channel[5].init(jyrina_channel_params,k_num_synth_params);
g_channel[6].init(hhdelay_channel_params,k_num_synth_params);
g_channel[7].init(jyrinadelay_channel_params,k_num_synth_params);
g_channel[8].init(sustbass_channel_params,k_num_synth_params);
g_channel[9].init(bass_channel_params,k_num_synth_params);
float division=955.5f; // supposed to be 960 but for some reason drifts a lot TODO: figure out why
float tempo=500000.0f;
//not sure if accurate, could drift/jitter
uint32_t srtick = static_cast<uint32_t>((AUDIO_FREQUENCY/1000000.0f)*(tempo/division));
//song data is split into inner and outer arrays with some decimation to save space without affecting timing
for(i = 0; (i < (buffer_length / sizeof(int16_t) / 2)); ++i)
{
unsigned ee = sizeof(g_song_data_outer) / sizeof(*g_song_data_outer) / 2;
for(j = 0; (j < ee); ++j)
{
const uint8_t *inner = g_song_data_inner + j * 3;
const uint16_t *outer = g_song_data_outer + j * 2;
if((outer[0] * OUTER_SONG_DATA_DIVISOR * srtick) == i)
{
switch(inner[0])
{
//since it turns out sonar is not using standard note off messages but uses 0 velocity note ons
//for note off messages, doing it like this
case k_event_num_standard_note_on:
if(0 < outer[1])
{
g_channel[inner[1]].noteOn(inner[2], static_cast<float>(outer[1]) * (1.0f / 127.0f));
}
else
{
g_channel[inner[1]].noteOff(inner[2]);
}
break;
case k_event_num_standard_pitch_bend:
g_channel[inner[1]].setParameter(k_pitchbend, static_cast<float>(outer[1]) * (1.0f / 16383.0f));
break;
case k_event_num_custom_nrpn:
g_channel[inner[1]].setParameter(inner[2], static_cast<float>(outer[1]) * (1.0f / 16383.0f));
break;
default:
break;
}
}
}
out_l=0.0f;
out_r=0.0f;
// to avoid overflowing or wrapping or just plain nasty clipping,
// running through tanh for tamer clipping.
for(k=0;k<NUM_CHANNELS;k++)
{
if(g_channel[k].getIsActive())
{
//last minute quick hack to lower the level of channels that were too hot to begin with
//should not have been necessary but with this we get acceptable audio levels without crazy distortion
//this was tuned by ear hence the magic numbers
g_channel[k].getSample(temp_l, temp_r);
out_l += (0.9f-(0.015f*(float)k))*ctanhf(0.69f*temp_l);
out_r += (0.9f-(0.015f*(float)k))*ctanhf(0.69f*temp_r);
}
}
// Scaling of 0.9f was too low. Rescales so that output maximum just caps and rounds into 0.999f.
out_l = ctanhf(5.9f * out_l);
out_r = ctanhf(5.9f * out_r);
#if 0
{
static float out_max = -FLT_MAX;
static float out_min = FLT_MAX;
float new_max = fmaxf(fmaxf(fabsf(out_l), fabsf(out_r)), out_max);
float new_min = fminf(fminf(fabsf(out_l), fabsf(out_r)), out_min);
if(new_max > out_max)
{
printf("new output maximum; %1.3f\n", new_max);
out_max = new_max;
}
if(new_min < out_min)
{
printf("new output minimum; %1.3f\n", new_min);
out_min = new_min;
}
}
#endif
// master out, scaling close to maximum dynamic range, should stay within range based on statistics
out[i * 2 + 0] = static_cast<int16_t>(out_l * 32760.f); //left channel;
out[i * 2 + 1] = static_cast<int16_t>(out_r * 32760.f); //right channel
}
}
void
docomplexf (void)
{
#ifndef NO_FLOAT
complex float ca, cb, cc;
float f1;
ca = 1.0 + 1.0 * I;
cb = 1.0 - 1.0 * I;
f1 = cabsf (ca);
fprintf (stdout, "cabsf : %f\n", f1);
cc = cacosf (ca);
fprintf (stdout, "cacosf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = cacoshf (ca);
fprintf (stdout, "cacoshf: %f %fi\n", crealf (cc),
cimagf (cc));
f1 = cargf (ca);
fprintf (stdout, "cargf : %f\n", f1);
cc = casinf (ca);
fprintf (stdout, "casinf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = casinhf (ca);
fprintf (stdout, "casinhf: %f %fi\n", crealf (cc),
cimagf (cc));
cc = catanf (ca);
fprintf (stdout, "catanf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = catanhf (ca);
fprintf (stdout, "catanhf: %f %fi\n", crealf (cc),
cimagf (cc));
cc = ccosf (ca);
fprintf (stdout, "ccosf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = ccoshf (ca);
fprintf (stdout, "ccoshf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = cexpf (ca);
fprintf (stdout, "cexpf : %f %fi\n", crealf (cc),
cimagf (cc));
f1 = cimagf (ca);
fprintf (stdout, "cimagf : %f\n", f1);
cc = clogf (ca);
fprintf (stdout, "clogf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = conjf (ca);
fprintf (stdout, "conjf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = cpowf (ca, cb);
fprintf (stdout, "cpowf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = cprojf (ca);
fprintf (stdout, "cprojf : %f %fi\n", crealf (cc),
cimagf (cc));
f1 = crealf (ca);
fprintf (stdout, "crealf : %f\n", f1);
cc = csinf (ca);
fprintf (stdout, "csinf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = csinhf (ca);
fprintf (stdout, "csinhf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = csqrtf (ca);
fprintf (stdout, "csqrtf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = ctanf (ca);
fprintf (stdout, "ctanf : %f %fi\n", crealf (cc),
cimagf (cc));
cc = ctanhf (ca);
fprintf (stdout, "ctanhf : %f %fi\n", crealf (cc),
cimagf (cc));
#endif
}
TEST(complex, ctanhf) {
ASSERT_EQ(0.0, ctanhf(0));
}
