diff --git a/pymc/gp/cov.py b/pymc/gp/cov.py
index 9b04e23977..0eddaea932 100644
--- a/pymc/gp/cov.py
+++ b/pymc/gp/cov.py
@@ -810,7 +810,7 @@ def power_spectral_density_approx(self, J: TensorLike) -> TensorVariable:
         a = 1 / pt.square(self.ls)
         c = pt.where(J > 0, 2, 1)
 
-        q2 = c * pt.iv(J, a) / pt.exp(a)
+        q2 = c * pt.ive(J, a)
 
         return q2