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

Re: [petsc-users] configure error

Thanks, I can not find the damn (PGI) c++ compiler (on SUMMIT) but this is for a Fortran code so I just set it to =0 and it seems to be chugging along. On Thu, May 3, 2018 at 3:14 PM, Smith, Barry F. wrote: > > Mark, > > You pass the C and C++ compiler names to

Re: [petsc-users] Could not determine how to create a shared library!

The MKL version you are trying to use does not match the environment setting on Theta. hongzh@thetalogin6:~/Projects/petsc (master)$ module list Currently Loaded Modulefiles: 1) modules/3.2.10.6 10) dmapp/7.1.1-6.0.4.0_46.2__gb8abda2.ari19)

Re: [petsc-users] configure error

Mark, You pass the C and C++ compiler names to ./configure with --with-cc=mpicc --with-fc=mpif90 but do not pass a C++ compiler hence it defaults to g++ which does not know about MPI. Barry > On May 3, 2018, at 2:02 PM, Mark Adams wrote: > > I am getting this

Re: [petsc-users] configure error

>>> Configure Options: --configModules=PETSc.Configure --optionsModule=config.compilerOptions --COPTFLAGS=-O2 --CXXOPTFLAGS=-O2 --FOPTFLAGS=-O2 --with-ssl=0 --with-batch=0 --prefix=/ccs/proj/env003/petscv3.9-opt64-summit-pgi --download-metis --with-hwloc=0 --download-parmetis

Re: [petsc-users] Could not determine how to create a shared library!

On Thu, May 3, 2018 at 11:50 AM, Zhang, Hong wrote: > Alternatively you can use --with-blaslapack-dir=/opt/ > intel/compilers_and_libraries/linux/mkl/lib/intel64 to let petsc pick the > right libs for you. > This used to work, but does not work any more. Thanks, Fande, >

Re: [petsc-users] Could not determine how to create a shared library!

Alternatively you can use --with-blaslapack-dir=/opt/intel/compilers_and_libraries/linux/mkl/lib/intel64 to let petsc pick the right libs for you. Hong (Mr.) On May 3, 2018, at 11:32 AM, Fande Kong > wrote: --with-blaslapack-lib=-mkl -L' +

Re: [petsc-users] Could not determine how to create a shared library!

--with-blaslapack-lib=-mkl -L' + os.environ['MKLROOT'] + '/lib/intel64 works. Fande, On Thu, May 3, 2018 at 10:09 AM, Satish Balay wrote: > Ok you are not 'building blaslapack' - but using mkl [as per > configure.log]. > > I'll have to check the issue. It might be something

Re: [petsc-users] Could not determine how to create a shared library!

Ok you are not 'building blaslapack' - but using mkl [as per configure.log]. I'll have to check the issue. It might be something to do with using mkl as a static library.. Hong might have some suggestions wrt theta builds. Satish On Thu, 3 May 2018, Satish Balay wrote: > Perhaps you should

Re: [petsc-users] Could not determine how to create a shared library!

Perhaps you should use MKL on theta? Again check config/examples/arch-cray-xc40-knl-opt.py Satish On Thu, 3 May 2018, Kong, Fande wrote: > Thanks, > > I get the PETSc complied, but theta does not like the shared lib, I think. > > I am switching back to a static lib. I ever successfully