diff options
Fix handling of outer boundary cellsstandalone
As we can only use a multiplicative mask to distinguish between cell
types and streaming memory offsets are statically resolved the buffers
have to provide a well-defined padding in both directions. Otherwise
undefined data is accessed which may distort results.
Diffstat (limited to 'standalone_cpp_codegen.py')
-rw-r--r-- | standalone_cpp_codegen.py | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/standalone_cpp_codegen.py b/standalone_cpp_codegen.py index 24a97af..e4dd0b7 100644 --- a/standalone_cpp_codegen.py +++ b/standalone_cpp_codegen.py @@ -12,7 +12,7 @@ lbm = LBM(D3Q19) moments = lbm.moments(optimize = True) collide = lbm.bgk(f_eq = lbm.equilibrium(), tau = 0.6, optimize = True) -geometry = Geometry(64, 64, 64) +geometry = Geometry(32, 32, 32) program_src = Template(filename = str(Path(__file__).parent/'template/standalone.mako')).render( descriptor = lbm.descriptor, @@ -26,7 +26,9 @@ program_src = Template(filename = str(Path(__file__).parent/'template/standalone collide_assignment = collide[1], float_type = 'double', - ccode = sympy.ccode + ccode = sympy.ccode, + + enable_omp_simd = True ) print(program_src) |