Ack.&nbsp; I may be able to get references to lu, lu_factor, et al, but they don't actually work with numpy arrays:<br><br>from scipy.linalg import lu,lu_factor,lu_solve<br>import scipy as S<br>A = S.rand(2,2)<br>lu(A)<br>Traceback (most recent call last):
<br>&nbsp; File &quot;&lt;input&gt;&quot;, line 1, in ?<br>&nbsp; File &quot;C:\Python24\Lib\site-packages\scipy\linalg\decomp.py&quot;, line 249, in lu<br>&nbsp;&nbsp;&nbsp; flu, = get_flinalg_funcs(('lu',),(a1,))<br>&nbsp; File &quot;C:\Python24\Lib\site-packages\scipy\linalg\flinalg.py&quot;, line 30, in get_flinalg_funcs
<br>&nbsp;&nbsp;&nbsp; t = arrays[i].dtypechar<br>AttributeError: 'numpy.ndarray' object has no attribute 'dtypechar'<br><br><br>Ok, so, once again, does anyone have an lu_factor / lu_solve implementation in python that I could borrow?<br>
<br>Apologies for the monologue.<br><br>--bb<br><br><br><div><span class="gmail_quote">On 2/20/06, <b class="gmail_sendername">Bill Baxter</b> &lt;<a href="mailto:wbaxter@gmail.com">wbaxter@gmail.com</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Upon further inspection I find that if I call 'from scipy import *' then <a href="http://linalg.lu" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">
linalg.lu</a> etc are defined.<br>But if I do anything else to import scipy like 'import scipy' or 'import scipy as S' or 'from scipy import linalg', then lu, cg etc are not defined.&nbsp; 
<br><br>Why is that?&nbsp; <br><br>I can get at them without importing * by doing 'from scipy.linalg import lu', but that's kind of odd to have to do that.<br><br>--bb<div><span class="e" id="q_109862c934b6c2be_1"><br><br><div>
<span class="gmail_quote">On 2/20/06, <b class="gmail_sendername">
Bill Baxter</b> &lt;<a href="mailto:wbaxter@gmail.com" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">wbaxter@gmail.com</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

This url <a href="http://www.rexx.com/%7Edkuhlman/scipy_course_01.html" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://www.rexx.com/~dkuhlman/scipy_course_01.html</a> seems to
keep turning up in my searches for numpy and scipy things,<br>
but many of the linalg operations it lists don't seem to exist in recent versions of numpy (or scipy).&nbsp; <br><br>Some of them are:<br>
<br>
* norm<br>* factorizations: lu, lu_factor, lu_solve, qr<br>* iterative solvers: cg, cgs, gmres etc.<br><br clear="all">Did these things used to exist in Numeric but they haven't been ported over?&nbsp; Will they be re-introduced sometime?
<br><br>In the short term, the one I'm after right now is LU decompose and solve functionality.&nbsp; Anyone have a numpy implementation?<br><span><br>--Bill Baxter<br>

</span></blockquote></div></span></div><span class="sg"></span></blockquote></div>