lal: fixed missing factor in mkl PSD

18fd13c1 · Duncan Macleod · 302fb352 · 18fd13c1
Verified Commit 18fd13c1 authored Jun 11, 2019 by Duncan Macleod
Hide whitespace changes
Inline Side-by-side

Showing with 1 addition and 0 deletions

lal/src/fft/IntelRealFFT_source.c lal/src/fft/IntelRealFFT_source.c +1 -0

No files found.
--- a/lal/src/fft/IntelRealFFT_source.c
+++ b/lal/src/fft/IntelRealFFT_source.c
@@ -295,6 +295,7 @@ int POWER_SPECTRUM_FUNCTION(REAL_VECTOR_TYPE * spec, const REAL_VECTOR_TYPE * da
        REAL_TYPE re = plan->tmp[2 * k - 1];
        REAL_TYPE im = plan->tmp[2 * k];
        spec->data[k] = re * re + im * im;
+        spec->data[k] *= 2.0;  /* accounts for negative frequency part */
    }

    /* Nyquist frequency */