ghost exchangers

examples/scalar_advection/scalar_advection.py

@@ -20,6 +20,8 @@ def compute(args):
     def init_velocity(data, coords):
         for i,d in enumerate(data):
             d[...] = args.velocity[::-1][i]
+        # data[0][...] = np.arange(data[0].shape[0])[:,None]
+        # data[1][...] = -np.arange(data[0].shape[0])[:,None]
 
     ## Function to compute initial scalar values
     def init_scalar(data, coords):
@@ -101,6 +103,8 @@ def compute(args):
     dfields = problem.input_discrete_fields
     dfields[velo].initialize(formula=init_velocity)
     dfields[scalar].initialize(formula=init_scalar)
+    # import sys
+    # sys.exit(1)
     
     # Create a simulation and solve the problem 
     # (do not forget to specify the dt parameter here)

hysop/constants.py

@@ -102,6 +102,8 @@ def implementation_to_backend(implementation):
 
 GhostOperation = EnumFactory.create('GhostOperation',
         ['EXCHANGE', 'ACCUMULATE'])
+ExchangeMethod = EnumFactory.create('ExchangeMethod',
+        ['ISEND_IRECV', 'NEIGHBOR_ALL_TO_ALL_V', 'NEIGHBOR_ALL_TO_ALL_W'])
 
 Precision = EnumFactory.create('Precision',
     ['DEFAULT', 'SAME', 'QUAD', 'LONG_DOUBLE', 'DOUBLE', 'FLOAT', 'HALF'])