accounted for case where inc_min=-inc_max in priors

kylemede · Sep 23, 2017 · b58e6aa · b58e6aa
1 parent c73598a
commit b58e6aa
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 8 deletions.
diff --git a/ExoSOFT/tools/priors.py b/ExoSOFT/tools/priors.py
@@ -230,13 +230,18 @@ def inc_prior_fn(self, inc):
             ret = self.uniform_fn(inc,7)
         else:
             if inc not in [0.0,90.0,180.0]:
-                mn = np.radians(self.mins_ary[7])
-                mx = np.radians(self.maxs_ary[7])
+                mn = self.mins_ary[7]
+                mx = self.maxs_ary[7]
+                # Account for case where min = -(max), as it causes error in denominator
+                if mn == -1*mx:
+                    mn = mn-0.1
+                mn_rad = np.radians(mn)
+                mx_rad = np.radians(mx)
                 inc_rad = np.radians(inc)
                 if (self.inc_prior == True) or (self.inc_prior == 'sin'):
-                    ret = np.abs(np.sin(inc_rad)) / np.abs(np.cos(mn)-np.cos(mx))
+                    ret = np.abs(np.sin(inc_rad)) / np.abs(np.cos(mn_rad)-np.cos(mx_rad))
                 elif self.inc_prior == 'cos':
-                    ret =  np.abs(np.cos(inc_rad)) / np.abs(np.cos(mn)-np.cos(mx))
+                    ret =  np.abs(np.cos(inc_rad)) / np.abs(np.cos(mn_rad)-np.cos(mx_rad)))
         #if ret==0: ret=-np.inf
         return ret
 

diff --git a/examples/priors.py b/examples/priors.py
@@ -230,13 +230,18 @@ def inc_prior_fn(self, inc):
             ret = self.uniform_fn(inc,7)
         else:
             if inc not in [0.0,90.0,180.0]:
-                mn = np.radians(self.mins_ary[7])
-                mx = np.radians(self.maxs_ary[7])
+                mn = self.mins_ary[7]
+                mx = self.maxs_ary[7]
+                # Account for case where min = -(max), as it causes error in denominator
+                if mn == -1*mx:
+                    mn = mn-0.1
+                mn_rad = np.radians(mn)
+                mx_rad = np.radians(mx)
                 inc_rad = np.radians(inc)
                 if (self.inc_prior == True) or (self.inc_prior == 'sin'):
-                    ret = np.abs(np.sin(inc_rad)) / np.abs(np.cos(mn)-np.cos(mx))
+                    ret = np.abs(np.sin(inc_rad)) / np.abs(np.cos(mn_rad)-np.cos(mx_rad))
                 elif self.inc_prior == 'cos':
-                    ret =  np.abs(np.cos(inc_rad)) / np.abs(np.cos(mn)-np.cos(mx))
+                    ret =  np.abs(np.cos(inc_rad)) / np.abs(np.cos(mn_rad)-np.cos(mx_rad))
         #if ret==0: ret=-np.inf
         return ret