blah .. for spline map

2026-07-09 20:30:17 +08:00 · 2015-09-25 23:58:08 -07:00
parent 23a3768b52
commit 2b17ad45ac
1 changed files with 8 additions and 1 deletions
@@ -991,4 +991,11 @@ class SplineMap(IdentityMap):
                    fderiv = (flinesa-flinesb)/(2*dy)                                
                    g3[:,i] = Utils.sdiag(alpha*(sig2-sig1)/(1.+(alpha*f)**2)/np.pi)*fderiv
        else :
-            raise(Exception("Not Implemented for Y and Z, your turn :)"))
+            raise(Exception("Not Implemented for Y and Z, your turn :)"))
+        return sp.csr_matrix(np.c_[g1,g2,g3])  
+
+
+
+
+
+