cinelerra-5.1/plugins/reverb/reverb.C

   1
   2 /*
   3  * CINELERRA
   4  * Copyright (C) 2017 Adam Williams <broadcast at earthling dot net>
   5  *
   6  * This program is free software; you can redistribute it and/or modify
   7  * it under the terms of the GNU General Public License as published by
   8  * the Free Software Foundation; either version 2 of the License, or
   9  * (at your option) any later version.
  10  *
  11  * This program is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14  * GNU General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU General Public License
  17  * along with this program; if not, write to the Free Software
  18  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  19  *
  20  */
  21
  22 #include "clip.h"
  23 #include "confirmsave.h"
  24 #include "bchash.h"
  25 #include "bcsignals.h"
  26 #include "errorbox.h"
  27 #include "filexml.h"
  28 #include "language.h"
  29 #include "reverb.h"
  30 #include "reverbwindow.h"
  31 #include "samples.h"
  32 #include "units.h"
  33
  34 #include "vframe.h"
  35
  36 #include <math.h>
  37 #include <string.h>
  38 #include <time.h>
  39 #include <unistd.h>
  40
  41
  42
  43 PluginClient* new_plugin(PluginServer *server)
  44 {
  45         return new Reverb(server);
  46 }
  47
  48
  49
  50 Reverb::Reverb(PluginServer *server)
  51  : PluginAClient(server)
  52 {
  53         srand(time(0));
  54         redo_buffers = 1;       // set to redo buffers before the first render
  55         dsp_in_length = 0;
  56         ref_channels = 0;
  57         ref_offsets = 0;
  58         ref_levels = 0;
  59         ref_lowpass = 0;
  60         dsp_in = 0;
  61         lowpass_in1 = 0;
  62         lowpass_in2 = 0;
  63         initialized = 0;
  64
  65 }
  66
  67 Reverb::~Reverb()
  68 {
  69
  70
  71         if(initialized)
  72         {
  73                 for(int i = 0; i < total_in_buffers; i++)
  74                 {
  75                         delete [] dsp_in[i];
  76                         delete [] ref_channels[i];
  77                         delete [] ref_offsets[i];
  78                         delete [] ref_levels[i];
  79                         delete [] ref_lowpass[i];
  80                         delete [] lowpass_in1[i];
  81                         delete [] lowpass_in2[i];
  82                 }
  83
  84                 delete [] dsp_in;
  85                 delete [] ref_channels;
  86                 delete [] ref_offsets;
  87                 delete [] ref_levels;
  88                 delete [] ref_lowpass;
  89                 delete [] lowpass_in1;
  90                 delete [] lowpass_in2;
  91
  92                 for(int i = 0; i < (smp + 1); i++)
  93                 {
  94                         delete engine[i];
  95                 }
  96                 delete [] engine;
  97                 initialized = 0;
  98         }
  99 }
 100
 101 const char* Reverb::plugin_title() { return N_("Reverb"); }
 102 int Reverb::is_realtime() { return 1; }
 103 int Reverb::is_multichannel() { return 1; }
 104 int Reverb::is_synthesis() { return 1; }
 105
 106 int Reverb::process_realtime(int64_t size,
 107         Samples **input_ptr,
 108         Samples **output_ptr)
 109 {
 110         int64_t new_dsp_length, i, j;
 111         main_in = new double*[total_in_buffers];
 112         main_out = new double*[total_in_buffers];
 113
 114         for(i = 0; i < total_in_buffers; i++)
 115         {
 116                 main_in[i] = input_ptr[i]->get_data();
 117                 main_out[i] = output_ptr[i]->get_data();
 118         }
 119
 120 //printf("Reverb::process_realtime 1\n");
 121         redo_buffers |= load_configuration();
 122
 123 //printf("Reverb::process_realtime 1\n");
 124         if(!config.ref_total)
 125         {
 126                 delete [] main_in;
 127                 delete [] main_out;
 128                 return 0;
 129         }
 130
 131
 132         if(!initialized)
 133         {
 134                 dsp_in = new double*[total_in_buffers];
 135                 ref_channels = new int64_t*[total_in_buffers];
 136                 ref_offsets = new int64_t*[total_in_buffers];
 137                 ref_levels = new double*[total_in_buffers];
 138                 ref_lowpass = new int64_t*[total_in_buffers];
 139                 lowpass_in1 = new double*[total_in_buffers];
 140                 lowpass_in2 = new double*[total_in_buffers];
 141
 142
 143
 144                 for(i = 0; i < total_in_buffers; i++)
 145                 {
 146                         dsp_in[i] = new double[1];
 147                         ref_channels[i] = new int64_t[1];
 148                         ref_offsets[i] = new int64_t[1];
 149                         ref_levels[i] = new double[1];
 150                         ref_lowpass[i] = new int64_t[1];
 151                         lowpass_in1[i] = new double[1];
 152                         lowpass_in2[i] = new double[1];
 153                 }
 154
 155                 engine = new ReverbEngine*[(smp + 1)];
 156                 for(i = 0; i < (smp + 1); i++)
 157                 {
 158                         engine[i] = new ReverbEngine(this);
 159 //printf("Reverb::start_realtime %d\n", Thread::calculate_realtime());
 160 // Realtime priority moved to sound driver
 161 //              engine[i]->set_realtime(realtime_priority);
 162                         engine[i]->start();
 163                 }
 164                 initialized = 1;
 165                 redo_buffers = 1;
 166         }
 167
 168         new_dsp_length = size +
 169                 ((int64_t)config.delay_init + config.ref_length) * project_sample_rate / 1000 + 1;
 170
 171         if(redo_buffers || new_dsp_length != dsp_in_length)
 172         {
 173                 for(i = 0; i < total_in_buffers; i++)
 174                 {
 175                         double *old_dsp = dsp_in[i];
 176                         double *new_dsp = new double[new_dsp_length];
 177                         for(j = 0; j < dsp_in_length && j < new_dsp_length; j++)
 178                                 new_dsp[j] = old_dsp[j];
 179
 180
 181                         for( ; j < new_dsp_length; j++) new_dsp[j] = 0;
 182                         delete [] old_dsp;
 183                         dsp_in[i] = new_dsp;
 184                 }
 185
 186                 dsp_in_length = new_dsp_length;
 187                 redo_buffers = 1;
 188         }
 189 //printf("Reverb::process_realtime 1\n");
 190
 191         if(redo_buffers)
 192         {
 193                 for(i = 0; i < total_in_buffers; i++)
 194                 {
 195                         delete [] ref_channels[i];
 196                         delete [] ref_offsets[i];
 197                         delete [] ref_lowpass[i];
 198                         delete [] ref_levels[i];
 199                         delete [] lowpass_in1[i];
 200                         delete [] lowpass_in2[i];
 201
 202                         ref_channels[i] = new int64_t[config.ref_total + 1];
 203                         ref_offsets[i] = new int64_t[config.ref_total + 1];
 204                         ref_lowpass[i] = new int64_t[config.ref_total + 1];
 205                         ref_levels[i] = new double[config.ref_total + 1];
 206                         lowpass_in1[i] = new double[config.ref_total + 1];
 207                         lowpass_in2[i] = new double[config.ref_total + 1];
 208
 209 // set channels
 210                         ref_channels[i][0] = i;         // primary noise
 211                         ref_channels[i][1] = i;         // first reflection
 212 // set offsets
 213                         ref_offsets[i][0] = 0;
 214                         ref_offsets[i][1] = config.delay_init * project_sample_rate / 1000;
 215 // set levels
 216                         ref_levels[i][0] = db.fromdb(config.level_init);
 217                         ref_levels[i][1] = db.fromdb(config.ref_level1);
 218 // set lowpass
 219                         ref_lowpass[i][0] = -1;     // ignore first noise
 220                         ref_lowpass[i][1] = config.lowpass1;
 221                         lowpass_in1[i][0] = 0;
 222                         lowpass_in2[i][0] = 0;
 223                         lowpass_in1[i][1] = 0;
 224                         lowpass_in2[i][1] = 0;
 225
 226                         int64_t ref_division = config.ref_length * project_sample_rate / 1000 / (config.ref_total + 1);
 227                         for(j = 2; j < config.ref_total + 1; j++)
 228                         {
 229 // set random channels for remaining reflections
 230                                 ref_channels[i][j] = rand() % total_in_buffers;
 231
 232 // set random offsets after first reflection
 233                                 ref_offsets[i][j] = ref_offsets[i][1];
 234                                 if( ref_division > 0 )
 235                                         ref_offsets[i][j] += ref_division * j - (rand() % ref_division) / 2;
 236
 237 // set changing levels
 238                                 ref_levels[i][j] = db.fromdb(config.ref_level1 + (config.ref_level2 - config.ref_level1) / (config.ref_total - 1) * (j - 2));
 239                                 //ref_levels[i][j] /= 100;
 240
 241 // set changing lowpass as linear
 242                                 ref_lowpass[i][j] = (int64_t)(config.lowpass1 + (double)(config.lowpass2 - config.lowpass1) / (config.ref_total - 1) * (j - 2));
 243                                 lowpass_in1[i][j] = 0;
 244                                 lowpass_in2[i][j] = 0;
 245                         }
 246                 }
 247
 248                 redo_buffers = 0;
 249         }
 250 //printf("Reverb::process_realtime 1\n");
 251
 252         for(i = 0; i < total_in_buffers; )
 253         {
 254                 for(j = 0; j < (smp + 1) && (i + j) < total_in_buffers; j++)
 255                 {
 256                         engine[j]->process_overlays(i + j, size);
 257                 }
 258
 259                 for(j = 0; j < (smp + 1) && i < total_in_buffers; j++, i++)
 260                 {
 261                         engine[j]->wait_process_overlays();
 262                 }
 263         }
 264 //printf("Reverb::process_realtime 2 %d %d\n", total_in_buffers, size);
 265
 266         for(i = 0; i < total_in_buffers; i++)
 267         {
 268                 double *current_out = main_out[i];
 269                 double *current_in = dsp_in[i];
 270
 271                 for(j = 0; j < size; j++) current_out[j] = current_in[j];
 272
 273                 int64_t k;
 274                 for(k = 0; j < dsp_in_length; j++, k++) current_in[k] = current_in[j];
 275
 276                 for(; k < dsp_in_length; k++) current_in[k] = 0;
 277         }
 278 //printf("Reverb::process_realtime 2 %d %d\n", total_in_buffers, size);
 279
 280
 281         delete [] main_in;
 282         delete [] main_out;
 283         return 0;
 284 }
 285
 286
 287 NEW_WINDOW_MACRO(Reverb, ReverbWindow)
 288
 289
 290 LOAD_CONFIGURATION_MACRO(Reverb, ReverbConfig)
 291
 292
 293 void Reverb::save_data(KeyFrame *keyframe)
 294 {
 295 //printf("Reverb::save_data 1\n");
 296         FileXML output;
 297 //printf("Reverb::save_data 1\n");
 298
 299 // cause xml file to store data directly in text
 300         output.set_shared_output(keyframe->xbuf);
 301 //printf("Reverb::save_data 1\n");
 302
 303         output.tag.set_title("REVERB");
 304         output.tag.set_property("LEVELINIT", config.level_init);
 305         output.tag.set_property("DELAY_INIT", config.delay_init);
 306         output.tag.set_property("REF_LEVEL1", config.ref_level1);
 307         output.tag.set_property("REF_LEVEL2", config.ref_level2);
 308         output.tag.set_property("REF_TOTAL", config.ref_total);
 309 //printf("Reverb::save_data 1\n");
 310         output.tag.set_property("REF_LENGTH", config.ref_length);
 311         output.tag.set_property("LOWPASS1", config.lowpass1);
 312         output.tag.set_property("LOWPASS2", config.lowpass2);
 313 //printf("Reverb::save_data config.ref_level2 %f\n", config.ref_level2);
 314         output.append_tag();
 315         output.tag.set_title("/REVERB");
 316         output.append_tag();
 317         output.append_newline();
 318         output.terminate_string();
 319 //printf("Reverb::save_data 2\n");
 320 }
 321
 322 void Reverb::read_data(KeyFrame *keyframe)
 323 {
 324         FileXML input;
 325 // cause xml file to read directly from text
 326         input.set_shared_input(keyframe->xbuf);
 327         int result = 0;
 328
 329         result = input.read_tag();
 330
 331         if(!result)
 332         {
 333                 if(input.tag.title_is("REVERB"))
 334                 {
 335                         config.level_init = input.tag.get_property("LEVELINIT", config.level_init);
 336                         config.delay_init = input.tag.get_property("DELAY_INIT", config.delay_init);
 337                         config.ref_level1 = input.tag.get_property("REF_LEVEL1", config.ref_level1);
 338                         config.ref_level2 = input.tag.get_property("REF_LEVEL2", config.ref_level2);
 339                         config.ref_total = input.tag.get_property("REF_TOTAL", config.ref_total);
 340                         config.ref_length = input.tag.get_property("REF_LENGTH", config.ref_length);
 341                         config.lowpass1 = input.tag.get_property("LOWPASS1", config.lowpass1);
 342                         config.lowpass2 = input.tag.get_property("LOWPASS2", config.lowpass2);
 343                 }
 344         }
 345
 346         config.boundaries();
 347 }
 348
 349 void Reverb::update_gui()
 350 {
 351         if(thread)
 352         {
 353                 if(load_configuration())
 354                 {
 355 //printf("Reverb::update_gui %d %d\n", __LINE__, config.ref_length);
 356                         thread->window->lock_window("Reverb::update_gui");
 357                         ((ReverbWindow*)thread->window)->level_init->update(config.level_init);
 358                         ((ReverbWindow*)thread->window)->delay_init->update(config.delay_init);
 359                         ((ReverbWindow*)thread->window)->ref_level1->update(config.ref_level1);
 360                         ((ReverbWindow*)thread->window)->ref_level2->update(config.ref_level2);
 361                         ((ReverbWindow*)thread->window)->ref_total->update(config.ref_total);
 362                         ((ReverbWindow*)thread->window)->ref_length->update(config.ref_length);
 363                         ((ReverbWindow*)thread->window)->lowpass1->update(config.lowpass1);
 364                         ((ReverbWindow*)thread->window)->lowpass2->update(config.lowpass2);
 365                         thread->window->unlock_window();
 366                 }
 367         }
 368 }
 369
 370
 371
 372
 373 // int Reverb::load_from_file(char *path)
 374 // {
 375 //      FILE *in;
 376 //      int result = 0;
 377 //      int length;
 378 //      char string[1024];
 379 //
 380 //      if(in = fopen(path, "rb"))
 381 //      {
 382 //              fseek(in, 0, SEEK_END);
 383 //              length = ftell(in);
 384 //              fseek(in, 0, SEEK_SET);
 385 //              int temp = fread(string, length, 1, in);
 386 //              fclose(in);
 387 // //           read_data(string);
 388 //      }
 389 //      else
 390 //      {
 391 //              perror("fopen:");
 392 // // failed
 393 //              ErrorBox errorbox("");
 394 //              char string[1024];
 395 //              sprintf(string, _("Couldn't open %s."), path);
 396 //              errorbox.create_objects(string);
 397 //              errorbox.run_window();
 398 //              result = 1;
 399 //      }
 400 //
 401 //      return result;
 402 // }
 403 //
 404 // int Reverb::save_to_file(char *path)
 405 // {
 406 //      FILE *out;
 407 //      int result = 0;
 408 //      char string[1024];
 409 //
 410 //      {
 411 // //           ConfirmSave confirm;
 412 // //           result = confirm.test_file("", path);
 413 //      }
 414 //
 415 //      if(!result)
 416 //      {
 417 //              if(out = fopen(path, "wb"))
 418 //              {
 419 // //                   save_data(string);
 420 //                      fwrite(string, strlen(string), 1, out);
 421 //                      fclose(out);
 422 //              }
 423 //              else
 424 //              {
 425 //                      result = 1;
 426 // // failed
 427 //                      ErrorBox errorbox("");
 428 //                      char string[1024];
 429 //                      sprintf(string, _("Couldn't save %s."), path);
 430 //                      errorbox.create_objects(string);
 431 //                      errorbox.run_window();
 432 //                      result = 1;
 433 //              }
 434 //      }
 435 //
 436 //      return result;
 437 // }
 438
 439 ReverbEngine::ReverbEngine(Reverb *plugin)
 440  : Thread(1, 0, 0)
 441 {
 442         this->plugin = plugin;
 443         completed = 0;
 444         input_lock.lock();
 445         output_lock.lock();
 446 }
 447
 448 ReverbEngine::~ReverbEngine()
 449 {
 450         completed = 1;
 451         input_lock.unlock();
 452         join();
 453 }
 454
 455 int ReverbEngine::process_overlays(int output_buffer, int64_t size)
 456 {
 457         this->output_buffer = output_buffer;
 458         this->size = size;
 459         input_lock.unlock();
 460         return 0;
 461 }
 462
 463 int ReverbEngine::wait_process_overlays()
 464 {
 465         output_lock.lock();
 466         return 0;
 467 }
 468
 469 int ReverbEngine::process_overlay(double *in, double *out, double &out1, double &out2, double level, int64_t lowpass, int64_t samplerate, int64_t size)
 470 {
 471 // Modern niquist frequency is 44khz but pot limit is 20khz so can't use
 472 // niquist
 473         if(lowpass == -1 || lowpass >= 20000)
 474         {
 475 // no lowpass filter
 476                 for(int i = 0; i < size; i++) out[i] += in[i] * level;
 477         }
 478         else
 479         {
 480                 double coef = 0.25 * 2.0 * M_PI * (double)lowpass / (double)plugin->project_sample_rate;
 481                 double a = coef * 0.25;
 482                 double b = coef * 0.50;
 483
 484                 for(int i = 0; i < size; i++)
 485                 {
 486                         out2 += a * (3 * out1 + in[i] - out2);
 487                         out2 += b * (out1 + in[i] - out2);
 488                         out2 += a * (out1 + 3 * in[i] - out2);
 489                         out2 += coef * (in[i] - out2);
 490                         out1 = in[i];
 491                         out[i] += out2 * level;
 492                 }
 493         }
 494         return 0;
 495 }
 496
 497 void ReverbEngine::run()
 498 {
 499         int j, i;
 500 //printf("ReverbEngine::run 1 %d\n", calculate_realtime());
 501         while(1)
 502         {
 503                 input_lock.lock();
 504                 if(completed) return;
 505
 506 // Process reverb
 507                 for(i = 0; i < plugin->total_in_buffers; i++)
 508                 {
 509                         for(j = 0; j < plugin->config.ref_total + 1; j++)
 510                         {
 511                                 if(plugin->ref_channels[i][j] == output_buffer)
 512                                         process_overlay(plugin->main_in[i],
 513                                                                 &(plugin->dsp_in[plugin->ref_channels[i][j]][plugin->ref_offsets[i][j]]),
 514                                                                 plugin->lowpass_in1[i][j],
 515                                                                 plugin->lowpass_in2[i][j],
 516                                                                 plugin->ref_levels[i][j],
 517                                                                 plugin->ref_lowpass[i][j],
 518                                                                 plugin->project_sample_rate,
 519                                                                 size);
 520                         }
 521                 }
 522
 523                 output_lock.unlock();
 524         }
 525 }
 526
 527
 528
 529
 530
 531
 532
 533 ReverbConfig::ReverbConfig()
 534 {
 535         level_init = 0;
 536         delay_init = 0;
 537         ref_level1 = -20;
 538         ref_level2 = INFINITYGAIN;
 539         ref_total = 100;
 540         ref_length = 600;
 541         lowpass1 = Freq::tofreq(TOTALFREQS);
 542         lowpass2 = Freq::tofreq(TOTALFREQS);
 543
 544 }
 545
 546 int ReverbConfig::equivalent(ReverbConfig &that)
 547 {
 548         return (EQUIV(level_init, that.level_init) &&
 549                 delay_init == that.delay_init &&
 550                 EQUIV(ref_level1, that.ref_level1) &&
 551                 EQUIV(ref_level2, that.ref_level2) &&
 552                 ref_total == that.ref_total &&
 553                 ref_length == that.ref_length &&
 554                 lowpass1 == that.lowpass1 &&
 555                 lowpass2 == that.lowpass2);
 556 }
 557
 558 void ReverbConfig::copy_from(ReverbConfig &that)
 559 {
 560         level_init = that.level_init;
 561         delay_init = that.delay_init;
 562         ref_level1 = that.ref_level1;
 563         ref_level2 = that.ref_level2;
 564         ref_total = that.ref_total;
 565         ref_length = that.ref_length;
 566         lowpass1 = that.lowpass1;
 567         lowpass2 = that.lowpass2;
 568 }
 569
 570 void ReverbConfig::interpolate(ReverbConfig &prev,
 571         ReverbConfig &next,
 572         int64_t prev_frame,
 573         int64_t next_frame,
 574         int64_t current_frame)
 575 {
 576         level_init = prev.level_init;
 577         delay_init = prev.delay_init;
 578         ref_level1 = prev.ref_level1;
 579         ref_level2 = prev.ref_level2;
 580         ref_total = prev.ref_total;
 581         ref_length = prev.ref_length;
 582         lowpass1 = prev.lowpass1;
 583         lowpass2 = prev.lowpass2;
 584 }
 585
 586 void ReverbConfig::boundaries()
 587 {
 588
 589         CLAMP(level_init, INFINITYGAIN, 0);
 590         CLAMP(delay_init, 0, MAX_DELAY_INIT);
 591         CLAMP(ref_level1, INFINITYGAIN, 0);
 592         CLAMP(ref_level2, INFINITYGAIN, 0);
 593         CLAMP(ref_total, MIN_REFLECTIONS, MAX_REFLECTIONS);
 594         CLAMP(ref_length, 0, MAX_REFLENGTH);
 595         CLAMP(lowpass1, 0, Freq::tofreq(TOTALFREQS));
 596         CLAMP(lowpass2, 0, Freq::tofreq(TOTALFREQS));
 597 }
 598
 599 void ReverbConfig::dump()
 600 {
 601         printf("ReverbConfig::dump %f %jd %f %f %jd %jd %jd %jd\n",
 602                 level_init, delay_init, ref_level1, ref_level2,
 603                 ref_total, ref_length, lowpass1, lowpass2);
 604 }
 605
 606