So for dense matrices, we assume column alignment when we pass in the pointer? That is<br>probably what was wrong. Crap it does say this at the top, but not at the argument. I will<br>augment the documentation.<br><br>  Matt<br>
<br><div class="gmail_quote">On Sun, Jan 17, 2010 at 10:42 AM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">On Fri, 15 Jan 2010 09:58:46 -0600, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
> SeqDense MatMult was inexplicably giving me the wrong answers. I<br>
> checked the BLAS call (and the underlying defines) and they also<br>
> looked right to me. I temporarily checked in the simple-minded loops<br>
> which do give me the correct answer.<br>
<br>
</div>The matrices are column-aligned.  The "simple-minded loops" you checked<br>
in (cc0b86d7fd3f) assume row-aligned matrices, and are thus delivering<br>
the wrong results for everyone that assembled their matrix correctly, so<br>
I'm reverting that change.<br>
<br>
How were you assembling this matrix?<br>
<font color="#888888"><br>
Jed<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>