DLASY2: INFO not set

Post here if you want to report a bug to the LAPACK team

DLASY2: INFO not set

Postby vchalupecky » Fri Jun 10, 2016 3:06 am

In dlasy2.f, it seems that INFO = 1 should be inserted before line 442.
vchalupecky
 
Posts: 6
Joined: Fri Jun 10, 2016 2:54 am

Re: DLASY2: INFO not set

Postby Julien Langou » Tue Jun 14, 2016 1:41 pm

Thanks Vladimir. The patch is below and was committed to SVN repository. Cheers, Julien.

Code: Select all
Index: dlasy2.f
===================================================================
--- dlasy2.f   (revision 1753)
+++ dlasy2.f   (working copy)
@@ -438,8 +438,10 @@
    80       CONTINUE
    90    CONTINUE
   100 CONTINUE
-      IF( ABS( T16( 4, 4 ) ).LT.SMIN )
-     $   T16( 4, 4 ) = SMIN
+      IF( ABS( T16( 4, 4 ) ).LT.SMIN ) THEN
+         INFO = 1
+         T16( 4, 4 ) = SMIN
+      END IF
       SCALE = ONE
       IF( ( EIGHT*SMLNUM )*ABS( BTMP( 1 ) ).GT.ABS( T16( 1, 1 ) ) .OR.
      $    ( EIGHT*SMLNUM )*ABS( BTMP( 2 ) ).GT.ABS( T16( 2, 2 ) ) .OR.
Index: slasy2.f
===================================================================
--- slasy2.f   (revision 1753)
+++ slasy2.f   (working copy)
@@ -438,8 +438,10 @@
    80       CONTINUE
    90    CONTINUE
   100 CONTINUE
-      IF( ABS( T16( 4, 4 ) ).LT.SMIN )
-     $   T16( 4, 4 ) = SMIN
+      IF( ABS( T16( 4, 4 ) ).LT.SMIN ) THEN
+         INFO = 1
+         T16( 4, 4 ) = SMIN
+      END IF
       SCALE = ONE
       IF( ( EIGHT*SMLNUM )*ABS( BTMP( 1 ) ).GT.ABS( T16( 1, 1 ) ) .OR.
      $    ( EIGHT*SMLNUM )*ABS( BTMP( 2 ) ).GT.ABS( T16( 2, 2 ) ) .OR.
Julien Langou
 
Posts: 821
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA


Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest