delay.h

Go to the documentation of this file.

/* Calf DSP Library
 * Reusable audio effect classes.
 *
 * Copyright (C) 2001-2007 Krzysztof Foltman
 *
 * This program is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General
 * Public License along with this program; if not, write to the
 * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, 
 * Boston, MA  02110-1301  USA
 */
#ifndef __CALF_DELAY_H
#define __CALF_DELAY_H
 
#include "primitives.h"
#include "buffer.h"
#include "onepole.h"
 
namespace dsp {

template<int N, class T> 
struct simple_delay {
    auto_buffer<N, T> data;
    int pos;
 
    simple_delay() {
        reset();
    }
    void reset() {
        pos = 0;
        for (int i=0; i<N; i++)
            zero(data[i]);
    }    

    inline void put(T idata) {
        data[pos] = idata;
        pos = wrap_around<N>(pos+1);
    }

    template<class U>
    inline void get(U &odata, int delay) {
        assert(delay >= 0 && delay < N);
        int ppos = wrap_around<N>(pos + N - delay);
        odata = data[ppos];
    }
    
    inline T process(T idata, int delay)
    {
        assert(delay >= 0 && delay < N);
        int ppos = wrap_around<N>(pos + N - delay);
        T odata = data[ppos];
        data[pos] = idata;
        pos = wrap_around<N>(pos+1);
        return odata;
    }

    template<class U>
    inline void get_interp(U &odata, int delay, float udelay) {
//        assert(delay >= 0 && delay <= N-1);
        int ppos = wrap_around<N>(pos + N - delay);
        int pppos = wrap_around<N>(ppos + N - 1);
        odata = lerp(data[ppos], data[pppos], udelay);
    }
    
    inline T get_interp_1616(unsigned int delay) {
        float udelay = (float)((delay & 0xFFFF) * (1.0 / 65536.0));
        delay = delay >> 16;
//        assert(delay >= 0 && delay < N-1);
        int ppos = wrap_around<N>(pos + N - delay);
        int pppos = wrap_around<N>(ppos + N - 1);
        return lerp(data[ppos], data[pppos], udelay);
    }
    
    inline T process_comb(T in, unsigned int delay, float fb)
    {
        T old, cur;
        get(old, delay);
        cur = in + fb*old;
        sanitize(cur);
        put(cur);
        return old;
    }
    
    inline T process_comb_lerp16(T in, unsigned int delay, float udelay, float fb)
    {
        T old, cur;
        get_interp(old, delay>>16, dsp::fract16(delay));
        cur = in + fb*old;
        sanitize(cur);
        put(cur);
        return old;
    }
    
    inline T process_allpass_comb(T in, unsigned int delay, float fb)
    {
        T old, cur;
        get(old, delay);
        cur = in + fb*old;
        sanitize(cur);
        put(cur);
        return old - fb * cur;
    }

    inline T process_allpass_comb_lerp16(T in, unsigned int delay, float fb)
    {
        T old, cur;
        get_interp(old, delay>>16, dsp::fract16(delay));
        cur = in + fb*old;
        sanitize(cur);
        put(cur);
        return old - fb * cur;
    }
};
 
};
 
#endif