On Sat, Jan 14, 2012 at 13:32, Junchao Zhang <junchao.zhang at gmail.com>wrote:
> My matrices are from Florida matrix collection. > These ASCII formats cannot be read efficiently in parallel. > I want to use them to test MatMult() > You just want to benchmark sparse matrix kernels? 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. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120114/243e66d6/attachment.htm>
