[Bug fortran/29621] lapack runs into infinite loop with -03

["fkar at nemesis-project dot org" <gcc-bugzilla at gcc dot gnu dot org>]

------- Comment #2 from fkar at nemesis-project dot org  2006-10-27 21:14 -------
I just tried to submit a reduced test case. Please reconsider this bug with
these two cases (one linking with a Fortran program and one with a C/C++ one):

===========================================================================
Test case #1 [Fortran with DOUBLE PRECISION]
===========================================================================

      PROGRAM TESTCASE
      IMPLICIT NONE

      DOUBLE PRECISION A(3,3)
      DOUBLE PRECISION X(3)
      DOUBLE PRECISION WORK(9)
      INTEGER INFO

      A(1,1)= 1.        
      A(2,1)= 0.        
      A(3,1)= 0.        
      A(1,2)= 0.
      A(2,2)= 1.
      A(3,2)= 0.
      A(1,3)= 0.
      A(2,3)= 0.
      A(3,3)= 1.

      CALL DSYEV('N','L',3,A,3,X,WORK,3*3,INFO);
      END PROGRAM

===========================================================================
Test case #2 [C/C++]
===========================================================================
 extern "C" void dsyev_(char* JOBZ,char* UPLO,int* N,double* A,int* LDA,
                        double* W,double* WORK,int* LWORK,int* INFO);

int main()
{
        int N=3;
        double x[3];
        double a[3*3];
        a[0] =  1; a[3] =   0; a[6] =  0;
        a[1] =  0; a[4] =   1; a[7] =  0;
        a[2] =  0; a[5] =   0; a[8] =  1;
        char JOBZ='N';
        char UPLO='L';
        int LDA=N;
        int LWORK=3*3;
        double WORK[LWORK];
        int INFO;
        dsyev_(&JOBZ,&UPLO,&N,&a[0],&LDA,&x[0],&WORK[0],&LWORK,&INFO);
        return 0;
}

The result is in both cases identical: an infinite loop.

Kind regards,
F.E. Karaoulanis (www.nemesis-project.org)


-- 

fkar at nemesis-project dot org changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
             Status|RESOLVED                    |UNCONFIRMED
         Resolution|INVALID                     |


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=29621