fixed scalar advection

examples/scalar_advection/scalar_advection.py

@@ -12,9 +12,6 @@ from hysop.methods import StrangOrder, Remesh, TimeIntegrator, \
 from hysop.numerics.splitting.strang import StrangSplitting
 from hysop.numerics.odesolvers.runge_kutta import Euler, RK2, RK3, RK4
 
-from hysop.operators import CustomSymbolicOperator
-
-
 IO.set_default_path('/tmp/scalar_advection')
 
 pi  = np.pi
@@ -65,17 +62,8 @@ def run(npts=64+1, cfl=0.5):
                     order=StrangOrder.STRANG_SECOND_ORDER)
     splitting.push_operators(advec)
     
-    from hysop.symbolic.relational import Assignment
-    custom_expr0 = Assignment(scalar.s(), 2*scalar.s()*scalar.s()*velo.s[0]()*velo.s[1]()*velo.s[2]())
-    custom_expr1 = Assignment(velo.s()[1], 12)
-    custom_expr = (custom_expr0, custom_expr1)
-    custom_op = CustomSymbolicOperator(expr=custom_expr,
-            name='custom_op', 
-            variables = {scalar: d3d, velo: d3d})
-    
     problem = Problem(method=method)
-    #problem.insert(splitting)
-    problem.insert(custom_op)
+    problem.insert(splitting)
     
     io_params = IOParams(filename='S0', frequency=1)
     problem.dump_inputs(fields=scalar, io_params=io_params)