31 Dec
2010
31 Dec
'10
12:12 p.m.
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