<div dir="ltr">Thanks Dave!  I'll provide feedback after testing...<div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 15, 2014 at 9:03 AM, otter <span dir="ltr"><<a href="mailto:otter@otter-rsch.com" target="_blank">otter@otter-rsch.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">After much pain I have produced a C++ version of the LU decomposition which<br>
is suitable for producing adjoint code for ADMB and perhaps cppad. (Don't know what adjoint<br>
code for cppad is as yet!)  This code is about 25 times faster than the current<br>
ADMB code for a 2,000 x 2,000 matrix and about 4 times slower than the Openblas<br>
code which contains optimized assembler and Fortrash.  Any suggestion for improvements<br>
would be welcome.  Per usual I will hold my breath.<br>
<br>
<br>
<br>_______________________________________________<br>
Developers mailing list<br>
<a href="mailto:Developers@admb-project.org">Developers@admb-project.org</a><br>
<a href="http://lists.admb-project.org/mailman/listinfo/developers" target="_blank">http://lists.admb-project.org/mailman/listinfo/developers</a><br>
<br></blockquote></div><br></div>