[Scipy-tickets] [SciPy] #1306: Allow for custom inner products in CG, ...

SciPy Trac scipy-tickets@scipy....
Thu Oct 14 10:50:18 CDT 2010


#1306: Allow for custom inner products in CG,...
---------------------------------+------------------------------------------
 Reporter:  nicki                |       Owner:  wnbell
     Type:  enhancement          |      Status:  new   
 Priority:  normal               |   Milestone:  0.9.0 
Component:  scipy.sparse.linalg  |     Version:  0.8.0 
 Keywords:                       |  
---------------------------------+------------------------------------------

Comment(by nicki):

 Somewhere in CG, once per step, you'll find

     Ax = A * x

 or similar. Indeed, you can use any linear operator here, so there's not a
 problem.

 At some _other_ points, inner products are calculated, e.g.,

     alpha = < z, p >

 And here's where the problem sits. Typicall you'll fine dot() or vdot()
 here, but one may want to use something else, e.g., vdot().real.

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


More information about the Scipy-tickets mailing list