Fix

Author: xfong

engine/rk23.go

@@ -104,7 +104,7 @@ func (rk *RK23) Step() {
 			rlerr = float64(cuda.MaxAbs(errnorm))
 			fail = 1
 		}
-		if fail == 0 || RelErr < 0.0 || rlerr < RelErr || Dt_si <= MinDt || FixDt != 0 { // mindt check to avoid infinite loop
+		if fail == 0 || RelErr <= 0.0 || rlerr < RelErr || Dt_si <= MinDt || FixDt != 0 { // mindt check to avoid infinite loop
 			// step OK
 			setLastErr(err)
 			setMaxTorque(k4)

engine/rk45dp.go

@@ -118,7 +118,7 @@ func (rk *RK45DP) Step() {
 			rlerr = float64(cuda.MaxAbs(errnorm))
 			fail = 1
 		}
-		if fail == 0 || RelErr < 0.0 || rlerr < RelErr || Dt_si <= MinDt || FixDt != 0 { // mindt check to avoid infinite loop
+		if fail == 0 || RelErr <= 0.0 || rlerr < RelErr || Dt_si <= MinDt || FixDt != 0 { // mindt check to avoid infinite loop
 			// step OK
 			setLastErr(err)
 			setMaxTorque(k7)

engine/rk56.go

@@ -128,7 +128,7 @@ func (rk *RK56) Step() {
 			rlerr = float64(cuda.MaxAbs(errnorm))
 			fail = 1
 		}
-		if fail == 0 || RelErr < 0.0 || rlerr < RelErr || Dt_si <= MinDt || FixDt != 0 { // mindt check to avoid infinite loop
+		if fail == 0 || RelErr <= 0.0 || rlerr < RelErr || Dt_si <= MinDt || FixDt != 0 { // mindt check to avoid infinite loop
 			// step OK
 			setLastErr(err)
 			setMaxTorque(k2)