Working version 1

Author: llaniewski

models/flow/d2q9/Dynamics.c.Rt

@@ -52,6 +52,10 @@ CudaDeviceFunction vector_t getU(){
         u.y += BC[1]*0.5;
     }
 <?R } ?>
+    if (NodeType == NODE_Wall) {
+        u.x = 0;
+        u.y = 0;
+    }
     u.x += GravitationX*0.5;
     u.y += GravitationY*0.5;
 	u.z = 0.0;
@@ -62,7 +66,7 @@ CudaDeviceFunction float2 Color() {
         float2 ret;
         vector_t u = getU();
         ret.x = sqrt(u.x*u.x + u.y*u.y);
-        if (NodeType == NODE_Wall){
+        if (NodeType == NODE_Solid){
                 ret.y = 0;
         } else {
                 ret.y = 1;

src/LatticeContainer.inc.cpp.Rt

@@ -379,7 +379,7 @@ CudaDeviceFunction void NodeToColor( int x, int y, int z, uchar4 *optr )
 if (ISFINITE(l)) {
 
     l = l * 111;
-    /*
+    
     if (               (l <-111)) {r = 255; g = 255; b = 255; }
     if ((l >= -111) && (l < -11)) {r = 255*(-l-11)/100; g = 255; b = 255; }
     if ((l >=  -11) && (l <  -1)) {r = 0; g = (255*(-l-1))/10; b = 255; }
@@ -388,13 +388,7 @@ if (ISFINITE(l)) {
     if ((l >=    1) && (l <  11)) {r = 255; g = 255*(l-1)/10; b = 0; }
     if ((l >=   11) && (l < 111)) {r = 255; g = 255; b = 255*(l-11)/100; }
     if ((l >=  111)             ) {r = 255; g = 255; b = 255; }
-    */
-    if (               (l < -11)) {r = 255; g = 0; b = 0; }
-    if ((l >=  -11) && (l <  -1)) {r = 0; g = 255*(1-(-l-1)/10); b = 0; }
-    if ((l >=   -1) && (l <   0)) {r = 255; g = 255; b = 255*(1+l); }
-    if ((l >=    0) && (l <   1)) {r = 255; g = 255*(1-l); b = 255; }
-    if ((l >=    1) && (l <  11)) {r = 255; g = 0; b = 255*(1-(l-1)/10); }
-    if ((l >=   11)             ) {r = 255; g = 0; b = 0; }
+    
     r=r*w;
     g=g*w + (1-w)*255;
     b=b*w;