Changeset 3335 for development/pyvolution-1d/quantity.py

Timestamp:

Jul 14, 2006, 2:11:20 PM (18 years ago)

Author:

jakeman

Message:

New file dam.py tests numerical solution against analytical solution of
Stoker 57. Numerical results match Stoker's roughly but oscillations
present. Most likely limiting is not working correctly. Q.limit does
have an effect but unsure whether this effect is the right one.

File:

• development/pyvolution-1d/quantity.py (modified) (5 diffs)

development/pyvolution-1d/quantity.py

@@ -134 +134 @@
                   Default is "vertices"
         """
-
-        print "In set_function value"
-        print self.domain.vertices
         
         if location == 'centroids':
-
-            print "In set at centroid"            
+            
             P = self.domain.centroids
             self.set_values(f(P), location)
@@ -146 +142 @@
             #Vertices
             P = self.domain.get_vertices()
-
-            print "P"
-            print P
-            print P[:,0]
-            print f(P[:,0])
-
-            print P[:,1]
-            print f(P[:,1])
+            
             for i in range(2):
-                self.vertex_values[:,i] = f(P[:,i])
-
-
-            
-
-            print "In set at vertices"
-            print self.vertex_values
-
-        print "Out set_function value"
-        print self.domain.vertices        
-
+                self.vertex_values[:,i] = f(P[:,i])        
 
     def set_array_values(self, values, location='vertices'):
@@ -564 +543 @@
         vertex values are updated
         """
-
+        print "in Q.limit"
         from Numeric import zeros, Float
 
@@ -706 +685 @@
     print Q2.centroid_values
     print Qc
-
+    raw_input('press_return')
+    
     g3 = newLinePlot('plot 3')
     linePlot(g3,Xc,Qc)
@@ -720 +700 @@
     raw_input('press return')
 
-
-
-
-
-
-