[Scipy-svn] r3848 - trunk/scipy/linsolve

scipy-svn@scip... scipy-svn@scip...
Thu Jan 17 09:41:47 CST 2008


Author: wnbell
Date: 2008-01-17 09:41:44 -0600 (Thu, 17 Jan 2008)
New Revision: 3848

Modified:
   trunk/scipy/linsolve/linsolve.py
Log:
fixed bug with umfpack and 2D RHS


Modified: trunk/scipy/linsolve/linsolve.py
===================================================================
--- trunk/scipy/linsolve/linsolve.py	2008-01-17 15:13:50 UTC (rev 3847)
+++ trunk/scipy/linsolve/linsolve.py	2008-01-17 15:41:44 UTC (rev 3848)
@@ -71,6 +71,8 @@
         if A.dtype.char not in 'dD':
             raise ValueError, "convert matrix data to double, please, using"\
                   " .astype(), or set linsolve.useUmfpack = False"
+        
+        b = asarray(b, dtype=A.dtype).reshape(-1)
 
         family = {'d' : 'di', 'D' : 'zi'}
         umf = umfpack.UmfpackContext( family[A.dtype.char] )



More information about the Scipy-svn mailing list