Fix logic for switching to Keplerian behavior in DoubleExpDisk

jobovy · Aug 20, 2019 · 1a1f9e0 · 1a1f9e0
1 parent 3ac459d
commit 1a1f9e0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/galpy/potential/DoubleExponentialDiskPotential.py b/galpy/potential/DoubleExponentialDiskPotential.py
@@ -135,7 +135,7 @@ def _evaluate(self,R,z,phi=0.,t=0.,dR=0,dphi=0):
                 floatIn= False
             out= nu.empty(len(R))
             indx= (R <= 6.)
-            if nu.sum(True^indx) > 1:
+            if nu.sum(True^indx) > 0:
                 out[True^indx]= self._kp(R[True^indx],z[True^indx])
             R4max= nu.copy(R)
             R4max[(R < 1.)]= 1.