fix iono warning

commaai · Aug 25, 2023 · ac5b966 · ac5b966
1 parent a028f6a
commit ac5b966
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/laika/iono.py b/laika/iono.py
@@ -12,7 +12,8 @@
 def get_alpha_beta(rcv_pos, el):
  geocentric_alt = np.linalg.norm(rcv_pos)
  alpha = np.pi/2 + el
- beta = np.arcsin(geocentric_alt*np.sin(alpha)/IPP_ALT)
+ arcsin_arg = geocentric_alt*np.sin(alpha)/IPP_ALT
+ beta = np.arcsin(np.clip(arcsin_arg, -1, 1))
  return alpha, beta
 
 def get_slant_delay(rcv_pos, az, el, sat_pos, time, freq, vertical_delay):