2 * Copyright (C) 2016 Damien Zammit <damien@zamaudio.com>
4 * This program is free software; you can redistribute it and/or
5 * modify it under the terms of the GNU General Public License
6 * as published by the Free Software Foundation; either version 2
7 * of the License, or (at your option) any later version.
9 * This program is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
21 #include <cairo/cairo.h>
22 #include "ardour/lv2_extensions.h"
25 #include "lv2/lv2plug.in/ns/lv2core/lv2.h"
27 #define ACOMP_URI "urn:ardour:a-comp"
28 #define ACOMP_STEREO_URI "urn:ardour:a-comp#stereo"
48 ACOMP_STEREO_ATTACK = 0,
52 ACOMP_STEREO_THRESHOLD,
56 ACOMP_STEREO_OUTLEVEL,
57 ACOMP_STEREO_SIDECHAIN,
91 LV2_Inline_Display_Image_Surface surf;
92 cairo_surface_t* display;
93 LV2_Inline_Display* queue_draw;
96 /* ports pointers are only valid during run so we'll
97 * have to cache them for the display, besides
98 * we do want to check for changes
110 instantiate(const LV2_Descriptor* descriptor,
112 const char* bundle_path,
113 const LV2_Feature* const* features)
115 AComp* acomp = (AComp*)calloc(1, sizeof(AComp));
117 for (int i=0; features[i]; ++i) {
119 if (!strcmp(features[i]->URI, LV2_INLINEDISPLAY__queue_draw)) {
120 acomp->queue_draw = (LV2_Inline_Display*) features[i]->data;
126 acomp->old_yl=acomp->old_y1=acomp->old_yg=0.f;
127 acomp->need_expose = true;
129 return (LV2_Handle)acomp;
134 connect_port_mono(LV2_Handle instance,
138 AComp* acomp = (AComp*)instance;
140 switch ((PortIndexMono)port) {
142 acomp->attack = (float*)data;
145 acomp->release = (float*)data;
148 acomp->knee = (float*)data;
151 acomp->ratio = (float*)data;
153 case ACOMP_THRESHOLD:
154 acomp->thresdb = (float*)data;
157 acomp->makeup = (float*)data;
160 acomp->gainr = (float*)data;
163 acomp->outlevel = (float*)data;
165 case ACOMP_SIDECHAIN:
166 acomp->sidechain = (float*)data;
169 acomp->input0 = (float*)data;
172 acomp->sc = (float*)data;
175 acomp->output0 = (float*)data;
181 connect_port_stereo(LV2_Handle instance,
185 AComp* acomp = (AComp*)instance;
187 switch ((PortIndexStereo)port) {
188 case ACOMP_STEREO_ATTACK:
189 acomp->attack = (float*)data;
191 case ACOMP_STEREO_RELEASE:
192 acomp->release = (float*)data;
194 case ACOMP_STEREO_KNEE:
195 acomp->knee = (float*)data;
197 case ACOMP_STEREO_RATIO:
198 acomp->ratio = (float*)data;
200 case ACOMP_STEREO_THRESHOLD:
201 acomp->thresdb = (float*)data;
203 case ACOMP_STEREO_MAKEUP:
204 acomp->makeup = (float*)data;
206 case ACOMP_STEREO_GAINR:
207 acomp->gainr = (float*)data;
209 case ACOMP_STEREO_OUTLEVEL:
210 acomp->outlevel = (float*)data;
212 case ACOMP_STEREO_SIDECHAIN:
213 acomp->sidechain = (float*)data;
215 case ACOMP_STEREO_INPUT0:
216 acomp->input0 = (float*)data;
218 case ACOMP_STEREO_INPUT1:
219 acomp->input1 = (float*)data;
221 case ACOMP_STEREO_SC:
222 acomp->sc = (float*)data;
224 case ACOMP_STEREO_OUTPUT0:
225 acomp->output0 = (float*)data;
227 case ACOMP_STEREO_OUTPUT1:
228 acomp->output1 = (float*)data;
233 // Force already-denormal float value to zero
235 sanitize_denormal(float value) {
236 if (!isnormal(value)) {
244 return (exp(gdb/20.f*log(10.f)));
249 return (20.f*log10(g));
253 activate(LV2_Handle instance)
255 AComp* acomp = (AComp*)instance;
257 *(acomp->gainr) = 0.0f;
258 *(acomp->outlevel) = -45.0f;
259 acomp->old_yl=acomp->old_y1=acomp->old_yg=0.f;
263 run_mono(LV2_Handle instance, uint32_t n_samples)
265 AComp* acomp = (AComp*)instance;
267 const float* const input = acomp->input0;
268 const float* const sc = acomp->sc;
269 float* const output = acomp->output0;
271 float srate = acomp->srate;
272 float width = (6.f * *(acomp->knee)) + 0.01;
274 float attack_coeff = exp(-1000.f/(*(acomp->attack) * srate));
275 float release_coeff = exp(-1000.f/(*(acomp->release) * srate));
278 float lgaininp = 0.f;
280 float Lxg, Lxl, Lyg, Lyl, Ly1;
281 int usesidechain = (*(acomp->sidechain) < 0.5) ? 0 : 1;
286 float ratio = *(acomp->ratio);
287 float thresdb = *(acomp->thresdb);
290 if (acomp->v_knee != *acomp->knee) {
291 acomp->v_knee = *acomp->knee;
292 acomp->need_expose = true;
295 if (acomp->v_ratio != *acomp->ratio) {
296 acomp->v_ratio = *acomp->ratio;
297 acomp->need_expose = true;
300 if (acomp->v_thresdb != *acomp->thresdb) {
301 acomp->v_thresdb = *acomp->thresdb;
302 acomp->need_expose = true;
308 for (i = 0; i < n_samples; i++) {
311 ingain = usesidechain ? fabs(sc0) : fabs(in0);
313 Lxg = (ingain==0.f) ? -160.f : to_dB(ingain);
314 Lxg = sanitize_denormal(Lxg);
316 Lyg = Lxg + (1.f/ratio-1.f)*(Lxg-thresdb+width/2.f)*(Lxg-thresdb+width/2.f)/(2.f*width);
318 if (2.f*(Lxg-thresdb) < -width) {
321 Lyg = thresdb + (Lxg-thresdb)/ratio;
322 Lyg = sanitize_denormal(Lyg);
327 acomp->old_y1 = sanitize_denormal(acomp->old_y1);
328 acomp->old_yl = sanitize_denormal(acomp->old_yl);
329 Ly1 = fmaxf(Lxl, release_coeff * acomp->old_y1+(1.f-release_coeff)*Lxl);
330 Lyl = attack_coeff * acomp->old_yl+(1.f-attack_coeff)*Ly1;
331 Ly1 = sanitize_denormal(Ly1);
332 Lyl = sanitize_denormal(Lyl);
335 Lgain = from_dB(cdb);
337 *(acomp->gainr) = Lyl;
339 if (ingain > in_peak) {
342 lgaininp = in0 * Lgain;
343 output[i] = lgaininp * from_dB(*(acomp->makeup));
345 max = (fabsf(output[i]) > max) ? fabsf(output[i]) : sanitize_denormal(max);
347 // TODO re-use local variables on stack
348 // store values back to acomp at the end of the inner-loop
354 *(acomp->outlevel) = (max < 0.0056f) ? -45.f : to_dB(max);
357 acomp->v_lvl += .1 * (in_peak - acomp->v_lvl); // crude LPF TODO use n_samples/rate TC
358 const float v_lvl_in = (acomp->v_lvl < 0.001f) ? -60.f : to_dB(acomp->v_lvl);
359 const float v_lvl_out = (max < 0.001f) ? -60.f : to_dB(max);
360 if (fabsf (acomp->v_lvl_out - v_lvl_out) >= 1 || fabsf (acomp->v_lvl_in - v_lvl_in) >= 1) {
362 acomp->need_expose = true;
363 acomp->v_lvl_in = v_lvl_in;
364 acomp->v_lvl_out = v_lvl_out - *acomp->makeup;
366 if (acomp->need_expose && acomp->queue_draw) {
367 acomp->need_expose = false;
368 acomp->queue_draw->queue_draw (acomp->queue_draw->handle);
374 run_stereo(LV2_Handle instance, uint32_t n_samples)
376 AComp* acomp = (AComp*)instance;
378 const float* const input0 = acomp->input0;
379 const float* const input1 = acomp->input1;
380 const float* const sc = acomp->sc;
381 float* const output0 = acomp->output0;
382 float* const output1 = acomp->output1;
384 float srate = acomp->srate;
385 float width = (6.f * *(acomp->knee)) + 0.01;
387 float attack_coeff = exp(-1000.f/(*(acomp->attack) * srate));
388 float release_coeff = exp(-1000.f/(*(acomp->release) * srate));
391 float lgaininp = 0.f;
392 float rgaininp = 0.f;
394 float Lxg, Lxl, Lyg, Lyl, Ly1;
395 int usesidechain = (*(acomp->sidechain) < 0.5) ? 0 : 1;
402 float ratio = *(acomp->ratio);
403 float thresdb = *(acomp->thresdb);
406 if (acomp->v_knee != *acomp->knee) {
407 acomp->v_knee = *acomp->knee;
408 acomp->need_expose = true;
411 if (acomp->v_ratio != *acomp->ratio) {
412 acomp->v_ratio = *acomp->ratio;
413 acomp->need_expose = true;
416 if (acomp->v_thresdb != *acomp->thresdb) {
417 acomp->v_thresdb = *acomp->thresdb;
418 acomp->need_expose = true;
424 for (i = 0; i < n_samples; i++) {
428 maxabslr = fmaxf(fabs(in0), fabs(in1));
429 ingain = usesidechain ? fabs(sc0) : maxabslr;
431 Lxg = (ingain==0.f) ? -160.f : to_dB(ingain);
432 Lxg = sanitize_denormal(Lxg);
434 Lyg = Lxg + (1.f/ratio-1.f)*(Lxg-thresdb+width/2.f)*(Lxg-thresdb+width/2.f)/(2.f*width);
436 if (2.f*(Lxg-thresdb) < -width) {
439 Lyg = thresdb + (Lxg-thresdb)/ratio;
440 Lyg = sanitize_denormal(Lyg);
445 acomp->old_y1 = sanitize_denormal(acomp->old_y1);
446 acomp->old_yl = sanitize_denormal(acomp->old_yl);
447 Ly1 = fmaxf(Lxl, release_coeff * acomp->old_y1+(1.f-release_coeff)*Lxl);
448 Lyl = attack_coeff * acomp->old_yl+(1.f-attack_coeff)*Ly1;
449 Ly1 = sanitize_denormal(Ly1);
450 Lyl = sanitize_denormal(Lyl);
453 Lgain = from_dB(cdb);
455 *(acomp->gainr) = Lyl;
457 if (ingain > in_peak) {
460 lgaininp = in0 * Lgain;
461 rgaininp = in1 * Lgain;
462 output0[i] = lgaininp * from_dB(*(acomp->makeup));
463 output1[i] = rgaininp * from_dB(*(acomp->makeup));
465 max = (fmaxf(fabs(output0[i]), fabs(output1[i])) > max) ? fmaxf(fabs(output0[i]), fabs(output1[i])) : sanitize_denormal(max);
467 // TODO re-use local variables on stack
468 // store values back to acomp at the end of the inner-loop
474 *(acomp->outlevel) = (max < 0.0056f) ? -45.f : to_dB(max);
477 acomp->v_lvl += .1 * (in_peak - acomp->v_lvl); // crude LPF TODO use n_samples/rate TC
478 const float v_lvl_in = (acomp->v_lvl < 0.001f) ? -60.f : to_dB(acomp->v_lvl);
479 const float v_lvl_out = (max < 0.001f) ? -60.f : to_dB(max);
480 if (fabsf (acomp->v_lvl_out - v_lvl_out) >= 1 || fabsf (acomp->v_lvl_in - v_lvl_in) >= 1) {
482 acomp->need_expose = true;
483 acomp->v_lvl_in = v_lvl_in;
484 acomp->v_lvl_out = v_lvl_out - *acomp->makeup;
486 if (acomp->need_expose && acomp->queue_draw) {
487 acomp->need_expose = false;
488 acomp->queue_draw->queue_draw (acomp->queue_draw->handle);
494 deactivate(LV2_Handle instance)
500 cleanup(LV2_Handle instance)
503 AComp* acomp = (AComp*)instance;
504 if (acomp->display) {
505 cairo_surface_destroy (acomp->display);
514 #define MIN(A,B) ((A) < (B)) ? (A) : (B)
519 comp_curve (AComp* self, float xg) {
520 const float knee = self->v_knee;
521 const float ratio = self->v_ratio;
522 const float thresdb = self->v_thresdb;
524 const float width = 6.f * knee + 0.01f;
527 if (2.f * (xg - thresdb) < -width) {
529 } else if (2.f * fabs (xg - thresdb) <= width) {
530 yg = xg + (1.f / ratio - 1.f ) * (xg - thresdb + width / 2.f) * (xg - thresdb + width / 2.f) / (2.f * width);
531 } else if (2.f * (xg - thresdb) > width) {
532 yg = thresdb + (xg - thresdb) / ratio;
537 static LV2_Inline_Display_Image_Surface *
538 render_inline (LV2_Handle instance, uint32_t w, uint32_t max_h)
540 AComp* self = (AComp*)instance;
541 uint32_t h = MIN (w, max_h);
543 if (!self->display || self->w != w || self->h != h) {
544 if (self->display) cairo_surface_destroy(self->display);
545 self->display = cairo_image_surface_create (CAIRO_FORMAT_ARGB32, w, h);
550 cairo_t* cr = cairo_create (self->display);
553 cairo_rectangle (cr, 0, 0, w, h);
554 cairo_set_source_rgba (cr, .2, .2, .2, 1.0);
557 cairo_set_line_width(cr, 1.0);
559 // draw grid 10dB steps
560 const double dash1[] = {1, 2};
561 const double dash2[] = {1, 3};
563 cairo_set_line_cap(cr, CAIRO_LINE_CAP_ROUND);
564 cairo_set_dash(cr, dash2, 2, 2);
565 cairo_set_source_rgba (cr, 0.5, 0.5, 0.5, 0.5);
567 for (uint32_t d = 1; d < 6; ++d) {
568 const float x = -.5 + floorf (w * (d * 10.f / 60.f));
569 const float y = -.5 + floorf (h * (d * 10.f / 60.f));
571 cairo_move_to (cr, x, 0);
572 cairo_line_to (cr, x, h);
575 cairo_move_to (cr, 0, y);
576 cairo_line_to (cr, w, y);
579 if (self->v_thresdb < 0) {
580 cairo_set_source_rgba (cr, 0.5, 0.5, 0.5, 1.0);
581 const float y = -.5 + floorf (h * (self->v_thresdb / -60.f));
582 cairo_set_dash(cr, dash1, 2, 2);
583 cairo_move_to (cr, 0, y);
584 cairo_line_to (cr, w, y);
586 cairo_move_to (cr, 0, h);
587 cairo_line_to (cr, w, 0);
594 cairo_set_source_rgba (cr, .8, .8, .8, 1.0);
595 cairo_move_to (cr, 0, h);
597 for (uint32_t x = 0; x < w; ++x) {
599 const float x_db = 60.f * (-1.f + x / (float)w);
600 const float y_db = comp_curve (self, x_db);
601 const float y = h * (y_db / -60.f);
602 cairo_line_to (cr, x, y);
604 cairo_stroke_preserve (cr);
606 cairo_line_to (cr, w, h);
607 cairo_close_path (cr);
610 // draw signal level & reduction/gradient
611 const float top = comp_curve (self, 0);
612 cairo_pattern_t* pat = cairo_pattern_create_linear (0.0, 0.0, 0.0, h);
613 if (top > self->v_thresdb) {
614 cairo_pattern_add_color_stop_rgba (pat, 0.0, 0.8, 0.1, 0.1, 0.5);
615 cairo_pattern_add_color_stop_rgba (pat, top / -60.f, 0.8, 0.1, 0.1, 0.5);
617 if (self->v_knee > 0) {
618 cairo_pattern_add_color_stop_rgba (pat, (self->v_thresdb / -60.f), 0.7, 0.7, 0.2, 0.5);
619 cairo_pattern_add_color_stop_rgba (pat, ((self->v_thresdb - self->v_knee) / -60.f), 0.5, 0.5, 0.5, 0.5);
621 cairo_pattern_add_color_stop_rgba (pat, (self->v_thresdb / -60.f), 0.7, 0.7, 0.2, 0.5);
622 cairo_pattern_add_color_stop_rgba (pat, ((self->v_thresdb - .01) / -60.f), 0.5, 0.5, 0.5, 0.5);
624 cairo_pattern_add_color_stop_rgba (pat, 1.0, 0.5, 0.5, 0.5, 0.5);
626 // maybe cut off at x-position?
627 const float x = w * (self->v_lvl_in + 60) / 60.f;
628 const float y = h * (self->v_lvl_out + 60) / 60.f;
629 cairo_rectangle (cr, 0, h - y, x, y);
630 if (self->v_ratio > 1.0) {
631 cairo_set_source (cr, pat);
633 cairo_set_source_rgba (cr, 0.5, 0.5, 0.5, 0.5);
637 cairo_pattern_destroy (pat); // TODO cache pattern
640 // create RGBA surface
642 cairo_surface_flush (self->display);
643 self->surf.width = cairo_image_surface_get_width (self->display);
644 self->surf.height = cairo_image_surface_get_height (self->display);
645 self->surf.stride = cairo_image_surface_get_stride (self->display);
646 self->surf.data = cairo_image_surface_get_data (self->display);
653 extension_data(const char* uri)
655 static const LV2_Inline_Display_Interface display = { render_inline };
656 if (!strcmp(uri, LV2_INLINEDISPLAY__interface)) {
662 static const LV2_Descriptor descriptor_mono = {
673 static const LV2_Descriptor descriptor_stereo = {
685 const LV2_Descriptor*
686 lv2_descriptor(uint32_t index)
690 return &descriptor_mono;
692 return &descriptor_stereo;