[petsc-users] Sparse matrix partitioning in PETSc

[Jed Brown]

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>