Hi,
May I know if there's a routine which can convert a matrix initially
created using MatCreateSeqAIJ to compressed row format (CSR) ? Then I
would like to view it in my visual fortran debugger. In other words, the
matrix is in standard array format, not Mat.
Thanks
Hi,
when KSPSolve is called, I got the following errors:
Detected zero pivot in LU factorization
Empty row in matrix!
There's definitely some errors in the code, but I just can't detect it.
PETSc does not tell me which row. How can I find out? Btw, I'm working
in visual fortran
Thank you
Hi,
I didn't touch the petscmat.h and I've verified that the 2 lines are
there. Could it be some compilation error of PETSc? I can send in the
configure.log if required.
thanks.
Satish Balay wrote:
I can't reproduce this error. Can you verify that finclude/petscmat.h
is unmodified?
It
You can use
MatGetRowIJ()
(http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatGetRowIJ.html)
and
MatGetArray()
(http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatGetArray.html#MatGetArray)
Hong
On Sat, 25 Aug 2007,
Your matrix might have an empty row.
If your matrix indeed has an empty row, you need set val=0
to the diagonal entry.
When petsc LU or ILU detects a zero, it
should give an error of type
[0]PETSC ERROR: Detected zero pivot in LU factorization
see