Re: [petsc-users] Question on incomplete factorization level and fill

Hi All, I just delete the .info file and it works without problem now. Thanks, Danyang On 17-05-24 06:32 PM, Hong wrote: Remove your option '-vecload_block_size 10'. Hong On Wed, May 24, 2017 at 3:06 PM, Danyang Su > wrote: Dear

Re: [petsc-users] Question on incomplete factorization level and fill

Remove your option '-vecload_block_size 10'. Hong On Wed, May 24, 2017 at 3:06 PM, Danyang Su wrote: > Dear Hong, > > I just tested with different number of processors for the same matrix. It > sometimes got "ERROR: Arguments are incompatible" for different number of >

Re: [petsc-users] Question on incomplete factorization level and fill

I don't think this has anything to do with the specific solver but is because you are loading both a vector and matrix from a file and when it uses the default parallel layout for each, because you have -matload_block_size 1 and -vecload_block_size 10 they do not get the same layout.

Re: [petsc-users] Question on incomplete factorization level and fill

Dear Hong, I just tested with different number of processors for the same matrix. It sometimes got "ERROR: Arguments are incompatible" for different number of processors. It works fine using 4, 8, or 24 processors, but failed with "ERROR: Arguments are incompatible" using 16 or 48

Re: [petsc-users] Question on incomplete factorization level and fill

Hi Hong, Awesome. Thanks for testing the case. I will try your options for the code and get back to you later. Regards, Danyang On 17-05-24 12:21 PM, Hong wrote: Danyang : I tested your data. Your matrices encountered zero pivots, e.g. petsc/src/ksp/ksp/examples/tutorials (master) $

Re: [petsc-users] Question on incomplete factorization level and fill

Danyang : I tested your data. Your matrices encountered zero pivots, e.g. petsc/src/ksp/ksp/examples/tutorials (master) $ mpiexec -n 24 ./ex10 -f0 a_react_in_2.bin -rhs b_react_in_2.bin -ksp_monitor -ksp_error_if_not_converged [15]PETSC ERROR: Zero pivot in LU factorization:

Re: [petsc-users] Question on incomplete factorization level and fill

Hi Matt, Yes. The matrix is 45x45 sparse. The hypre takes hundreds of iterates, not for all but in most of the timesteps. The matrix is not well conditioned, with nonzero entries range from 1.0e-29 to 1.0e2. I also made double check if there is anything wrong in the parallel version,

Re: [petsc-users] Question on incomplete factorization level and fill

On Wed, May 24, 2017 at 12:50 PM, Danyang Su wrote: > Hi Matthew and Barry, > > Thanks for the quick response. > > I also tried superlu and mumps, both work but it is about four times > slower than ILU(dt) prec through hypre, with 24 processors I have tested. > You mean the

Re: [petsc-users] Question on incomplete factorization level and fill

Hi Matthew and Barry, Thanks for the quick response. I also tried superlu and mumps, both work but it is about four times slower than ILU(dt) prec through hypre, with 24 processors I have tested. When I look into the convergence information, the method using ILU(dt) still takes 200 to 3000

Re: [petsc-users] Question on incomplete factorization level and fill

> On May 24, 2017, at 2:21 AM, Danyang Su wrote: > > Dear All, > > I use PCFactorSetLevels for ILU and PCFactorSetFill for other preconditioning > in my code to help solve the problems that the default option is hard to > solve. However, I found the latter one,

Re: [petsc-users] Question on incomplete factorization level and fill

On Wed, May 24, 2017 at 2:21 AM, Danyang Su wrote: > Dear All, > > I use PCFactorSetLevels for ILU and PCFactorSetFill for other > preconditioning in my code to help solve the problems that the default > option is hard to solve. However, I found the latter one,