Changeset 5738 for anuga_work/development/anuga_1d/quantity.py
- Timestamp:
- Sep 5, 2008, 3:06:01 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/anuga_1d/quantity.py
r5536 r5738 51 51 #Allocate space for other quantities 52 52 self.centroid_values = zeros(N, Float) 53 self.centroid_backup_values = zeros(N, Float) 53 54 #self.edge_values = zeros((N, 2), Float) 54 55 #edge values are values of the ends of each interval … … 854 855 else: 855 856 qv[k,i] = qc[k] 856 857 858 def backup_centroid_values(self): 859 # Call correct module function 860 # (either from this module or C-extension) 861 #backup_centroid_values(self) 862 N = self.domain.number_of_elements 863 for i in range(self.N): 864 quantity.centroid_backup_values[i] = quantity.centroid_values[i] 865 866 def saxpy_centroid_values(self,a,b): 867 # Call correct module function 868 # (either from this module or C-extension) 869 #saxpy_centroid_values(self,a,b) 870 N = self.domain.number_of_elements 871 for i in range(self.N): 872 quantity.centroid_backup_values[i] = a*quantity.centroid_values[i] + b*quantity.centroid_backup_values[i] 857 873 858 874 class Conserved_quantity(Quantity): … … 869 885 print "Use Quantity instead of Conserved_quantity" 870 886 871 887 872 888 873 889 def newLinePlot(title='Simple Plot'):
Note: See TracChangeset
for help on using the changeset viewer.