Changeset 405 for inundation/ga/storm_surge/pyvolution/least_squares.py

Timestamp:

Oct 15, 2004, 9:12:12 AM (20 years ago)

Author:

duncan

Message:

if alpha is 0, don't build the D matrix

File:

• inundation/ga/storm_surge/pyvolution/least_squares.py (modified) (2 diffs)

inundation/ga/storm_surge/pyvolution/least_squares.py

@@ -64 +64 @@
     
     # load in the .xya file
-    print point_file,mesh_file,alpha
-
     
     point_dict = load_xya_file(point_file)
@@ -178 +176 @@
         """Build final coefficient matrix
         """
-
-        self.build_smoothing_matrix_D()
+        if self.alpha <> 0:
+            self.build_smoothing_matrix_D()
+        
         
         if point_coordinates:
             self.build_interpolation_matrix_A(point_coordinates)
             AtA = dot(self.At, self.A)
-            self.B = AtA + self.alpha*self.D
+            
+            if self.alpha <> 0:
+                self.B = AtA + self.alpha*self.D
+            else:
+                self.B = AtA