blas.f90

 
module blas
    implicit none
 
    interface
        subroutine dgemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, &
            c, ldc)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: transa, transb
            integer(int32), intent(in) :: m, n, k, lda, ldb, ldc
            real(real64), intent(in) :: alpha, beta, a(lda,*), b(ldb,*)
            real(real64), intent(inout) :: c(ldc,*)
        end subroutine
 
        subroutine dgemv(trans, m, n, alpha, a, lda, x, incx, beta, y, incy)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: trans
            integer(int32), intent(in) :: m, n, lda, incx, incy
            real(real64), intent(in) :: alpha, beta, a(lda,*), x(*)
            real(real64), intent(inout) :: y(*)
        end subroutine
 
        subroutine zgemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, &
            c, ldc)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: transa, transb
            integer(int32), intent(in) :: m, n, k, lda, ldb, ldc
            complex(real64), intent(in) :: alpha, beta, a(lda,*), b(ldb,*)
            complex(real64), intent(inout) :: c(ldc,*)
        end subroutine
 
        subroutine zgemv(trans, m, n, alpha, a, lda, x, incx, beta, y, incy)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: trans
            integer(int32), intent(in) :: m, n, lda, incx, incy
            complex(real64), intent(in) :: alpha, beta, a(lda,*), x(*)
            complex(real64), intent(inout) :: y(*)
        end subroutine
 
        subroutine dtrsm(side, uplo, transa, diag, m, n, alpha, a, lda, b, ldb)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: side, uplo, transa, diag
            integer(int32), intent(in) :: m, n, lda, ldb
            real(real64), intent(in) :: alpha, a(lda,*)
            real(real64), intent(inout) :: b(ldb,*)
        end subroutine
 
        subroutine ztrsm(side, uplo, transa, diag, m, n, alpha, a, lda, b, ldb)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: side, uplo, transa, diag
            integer(int32), intent(in) :: m, n, lda, ldb
            complex(real64), intent(in) :: alpha, a(lda,*)
            complex(real64), intent(inout) :: b(ldb,*)
        end subroutine
 
        subroutine dtrsv(uplo, trans, diag, n, a, lda, x, incx)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: uplo, trans, diag
            integer(int32), intent(in) :: n, lda, incx
            real(real64), intent(in) :: a(lda,*)
            real(real64), intent(inout) :: x(*)
        end subroutine
 
        subroutine ztrsv(uplo, trans, diag, n, a, lda, x, incx)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: uplo, trans, diag
            integer(int32), intent(in) :: n, lda, incx
            complex(real64), intent(in) :: a(lda,*)
            complex(real64), intent(inout) :: x(*)
        end subroutine
 
        subroutine dscal(n, da, dx, incx)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, incx
            real(real64), intent(in) :: da
            real(real64), intent(inout) :: dx(*)
        end subroutine
 
        subroutine zscal(n, za, zx, incx)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, incx
            complex(real64), intent(in) :: za
            complex(real64), intent(inout) :: zx(*)
        end subroutine
 
        subroutine zdscal(n, da, zx, incx)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, incx
            real(real64), intent(in) :: da
            complex(real64), intent(inout) :: zx(*)
        end subroutine
 
        subroutine dgbmv(trans, m, n, kl, ku, alpha, a, lda, x, incx, beta, y, &
            incy)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: trans
            integer(int32), intent(in) :: m, n, kl, ku, lda, incx, incy
            real(real64), intent(in) :: alpha, beta, a(lda,*), x(*)
            real(real64), intent(inout) :: y(*)
        end subroutine
 
        subroutine zgbmv(trans, m, n, kl, ku, alpha, a, lda, x, incx, beta, y, &
            incy)
            use iso_fortran_env, only : int32, real64
            character, intent(in) :: trans
            integer(int32), intent(in) :: m, n, kl, ku, lda, incx, incy
            complex(real64), intent(in) :: alpha, beta, a(lda,*), x(*)
            complex(real64), intent(inout) :: y(*)
        end subroutine
 
        function ddot(n, dx, incx, dy, incy) result(rst)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, incx, incy
            real(real64), intent(in) :: dx(*), dy(*)
            real(real64) :: rst
        end function
    end interface
end module