[petsc-users] pastix solver break at pastix_checkMatrix

[Barry Smith]

  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>