Don't use array omega

jobovy · Nov 17, 2015 · 256c289 · 256c289
1 parent 15e1ff7
commit 256c289
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/galpy/df_src/streamgapdf.py b/galpy/df_src/streamgapdf.py
@@ -158,7 +158,7 @@ def _rewind_angle_impact(self,dangle):
             else: self._caustic_phase= True
         if not self._caustic_phase: # Single-valued everywhere
             guess= dangle
-            dguess= self._dOparMax*self._timpact*1.2
+            dguess= numpy.fabs(self._dOparMax[0])*self._timpact*1.2
             out= optimize.brentq(lambda da: dangle
                                  -self._kick_interpdOpar(da)*self._timpact
                                  -da,