fill zero velocity values for the solid nodes

2017-11-23 16:07:24 +11:00 · 2017-11-23 16:07:24 +11:00 · 0efeda94cc
commit 0efeda94cc
parent 9c7605403f
1 changed files with 5 additions and 0 deletions
--- a/cpu/D3Q19.cpp
+++ b/cpu/D3Q19.cpp
@ -714,6 +714,11 @@ extern "C" void ScaLBL_D3Q19_Velocity(char *ID, double *disteven, double *distod
 				distodd[q*N+n] = -1.0;
 			}
 			disteven[9*N+n] = -1.0;
+            
+            //For ID[n]<0 - solid nodes
+			vel[n] = 0.0;
+			vel[N+n] = 0.0;
+			vel[2*N+n] = 0.0;
 		}
 	}
 }