From 46dd5f19538ad06f1b799b14912b91c95e582359 Mon Sep 17 00:00:00 2001
From: Ollie Farley <o.j.d.farley@durham.ac.uk>
Date: Thu, 2 Apr 2026 15:36:36 +0100
Subject: [PATCH] fixed erroneous float32 conversion in phase_covariance which
 was causing numerical issues

---
 aotools/turbulence/turb.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/aotools/turbulence/turb.py b/aotools/turbulence/turb.py
index 9ca662d..0db1a00 100644
--- a/aotools/turbulence/turb.py
+++ b/aotools/turbulence/turb.py
@@ -21,7 +21,7 @@ def phase_covariance(r, r0, L0):
         L0 (float): Outer scale of turbulence in metres
     """
     # Make sure everything is a float to avoid nasty surprises in division!
-    r = numpy.float32(r)
+    r = numpy.float64(r)
     r0 = float(r0)
     L0 = float(L0)