Merge pull request #455 from timleslie/fix-rotmg-bug

[BLAS] Fix missing logic path in xrotmg.f

Author: langou

BLAS/SRC/drotmg.f

@@ -152,6 +152,19 @@ SUBROUTINE DROTMG(DD1,DD2,DX1,DY1,DPARAM)
              DD1 = DD1/DU
              DD2 = DD2/DU
              DX1 = DX1*DU
+           ELSE
+*            This code path if here for safety. We do not expect this
+*            condition to ever hold except in edge cases with rounding
+*            errors. See DOI: 10.1145/355841.355847
+             DFLAG = -ONE
+             DH11 = ZERO
+             DH12 = ZERO
+             DH21 = ZERO
+             DH22 = ZERO
+*
+             DD1 = ZERO
+             DD2 = ZERO
+             DX1 = ZERO
            END IF
          ELSE
 

BLAS/SRC/srotmg.f

@@ -152,6 +152,19 @@ SUBROUTINE SROTMG(SD1,SD2,SX1,SY1,SPARAM)
              SD1 = SD1/SU
              SD2 = SD2/SU
              SX1 = SX1*SU
+           ELSE
+*            This code path if here for safety. We do not expect this
+*            condition to ever hold except in edge cases with rounding
+*            errors. See DOI: 10.1145/355841.355847
+             SFLAG = -ONE
+             SH11 = ZERO
+             SH12 = ZERO
+             SH21 = ZERO
+             SH22 = ZERO
+*
+             SD1 = ZERO
+             SD2 = ZERO
+             SX1 = ZERO
            END IF
          ELSE