From 11efd9296e26dc4db68a3081c6f938d7c5f7fa30 Mon Sep 17 00:00:00 2001
From: sjvenditto <sjvenditto@gmail.com>
Date: Tue, 21 Jan 2025 16:12:28 -0500
Subject: [PATCH] bug fix

---
 pynapple/process/spectrum.py | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/pynapple/process/spectrum.py b/pynapple/process/spectrum.py
index 723a4e29..9c7b98f9 100644
--- a/pynapple/process/spectrum.py
+++ b/pynapple/process/spectrum.py
@@ -160,7 +160,9 @@ def compute_power_spectral_density(sig, fs=None, ep=None, full_range=False, n=No
     # subtract from the nyquist frequency to adjust for floating point error in np.fft.fftfreq
     # nyquist freq may occur at negative end of frequencies if N is even
     doubled_freqs = (
-        (fft.index != 0) & (fft.index < (fs / 2 - 1e-6)) & (fft.index > (fs / 2 + 1e-6))
+        (fft.index != 0)  # not 0
+        & (fft.index < (fs / 2 - 1e-6))  # less than positive nyquist freq
+        & (fft.index > (-fs / 2 + 1e-6))  # greater than negative nyquist freq
     )
     psd[doubled_freqs] *= 2