search for: nfft

Hello Jean-Marc, Below are the results that show test_unit_dft passes, but test_unit_mdct fails (only for nfft=480, 960, 1920) Note: Tested on BeagleboneBlack(Cortex-A8) fixed point on branch [1] ./test_unit_dft nfft=32 inverse=0,snr = 88.394372 nfft=32 inverse=1,snr = 93.896470 nfft=128 inverse=0,snr = 89.185895 nfft=128 inverse=1,snr = 93.537021 nfft=256 inverse=0,snr = 88.353151 nfft=256 inverse=1,snr =...

Hello Jean-Marc, Yep, that was it.. with your patch, test_unit_mdct passes for all nfft. So, what you do you suggest the next step here is? Regards, Vish On 8 May 2015 at 12:30, Jean-Marc Valin <jmvalin at jmvalin.ca> wrote: > Hi, > > Can you apply this change to the MDCT test and run it again. See if more > (all) sizes pass. Given the results, I strongly suspect...

...but running make on revision e3187444692195957eb66989622c7b1ad8448b06 fails one of the tests when using fixed point configuration (floating point is ok) on my linux x86. Note that libopus1.1, as extracted from the tar ball, is OK. Specifically, the tests that fail are in celt/tests/test_unit_mdct: nfft=32 inverse=0,snr = 85.341197 nfft=32 inverse=1,snr = 96.285703 nfft=256 inverse=0,snr = 86.726515 nfft=256 inverse=1,snr = 87.308386 nfft=512 inverse=0,snr = 82.577935 nfft=512 inverse=1,snr = 87.103849 nfft=1024 inverse=0,snr = 84.950837 nfft=1024 inverse=1,snr = 87.105776 nfft=2048 inverse=0,snr...

Hello Jean-Marc, **Resending.. not sure why subject got removed earlier** Below are the results that show test_unit_dft passes, but test_unit_mdct fails (only for nfft=480, 960, 1920) Note: Tested on BeagleboneBlack(Cortex-A8) fixed point on branch [1] ./test_unit_dft nfft=32 inverse=0,snr = 88.394372 nfft=32 inverse=1,snr = 93.896470 nfft=128 inverse=0,snr = 89.185895 nfft=128 inverse=1,snr = 93.537021 nfft=256 inverse=0,snr = 88.353151 nfft=256 inverse=1,snr =...

...MDCT test and run it again. See if more (all) sizes pass. Given the results, I strongly suspect an overflow. Jean-Marc On 08/05/15 01:21 PM, Viswanath Puttagunta wrote: > Hello Jean-Marc, > > Below are the results that show test_unit_dft passes, but > test_unit_mdct fails (only for nfft=480, 960, 1920) > Note: Tested on BeagleboneBlack(Cortex-A8) fixed point on branch [1] > > ./test_unit_dft > nfft=32 inverse=0,snr = 88.394372 > nfft=32 inverse=1,snr = 93.896470 > nfft=128 inverse=0,snr = 89.185895 > nfft=128 inverse=1,snr = 93.537021 > nfft=256 inverse=0,...

...revision e3187444692195957eb66989622c7b1ad8448b06 fails one of the tests when using fixed point configuration (floating point is ok) on my linux x86. > Note that libopus1.1, as extracted from the tar ball, is OK. > > Specifically, the tests that fail are in celt/tests/test_unit_mdct: > nfft=32 inverse=0,snr = 85.341197 > nfft=32 inverse=1,snr = 96.285703 > nfft=256 inverse=0,snr = 86.726515 > nfft=256 inverse=1,snr = 87.308386 > nfft=512 inverse=0,snr = 82.577935 > nfft=512 inverse=1,snr = 87.103849 > nfft=1024 inverse=0,snr = 84.950837 > nfft=1024 inverse=1,snr =...

...g CWRS with N=14, K=6... Testing CWRS with N=14, K=7... Testing CWRS with N=14, K=8... Testing CWRS with N=14, K=9... Testing CWRS with N=14, K=10... Testing CWRS with N=14, K=11... Testing CWRS with N=14, K=12... Testing CWRS with N=14, K=13... PASS: cwrs32-test nfft=32 inverse=0,snr = 150.247514 nfft=32 inverse=1,snr = 149.620755 nfft=128 inverse=0,snr = 145.260928 nfft=128 inverse=1,snr = 144.515158 nfft=256 inverse=0,snr = 144.879667 nfft=256 inverse=1,snr = 144.452727 nfft=36 inverse=0,snr = 144.039051 nfft=36 inverse=1,snr = 144...

I have an irregular time series in a Zoo object, and I've been unable to find any way to do an FFT on it. More precisely, I'd like to do an NFFT (non-equispaced / non-uniform time FFT) on the data. The data is timestamped samples from a cheap self-logging accelerometer. The data is weakly regular, with the following characteristics: - short gaps every ~20ms - large gaps every ~200ms - jitter/noise in the timestamp The gaps cover ~10%...

..._fft_cfg_float32_t +#define NE10_FFT_STATE_TYPE_T ne10_fft_state_float32_t +#define NE10_FFT_DESTROY_C2C_TYPE ne10_fft_destroy_c2c_float32 +#define NE10_FFT_CPX_TYPE_T ne10_fft_cpx_float32_t +#define NE10_FFT_C2C_1D_TYPE_NEON ne10_fft_c2c_1d_float32_neon +#else +#define NE10_FFT_ALLOC_C2C_TYPE_NEON(nfft) ne10_fft_alloc_c2c_int32_neon(nfft) +#define NE10_FFT_CFG_TYPE_T ne10_fft_cfg_int32_t +#define NE10_FFT_STATE_TYPE_T ne10_fft_state_int32_t +#define NE10_FFT_DESTROY_C2C_TYPE ne10_fft_destroy_c2c_int32 +#define NE10_FFT_DESTROY_C2C_TYPE ne10_fft_destroy_c2c_int32 +#define NE10_FFT_CPX_TYPE_T ne10_...

..."config.h" +#endif +#endif + +#include <arm_neon.h> +#include <NE10_init.h> +#include <NE10_dsp.h> +#include "../kiss_fft.h" +#include "stack_alloc.h" +#include "os_support.h" +#include "stack_alloc.h" + +#ifdef CUSTOM_MODES + +/* nfft lengths in NE10 that support scaled fft */ +#define NE10_FFTSCALED_SUPPORT_MAX 4 +static const int ne10_fft_scaled_support[NE10_FFTSCALED_SUPPORT_MAX] = { + 480, 240, 120, 60 +}; + +int opus_fft_alloc_arm_float_neon(kiss_fft_state *st) +{ + int i; + size_t memneeded = sizeof(struct arch_fft_s...

..."config.h" +#endif +#endif + +#include <arm_neon.h> +#include <NE10_init.h> +#include <NE10_dsp.h> +#include "../kiss_fft.h" +#include "stack_alloc.h" +#include "os_support.h" +#include "stack_alloc.h" + +#ifdef CUSTOM_MODES + +/* nfft lengths in NE10 that support scaled fft */ +#define NE10_FFTSCALED_SUPPORT_MAX 4 +static const int ne10_fft_scaled_support[NE10_FFTSCALED_SUPPORT_MAX] = { + 480, 240, 120, 60 +}; + +int opus_fft_alloc_arm_float_neon(kiss_fft_state *st) +{ + int i; + size_t memneeded = sizeof(struct arch_fft_s...

..."config.h" +#endif +#endif + +#include <arm_neon.h> +#include <NE10_init.h> +#include <NE10_dsp.h> +#include "../kiss_fft.h" +#include "stack_alloc.h" +#include "os_support.h" +#include "stack_alloc.h" + +#ifdef CUSTOM_MODES + +/* nfft lengths in NE10 that support scaled fft */ +#define NE10_FFTSCALED_SUPPORT_MAX 4 +static const int ne10_fft_scaled_support[NE10_FFTSCALED_SUPPORT_MAX] = { + 480, 240, 120, 60 +}; + +int opus_fft_alloc_arm_float_neon(kiss_fft_state *st) +{ + int i; + size_t memneeded = sizeof(struct arch_fft_s...

Changes from RFC PATCH v3 - Just rebased on tip - For all else, please see notes from RFC PATCH v3 at http://lists.xiph.org/pipermail/opus/2015-March/002902.html - latest wip opus tree/branch https://git.linaro.org/people/viswanath.puttagunta/Ne10.git branch: rfcv4_final_fft_ne10 Viswanath Puttagunta (1): armv7(float): Optimize encode usecase using NE10 library Makefile.am

Changes from RFC PATCH v1: - passing arch parameter explicitly - reduced stack usage by ~3.5K by using scaled NE10 fft version - moved all optimization array functions to arm_celt_map.c - Other cleanups pointed out by Timothy Phil, As you mentioned earlier, could you please address all compile and linker errors/warnings coming out of Ne10 library? You can find my working Ne10 repo at [1] You

Changes from RFC PATCH v2 - fixed compile issue when just compiling for --enable-intrinsics for ARMv7 without NE10 - Notes for NE10: - All compile/link warnings are now in upstream NE10 - Only patch pending upstream in NE10 is the one that needs to add -funsafe-math-optimizations for ARMv7 targets. - Phil Wang @ ARM is working on getting this fixed. - Note that even without

I'm trying to get these cleaned up and landed, but I'm running into some trouble with this patch. Using commit a08b29d88e3c (July 21) of Ne10, I'm seeing test failures for 60-point FFTs: nfft=60 inverse=0,snr = -3.312408 ** poor snr: -3.312408 ** nfft=60 inverse=1,snr = -16.079597 ** poor snr: -16.079597 ** All other sizes tested appear to work fine (84 to 140 dB of SNR). This doesn't match the results you reported on this list in May, where the FFTs appeared to pass tests, but...

...neon(const kiss_fft_state *st, + const kiss_fft_cpx *fin, + kiss_fft_cpx *fout) +{ + ne10_fft_state_float32_t state; + ne10_fft_cfg_float32_t cfg = &state; + VARDECL(ne10_fft_cpx_float32_t, buffer); + SAVE_STACK; + ALLOC(buffer, st->nfft, ne10_fft_cpx_float32_t); + + if (!st->arch_fft->is_supported) { + /* This nfft length (scaled fft) not supported in NE10 */ + opus_ifft_c(st, fin, fout); + } + else { + memcpy((void *)cfg, st->arch_fft->priv, sizeof(ne10_fft_state_float32_t)); + state.buffer =...

...neon(const kiss_fft_state *st, + const kiss_fft_cpx *fin, + kiss_fft_cpx *fout) +{ + ne10_fft_state_float32_t state; + ne10_fft_cfg_float32_t cfg = &state; + VARDECL(ne10_fft_cpx_float32_t, buffer); + SAVE_STACK; + ALLOC(buffer, st->nfft, ne10_fft_cpx_float32_t); + + if (!st->arch_fft->is_supported) { + /* This nfft length (scaled fft) not supported in NE10 */ + opus_ifft_c(st, fin, fout); + } + else { + memcpy((void *)cfg, st->arch_fft->priv, sizeof(ne10_fft_state_float32_t)); + state.buffer =...

Signed-off-by: Viswanath Puttagunta <viswanath.puttagunta at linaro.org> --- celt/tests/test_unit_dft.c | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/celt/tests/test_unit_dft.c b/celt/tests/test_unit_dft.c index 5ea10fb..111c249 100644 --- a/celt/tests/test_unit_dft.c +++ b/celt/tests/test_unit_dft.c @@ -177,8 +177,14 @@ int main(int argc,char ** argv)

Signed-off-by: Viswanath Puttagunta <viswanath.puttagunta at linaro.org> --- celt/tests/test_unit_dft.c | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/celt/tests/test_unit_dft.c b/celt/tests/test_unit_dft.c index e17e26f..f08eb65 100644 --- a/celt/tests/test_unit_dft.c +++ b/celt/tests/test_unit_dft.c @@ -177,8 +177,14 @@ int main(int argc,char ** argv)