Changeset 7868 for anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py

Timestamp:

Jun 22, 2010, 5:30:32 PM (13 years ago)

Author:

steve

Message:

Added in some more c based limiters

File:

• anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py (modified) (3 diffs)

anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py

@@ -17 +17 @@
                  coordinates,
                  boundary = None,
-                 conserved_quantities = None,
-                 evolved_quantities = None,
-                 other_quantities = None,
+                 conserved_quantities = [],
+                 evolved_quantities = [],
+                 other_quantities = [],
                  tagged_elements = None):
         """
@@ -35 +35 @@
         self.number_of_elements = N = len(self.coordinates)-1
 
-        self.beta = 1.0
-        self.set_limiter("minmod_kurganov")
         self.set_CFL(CFL)
         self.set_timestepping_method(timestepping_method)
@@ -869 +867 @@
     
 
+    def get_beta(self):
+
+        warn('limiter parameter beta associated with quantity not domain')
+
+    def set_beta(self,beta):
+        """Set the same limiter beta parameter to all evolving quantities
+        """
+
+        for name in self.evolved_quantities:
+            Q = self.quantities[name]
+            Q.set_beta(beta)
+
+
     def get_limiter(self):
-        return self.limiter
+
+        warn('limiter associated with quantity not domain')
 
     def set_limiter(self,limiter):
 
-        possible_limiters = \
-        ['pyvolution', 'minmod_steve', 'minmod', 'minmod_kurganov', 'superbee', 'vanleer', 'vanalbada']
-
-        if limiter in possible_limiters:
-            self.limiter = limiter
-            return
-
-        msg = '%s is an incorrect limiter type.\n'% limiter
-        msg += 'Possible types are: '+ ", ".join(["%s" % el for el in possible_limiters])
-        raise Exception, msg
+        for name in self.evolved_quantities:
+            Q = self.quantities[name]
+            Q.set_limiter(limiter)
+