Using faster angular distance

nyctophile-aayu · Apr 25, 2018 · d426c73 · d426c73
1 parent 7ea9583
commit d426c73
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/astromodels/functions/functions_2D.py b/astromodels/functions/functions_2D.py
@@ -1,5 +1,4 @@
 import numpy as np
-from scipy.integrate import quad
 from astropy.coordinates import SkyCoord, ICRS, BaseCoordinateFrame
 from astropy.io import fits
 from astropy import wcs
@@ -147,8 +146,9 @@ def evaluate(self, x, y, lon0, lat0, sigma):
 
         angsep = angular_distance(lon0, lat0, lon, lat)
 
-        return np.power(180 / np.pi, 2) * 1. / (2 * np.pi * sigma ** 2) * np.exp(
-            -0.5 * np.power(angsep, 2) / sigma ** 2)
+        s2 = sigma**2
+
+        return (180 / np.pi)**2 * 1 / (2.0 * np.pi * s2) * np.exp(-0.5 * angsep**2/s2)
 
     def get_boundaries(self):