diff --git a/lbmpy/boundaries/boundaryhandling.py b/lbmpy/boundaries/boundaryhandling.py index b840c6c708c0fd0055bb05798160d4392bcf3fce..fc1fae4a5e04d09ebdb0a671f00a1294e713a839 100644 --- a/lbmpy/boundaries/boundaryhandling.py +++ b/lbmpy/boundaries/boundaryhandling.py @@ -52,11 +52,12 @@ class LatticeBoltzmannBoundaryHandling(BoundaryHandling): if boundary_obj not in self._boundary_object_to_boundary_info: sym_index_field = Field.create_generic('indexField', spatial_dimensions=1, dtype=numpy_data_type_for_boundary_object(boundary_obj, self.dim)) - kernels = [self._create_boundary_kernel( - self._data_handling.fields[self._field_name], sym_index_field, boundary_obj, Timestep.EVEN).compile(), - self._create_boundary_kernel( - self._data_handling.fields[self._field_name], sym_index_field, boundary_obj, - Timestep.ODD).compile()] + + ast_even = self._create_boundary_kernel(self._data_handling.fields[self._field_name], sym_index_field, + boundary_obj, Timestep.EVEN) + ast_odd = self._create_boundary_kernel(self._data_handling.fields[self._field_name], sym_index_field, + boundary_obj, Timestep.ODD) + kernels = [ast_even.compile(), ast_odd.compile()] if flag is None: flag = self.flag_interface.reserve_next_flag() boundary_info = self.InplaceStreamingBoundaryInfo(self, boundary_obj, flag, kernels)