Re: [petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

That version is 4-year old, I recommend you to upgrade to the latest version. We have made the change suggested by Barry. It will be included in 3.9.2. Jose > El 4 may 2018, a las 5:45, Harshad Sahasrabudhe > escribió: > > I think the residual of the first

Re: [petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

> > I think the residual of the first unconverged eigenvalue is set to INF on > purpose and not calculated in the iteration when a converged eigenvalue is > found. Dividing the INF residual is a mistake: Found the line where d->nR[i] is set to infinity:

Re: [petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

> > Jumps like this usually indicate something is going very poorly in the > algorithm. I hope the algorithmic experts have a chance to look at your > case. I think the residual of the first unconverged eigenvalue is set to INF on purpose and not calculated in the iteration when a converged

Re: [petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

> On May 3, 2018, at 10:24 PM, Harshad Sahasrabudhe > wrote: > > Hi Barry, > > There's an overflow in the division: > > Program received signal SIGFPE, Arithmetic exception. > 0x2aaab377ea26 in dvd_improvex_jd_lit_const_0 (d=0x1d29078, i=0, >

Re: [petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

Hi Barry, There's an overflow in the division: Program received signal SIGFPE, Arithmetic exception. 0x2aaab377ea26 in dvd_improvex_jd_lit_const_0 (d=0x1d29078, i=0, theta=0x1f396f8, thetai=0x1f39718, maxits=0x7fff816c, tol=0x7fff8140) at

Re: [petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

at /depot/ncn/apps/conte/conte-gcc-petsc35-dbg/libs/slepc/build-real/src/eps/impls/davidson/common/dvd_improvex.c:1112 1112 if (d->nR[i]/a < data->fix) { Likely the problem is due to the division by a when a is zero. Perhaps the code needs above a check that a is not zero. Or rewrite the

[petsc-users] [SLEPc] SIGFPE Arithmetic exception in EPSGD

Hello, I am solving for the lowest eigenvalues and eigenvectors of symmetric positive definite matrices in the generalized eigenvalue problem. I am using the GD solver with the default settings of PCBJACOBI. When I run a standalone executable on 16 processes which loads the matrices from a file