Should have mentioned -- I was using numpy 0.9.4 / scipy 0.4.4.<br>Looks like it works in numpy 0.9.5 / scipy 0.4.6<br><br>But matplotlib, which I also need, hasn't been updated for numpy 0.9.5 yet. :-(<br><br>It's also still pretty weird to me that you have to do &quot;from 
scipy.linalg import lu&quot; specifically.&nbsp; And then after doing that one import, then all the other scipy.linalg.* functions magically spring into existence too.&nbsp; Is that sort of hing expected behavior from Python imports?
<br><br>&gt;&gt;&gt; import numpy as N<br>&gt;&gt;&gt; import scipy as S<br>&gt;&gt;&gt; <a href="http://S.linalg.lu">S.linalg.lu</a><br>Traceback (most recent call last):<br>&nbsp; File &quot;&lt;input&gt;&quot;, line 1, in ?
<br>AttributeError: 'module' object has no attribute 'lu'<br>&gt;&gt;&gt; from scipy.linalg import lu<br>&gt;&gt;&gt; S.linalg.lu(N.rand(2,2))<br>(array([[ 0.,&nbsp; 1.],<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [ 1.,&nbsp; 0.]]), array([[ 1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ,&nbsp; 0.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ],
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [ 0.18553085,&nbsp; 1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ]]), array([[ 0.71732168,&nbsp; 0.48540043],<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [ 0.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ,&nbsp; 0.61379118]]))<br>&gt;&gt;&gt; (N.__version__, S.__version__)<br>('0.9.5', '0.4.6')<br><br>--bb<br><br><div><span class="gmail_quote">
On 2/20/06, <b class="gmail_sendername">Nils Wagner</b> &lt;<a href="mailto:nwagner@mecha.uni-stuttgart.de">nwagner@mecha.uni-stuttgart.de</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;">
Bill Baxter wrote:<br>&gt; Ack.&nbsp;&nbsp;I may be able to get references to lu, lu_factor, et al, but<br>&gt; they don't actually work with numpy arrays:<br>&gt;<br>&gt; from scipy.linalg import lu,lu_factor,lu_solve<br>&gt; import scipy as S
<br>&gt; A = S.rand(2,2)<br>&gt; lu(A)<br>&gt; Traceback (most recent call last):<br>&gt;&nbsp;&nbsp; File &quot;&lt;input&gt;&quot;, line 1, in ?<br>&gt;&nbsp;&nbsp; File &quot;C:\Python24\Lib\site-packages\scipy\linalg\decomp.py&quot;, line
<br>&gt; 249, in lu<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; flu, = get_flinalg_funcs(('lu',),(a1,))<br>&gt;&nbsp;&nbsp; File &quot;C:\Python24\Lib\site-packages\scipy\linalg\flinalg.py&quot;, line<br>&gt; 30, in get_flinalg_funcs<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; t = arrays[i].dtypechar
<br>&gt; AttributeError: 'numpy.ndarray' object has no attribute 'dtypechar'<br>&gt;<br>&gt;<br>&gt; Ok, so, once again, does anyone have an lu_factor / lu_solve<br>&gt; implementation in python that I could borrow?<br>&gt;
<br>&gt; Apologies for the monologue.<br>&gt;<br>&gt; --bb<br>&gt;<br>&gt;<br>&gt; On 2/20/06, *Bill Baxter* &lt;<a href="mailto:wbaxter@gmail.com">wbaxter@gmail.com</a><br>&gt; &lt;mailto:<a href="mailto:wbaxter@gmail.com">
wbaxter@gmail.com</a>&gt;&gt; wrote:<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; Upon further inspection I find that if I call 'from scipy import<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; *' then <a href="http://linalg.lu">linalg.lu</a> &lt;<a href="http://linalg.lu">http://linalg.lu
</a>&gt; etc are defined.<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; But if I do anything else to import scipy like 'import scipy' or<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 'import scipy as S' or 'from scipy import linalg', then lu, cg etc<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; are not defined.<br>&gt;<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp; Why is that?<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; I can get at them without importing * by doing 'from scipy.linalg<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; import lu', but that's kind of odd to have to do that.<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; --bb<br>&gt;<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; On 2/20/06, * Bill Baxter* &lt;
<a href="mailto:wbaxter@gmail.com">wbaxter@gmail.com</a><br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; &lt;mailto:<a href="mailto:wbaxter@gmail.com">wbaxter@gmail.com</a>&gt;&gt; wrote:<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; This url <a href="http://www.rexx.com/~dkuhlman/scipy_course_01.html">
http://www.rexx.com/~dkuhlman/scipy_course_01.html</a><br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;<a href="http://www.rexx.com/%7Edkuhlman/scipy_course_01.html">http://www.rexx.com/%7Edkuhlman/scipy_course_01.html</a>&gt; seems<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; to keep turning up in my searches for numpy and scipy things,
<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; but many of the linalg operations it lists don't seem to exist<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; in recent versions of numpy (or scipy).<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Some of them are:<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; * norm<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; * factorizations: lu, lu_factor, lu_solve, qr
<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; * iterative solvers: cg, cgs, gmres etc.<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Did these things used to exist in Numeric but they haven't<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; been ported over?&nbsp;&nbsp;Will they be re-introduced sometime?<br>&gt;<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; In the short term, the one I'm after right now is LU decompose<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; and solve functionality.&nbsp;&nbsp;Anyone have a numpy implementation?<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; --Bill Baxter<br>&gt;<br>No problem here.<br>
<br>&gt;&gt;&gt; from scipy.linalg import lu,lu_factor,lu_solve<br>&gt;&gt;&gt; import scipy as S<br>&gt;&gt;&gt; A = S.rand(2,2)<br>&gt;&gt;&gt; lu(A)<br>(array([[ 0.,&nbsp;&nbsp;1.],<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [ 1.,&nbsp;&nbsp;0.]]), array([[ 1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;,&nbsp;&nbsp;0.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;],
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [ 0.81367315,&nbsp;&nbsp;1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;]]), array([[ 0.49886054,&nbsp;&nbsp;0.57065709],<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [ 0.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;, -0.30862809]]))<br>&gt;&gt;&gt; S.__version__<br>'0.4.7.1614'<br><br><br>Nils<br><br></blockquote></div><br><br>