Add missing subroutine parameter checks for linear solver

Author: sh-zheng

TESTING/LIN/derrky.f

@@ -84,7 +84,7 @@ SUBROUTINE DERRKY( PATH, NUNIT )
       EXTERNAL           LSAMEN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ALAESM, CHKXER, DKYTRI2X, DKYTF2
+      EXTERNAL           ALAESM, CHKXER, DKYTRI2X, DKYTF2,
      $                   DKYTRF, DKYTRI, DKYTRS, DKYTRI2
 *     ..
 *     .. Scalars in Common ..

TESTING/LIN/derrvx.f

@@ -89,7 +89,7 @@ SUBROUTINE DERRVX( PATH, NUNIT )
      $                   DGTSVX, DPBSV, DPBSVX, DPOSV, DPOSVX, DPPSV,
      $                   DPPSVX, DPTSV, DPTSVX, DSPSV, DSPSVX, DSYSV,
      $                   DSYSV_AA, DSYSV_RK, DSYSV_ROOK, DSYSVX,
-     $                   DSYSV_AA_2STAGE
+     $                   DSYSV_AA_2STAGE, DKYSV
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
@@ -634,6 +634,33 @@ SUBROUTINE DERRVX( PATH, NUNIT )
          CALL DSYSVX( 'N', 'U', 2, 0, A, 2, AF, 2, IP, B, 2, X, 2,
      $                RCOND, R1, R2, W, 3, IW, INFO )
          CALL CHKXER( 'DSYSVX', INFOT, NOUT, LERR, OK )
+*
+      ELSE IF( LSAMEN( 2, C2, 'KY' ) ) THEN
+*
+*        DKYSV
+*
+         SRNAMT = 'DKYSV '
+         INFOT = 1
+         CALL DKYSV( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL DKYSV( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL DKYSV( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL DKYSV( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
+         CALL CHKXER( 'DKYSV_ROOK', INFOT, NOUT, LERR, OK )
+         INFOT = 8
+         CALL DKYSV( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 10
+         CALL DKYSV( 'U', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
+         CALL CHKXER( 'DKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 10
+         CALL DKYSV( 'U', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
+         CALL CHKXER( 'DKYSV ', INFOT, NOUT, LERR, OK )
 *
       ELSE IF( LSAMEN( 2, C2, 'SR' ) ) THEN
 *

TESTING/LIN/serrky.f

@@ -84,7 +84,7 @@ SUBROUTINE SERRKY( PATH, NUNIT )
       EXTERNAL           LSAMEN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ALAESM, CHKXER, SKYTRI2X, SKYTF2
+      EXTERNAL           ALAESM, CHKXER, SKYTRI2X, SKYTF2,
      $                   SKYTRF, SKYTRI, SKYTRS, SKYTRI2
 *     ..
 *     .. Scalars in Common ..

TESTING/LIN/serrvx.f

@@ -89,7 +89,7 @@ SUBROUTINE SERRVX( PATH, NUNIT )
      $                   SGTSVX, SPBSV, SPBSVX, SPOSV, SPOSVX, SPPSV,
      $                   SPPSVX, SPTSV, SPTSVX, SSPSV, SSPSVX, SSYSV,
      $                   SSYSV_AA, SSYSV_RK, SSYSV_ROOK, SSYSVX,
-     $                   SSYSV_AA_2STAGE
+     $                   SSYSV_AA_2STAGE, SKYSV
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
@@ -635,6 +635,32 @@ SUBROUTINE SERRVX( PATH, NUNIT )
      $                RCOND, R1, R2, W, 3, IW, INFO )
          CALL CHKXER( 'SSYSVX', INFOT, NOUT, LERR, OK )
 *
+      ELSE IF( LSAMEN( 2, C2, 'KY' ) ) THEN
+*
+*        SKYSV
+*
+         SRNAMT = 'SKYSV '
+         INFOT = 1
+         CALL SKYSV( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL SKYSV( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL SKYSV( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL SKYSV( 'U', 2, 0, A, 1, IP, B, 2, W, 1, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 8
+         CALL SKYSV( 'U', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 10
+         CALL SKYSV( 'U', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
+         INFOT = 10
+         CALL SKYSV( 'U', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
+         CALL CHKXER( 'SKYSV ', INFOT, NOUT, LERR, OK )
 *
       ELSE IF( LSAMEN( 2, C2, 'SR' ) ) THEN
 *