I agree. It would be a good idea to have matrices out of numpy as a standalone package.<div><br></div><div>Indeed, having matrices in the numpy core comes at a pedagogical cost. Newcomers (as I once was) do not know which to use. Matrix or array? It turns out that the vast majority of numpy/scipy modules use arrays, so arrays is the preferred way to go.</div>
<div><br></div><div>It would thus be clearer to have arrays in numpy and matrices available as an external package.</div><div><br></div><div>Besides, I think matrices can be pretty tricky when used for teaching. For instance, you have to explain that all the operators work component-wise, except the multiplication! Another caveat is that since matrices are always 2x2, the "scalar product" of two column vectors computed as " x.T * y" will not be a scalar, but a 2x2 matrix. There is also the fact that you must cast all your vectors to column/raw matrices (as in matlab). For all these reasons, I prefer to use arrays and dot for teaching, and I have never had any complaints.</div>
<div><br></div><div>== Olivier</div><div><br><div class="gmail_quote">2009/6/4 Tommy Grav <span dir="ltr"><<a href="mailto:tgrav@mac.com">tgrav@mac.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
On Jun 4, 2009, at 5:41 PM, Alan G Isaac wrote:<br>
> On 6/4/2009 5:27 PM Tommy Grav apparently wrote:<br>
>> Or the core development team split the matrices out of numpy and<br>
>> make it<br>
>> as separate package that the people that use them could pick up and<br>
>> run with.<br>
><br>
><br>
> This too would be a mistake, I believe.<br>
> But it depends on whether a goal is to<br>
> have more people use NumPy. ¬†I believe<br>
> the community will gain from growth.<br>
><br>
</div><div class="im">> In sum, my argument is this:<br>
> Keeping a matrix object in NumPy has substantial<br>
> benefits in encouraging growth of the NumPy<br>
> community, and as far as I can tell, it is<br>
> imposing few costs. ¬†Therefore I think there is<br>
> a very substantial burden on people who propose<br>
> removing the matrix object to demonstrate<br>
> just how the NumPy community will benefit from<br>
> this change.<br>
<br>
</div>This is a perfectly valid argument. I am actually quite happy with the<br>
numpy package as it is (I work in astronomy), I was just pointing out<br>
that if there are few of the core numpy people interested in maintaing<br>
or upgrading the matrix class one solution might be to make it a<br>
scipy-like package that easily can be installed on top of numpy, but<br>
where the code base might be more accessible to those that are<br>
interested in matrices, but feel that numpy is a daunting beast to<br>
tackle.<br>
Some sense of ownership of a matrixpy package might encourage more<br>
people to contribute.<br>
<br>
Just an idea ;-)<br>
<font color="#888888">Tommy<br>
</font><div><div></div><div class="h5"><br>
_______________________________________________<br>
Numpy-discussion mailing list<br>
<a href="mailto:Numpy-discussion@scipy.org">Numpy-discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</div></div></blockquote></div><br></div>