1D fields quickfix

hysop/backend/fortran/fftw.py

@@ -36,7 +36,7 @@ class FortranFFTWOperator(ComputationalGraphOperator):
         self.domain   = domain
         self.topology = topology
 
-    def initialize(self,topgraph_method):
+    def initialize(self,topgraph_method=None):
         super(FortranFFTWOperator,self).initialize(topgraph_method)
     
     @debug

hysop/fields/discrete.py

@@ -86,6 +86,8 @@ class DiscreteField(object):
                      for _ in xrange(self.nb_components)])
 
     def set_data(self, new_data):
+        if (self.nb_components==1) and not isinstance(new_data,(list,tuple)):
+            new_data = [new_data]
         check_instance(new_data, (list,tuple))
         new_data = list(new_data)
         backend = self.topology.backend()

hysop/operator/analytic.py

@@ -78,7 +78,6 @@ class Analytic(ComputationalGraphOperator):
         assert simulation is not None, \
             "Missing simulation value for computation."
         for (var,topo) in self.variables.iteritems():
-            print var.name, 'INITIALIZE'
             var.initialize(time=simulation.time, topo=topo)
     
     @classmethod