Changeset 13482 for issm/trunk-jpl/src/m/miscellaneous/MatlabFuncs.py

Timestamp:

09/28/12 16:04:47 (12 years ago)

Author:

jschierm

Message:

NEW: New python functions in support of test514.py (which doesn't yet work).

File:

• issm/trunk-jpl/src/m/miscellaneous/MatlabFuncs.py (modified) (1 diff)

issm/trunk-jpl/src/m/miscellaneous/MatlabFuncs.py

@@ -74 +74 @@
                 raise TypeError("MatlabFunc.det only implemented for shape (2, 2), not for shape %s." % str(a.shape))
 
+def sparse(ivec,jvec,svec,m=0,n=0,nzmax=0):
+        import numpy
+
+        if not m:
+                m=numpy.max(ivec)
+        if not n:
+                n=numpy.max(jvec)
+
+        a=numpy.zeros((m,n))
+
+        for i,j,s in zip(ivec.reshape(-1,order='F'),jvec.reshape(-1,order='F'),svec.reshape(-1,order='F')):
+                a[i-1,j-1]+=s
+
+        return a
+