>
> 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:
src/eps/impls/davidson/common/dvd_updatev
>
> 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 eig
> 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,
> theta=0x1f396f8, thetai=0x1f39718, maxits=0x
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
/depot/ncn/apps/conte/conte-gcc-petsc35-db
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
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
an
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 ./configure with
> --with-cc
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) PrgEnv-intel/6.0
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 configure error on t
>>>
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 --download-cma
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,
>
> Hong (Mr.)
>
>
>
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
mailto:fdkong...@gmail.com>> wrote:
--with-blaslapack-lib=-mkl -L' + os.environ['MKLROOT'] + '/
--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 to do with using
>
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 use
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 buil
15 matches
Mail list logo
