test

examples/scalar_advection/scalar_advection.py

@@ -35,7 +35,7 @@ def run(npts=64+1, cfl=0.5):
     d3d=(npts,)*dim
     box = Box(length=(2*pi,)*dim)
 
-    impl = Implementation.OPENCL_CODEGEN
+    impl = Implementation.PYTHON#OPENCL_CODEGEN
 
     if impl is Implementation.OPENCL_CODEGEN:
         autotuner_config = OpenClKernelAutotunerConfig(
@@ -78,7 +78,7 @@ def run(npts=64+1, cfl=0.5):
     
     dx = dfields[scalar].space_step.min()
     dt = cfl*dx 
-    simu = Simulation(start=0.0, end=pi, dt0=dt)
+    simu = Simulation(start=0.0, end=pi, dt0=0.01)
     
     problem.solve(simu)
     problem.finalize()

hysop/backend/host/python/operator/directional/advection_dir.py

@@ -75,8 +75,6 @@ class PythonDirectionalAdvection(DirectionalAdvectionBase, HostDirectionalOperat
         velo_mesh = self.dvelocity.mesh
         velo_mesh_iterator = velo_mesh.build_compute_mesh_iterator(cr)
         X0     = velo_mesh.local_compute_coords[-1]
-        print 'X0'
-        print X0
         dx     = velo_mesh.space_step[-1]
         inv_dx = (1.0/dx)
         velo_compute_view = velo_mesh.local_compute_slices
@@ -129,16 +127,10 @@ class PythonDirectionalAdvection(DirectionalAdvectionBase, HostDirectionalOperat
 
         dt  = simulation.time_step
         dt *= self.dt_coeff
-        print 'DT'
-        print float(dt)
         
         self._compute_advection(dt)
-        print 'POSITION'
-        print self.dpos
 
         self._compute_remesh(dt)
-        print 'VORTICITY OUT'
-        self.dadvected_fields_out.values()[0].print_with_ghosts()
 
     def _interp_velocity(self, Vin, Vout, dX, I, Ig, lidx, ridx, inv_dx, is_periodic):
         """
@@ -374,8 +366,6 @@ class PythonDirectionalAdvection(DirectionalAdvectionBase, HostDirectionalOperat
         # exchange and sum remeshed ghosts
         ghost_exchangers = self._out_ghost_exchangers
         for sout in dsoutputs.values():
-            print 'VORTICITY PRE-REDISTRIBUTE'
-            print sout[0]
             buffers      = sout.buffers
             outer_ghosts = self._outer_ghosts[sout]
             inner_ghosts = self._inner_ghosts[sout]