[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octal-dev] SoX flanger
|
From: |
Avelino |
|
Subject: |
[Octal-dev] SoX flanger |
|
Date: |
Mon, 26 Jun 2000 19:31:38 +0100 (BST) |
Sox flanger sounds great!
I have test it with this configuration:
dtosquare --> svfilter --> soxflanger --> delay --> master
... and sounds very nice ;-)
Bye.
-------------------------------------------------------------------------------
Avelino Herrera Morales. E-Mail: address@hidden
address@hidden
HTTP: http://www.avelino.turincon.com
http://www.geocities.com/avelinoherrera/
Powered by SunOS 5.7
Centro Superior de Informatica
UNIVERSIDAD DE LA LAGUNA.
-------------------------------------------------------------------------------
Received: from mta5.snfc21.pbi.net (mta5.snfc21.pbi.net [206.13.28.241])
by mescaline.gnu.org (8.9.1a/8.9.1) with ESMTP id XAA25210
for <address@hidden>; Mon, 26 Jun 2000 23:17:06 -0400
Received: from pacbell.net ([206.170.64.80])
by mta5.snfc21.pbi.net (Sun Internet Mail Server sims.3.5.2000.01.05.12.18.p9)
with ESMTP id <address@hidden> for address@hidden; Mon,
26 Jun 2000 18:33:07 -0700 (PDT)
Date: Mon, 26 Jun 2000 18:31:06 -0700
From: Neil Nelson <address@hidden>
To: address@hidden
Message-id: <address@hidden>
MIME-version: 1.0
X-Mailer: Mozilla 4.72 [en] (X11; U; Linux 2.2.14-5.0 i686)
Content-type: multipart/alternative;
boundary="------------24819DF4628C2DE0E67AAF2F"
X-Accept-Language: en
Subject: [Octal-dev] soxchorus.c - SoX chorus.c converted to Octal
Sender: address@hidden
Errors-To: address@hidden
X-BeenThere: address@hidden
X-Mailman-Version: 1.2 (beta 1)
Precedence: bulk
Reply-To: address@hidden
List-Id: OCTAL tracker http://www.gnu.org/software/octal <octal-dev.gnu.org>
--------------24819DF4628C2DE0E67AAF2F
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit
After working most of the day to eliminate clicks on some settings
associated with
the cycle speed of a voice and not being successful, I have released
this version
so that anyone else may work on the problem. I had this similar problem
with the
soxflanger.c but was able to overcome it. We could use an ability to
write the
final audio output to disk in, say, the signed, 16-bit, little-endian
audio format
seen in output.c.
Additionally, I have made a number of changes to flanger.c but will wait
to release
that code to avoid posting too many copies quickly.
/* SOXCHORUS.C
*
* This code has been adapted from the SoX chorus.c code:
* August 24, 1998
*
* Copyright (C) 1998 Juergen Mueller And Sundry Contributors
* This source code is freely redistributable and may be used for
* any purpose. This copyright notice must be maintained.
* Juergen Mueller And Sundry Contributors are not responsible for
* the consequences of using this software.
*
* Conversion to Octal made by:
*
* Copyright 2000 Neil Nelson <address@hidden> June 26, 2000.
*
* This software is distributed under the terms of the
* GNU General Public License (GPL). Read the included file
* COPYING for more information.
*
*
* Chorus effect.
*
* Flow diagram scheme for n delays ( 1 <= n <= max_chorus ):
*
* * gain-in ___
* ibuff -----+--------------------------------------------->| |
* | _________ | |
* | | | * decay 1 | |
* +---->| delay 1 |----------------------------->| |
* | |_________| | |
* | /|\ | |
* : | | |
* : +-----------------+ +--------------+ | + |
* : | Delay control 1 |<--| mod. speed 1 | | |
* : +-----------------+ +--------------+ | |
* | _________ | |
* | | | * decay n | |
* +---->| delay n |----------------------------->| |
* |_________| | |
* /|\ |___|
* | |
* +-----------------+ +--------------+ | *
gain-out
* | Delay control n |<--| mod. speed n | |
* +-----------------+ +--------------+
+----->obuff
*
*
* The delay i is controled by a sine or triangle modulation i ( 1 <= i
<= n).
*
* Usage: // SoX command line format.
* chorus gain-in gain-out delay-1 decay-1 speed-1 depth-1 -s1|t1 [
* delay-2 decay-2 speed-2 depth-2 -s2|-t2 ... ]
*
* Where:
* gain-in, decay-1 ... decay-n : 0.0 ... 1.0 volume
* gain-out : 0.0 ... volume
* delay-1 ... delay-n : 20.0 ... 100.0 msec
* speed-1 ... speed-n : 0.1 ... 5.0 Hz modulation 1 ... n
* depth-1 ... depth-n : 0.0 ... 10.0 msec modulated delay 1 ... n
* -s1 ... -sn : modulation by sine 1 ... n
* -t1 ... -tn : modulation by triangle 1 ... n
*
* Note:
* when decay is close to 1.0, the samples can begin clipping and the
output
* can saturate!
*
* Hint:
* 1 / out-gain < gain-in ( 1 + decay-1 + ... + decay-n )
*
* Potential Improvements:
* (1) Except for the expected fixed number of sxc_params, this program
* could handle an arbitrary number of chorus voices. What trouble
* would it be to provide for an arbitrary number of sxc_params?
*
*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "util.h"
#include "machine.h"
#define PI2 (6.28318530717958465692)
#define MOD_SINE 0
#define MOD_TRIANGLE 1
#define MOD_SILENT 2
#define DYN_RANGE 10
#define max_chorus 3
/* The parameters are: */
enum {
ix_gain_in, ix_gain_out,
ix_sxcdecay, ix_sxcdelay, ix_speed, ix_depth, ix_mod_type,
ix_sxcdecay_1, ix_sxcdelay_1, ix_speed_1, ix_depth_1, ix_mod_type_1,
ix_sxcdecay_2, ix_sxcdelay_2, ix_speed_2, ix_depth_2, ix_mod_type_2
} param_index;
param_spec sxc_params[] = {
/* Param 1: gain_in. */
{
small,
slider,
"Gain In",
"Gain In",
0, /* 0 gain */
255, /* No attenuation */
200 /* Approx. 80% of dynamic range */
},
/* Param 2: gain_out. */
{
small,
slider,
"Gain Out",
"Gain Out",
0, /* 0 gain */
255, /* No attenuation */
255
},
/* Param 3: sxcdecay[0]. */
{
small,
slider,
"Decay 0",
"Decay of 0th chorus voice.",
0,
255, /* 1 */
200 /* Approx. 80% of dynamic range */
},
/* Param 4: sxcdelay[0]. */
{
small,
slider,
"Delay 0",
"Delay: 20 to 100 msec., for 0th voice.",
0, /* Minimum 20 msec. delay */
255, /* Maximum 100 msec. delay */
128 /* Approx. 60 msec. delay */
},
/* Param 5: speed[0]. */
{
small,
slider,
"Speed 0",
"Modulation Speed: 0.1 to 5 Hz., for 0th voice.",
0, /* 0.1 Hz */
255, /* 5 Hz */
128 /* Approx. 2.5 Hz */
},
/* Param 6: depth[0]. */
{
small,
slider,
"Depth 0",
"Depth Modulation Speed: 0 to 10 msec., for 0th voice.",
0, /* 0 msec. depth */
255, /* Maximum 10 msec. depth */
128 /* Approx. 5 msec. depth */
},
/* Param 7: mod_type[0]. */
{
small,
slider,
"Modulation 0"
"Modulation Type: 0=sine (default), 1=triangle, 2=voice off, for 0th
voice.",
0, /* Sine */
2, /* Voice off */
0 /* Sine by default */
},
/* Param 8: sxcdecay[1]). */
{
small,
slider,
"Decay 1",
"Decay of 1st chorus voice.",
0,
255, /* 1 */
200 /* Approx. 80% of dynamic range */
},
/* Param 9: sxcdelay[1]. */
{
small,
slider,
"Delay 1",
"Delay: 20 to 100 msec., for 1st voice.",
0, /* Minimum 20 msec. delay */
255, /* Maximum 100 msec. delay */
128 /* Approx. 60 msec. delay */
},
/* Param 10: speed[1]. */
{
small,
slider,
"Speed 1",
"Modulation Speed: 0.1 to 5 Hz., for 1st voice.",
0, /* 0.1 Hz */
255, /* 5 Hz */
128 /* Approx. 2.5 Hz */
},
/* Param 11: depth[1]. */
{
small,
slider,
"Depth 1",
"Depth Modulation Speed: 0 to 10 msec., for 1st voice.",
0, /* 0 msec. depth */
255, /* Maximum 10 msec. depth */
128 /* Approx. 5 msec. depth */
},
/* Param 12: mod_type[1]. */
{
small,
slider,
"Modulation 1"
"Modulation Type: 0=sine (default), 1=triangle, 2=voice off, for 1st
voice.",
0, /* Sine */
2, /* Voice off */
0 /* Sine by default */
},
/* Param 13: sxcdecay[2]. */
{
small,
slider,
"Decay 2",
"Decay of 2nd chorus voice.",
0,
255, /* 1 */
200 /* Approx. 80% of dynamic range */
},
/* Param 14: sxcdelay[2]. */
{
small,
slider,
"Delay 2",
"Delay: 20 to 100 msec., for 2nd voice.",
0, /* Minimum 20 msec. delay */
255, /* Maximum 100 msec. delay */
128 /* Approx. 60 msec. delay */
},
/* Param 15: speed[2]. */
{
small,
slider,
"Speed 2",
"Modulation Speed: 0.1 to 5 Hz., for 2nd voice.",
0, /* 0.1 Hz */
255, /* 5 Hz */
128 /* Approx. 2.5 Hz */
},
/* Param 16: depth[2]. */
{
small,
slider,
"Depth 2",
"Depth Modulation Speed: 0 to 10 msec., for 2nd voice.",
0, /* 0 msec. depth */
255, /* Maximum 10 msec. depth */
128 /* Approx. 5 msec. depth */
},
/* Param 17: mod_type[2]. */
{
small,
slider,
"Modulation 2"
"Modulation Type: 0=sine (default), 1=triangle, 2=voice off, for 2nd
voice.",
0, /* Sine */
2, /* Voice off */
0 /* Sine by default */
},
};
/* State of the soxchorus. */
typedef struct {
float in_gain, out_gain;
float sxcdecay[max_chorus];
float sxcdelay[max_chorus];
float speed[max_chorus];
float depth[max_chorus];
int mod_type[max_chorus];
long depth_smpls[max_chorus];
long ttlchrssmpls[max_chorus];
int offset[max_chorus];
int phase[max_chorus];
long phz_width[max_chorus];
samp* chorusbuf;
float maxchrsdelay;
long maxttlchrssmpls;
long maxchrssmpls;
int chorus_cntr;
int chorus_voices;
double unit_vol_scale;
double exp_fract_vol;
// float clip; // Put in after clipping issue decided.
} sxc_state;
int ox_init(machine_type *t) {
t->long_name = "soxchorus (Mueller, Nelson, et al.)";
t->short_name = "soxchorus";
t->max_tracks = 1;
t->input_channels = 1;
t->output_channels = 1;
t->num_params = 16;
t->param_specs = sxc_params;
return 1;
}
void ox_create(machine *m) {
sxc_state *s;
int i;
s = (sxc_state*) malloc(sizeof(sxc_state));
s->chorus_voices = max_chorus;
s->unit_vol_scale = exp((double)DYN_RANGE / 3.0) - 1.0;
s->exp_fract_vol = ((double)DYN_RANGE / 3) / 255.0;
/* Set defaults */
// s->in_gain = 128.0/255.0; /* Approx. 0.5 (old) */
s->in_gain = ((float)exp(200.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->out_gain = 1.0; /* No attenuation */
s->maxttlchrssmpls = 0;
s->sxcdecay[0] = ((float)exp(200.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxcdelay[0] = 128.0/255.0 * 80.0 + 20.0; /* Approx. 60 msec. */
s->speed[0] = 128.0/255.0 * 4.9 + 0.1; /* Approx. 2.5 Hz */
s->depth[0] = 128.0/255.0 * 10.0; /* Approx. 5 msec. */
s->mod_type[0] = 0; /* Sine */
s->sxcdecay[1] = ((float)exp(200.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxcdelay[1] = 0.0/255.0 * 80.0 + 20.0; /* Approx. 20 msec. */
s->speed[1] = 255.0/255.0 * 4.9 + 0.1; /* 5 Hz */
s->depth[1] = 32.0/255.0 * 10.0; /* Approx. 1.7 msec. */
s->mod_type[1] = 1; /* Triangle */
s->sxcdecay[2] = ((float)exp(200.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxcdelay[2] = 255.0/255.0 * 80.0 + 20.0; /* Approx. 100 msec. */
s->speed[2] = 25.0/255.0 * 4.9 + 0.1; /* Approx. 0.5 Hz */
s->depth[2] = 255.0/255.0 * 10.0; /* Approx. 10 msec. */
s->mod_type[2] = 0; /* Sine */
for(i = 0; i < s->chorus_voices; i++) {
// s->sxcdecay[i] = 128.0/255.0; /* Approx. 0.5 (old)
*/
// s->sxcdecay[i] = ((float)exp(200.0 * s->exp_fract_vol) - 1.0)
// / s->unit_vol_scale;
// s->sxcdelay[i] = 128.0/255.0 * 80.0 + 20.0; /* Approx. 60 msec.
*/
// s->speed[i] = 128.0/255.0 * 4.9 + 0.1; /* Approx. 2.5 Hz */
// s->depth[i] = 128.0/255.0 * 10.0; /* Approx. 5 msec.
*/
// s->mod_type[i] = 0; /* Sine */
s->depth_smpls[i] = (long)(s->depth[i] * ((float)OX_SAMPLING_RATE)
/ 1000.0
+ 0.5);
s->ttlchrssmpls[i] = (long)((s->sxcdelay[i] + s->depth[i])
* ((float)OX_SAMPLING_RATE) / 1000.0 + 0.5);
if (s->ttlchrssmpls[i] > s->maxttlchrssmpls)
s->maxttlchrssmpls = s->ttlchrssmpls[i];
// s->maxttlchrssmpls = s->ttlchrssmpls[i] + 11;
if (s->mod_type[i] == 0)
s->offset[i] = (int)((s->ttlchrssmpls[i] -1) -
s->depth_smpls[i]);
else if (s->mod_type[i] == 1)
s->offset[i] = (int)((s->ttlchrssmpls[i] -1) - 2 *
s->depth_smpls[i]);
else // if (s->mod_type[i] == 2)
s->offset[i] = 0;
s->phz_width[i] = (long)(((float)OX_SAMPLING_RATE) / s->speed[i] +
0.5);
s->phase[i] = 0;
}
s->maxchrsdelay = 120.0; /* 120 msecs, only 110 msecs obtainable in
program */
s->maxchrssmpls = (long)(s->maxchrsdelay * ((float)OX_SAMPLING_RATE)
/ 1000.0 + 0.5);
// s->chorus_cntr = 1; /* Start at 1 to account for add of 1 in phase
shift below */
s->chorus_cntr = 0;
// s->clip = ?; // Put in after clipping issue decided.
s->chorusbuf = (samp*) malloc(s->maxchrssmpls * sizeof(samp));
// memset(s->chorusbuf, 0, s->maxchrssmpls * sizeof(samp)); //
Unnecessary
m->state = (void *) s;
return;
}
void ox_destroy(machine *m) {
free( ((sxc_state*)(m->state)) -> chorusbuf);
free(m->state);
m->state = NULL;
return;
}
void ox_update(machine *m) {
sxc_state *s = (sxc_state *) m->state;
param temp;
int i, y;
temp = m->params[0][ix_gain_in];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
s->in_gain = ((float)temp) / 255.0;
}
temp = m->params[0][ix_gain_out];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
s->out_gain = ((float)temp) / 255.0;
}
s->maxttlchrssmpls = 0;
for(i = 0; i < s->chorus_voices; i++) {
temp = m->params[0][ix_sxcdecay + i * 5];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
// s->sxcdecay[i] = ((float)temp) / 255.0; // Old, linear.
s->sxcdecay[i] = ((float)exp((double)temp * s->exp_fract_vol) -
1.0)
/ s->unit_vol_scale;
}
temp = m->params[0][ix_sxcdelay + i * 5];
if (temp != nochange) {
if (temp < 1) temp = 1;
else if (temp > 255) temp = 255;
s->sxcdelay[i] = ((float)temp) / 255.0 * 5.0;
s->ttlchrssmpls[i] = (long)(s->sxcdelay[i] *
((float)OX_SAMPLING_RATE)
/ 1000.0 + 0.5);
}
temp = m->params[0][ix_speed + i * 5];
if (temp != nochange) {
if (temp < 10) temp = 10;
else if (temp > 255) temp = 255;
s->speed[i] = ((float)temp) / 255.0 * 2.0;
s->phz_width[i] = (long)(((float)OX_SAMPLING_RATE) /
s->speed[i] + 0.5);
}
temp = m->params[0][ix_depth + i * 5];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
s->depth[i] = ((float)temp) / 255.0 * 10.0;
s->depth_smpls[i] = (long)(s->depth[i] *
((float)OX_SAMPLING_RATE) / 1000.0
+ 0.5);
}
temp = m->params[0][ix_mod_type + i * 5];
if (temp != nochange) {
if (temp > 2 || temp < 0) temp = 2; // If out-of-range, turn
off voice.
s->mod_type[i] = temp;
}
s->ttlchrssmpls[i] = (long)((s->sxcdelay[i] + s->depth[i])
* ((float)OX_SAMPLING_RATE) / 1000.0 + 0.5);
if (s->ttlchrssmpls[i] > s->maxttlchrssmpls)
s->maxttlchrssmpls = s->ttlchrssmpls[i];
// s->maxttlchrssmpls = s->ttlchrssmpls[i] + 11;
if (s->mod_type[i] == 0)
s->offset[i] = (int)((s->ttlchrssmpls[i] - 1) -
s->depth_smpls[i]);
else if (s->mod_type[i] == 1)
s->offset[i] = (int)((s->ttlchrssmpls[i] - 1) - 2 *
s->depth_smpls[i]);
else // if (s->mod_type[i] == 2)
s->offset[i] = 0;
}
return;
}
const char *ox_desc(int which_param, param value) {
static char temp_string[80];
float user_value;
int i, percent;
sprintf(temp_string, "ERROR");
if (which_param == ix_gain_in || which_param == ix_gain_out) {
if (value < 0) sprintf(temp_string, "Gain must not be less than
0.");
else {
user_value = ((float)value) / 255.0;
percent = (int)(user_value * 100.0 + 0.5);
if (user_value > 1.0)
sprintf(temp_string, "Gain must be not greater than 100%.");
else sprintf(temp_string, "%3.0f%%", percent);
}
}
else {
for(i = 0; i < max_chorus; i++) {
if (which_param == (ix_sxcdecay + i * 5)) {
if (value < 0) sprintf(temp_string, "%s%d%s", "Decay ", i,
" must not be less than 0.");
else {
user_value = (float)value / 255.0;
percent = (int)(user_value * 100.0 + 0.5);
if (user_value > 1.0)
sprintf(temp_string, "%s%d%s", "Decay ", i,
"must be not greater than 100%.");
else sprintf(temp_string, "%3.0f%%", percent);
}
}
else if (which_param == (ix_sxcdelay + i * 5)) {
if (value < 0) sprintf(temp_string, "%s%d%s", "Delay ", i,
" must not be less than 0 (20
msec.");
else {
user_value = (float)value / 255.0 * 80.0 + 20.0;
if (user_value > 100.0)
sprintf(temp_string, "%s%d%s", "Delay ", i,
" must be not greater than 100 msec.");
else sprintf(temp_string, "%3.0f", user_value);
}
}
else if (which_param == ix_speed + i * 5) {
if (value < 0) sprintf(temp_string, "%s%d%s", "Speed ", i,
" must not be less than 0 (0.1
Hz).");
else {
user_value = ((float) value) / 255.0 * 5.0;
if (user_value > 5.0)
sprintf(temp_string, "%s%d%s", "Speed ", i,
" must be not greater than 5 Hz.");
else sprintf(temp_string, "%3.2f", user_value);
}
}
else if (which_param == ix_depth + i * 5) {
if (value < 0) sprintf(temp_string, "%s%d%s", "Depth ", i,
" must not be less than 0 msec.");
else {
user_value = ((float) value) / 255.0 * 10.0;
if (user_value > 10.0)
sprintf(temp_string, "%s%d%s", "Depth ", i,
" must be not greater than 10 msec.");
else sprintf(temp_string, "%3.1f", user_value);
}
}
else if (which_param == ix_mod_type + i * 5) {
switch (value) {
case 0 : sprintf(temp_string, "Sine"); break;
case 1 : sprintf(temp_string, "Triangle"); break;
default : sprintf(temp_string, "%s%d%s", "Error:
Modulation", i,
"type must be 0 = sine, 1 = triangle, or 2 =
voice off.");
}
}
}
}
/* SoX: Clipping - Be nice and check the hint with warning, if...
if ( flanger->in_gain * ( 1.0 + flanger->sxcdecay ) > 1.0 /
flanger->out_gain )
warn("flanger: warning >>> gain-out can cause saturation or
clipping of output <<<");
*/
}
int ox_work(machine *m, int block_size) {
sxc_state *s = (sxc_state *) m->state;
int i, j, x_phz, ytemp;
float x_flt, ingain, outgain; // dcy, phzdepth, phz2depth, phz2width;
samp* mlout;
// samp* mrout;
samp* mlin;
samp* chrsbuf;
mlout = m->lout;
// mrout = m->rout;
mlin = m->lin;
ingain = s->in_gain;
outgain = s->out_gain;
chrsbuf = s->chorusbuf;
for(i = 0; i < block_size; i++) {
/* Compute output first */
mlout[i] = mlin[i] * ingain;
for ( j = 0; j < s->chorus_voices; j++ ) {
if (s->mod_type[j] != MOD_SILENT) { // Start no indent
if (s->mod_type[j] == MOD_SINE) {
// Get sine phase shift. 1 added to correct noise at 0.
x_flt = (float)
sin((double)s->phase[j]/(double)s->phz_width[j] * PI2);
x_phz = s->offset[j] + (int) (x_flt *
(double)s->depth_smpls[j] + 0.5);
// + 1;
}
else { // s->mod_type[j] == MOD_TRIANGLE
// Compute triangle - ramp up and down.
if (s->phase[j] < (s->phz_width[j] / 2)) { // Ramp up
x_flt = (float)s->phase[j] * 2.0 /
(float)s->phz_width[j];
x_phz = s->offset[j] + (int) (x_flt * 2.0 *
(float)s->depth_smpls[j]
+ 0.5);
}
else { // Ramp down
x_flt = (float)(s->phz_width[j] - s->phase[j]) * 2.0
/ (float)s->phz_width[j];
x_phz = s->offset[j] + (int) (x_flt * 2.0 *
(float)s->depth_smpls[j]
+ 0.5);
}
}
mlout[i] += chrsbuf[(s->maxttlchrssmpls + s->chorus_cntr
- x_phz) % s->maxttlchrssmpls] *
s->sxcdecay[j];
} // End no indent
}
/* Adjust the output volume. */
mlout[i] *= outgain;
/* Check for clipping. Not implemented. */
// if (mlout[i] > hig_clip) mlout[i] = hig_clip;
// else if (mlout[i] < low_clip) mlout[i] = low_clip;
/* Mix decay of delay and input */
chrsbuf[s->chorus_cntr] = mlin[i];
s->chorus_cntr = ( s->chorus_cntr + 1 ) % s->maxttlchrssmpls;
for ( j = 0; j < s->chorus_voices; j++ )
s->phase[j] = ( s->phase[j] + 1 ) % s->phz_width[j];
}
/* Save speedup variables to state. */
return 1;
}
void ox_track(machine *m, int change) {
}
--------------24819DF4628C2DE0E67AAF2F
Content-Type: text/html; charset=us-ascii
Content-Transfer-Encoding: 7bit
<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
After working most of the day to eliminate clicks on some settings associated
with
<br>the cycle speed of a voice and not being successful, I have released
this version
<br>so that anyone else may work on the problem. I had this similar
problem with the
<br>soxflanger.c but was able to overcome it. We could use an ability
to write the
<br>final audio output to disk in, say, the signed, 16-bit, little-endian
audio format
<br>seen in output.c.
<p>Additionally, I have made a number of changes to flanger.c but will
wait to release
<br>that code to avoid posting too many copies quickly.<tt></tt>
<p><tt>/* SOXCHORUS.C</tt>
<br><tt>*</tt>
<br><tt>* This code has been adapted from the SoX chorus.c code:</tt>
<br><tt>* August 24, 1998</tt>
<br><tt>*</tt>
<br><tt>* Copyright (C) 1998 Juergen Mueller And Sundry Contributors</tt>
<br><tt>* This source code is freely redistributable and may be used for</tt>
<br><tt>* any purpose. This copyright notice must be maintained.</tt>
<br><tt>* Juergen Mueller And Sundry Contributors are not responsible for</tt>
<br><tt>* the consequences of using this software.</tt>
<br><tt>*</tt>
<br><tt>* Conversion to Octal made by:</tt>
<br><tt>*</tt>
<br><tt>* Copyright 2000 Neil Nelson <address@hidden> June 26,
2000.</tt>
<br><tt>*</tt>
<br><tt>* This software is distributed under the terms of the</tt>
<br><tt>* GNU General Public License (GPL). Read the included file</tt>
<br><tt>* COPYING for more information.</tt>
<br><tt>*</tt>
<br><tt>*</tt>
<br><tt>* Chorus effect.</tt>
<br><tt>*</tt>
<br><tt>* Flow diagram scheme for n delays ( 1 <= n <= max_chorus
):</tt>
<br><tt>*</tt>
<br><tt>* *
gain-in
___</tt>
<br><tt>* ibuff
-----+--------------------------------------------->|
|</tt>
<br><tt>*
|
_________
| |</tt>
<br><tt>*
| |
|
* decay 1 | |</tt>
<br><tt>*
+---->| delay 1 |----------------------------->| |</tt>
<br><tt>*
|
|_________|
| |</tt>
<br><tt>*
|
/|\
| |</tt>
<br><tt>*
:
|
| |</tt>
<br><tt>*
: +-----------------+
+--------------+
| + |</tt>
<br><tt>*
: | Delay control 1 |<--| mod. speed 1 |
| |</tt>
<br><tt>*
: +-----------------+
+--------------+
| |</tt>
<br><tt>*
|
_________
| |</tt>
<br><tt>*
| |
|
* decay n | |</tt>
<br><tt>*
+---->| delay n |----------------------------->| |</tt>
<br><tt>*
|_________|
| |</tt>
<br><tt>*
/|\
|___|</tt>
<br><tt>*
|
|</tt>
<br><tt>*
+-----------------+
+--------------+
| * gain-out</tt>
<br><tt>*
| Delay control n |<--| mod. speed n
|
|</tt>
<br><tt>*
+-----------------+
+--------------+
+----->obuff</tt>
<br><tt>*</tt>
<br><tt>*</tt>
<br><tt>* The delay i is controled by a sine or triangle modulation i (
1 <= i <= n).</tt>
<br><tt>*</tt>
<br><tt>* Usage: // SoX command line format.</tt>
<br><tt>* chorus gain-in gain-out delay-1 decay-1 speed-1 depth-1
-s1|t1 [</tt>
<br><tt>* delay-2 decay-2 speed-2 depth-2
-s2|-t2 ... ]</tt>
<br><tt>*</tt>
<br><tt>* Where:</tt>
<br><tt>* gain-in, decay-1 ... decay-n : 0.0 ...
1.0
volume</tt>
<br><tt>* gain-out : 0.0 ...
volume</tt>
<br><tt>* delay-1 ... delay-n : 20.0 ... 100.0 msec</tt>
<br><tt>* speed-1 ... speed-n : 0.1 ... 5.0
Hz
modulation 1 ... n</tt>
<br><tt>* depth-1 ... depth-n : 0.0 ... 10.0
msec
modulated delay 1 ... n</tt>
<br><tt>* -s1 ... -sn : modulation by sine 1 ... n</tt>
<br><tt>* -t1 ... -tn : modulation by triangle 1 ... n</tt>
<br><tt>*</tt>
<br><tt>* Note:</tt>
<br><tt>* when decay is close to 1.0, the samples can begin
clipping and the output</tt>
<br><tt>* can saturate!</tt>
<br><tt>*</tt>
<br><tt>* Hint:</tt>
<br><tt>* 1 / out-gain < gain-in ( 1 + decay-1 + ... + decay-n
)</tt>
<br><tt>*</tt>
<br><tt>* Potential Improvements:</tt>
<br><tt>* (1) Except for the expected fixed number of sxc_params,
this program</tt>
<br><tt>* could handle an arbitrary number
of chorus voices. What trouble</tt>
<br><tt>* would it be to provide for an arbitrary
number of sxc_params?</tt>
<br><tt>*</tt>
<br><tt>*/</tt><tt></tt>
<p><tt>#include <stdio.h></tt>
<br><tt>#include <stdlib.h></tt>
<br><tt>#include <math.h></tt>
<br><tt>#include "util.h"</tt>
<br><tt>#include "machine.h"</tt><tt></tt>
<p><tt>#define PI2 (6.28318530717958465692)</tt>
<br><tt>#define MOD_SINE 0</tt>
<br><tt>#define MOD_TRIANGLE 1</tt>
<br><tt>#define MOD_SILENT 2</tt><tt></tt>
<p><tt>#define DYN_RANGE 10</tt>
<br><tt>#define max_chorus 3</tt><tt></tt>
<p><tt>/* The parameters are: */</tt>
<br><tt>enum {</tt>
<br><tt> ix_gain_in, ix_gain_out,</tt>
<br><tt> ix_sxcdecay, ix_sxcdelay, ix_speed, ix_depth,
ix_mod_type,</tt>
<br><tt> ix_sxcdecay_1, ix_sxcdelay_1, ix_speed_1, ix_depth_1,
ix_mod_type_1,</tt>
<br><tt> ix_sxcdecay_2, ix_sxcdelay_2, ix_speed_2, ix_depth_2,
ix_mod_type_2</tt>
<br><tt>} param_index;</tt><tt></tt>
<p><tt>param_spec sxc_params[] = {</tt>
<br><tt> /* Param 1: gain_in. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Gain In",</tt>
<br><tt> "Gain In",</tt>
<br><tt>
0,
/* 0 gain */</tt>
<br><tt> 255,
/* No attenuation */</tt>
<br><tt> 200
/* Approx. 80% of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 2: gain_out. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Gain Out",</tt>
<br><tt> "Gain Out",</tt>
<br><tt>
0,
/* 0 gain */</tt>
<br><tt> 255,
/* No attenuation */</tt>
<br><tt> 255</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 3: sxcdecay[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 0",</tt>
<br><tt> "Decay of 0th chorus voice.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 200
/* Approx. 80% of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 4: sxcdelay[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 0",</tt>
<br><tt> "Delay: 20 to 100 msec., for 0th voice.",</tt>
<br><tt>
0,
/* Minimum 20 msec. delay */</tt>
<br><tt> 255,
/* Maximum 100 msec. delay */</tt>
<br><tt> 128
/* Approx. 60 msec. delay */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 5: speed[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Speed 0",</tt>
<br><tt> "Modulation Speed: 0.1 to 5 Hz., for 0th voice.",</tt>
<br><tt>
0,
/* 0.1 Hz */</tt>
<br><tt> 255,
/* 5 Hz */</tt>
<br><tt> 128
/* Approx. 2.5 Hz */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 6: depth[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Depth 0",</tt>
<br><tt> "Depth Modulation Speed: 0 to 10 msec., for 0th
voice.",</tt>
<br><tt>
0,
/* 0 msec. depth */</tt>
<br><tt> 255,
/* Maximum 10 msec. depth */</tt>
<br><tt> 128
/* Approx. 5 msec. depth */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 7: mod_type[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Modulation 0"</tt>
<br><tt> "Modulation Type: 0=sine (default), 1=triangle, 2=voice
off, for 0th voice.",</tt>
<br><tt>
0,
/* Sine */</tt>
<br><tt>
2,
/* Voice off */</tt>
<br><tt>
0
/* Sine by default */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 8: sxcdecay[1]). */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 1",</tt>
<br><tt> "Decay of 1st chorus voice.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 200
/* Approx. 80% of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 9: sxcdelay[1]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 1",</tt>
<br><tt> "Delay: 20 to 100 msec., for 1st voice.",</tt>
<br><tt>
0,
/* Minimum 20 msec. delay */</tt>
<br><tt> 255,
/* Maximum 100 msec. delay */</tt>
<br><tt> 128
/* Approx. 60 msec. delay */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 10: speed[1]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Speed 1",</tt>
<br><tt> "Modulation Speed: 0.1 to 5 Hz., for 1st voice.",</tt>
<br><tt>
0,
/* 0.1 Hz */</tt>
<br><tt> 255,
/* 5 Hz */</tt>
<br><tt> 128
/* Approx. 2.5 Hz */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 11: depth[1]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Depth 1",</tt>
<br><tt> "Depth Modulation Speed: 0 to 10 msec., for 1st
voice.",</tt>
<br><tt>
0,
/* 0 msec. depth */</tt>
<br><tt> 255,
/* Maximum 10 msec. depth */</tt>
<br><tt> 128
/* Approx. 5 msec. depth */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 12: mod_type[1]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Modulation 1"</tt>
<br><tt> "Modulation Type: 0=sine (default), 1=triangle, 2=voice
off, for 1st voice.",</tt>
<br><tt>
0,
/* Sine */</tt>
<br><tt>
2,
/* Voice off */</tt>
<br><tt>
0
/* Sine by default */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 13: sxcdecay[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 2",</tt>
<br><tt> "Decay of 2nd chorus voice.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 200
/* Approx. 80% of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 14: sxcdelay[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 2",</tt>
<br><tt> "Delay: 20 to 100 msec., for 2nd voice.",</tt>
<br><tt>
0,
/* Minimum 20 msec. delay */</tt>
<br><tt> 255,
/* Maximum 100 msec. delay */</tt>
<br><tt> 128
/* Approx. 60 msec. delay */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 15: speed[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Speed 2",</tt>
<br><tt> "Modulation Speed: 0.1 to 5 Hz., for 2nd voice.",</tt>
<br><tt>
0,
/* 0.1 Hz */</tt>
<br><tt> 255,
/* 5 Hz */</tt>
<br><tt> 128
/* Approx. 2.5 Hz */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 16: depth[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Depth 2",</tt>
<br><tt> "Depth Modulation Speed: 0 to 10 msec., for 2nd
voice.",</tt>
<br><tt>
0,
/* 0 msec. depth */</tt>
<br><tt> 255,
/* Maximum 10 msec. depth */</tt>
<br><tt> 128
/* Approx. 5 msec. depth */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 17: mod_type[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Modulation 2"</tt>
<br><tt> "Modulation Type: 0=sine (default), 1=triangle, 2=voice
off, for 2nd voice.",</tt>
<br><tt>
0,
/* Sine */</tt>
<br><tt>
2,
/* Voice off */</tt>
<br><tt>
0
/* Sine by default */</tt>
<br><tt> },</tt>
<br><tt>};</tt><tt></tt>
<p><tt>/* State of the soxchorus. */</tt><tt></tt>
<p><tt>typedef struct {</tt>
<br><tt> float in_gain, out_gain;</tt>
<br><tt> float sxcdecay[max_chorus];</tt>
<br><tt> float sxcdelay[max_chorus];</tt>
<br><tt> float speed[max_chorus];</tt>
<br><tt> float depth[max_chorus];</tt>
<br><tt> int mod_type[max_chorus];</tt>
<br><tt> long depth_smpls[max_chorus];</tt>
<br><tt> long ttlchrssmpls[max_chorus];</tt>
<br><tt> int offset[max_chorus];</tt>
<br><tt> int phase[max_chorus];</tt>
<br><tt> long phz_width[max_chorus];</tt>
<br><tt> samp* chorusbuf;</tt>
<br><tt> float maxchrsdelay;</tt>
<br><tt> long maxttlchrssmpls;</tt>
<br><tt> long maxchrssmpls;</tt>
<br><tt> int chorus_cntr;</tt>
<br><tt> int chorus_voices;</tt>
<br><tt> double unit_vol_scale;</tt>
<br><tt> double exp_fract_vol;</tt>
<br><tt>// float clip; // Put in after clipping
issue decided.</tt>
<br><tt>} sxc_state;</tt><tt></tt>
<p><tt>int ox_init(machine_type *t) {</tt>
<br><tt> t->long_name = "soxchorus (Mueller, Nelson, et al.)";</tt>
<br><tt> t->short_name = "soxchorus";</tt>
<br><tt> t->max_tracks = 1;</tt>
<br><tt> t->input_channels = 1;</tt>
<br><tt> t->output_channels = 1;</tt>
<br><tt> t->num_params = 16;</tt>
<br><tt> t->param_specs = sxc_params;</tt><tt></tt>
<p><tt> return 1;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_create(machine *m) {</tt>
<br><tt> sxc_state *s;</tt>
<br><tt> int i;</tt><tt></tt>
<p><tt> s = (sxc_state*) malloc(sizeof(sxc_state));</tt><tt></tt>
<p><tt> s->chorus_voices = max_chorus;</tt><tt></tt>
<p><tt> s->unit_vol_scale = exp((double)DYN_RANGE / 3.0) -
1.0;</tt>
<br><tt> s->exp_fract_vol = ((double)DYN_RANGE / 3) /
255.0;</tt><tt></tt>
<p><tt> /* Set defaults */</tt>
<br><tt>// s->in_gain = 128.0/255.0;
/* Approx. 0.5 (old) */</tt>
<br><tt> s->in_gain = ((float)exp(200.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->out_gain =
1.0;
/* No attenuation */</tt><tt></tt>
<p><tt> s->maxttlchrssmpls = 0;</tt><tt></tt>
<p><tt> s->sxcdecay[0] = ((float)exp(200.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxcdelay[0] = 128.0/255.0 * 80.0 + 20.0; /* Approx.
60 msec. */</tt>
<br><tt> s->speed[0] = 128.0/255.0 * 4.9
+ 0.1; /* Approx. 2.5 Hz */</tt>
<br><tt> s->depth[0] = 128.0/255.0 *
10.0;
/* Approx. 5 msec. */</tt>
<br><tt> s->mod_type[0] =
0;
/* Sine */</tt><tt></tt>
<p><tt> s->sxcdecay[1] = ((float)exp(200.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxcdelay[1] = 0.0/255.0 * 80.0 + 20.0;
/* Approx. 20 msec. */</tt>
<br><tt> s->speed[1] = 255.0/255.0 * 4.9
+ 0.1; /* 5 Hz */</tt>
<br><tt> s->depth[1] = 32.0/255.0 *
10.0;
/* Approx. 1.7 msec. */</tt>
<br><tt> s->mod_type[1] =
1;
/* Triangle */</tt><tt></tt>
<p><tt> s->sxcdecay[2] = ((float)exp(200.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxcdelay[2] = 255.0/255.0 * 80.0 + 20.0; /* Approx.
100 msec. */</tt>
<br><tt> s->speed[2] = 25.0/255.0 * 4.9 +
0.1; /* Approx. 0.5 Hz */</tt>
<br><tt> s->depth[2] = 255.0/255.0 *
10.0;
/* Approx. 10 msec. */</tt>
<br><tt> s->mod_type[2] =
0;
/* Sine */</tt><tt></tt>
<p><tt> for(i = 0; i < s->chorus_voices; i++) {</tt>
<br><tt>// s->sxcdecay[i] =
128.0/255.0;
/* Approx. 0.5 (old) */</tt>
<br><tt>// s->sxcdecay[i] = ((float)exp(200.0
* s->exp_fract_vol) - 1.0)</tt>
<br><tt>//
/ s->unit_vol_scale;</tt>
<br><tt>// s->sxcdelay[i] = 128.0/255.0 *
80.0 + 20.0; /* Approx. 60 msec. */</tt>
<br><tt>// s->speed[i]
= 128.0/255.0 * 4.9 + 0.1; /* Approx. 2.5 Hz */</tt>
<br><tt>// s->depth[i]
= 128.0/255.0 * 10.0; /* Approx.
5 msec. */</tt>
<br><tt>// s->mod_type[i] =
0;
/* Sine */</tt><tt></tt>
<p><tt> s->depth_smpls[i] = (long)(s->depth[i]
* ((float)OX_SAMPLING_RATE) / 1000.0</tt>
<br><tt>
+ 0.5);</tt>
<br><tt> s->ttlchrssmpls[i] =
(long)((s->sxcdelay[i]
+ s->depth[i])</tt>
<br><tt>
* ((float)OX_SAMPLING_RATE) / 1000.0 + 0.5);</tt><tt></tt>
<p><tt> if (s->ttlchrssmpls[i] >
s->maxttlchrssmpls)</tt>
<br><tt> s->maxttlchrssmpls
= s->ttlchrssmpls[i];</tt>
<br><tt>// s->maxttlchrssmpls
= s->ttlchrssmpls[i] + 11;</tt><tt></tt>
<p><tt> if (s->mod_type[i] == 0)</tt>
<br><tt> s->offset[i] =
(int)((s->ttlchrssmpls[i] -1) - s->depth_smpls[i]);</tt>
<br><tt> else if (s->mod_type[i] == 1)</tt>
<br><tt> s->offset[i] =
(int)((s->ttlchrssmpls[i] -1) - 2 * s->depth_smpls[i]);</tt>
<br><tt> else // if (s->mod_type[i] == 2)</tt>
<br><tt> s->offset[i] =
0;</tt><tt></tt>
<p><tt> s->phz_width[i] =
(long)(((float)OX_SAMPLING_RATE)
/ s->speed[i] + 0.5);</tt>
<br><tt> s->phase[i] = 0;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> s->maxchrsdelay = 120.0; /* 120 msecs, only 110
msecs obtainable in program */</tt>
<br><tt> s->maxchrssmpls = (long)(s->maxchrsdelay *
((float)OX_SAMPLING_RATE)</tt>
<br><tt>
/ 1000.0 + 0.5);</tt><tt></tt>
<p><tt>// s->chorus_cntr = 1; /* Start at 1 to account for add of 1 in
phase shift below */</tt>
<br><tt> s->chorus_cntr = 0;</tt><tt></tt>
<p><tt>// s->clip = ?; // Put in after clipping issue
decided.</tt><tt></tt>
<p><tt> s->chorusbuf = (samp*) malloc(s->maxchrssmpls *
sizeof(samp));</tt>
<br><tt>// memset(s->chorusbuf, 0, s->maxchrssmpls * sizeof(samp));
// Unnecessary</tt><tt></tt>
<p><tt> m->state = (void *) s;</tt>
<br><tt> return;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_destroy(machine *m) {</tt>
<br><tt> free( ((sxc_state*)(m->state)) -> chorusbuf);</tt>
<br><tt> free(m->state);</tt>
<br><tt> m->state = NULL;</tt>
<br><tt> return;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_update(machine *m) {</tt>
<br><tt> sxc_state *s = (sxc_state *) m->state;</tt>
<br><tt> param temp;</tt>
<br><tt> int i, y;</tt><tt></tt>
<p><tt> temp = m->params[0][ix_gain_in];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt> if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp > 255) temp = 255;</tt>
<br><tt> s->in_gain = ((float)temp) / 255.0;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_gain_out];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt> if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp > 255) temp = 255;</tt>
<br><tt> s->out_gain = ((float)temp) / 255.0;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> s->maxttlchrssmpls = 0;</tt><tt></tt>
<p><tt> for(i = 0; i < s->chorus_voices; i++) {</tt>
<br><tt> temp = m->params[0][ix_sxcdecay
+ i * 5];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt>
if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp
> 255) temp = 255;</tt><tt></tt>
<p><tt>// s->sxcdecay[i]
= ((float)temp) / 255.0; // Old, linear.</tt>
<br><tt> s->sxcdecay[i]
= ((float)exp((double)temp * s->exp_fract_vol) - 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_sxcdelay +
i * 5];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt>
if
(temp < 1) temp = 1;</tt>
<br><tt> else if (temp
> 255) temp = 255;</tt><tt></tt>
<p><tt> s->sxcdelay[i]
= ((float)temp) / 255.0 * 5.0;</tt>
<br><tt> s->ttlchrssmpls[i]
= (long)(s->sxcdelay[i] * ((float)OX_SAMPLING_RATE)</tt>
<br><tt>
/ 1000.0 + 0.5);</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_speed + i
* 5];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt>
if
(temp < 10) temp = 10;</tt>
<br><tt> else if (temp
> 255) temp = 255;</tt><tt></tt>
<p><tt> s->speed[i] =
((float)temp)
/ 255.0 * 2.0;</tt>
<br><tt> s->phz_width[i]
= (long)(((float)OX_SAMPLING_RATE) / s->speed[i] + 0.5);</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_depth + i
* 5];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt>
if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp
> 255) temp = 255;</tt><tt></tt>
<p><tt> s->depth[i] =
((float)temp)
/ 255.0 * 10.0;</tt>
<br><tt> s->depth_smpls[i]
= (long)(s->depth[i] * ((float)OX_SAMPLING_RATE) / 1000.0</tt>
<br><tt>
+ 0.5);</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_mod_type +
i * 5];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt> if (temp > 2 ||
temp < 0) temp = 2; // If out-of-range, turn off voice.</tt>
<br><tt> s->mod_type[i]
= temp;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> s->ttlchrssmpls[i] =
(long)((s->sxcdelay[i]
+ s->depth[i])</tt>
<br><tt>
* ((float)OX_SAMPLING_RATE) / 1000.0 + 0.5);</tt><tt></tt>
<p><tt> if (s->ttlchrssmpls[i] >
s->maxttlchrssmpls)</tt>
<br><tt> s->maxttlchrssmpls
= s->ttlchrssmpls[i];</tt>
<br><tt>// s->maxttlchrssmpls
= s->ttlchrssmpls[i] + 11;</tt><tt></tt>
<p><tt> if (s->mod_type[i] == 0)</tt>
<br><tt> s->offset[i] =
(int)((s->ttlchrssmpls[i] - 1) - s->depth_smpls[i]);</tt>
<br><tt> else if (s->mod_type[i] == 1)</tt>
<br><tt> s->offset[i] =
(int)((s->ttlchrssmpls[i] - 1) - 2 * s->depth_smpls[i]);</tt>
<br><tt> else // if (s->mod_type[i] == 2)</tt>
<br><tt> s->offset[i] =
0;</tt><tt></tt>
<p><tt> }</tt><tt></tt>
<p><tt> return;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>const char *ox_desc(int which_param, param value) {</tt>
<br><tt> static char temp_string[80];</tt>
<br><tt> float user_value;</tt>
<br><tt> int i, percent;</tt><tt></tt>
<p><tt> sprintf(temp_string, "ERROR");</tt><tt></tt>
<p><tt> if (which_param == ix_gain_in || which_param == ix_gain_out)
{</tt>
<br><tt> if (value < 0) sprintf(temp_string,
"Gain must not be less than 0.");</tt>
<br><tt> else {</tt>
<br><tt> user_value =
((float)value)
/ 255.0;</tt>
<br><tt> percent =
(int)(user_value
* 100.0 + 0.5);</tt>
<br><tt> if (user_value
> 1.0)</tt>
<br><tt>
sprintf(temp_string, "Gain must be not greater than 100%.");</tt><tt></tt>
<p><tt> else
sprintf(temp_string,
"%3.0f%%", percent);</tt>
<br><tt> }</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else {</tt><tt></tt>
<p><tt> for(i = 0; i < max_chorus; i++)
{</tt>
<br><tt> if (which_param
== (ix_sxcdecay + i * 5)) {</tt>
<br><tt>
if (value < 0) sprintf(temp_string, "%s%d%s", "Decay ", i,</tt>
<br><tt>
" must not be less than 0.");</tt>
<br><tt>
else {</tt>
<br><tt>
user_value = (float)value / 255.0;</tt>
<br><tt>
percent = (int)(user_value * 100.0 + 0.5);</tt>
<br><tt>
if (user_value > 1.0)</tt>
<br><tt>
sprintf(temp_string, "%s%d%s", "Decay ", i,</tt>
<br><tt>
"must be not greater than 100%.");</tt><tt></tt>
<p><tt>
else sprintf(temp_string, "%3.0f%%", percent);</tt>
<br><tt>
}</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else if (which_param
== (ix_sxcdelay + i * 5)) {</tt>
<br><tt>
if (value < 0) sprintf(temp_string, "%s%d%s", "Delay ", i,</tt>
<br><tt>
" must not be less than 0 (20 msec.");</tt>
<br><tt>
else {</tt>
<br><tt>
user_value = (float)value / 255.0 * 80.0 + 20.0;</tt>
<br><tt>
if (user_value > 100.0)</tt>
<br><tt>
sprintf(temp_string, "%s%d%s", "Delay ", i,</tt>
<br><tt>
" must be not greater than 100 msec.");</tt><tt></tt>
<p><tt>
else sprintf(temp_string, "%3.0f", user_value);</tt>
<br><tt>
}</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else if (which_param
== ix_speed + i * 5) {</tt>
<br><tt>
if (value < 0) sprintf(temp_string, "%s%d%s", "Speed ", i,</tt>
<br><tt>
" must not be less than 0 (0.1 Hz).");</tt>
<br><tt>
else {</tt>
<br><tt>
user_value = ((float) value) / 255.0 * 5.0;</tt>
<br><tt>
if (user_value > 5.0)</tt>
<br><tt>
sprintf(temp_string, "%s%d%s", "Speed ", i,</tt>
<br><tt>
" must be not greater than 5 Hz.");</tt><tt></tt>
<p><tt>
else sprintf(temp_string, "%3.2f", user_value);</tt>
<br><tt>
}</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else if (which_param
== ix_depth + i * 5) {</tt>
<br><tt>
if (value < 0) sprintf(temp_string, "%s%d%s", "Depth ", i,</tt>
<br><tt>
" must not be less than 0 msec.");</tt>
<br><tt>
else {</tt>
<br><tt>
user_value = ((float) value) / 255.0 * 10.0;</tt>
<br><tt>
if (user_value > 10.0)</tt>
<br><tt>
sprintf(temp_string, "%s%d%s", "Depth ", i,</tt>
<br><tt>
" must be not greater than 10 msec.");</tt><tt></tt>
<p><tt>
else sprintf(temp_string, "%3.1f", user_value);</tt>
<br><tt>
}</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else if (which_param
== ix_mod_type + i * 5) {</tt>
<br><tt>
switch (value) {</tt>
<br><tt>
case 0 : sprintf(temp_string, "Sine"); break;</tt>
<br><tt>
case 1 : sprintf(temp_string, "Triangle"); break;</tt>
<br><tt>
default : sprintf(temp_string, "%s%d%s", "Error: Modulation", i,</tt>
<br><tt>
"type must be 0 = sine, 1 = triangle, or 2 = voice off.");</tt>
<br><tt>
}</tt>
<br><tt> }</tt>
<br><tt> }</tt>
<br><tt> }</tt><tt></tt>
<p><tt>/* SoX: Clipping - Be nice and check the hint with warning, if...</tt>
<br><tt> if ( flanger->in_gain * ( 1.0 + flanger->sxcdecay
) > 1.0 / flanger->out_gain )</tt>
<br><tt> warn("flanger: warning >>> gain-out
can cause saturation or clipping of output <<<");</tt>
<br><tt>*/</tt><tt></tt>
<p><tt>}</tt><tt></tt>
<p><tt>int ox_work(machine *m, int block_size) {</tt>
<br><tt> sxc_state *s = (sxc_state *) m->state;</tt>
<br><tt> int i, j, x_phz, ytemp;</tt>
<br><tt> float x_flt, ingain, outgain; // dcy, phzdepth, phz2depth,
phz2width;</tt><tt></tt>
<p><tt> samp* mlout;</tt>
<br><tt>// samp* mrout;</tt>
<br><tt> samp* mlin;</tt>
<br><tt> samp* chrsbuf;</tt><tt></tt>
<p><tt> mlout = m->lout;</tt>
<br><tt>// mrout = m->rout;</tt>
<br><tt> mlin = m->lin;</tt><tt></tt>
<p><tt> ingain = s->in_gain;</tt>
<br><tt> outgain = s->out_gain;</tt>
<br><tt> chrsbuf = s->chorusbuf;</tt><tt></tt>
<p><tt> for(i = 0; i < block_size; i++) {</tt><tt></tt>
<p><tt> /* Compute output first */</tt>
<br><tt> mlout[i] = mlin[i] *
ingain;</tt><tt></tt>
<p><tt> for ( j = 0; j < s->chorus_voices;
j++ ) {</tt><tt></tt>
<p><tt> if (s->mod_type[j]
!= MOD_SILENT) { // Start no indent</tt>
<br><tt> if (s->mod_type[j]
== MOD_SINE) {</tt><tt></tt>
<p><tt>
// Get sine phase shift. 1 added to correct noise at 0.</tt>
<br><tt>
x_flt = (float) sin((double)s->phase[j]/(double)s->phz_width[j] * PI2);</tt>
<br><tt>
x_phz = s->offset[j] + (int) (x_flt * (double)s->depth_smpls[j] + 0.5);</tt>
<br><tt>//
+ 1;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else { // s->mod_type[j]
== MOD_TRIANGLE</tt><tt></tt>
<p><tt>
// Compute triangle - ramp up and down.</tt>
<br><tt>
if (s->phase[j] < (s->phz_width[j] / 2)) { // Ramp up</tt>
<br><tt>
x_flt = (float)s->phase[j] * 2.0 / (float)s->phz_width[j];</tt>
<br><tt>
x_phz = s->offset[j] + (int) (x_flt * 2.0 * (float)s->depth_smpls[j]</tt>
<br><tt>
+ 0.5);</tt>
<br><tt>
}</tt><tt></tt>
<p><tt>
else { // Ramp down</tt>
<br><tt>
x_flt = (float)(s->phz_width[j] - s->phase[j]) * 2.0</tt>
<br><tt>
/ (float)s->phz_width[j];</tt>
<br><tt>
x_phz = s->offset[j] + (int) (x_flt * 2.0 * (float)s->depth_smpls[j]</tt>
<br><tt>
+ 0.5);</tt>
<br><tt>
}</tt>
<br><tt> }</tt><tt></tt>
<p><tt> mlout[i] +=
chrsbuf[(s->maxttlchrssmpls
+ s->chorus_cntr</tt>
<br><tt>
- x_phz) % s->maxttlchrssmpls] * s->sxcdecay[j];</tt>
<br><tt> } // End no indent</tt>
<br><tt> }</tt><tt></tt>
<p><tt> /* Adjust the output volume. */</tt>
<br><tt> mlout[i] *= outgain;</tt><tt></tt>
<p><tt> /* Check for clipping. Not
implemented. */</tt>
<br><tt>// if
(mlout[i] > hig_clip) mlout[i] = hig_clip;</tt>
<br><tt>// else if (mlout[i] < low_clip)
mlout[i] = low_clip;</tt><tt></tt>
<p><tt> /* Mix decay of delay and input */</tt>
<br><tt> chrsbuf[s->chorus_cntr] =
mlin[i];</tt><tt></tt>
<p><tt> s->chorus_cntr = ( s->chorus_cntr
+ 1 ) % s->maxttlchrssmpls;</tt><tt></tt>
<p><tt> for ( j = 0; j < s->chorus_voices;
j++ )</tt>
<br><tt> s->phase[j] =
( s->phase[j] + 1 ) % s->phz_width[j];</tt>
<br><tt> }</tt><tt></tt>
<p><tt> /* Save speedup variables to state. */</tt><tt></tt>
<p><tt> return 1;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_track(machine *m, int change) {</tt>
<br><tt>}</tt>
<br><tt></tt> </html>
--------------24819DF4628C2DE0E67AAF2F--
Received: from mta6.snfc21.pbi.net (mta6.snfc21.pbi.net [206.13.28.240])
by mescaline.gnu.org (8.9.1a/8.9.1) with ESMTP id OAA02904
for <address@hidden>; Tue, 27 Jun 2000 14:21:37 -0400
Received: from pacbell.net ([207.214.184.75])
by mta6.snfc21.pbi.net (Sun Internet Mail Server sims.3.5.2000.01.05.12.18.p9)
with ESMTP id <address@hidden> for address@hidden; Tue,
27 Jun 2000 10:24:28 -0700 (PDT)
Date: Tue, 27 Jun 2000 10:22:33 -0700
From: Neil Nelson <address@hidden>
To: address@hidden
Message-id: <address@hidden>
MIME-version: 1.0
X-Mailer: Mozilla 4.72 [en] (X11; U; Linux 2.2.14-5.0 i686)
Content-type: multipart/alternative;
boundary="------------801DEC9FA1F1642AE599CEF4"
X-Accept-Language: en
Subject: [Octal-dev] soxecho.c - SoX echo.c converted to Octal
Sender: address@hidden
Errors-To: address@hidden
X-BeenThere: address@hidden
X-Mailman-Version: 1.2 (beta 1)
Precedence: bulk
Reply-To: address@hidden
List-Id: OCTAL tracker http://www.gnu.org/software/octal <octal-dev.gnu.org>
--------------801DEC9FA1F1642AE599CEF4
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit
No clicks in this one. Contains log volume (gain) controls.
/* SOXECHO.C
*
* This code has been adapted from the SoX echo.c code:
*
* August 24, 1998
* Copyright (C) 1998 Juergen Mueller And Sundry Contributors
* This source code is freely redistributable and may be used for
* any purpose. This copyright notice must be maintained.
* Juergen Mueller And Sundry Contributors are not responsible for
* the consequences of using this software.
*
* Conversion to Octal made by:
*
* Copyright 2000 Neil Nelson <address@hidden> June 27, 2000.
*
* This software is distributed under the terms of the
* GNU General Public License (GPL). Read the included file
* COPYING for more information.
*
*
* From SoX:
* This is the "echo.c" while the old "echo.c" from version 12 moves to
* "reverb.c" satisfying the defintions made in the Guitar FX FAQ.
*
*
* Echo effect for dsp.
*
* Flow diagram scheme for n delays ( 1 <= n <= MAX_ECHOS ):
*
* * gain-in ___
* ibuff -----------+------------------------------------------>| |
* | _________ | |
* | | | * decay 1 | |
* +----->| delay 1 |------------------------->| |
* | |_________| | |
* | _________ | + |
* | | | * decay 2 | |
* +---------->| delay 2 |-------------------->| |
* | |_________| | |
* : _________ | |
* | | | * decay n | |
* +--------------->| delay n |--------------->|___|
* |_________| |
* | *
gain-out
* |
*
+----->obuff
*
* Usage:
* echo gain-in gain-out delay-1 decay-1 [delay-2 decay-2 ... delay-n
decay-n]
*
* Where:
* gain-in, decay-1 ... decay-n : 0.0 ... 1.0 volume
* gain-out : 0.0 ... volume
* delay-1 ... delay-n : > 0.0 msec
*
* Note:
* when decay is close to 1.0, the samples can begin clipping and the
output
* can saturate!
*
* Hint:
* 1 / out-gain > gain-in ( 1 + decay-1 + ... + decay-n )
*
*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "util.h"
#include "machine.h"
#define DYN_RANGE 10
#define MAX_ECHOS 7
#define MAX_DELAY 500 // msec
// #define DELAY_LOG_RANGE 2 // Higher number skews toward selecting
shorter delays.
/* The parameters are: gain_in, gain_out, decay, delay, ... 6 more
times */
enum {
ix_gain_in, ix_gain_out,
ix_sxedecay, ix_sxedelay,
ix_sxedecay_1, ix_sxedelay_1,
ix_sxedecay_2, ix_sxedelay_2,
ix_sxedecay_3, ix_sxedelay_3,
ix_sxedecay_4, ix_sxedelay_4,
ix_sxedecay_5, ix_sxedelay_5,
ix_sxedecay_6, ix_sxedelay_6
} param_index;
param_spec sxe_params[] = {
/* Param 1: gain_in. */
{
small,
slider,
"Gain In",
"Gain In",
0, /* 0 gain */
255, /* No attenuation */
70 /* Approx. 70/255 of dynamic range */
},
/* Param 2: gain_out. */
{
small,
slider,
"Gain Out",
"Gain Out",
0, /* 0 gain */
255, /* No attenuation */
100 /* Approx. 100/255 of dynamic range */
},
/* Param 3: sxedecay[0]. */
{
small,
slider,
"Decay 0",
"Decay of Delay 0.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 4: sxedelay[0]. */
{
small,
slider,
"Delay 0",
"Delay 0 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
35 /* 35/255 of MAX_DELAY msec */
},
/* Param 5: sxedecay[1]. */
{
small,
slider,
"Decay 1",
"Decay of Delay 1.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 6: sxedelay[1]. */
{
small,
slider,
"Delay 1",
"Delay 1 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
70 /* 70/255 of MAX_DELAY msec */
},
/* Param 7: sxedecay[2]. */
{
small,
slider,
"Decay 2",
"Decay of Delay 2.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 8: sxedelay[2]. */
{
small,
slider,
"Delay 2",
"Delay 2 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
105 /* 105/255 of MAX_DELAY msec */
},
/* Param 9: sxedecay[3]. */
{
small,
slider,
"Decay 3",
"Decay of Delay 3.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 10: sxedelay[3]. */
{
small,
slider,
"Delay 3",
"Delay 3 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
140 /* 140/255 of MAX_DELAY msec */
},
/* Param 11: sxedecay[4]. */
{
small,
slider,
"Decay 4",
"Decay of Delay 4.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 12: sxedelay[4]. */
{
small,
slider,
"Delay 4",
"Delay 4 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
175 /* 175/255 of MAX_DELAY msec */
},
/* Param 13: sxedecay[5]. */
{
small,
slider,
"Decay 5",
"Decay of Delay 5.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 14: sxedelay[5]. */
{
small,
slider,
"Delay 5",
"Delay 5 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
210 /* 210/255 of MAX_DELAY msec */
},
/* Param 15: sxedecay[6]. */
{
small,
slider,
"Decay 6",
"Decay of Delay 6.",
0,
255, /* 1 */
220 /* 220/255 of dynamic range */
},
/* Param 16: sxedelay[6]. */
{
small,
slider,
"Delay 6",
"Delay 6 time: 0 to maximum.",
0, /* 0 msec. delay */
255, /* Highest delay of range */
255
},
};
/* State of the soxflanger. */
typedef struct {
float in_gain, out_gain;
float sxedecay[MAX_ECHOS];
float sxedelay[MAX_ECHOS];
long sxesamples[MAX_ECHOS];
samp* delaybuf;
long maxsxesamples;
float maxanydelay;
long maxanysmpls;
int echo_cntr;
int num_delays;
double unit_vol_scale;
double exp_fract_vol;
// double unit_delay_scale;
// double exp_fract_delay;
// float clip; // Put in after clipping issue decided.
} sxe_state;
int ox_init(machine_type *t) {
t->long_name = "soxecho (Mueller, Nelson, et al.)";
t->short_name = "soxecho";
t->max_tracks = 1;
t->input_channels = 1;
t->output_channels = 1;
t->num_params = 15;
t->param_specs = sxe_params;
return 1;
}
void ox_create(machine *m) {
sxe_state *s;
int i;
s = (sxe_state*) malloc(sizeof(sxe_state));
m->state = (void *) s;
s->num_delays = MAX_ECHOS;
s->unit_vol_scale = exp((double)DYN_RANGE / 3.0) - 1.0;
s->exp_fract_vol = ((double)DYN_RANGE / 3) / 255.0;
// s->unit_delay_scale = (exp((double)DELAY_LOG_RANGE) - 1.0) /
(double)MAX_DELAY;
// s->exp_fract_delay = ((double)DELAY_LOG_RANGE) / 255.0;
/* Set defaults */
// s->in_gain = 128.0/255.0; /* Approx. 0.5 (old) */
s->in_gain = ((float)exp(70.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
// s->out_gain = 1.0; /* No attenuation */
s->out_gain = ((float)exp(100.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedecay[0] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[0] = 35.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[0] = ((float)exp(35.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxedecay[1] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[1] = 70.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[1] = ((float)exp(70.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxedecay[2] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[2] = 105.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[2] = ((float)exp(105.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxedecay[3] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[3] = 140.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[3] = ((float)exp(140.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxedecay[4] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[4] = 175.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[4] = ((float)exp(175.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxedecay[5] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[5] = 210.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[5] = ((float)exp(210.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxedecay[6] = ((float)exp(220.0 * s->exp_fract_vol) - 1.0)
/ s->unit_vol_scale;
s->sxedelay[6] = 255.0/255.0 * (float)MAX_DELAY;
// s->sxedelay[6] = ((float)exp(255.0 * s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->maxsxesamples = 0; // Accumulated in following loop.
for(i = 0; i < s->num_delays; i++) {
s->sxesamples[i] = (long)(s->sxedelay[i] *
((float)OX_SAMPLING_RATE) / 1000.0
+ 0.5);
if (s->sxesamples[i] > s->maxsxesamples)
s->maxsxesamples = s->sxesamples[i];
}
s->maxanydelay = MAX_DELAY;
s->maxanysmpls = (long)(s->maxanydelay * ((float)OX_SAMPLING_RATE) /
1000.0
+ 10.0); // Add an extra 10 to avoid possible exp
overflow.
if (s->maxsxesamples > s->maxanysmpls)
printf("%s\n",
"soxecho create error: delay sample requirement exceeds
allocation.");
s->echo_cntr = 0;
// s->clip = ?; // Put in after clipping issue decided.
s->delaybuf = (samp*) malloc(s->maxanysmpls * sizeof(samp));
// memset(s->delaybuf, 0, s->maxanysmpls * sizeof(samp)); //
Unnecessary
// for (i=0, i < s->maxanysmpls, i++) s->delaybuf[i] = 0; // Speed
this up.
return;
}
void ox_destroy(machine *m) {
free( ((sxe_state*)(m->state)) -> delaybuf);
free(m->state);
m->state = NULL;
return;
}
void ox_update(machine *m) {
sxe_state *s = (sxe_state *) m->state;
param temp;
int i;
temp = m->params[0][ix_gain_in];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
s->in_gain = ((float)temp) / 255.0;
}
temp = m->params[0][ix_gain_out];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
s->out_gain = ((float)temp) / 255.0;
}
s->maxsxesamples = 0; // Accumulated in following loop.
for(i = 0; i < s->num_delays; i++) {
temp = m->params[0][ix_sxedecay + i * 2];
if (temp != nochange) {
if (temp < 0) temp = 0;
else if (temp > 255) temp = 255;
// s->sxedecay[i] = ((float)temp) / 255.0; // Old, linear.
s->sxedecay[i] = ((float)exp((double)temp * s->exp_fract_vol) -
1.0)
/ s->unit_vol_scale;
}
temp = m->params[0][ix_sxedelay + i * 2];
if (temp != nochange) {
if (temp < 1) temp = 0;
else if (temp > 255) temp = 255;
s->sxedelay[i] = (float)temp/255.0 * (float)MAX_DELAY;
// s->sxedelay[i] = ((float)exp((double)temp *
s->exp_fract_delay) - 1.0)
// / s->unit_delay_scale;
s->sxesamples[i] = (long)(s->sxedelay[i] *
((float)OX_SAMPLING_RATE)
/ 1000.0 + 0.5);
}
if (s->sxesamples[i] > s->maxsxesamples) s->maxsxesamples =
s->sxesamples[i];
}
return;
}
const char *ox_desc(int which_param, param value) {
static char temp_string[80];
float user_value;
int i, percent;
sprintf(temp_string, "ERROR");
if (which_param == ix_gain_in || which_param == ix_gain_out) {
if (value < 0) sprintf(temp_string, "Gain must not be less than
0.");
else {
user_value = ((float)value) / 255.0;
percent = (int)(user_value * 100.0 + 0.5);
if (user_value > 1.0)
sprintf(temp_string, "Gain must be not greater than 100%.");
else sprintf(temp_string, "%3.0f%%", percent);
}
}
else {
for(i = 0; i < MAX_ECHOS; i++) {
if (which_param == (ix_sxedecay + i * 2)) {
if (value < 0) sprintf(temp_string, "%s%d%s", "Decay ", i,
" must not be less than 0. 0 removes this
delay");
else {
user_value = (float)value / 255.0;
percent = (int)(user_value * 100.0 + 0.5);
if (user_value > 1.0)
sprintf(temp_string, "%s%d%s", "Decay ", i,
"must be not greater than 100%.");
else sprintf(temp_string, "%3.0f%%", percent);
}
}
else if (which_param == (ix_sxedelay + i * 2)) {
if (value < 0) sprintf(temp_string, "%s%d%s", "Delay ", i,
" must not be less than 0.");
else {
user_value = (float)value / 255.0 * (float)MAX_DELAY;
if (user_value > (float)MAX_DELAY)
sprintf(temp_string, "%s%d%s%d%s", "Delay ", i,
" must be not greater than ", MAX_DELAY, "
msec.");
else sprintf(temp_string, "%3.0f", user_value);
}
}
}
}
/*
*** Be nice and check the hint with warning, if...
sum_in_volume = 1.0;
for ( i = 0; i < echo->num_delays; i++ )
sum_in_volume += echo->decay[i];
if ( sum_in_volume * echo->in_gain > 1.0 / echo->out_gain )
warn("echo: warning >>> gain-out can cause saturation of output <<<");
*/
}
int ox_work(machine *m, int block_size) {
sxe_state *s = (sxe_state *) m->state;
samp* mlout;
// samp* mrout;
samp* mlin;
samp* dlybf;
int i, j, dlyctr, numdlys;
long mxsmpls;
long* sxsmpls;
float ingain, outgain;
float* sxdecy;
// float hig_clip, low_clip;
// hig_clip = s->clip;
// low_clip = 0.0 - s->clip;
/* Speed-up variables */
mlout = m->lout;
// mrout = m->rout;
mlin = m->lin;
ingain = s->in_gain;
outgain = s->out_gain;
dlybf = s->delaybuf;
dlyctr = s->echo_cntr;
mxsmpls = s->maxsxesamples;
sxsmpls = (long*)s->sxesamples;
sxdecy = (float*)s->sxedecay;
numdlys = s->num_delays;
for(i = 0; i < block_size; i++) {
/* Compute output first */
mlout[i] = mlin[i] * ingain;
for (j = 0; j < numdlys; j++) {
if (sxdecy[j])
mlout[i] += dlybf[(mxsmpls + dlyctr - sxsmpls[j]) % mxsmpls]
* sxdecy[j];
}
/* Adjust the output volume */
mlout[i] *= outgain;
/* Check for clipping. Not implemented. */
// if (mlout[i] > hig_clip) mlout[i] = hig_clip;
// else if (mlout[i] < low_clip) mlout[i] = low_clip;
/* Store input in delay buffer */
dlybf[dlyctr] = mlin[i];
/* Adjust the counter */
dlyctr = ( dlyctr + 1 ) % mxsmpls;
}
/* Save speed-up variables to state. */
s->echo_cntr = dlyctr;
return 1;
}
void ox_track(machine *m, int change) {
}
--------------801DEC9FA1F1642AE599CEF4
Content-Type: text/html; charset=us-ascii
Content-Transfer-Encoding: 7bit
<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
No clicks in this one. Contains log volume (gain) controls.<tt></tt>
<p><tt>/* SOXECHO.C</tt>
<br><tt>*</tt>
<br><tt>* This code has been adapted from the SoX echo.c code:</tt>
<br><tt>*</tt>
<br><tt>* August 24, 1998</tt>
<br><tt>* Copyright (C) 1998 Juergen Mueller And Sundry Contributors</tt>
<br><tt>* This source code is freely redistributable and may be used for</tt>
<br><tt>* any purpose. This copyright notice must be maintained.</tt>
<br><tt>* Juergen Mueller And Sundry Contributors are not responsible for</tt>
<br><tt>* the consequences of using this software.</tt>
<br><tt>*</tt>
<br><tt>* Conversion to Octal made by:</tt>
<br><tt>*</tt>
<br><tt>* Copyright 2000 Neil Nelson <address@hidden> June 27,
2000.</tt>
<br><tt>*</tt>
<br><tt>* This software is distributed under the terms of the</tt>
<br><tt>* GNU General Public License (GPL). Read the included file</tt>
<br><tt>* COPYING for more information.</tt>
<br><tt>*</tt>
<br><tt>*</tt>
<br><tt>* From SoX:</tt>
<br><tt>* This is the "echo.c" while the old "echo.c" from version 12 moves
to</tt>
<br><tt>* "reverb.c" satisfying the defintions made in the Guitar FX FAQ.</tt>
<br><tt>*</tt>
<br><tt>*</tt>
<br><tt>* Echo effect for dsp.</tt>
<br><tt>*</tt>
<br><tt>* Flow diagram scheme for n delays ( 1 <= n <= MAX_ECHOS
):</tt>
<br><tt>*</tt>
<br><tt>* *
gain-in
___</tt>
<br><tt>* ibuff
-----------+------------------------------------------>|
|</tt>
<br><tt>*
|
_________
| |</tt>
<br><tt>*
|
|
|
* decay 1 | |</tt>
<br><tt>*
+----->| delay 1 |------------------------->| |</tt>
<br><tt>*
|
|_________|
| |</tt>
<br><tt>*
|
_________
| + |</tt>
<br><tt>*
|
|
| * decay 2
| |</tt>
<br><tt>*
+---------->| delay 2 |-------------------->| |</tt>
<br><tt>*
|
|_________|
| |</tt>
<br><tt>*
:
_________
| |</tt>
<br><tt>*
|
|
|
* decay n | |</tt>
<br><tt>*
+--------------->| delay n |--------------->|___|</tt>
<br><tt>*
|_________|
|</tt>
<br><tt>*
| * gain-out</tt>
<br><tt>*
|</tt>
<br><tt>*
+----->obuff</tt>
<br><tt>*</tt>
<br><tt>* Usage:</tt>
<br><tt>* echo gain-in gain-out delay-1 decay-1 [delay-2 decay-2
... delay-n decay-n]</tt>
<br><tt>*</tt>
<br><tt>* Where:</tt>
<br><tt>* gain-in, decay-1 ... decay-n : 0.0 ...
1.0
volume</tt>
<br><tt>* gain-out : 0.0 ...
volume</tt>
<br><tt>* delay-1 ... delay-n : > 0.0 msec</tt>
<br><tt>*</tt>
<br><tt>* Note:</tt>
<br><tt>* when decay is close to 1.0, the samples can begin
clipping and the output</tt>
<br><tt>* can saturate!</tt>
<br><tt>*</tt>
<br><tt>* Hint:</tt>
<br><tt>* 1 / out-gain > gain-in ( 1 + decay-1 + ... + decay-n
)</tt>
<br><tt>*</tt>
<br><tt>*/</tt><tt></tt>
<p><tt>#include <stdio.h></tt>
<br><tt>#include <stdlib.h></tt>
<br><tt>#include <math.h></tt>
<br><tt>#include "util.h"</tt>
<br><tt>#include "machine.h"</tt><tt></tt>
<p><tt>#define DYN_RANGE 10</tt>
<br><tt>#define MAX_ECHOS 7</tt>
<br><tt>#define MAX_DELAY 500 // msec</tt>
<br><tt>// #define DELAY_LOG_RANGE 2 // Higher number skews toward selecting
shorter delays.</tt><tt></tt>
<p><tt>/* The parameters are: gain_in, gain_out, decay, delay,
... 6 more times */</tt>
<br><tt>enum {</tt>
<br><tt> ix_gain_in, ix_gain_out,</tt>
<br><tt> ix_sxedecay, ix_sxedelay,</tt>
<br><tt> ix_sxedecay_1, ix_sxedelay_1,</tt>
<br><tt> ix_sxedecay_2, ix_sxedelay_2,</tt>
<br><tt> ix_sxedecay_3, ix_sxedelay_3,</tt>
<br><tt> ix_sxedecay_4, ix_sxedelay_4,</tt>
<br><tt> ix_sxedecay_5, ix_sxedelay_5,</tt>
<br><tt> ix_sxedecay_6, ix_sxedelay_6</tt>
<br><tt>} param_index;</tt><tt></tt>
<p><tt>param_spec sxe_params[] = {</tt>
<br><tt> /* Param 1: gain_in. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Gain In",</tt>
<br><tt> "Gain In",</tt>
<br><tt>
0,
/* 0 gain */</tt>
<br><tt> 255,
/* No attenuation */</tt>
<br><tt>
70
/* Approx. 70/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 2: gain_out. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Gain Out",</tt>
<br><tt> "Gain Out",</tt>
<br><tt>
0,
/* 0 gain */</tt>
<br><tt> 255,
/* No attenuation */</tt>
<br><tt> 100
/* Approx. 100/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 3: sxedecay[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 0",</tt>
<br><tt> "Decay of Delay 0.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 4: sxedelay[0]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 0",</tt>
<br><tt> "Delay 0 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt>
35
/* 35/255 of MAX_DELAY msec */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 5: sxedecay[1]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 1",</tt>
<br><tt> "Decay of Delay 1.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 6: sxedelay[1]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 1",</tt>
<br><tt> "Delay 1 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt>
70
/* 70/255 of MAX_DELAY msec */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 7: sxedecay[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 2",</tt>
<br><tt> "Decay of Delay 2.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 8: sxedelay[2]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 2",</tt>
<br><tt> "Delay 2 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt> 105
/* 105/255 of MAX_DELAY msec */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 9: sxedecay[3]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 3",</tt>
<br><tt> "Decay of Delay 3.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 10: sxedelay[3]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 3",</tt>
<br><tt> "Delay 3 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt> 140
/* 140/255 of MAX_DELAY msec */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 11: sxedecay[4]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 4",</tt>
<br><tt> "Decay of Delay 4.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 12: sxedelay[4]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 4",</tt>
<br><tt> "Delay 4 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt> 175
/* 175/255 of MAX_DELAY msec */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 13: sxedecay[5]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 5",</tt>
<br><tt> "Decay of Delay 5.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 14: sxedelay[5]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 5",</tt>
<br><tt> "Delay 5 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt> 210
/* 210/255 of MAX_DELAY msec */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 15: sxedecay[6]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Decay 6",</tt>
<br><tt> "Decay of Delay 6.",</tt>
<br><tt> 0,</tt>
<br><tt> 255,
/* 1 */</tt>
<br><tt> 220
/* 220/255 of dynamic range */</tt>
<br><tt> },</tt><tt></tt>
<p><tt> /* Param 16: sxedelay[6]. */</tt>
<br><tt> {</tt>
<br><tt> small,</tt>
<br><tt> slider,</tt>
<br><tt> "Delay 6",</tt>
<br><tt> "Delay 6 time: 0 to maximum.",</tt>
<br><tt>
0,
/* 0 msec. delay */</tt>
<br><tt> 255,
/* Highest delay of range */</tt>
<br><tt> 255</tt>
<br><tt> },</tt>
<br><tt>};</tt><tt></tt>
<p><tt>/* State of the soxflanger. */</tt><tt></tt>
<p><tt>typedef struct {</tt>
<br><tt> float in_gain, out_gain;</tt>
<br><tt> float sxedecay[MAX_ECHOS];</tt>
<br><tt> float sxedelay[MAX_ECHOS];</tt>
<br><tt> long sxesamples[MAX_ECHOS];</tt>
<br><tt> samp* delaybuf;</tt>
<br><tt> long maxsxesamples;</tt>
<br><tt> float maxanydelay;</tt>
<br><tt> long maxanysmpls;</tt>
<br><tt> int echo_cntr;</tt>
<br><tt> int num_delays;</tt>
<br><tt> double unit_vol_scale;</tt>
<br><tt> double exp_fract_vol;</tt>
<br><tt>// double unit_delay_scale;</tt>
<br><tt>// double exp_fract_delay;</tt>
<br><tt>// float clip; // Put in after clipping
issue decided.</tt>
<br><tt>} sxe_state;</tt><tt></tt>
<p><tt>int ox_init(machine_type *t) {</tt>
<br><tt> t->long_name = "soxecho (Mueller, Nelson, et al.)";</tt>
<br><tt> t->short_name = "soxecho";</tt>
<br><tt> t->max_tracks = 1;</tt>
<br><tt> t->input_channels = 1;</tt>
<br><tt> t->output_channels = 1;</tt>
<br><tt> t->num_params = 15;</tt>
<br><tt> t->param_specs = sxe_params;</tt><tt></tt>
<p><tt> return 1;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_create(machine *m) {</tt>
<br><tt> sxe_state *s;</tt>
<br><tt> int i;</tt><tt></tt>
<p><tt> s = (sxe_state*) malloc(sizeof(sxe_state));</tt>
<br><tt> m->state = (void *) s;</tt><tt></tt>
<p><tt> s->num_delays = MAX_ECHOS;</tt><tt></tt>
<p><tt> s->unit_vol_scale = exp((double)DYN_RANGE / 3.0) -
1.0;</tt>
<br><tt> s->exp_fract_vol = ((double)DYN_RANGE / 3) /
255.0;</tt><tt></tt>
<p><tt>// s->unit_delay_scale = (exp((double)DELAY_LOG_RANGE)
- 1.0) / (double)MAX_DELAY;</tt>
<br><tt>// s->exp_fract_delay = ((double)DELAY_LOG_RANGE) /
255.0;</tt><tt></tt>
<p><tt> /* Set defaults */</tt>
<br><tt>// s->in_gain = 128.0/255.0;
/* Approx. 0.5 (old) */</tt>
<br><tt> s->in_gain = ((float)exp(70.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt>// s->out_gain =
1.0;
/* No attenuation */</tt>
<br><tt> s->out_gain = ((float)exp(100.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt><tt></tt>
<p><tt> s->sxedecay[0] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[0] = 35.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[0] = ((float)exp(35.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxedecay[1] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[1] = 70.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[1] = ((float)exp(70.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxedecay[2] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[2] = 105.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[2] = ((float)exp(105.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxedecay[3] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[3] = 140.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[3] = ((float)exp(140.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxedecay[4] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[4] = 175.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[4] = ((float)exp(175.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxedecay[5] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[5] = 210.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[5] = ((float)exp(210.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxedecay[6] = ((float)exp(220.0 * s->exp_fract_vol)
- 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> s->sxedelay[6] = 255.0/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[6] = ((float)exp(255.0 * s->exp_fract_delay)
- 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->maxsxesamples = 0; // Accumulated in following
loop.</tt><tt></tt>
<p><tt> for(i = 0; i < s->num_delays; i++) {</tt><tt></tt>
<p><tt> s->sxesamples[i] = (long)(s->sxedelay[i]
* ((float)OX_SAMPLING_RATE) / 1000.0</tt>
<br><tt>
+ 0.5);</tt><tt></tt>
<p><tt> if (s->sxesamples[i] >
s->maxsxesamples)</tt>
<br><tt> s->maxsxesamples
= s->sxesamples[i];</tt>
<br><tt> }</tt><tt></tt>
<p><tt> s->maxanydelay = MAX_DELAY;</tt>
<br><tt> s->maxanysmpls = (long)(s->maxanydelay *
((float)OX_SAMPLING_RATE)
/ 1000.0</tt>
<br><tt>
+ 10.0); // Add an extra 10 to avoid possible exp overflow.</tt><tt></tt>
<p><tt> if (s->maxsxesamples > s->maxanysmpls)</tt>
<br><tt> printf("%s\n",</tt>
<br><tt>
"soxecho create error: delay sample requirement exceeds
allocation.");</tt><tt></tt>
<p><tt> s->echo_cntr = 0;</tt><tt></tt>
<p><tt>// s->clip = ?; // Put in after clipping issue
decided.</tt><tt></tt>
<p><tt> s->delaybuf = (samp*) malloc(s->maxanysmpls *
sizeof(samp));</tt><tt></tt>
<p><tt>// memset(s->delaybuf, 0, s->maxanysmpls * sizeof(samp));
// Unnecessary</tt>
<br><tt>// for (i=0, i < s->maxanysmpls, i++) s->delaybuf[i]
= 0; // Speed this up.</tt><tt></tt>
<p><tt> return;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_destroy(machine *m) {</tt>
<br><tt> free( ((sxe_state*)(m->state)) -> delaybuf);</tt>
<br><tt> free(m->state);</tt>
<br><tt> m->state = NULL;</tt>
<br><tt> return;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_update(machine *m) {</tt>
<br><tt> sxe_state *s = (sxe_state *) m->state;</tt>
<br><tt> param temp;</tt>
<br><tt> int i;</tt><tt></tt>
<p><tt> temp = m->params[0][ix_gain_in];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt> if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp > 255) temp = 255;</tt>
<br><tt> s->in_gain = ((float)temp) / 255.0;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_gain_out];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt> if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp > 255) temp = 255;</tt>
<br><tt> s->out_gain = ((float)temp) / 255.0;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> s->maxsxesamples = 0; // Accumulated in following
loop.</tt><tt></tt>
<p><tt> for(i = 0; i < s->num_delays; i++) {</tt>
<br><tt> temp = m->params[0][ix_sxedecay
+ i * 2];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt>
if
(temp < 0) temp = 0;</tt>
<br><tt> else if (temp
> 255) temp = 255;</tt><tt></tt>
<p><tt>// s->sxedecay[i]
= ((float)temp) / 255.0; // Old, linear.</tt>
<br><tt> s->sxedecay[i]
= ((float)exp((double)temp * s->exp_fract_vol) - 1.0)</tt>
<br><tt>
/ s->unit_vol_scale;</tt>
<br><tt> }</tt><tt></tt>
<p><tt> temp = m->params[0][ix_sxedelay +
i * 2];</tt>
<br><tt> if (temp != nochange) {</tt>
<br><tt>
if
(temp < 1) temp = 0;</tt>
<br><tt> else if (temp
> 255) temp = 255;</tt><tt></tt>
<p><tt> s->sxedelay[i]
= (float)temp/255.0 * (float)MAX_DELAY;</tt>
<br><tt>// s->sxedelay[i]
= ((float)exp((double)temp * s->exp_fract_delay) - 1.0)</tt>
<br><tt>//
/ s->unit_delay_scale;</tt><tt></tt>
<p><tt> s->sxesamples[i]
= (long)(s->sxedelay[i] * ((float)OX_SAMPLING_RATE)</tt>
<br><tt>
/ 1000.0 + 0.5);</tt>
<br><tt> }</tt><tt></tt>
<p><tt> if (s->sxesamples[i] > s->maxsxesamples)
s->maxsxesamples = s->sxesamples[i];</tt>
<br><tt> }</tt><tt></tt>
<p><tt> return;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>const char *ox_desc(int which_param, param value) {</tt>
<br><tt> static char temp_string[80];</tt>
<br><tt> float user_value;</tt>
<br><tt> int i, percent;</tt><tt></tt>
<p><tt> sprintf(temp_string, "ERROR");</tt><tt></tt>
<p><tt> if (which_param == ix_gain_in || which_param == ix_gain_out)
{</tt>
<br><tt> if (value < 0) sprintf(temp_string,
"Gain must not be less than 0.");</tt>
<br><tt> else {</tt>
<br><tt> user_value =
((float)value)
/ 255.0;</tt>
<br><tt> percent =
(int)(user_value
* 100.0 + 0.5);</tt>
<br><tt> if (user_value
> 1.0)</tt>
<br><tt>
sprintf(temp_string, "Gain must be not greater than 100%.");</tt><tt></tt>
<p><tt> else
sprintf(temp_string,
"%3.0f%%", percent);</tt>
<br><tt> }</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else {</tt><tt></tt>
<p><tt> for(i = 0; i < MAX_ECHOS; i++)
{</tt>
<br><tt> if (which_param
== (ix_sxedecay + i * 2)) {</tt>
<br><tt>
if (value < 0) sprintf(temp_string, "%s%d%s", "Decay ", i,</tt>
<br><tt>
" must not be less than 0. 0 removes this delay");</tt>
<br><tt>
else {</tt>
<br><tt>
user_value = (float)value / 255.0;</tt>
<br><tt>
percent = (int)(user_value * 100.0 + 0.5);</tt>
<br><tt>
if (user_value > 1.0)</tt>
<br><tt>
sprintf(temp_string, "%s%d%s", "Decay ", i,</tt>
<br><tt>
"must be not greater than 100%.");</tt><tt></tt>
<p><tt>
else sprintf(temp_string, "%3.0f%%", percent);</tt>
<br><tt>
}</tt>
<br><tt> }</tt><tt></tt>
<p><tt> else if (which_param
== (ix_sxedelay + i * 2)) {</tt>
<br><tt>
if (value < 0) sprintf(temp_string, "%s%d%s", "Delay ", i,</tt>
<br><tt>
" must not be less than 0.");</tt>
<br><tt>
else {</tt>
<br><tt>
user_value = (float)value / 255.0 * (float)MAX_DELAY;</tt>
<br><tt>
if (user_value > (float)MAX_DELAY)</tt>
<br><tt>
sprintf(temp_string, "%s%d%s%d%s", "Delay ", i,</tt>
<br><tt>
" must be not greater than ", MAX_DELAY, " msec.");</tt><tt></tt>
<p><tt>
else sprintf(temp_string, "%3.0f", user_value);</tt>
<br><tt>
}</tt>
<br><tt> }</tt>
<br><tt> }</tt>
<br><tt> }</tt>
<br><tt>/*</tt>
<br><tt> *** Be nice and check the hint with warning, if...</tt>
<br><tt> sum_in_volume = 1.0;</tt>
<br><tt> for ( i = 0; i < echo->num_delays; i++ )</tt>
<br><tt> sum_in_volume += echo->decay[i];</tt>
<br><tt> if ( sum_in_volume * echo->in_gain > 1.0 / echo->out_gain
)</tt>
<br><tt> warn("echo: warning >>> gain-out can cause saturation of
output <<<");</tt>
<br><tt>*/</tt><tt></tt>
<p><tt>}</tt><tt></tt>
<p><tt>int ox_work(machine *m, int block_size) {</tt>
<br><tt> sxe_state *s = (sxe_state *) m->state;</tt><tt></tt>
<p><tt> samp* mlout;</tt>
<br><tt>// samp* mrout;</tt>
<br><tt> samp* mlin;</tt>
<br><tt> samp* dlybf;</tt>
<br><tt> int i, j, dlyctr, numdlys;</tt>
<br><tt> long mxsmpls;</tt>
<br><tt> long* sxsmpls;</tt>
<br><tt> float ingain, outgain;</tt>
<br><tt> float* sxdecy;</tt><tt></tt>
<p><tt>// float hig_clip, low_clip;</tt><tt></tt>
<p><tt>// hig_clip = s->clip;</tt>
<br><tt>// low_clip = 0.0 - s->clip;</tt><tt></tt>
<p><tt> /* Speed-up variables */</tt>
<br><tt> mlout = m->lout;</tt>
<br><tt>// mrout = m->rout;</tt>
<br><tt> mlin = m->lin;</tt><tt></tt>
<p><tt> ingain = s->in_gain;</tt>
<br><tt> outgain = s->out_gain;</tt><tt></tt>
<p><tt> dlybf = s->delaybuf;</tt>
<br><tt> dlyctr = s->echo_cntr;</tt>
<br><tt> mxsmpls = s->maxsxesamples;</tt>
<br><tt> sxsmpls = (long*)s->sxesamples;</tt>
<br><tt> sxdecy = (float*)s->sxedecay;</tt>
<br><tt> numdlys = s->num_delays;</tt><tt></tt>
<p><tt> for(i = 0; i < block_size; i++) {</tt><tt></tt>
<p><tt> /* Compute output first */</tt>
<br><tt> mlout[i] = mlin[i] *
ingain;</tt><tt></tt>
<p><tt> for (j = 0; j < numdlys; j++)
{</tt>
<br><tt> if (sxdecy[j])</tt>
<br><tt>
mlout[i] += dlybf[(mxsmpls + dlyctr - sxsmpls[j]) % mxsmpls] * sxdecy[j];</tt>
<br><tt> }</tt><tt></tt>
<p><tt> /* Adjust the output volume */</tt>
<br><tt> mlout[i] *= outgain;</tt><tt></tt>
<p><tt> /* Check for clipping. Not
implemented. */</tt>
<br><tt>// if
(mlout[i] > hig_clip) mlout[i] = hig_clip;</tt>
<br><tt>// else if (mlout[i] < low_clip)
mlout[i] = low_clip;</tt><tt></tt>
<p><tt> /* Store input in delay buffer */</tt>
<br><tt> dlybf[dlyctr] = mlin[i];</tt><tt></tt>
<p><tt> /* Adjust the counter */</tt>
<br><tt> dlyctr = ( dlyctr + 1 ) %
mxsmpls;</tt><tt></tt>
<p><tt> }</tt><tt></tt>
<p><tt> /* Save speed-up variables to state. */</tt>
<br><tt> s->echo_cntr = dlyctr;</tt><tt></tt>
<p><tt> return 1;</tt>
<br><tt>}</tt><tt></tt>
<p><tt>void ox_track(machine *m, int change) {</tt>
<br><tt>}</tt>
<br><tt></tt> </html>
--------------801DEC9FA1F1642AE599CEF4--
| [Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Octal-dev] SoX flanger,
Avelino <=