Re: [petsc-users] Matrix is missing diagonal entry

2016-10-18 Thread Barry Smith
You could also fork off of master and find all the places where the code depends on diagonal entries existing and modify the code in each of these places to handle the missing diagonal. Then run all the test suite and add some new tests that explicitly test this functionality and then make

Re: [petsc-users] Matrix is missing diagonal entry

2016-10-18 Thread Kong, Fande
Thanks, Hong and Jed. I am going to explicitly add a few zeros into the matrix. Regards, Fande, On Tue, Oct 18, 2016 at 9:46 AM, Hong wrote: > You need set 0.0 to the diagonals. > Diagonal storage is used in PETSc library. > > Hong > > > On Tue, Oct 18, 2016 at 10:11 AM,

Re: [petsc-users] Matrix is missing diagonal entry

2016-10-18 Thread Hong
You need set 0.0 to the diagonals. Diagonal storage is used in PETSc library. Hong On Tue, Oct 18, 2016 at 10:11 AM, Kong, Fande wrote: > Hi Developers, > > Any reason to force users provide a matrix which does not miss any > diagonal entries when using a LU-type solver? >

Re: [petsc-users] Matrix is missing diagonal entry

2016-10-18 Thread Jed Brown
"Kong, Fande" writes: > Hi Developers, > > Any reason to force users provide a matrix which does not miss any diagonal > entries when using a LU-type solver? Automatically adding the entries for incomplete factorization is a slight increase in complexity and means that