Page 1 of 1

DLASY2: INFO not set

PostPosted: Fri Jun 10, 2016 3:06 am
by vchalupecky
In dlasy2.f, it seems that INFO = 1 should be inserted before line 442.

Re: DLASY2: INFO not set

PostPosted: Tue Jun 14, 2016 1:41 pm
by Julien Langou
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.