fix data type
authorMartin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Fri, 29 Jan 2021 09:45:36 +0000 (10:45 +0100)
committerGitHub <noreply@github.com>
Fri, 29 Jan 2021 09:45:36 +0000 (10:45 +0100)
lapack-netlib/SRC/chgeqz.f

index 0d37879..4725e71 100644 (file)
                TEMP = MAX( TEMP, ABS1( X ) )
                Y = TEMP*SQRT( ( X / TEMP )**2+( CTEMP / TEMP )**2 )
                IF( TEMP2.GT.ZERO ) THEN
-                  IF( DBLE( X / TEMP2 )*DBLE( Y )+
-     $                DIMAG( X / TEMP2 )*DIMAG( Y ).LT.ZERO )Y = -Y
+                  IF( REAL( X / TEMP2 )*REAL( Y )+
+     $                AIMAG( X / TEMP2 )*AIMAG( Y ).LT.ZERO )Y = -Y
                END IF
                SHIFT = SHIFT - CTEMP*CLADIV( CTEMP, ( X+Y ) )
             END IF