debug traces

examples/scalar_advection/scalar_advection.py

@@ -34,8 +34,8 @@ def run(npts=64+1, cfl=0.5):
     d3d=(npts,)*dim
     box = Box(length=(2*pi,)*dim)
     
-    # impl = Implementation.PYTHON
-    impl = Implementation.OPENCL_CODEGEN
+    impl = Implementation.PYTHON
+    # impl = Implementation.OPENCL_CODEGEN
 
     if impl is Implementation.OPENCL_CODEGEN:
         autotuner_config = OpenClKernelAutotunerConfig(
@@ -77,6 +77,7 @@ def run(npts=64+1, cfl=0.5):
     
     dx = dfields[scalar].space_step.min()
     dt = cfl*dx 
+    print 'DT = {}'.format(dt)
     simu = Simulation(start=0.0, end=pi, dt0=0.01)
     
     problem.solve(simu)

hysop/__init__.py

@@ -27,7 +27,7 @@ __KERNEL_DEBUG__   = True
 __ENABLE_LONG_TESTS__ = "OFF" is "ON"
 
 # OpenCL
-__DEFAULT_PLATFORM_ID__ = 1
+__DEFAULT_PLATFORM_ID__ = 0
 __DEFAULT_DEVICE_ID__   = 0
 
 

hysop/backend/device/opencl/operator/directional/advection_dir.py

@@ -242,14 +242,14 @@ class OpenClDirectionalAdvection(DirectionalAdvectionBase, OpenClDirectionalOper
         queue = self.cl_env.default_queue
         dt  = self.precision(simulation.time_step * self.dt_coeff)
         
-        print 'DT= {}'.format(dt)
+        # print 'DT= {}'.format(dt)
         evt = self.advection_kernel_launcher(queue=queue, dt=dt)
-        print 'P'
-        self.dposition.print_with_ghosts()
+        # print 'P'
+        # self.dposition.print_with_ghosts()
         evt = self.remesh_kernel_launcher(queue=queue)
-        print 'S (before accumulation)'
-        self.dadvected_fields_out.values()[0].print_with_ghosts()
+        # print 'S (before accumulation)'
+        # self.dadvected_fields_out.values()[0].print_with_ghosts()
         evt = self.accumulate_and_exchange(queue=queue)
-        print 'S (after accumulation)'
-        self.dadvected_fields_out.values()[0].print_with_ghosts()
+        # print 'S (after accumulation)'
+        # self.dadvected_fields_out.values()[0].print_with_ghosts()
         evt.wait()