Changeset 599
- Timestamp:
- Nov 19, 2004, 12:13:15 PM (20 years ago)
- Location:
- inundation/ga/storm_surge/pyvolution
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/pyvolution/cg_solve.py
r477 r599 58 58 i = i+1 59 59 if i%iprint == 0 : 60 print 'i = %g rTr = %20.15e'% (i,rTr) 60 pass 61 #FIXME: Should depend on verbosity 62 #print 'i = %g rTr = %20.15e'% (i,rTr) 61 63 62 64 #FIXME: Should this raise an exception? -
inundation/ga/storm_surge/pyvolution/sparse_ext.c
r597 r599 52 52 *y; //Return vector array 53 53 54 double *xdata; 55 54 56 55 int dimensions[1], M, err; 57 int xdim, i;58 56 59 57 // Convert Python arguments to C -
inundation/ga/storm_surge/pyvolution/test_cg_solve.py
r596 r599 103 103 104 104 # Convert to csr format 105 print 'start covert'105 #print 'start covert' 106 106 A = Sparse_CSR(A) 107 print 'finish covert'107 #print 'finish covert' 108 108 b = A*xe 109 109 x = conjugate_gradient(A,b,b,iprint=20) -
inundation/ga/storm_surge/pyvolution/test_sparse.py
r596 r599 181 181 A[2,0] = 5 182 182 183 print ' '184 print A.todense()183 #print ' ' 184 #print A.todense() 185 185 186 186 B = Sparse_CSR(A) 187 187 188 print B.todense()188 #print B.todense() 189 189 190 190 C = [1, 2, 3] … … 194 194 C2 = [[1,2],[2,4],[3,6]] 195 195 196 print B*C2196 #print B*C2 197 197 198 198 assert allclose(B*C2, [[15.0, 30.0],[10.0, 20.0],[8.0, 16.0],[0.0, 0.0]])
Note: See TracChangeset
for help on using the changeset viewer.