From 0822625cf870ffaa6851cc6e60dd28d5f8e662fc Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Mon, 10 Aug 2020 20:20:55 +0200
Subject: [PATCH] fixed a NAN in rots::formula_exp

---
 nonisotropic.cpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/nonisotropic.cpp b/nonisotropic.cpp
index 2aec2787..fe3cc6e4 100644
--- a/nonisotropic.cpp
+++ b/nonisotropic.cpp
@@ -2229,7 +2229,7 @@ EX namespace rots {
     if(vel[2] < 0) len = -len;
     
     ld z_part = vel[2]/len;
-    ld x_part = sqrt(1 - z_part * z_part);
+    ld x_part = sqrt(max<ld>(1 - z_part * z_part, 0));
   
     ld SV = stretch::not_squared();
     
@@ -2241,6 +2241,8 @@ EX namespace rots {
     if(sl2 && rparam > 1) {  
       ld cr = 1 / sqrt(rparam*rparam - 1); // *i
       ld sr = rparam * cr;                 // *i
+      
+      if(z_part == 0) cr = 0, sr = 1;
   
       ld z = cr * (K - 1/SV/SV); // *i