<div class="gmail_quote">On Sat, Jan 14, 2012 at 13:32, Junchao Zhang <span dir="ltr"><<a href="mailto:junchao.zhang@gmail.com">junchao.zhang@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":362">My matrices are from Florida matrix collection.</div></blockquote><div><br></div><div>These ASCII formats cannot be read efficiently in parallel.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":362"> I want to use them to test MatMult()</div></blockquote><div><br></div><div>You just want to benchmark sparse matrix kernels?</div><div><br></div><div>Some might consider it to be cheating, but you should consider also reordering the matrices using MatGetOrdering(mat,MATORDERINGRCM,&rperm,&cperm). This will tend to improve cache reuse and can offer a large (e.g. 2x depending on the matrix) speedup relative to the original ordering.</div>
</div>