Update adaptive.h

Author: TimGoTheCreator

src/dt/adaptive.h

@@ -16,13 +16,15 @@ real computeAdaptiveDt(const std::vector<Particle>& p, real base_dt)
             maxSpeed = speed;
     }
 
+    maxSpeed = std::min(maxSpeed, real(1e4));
+
     // If everything is basically stationary, use base dt
     if (maxSpeed < real(1e-8))
         return base_dt;
 
     // Smaller dt when speeds are high
     real dt = base_dt / (1 + maxSpeed);
-
+     
     // Clamp dt to a reasonable range
     dt = std::max(dt, base_dt * real(0.01)); // never smaller than 1% of base
     dt = std::min(dt, base_dt * real(1.0));  // never larger than base