combine_channels.h

Go to the documentation of this file.

/******************************************************************************
 Copyright (c) 2012-2016 Institut für Nachrichtentechnik, Universität Rostock
 Copyright (c) 2006-2012 Quality & Usability Lab
                         Deutsche Telekom Laboratories, TU Berlin
 
 Permission is hereby granted, free of charge, to any person obtaining a copy
 of this software and associated documentation files (the "Software"), to deal
 in the Software without restriction, including without limitation the rights
 to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 copies of the Software, and to permit persons to whom the Software is
 furnished to do so, subject to the following conditions:
 
 The above copyright notice and this permission notice shall be included in
 all copies or substantial portions of the Software.
 
 THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
*******************************************************************************/
 
// https://AudioProcessingFramework.github.io/
 
 
#ifndef APF_COMBINE_CHANNELS_H
#define APF_COMBINE_CHANNELS_H
 
#include <vector>
#include <cassert>  // for assert()
#include <stdexcept>  // for std::logic_error
#include <algorithm>  // for std::transform(), std::copy(), std::fill()
 
#include <functional>  // for std::bind()
#include <type_traits>  // for std::remove_reference
 
#include "apf/iterator.h" // for *_iterator, make_*_iterator(), cast_proxy_const
#include "apf/misc.h"  // for CRTP
 
namespace apf
{
 
namespace CombineChannelsResult
{
  enum type
  {
    nothing  = 0,
    constant = 1,
    change   = 2,
    fade_in  = 3,
    fade_out = 4
  };
}
 
template<typename Derived, typename ListProxy, typename Out>
class CombineChannelsBase : public CRTP<Derived>
{
  protected:
    using T = typename std::iterator_traits<typename std::remove_reference<
      ListProxy>::type::value_type::iterator>::value_type;
 
  public:
    template<typename L>
    CombineChannelsBase(L& in, Out& out)
      : _in(in)
      , _out(out)
    {}
 
    template<typename F>
    void process(F f)
    {
      // We pass f by value because this is common in STL-like algorithms.
      // After select() is called, it is passed to case_one() and case_two() as
      // non-const reference to avoid a further copy.
 
      _accumulate = false;
 
      this->derived().before_the_loop();
 
      for (auto& item: _in)
      {
        using namespace CombineChannelsResult;
 
        switch (_selection = f.select(item))
        {
          case nothing:
            continue;  // jump to next list item
 
          case constant:
            this->derived().case_one(item, f);
            break;
 
          case change:
          case fade_in:
          case fade_out:
            this->derived().case_two(item, f);
            break;
 
          default:
            throw std::runtime_error("Predicate must return 0, 1 or 2!");
        }
      }
 
      this->derived().after_the_loop();
 
      if (!_accumulate)
      {
        std::fill(_out.begin(), _out.end(), T());
      }
    }
 
    void before_the_loop() {}
 
    template<typename ItemType, typename F>
    void case_one(const ItemType&, F&)
    {
      throw std::logic_error("CombineChannelsBase: case 1 not implemented!");
    }
 
    template<typename ItemType, typename F>
    void case_two(const ItemType&, F&)
    {
      throw std::logic_error("CombineChannelsBase: case 2 not implemented!");
    }
 
    void after_the_loop() {}
 
  private:
    ListProxy _in;
 
  protected:
    template<typename ItemType>
    void _case_one_copy(const ItemType& item)
    {
      if (_accumulate)
      {
        std::copy(item.begin(), item.end()
            , make_accumulating_iterator(_out.begin()));
      }
      else
      {
        std::copy(item.begin(), item.end(), _out.begin());
        _accumulate = true;
      }
    }
 
    template<typename ItemType, typename FunctionType>
    void _case_one_transform(const ItemType& item, FunctionType& f)
    {
      if (_accumulate)
      {
        std::transform(item.begin(), item.end()
            , make_accumulating_iterator(_out.begin()), f);
      }
      else
      {
        std::transform(item.begin(), item.end(), _out.begin(), f);
        _accumulate = true;
      }
    }
 
    Out& _out;
    CombineChannelsResult::type _selection;
    bool _accumulate;
};
 
template<typename L, typename Out>
class CombineChannelsCopy : public CombineChannelsBase<
                                            CombineChannelsCopy<L, Out>, L, Out>
{
  private:
    using _base = CombineChannelsBase<CombineChannelsCopy<L, Out>, L, Out>;
 
  public:
    CombineChannelsCopy(const L& in, Out& out) : _base(in, out) {}
 
    template<typename ItemType, typename F>
    void case_one(const ItemType& item, F&)
    {
      this->_case_one_copy(item);
    }
 
    // Case 2 is not implemented and shall not be used!
};
 
template<typename L, typename Out>
class CombineChannels: public CombineChannelsBase<
                                                CombineChannels<L, Out>, L, Out>
{
  private:
    using _base = CombineChannelsBase<CombineChannels<L, Out>, L, Out>;
 
  public:
    CombineChannels(const L& in, Out& out) : _base(in, out) {}
 
    template<typename ItemType, typename F>
    void case_one(const ItemType& item, F& f)
    {
      this->_case_one_transform(item, f);
    }
 
    // Case 2 is not implemented and shall not be used!
};
 
template<typename L, typename Out>
class CombineChannelsInterpolation: public CombineChannelsBase<
                                   CombineChannelsInterpolation<L, Out>, L, Out>
{
  private:
    using _base
      = CombineChannelsBase<CombineChannelsInterpolation<L, Out>, L, Out>;
    using typename _base::T;
    using _base::_selection;
    using _base::_accumulate;
    using _base::_out;
 
  public:
    CombineChannelsInterpolation(const L& in, Out& out) : _base(in, out) {}
 
    template<typename ItemType, typename F>
    void case_one(const ItemType& item, F& f)
    {
      this->_case_one_transform(item, f);
    }
 
    template<typename ItemType, typename F>
    void case_two(const ItemType& item, F& f)
    {
      assert(_selection == CombineChannelsResult::change);
 
      if (_accumulate)
      {
        std::transform(item.begin(), item.end(), index_iterator<T>()
            , make_accumulating_iterator(_out.begin()), f);
      }
      else
      {
        std::transform(item.begin(), item.end(), index_iterator<T>()
            , _out.begin(), f);
        _accumulate = true;
      }
    }
};
 
struct fade_out_tag {};
 
template<typename Derived, typename L, typename Out, typename Crossfade>
class CombineChannelsCrossfadeBase : public CombineChannelsBase<Derived, L, Out>
{
  private:
    using _base = CombineChannelsBase<Derived, L, Out>;
    using typename _base::T;
    using _base::_accumulate;
    using _base::_out;
 
  public:
    CombineChannelsCrossfadeBase(const L& in, Out& out, const Crossfade& fade)
      : _base(in, out)
      , _fade_out_buffer(fade.size())
      , _fade_in_buffer(fade.size())
      , _crossfade_data(fade)
    {}
 
    void before_the_loop()
    {
      _accumulate_fade_in = _accumulate_fade_out = false;
    }
 
    void after_the_loop()
    {
      if (_accumulate_fade_out)
      {
        if (_accumulate)
        {
          std::transform(_fade_out_buffer.begin(), _fade_out_buffer.end()
              , _crossfade_data.fade_out_begin()
              , make_accumulating_iterator(_out.begin())
              , std::multiplies<T>());
        }
        else
        {
          std::transform(_fade_out_buffer.begin(), _fade_out_buffer.end()
              , _crossfade_data.fade_out_begin()
              , _out.begin()
              , std::multiplies<T>());
          _accumulate = true;
        }
      }
      if (_accumulate_fade_in)
      {
        if (_accumulate)
        {
          std::transform(_fade_in_buffer.begin(), _fade_in_buffer.end()
              , _crossfade_data.fade_in_begin()
              , make_accumulating_iterator(_out.begin())
              , std::multiplies<T>());
        }
        else
        {
          std::transform(_fade_in_buffer.begin(), _fade_in_buffer.end()
              , _crossfade_data.fade_in_begin()
              , _out.begin()
              , std::multiplies<T>());
          _accumulate = true;
        }
      }
    }
 
  protected:
    bool _accumulate_fade_in, _accumulate_fade_out;
    std::vector<T> _fade_out_buffer, _fade_in_buffer;
 
  private:
    const Crossfade& _crossfade_data;
};
 
template<typename L, typename Out, typename Crossfade>
class CombineChannelsCrossfadeCopy : public CombineChannelsCrossfadeBase<
             CombineChannelsCrossfadeCopy<L, Out, Crossfade>, L, Out, Crossfade>
{
  private:
    using _base = CombineChannelsCrossfadeBase<CombineChannelsCrossfadeCopy<
      L, Out, Crossfade>, L, Out, Crossfade>;
 
    using _base::_fade_out_buffer;
    using _base::_fade_in_buffer;
    using _base::_accumulate_fade_in;
    using _base::_accumulate_fade_out;
    using _base::_selection;
 
  public:
    CombineChannelsCrossfadeCopy(const L& in, Out& out, const Crossfade& fade)
      : _base(in, out, fade)
    {}
 
    template<typename ItemType, typename F>
    void case_one(const ItemType& item, F&)
    {
      this->_case_one_copy(item);
    }
 
    template<typename ItemType, typename F>
    void case_two(ItemType& item, F& f)
    {
      if (_selection != CombineChannelsResult::fade_in)
      {
        if (_accumulate_fade_out)
        {
          std::copy(item.begin(), item.end()
              , make_accumulating_iterator(_fade_out_buffer.begin()));
        }
        else
        {
          std::copy(item.begin(), item.end(), _fade_out_buffer.begin());
          _accumulate_fade_out = true;
        }
      }
      if (_selection != CombineChannelsResult::fade_out)
      {
        f.update();
 
        if (_accumulate_fade_in)
        {
          std::copy(item.begin(), item.end()
              , make_accumulating_iterator(_fade_in_buffer.begin()));
        }
        else
        {
          std::copy(item.begin(), item.end(), _fade_in_buffer.begin());
          _accumulate_fade_in = true;
        }
      }
    }
};
 
template<typename L, typename Out, typename Crossfade>
class CombineChannelsCrossfade : public CombineChannelsCrossfadeBase<
                 CombineChannelsCrossfade<L, Out, Crossfade>, L, Out, Crossfade>
{
  private:
    using _base = CombineChannelsCrossfadeBase<CombineChannelsCrossfade<
      L, Out, Crossfade>, L, Out, Crossfade>;
    using _base::_selection;
    using _base::_accumulate_fade_in;
    using _base::_accumulate_fade_out;
 
  public:
    CombineChannelsCrossfade(const L& in, Out& out, const Crossfade& fade)
      : _base(in, out, fade)
    {}
 
    template<typename ItemType, typename F>
    void case_one(const ItemType& item, F& f)
    {
      this->_case_one_transform(item, f);
    }
 
    template<typename ItemType, typename F>
    void case_two(ItemType& item, F& f)
    {
      if (_selection != CombineChannelsResult::fade_in)
      {
        if (_accumulate_fade_out)
        {
          std::transform(item.begin(), item.end()
              , make_accumulating_iterator(this->_fade_out_buffer.begin())
              , std::bind(f, std::placeholders::_1, fade_out_tag()));
        }
        else
        {
          std::transform(item.begin(), item.end()
              , this->_fade_out_buffer.begin()
              , std::bind(f, std::placeholders::_1, fade_out_tag()));
          _accumulate_fade_out = true;
        }
      }
      if (_selection != CombineChannelsResult::fade_out)
      {
        f.update();
 
        if (_accumulate_fade_in)
        {
          std::transform(item.begin(), item.end()
              , make_accumulating_iterator(this->_fade_in_buffer.begin()), f);
        }
        else
        {
          std::transform(item.begin(), item.end()
              , this->_fade_in_buffer.begin(), f);
          _accumulate_fade_in = true;
        }
      }
    }
};
 
template<typename T>
class raised_cosine_fade
{
  private:
    using iterator_type
      = transform_iterator<index_iterator<T>, math::raised_cosine<T>>;
 
  public:
    using iterator = typename std::vector<T>::const_iterator;
    using reverse_iterator = typename std::vector<T>::const_reverse_iterator;
 
    raised_cosine_fade(size_t block_size)
      : _crossfade_data(
          iterator_type(index_iterator<T>()
            , math::raised_cosine<T>(static_cast<T>(2 * block_size))),
          // block_size + 1 because we also use it in reverse order
          iterator_type(index_iterator<T>(static_cast<T>(block_size + 1))))
      , _size(block_size)
    {}
 
    iterator fade_out_begin() const { return _crossfade_data.begin(); }
    reverse_iterator fade_in_begin() const { return _crossfade_data.rbegin(); }
    size_t size() const { return _size; }
 
  private:
    const std::vector<T> _crossfade_data;
    const size_t _size;
};
 
}  // namespace apf
 
#endif