! Replaces DSPR2  8 procedures

    ! .       alpha    beta  lower  c

    ! 201     any     any     L    c < alpha a b(T) + alpha b a(T) + beta c
    ! 202     any     any          c < alpha a b(T) + alpha b a(T) + beta c
    ! 203     any     0.0     L    c < alpha a b(T) + alpha b a(T)
    ! 204     any     0.0          c < alpha a b(T) + alpha b a(T)
    ! 205     1.0     any     L    c < a b(T) + b a(T) + beta c
    ! 206     1.0     any          c < a b(T) + b a(T) + beta c
    ! 207     1.0     0.0     L    c < a b(T) + b a(T)
    ! 208     1.0     0.0          c < a b(T) + b a(T)


    MODULE dspr2_01

    CONTAINS
      SUBROUTINE dspr2_201(alpha,a,b,beta,c,lower)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        USE blas_keywords, ONLY : blas_lower_type
        ! .. Parameters ..
        REAL (wp), PARAMETER :: one = 1.0_wp
        ! .. Scalar Arguments ..
        REAL (wp), INTENT (IN) :: alpha, beta
        TYPE (blas_lower_type), INTENT (IN) :: lower
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        IF (beta/=one) c = beta*c
        CALL dspr2('L',n,alpha,a,1,b,1,c)
      END SUBROUTINE dspr2_201
    END MODULE dspr2_01

    MODULE dspr2_02

    CONTAINS
      SUBROUTINE dspr2_202(alpha,a,b,beta,c)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        ! .. Parameters ..
        REAL (wp), PARAMETER :: one = 1.0_wp
        ! .. Scalar Arguments ..
        REAL (wp), INTENT (IN) :: alpha, beta
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        IF (beta/=one) c = beta*c
        CALL dspr2('U',n,alpha,a,1,b,1,c)
      END SUBROUTINE dspr2_202
    END MODULE dspr2_02

    MODULE dspr2_03

    CONTAINS
      SUBROUTINE dspr2_203(alpha,a,b,c,lower)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        USE blas_keywords, ONLY : blas_lower_type
        ! .. Parameters ..
        REAL (wp), PARAMETER :: zero = 0.0_wp
        ! .. Scalar Arguments ..
        REAL (wp), INTENT (IN) :: alpha
        TYPE (blas_lower_type), INTENT (IN) :: lower
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        c = zero
        CALL dspr2('L',n,alpha,a,1,b,1,c)
      END SUBROUTINE dspr2_203
    END MODULE dspr2_03

    MODULE dspr2_04

    CONTAINS
      SUBROUTINE dspr2_204(alpha,a,b,c)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        ! .. Parameters ..
        REAL (wp), PARAMETER :: zero = 0.0_wp
        ! .. Scalar Arguments ..
        REAL (wp), INTENT (IN) :: alpha
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        c = zero
        CALL dspr2('U',n,alpha,a,1,b,1,c)
      END SUBROUTINE dspr2_204
    END MODULE dspr2_04

    MODULE dspr2_05

    CONTAINS
      SUBROUTINE dspr2_205(a,b,beta,c,lower)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        USE blas_keywords, ONLY : blas_lower_type
        ! .. Parameters ..
        REAL (wp), PARAMETER :: one = 1.0_wp
        ! .. Scalar Arguments ..
        REAL (wp), INTENT (IN) :: beta
        TYPE (blas_lower_type), INTENT (IN) :: lower
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        IF (beta/=one) c = beta*c
        CALL dspr2('L',n,one,a,1,b,1,c)
      END SUBROUTINE dspr2_205
    END MODULE dspr2_05

    MODULE dspr2_06

    CONTAINS
      SUBROUTINE dspr2_206(a,b,beta,c)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        ! .. Parameters ..
        REAL (wp), PARAMETER :: one = 1.0_wp
        ! .. Scalar Arguments ..
        REAL (wp), INTENT (IN) :: beta
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        IF (beta/=one) c = beta*c
        CALL dspr2('U',n,one,a,1,b,1,c)
      END SUBROUTINE dspr2_206
    END MODULE dspr2_06

    MODULE dspr2_07

    CONTAINS
      SUBROUTINE dspr2_207(a,b,c,lower)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        USE blas_keywords, ONLY : blas_lower_type
        ! .. Parameters ..
        REAL (wp), PARAMETER :: zero = 0.0_wp
        REAL (wp), PARAMETER :: one = 1.0_wp
        ! .. Scalar Arguments ..
        TYPE (blas_lower_type), INTENT (IN) :: lower
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        c = zero
        CALL dspr2('L',n,one,a,1,b,1,c)
      END SUBROUTINE dspr2_207
    END MODULE dspr2_07

    MODULE dspr2_08

    CONTAINS
      SUBROUTINE dspr2_208(a,b,c)
        ! .. Use Statements ..
        USE blas_kinds, ONLY : wp => dp
        ! .. Parameters ..
        REAL (wp), PARAMETER :: zero = 0.0_wp
        REAL (wp), PARAMETER :: one = 1.0_wp
        ! .. Array Arguments ..
        REAL (wp), INTENT (IN) :: a(:), b(:)
        REAL (wp), INTENT (INOUT) :: c(:)
        ! .. Local Scalar ..
        INTEGER :: n
        ! .. External Procedures ..
        EXTERNAL dspr2

        n = SIZE(a)
        c = zero
        CALL dspr2('U',n,one,a,1,b,1,c)
      END SUBROUTINE dspr2_208
    END MODULE dspr2_08