------- Comment #8 from nuno dot bandeira at ist dot utl dot pt 2006-04-07 00:32 ------- Subject: Re: -ffast-math crash
This is what I get from the compiler: $ gfortran -I../SOURCE -c -fdefault-real-8 -O3 -fforce-addr -march=pentium4 -fcray-pointer -ffast-math ./graden.f -o ./graden.o ./graden.f: In function 'graden': ./graden.f:2: internal compiler error: in find_lattice_value, at tree-complex.c:133 Please submit a full bug report, with preprocessed source if appropriate. See <URL:http://gcc.gnu.org/bugs.html> for instructions. Attached are also the include files for the compilation. Thank you all for such a quick reply ! The code is not mine by the way hence my ignorance... C ================================================================== C == DYNAMIC ALLOCATION OF PERMANENT ARRAYS == C ================================================================== C == TSHELS = FALSE if all TSHEL(IS) false == C ==--------------------------------------------------------------== LOGICAL TSHELS,TSHEL(MAXSP) COMMON/CPPSW/ TSHELS,TSHEL C ================================================================== C == INYH(3,NHG) coordinates +NHI (I=1,2,3) for G-vectors == C == IGL(NHG) the index of the shell for each G-vector == C == NZH(NHG) index in PSI or RHO array (IG1>=0) == C == INDZ(NHG) index in PSI or RHO array (IG1<=0) == C == ISPTR(NHGL+1) last index IG for the shell == C == INDZS(NGW) index for G-compon. of wavefunction in PSI (I1<0)== C == NZH(NGW) index for G-compon. of wavefunction in PSI (I1>0)== C ==--------------------------------------------------------------== INTEGER INYH(3,NHG),IGL(NHG),NZH(NHG),INDZ(NHG), & ISPTR(*),INDZS(NGW),NZHS(NGW) POINTER (IP_INYH,INYH), & (IP_IGL,IGL), & (IP_NZH,NZH), & (IP_INDZ,INDZ), & (IP_ISPTR,ISPTR), & (IP_INDZS,INDZS), & (IP_NZHS,NZHS) COMMON/CPPTI/ IP_INYH,IP_IGL,IP_NZH,IP_INDZ,IP_ISPTR,IP_INDZS, & IP_NZHS C ================================================================== C == HG(1:NHG) Square norm of G-vectors == C == GK(1:3,1:NHG) Components of G-vectors == C == GL(1:NHGL) Square norm of G-vectors for each shell == C == VPS(1:NSP,1:NHG) Local pseudopotential per species in G space== C == RHOPS: smeared ionic charge density in G space == C == ionic point charges replaced by Gaussian charge == C == distributions (see ener.inc) calculated in PUTPS == C == TWNL(1:NGW,1:NGH(IS),1:NSP) Non-Local projectors array == C == for each G-components (Kleinman-Bylander form) == C == USED BY VANDERBILT PSEUDOPOTENTIALS == C == QRAD == C == YLMB(NHGK,LPMAX,NKPNT) Initialized in PUTWNL == C ==--------------------------------------------------------------== REAL*8 HG(NHG),GK(3,NHG),GL(*), & VPS(NSX,NHG),RHOPS(NSX,NHG), & TWNL,QRAD,TWNLS,YLMB POINTER (IP_HG,HG), & (IP_GK,GK), & (IP_GL,GL), & (IP_VPS,VPS), & (IP_RHOPS,RHOPS), & (IP_TWNL,TWNL), & (IP_QRAD,QRAD), & (IP_TWNLS,TWNLS), & (IP_YLMB,YLMB) COMMON/CPPTR/ IP_HG,IP_GK,IP_GL,IP_VPS,IP_RHOPS, & IP_TWNL,IP_TWNLS,IP_QRAD,IP_YLMB C ================================================================== C == Dimension of HGPOT (for isolated system -- HIP) == C ==--------------------------------------------------------------== INTEGER NR1H,NR2H,NR3H,NR3PL COMMON/ISOSI/ NR1H,NR2H,NR3H,NR3PL C ================================================================== REAL*8 HGPOT(NR1S+1,NR2S+1,NR3PL),HIPZ(NHG) COMPLEX*16 SCG(NHG),SCGX(NHG) POINTER (IP_HGPOT,HGPOT),(IP_HIPZ,HIPZ) POINTER (IP_SCG,SCG),(IP_SCGX,SCGX) COMMON/ISOSP/ IP_HGPOT,IP_HIPZ,IP_SCG,IP_SCGX C ================================================================== C ================================================================== C == INCLUDE FILE OF NUMERIC CONSTANTS (SEE SETCNST.F) == C ==--------------------------------------------------------------== C == UIMAG = DCMPLX(0.D0,1.D0) == C == PI PI NUMBER == C == FPI 4*PI == C == RY RYDBERG IN ELECTRON-VOLT == C == FACTEM 1 HARTREE IN KELVIN == C == SCMASS PROTON MASS (IN ATOMIC UNITS = 1822*ELECTRON MASS) == C == FBOHR ANGSTROM TO ATOMIC UNITS (BOHR) == C == AU_K ATOMIC UNITS TO KBAR == C == KB_AU KBAR TO ATOMIC UNITS == C ================================================================== COMPLEX*16 UIMAG REAL*8 PI,FPI,RY,FACTEM,SCMASS,FBOHR,AU_KB,KB_AU COMMON/CNST/ UIMAG,PI,FPI,RY,FACTEM,SCMASS,FBOHR, * AU_KB,KB_AU C ================================================================== C ================================================================== INTEGER MG,MS,MZ,NAUX1,NAUX2 REAL*8 AUX1,AUX2 COMPLEX*16 XF,YF,ZF POINTER (IP_AUX1,AUX1), (IP_AUX2,AUX2), * (IP_MG,MG), (IP_MS,MS), (IP_MZ,MZ), * (IP_XF,XF), (IP_YF,YF), (IP_ZF,ZF) COMMON/CFFT_X/ IP_AUX1,IP_AUX2, * IP_MG,IP_MS,IP_MZ, * IP_XF,IP_YF,IP_ZF COMMON/CFFT_Y/ NAUX1,NAUX2 C ================================================================== INTEGER MAXFFT,NR1M,NR2M,NR3M,KR1M,KR2M,KR3M,NHRM,NGRM, * KR2MIN,KR2MAX,KR3MIN,KR3MAX, * MXY(MAXCPU),MAXRPT COMMON/DFFT_X/ MAXFFT,NR1M,NR2M,NR3M,KR1M,KR2M,KR3M,NHRM,NGRM, * KR2MIN,KR2MAX,KR3MIN,KR3MAX, * MXY,MAXRPT C ================================================================== C == GATHER/SCATTER ARRAYS == C ================================================================== INTEGER MSP POINTER (IP_MSP,MSP) COMMON/GSAR/ IP_MSP C ================================================================== C NEW GENERAL PARALLEL FFT CODE C ================================================================== INTEGER MSRAYS,MFRAYS,LLR1 INTEGER QR1S,QR2S,QR3S,QR1,QR2,QR3 INTEGER LR1S,LR2S,LR3S,LR1,LR2,LR3 INTEGER QR2MAX,QR2MIN,QR3MAX,QR3MIN INTEGER LSRM,LFRM,LR1M,LMSQ,MAXFFTN INTEGER JGW,JGWS,JHG,JHGS COMMON/FFTN1/ MSRAYS,MFRAYS,LLR1,QR1S,QR2S,QR3S,QR1,QR2,QR3, * LR1S,LR2S,LR3S,LR1,LR2,LR3, * QR2MAX,QR2MIN,QR3MAX,QR3MIN,LSRM,LFRM,LMSQ,MAXFFTN, * LR1M,JGW,JGWS,JHG,JHGS INTEGER LRXPL(0:MAXCPU,2) INTEGER SP5(0:MAXCPU),SP8(0:MAXCPU),SP9(0:MAXCPU) COMMON/FFTN2/ LRXPL,SP5,SP8,SP9 INTEGER MSQS,MSQF POINTER (IP_MSQS,MSQS),(IP_MSQF,MSQF) INTEGER NZFF(*),NZFS(*),INZF(*),INZS(*),INZH(3,*) POINTER (IP_NZFF,NZFF),(IP_NZFS,NZFS) POINTER (IP_INZF,INZF),(IP_INZS,INZS),(IP_INZH,INZH) COMMON/FFTN3/ IP_MSQS,IP_MSQF,IP_NZFF,IP_NZFS,IP_INZF,IP_INZS, * IP_INZH C ================================================================== C POOL C ================================================================== INTEGER FFTPOOLSIZE PARAMETER (FFTPOOLSIZE=3) INTEGER FFTPOOL COMMON/FFTP0/ FFTPOOL INTEGER FPOOLV(28,FFTPOOLSIZE) COMMON/FFTP1/ FPOOLV INTEGER LRXPOOL(0:MAXCPU,2,FFTPOOLSIZE) INTEGER SPM(9,0:MAXCPU,FFTPOOLSIZE) COMMON/FFTP2/ LRXPOOL,SPM INTEGER LMSQMAX,LNZF,LNZS COMMON/FFTP3/ LMSQMAX,LNZF,LNZS INTEGER MSQSPOOL,MSQFPOOL POINTER (IP_MSQSPOOL,MSQSPOOL),(IP_MSQFPOOL,MSQFPOOL) COMMON/FFTP4/ IP_MSQSPOOL,IP_MSQFPOOL INTEGER NZFFP(LNZF,*),NZFSP(LNZS,*), * INZFP(LNZF,*),INZSP(LNZS,*),INZHP(3,LNZF,*) POINTER (IP_NZFFP,NZFFP),(IP_NZFSP,NZFSP) POINTER (IP_INZFP,INZFP),(IP_INZSP,INZSP),(IP_INZHP,INZHP) COMMON/FFTP5/ IP_NZFFP,IP_NZFSP,IP_INZFP,IP_INZSP,IP_INZHP C ================================================================== -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=27069
