Only evaluate Keplerian part of DoubleExponentialDiskPotential when t…

…here are points given in the Keplerian part
jobovy · Aug 20, 2019 · 3ac459d · 3ac459d
1 parent f709d04
commit 3ac459d
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/galpy/potential/DoubleExponentialDiskPotential.py b/galpy/potential/DoubleExponentialDiskPotential.py
@@ -135,7 +135,8 @@ def _evaluate(self,R,z,phi=0.,t=0.,dR=0,dphi=0):
                 floatIn= False
             out= nu.empty(len(R))
             indx= (R <= 6.)
-            out[True^indx]= self._kp(R[True^indx],z[True^indx])
+            if nu.sum(True^indx) > 1:
+                out[True^indx]= self._kp(R[True^indx],z[True^indx])
             R4max= nu.copy(R)
             R4max[(R < 1.)]= 1.
             kmax= self._kmaxFac*self._beta