MSVC uses a slightly smaller value for HUGE
[ardour.git] / libs / backends / dummy / dummy_audiobackend.cc
index 91efd79a495f2226b120374f2e79c60da32cf386..94cd5f440aa157c9af92952f3fe041ef8fdc9b30 100644 (file)
@@ -17,6 +17,7 @@
  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  */
 
+#include <math.h>
 #include <sys/time.h>
 #include <regex.h>
 #include <stdlib.h>
@@ -32,8 +33,9 @@
 #include "dummy_midi_seq.h"
 
 #include "pbd/error.h"
+#include "pbd/compose.h"
 #include "ardour/port_manager.h"
-#include "i18n.h"
+#include "pbd/i18n.h"
 
 using namespace ARDOUR;
 
@@ -112,7 +114,10 @@ DummyAudioBackend::enumerate_devices () const
 {
        if (_device_status.empty()) {
                _device_status.push_back (DeviceStatus (_("Silence"), true));
+               _device_status.push_back (DeviceStatus (_("DC -6dBFS (+.5)"), true));
+               _device_status.push_back (DeviceStatus (_("Demolition"), true));
                _device_status.push_back (DeviceStatus (_("Sine Wave"), true));
+               _device_status.push_back (DeviceStatus (_("Sine Wave 1K, 1/3 Oct"), true));
                _device_status.push_back (DeviceStatus (_("Square Wave"), true));
                _device_status.push_back (DeviceStatus (_("Impulses"), true));
                _device_status.push_back (DeviceStatus (_("Uniform White Noise"), true));
@@ -123,6 +128,7 @@ DummyAudioBackend::enumerate_devices () const
                _device_status.push_back (DeviceStatus (_("Sine Sweep Swell"), true));
                _device_status.push_back (DeviceStatus (_("Square Sweep"), true));
                _device_status.push_back (DeviceStatus (_("Square Sweep Swell"), true));
+               _device_status.push_back (DeviceStatus (_("LTC"), true));
                _device_status.push_back (DeviceStatus (_("Loopback"), true));
        }
        return _device_status;
@@ -178,7 +184,7 @@ DummyAudioBackend::available_output_channel_count (const std::string&) const
 bool
 DummyAudioBackend::can_change_sample_rate_when_running () const
 {
-       return true;
+       return false;
 }
 
 bool
@@ -451,8 +457,6 @@ DummyAudioBackend::_start (bool /*for_latency_measurement*/)
        engine.sample_rate_change (_samplerate);
        engine.buffer_size_change (_samples_per_period);
 
-       _dsp_load_calc.set_max_time (_samplerate, _samples_per_period);
-
        if (engine.reestablish_ports ()) {
                PBD::error << _("DummyAudioBackend: Could not re-establish ports.") << endmsg;
                stop ();
@@ -634,7 +638,7 @@ DummyAudioBackend::my_name () const
 bool
 DummyAudioBackend::available () const
 {
-       return true;
+       return _running;
 }
 
 uint32_t
@@ -727,7 +731,7 @@ DummyAudioBackend::get_ports (
                }
        }
 
-       for (PortIndex::iterator i = _ports.begin (); i != _ports.end (); ++i) {
+       for (PortIndex::const_iterator i = _ports.begin (); i != _ports.end (); ++i) {
                DummyPort* port = *i;
                if ((port->type () == type) && flags == (port->flags () & flags)) {
                        if (!use_regexp || !regexec (&port_regex, port->name ().c_str (), 0, NULL, 0)) {
@@ -805,7 +809,7 @@ DummyAudioBackend::unregister_port (PortEngine::PortHandle port_handle)
                return;
        }
        DummyPort* port = static_cast<DummyPort*>(port_handle);
-       PortIndex::iterator i = _ports.find (static_cast<DummyPort*>(port_handle));
+       PortIndex::iterator i = std::find (_ports.begin(), _ports.end(), static_cast<DummyPort*>(port_handle));
        if (i == _ports.end ()) {
                PBD::error << _("DummyBackend::unregister_port: Failed to find port") << endmsg;
                return;
@@ -831,6 +835,8 @@ DummyAudioBackend::register_system_ports()
                gt = DummyAudioPort::PonyNoise;
        } else if (_device == _("Sine Wave")) {
                gt = DummyAudioPort::SineWave;
+       } else if (_device == _("Sine Wave 1K, 1/3 Oct")) {
+               gt = DummyAudioPort::SineWaveOctaves;
        } else if (_device == _("Square Wave")) {
                gt = DummyAudioPort::SquareWave;
        } else if (_device == _("Impulses")) {
@@ -843,8 +849,14 @@ DummyAudioBackend::register_system_ports()
                gt = DummyAudioPort::SquareSweep;
        } else if (_device == _("Square Sweep Swell")) {
                gt = DummyAudioPort::SquareSweepSwell;
+       } else if (_device == _("LTC")) {
+               gt = DummyAudioPort::LTC;
        } else if (_device == _("Loopback")) {
                gt = DummyAudioPort::Loopback;
+       } else if (_device == _("Demolition")) {
+               gt = DummyAudioPort::Demolition;
+       } else if (_device == _("DC -6dBFS (+.5)")) {
+               gt = DummyAudioPort::DC05;
        } else {
                gt = DummyAudioPort::Silence;
        }
@@ -868,7 +880,10 @@ DummyAudioBackend::register_system_ports()
                if (!p) return -1;
                set_latency_range (p, false, lr);
                _system_inputs.push_back (static_cast<DummyAudioPort*>(p));
-               static_cast<DummyAudioPort*>(p)->setup_generator (gt, _samplerate);
+               std::string name = static_cast<DummyAudioPort*>(p)->setup_generator (gt, _samplerate, i - 1, a_ins);
+               if (!name.empty ()) {
+                       static_cast<DummyAudioPort*>(p)->set_pretty_name (name);
+               }
        }
 
        lr.min = lr.max = _systemic_output_latency;
@@ -891,8 +906,10 @@ DummyAudioBackend::register_system_ports()
                set_latency_range (p, false, lr);
                _system_midi_in.push_back (static_cast<DummyMidiPort*>(p));
                if (_midi_mode == MidiGenerator) {
-                       static_cast<DummyMidiPort*>(p)->setup_generator (i % NUM_MIDI_EVENT_GENERATORS, _samplerate);
-                       static_cast<DummyMidiPort*>(p)->set_pretty_name (DummyMidiData::sequence_names[i % NUM_MIDI_EVENT_GENERATORS]);
+                       std::string name = static_cast<DummyMidiPort*>(p)->setup_generator (i % NUM_MIDI_EVENT_GENERATORS, _samplerate);
+                       if (!name.empty ()) {
+                               static_cast<DummyMidiPort*>(p)->set_pretty_name (name);
+                       }
                }
        }
 
@@ -1065,7 +1082,7 @@ DummyAudioBackend::get_connections (PortEngine::PortHandle port, std::vector<std
 int
 DummyAudioBackend::midi_event_get (
                pframes_t& timestamp,
-               size_t& size, uint8_t** buf, void* port_buffer,
+               size_t& size, uint8_t const** buf, void* port_buffer,
                uint32_t event_index)
 {
        assert (buf && port_buffer);
@@ -1227,7 +1244,7 @@ DummyAudioBackend::n_physical_outputs () const
 {
        int n_midi = 0;
        int n_audio = 0;
-       for (PortIndex::iterator i = _ports.begin (); i != _ports.end (); ++i) {
+       for (PortIndex::const_iterator i = _ports.begin (); i != _ports.end (); ++i) {
                DummyPort* port = *i;
                if (port->is_output () && port->is_physical ()) {
                        switch (port->type ()) {
@@ -1248,7 +1265,7 @@ DummyAudioBackend::n_physical_inputs () const
 {
        int n_midi = 0;
        int n_audio = 0;
-       for (PortIndex::iterator i = _ports.begin (); i != _ports.end (); ++i) {
+       for (PortIndex::const_iterator i = _ports.begin (); i != _ports.end (); ++i) {
                DummyPort* port = *i;
                if (port->is_input () && port->is_physical ()) {
                        switch (port->type ()) {
@@ -1288,6 +1305,7 @@ DummyAudioBackend::main_process_thread ()
        int64_t clock1;
        clock1 = -1;
        while (_running) {
+               const size_t samples_per_period = _samples_per_period;
 
                if (_freewheeling != _freewheel) {
                        _freewheel = _freewheeling;
@@ -1302,17 +1320,17 @@ DummyAudioBackend::main_process_thread ()
                        (*it)->next_period();
                }
 
-               if (engine.process_callback (_samples_per_period)) {
+               if (engine.process_callback (samples_per_period)) {
                        return 0;
                }
-               _processed_samples += _samples_per_period;
+               _processed_samples += samples_per_period;
 
                if (_device == _("Loopback") && _midi_mode != MidiToAudio) {
                        int opn = 0;
                        int opc = _system_outputs.size();
                        for (std::vector<DummyAudioPort*>::const_iterator it = _system_inputs.begin (); it != _system_inputs.end (); ++it, ++opn) {
                                DummyAudioPort* op = _system_outputs[(opn % opc)];
-                               (*it)->fill_wavetable ((const float*)op->get_buffer (_samples_per_period), _samples_per_period);
+                               (*it)->fill_wavetable ((const float*)op->get_buffer (samples_per_period), samples_per_period);
                        }
                }
 
@@ -1331,11 +1349,12 @@ DummyAudioBackend::main_process_thread ()
                        for (std::vector<DummyAudioPort*>::const_iterator it = _system_inputs.begin (); it != _system_inputs.end (); ++it, ++opn) {
                                DummyMidiPort* op = _system_midi_out[(opn % opc)];
                                op->get_buffer(0); // mix-down
-                               (*it)->midi_to_wavetable (op->const_buffer(), _samples_per_period);
+                               (*it)->midi_to_wavetable (op->const_buffer(), samples_per_period);
                        }
                }
 
                if (!_freewheel) {
+                       _dsp_load_calc.set_max_time (_samplerate, samples_per_period);
                        _dsp_load_calc.set_start_timestamp_us (clock1);
                        _dsp_load_calc.set_stop_timestamp_us (_x_get_monotonic_usec());
                        _dsp_load = _dsp_load_calc.get_dsp_load_unbound ();
@@ -1640,17 +1659,57 @@ DummyAudioPort::DummyAudioPort (DummyAudioBackend &b, const std::string& name, P
        , _gen_count2 (0)
        , _pass (false)
        , _rn1 (0)
+       , _ltc (0)
+       , _ltcbuf (0)
 {
        memset (_buffer, 0, sizeof (_buffer));
 }
 
 DummyAudioPort::~DummyAudioPort () {
        free(_wavetable);
+       ltc_encoder_free (_ltc);
+       delete _ltcbuf;
        _wavetable = 0;
+       _ltc = 0;
+       _ltcbuf = 0;
+}
+
+static std::string format_hz (float freq) {
+       std::stringstream ss;
+       if (freq >= 10000) {
+               ss <<  std::setprecision (1) << std::fixed << freq / 1000 << "kHz";
+       } else if (freq >= 1000) {
+               ss <<  std::setprecision (2) << std::fixed << freq / 1000 << "kHz";
+       } else {
+               ss <<  std::setprecision (1) << std::fixed << freq << "Hz";
+       }
+       return ss.str ();
+}
+
+static size_t fit_wave (float freq, float rate, float precision = 0.001) {
+       const size_t max_mult = floor (freq * rate);
+       float minErr = 2;
+       size_t fact = 1;
+       for (size_t i = 1; i < max_mult; ++i) {
+               const float isc = rate * (float)i / freq; // ideal sample count
+               const float rsc = rintf (isc); // rounded sample count
+               const float err = fabsf (isc - rsc);
+               if (err < minErr) {
+                       minErr = err;
+                       fact = i;
+               }
+               if (err < precision) {
+                       break;
+               }
+       }
+       //printf(" FIT %8.1f Hz / %8.1f Hz * %ld = %.0f (err: %e)\n", freq, rate, fact, fact * rate / freq, minErr);
+       return fact;
 }
 
-void DummyAudioPort::setup_generator (GeneratorType const g, float const samplerate)
+std::string
+DummyAudioPort::setup_generator (GeneratorType const g, float const samplerate, int c, int total)
 {
+       std::string name;
        DummyPort::setup_random_number_generator();
        _gen_type = g;
 
@@ -1659,16 +1718,37 @@ void DummyAudioPort::setup_generator (GeneratorType const g, float const sampler
                case PonyNoise:
                case UniformWhiteNoise:
                case GaussianWhiteNoise:
+               case DC05:
                case Silence:
                        break;
+               case Demolition:
+                       _gen_period = 3 * samplerate;
+                       break;
                case KronekerDelta:
                        _gen_period = (5 + randi() % (int)(samplerate / 20.f));
+                       name = "Delta " + format_hz (samplerate / _gen_period);
                        break;
                case SquareWave:
                        _gen_period = (5 + randi() % (int)(samplerate / 20.f)) & ~1;
+                       name = "Square " + format_hz (samplerate / _gen_period);
+                       break;
+               case SineWaveOctaves:
+                       {
+                               const int x = c - floor (((float)total / 2));
+                               float f = powf (2.f, x / 3.f) * 1000.f;
+                               f = std::max (10.f, std::min (samplerate *.5f, f));
+                               const size_t mult = fit_wave (f, samplerate);
+                               _gen_period = rintf ((float)mult * samplerate / f);
+                               name = "Sine " + format_hz (samplerate * mult / (float)_gen_period);
+                               _wavetable = (Sample*) malloc (_gen_period * sizeof(Sample));
+                               for (uint32_t i = 0 ; i < _gen_period; ++i) {
+                                       _wavetable[i] = .12589f * sinf(2.0f * M_PI * (float)mult * (float)i / (float)(_gen_period)); // -18dBFS
+                               }
+                       }
                        break;
                case SineWave:
                        _gen_period = 5 + randi() % (int)(samplerate / 20.f);
+                       name = "Sine " + format_hz (samplerate / _gen_period);
                        _wavetable = (Sample*) malloc (_gen_period * sizeof(Sample));
                        for (uint32_t i = 0 ; i < _gen_period; ++i) {
                                _wavetable[i] = .12589f * sinf(2.0f * M_PI * (float)i / (float)_gen_period); // -18dBFS
@@ -1723,10 +1803,49 @@ void DummyAudioPort::setup_generator (GeneratorType const g, float const sampler
                                }
                        }
                        break;
+               case LTC:
+                       switch (c % 4) {
+                               case 0:
+                                       _ltc = ltc_encoder_create (samplerate, 25, LTC_TV_625_50, 0);
+                                       name = "LTC25";
+                                       break;
+                               case 1:
+                                       _ltc = ltc_encoder_create (samplerate, 30, LTC_TV_1125_60, 0);
+                                       name = "LTC30";
+                                       break;
+                               case 2:
+                                       _ltc = ltc_encoder_create (samplerate, 30001.f / 1001.f, LTC_TV_525_60, 0);
+                                       name = "LTC29df";
+                                       break;
+                               case 3:
+                                       _ltc = ltc_encoder_create (samplerate, 24, LTC_TV_FILM_24, 0);
+                                       name = "LTC24";
+                                       break;
+                       }
+                       _ltc_spd = 1.0;
+                       _ltc_rand = floor((float)c / 4) * .001f;
+                       if (c < 4) {
+                                       name += " (locked)";
+                       } else {
+                                       name += " (varspd)";
+                       }
+                       SMPTETimecode tc;
+                       tc.years = 0;
+                       tc.months = 0;
+                       tc.days = 0;
+                       tc.hours = (3 * (c / 4)) % 24; // XXX
+                       tc.mins = 0;
+                       tc.secs = 0;
+                       tc.frame = 0;
+                       ltc_encoder_set_timecode (_ltc, &tc);
+                                       name += string_compose ("@%1h", (int)tc.hours);
+                       _ltcbuf = new RingBuffer<Sample> (std::max (DummyAudioBackend::max_buffer_size() * 2.f, samplerate));
+                       break;
                case Loopback:
                        _wavetable = (Sample*) malloc (DummyAudioBackend::max_buffer_size() * sizeof(Sample));
                        break;
        }
+       return name;
 }
 
 void DummyAudioPort::midi_to_wavetable (DummyMidiBuffer const * const src, size_t n_samples)
@@ -1781,6 +1900,23 @@ float DummyAudioPort::grandf ()
        return r * x1;
 }
 
+/* inspired by jack-demolition by Steve Harris */
+static const float _demolition[] = {
+        0.0f,             /* special case - 0dbFS white noise */
+        0.0f,             /* zero, may cause denomrals following a signal */
+        0.73 / 1e45,      /* very small - should be denormal when floated */
+        3.7f,             /* arbitrary number > 0dBFS */
+       -4.3f,             /* arbitrary negative number > 0dBFS */
+        4294967395.0f,    /* 2^16 + 100 */
+       -4294967395.0f,
+        3.402823466e+38F, /* HUGE, HUGEVALF, non-inf number */
+        INFINITY,         /* +inf */
+       -INFINITY,         /* -inf */
+       -NAN,              /* -nan */
+        NAN,              /*  nan */
+        0.0f,             /* some silence to check for recovery */
+};
+
 void DummyAudioPort::generate (const pframes_t n_samples)
 {
        Glib::Threads::Mutex::Lock lm (generator_lock);
@@ -1792,6 +1928,30 @@ void DummyAudioPort::generate (const pframes_t n_samples)
                case Silence:
                        memset (_buffer, 0, n_samples * sizeof (Sample));
                        break;
+               case DC05:
+                       for (pframes_t i = 0 ; i < n_samples; ++i) {
+                               _buffer[i] = 0.5f;
+                       }
+                       break;
+               case Demolition:
+                       switch (_gen_count2) {
+                               case 0: // noise
+                                       for (pframes_t i = 0 ; i < n_samples; ++i) {
+                                               _buffer[i] = randf();
+                                       }
+                                       break;
+                               default:
+                                       for (pframes_t i = 0 ; i < n_samples; ++i) {
+                                               _buffer[i] = _demolition [_gen_count2];
+                                       }
+                                       break;
+                       }
+                       _gen_offset += n_samples;
+                       if (_gen_offset > _gen_period) {
+                               _gen_offset = 0;
+                               _gen_count2 = (_gen_count2 + 1) % (sizeof (_demolition) / sizeof (float));
+                       }
+                       break;
                case SquareWave:
                        assert(_gen_period > 0);
                        for (pframes_t i = 0 ; i < n_samples; ++i) {
@@ -1827,8 +1987,10 @@ void DummyAudioPort::generate (const pframes_t n_samples)
                        }
                        break;
                case Loopback:
-                       _gen_period = n_samples; // XXX DummyBackend::_samples_per_period;
+                       memcpy((void*)_buffer, (void*)_wavetable, n_samples * sizeof(Sample));
+                       break;
                case SineWave:
+               case SineWaveOctaves:
                case SineSweep:
                case SquareSweep:
                        assert(_wavetable && _gen_period > 0);
@@ -1883,6 +2045,28 @@ void DummyAudioPort::generate (const pframes_t n_samples)
                                _buffer[i] = _b0 + _b1 + _b2 + white * 0.1848f;
                        }
                        break;
+               case LTC:
+                       while (_ltcbuf->read_space () < n_samples) {
+                               // we should pre-allocate (or add a zero-copy libltc API), whatever.
+                               ltcsnd_sample_t* enc_buf = (ltcsnd_sample_t*) malloc (ltc_encoder_get_buffersize (_ltc) * sizeof (ltcsnd_sample_t));
+                               for (int byteCnt = 0; byteCnt < 10; byteCnt++) {
+                                       if (_ltc_rand != 0.f) {
+                                               _ltc_spd += randf () * _ltc_rand;
+                                               _ltc_spd = std::min (1.5f, std::max (0.5f, _ltc_spd));
+                                       }
+                                       ltc_encoder_encode_byte (_ltc, byteCnt, _ltc_spd);
+                                       const int len = ltc_encoder_get_buffer (_ltc, enc_buf);
+                                       for (int i = 0; i < len; ++i) {
+                                               const float v1 = enc_buf[i] - 128;
+                                               Sample v = v1 * 0.002;
+                                               _ltcbuf->write (&v, 1);
+                                       }
+                               }
+                               ltc_encoder_inc_timecode (_ltc);
+                               free (enc_buf);
+                       }
+                       _ltcbuf->read (_buffer, n_samples);
+                       break;
        }
        _gen_cycle = true;
 }
@@ -1952,13 +2136,15 @@ void DummyMidiPort::set_loopback (DummyMidiBuffer const * const src)
        }
 }
 
-void DummyMidiPort::setup_generator (int seq_id, const float sr)
+std::string
+DummyMidiPort::setup_generator (int seq_id, const float sr)
 {
        DummyPort::setup_random_number_generator();
        _midi_seq_dat = DummyMidiData::sequences[seq_id % NUM_MIDI_EVENT_GENERATORS];
        _midi_seq_spb = sr * .5f; // 120 BPM, beat_time 1.0 per beat.
        _midi_seq_pos = 0;
        _midi_seq_time = 0;
+       return DummyMidiData::sequence_names[seq_id];
 }
 
 void DummyMidiPort::midi_generate (const pframes_t n_samples)
@@ -2019,7 +2205,7 @@ void* DummyMidiPort::get_buffer (pframes_t n_samples)
                                _buffer.push_back (boost::shared_ptr<DummyMidiEvent>(new DummyMidiEvent (**it)));
                        }
                }
-               std::sort (_buffer.begin (), _buffer.end (), MidiEventSorter());
+               std::stable_sort (_buffer.begin (), _buffer.end (), MidiEventSorter());
        } else if (is_output () && is_physical () && is_terminal()) {
                if (!_gen_cycle) {
                        midi_generate(n_samples);