Sorry. Yes there was another bug. What caused both of these problems is that Pastix requires a symmetric nonzero structure and the interface made some assumptions that the PETSc matrix had a symmetric nonzero structure which is usually true, hence it did not crash for most matrices people use.
I've attached another copy of pastix.c follow the same procedure again. Barry On Dec 30, 2010, at 10:47 PM, Gong Ding wrote: > Dear Barry > Sorry, I should complement that > when I set DIFFERENT_NONZERO_PATTERN to the MatStructure, pastix works well. > But it crash for SAME_NONZERO_PATTERN. > However, I always use SAME_NONZERO_PATTERN before, which works for ksp > solvers, mumps, superlu_dist, etc. > > Yours > Gong Ding > > -------------- next part -------------- A non-text attachment was scrubbed... Name: pastix.c Type: application/octet-stream Size: 28044 bytes Desc: not available URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101231/54d5ca39/attachment.obj>
