[Scipy-svn] r3848 - trunk/scipy/linsolve
scipy-svn at scipy.org
scipy-svn at scipy.org
Thu Jan 17 10:41:47 EST 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