Routine to convert matrix to CSR format

2007-08-25 Thread Ben Tay
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

Best way to find empty row in matrix

2007-08-25 Thread Ben Tay
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

Error during example compilation of test examples

2007-08-25 Thread Ben Tay
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

Routine to convert matrix to CSR format

2007-08-25 Thread Hong Zhang
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,

Best way to find empty row in matrix

2007-08-25 Thread Hong Zhang
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