MPI error in cluster SCF eigenvaluesolver (OpenMX3.7) |
- Date: 2014/05/26 23:24
- Name: Jae Nyeong Kim
<algetnya@postech.ac.kr>
- Dear Professor,
I have compiled OpenMx with the following option with out error.
MKLROOT=/opt/intel/Compiler/11.1/059 CC = mpicc -O3 -xHOST -openmp -ip -no-prec-div -mkl=parallel -I$(MKLROOT)/mkl/include/fftw -shared-intel FC = mpif90 -O3 -xHOST -openmp -ip -no-prec-div -mkl=parallel -I$(MKLROOT)/mkl/include/fftw -shared-intel LIB = -L$(MKLROOT)/lib/intel64 -lifcore -limf
However, in the test calculation. I have error message like
./openmx Methane.dat
The number of threads in each node for OpenMP parallelization is 1.
******************************************************* ******************************************************* Welcome to OpenMX Ver. 3.7.8 Copyright (C), 2002-2013, T. Ozaki OpenMX comes with ABSOLUTELY NO WARRANTY. This is free software, and you are welcome to redistribute it under the constitution of the GNU-GPL. ******************************************************* *******************************************************
<Input_std> Your input file was normally read. <Input_std> The system includes 2 species and 5 atoms.
******************************************************* PAO and VPS *******************************************************
<SetPara_DFT> PAOs of species H were normally found. <SetPara_DFT> PAOs of species C were normally found. <SetPara_DFT> VPSs of species H were normally found. H_PBE13.vps is j-dependent. In case of scf.SpinOrbit.Coupling=off, j-dependent pseudo potentials are averaged by j-degeneracy, which corresponds to a scalar relativistic treatment. <SetPara_DFT> VPSs of species C were normally found. C_PBE13.vps is j-dependent. In case of scf.SpinOrbit.Coupling=off, j-dependent pseudo potentials are averaged by j-degeneracy, which corresponds to a scalar relativistic treatment.
******************************************************* Fourier transform of PAO and projectors of VNL *******************************************************
<FT_PAO> Fourier transform of pseudo atomic orbitals <FT_NLP> Fourier transform of non-local projectors <FT_ProExpn_VNA> Fourier transform of VNA separable projectors <FT_VNA> Fourier transform of VNA potentials <FT_ProductPAO> Fourier transform of product of PAOs
******************************************************* Allocation of atoms to proccesors at MD_iter= 1 *******************************************************
proc = 0 # of atoms= 5 estimated weight= 5.00000
******************************************************* Analysis of neigbors and setting of grids *******************************************************
TFNAN= 20 Average FNAN= 4.00000 TSNAN= 0 Average SNAN= 0.00000 <truncation> CpyCell= 1 ct_AN= 1 FNAN SNAN 4 0 <truncation> CpyCell= 1 ct_AN= 2 FNAN SNAN 4 0 <truncation> CpyCell= 1 ct_AN= 3 FNAN SNAN 4 0 <truncation> CpyCell= 1 ct_AN= 4 FNAN SNAN 4 0 <truncation> CpyCell= 1 ct_AN= 5 FNAN SNAN 4 0 TFNAN= 20 Average FNAN= 4.00000 TSNAN= 0 Average SNAN= 0.00000 <truncation> CpyCell= 2 ct_AN= 1 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 2 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 3 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 4 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 5 FNAN SNAN 4 0 TFNAN= 20 Average FNAN= 4.00000 TSNAN= 0 Average SNAN= 0.00000 <truncation> CpyCell= 2 ct_AN= 1 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 2 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 3 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 4 FNAN SNAN 4 0 <truncation> CpyCell= 2 ct_AN= 5 FNAN SNAN 4 0 <Check_System> The system is molecule. lattice vectors (bohr) A = 18.897259885789, 0.000000000000, 0.000000000000 B = 0.000000000000, 18.897259885789, 0.000000000000 C = 0.000000000000, 0.000000000000, 18.897259885789 reciprocal lattice vectors (bohr^-1) RA = 0.332491871581, 0.000000000000, 0.000000000000 RB = 0.000000000000, 0.332491871581, 0.000000000000 RB = 0.000000000000, 0.000000000000, 0.332491871581 Grid_Origin -9.300995100037 -9.300995100037 -9.300995100037 Cell_Volume = 6748.333037104149 (Bohr^3) GridVol = 0.025742847584 (Bohr^3) Grid_Origin -9.300995100037 -9.300995100037 -9.300995100037 Cell_Volume = 6748.333037104149 (Bohr^3) GridVol = 0.025742847584 (Bohr^3) <UCell_Box> Info. of cutoff energy and num. of grids lattice vectors (bohr) A = 18.897259885789, 0.000000000000, 0.000000000000 B = 0.000000000000, 18.897259885789, 0.000000000000 C = 0.000000000000, 0.000000000000, 18.897259885789 reciprocal lattice vectors (bohr^-1) RA = 0.332491871581, 0.000000000000, 0.000000000000 RB = 0.000000000000, 0.332491871581, 0.000000000000 RB = 0.000000000000, 0.000000000000, 0.332491871581 Required cutoff energy (Ryd) for 3D-grids = 120.0000 Used cutoff energy (Ryd) for 3D-grids = 113.2041, 113.2041, 113.2041 Num. of grids of a-, b-, and c-axes = 64, 64, 64 Grid_Origin -9.300995100037 -9.300995100037 -9.300995100037 Cell_Volume = 6748.333037104149 (Bohr^3) GridVol = 0.025742847584 (Bohr^3) Cell vectors (bohr) of the grid cell (gtv) gtv_a = 0.295269685715, 0.000000000000, 0.000000000000 gtv_b = 0.000000000000, 0.295269685715, 0.000000000000 gtv_c = 0.000000000000, 0.000000000000, 0.295269685715 |gtv_a| = 0.295269685715 |gtv_b| = 0.295269685715 |gtv_c| = 0.295269685715 Num. of grids overlapping with atom 1 = 20336 Num. of grids overlapping with atom 2 = 20346 Num. of grids overlapping with atom 3 = 20346 Num. of grids overlapping with atom 4 = 20346 Num. of grids overlapping with atom 5 = 20346
******************************************************* SCF calculation at MD = 1 *******************************************************
<MD= 1> Calculation of the overlap matrix <MD= 1> Calculation of the nonlocal matrix <MD= 1> Calculation of the VNA projector matrix
******************* MD= 1 SCF= 1 ******************* <Cluster> Solving the eigenvalue problem... Assertion failed in file ch3u_buffer.c at line 77: FALSE memcpy argument memory ranges overlap, dst_=0x251bc50 src_=0x251bc50 len_=512
internal ABORT - process 0
As I searched the internet and got to know that this is error from MPI. (I could not understand meaning fully)
Could I know how to fix this?
I am using mpich2-1.4.1p1 version and mkl library 11.1/059 version.
I also have openmpi, but I think that mkl library is compiled with mpich2.
| |