[Scipy-tickets] [SciPy] #1770: root() doesn't work well with matrices

SciPy Trac scipy-tickets@scipy....
Thu Nov 8 12:29:48 CST 2012


#1770: root() doesn't work well with matrices
----------------------------+-----------------------------------------------
 Reporter:  pv              |       Owner:  somebody   
     Type:  defect          |      Status:  new        
 Priority:  normal          |   Milestone:  Unscheduled
Component:  scipy.optimize  |     Version:  0.11.0     
 Keywords:                  |  
----------------------------+-----------------------------------------------
 {{{
 from scipy.optimize import root
 import numpy as np

 x0 = np.matrix([0, 0])
 print root(lambda z: z - 1, x0)
 }}}
 The problem is that minpack.py:179 does
 {{{
     x0 = array(x0, ndmin=1)
     n = len(x0)
     if type(args) != type(()):
         args = (args,)
     _check_func('fsolve', 'func', func, x0, args, n, (n,))
 }}}
 which will not work for 2-D arrays.

 Obviously, tests are missing... There are possibly similar problems with
 the optimization routines, should be also checked.

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1770>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list