aboutsummaryrefslogtreecommitdiffstats
path: root/src/waveguide_mesh_2670.c
diff options
context:
space:
mode:
authorDavid Robillard <d@drobilla.net>2012-08-19 19:24:44 +0000
committerDavid Robillard <d@drobilla.net>2012-08-19 19:24:44 +0000
commit6a651aca0de7262e805ff8edb12cdc4692009997 (patch)
treef14451d008ea578c162df85d2e2605a227bf4bf6 /src/waveguide_mesh_2670.c
parent5de48cab1d561e6b37502fe823922b94d1489eda (diff)
Move omins to plugins directory.
git-svn-id: http://svn.drobilla.net/lad/trunk/plugins/omins.lv2@4725 a436a847-0d15-0410-975c-d299462d15a1
Diffstat (limited to 'src/waveguide_mesh_2670.c')
-rw-r--r--src/waveguide_mesh_2670.c355
1 files changed, 0 insertions, 355 deletions
diff --git a/src/waveguide_mesh_2670.c b/src/waveguide_mesh_2670.c
deleted file mode 100644
index d767dce..0000000
--- a/src/waveguide_mesh_2670.c
+++ /dev/null
@@ -1,355 +0,0 @@
-/* This file is an audio plugin.
- * Copyright 2005 Loki Davison.
- *
- * Based on code by Brook Eaton and a couple of papers...
- *
- * Implements a Waveguide Mesh drum. FIXME to be extended, to have rimguides, power normalisation and all
- * manner of other goodies.
- *
- * Tension is well, drum tension
- * power is how hard you hit it.
- *
- *
- * This plugin is free software; you can redistribute it and/or modify it under
- * the terms of the GNU General Public License as published by the Free
- * Software Foundation; either version 2 of the License, or (at your option)
- * any later version.
- *
- * This plugin is distributed in the hope that it will be useful, but WITHOUT
- * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
- * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for
- * details.
- *
- * You should have received a copy of the GNU General Public License along
- * with this program; if not, write to the Free Software Foundation, Inc.,
- * 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
- */
-
-#define _XOPEN_SOURCE 500 /* strdup */
-#include <stdlib.h>
-#include <string.h>
-#include <math.h>
-
-#include "ladspa.h"
-
-#define MESH_BASE_ID 2670
-
-#define MESH_NUM_PORTS 6
-
-/* Port Numbers */
-#define MESH_INPUT1 0
-#define MESH_OUTPUT 1
-#define MESH_TENSION 2
-#define MESH_POWER 3
-#define MESH_EX_X 4
-#define MESH_EX_Y 5
-
-#define LENGTH 8 // must be divisible by 4!!
-#define WIDTH 8
-#define SIZE LENGTH*WIDTH // Size of mesh
-#define INITIAL 0 // initial values stored at junctions
-#define LOSS 0.2 // loss in wave propagation
-#define INIT_DELTA 6 // initial values
-#define INIT_T 0.1 // for the impedances
-#define INIT_GAMMA 8 //
-#define PORTS 8 //for the initialization of junction velocities from excitation
-
-// 2D array of junctions.
-// The important function of the junction is to store
-// the velocities of travelling waves so that adjacent junction's
-// velocities can be calculated.
-typedef struct _junction
-{
- LADSPA_Data v_junction; // junction velocity
- LADSPA_Data n_junction; // velocity heading north into junction
- LADSPA_Data s_junction; // velocity heading south into junction
- LADSPA_Data e_junction; // velocity heading east into junction
- LADSPA_Data w_junction; // velocity heading west into junction
- LADSPA_Data c_junction; // velocity heading back into the junction (delayed)
- LADSPA_Data s_temp; // these two 'temp' values are used because calculation
- LADSPA_Data e_temp; // of the mesh updates s/e values before they are used
-} _junction; // to calculate north and south velocities!
-
-/* All state information for plugin */
-
-typedef struct
-{
- /* Ports */
- LADSPA_Data* input1;
- LADSPA_Data* output;
- LADSPA_Data* tension;
- LADSPA_Data* power;
- LADSPA_Data* ex_x;
- LADSPA_Data* ex_y;
-
- /* vars */
- _junction mesh[LENGTH][WIDTH];
- LADSPA_Data last_trigger;
-
-} WgMesh;
-
-/* Construct a new plugin instance */
-LADSPA_Handle
-wgmesh_instantiate(const LADSPA_Descriptor* descriptor,
- unsigned long srate)
-{
- WgMesh * plugin = (WgMesh *) malloc (sizeof (WgMesh));
- int i, j;
-
- //// Init Mesh
- for (i=0; i<LENGTH; i++) {
- for (j=0; j<WIDTH; j++) {
- plugin->mesh[i][j].v_junction = INITIAL;
- plugin->mesh[i][j].n_junction = INITIAL;
- plugin->mesh[i][j].s_junction = INITIAL;
- plugin->mesh[i][j].e_junction = INITIAL;
- plugin->mesh[i][j].w_junction = INITIAL;
- plugin->mesh[i][j].c_junction = INITIAL;
- plugin->mesh[i][j].s_temp = INITIAL;
- plugin->mesh[i][j].e_temp = INITIAL;
- }
- }
- plugin->last_trigger = 0.0;
- return (LADSPA_Handle)plugin;
-}
-
-/* Connect a port to a data location */
-void
-wgmesh_connect_port(LADSPA_Handle instance,
- unsigned long port,
- LADSPA_Data* location)
-{
- WgMesh* plugin;
-
- plugin = (WgMesh*)instance;
- switch (port) {
- case MESH_INPUT1:
- plugin->input1 = location;
- break;
- case MESH_OUTPUT:
- plugin->output = location;
- break;
- case MESH_TENSION:
- plugin->tension = location;
- break;
- case MESH_POWER:
- plugin->power = location;
- break;
- case MESH_EX_X:
- plugin->ex_x = location;
- break;
- case MESH_EX_Y:
- plugin->ex_y = location;
- break;
- }
-}
-
-inline static void excite_mesh(WgMesh* plugin, LADSPA_Data power, LADSPA_Data ex_x, LADSPA_Data ex_y)
-{
- int i=ex_x,j=ex_y;
- LADSPA_Data temp;
- LADSPA_Data Yj;
-
- Yj = 2*(INIT_DELTA*INIT_DELTA/((INIT_T*INIT_T)*(INIT_GAMMA*INIT_GAMMA))); // junction admittance
- temp = power*2/(LENGTH+WIDTH);
- plugin->mesh[i][j].v_junction = plugin->mesh[i][j].v_junction + temp;
- plugin->mesh[i][j].n_junction = plugin->mesh[i][j].n_junction + Yj*temp/PORTS;
- // All velocities leaving the junction are equal to
- plugin->mesh[i][j].s_junction = plugin->mesh[i][j].s_junction + Yj*temp/PORTS;
- // the total velocity in the junction * the admittance
- plugin->mesh[i][j].e_junction = plugin->mesh[i][j].e_junction + Yj*temp/PORTS;
- // divided by the number of outgoing ports.
- plugin->mesh[i][j].w_junction = plugin->mesh[i][j].w_junction + Yj*temp/PORTS;
- //mesh[i][j].c_junction = 0;
-}
-
-void
-wgmesh_run_cr(LADSPA_Handle instance, unsigned long nframes)
-{
- WgMesh* plugin = (WgMesh*)instance;
- LADSPA_Data tension = *(plugin->tension);
- LADSPA_Data ex_x = *(plugin->ex_x);
- LADSPA_Data ex_y = *(plugin->ex_y);
- LADSPA_Data* input = plugin->input1;
- LADSPA_Data* out = plugin->output;
- LADSPA_Data* power = plugin->power;
- LADSPA_Data last_trigger = plugin->last_trigger;
- size_t i,j,k;
- LADSPA_Data filt, trg, oldfilt;
- LADSPA_Data Yc,Yj,tempN,tempS,tempE,tempW;
-
- // Set input variables //
- oldfilt = plugin->mesh[LENGTH-LENGTH/4][WIDTH-WIDTH/4].v_junction;
-
- for (k=0; k<nframes; k++) {
- if (tension==0)
- tension=0.0001;
- trg = input[k];
-
- if (trg > 0.0f && !(last_trigger > 0.0f))
- {
- //printf("got trigger, exciting mesh, %f \n", tension);
- excite_mesh(plugin, power[k], ex_x, ex_y);
- }
-
- //junction admitance
- Yj = 2*INIT_DELTA*INIT_DELTA/(( (tension)*((tension) )*(INIT_GAMMA*INIT_GAMMA)));
- Yc = Yj-4; // junction admittance (left shift is for multiply by 2!)
- //plugin->v_power = power[k];
-
- for (i=1; i<LENGTH-1; i++) {
- // INNER MESH //
- for (j=1; j<WIDTH-1; j++) { // to multiply by 2 - simply shift to the left by 1!
- plugin->mesh[i][j].v_junction = 2.0*(plugin->mesh[i][j].n_junction + plugin->mesh[i][j].s_junction
- + plugin->mesh[i][j].e_junction + plugin->mesh[i][j].w_junction + Yc*plugin->mesh[i][j].c_junction)/Yj;
-
- plugin->mesh[i][j+1].s_junction = plugin->mesh[i][j].v_junction - plugin->mesh[i][j].n_junction;
- plugin->mesh[i][j-1].n_junction = plugin->mesh[i][j].v_junction - plugin->mesh[i][j].s_temp;
-
- plugin->mesh[i+1][j].e_junction = plugin->mesh[i][j].v_junction - plugin->mesh[i][j].w_junction;
- plugin->mesh[i-1][j].w_junction = plugin->mesh[i][j].v_junction - plugin->mesh[i][j].e_temp;
-
- plugin->mesh[i][j].c_junction = plugin->mesh[i][j].v_junction - plugin->mesh[i][j].c_junction;
-
- plugin->mesh[i][j].s_temp = plugin->mesh[i][j].s_junction; //
- plugin->mesh[i][j].e_temp = plugin->mesh[i][j].e_junction; // update current values in the temp slots!
- }
- // BOUNDARY //
- tempS = plugin->mesh[i][0].s_junction;
- plugin->mesh[i][0].s_junction = -plugin->mesh[i][0].n_junction;
- plugin->mesh[i][1].s_junction = plugin->mesh[i][1].s_temp = tempS;
- tempN = plugin->mesh[i][WIDTH-1].n_junction;
- plugin->mesh[i][WIDTH-1].n_junction = -plugin->mesh[i][WIDTH-1].s_junction;
- plugin->mesh[i][WIDTH-2].n_junction = tempN;
- // 'i's in the neplugint few lines are really 'j's!
- tempE = plugin->mesh[0][i].e_junction;
- plugin->mesh[0][i].e_junction = -plugin->mesh[0][i].w_junction;
- plugin->mesh[1][i].e_junction = plugin->mesh[1][i].e_temp = tempE;
- tempW = plugin->mesh[WIDTH-1][i].w_junction;
- plugin->mesh[WIDTH-1][i].w_junction = -plugin->mesh[WIDTH-1][i].e_junction;
- plugin->mesh[WIDTH-2][i].w_junction = tempW;
- }
-
- filt = LOSS*(plugin->mesh[LENGTH-LENGTH/4][WIDTH-WIDTH/4].v_junction + oldfilt);
- oldfilt = plugin->mesh[LENGTH-LENGTH/4][WIDTH-WIDTH/4].v_junction;
- plugin->mesh[LENGTH-LENGTH/4][WIDTH-WIDTH/4].v_junction = filt;
-
- out[k] = plugin->mesh[WIDTH/4][WIDTH/4-1].v_junction;
- last_trigger = trg;
- }
- plugin->last_trigger = last_trigger;
-
-}
-
-void
-wgmesh_cleanup(LADSPA_Handle instance)
-{
- free(instance);
-}
-
-LADSPA_Descriptor* wg_mesh_cr_desc = NULL;
-
-/* Called automatically when the plugin library is first loaded. */
-void
-_init()
-{
- char** port_names;
- LADSPA_PortDescriptor* port_descriptors;
- LADSPA_PortRangeHint* port_range_hints;
-
- wg_mesh_cr_desc = (LADSPA_Descriptor*)malloc(sizeof(LADSPA_Descriptor));
-
- if (wg_mesh_cr_desc) {
-
- wg_mesh_cr_desc->UniqueID = MESH_BASE_ID;
- wg_mesh_cr_desc->Label = strdup("wg_mesh_cr");
- wg_mesh_cr_desc->Properties = LADSPA_PROPERTY_HARD_RT_CAPABLE;
- wg_mesh_cr_desc->Name = strdup("Simple waveguide mesh (CR Controls)");
- wg_mesh_cr_desc->Maker = strdup("Loki Davison");
- wg_mesh_cr_desc->Copyright = strdup("GPL");
- wg_mesh_cr_desc->PortCount = MESH_NUM_PORTS;
- port_descriptors = (LADSPA_PortDescriptor*)calloc(MESH_NUM_PORTS, sizeof(LADSPA_PortDescriptor));
- wg_mesh_cr_desc->PortDescriptors = (const LADSPA_PortDescriptor*)port_descriptors;
- port_descriptors[MESH_TENSION] = LADSPA_PORT_INPUT | LADSPA_PORT_CONTROL;
- port_descriptors[MESH_POWER] = LADSPA_PORT_INPUT | LADSPA_PORT_AUDIO;
- port_descriptors[MESH_INPUT1] = LADSPA_PORT_INPUT | LADSPA_PORT_AUDIO;
- port_descriptors[MESH_OUTPUT] = LADSPA_PORT_OUTPUT | LADSPA_PORT_AUDIO;
- port_descriptors[MESH_EX_X] = LADSPA_PORT_INPUT | LADSPA_PORT_CONTROL;
- port_descriptors[MESH_EX_Y] = LADSPA_PORT_INPUT | LADSPA_PORT_CONTROL;
- port_names = (char**)calloc(MESH_NUM_PORTS, sizeof(char*));
- wg_mesh_cr_desc->PortNames = (const char**)port_names;
- port_names[MESH_TENSION] = strdup("Tension");
- port_names[MESH_POWER] = strdup("Power");
- port_names[MESH_INPUT1] = strdup("Trigger");
- port_names[MESH_OUTPUT] = strdup("Output");
- port_names[MESH_EX_X] = strdup("Excitation X");
- port_names[MESH_EX_Y] = strdup("Excitation Y");
- port_range_hints = ((LADSPA_PortRangeHint *)
- calloc(MESH_NUM_PORTS, sizeof(LADSPA_PortRangeHint)));
- wg_mesh_cr_desc->PortRangeHints = (const LADSPA_PortRangeHint*)port_range_hints;
- port_range_hints[MESH_TENSION].HintDescriptor = LADSPA_HINT_BOUNDED_BELOW | LADSPA_HINT_BOUNDED_ABOVE
-| LADSPA_HINT_DEFAULT_MIDDLE;
- port_range_hints[MESH_TENSION].LowerBound = 0.0001f;
- port_range_hints[MESH_TENSION].UpperBound = 0.22f;
- port_range_hints[MESH_POWER].HintDescriptor = LADSPA_HINT_DEFAULT_1 | LADSPA_HINT_BOUNDED_BELOW;
- port_range_hints[MESH_POWER].LowerBound = 0.000f;
- port_range_hints[MESH_POWER].UpperBound = 20.000f;
-
- port_range_hints[MESH_EX_X].HintDescriptor = LADSPA_HINT_DEFAULT_1 | LADSPA_HINT_BOUNDED_BELOW | LADSPA_HINT_BOUNDED_ABOVE | LADSPA_HINT_INTEGER;
- port_range_hints[MESH_EX_X].LowerBound = 0.950f; //1
- port_range_hints[MESH_EX_X].UpperBound = LENGTH-0.99;//length-1 ish
-
- port_range_hints[MESH_EX_Y].HintDescriptor = LADSPA_HINT_DEFAULT_1 | LADSPA_HINT_BOUNDED_BELOW | LADSPA_HINT_BOUNDED_ABOVE | LADSPA_HINT_INTEGER;
- port_range_hints[MESH_EX_Y].LowerBound = 0.950f; //1
- port_range_hints[MESH_EX_Y].UpperBound = LENGTH-0.99;//length-1 ish
-
- port_range_hints[MESH_INPUT1].HintDescriptor = 0;
- port_range_hints[MESH_OUTPUT].HintDescriptor = 0;
- wg_mesh_cr_desc->instantiate = wgmesh_instantiate;
- wg_mesh_cr_desc->connect_port = wgmesh_connect_port;
- wg_mesh_cr_desc->activate = NULL;
- wg_mesh_cr_desc->run = wgmesh_run_cr;
- wg_mesh_cr_desc->run_adding = NULL;
- wg_mesh_cr_desc->set_run_adding_gain = NULL;
- wg_mesh_cr_desc->deactivate = NULL;
- wg_mesh_cr_desc->cleanup = wgmesh_cleanup;
- }
-}
-
-void
-wgmesh_delete_descriptor(LADSPA_Descriptor* psDescriptor)
-{
- unsigned long lIndex;
- if (psDescriptor) {
- free((char*)psDescriptor->Label);
- free((char*)psDescriptor->Name);
- free((char*)psDescriptor->Maker);
- free((char*)psDescriptor->Copyright);
- free((LADSPA_PortDescriptor *)psDescriptor->PortDescriptors);
- for (lIndex = 0; lIndex < psDescriptor->PortCount; lIndex++)
- free((char*)(psDescriptor->PortNames[lIndex]));
- free((char**)psDescriptor->PortNames);
- free((LADSPA_PortRangeHint *)psDescriptor->PortRangeHints);
- free(psDescriptor);
- }
-}
-
-/* Called automatically when the library is unloaded. */
-void
-_fini()
-{
- wgmesh_delete_descriptor(wg_mesh_cr_desc);
-}
-
-/* Return a descriptor of the requested plugin type. */
-const LADSPA_Descriptor*
-ladspa_descriptor(unsigned long Index)
-{
- switch (Index) {
- case 0:
- return wg_mesh_cr_desc;
- default:
- return NULL;
- }
-}
-