Updated D3Q19 UnpackDist() routine so that values are not unpacked into the solid

2014-08-10 10:19:32 -04:00
parent 65c0f04d3e
commit cc6ec0ca2f
2 changed files with 3 additions and 2 deletions
--- a/cpu/D3Q19.cpp
+++ b/cpu/D3Q19.cpp
@@ -48,7 +48,8 @@ extern "C" void UnpackDist(int q, int Cqx, int Cqy, int Cqz, int *list,  int sta
 	//		idx = count;
 	//	}
 //		list[idx] = nn;
-		dist[q*N+nn] = recvbuf[start+idx];
+		// Don't unpack distributions into the solid phase
+		if (dist(q*N+nn) > 0.0)	dist[q*N+nn] = recvbuf[start+idx];
 	}
 }