- const char * geoTypeStr[] = { "box", "channel", "pipe", "blocks" };
+ const char * geoTypeStr[] = { "box", "channel", "pipe", "blocks", "fluid" };
printf("# geometry: %d x %d x %d nodes, type %d %s\n", dims[0], dims[1], dims[2], type, geoTypeStr[type]);
printf("# geometry: %d x %d x %d nodes, type %d %s\n", dims[0], dims[1], dims[2], type, geoTypeStr[type]);
- if (blockSize < 0 || blockSize > dimMin / 2) {
- printf("ERROR: block size for geometry must be > 0 and smaller than half of the smalest dimension.\n");
- // TODO: find a better solution for handling errors in here.
- Verify(0);
- }
+ if (blockSize < 0 || blockSize > dimMin / 2) {
+ printf("ERROR: block size for geometry must be > 0 and smaller than half of the smalest dimension.\n");
+ // TODO: find a better solution for handling errors in here.
+ Verify(0);
+ }
- // Number of blocks in x, y, and z direction.
- int nbx = blockSize, nby = blockSize, nbz = blockSize;
+ // Number of blocks in x, y, and z direction.
+ int nbx = blockSize, nby = blockSize, nbz = blockSize;
- if ((x % (2 * nbx)) >= nbx) {
- lattice[L_INDEX_4(dims, x, y, z)] = LAT_CELL_OBSTACLE;
+ if ((x % (2 * nbx)) >= nbx) {
+ lattice[L_INDEX_4(dims, x, y, z)] = LAT_CELL_OBSTACLE;
+ }