Hi,<br><br><div class="gmail_quote">On Tue, Mar 27, 2012 at 9:48 PM, Ryan Krauss <span dir="ltr">&lt;<a href="mailto:ryanlists@gmail.com">ryanlists@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I am loosing my mind while trying to debug some code.  I am trying to<br>
find the cause of some differences between numpy analysis and analysis<br>
done first in maxima and then converted to python code.  The maxima<br>
approach is more difficult to do, but seems to lead to the correct<br>
answers.  The core issue seems to be one dot product of a 2x2 and a<br>
2x1 that are both complex numbers:<br>
<br>
here is the 2x2:<br>
<br>
ipdb&gt; submatinv<br>
array([[-0.22740265-1.63857451j, -0.07740957-0.55847886j],<br>
       [-3.20602957-4.93959054j, -0.36746252-1.68352465j]])<br>
<br>
here is the 2x1:<br>
<br>
ipdb&gt; augcol<br>
array([[ -3.74729148e-05-0.0005937j ],<br>
       [  7.96025801e-04+0.01137658j]])<br>
<br>
verifying their shapes and data types:<br>
<br>
ipdb&gt; submatinv.shape<br>
(2, 2)<br>
ipdb&gt; submatinv.dtype<br>
dtype(&#39;complex128&#39;)<br>
ipdb&gt; augcol.shape<br>
(2, 1)<br>
ipdb&gt; augcol.dtype<br>
dtype(&#39;complex128&#39;)<br>
<br>
I need to compute this result:<br>
<br>
ipdb&gt; -1*numpy.dot(submatinv,augcol)<br>
array([[  5.30985737e-05+0.00038316j],<br>
       [  1.72370377e-04+0.00115503j]])<br>
<br>
If I hard code how to do the matrix multiplication, I get the correct<br>
answer (it agrees with Maxima):<br>
<br>
For the first element:<br>
ipdb&gt; -1*(submatinv[0,0]*augcol[0,0]+submatinv[0,1]*augcol[1,0])<br>
(-0.005327660633034575+0.0011288088216130766j)<br>
<br>
For the second<br>
ipdb&gt; -1*(submatinv[1,0]*augcol[0,0]+submatinv[1,1]*augcol[1,0])<br>
(-0.016047752110848554+0.003432076134378004j)<br>
<br>
What is the dot product doing if it isn&#39;t dotting row by column?<br>
<br>
I am not seeing something.<br></blockquote><div>FWIIWO, I can&#39;t either reproduce your results:</div><div><div><font face="&#39;courier new&#39;, monospace">In []: sys.version</font></div><div><font face="&#39;courier new&#39;, monospace">Out[]: &#39;2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]&#39;</font></div>
<div><font face="&#39;courier new&#39;, monospace">In []: np.version.version</font></div><div><font face="&#39;courier new&#39;, monospace">Out[]: &#39;1.6.0&#39;</font></div><div><font face="&#39;courier new&#39;, monospace"><br>
</font></div><div><font face="&#39;courier new&#39;, monospace">In []: s.round(7)</font></div><div><font face="&#39;courier new&#39;, monospace">Out[]: </font></div><div><font face="&#39;courier new&#39;, monospace">array([[-0.2274026-1.6385745j, -0.0774096-0.5584789j],</font></div>
<div><font face="&#39;courier new&#39;, monospace">       [-3.2060296-4.9395905j, -0.3674625-1.6835246j]])</font></div><div><font face="&#39;courier new&#39;, monospace">In []: a.round(7)</font></div><div><font face="&#39;courier new&#39;, monospace">Out[]: </font></div>
<div><font face="&#39;courier new&#39;, monospace">array([[ -3.75000000e-05-0.0005937j],</font></div><div><font face="&#39;courier new&#39;, monospace">       [  7.96000000e-04+0.0113766j]])</font></div><div><font face="&#39;courier new&#39;, monospace"><br>
</font></div><div><span style="font-family:&#39;courier new&#39;,monospace">In []: -1* dot(s, a).round(7)</span><br></div><div><font face="&#39;courier new&#39;, monospace">Out[]: </font></div><div><font face="&#39;courier new&#39;, monospace">array([[-0.0053277+0.0011288j],</font></div>
<div><font face="&#39;courier new&#39;, monospace">       [-0.0160477+0.0034321j]])</font></div></div><div><br></div><div>Regards,</div><div>-eat </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<br>
Thanks,<br>
<br>
Ryan<br>
_______________________________________________<br>
SciPy-User mailing list<br>
<a href="mailto:SciPy-User@scipy.org">SciPy-User@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/scipy-user" target="_blank">http://mail.scipy.org/mailman/listinfo/scipy-user</a><br>
</blockquote></div><br>