diff options
author | David Doukhan <rudk@users.sourceforge.net> | 2011-02-07 13:09:27 +0000 |
---|---|---|
committer | David Doukhan <rudk@users.sourceforge.net> | 2011-02-07 13:09:27 +0000 |
commit | 45dfb5debd4a0d0bc79ad8b6691576e8462d1396 (patch) | |
tree | 7b4d1736586da6ac25a87d5135a357bab291648e /src/spectral_hrtfcont.cpp | |
parent | d69fd2f6b4cbbc538113e107c5d52a5d27d197bd (diff) |
cw_binaural~ code and examples!
Makefile still missing
svn path=/trunk/externals/ddoukhan/cw_binaural~/; revision=14855
Diffstat (limited to 'src/spectral_hrtfcont.cpp')
-rw-r--r-- | src/spectral_hrtfcont.cpp | 67 |
1 files changed, 67 insertions, 0 deletions
diff --git a/src/spectral_hrtfcont.cpp b/src/spectral_hrtfcont.cpp new file mode 100644 index 0000000..0710506 --- /dev/null +++ b/src/spectral_hrtfcont.cpp @@ -0,0 +1,67 @@ +/* + cw_binaural~: a binaural synthesis external for pure data + by David Doukhan - david.doukhan@gmail.com - http://www.limsi.fr/Individu/doukhan + and Anne Sedes - sedes.anne@gmail.com + Copyright (C) 2009-2011 David Doukhan and Anne Sedes + + For more details, see CW_binaural~, a binaural synthesis external for Pure Data + David Doukhan and Anne Sedes, PDCON09 + + + This program is free software: you can redistribute it and/or modify + it under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 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 General Public License for more details. + + You should have received a copy of the GNU General Public License + along with this program. If not, see <http://www.gnu.org/licenses/>. +*/ + + +#include <m_pd.h> //needed for fft implementation +#include "spectral_hrtfcont.hpp" +#include "flyweight_ir_factory.hpp" + +SpectralHrtfCont::SpectralHrtfCont(const ir_key& k): + HrtfCont(k) +{ + ir_key k2 = k; + + // length is twice the impulse response considered + _ir_length = k.length * 2; + // FIXME: IR size verification ??? => power of 2 + + // get the container storing the coresponding temporal impulse response + k2.spectral = false; + HrtfCont* hc = FlyweightIrFactory::instance()->hrtf_set_get(k2); + + // iterate on the container storing the temporal impulse response + for (angle1_cit ie = hc->map_get()->begin(); ie != hc->map_get()->end(); ie++) + for (angle2_cit ia2 = ie->second.begin(); ia2 != ie->second.end(); ia2++) + { + // current azimuth and elevation + const float el = ie->first; + const float az = ia2->first; + + // allocate buffers + float *lbuf = _m[el][az].lbuf = new float[_ir_length]; + float *rbuf = _m[el][az].rbuf = new float[_ir_length]; + + // fill buffers + for (size_t i = 0; i < k.length; ++i) + { + lbuf[i] = ia2->second.lbuf[i]; + rbuf[i] = ia2->second.rbuf[i]; + lbuf[i+k.length] = rbuf[i+k.length] = 0; + } + + // apply fft to the buffers + mayer_realfft(_ir_length, lbuf); + mayer_realfft(_ir_length, rbuf); + } +} |