/*
  An LV2 plugin to generate a bandlimited variable pulse waveform.
  Copyright 2011 David Robillard
  Copyright 2002 Mike Rawes

  This 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 3 of the License, or
  (at your option) any later version.

  This software 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 more details.

  You should have received a copy of the GNU General Public License
  along with this software.  If not, see <http://www.gnu.org/licenses/>.
*/

#include <stdlib.h>
#include "lv2/lv2plug.in/ns/ext/morph/morph.h"
#include "lv2/lv2plug.in/ns/lv2core/lv2.h"
#include "uris.h"
#include "wavedata.h"

#define PULSE_FREQUENCY  0
#define PULSE_PULSEWIDTH 1
#define PULSE_OUTPUT     2

typedef struct {
	float*   frequency;
	float*   pulsewidth;
	float*   output;
	float    phase;
	uint32_t frequency_is_cv;
	uint32_t pulsewidth_is_cv;
	Wavedata wdat;
	URIs     uris;
} Pulse;

static void
connect_port(LV2_Handle instance,
             uint32_t   port,
             void*      data)
{
	Pulse* plugin = (Pulse*)instance;

	switch (port) {
	case PULSE_FREQUENCY:
		plugin->frequency = (float*)data;
		break;
	case PULSE_PULSEWIDTH:
		plugin->pulsewidth = (float*)data;
		break;
	case PULSE_OUTPUT:
		plugin->output = (float*)data;
		break;
	}
}

static LV2_Morph_Status
morph_port(LV2_Handle                      instance,
           uint32_t                        port,
           LV2_URID                        type,
           const LV2_Morph_Property*const* properties)
{
	Pulse* plugin = (Pulse*)instance;

	if (type != plugin->uris.lv2_ControlPort &&
	    type != plugin->uris.lv2_CVPort) {
		return LV2_MORPH_ERR_BAD_TYPE;
	}

	switch (port) {
	case PULSE_FREQUENCY:
		plugin->frequency_is_cv = (type == plugin->uris.lv2_CVPort);
		break;
	case PULSE_PULSEWIDTH:
		plugin->pulsewidth_is_cv = (type == plugin->uris.lv2_CVPort);
		break;
	default:
		return LV2_MORPH_ERR_BAD_PORT;
	}

	return LV2_MORPH_SUCCESS;
}

static LV2_Handle
instantiate(const LV2_Descriptor*     descriptor,
            double                    sample_rate,
            const char*               bundle_path,
            const LV2_Feature* const* features)
{
	Pulse* plugin = (Pulse*)malloc(sizeof(Pulse));
	if (!plugin) {
		return NULL;
	}

	if (wavedata_load(&plugin->wdat, bundle_path, "sawtooth_data",
	                  BLOP_DLSYM_SAWTOOTH, sample_rate)) {
		free(plugin);
		return 0;
	}

	plugin->frequency_is_cv = 0;
	plugin->pulsewidth_is_cv = 0;
	map_uris(&plugin->uris, features);
	wavedata_get_table(&plugin->wdat, 440.0);

	return (LV2_Handle)plugin;
}

static void
cleanup(LV2_Handle instance)
{
	Pulse* plugin = (Pulse*)instance;

	wavedata_unload(&plugin->wdat);
	free(instance);
}

static void
activate(LV2_Handle instance)
{
	Pulse* plugin = (Pulse*)instance;

	plugin->phase = 0.0f;
}

static void
run(LV2_Handle instance,
    uint32_t   sample_count)
{
	Pulse* plugin = (Pulse*)instance;

	/* Frequency (array of float of length sample_count) */
	const float* frequency = plugin->frequency;

	/* Pulse Width (array of float of length sample_count) */
	const float* pulsewidth = plugin->pulsewidth;

	/* Output (pointer to float value) */
	float* output = plugin->output;

	/* Instance data */
	Wavedata* wdat  = &plugin->wdat;
	float     phase = plugin->phase;

	float last_pwidth = pulsewidth[0];
	float pwidth      = f_clip(last_pwidth, 0.0f, 1.0f);
	float dc_shift    = 1.0 - (2.0 * pwidth);
	float phase_shift = pwidth * wdat->sample_rate;

	for (uint32_t s = 0; s < sample_count; ++s) {
		const float freq = frequency[s * plugin->frequency_is_cv];
		if (freq != wdat->frequency) {
			/* Frequency changed, look up table to play */
			wavedata_get_table(wdat, freq);
		}

		const float this_pwidth = pulsewidth[s * plugin->pulsewidth_is_cv];
		if (this_pwidth != last_pwidth) {
			/* Pulsewidth changed, recalculate */
			last_pwidth = this_pwidth;
			pwidth      = f_clip(this_pwidth, 0.0f, 1.0f);
			dc_shift    = 1.0f - (2.0f * pwidth);
			phase_shift = pwidth * wdat->sample_rate;
		}

		/* Get samples from sawtooth and phase shifted inverted sawtooth,
		   with approriate DC offset */
		output[s] = wavedata_get_sample(wdat, phase)
		    - wavedata_get_sample(wdat, phase + phase_shift)
			+ dc_shift;

		/* Update phase, wrapping if necessary */
		phase += wdat->frequency;
		if (phase < 0.0f) {
			phase += wdat->sample_rate;
		} else if (phase > wdat->sample_rate) {
			phase -= wdat->sample_rate;
		}
	}
	plugin->phase = phase;
}

static const void*
extension_data(const char* uri)
{
	static const LV2_Morph_Interface morph = { morph_port, NULL };
	if (!strcmp(uri, LV2_MORPH__interface)) {
		return &morph;
	}
	return NULL;
}

static const LV2_Descriptor descriptor = {
	"http://drobilla.net/plugins/blop/pulse",
	instantiate,
	connect_port,
	activate,
	run,
	NULL,
	cleanup,
	extension_data,
};

LV2_SYMBOL_EXPORT const LV2_Descriptor*
lv2_descriptor(uint32_t index)
{
	switch (index) {
	case 0:  return &descriptor;
	default: return NULL;
	}
}