Abinit +intel compilators wierd issue while opt geom [SOLVED]
Posted: Wed Jul 09, 2014 4:35 pm
I have a vary strange situation. On the same machine (2 xenon e5-2340 processor based) I've compiled abinit 7.8.1 with both intel and gcc compilators. The configuration is as follows:
For intel (intel mpi used):
For gcc (mpich used):
And while optimizing the structure with abinit compiled by intel compilators it fails at some step, while optimizing with gcc abinit it finishes smooth. Here is the input I have:
ABINIT fails with warning
Moreover if I use abinit compiled with intel compilator on another machine (cluster) it whent smooth also. The optimization process goes absolutle the same way up to 4th step 5th step differs and the 6th fails.
There is somthing different with reduced forces they all goes to zero in basis at the fifth (pre fails) step of optimisation.
While in previous step its of order -4.30318458033374E-33
Maybe there is somekind of precision issure? Is there paramiter I should pass to intel compilator to rise precision of calculations.
In cluster (Xeon Processor X5550) where abinit compiled by intel the compilator set is 2013 upd 1. I was trying the set in my machine (Xenon E5-2340) to compile good abinit binary with no succes, but I can use the binary compiled in cluster on my E5-2340 based computer with full success. I've tried to use mpich and openmpi with no success also. Maybe I should add some additional flags to intel compilator?
For intel (intel mpi used):
Code: Select all
export FC=$MPICH/bin/mpiifort
export CC=$MPICH/bin/mpiicc
export CXX=$MPICH/bin/mpiicpc
./configure --prefix=/opt/abinit-$VER --enable-mpi --enable-mpi-io \
--with-linalg-incs="-I$INTEL_COMP/mkl/include" \
--with-fft-flavor="fftw3-mkl" --with-fft-incs="-I$INTEL_COMP/mkl/include/fftw" \
--with-fft-libs="-L$INTEL_COMP/mkl/lib/intel64 -lfftw3xf_intel" \
--with-linalg-flavor="mkl" \
--with-linalg-libs="-L$INTEL_COMP/mkl/lib/intel64 -Wl,--start-group -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -Wl,--end-group -liomp5 -lpthread -limf -lsvml -lirc" \
--with-linalg-incs="-I$INTEL_COMP/mkl/include" \
--with-dft-flavor="atompaw+bigdft+libxc+wannier90"
For gcc (mpich used):
Code: Select all
./configure --prefix=/opt/abinit-$VER --with-mpich=$MPICH2 \
--with-dft-flavor="atompaw+bigdft+libxc+wannier90"
And while optimizing the structure with abinit compiled by intel compilators it fails at some step, while optimizing with gcc abinit it finishes smooth. Here is the input I have:
Code: Select all
acell 5.88 5.88 57.1
angdeg 90.000 90.000 120.000
rprim sqrt(0.75) 0.5 0.0
-sqrt(0.75) 0.5 0.0
0.0 0.0 1.0
ntypat 3
znucl 31 7 13
natom 24
typat 1 2 1 2 1 2 1 2 1 2 3 2 3 2 3 2 3 2 3 2 3 2 3 2
xred
1/3 2/3 0.0000000000
1/3 2/3 0.0637500000
2/3 1/3 0.0841584158
2/3 1/3 0.1479084158
1/3 2/3 0.1683168317
1/3 2/3 0.2320668317
2/3 1/3 0.2524752475
2/3 1/3 0.3162252475
1/3 2/3 0.3366336634
1/3 2/3 0.4003836634
2/3 1/3 0.4207920792
2/3 1/3 0.4845420792
1/3 2/3 0.5049504950
1/3 2/3 0.5687004950
2/3 1/3 0.5923283505
2/3 1/3 0.6536028603
1/3 2/3 0.6732187925
1/3 2/3 0.7344933023
2/3 1/3 0.7541092345
2/3 1/3 0.8153837443
1/3 2/3 0.8349996764
1/3 2/3 0.8962741862
2/3 1/3 0.9158901184
2/3 1/3 0.9771646282
ixc 7 # Teter XC functional
ecut 25 # Kinetic energy cut-off (in Ha)
ecutsm 0.5
dilatmx 1.25
kptopt 1
ngkpt 6 6 4
nshiftk 1
shiftk 0.0 0.0 0.5
iscf 5 # Pulay mixing (default for norm-conserving PSP)
nstep 150 # Modify default
ionmov 3
ntime 150
toldff 1.0d-4
tolmxf 1.0d-3
optforces 1
optcell 0
use_gpu_cuda 0
AUTOPARAL=0
ABINIT fails with warning
Potential-based CG line minimization not converged after13 restarts.
Action : read the eventual warnings about lack of convergence.
Some might be relevant. Otherwise, raise nband. Returning
src_file: scfcge.F90
src_line: 375
Moreover if I use abinit compiled with intel compilator on another machine (cluster) it whent smooth also. The optimization process goes absolutle the same way up to 4th step 5th step differs and the 6th fails.
There is somthing different with reduced forces they all goes to zero in basis at the fifth (pre fails) step of optimisation.
Code: Select all
Cartesian forces (fcart) [Ha/bohr]; max,rms= 9.03107E-03 2.51165E-03 (free atoms)
0.00000000000000E+00 0.00000000000000E+00 1.54667733343065E-03
0.00000000000000E+00 0.00000000000000E+00 -2.06600109156256E-03
0.00000000000000E+00 0.00000000000000E+00 -9.83194404385286E-04
0.00000000000000E+00 0.00000000000000E+00 -4.16727320416967E-03
0.00000000000000E+00 0.00000000000000E+00 1.68339746419727E-03
0.00000000000000E+00 0.00000000000000E+00 -1.72370601891689E-03
0.00000000000000E+00 0.00000000000000E+00 5.67865482140431E-03
0.00000000000000E+00 0.00000000000000E+00 4.59214632078102E-03
0.00000000000000E+00 0.00000000000000E+00 4.48469221935451E-03
0.00000000000000E+00 0.00000000000000E+00 -1.97601745013876E-03
0.00000000000000E+00 0.00000000000000E+00 8.04770880797129E-03
0.00000000000000E+00 0.00000000000000E+00 6.86837849770853E-03
0.00000000000000E+00 0.00000000000000E+00 -1.29830590065030E-03
0.00000000000000E+00 0.00000000000000E+00 -9.03107255582250E-03
0.00000000000000E+00 0.00000000000000E+00 8.13504889824362E-03
0.00000000000000E+00 0.00000000000000E+00 -2.12193114733188E-03
0.00000000000000E+00 0.00000000000000E+00 -5.56983517030434E-03
0.00000000000000E+00 0.00000000000000E+00 -4.53171312473877E-03
0.00000000000000E+00 0.00000000000000E+00 5.62000940042989E-04
0.00000000000000E+00 0.00000000000000E+00 -1.60093617041179E-03
0.00000000000000E+00 0.00000000000000E+00 9.81585503490430E-04
0.00000000000000E+00 0.00000000000000E+00 -2.38942398773797E-03
0.00000000000000E+00 0.00000000000000E+00 -8.33702492019031E-04
0.00000000000000E+00 0.00000000000000E+00 -4.28717808843490E-03
While in previous step its of order -4.30318458033374E-33
Maybe there is somekind of precision issure? Is there paramiter I should pass to intel compilator to rise precision of calculations.
In cluster (Xeon Processor X5550) where abinit compiled by intel the compilator set is 2013 upd 1. I was trying the set in my machine (Xenon E5-2340) to compile good abinit binary with no succes, but I can use the binary compiled in cluster on my E5-2340 based computer with full success. I've tried to use mpich and openmpi with no success also. Maybe I should add some additional flags to intel compilator?