libs/evoral/src/Sequence.cpp

   1 /* This file is part of Evoral.
   2  * Copyright (C) 2008 Dave Robillard <http://drobilla.net>
   3  * Copyright (C) 2000-2008 Paul Davis
   4  *
   5  * Evoral is free software; you can redistribute it and/or modify it under the
   6  * terms of the GNU General Public License as published by the Free Software
   7  * Foundation; either version 2 of the License, or (at your option) any later
   8  * version.
   9  *
  10  * Evoral is distributed in the hope that it will be useful, but WITHOUT ANY
  11  * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  12  * FOR A PARTICULAR PURPOSE.  See the GNU General Public License for details.
  13  *
  14  * You should have received a copy of the GNU General Public License along
  15  * with this program; if not, write to the Free Software Foundation, Inc.,
  16  * 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  17  */
  18
  19 #define __STDC_LIMIT_MACROS 1
  20 #include <algorithm>
  21 #include <cmath>
  22 #include <iostream>
  23 #include <limits>
  24 #include <stdexcept>
  25 #include <stdint.h>
  26 #include <cstdio>
  27 #include "evoral/Control.hpp"
  28 #include "evoral/ControlList.hpp"
  29 #include "evoral/ControlSet.hpp"
  30 #include "evoral/EventSink.hpp"
  31 #include "evoral/MIDIParameters.hpp"
  32 #include "evoral/Sequence.hpp"
  33 #include "evoral/TypeMap.hpp"
  34 #include "evoral/midi_util.h"
  35
  36 // #define DEBUG_SEQUENCE 1
  37 #ifdef DEBUG_SEQUENCE
  38         #include <boost/format.hpp>
  39         using boost::format;
  40         #define DUMP(x) cerr << (x);
  41 #else
  42         #define DUMP(x)
  43 #endif
  44
  45 using namespace std;
  46
  47 namespace Evoral {
  48
  49 // Read iterator (const_iterator)
  50
  51 template<typename Time>
  52 Sequence<Time>::const_iterator::const_iterator()
  53         : _seq(NULL)
  54         , _is_end(true)
  55         , _control_iter(_control_iters.end())
  56 {
  57         _event = boost::shared_ptr< Event<Time> >(new Event<Time>());
  58 }
  59
  60 template<typename Time>
  61 Sequence<Time>::const_iterator::const_iterator(const Sequence<Time>& seq, Time t)
  62         : _seq(&seq)
  63         , _type(NIL)
  64         , _is_end((t == DBL_MAX) || seq.empty())
  65         , _note_iter(seq.notes().end())
  66         , _sysex_iter(seq.sysexes().end())
  67         , _control_iter(_control_iters.end())
  68 {
  69         DUMP(format("Created Iterator @ %1% (is end: %2%)\n)") % t % _is_end);
  70
  71         if (!_is_end) {
  72                 _lock = seq.read_lock();
  73         } else {
  74                 return;
  75         }
  76
  77         typename Sequence<Time>::ReadLock lock(seq.read_lock());
  78
  79         // Find first note which begins at or after t
  80         _note_iter = seq.note_lower_bound(t);
  81
  82         // Find first sysex event at or after t
  83         for (typename Sequence<Time>::SysExes::const_iterator i = seq.sysexes().begin();
  84                         i != seq.sysexes().end(); ++i) {
  85                 if ((*i)->time() >= t) {
  86                         _sysex_iter = i;
  87                         break;
  88                 }
  89         }
  90         assert(_sysex_iter == seq.sysexes().end() || (*_sysex_iter)->time() >= t);
  91
  92         // Find first control event after t
  93         ControlIterator earliest_control(boost::shared_ptr<ControlList>(), DBL_MAX, 0.0);
  94         _control_iters.reserve(seq._controls.size());
  95         bool   found                  = false;
  96         size_t earliest_control_index = 0;
  97         for (Controls::const_iterator i = seq._controls.begin(); i != seq._controls.end(); ++i) {
  98                 DUMP(format("Iterator: control: %1%\n") % seq._type_map.to_symbol(i->first));
  99                 double x, y;
 100                 bool ret = i->second->list()->rt_safe_earliest_event_unlocked(t, DBL_MAX, x, y, true);
 101                 if (!ret) {
 102                         DUMP(format("Iterator: CC %1% (size %2%) has no events past %3%\n")
 103                                         % i->first.id() % i->second->list()->size() % t);
 104                         continue;
 105                 }
 106
 107                 assert(x >= 0);
 108
 109                 if (y < i->first.min() || y > i->first.max()) {
 110                         cerr << "ERROR: Controller value " << y
 111                                 << " out of range [" << i->first.min() << "," << i->first.max()
 112                                 << "], event ignored" << endl;
 113                         continue;
 114                 }
 115
 116                 DUMP(format("Iterator: CC %1% added (%2%, %3%)\n") % i->first.id() % x % y);
 117
 118                 const ControlIterator new_iter(i->second->list(), x, y);
 119                 _control_iters.push_back(new_iter);
 120
 121                 // Found a new earliest_control
 122                 if (x < earliest_control.x) {
 123                         earliest_control = new_iter;
 124                         earliest_control_index = _control_iters.size() - 1;
 125                         found = true;
 126                 }
 127         }
 128
 129         if (found) {
 130                 _control_iter = _control_iters.begin() + earliest_control_index;
 131                 assert(_control_iter != _control_iters.end());
 132         } else {
 133                 _control_iter = _control_iters.end();
 134         }
 135
 136         // Now find the earliest event overall and point to it
 137         Time earliest_t = t;
 138
 139         if (_note_iter != seq.notes().end()) {
 140                 _type = NOTE_ON;
 141                 earliest_t = (*_note_iter)->time();
 142         }
 143
 144         if (_sysex_iter != seq.sysexes().end()
 145                         && ((*_sysex_iter)->time() < earliest_t || _type == NIL)) {
 146                 _type = SYSEX;
 147                 earliest_t = (*_sysex_iter)->time();
 148         }
 149
 150         if (_control_iter != _control_iters.end()
 151                         && earliest_control.list && earliest_control.x >= t
 152                         && (earliest_control.x < earliest_t || _type == NIL)) {
 153                 _type = CONTROL;
 154                 earliest_t = earliest_control.x;
 155         }
 156
 157         switch (_type) {
 158         case NOTE_ON:
 159                 DUMP(format("Starting at note on event @ %1%\n") % earliest_t);
 160                 _event = boost::shared_ptr< Event<Time> >(
 161                                 new Event<Time>((*_note_iter)->on_event(), true));
 162                 _active_notes.push(*_note_iter);
 163                 break;
 164         case SYSEX:
 165                 DUMP(format("Starting at sysex event @ %1%\n") % earliest_t);
 166                 _event = boost::shared_ptr< Event<Time> >(
 167                                 new Event<Time>(*(*_sysex_iter), true));
 168                 break;
 169         case CONTROL:
 170                 DUMP(format("Starting at control event @ %1%\n") % earliest_t);
 171                 seq.control_to_midi_event(_event, earliest_control);
 172                 break;
 173         default:
 174                 break;
 175         }
 176
 177         if (_type == NIL || !_event || _event->size() == 0) {
 178                 DUMP(format("Starting at end @ %1%\n") % t);
 179                 _type   = NIL;
 180                 _is_end = true;
 181         } else {
 182                 DUMP(printf("New iterator = 0x%x : 0x%x @ %f\n",
 183                             (int)_event->event_type(),
 184                             (int)((MIDIEvent<Time>*)_event.get())->type(),
 185                             _event->time()));
 186                 assert(midi_event_is_valid(_event->buffer(), _event->size()));
 187         }
 188 }
 189
 190 template<typename Time>
 191 Sequence<Time>::const_iterator::~const_iterator()
 192 {
 193 }
 194
 195 template<typename Time>
 196 void
 197 Sequence<Time>::const_iterator::invalidate()
 198 {
 199         while (!_active_notes.empty()) {
 200                 _active_notes.pop();
 201         }
 202         _type = NIL;
 203         _is_end = true;
 204         if (_seq) {
 205                 _note_iter = _seq->notes().end();
 206                 _sysex_iter = _seq->sysexes().end();
 207         }
 208         _control_iter = _control_iters.end();
 209         _lock.reset();
 210 }
 211
 212 template<typename Time>
 213 const typename Sequence<Time>::const_iterator&
 214 Sequence<Time>::const_iterator::operator++()
 215 {
 216         if (_is_end) {
 217                 throw std::logic_error("Attempt to iterate past end of Sequence");
 218         }
 219
 220         DUMP("Sequence::const_iterator++\n");
 221         assert(_event && _event->buffer() && _event->size() > 0);
 222
 223         const MIDIEvent<Time>& ev = *((MIDIEvent<Time>*)_event.get());
 224
 225         if (!(     ev.is_note()
 226                         || ev.is_cc()
 227                         || ev.is_pgm_change()
 228                         || ev.is_pitch_bender()
 229                         || ev.is_channel_pressure()
 230                         || ev.is_sysex()) ) {
 231                 cerr << "WARNING: Unknown event (type " << _type << "): " << hex
 232                         << int(ev.buffer()[0]) << int(ev.buffer()[1]) << int(ev.buffer()[2]) << endl;
 233         }
 234
 235         double x   = 0.0;
 236         double y   = 0.0;
 237         bool   ret = false;
 238
 239         // Increment past current event
 240         switch (_type) {
 241         case NOTE_ON:
 242                 ++_note_iter;
 243                 break;
 244         case NOTE_OFF:
 245                 break;
 246         case CONTROL:
 247                 // Increment current controller iterator
 248                 ret = _control_iter->list->rt_safe_earliest_event_unlocked(
 249                                 _control_iter->x, DBL_MAX, x, y, false);
 250                 assert(!ret || x > _control_iter->x);
 251                 if (ret) {
 252                         _control_iter->x = x;
 253                         _control_iter->y = y;
 254                 } else {
 255                         _control_iter->list.reset();
 256                         _control_iter->x = DBL_MAX;
 257                         _control_iter->y = DBL_MAX;
 258                 }
 259
 260                 // Find the controller with the next earliest event time
 261                 _control_iter = _control_iters.begin();
 262                 for (ControlIterators::iterator i = _control_iters.begin();
 263                                 i != _control_iters.end(); ++i) {
 264                         if (i->x < _control_iter->x) {
 265                                 _control_iter = i;
 266                         }
 267                 }
 268                 break;
 269         case SYSEX:
 270                 ++_sysex_iter;
 271                 break;
 272         default:
 273                 assert(false);
 274         }
 275
 276         // Now find the earliest event overall and point to it
 277         _type = NIL;
 278         Time earliest_t = std::numeric_limits<Time>::max();
 279
 280         // Next earliest note on
 281         if (_note_iter != _seq->notes().end()) {
 282                 _type = NOTE_ON;
 283                 earliest_t = (*_note_iter)->time();
 284         }
 285
 286         // Use the next note off iff it's earlier or the same time as the note on
 287         if (!_seq->percussive() && (!_active_notes.empty())) {
 288                 if (_type == NIL || _active_notes.top()->end_time() <= earliest_t) {
 289                         _type = NOTE_OFF;
 290                         earliest_t = _active_notes.top()->end_time();
 291                 }
 292         }
 293
 294         // Use the next earliest controller iff it's earlier than the note event
 295         if (_control_iter != _control_iters.end() && _control_iter->x != DBL_MAX) {
 296                 if (_type == NIL || _control_iter->x < earliest_t) {
 297                         _type = CONTROL;
 298                         earliest_t = _control_iter->x;
 299                 }
 300         }
 301
 302         // Use the next earliest SysEx iff it's earlier than the controller
 303         if (_sysex_iter != _seq->sysexes().end()) {
 304                 if (_type == NIL || (*_sysex_iter)->time() < earliest_t) {
 305                         _type = SYSEX;
 306                         earliest_t = (*_sysex_iter)->time();
 307                 }
 308         }
 309
 310         // Set event to reflect new position
 311         switch (_type) {
 312         case NOTE_ON:
 313                 DUMP("iterator = note on\n");
 314                 *_event = (*_note_iter)->on_event();
 315                 _active_notes.push(*_note_iter);
 316                 break;
 317         case NOTE_OFF:
 318                 DUMP("iterator = note off\n");
 319                 assert(!_active_notes.empty());
 320                 *_event = _active_notes.top()->off_event();
 321                 _active_notes.pop();
 322                 break;
 323         case CONTROL:
 324                 DUMP("iterator = control\n");
 325                 _seq->control_to_midi_event(_event, *_control_iter);
 326                 break;
 327         case SYSEX:
 328                 DUMP("iterator = sysex\n");
 329                 *_event = *(*_sysex_iter);
 330                 break;
 331         default:
 332                 DUMP("iterator = end\n");
 333                 _is_end = true;
 334         }
 335
 336         assert(_is_end || (_event->size() > 0 && _event->buffer() && _event->buffer()[0] != '\0'));
 337
 338         return *this;
 339 }
 340
 341 template<typename Time>
 342 bool
 343 Sequence<Time>::const_iterator::operator==(const const_iterator& other) const
 344 {
 345         if (_seq != other._seq) {
 346                 return false;
 347         } else if (_is_end || other._is_end) {
 348                 return (_is_end == other._is_end);
 349         } else if (_type != other._type) {
 350                 return false;
 351         } else {
 352                 return (_event == other._event);
 353         }
 354 }
 355
 356 template<typename Time>
 357 typename Sequence<Time>::const_iterator&
 358 Sequence<Time>::const_iterator::operator=(const const_iterator& other)
 359 {
 360         _seq           = other._seq;
 361         _event         = other._event;
 362         _active_notes  = other._active_notes;
 363         _type          = other._type;
 364         _is_end        = other._is_end;
 365         _note_iter     = other._note_iter;
 366         _sysex_iter    = other._sysex_iter;
 367         _control_iters = other._control_iters;
 368
 369         if (other._lock)
 370                 _lock = _seq->read_lock();
 371         else
 372                 _lock.reset();
 373
 374         if (other._control_iter == other._control_iters.end()) {
 375                 _control_iter = _control_iters.end();
 376         } else {
 377                 const size_t index = other._control_iter - other._control_iters.begin();
 378                 _control_iter  = _control_iters.begin() + index;
 379         }
 380
 381         return *this;
 382 }
 383
 384 // Sequence
 385
 386 template<typename Time>
 387 Sequence<Time>::Sequence(const TypeMap& type_map)
 388         : _edited(false)
 389         , _type_map(type_map)
 390         , _writing(false)
 391         , _end_iter(*this, DBL_MAX)
 392         , _percussive(false)
 393         , _lowest_note(127)
 394         , _highest_note(0)
 395 {
 396         DUMP(format("Sequence constructed: %1%\n") % this);
 397         assert(_end_iter._is_end);
 398         assert( ! _end_iter._lock);
 399 }
 400
 401 template<typename Time>
 402 Sequence<Time>::Sequence(const Sequence<Time>& other)
 403         : ControlSet (other)
 404         , _edited(false)
 405         , _type_map(other._type_map)
 406         , _writing(false)
 407         , _end_iter(*this, DBL_MAX)
 408         , _percussive(other._percussive)
 409         , _lowest_note(other._lowest_note)
 410         , _highest_note(other._highest_note)
 411 {
 412         for (typename Notes::const_iterator i = other._notes.begin(); i != other._notes.end(); ++i) {
 413                 boost::shared_ptr<Note<Time> > n (new Note<Time> (**i));
 414                 _notes.insert (n);
 415         }
 416
 417         for (typename SysExes::const_iterator i = other._sysexes.begin(); i != other._sysexes.end(); ++i) {
 418                 boost::shared_ptr<Event<Time> > n (new Event<Time> (**i, true));
 419                 _sysexes.push_back (n);
 420         }
 421
 422         DUMP(format("Sequence copied: %1%\n") % this);
 423         assert(_end_iter._is_end);
 424         assert(! _end_iter._lock);
 425 }
 426
 427 /** Write the controller event pointed to by \a iter to \a ev.
 428  * The buffer of \a ev will be allocated or resized as necessary.
 429  * The event_type of \a ev should be set to the expected output type.
 430  * \return true on success
 431  */
 432 template<typename Time>
 433 bool
 434 Sequence<Time>::control_to_midi_event(
 435                 boost::shared_ptr< Event<Time> >& ev,
 436                 const ControlIterator&            iter) const
 437 {
 438         assert(iter.list.get());
 439         const uint32_t event_type = iter.list->parameter().type();
 440
 441         // initialize the event pointer with a new event, if necessary
 442         if (!ev) {
 443                 ev = boost::shared_ptr< Event<Time> >(new Event<Time>(event_type, 0, 3, NULL, true));
 444         }
 445
 446         uint8_t midi_type = _type_map.parameter_midi_type(iter.list->parameter());
 447         ev->set_event_type(_type_map.midi_event_type(midi_type));
 448         switch (midi_type) {
 449         case MIDI_CMD_CONTROL:
 450                 assert(iter.list.get());
 451                 assert(iter.list->parameter().channel() < 16);
 452                 assert(iter.list->parameter().id() <= INT8_MAX);
 453                 assert(iter.y <= INT8_MAX);
 454
 455                 ev->time() = iter.x;
 456                 ev->realloc(3);
 457                 ev->buffer()[0] = MIDI_CMD_CONTROL + iter.list->parameter().channel();
 458                 ev->buffer()[1] = (uint8_t)iter.list->parameter().id();
 459                 ev->buffer()[2] = (uint8_t)iter.y;
 460                 break;
 461
 462         case MIDI_CMD_PGM_CHANGE:
 463                 assert(iter.list.get());
 464                 assert(iter.list->parameter().channel() < 16);
 465                 assert(iter.y <= INT8_MAX);
 466
 467                 ev->time() = iter.x;
 468                 ev->realloc(2);
 469                 ev->buffer()[0] = MIDI_CMD_PGM_CHANGE + iter.list->parameter().channel();
 470                 ev->buffer()[1] = (uint8_t)iter.y;
 471                 break;
 472
 473         case MIDI_CMD_BENDER:
 474                 assert(iter.list.get());
 475                 assert(iter.list->parameter().channel() < 16);
 476                 assert(iter.y < (1<<14));
 477
 478                 ev->time() = iter.x;
 479                 ev->realloc(3);
 480                 ev->buffer()[0] = MIDI_CMD_BENDER + iter.list->parameter().channel();
 481                 ev->buffer()[1] = uint16_t(iter.y) & 0x7F; // LSB
 482                 ev->buffer()[2] = (uint16_t(iter.y) >> 7) & 0x7F; // MSB
 483                 break;
 484
 485         case MIDI_CMD_CHANNEL_PRESSURE:
 486                 assert(iter.list.get());
 487                 assert(iter.list->parameter().channel() < 16);
 488                 assert(iter.y <= INT8_MAX);
 489
 490                 ev->time() = iter.x;
 491                 ev->realloc(2);
 492                 ev->buffer()[0] = MIDI_CMD_CHANNEL_PRESSURE + iter.list->parameter().channel();
 493                 ev->buffer()[1] = (uint8_t)iter.y;
 494                 break;
 495
 496         default:
 497                 return false;
 498         }
 499
 500         return true;
 501 }
 502
 503 /** Clear all events from the model.
 504  */
 505 template<typename Time>
 506 void
 507 Sequence<Time>::clear()
 508 {
 509         WriteLock lock(write_lock());
 510         _notes.clear();
 511         for (Controls::iterator li = _controls.begin(); li != _controls.end(); ++li)
 512                 li->second->list()->clear();
 513 }
 514
 515 /** Begin a write of events to the model.
 516  *
 517  * If \a mode is Sustained, complete notes with length are constructed as note
 518  * on/off events are received.  Otherwise (Percussive), only note on events are
 519  * stored; note off events are discarded entirely and all contained notes will
 520  * have length 0.
 521  */
 522 template<typename Time>
 523 void
 524 Sequence<Time>::start_write()
 525 {
 526         DUMP(format("%1% : start_write (percussive = %2%)\n") % this % _percussive);
 527         WriteLock lock(write_lock());
 528         _writing = true;
 529         for (int i = 0; i < 16; ++i) {
 530                 _write_notes[i].clear();
 531         }
 532         _dirty_controls.clear();
 533 }
 534
 535 /** Finish a write of events to the model.
 536  *
 537  * If \a delete_stuck is true and the current mode is Sustained, note on events
 538  * that were never resolved with a corresonding note off will be deleted.
 539  * Otherwise they will remain as notes with length 0.
 540  */
 541 template<typename Time>
 542 void
 543 Sequence<Time>::end_write(bool delete_stuck)
 544 {
 545         WriteLock lock(write_lock());
 546
 547         if (!_writing) {
 548                 return;
 549         }
 550
 551         DUMP(format("%1% : end_write (%2% notes)\n") % this % _notes.size());
 552
 553         if (!_percussive && delete_stuck) {
 554                 for (typename Notes::iterator n = _notes.begin(); n != _notes.end() ;) {
 555                         typename Notes::iterator next = n;
 556                         ++next;
 557                         if ((*n)->length() == 0) {
 558                                 cerr << "WARNING: Stuck note lost: " << (*n)->note() << endl;
 559                                 _notes.erase(n);
 560                         }
 561                         n = next;
 562                 }
 563         }
 564
 565         for (int i = 0; i < 16; ++i) {
 566                 if (!_write_notes[i].empty()) {
 567                         cerr << "WARNING: Sequence<Time>::end_write: Channel " << i << " has "
 568                                         << _write_notes[i].size() << " stuck notes" << endl;
 569                 }
 570                 _write_notes[i].clear();
 571         }
 572
 573         for (ControlLists::const_iterator i = _dirty_controls.begin(); i != _dirty_controls.end(); ++i) {
 574                 (*i)->mark_dirty();
 575         }
 576
 577         _writing = false;
 578 }
 579
 580 /** Append \a ev to model.  NOT realtime safe.
 581  *
 582  * Timestamps of events in \a buf are expected to be relative to
 583  * the start of this model (t=0) and MUST be monotonically increasing
 584  * and MUST be >= the latest event currently in the model.
 585  */
 586 template<typename Time>
 587 void
 588 Sequence<Time>::append(const Event<Time>& event)
 589 {
 590         WriteLock lock(write_lock());
 591         _edited = true;
 592
 593         const MIDIEvent<Time>& ev = (const MIDIEvent<Time>&)event;
 594
 595         assert(_notes.empty() || ev.time() >= (*_notes.rbegin())->time());
 596         assert(_writing);
 597
 598         if (!midi_event_is_valid(ev.buffer(), ev.size())) {
 599                 cerr << "WARNING: Sequence ignoring illegal MIDI event" << endl;
 600                 return;
 601         }
 602
 603         if (ev.is_note_on()) {
 604                 append_note_on_unlocked(ev.channel(), ev.time(), ev.note(), ev.velocity());
 605         } else if (ev.is_note_off()) {
 606                 append_note_off_unlocked(ev.channel(), ev.time(), ev.note());
 607         } else if (ev.is_sysex()) {
 608                 append_sysex_unlocked(ev);
 609         } else if (!_type_map.type_is_midi(ev.event_type())) {
 610                 printf("WARNING: Sequence: Unknown event type %X: ", ev.event_type());
 611                 for (size_t i=0; i < ev.size(); ++i) {
 612                         printf("%X ", ev.buffer()[i]);
 613                 }
 614                 printf("\n");
 615         } else if (ev.is_cc()) {
 616                 append_control_unlocked(
 617                                 Evoral::MIDI::ContinuousController(ev.event_type(), ev.channel(), ev.cc_number()),
 618                                 ev.time(), ev.cc_value());
 619         } else if (ev.is_pgm_change()) {
 620                 append_control_unlocked(
 621                                 Evoral::MIDI::ProgramChange(ev.event_type(), ev.channel()),
 622                                 ev.time(), ev.pgm_number());
 623         } else if (ev.is_pitch_bender()) {
 624                 append_control_unlocked(
 625                                 Evoral::MIDI::PitchBender(ev.event_type(), ev.channel()),
 626                                 ev.time(), double(  (0x7F & ev.pitch_bender_msb()) << 7
 627                                         | (0x7F & ev.pitch_bender_lsb()) ));
 628         } else if (ev.is_channel_pressure()) {
 629                 append_control_unlocked(
 630                                 Evoral::MIDI::ChannelPressure(ev.event_type(), ev.channel()),
 631                                 ev.time(), ev.channel_pressure());
 632         } else {
 633                 printf("WARNING: Sequence: Unknown MIDI event type %X\n", ev.type());
 634         }
 635 }
 636
 637 template<typename Time>
 638 void
 639 Sequence<Time>::append_note_on_unlocked(uint8_t chan, Time time, uint8_t note_num, uint8_t velocity)
 640 {
 641         DUMP(format("%1% c=%2% note %3% on @ %4% v=%5%\n")
 642                         % this % (int)chan % (int)note_num % time % (int)velocity);
 643         assert(note_num <= 127);
 644         assert(chan < 16);
 645         assert(_writing);
 646         _edited = true;
 647
 648         if (velocity == 0) {
 649                 append_note_off_unlocked(chan, time, note_num);
 650                 return;
 651         }
 652
 653         if (note_num < _lowest_note)
 654                 _lowest_note = note_num;
 655         if (note_num > _highest_note)
 656                 _highest_note = note_num;
 657
 658         boost::shared_ptr< Note<Time> > new_note(new Note<Time>(chan, time, 0, note_num, velocity));
 659         _notes.insert(new_note);
 660         if (!_percussive) {
 661                 DUMP(format("Sustained: Appending active note on %1% channel %2%\n")
 662                                 % (unsigned)(uint8_t)note_num % chan);
 663                 _write_notes[chan].insert(new_note);
 664         } else {
 665                 DUMP("Percussive: NOT appending active note on\n");
 666         }
 667 }
 668
 669 template<typename Time>
 670 void
 671 Sequence<Time>::append_note_off_unlocked(uint8_t chan, Time time, uint8_t note_num)
 672 {
 673         DUMP(format("%1% c=%2% note %3% off @ %4%\n")
 674                         % this % (int)chan % (int)note_num % time);
 675         assert(note_num <= 127);
 676         assert(chan < 16);
 677         assert(_writing);
 678         _edited = true;
 679
 680         if (_percussive) {
 681                 DUMP("Sequence Ignoring note off (percussive mode)\n");
 682                 return;
 683         }
 684
 685         // TODO: support note off velocity
 686
 687         bool resolved = false;
 688         for (typename WriteNotes::iterator n = _write_notes[chan].begin();
 689                         n != _write_notes[chan].end(); ++n) {
 690                 boost::shared_ptr< Note<Time> > note = *n;
 691                 if (note->note() == note_num) {
 692                         assert(time >= note->time());
 693                         note->set_length(time - note->time());
 694                         _write_notes[chan].erase(n);
 695                         DUMP(format("resolved note, length: %1%\n") % note->length());
 696                         resolved = true;
 697                         break;
 698                 }
 699         }
 700
 701         if (!resolved) {
 702                 cerr << this << " spurious note off chan " << (int)chan
 703                                 << ", note " << (int)note_num << " @ " << time << endl;
 704         }
 705 }
 706
 707 template<typename Time>
 708 void
 709 Sequence<Time>::append_control_unlocked(const Parameter& param, Time time, double value)
 710 {
 711         DUMP(format("%1% %2% @ %3%\t=\t%4% # controls: %5%\n")
 712                         % this % _type_map.to_symbol(param) % time % value % _controls.size());
 713         boost::shared_ptr<Control> c = control(param, true);
 714         c->list()->rt_add(time, value);
 715 }
 716
 717 template<typename Time>
 718 void
 719 Sequence<Time>::append_sysex_unlocked(const MIDIEvent<Time>& ev)
 720 {
 721         #ifdef DEBUG_SEQUENCE
 722         cerr << this << " SysEx @ " << ev.time() << " \t= \t [ " << hex;
 723         for (size_t i=0; i < ev.size(); ++i) {
 724                 cerr << int(ev.buffer()[i]) << " ";
 725         } cerr << "]" << endl;
 726         #endif
 727
 728         boost::shared_ptr<MIDIEvent<Time> > event(new MIDIEvent<Time>(ev, true));
 729         _sysexes.push_back(event);
 730 }
 731
 732 template<typename Time>
 733 bool
 734 Sequence<Time>::contains(const boost::shared_ptr< Note<Time> > note) const
 735 {
 736         ReadLock lock (read_lock());
 737
 738         for (typename Sequence<Time>::Notes::const_iterator i = note_lower_bound(note->time());
 739                         i != _notes.end() && (*i)->time() == note->time(); ++i) {
 740                 if (*i == note) {
 741                         cerr << "Existing note matches: " << *i << endl;
 742                         return true;
 743                 }
 744         }
 745         cerr << "No matching note for " << note << endl;
 746         return false;
 747 }
 748
 749 template<typename Time>
 750 bool
 751 Sequence<Time>::add_note_unlocked(const boost::shared_ptr< Note<Time> > note)
 752 {
 753         DUMP(format("%1% add note %2% @ %3%\n") % this % (int)note->note() % note->time());
 754
 755         for (typename Sequence<Time>::Notes::iterator i = note_lower_bound(note->time());
 756                         i != _notes.end() && (*i)->time() == note->time(); ++i) {
 757                 if (*i == note) {
 758                         return false;
 759                 }
 760         }
 761
 762         _edited = true;
 763         _notes.insert(note);
 764         return true;
 765 }
 766
 767 template<typename Time>
 768 void
 769 Sequence<Time>::remove_note_unlocked(const boost::shared_ptr< const Note<Time> > note)
 770 {
 771         _edited = true;
 772         DUMP(format("%1% remove note %2% @ %3%\n") % this % (int)note->note() % note->time());
 773         for (typename Sequence<Time>::Notes::iterator i = note_lower_bound(note->time());
 774                         i != _notes.end() && (*i)->time() == note->time(); ++i) {
 775                 if (*i == note) {
 776                         _notes.erase(i);
 777                 }
 778         }
 779         cerr << "Unable to find note to erase" << endl;
 780 }
 781
 782 template<typename Time>
 783 void
 784 Sequence<Time>::set_notes (const Sequence<Time>::Notes& n)
 785 {
 786         _notes = n;
 787 }
 788
 789 /** Return the earliest note with time >= t */
 790 template<typename Time>
 791 typename Sequence<Time>::Notes::const_iterator
 792 Sequence<Time>::note_lower_bound (Time t) const
 793 {
 794         boost::shared_ptr< Note<Time> > search_note(new Note<Time>(0, t, 0, 0, 0));
 795         typename Sequence<Time>::Notes::const_iterator i = _notes.lower_bound(search_note);
 796         assert(i == _notes.end() || (*i)->time() >= t);
 797         return i;
 798 }
 799
 800 template<typename Time>
 801 void
 802 Sequence<Time>::get_notes (Notes& n, NoteOperator op, uint8_t val, int chan_mask) const
 803 {
 804         ReadLock lock (read_lock());
 805
 806         for (typename Notes::const_iterator i = _notes.begin(); i != _notes.end(); ++i) {
 807
 808                 if (chan_mask != 0 && !((1<<(*i)->channel()) & chan_mask)) {
 809                         continue;
 810                 }
 811
 812                 switch (op) {
 813                 case PitchEqual:
 814                         if ((*i)->note() == val) {
 815                                 n.insert (*i);
 816                         }
 817                         break;
 818                 case PitchLessThan:
 819                         if ((*i)->note() < val) {
 820                                 n.insert (*i);
 821                         }
 822                         break;
 823                 case PitchLessThanOrEqual:
 824                         if ((*i)->note() <= val) {
 825                                 n.insert (*i);
 826                         }
 827                         break;
 828                 case PitchGreater:
 829                         if ((*i)->note() > val) {
 830                                 n.insert (*i);
 831                         }
 832                         break;
 833                 case PitchGreaterThanOrEqual:
 834                         if ((*i)->note() >= val) {
 835                                 n.insert (*i);
 836                         }
 837                         break;
 838                 case VelocityEqual:
 839                         if ((*i)->velocity() == val) {
 840                                 n.insert (*i);
 841                         }
 842                         break;
 843                 case VelocityLessThan:
 844                         if ((*i)->velocity() < val) {
 845                                 n.insert (*i);
 846                         }
 847                         break;
 848                 case VelocityLessThanOrEqual:
 849                         if ((*i)->velocity() <= val) {
 850                                 n.insert (*i);
 851                         }
 852                         break;
 853                 case VelocityGreater:
 854                         if ((*i)->velocity() > val) {
 855                                 n.insert (*i);
 856                         }
 857                         break;
 858                 case VelocityGreaterThanOrEqual:
 859                         if ((*i)->velocity() >= val) {
 860                                 n.insert (*i);
 861                         }
 862                         break;
 863                 }
 864         }
 865 }
 866
 867 template class Sequence<Evoral::MusicalTime>;
 868
 869 } // namespace Evoral
 870