From 3beb37b6854baa9f669bef0c7d8ffddde8c02aa2 Mon Sep 17 00:00:00 2001
From: Cyprien <cyprien.bosserelle@gmail.com>
Date: Wed, 12 Jul 2023 03:13:31 +1200
Subject: [PATCH] Add velocity threshold

---
 Flow_kernel.cu | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Flow_kernel.cu b/Flow_kernel.cu
index 6cdf437..d5b8995 100644
--- a/Flow_kernel.cu
+++ b/Flow_kernel.cu
@@ -1538,7 +1538,7 @@ __global__ void eulerustep(int nx, int ny, DECNUM dx, DECNUM dt, DECNUM g, DECNU
 		}
 		if (ix > 0)
 		{
-			uu[i] = uui[tx][ty];
+			uu[i] = min(uui[tx][ty],16.0f);
 
 		}
 	}
@@ -1594,7 +1594,7 @@ __global__ void eulervstep(int nx, int ny, DECNUM dx, DECNUM dt, DECNUM g, DECNU
 		}
 		if (ix > 0)// && iy>0 && iy<ny)
 		{
-			vv[i] = vvi[tx][ty];
+			vv[i] = min(vvi[tx][ty],16.0f);
 		}//vdvdy[i]=tauby;
 
 	}