[Bug middle-end/27069] -ffast-math crash
nuno dot bandeira at ist dot utl dot pt
gcc-bugzilla@gcc.gnu.org
Fri Apr 7 00:32:00 GMT 2006
------- 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
More information about the Gcc-bugs
mailing list