Changeset 447

Timestamp:

Oct 26, 2004, 11:55:30 AM (20 years ago)

Author:

ole

Message:

sparse test

File:

• inundation/ga/storm_surge/pyvolution/sparse.py (modified) (4 diffs)

inundation/ga/storm_surge/pyvolution/sparse.py

@@ -81 +81 @@
 
         else:
-            print 'Numeric matrix not yet implemented'
+            raise 'Numeric matrix not yet implemented'
 
         return R
@@ -88 +88 @@
 if __name__ == '__main__':
 
-    from Numeric import allclose
+    from Numeric import allclose, array, Float 
     
     A = Sparse(3,3)
@@ -111 +111 @@
     A[2,2] = 1
 
+    print A
     print A.todense()
 
@@ -121 +122 @@
     assert allclose(u, [6,14,4])
 
+    #Right hand side column
+    v = array([[2,4],[3,4],[4,4]])
+
+    u = A*v[:,0]
+    assert allclose(u, [6,14,4])
+
+    #u = A*v[:,1]
+    #print u