Changeset 606
 Timestamp:
 Nov 22, 2004, 12:18:16 AM (20 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

inundation/ga/storm_surge/pyvolution/least_squares.py
r595 r606 30 30 from Numeric import zeros, array, Float, Int, dot, transpose 31 31 from LinearAlgebra import solve_linear_equations 32 from sparse import Sparse 32 from sparse import Sparse, Sparse_CSR 33 33 from cg_solve import conjugate_gradient, VectorShapeError 34 34 … … 203 203 self.B = self.AtA 204 204 205 205 #Convert self.B matrix to CSR format for faster matrix vector 206 self.B = Sparse_CSR(self.B) 206 207 207 208 def build_interpolation_matrix_A(self, point_coordinates): … … 515 516 516 517 518 517 519 return conjugate_gradient(self.B, Atz, Atz,imax=2*len(Atz) ) 518 520 #FIXME: Should we store the result here for later use? (ON)
Note: See TracChangeset
for help on using the changeset viewer.