No longer allow infinite velocities

Author: Ruochun

src/DEM/kT.cpp

@@ -134,7 +134,7 @@ inline void DEMKinematicThread::unpackMyBuffer() {
                         solverScratchSpace);
     // Get the reduced maxVel value
     stateParams.maxVel.toHost();
-    if (*(stateParams.maxVel) > simParams->errOutVel) {
+    if (*(stateParams.maxVel) > simParams->errOutVel || (!std::isfinite(*(stateParams.maxVel)))) {
         DEME_ERROR(
             "System max velocity is %.7g, exceeded max allowance (%.7g).\nIf this velocity is not abnormal and you "
             "want to increase this allowance, use SetErrorOutVelocity before initializing simulation.\nOtherwise, the "

src/kernel/DEMMiscKernels.cu

@@ -42,7 +42,7 @@ __global__ void computeMarginFromAbsv(deme::DEMSimParams* simParams,
     if (ownerID < n) {
         float absv = granData->marginSize[ownerID];
         unsigned int my_family = granData->familyID[ownerID];
-        if (absv > simParams->approxMaxVel) {
+        if (absv > simParams->approxMaxVel || (!isfinite(absv))) {
             absv = simParams->approxMaxVel;
         }
         // User-specified extra margin also needs to be added here. This marginSize is used for bin--sph or bin--tri