[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