fixing another bug in contact points
@@ -416,7 +416,7 @@ class KeplerianOrbit(object):
mask = tt.any(tt.and_(dt >= t_start, dt <= t_end), axis=-1)
result = ifelse(tt.and_(tt.all(tt.eq(M_contact[2], 0)),
- tt.gt(t_end, t_start)),
+ tt.all(tt.gt(t_end, t_start))),
tt.arange(t.size)[mask],
tt.arange(t.size))
