[Scipy-tickets] [SciPy] #1387: Replace dot(linalg.inv(C), an) by solve(C, an) in scipy.misc.pade

SciPy Trac scipy-tickets@scipy....
Wed Feb 16 15:51:48 CST 2011


#1387: Replace dot(linalg.inv(C),an) by solve(C,an) in scipy.misc.pade
-------------------------+--------------------------------------------------
 Reporter:  robinro      |       Owner:  somebody   
     Type:  enhancement  |      Status:  new        
 Priority:  low          |   Milestone:  Unscheduled
Component:  scipy.misc   |     Version:  devel      
 Keywords:  pade         |  
-------------------------+--------------------------------------------------
 In the computation of the pade coefficients, in scipy/misc/common.py line
 336 is
 {{{
 pq = dot(linalg.inv(C),an)
 }}}
 which should be
 {{{
 solve(C,an)
 }}}
 for numerical stability and speed.
 It shouldn't be to big of an issue in most usecases, but dot(inv(..)) just
 looks wrong.

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


More information about the Scipy-tickets mailing list