sig
  type order = Blas.order = RowMajor | ColMajor
  type transpose = Blas.transpose = NoTrans | Trans | ConjTrans
  type uplo = Blas.uplo = Upper | Lower
  type diag = Blas.diag = NonUnit | Unit
  type side = Blas.side = Left | Right
  external dot : Vector_flat.vector -> Vector_flat.vector -> float
    = "ml_gsl_blas_ddot"
  external nrm2 : Vector_flat.vector -> float = "ml_gsl_blas_dnrm2"
  external asum : Vector_flat.vector -> float = "ml_gsl_blas_dasum"
  external iamax : Vector_flat.vector -> int = "ml_gsl_blas_idamax"
  external swap : Vector_flat.vector -> Vector_flat.vector -> unit
    = "ml_gsl_blas_dswap"
  external copy : Vector_flat.vector -> Vector_flat.vector -> unit
    = "ml_gsl_blas_dcopy"
  external axpy : float -> Vector_flat.vector -> Vector_flat.vector -> unit
    = "ml_gsl_blas_daxpy"
  external rot :
    Vector_flat.vector -> Vector_flat.vector -> float -> float -> unit
    = "ml_gsl_blas_drot"
  external scal : float -> Vector_flat.vector -> unit = "ml_gsl_blas_dscal"
  external gemv :
    Blas_flat.transpose ->
    alpha:float ->
    a:Matrix_flat.matrix ->
    x:Vector_flat.vector -> beta:float -> y:Vector_flat.vector -> unit
    = "ml_gsl_blas_dgemv_bc" "ml_gsl_blas_dgemv"
  external trmv :
    Blas_flat.uplo ->
    Blas_flat.transpose ->
    Blas_flat.diag -> a:Matrix_flat.matrix -> x:Vector_flat.vector -> unit
    = "ml_gsl_blas_dtrmv"
  external trsv :
    Blas_flat.uplo ->
    Blas_flat.transpose ->
    Blas_flat.diag -> a:Matrix_flat.matrix -> x:Vector_flat.vector -> unit
    = "ml_gsl_blas_dtrsv"
  external symv :
    Blas_flat.uplo ->
    alpha:float ->
    a:Matrix_flat.matrix ->
    x:Vector_flat.vector -> beta:float -> y:Vector_flat.vector -> unit
    = "ml_gsl_blas_dsymv_bc" "ml_gsl_blas_dsymv"
  external dger :
    alpha:float ->
    x:Vector_flat.vector ->
    y:Vector_flat.vector -> a:Matrix_flat.matrix -> unit = "ml_gsl_blas_dger"
  external syr :
    Blas_flat.uplo ->
    alpha:float -> x:Vector_flat.vector -> a:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dsyr"
  external syr2 :
    Blas_flat.uplo ->
    alpha:float ->
    x:Vector_flat.vector ->
    y:Vector_flat.vector -> a:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dsyr2"
  external gemm :
    ta:Blas_flat.transpose ->
    tb:Blas_flat.transpose ->
    alpha:float ->
    a:Matrix_flat.matrix ->
    b:Matrix_flat.matrix -> beta:float -> c:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dgemm_bc" "ml_gsl_blas_dgemm"
  external symm :
    Blas_flat.side ->
    Blas_flat.uplo ->
    alpha:float ->
    a:Matrix_flat.matrix ->
    b:Matrix_flat.matrix -> beta:float -> c:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dsymm_bc" "ml_gsl_blas_dsymm"
  external trmm :
    Blas_flat.side ->
    Blas_flat.uplo ->
    Blas_flat.transpose ->
    Blas_flat.diag ->
    alpha:float -> a:Matrix_flat.matrix -> b:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dtrmm_bc" "ml_gsl_blas_dtrmm"
  external trsm :
    Blas_flat.side ->
    Blas_flat.uplo ->
    Blas_flat.transpose ->
    Blas_flat.diag ->
    alpha:float -> a:Matrix_flat.matrix -> b:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dtrsm_bc" "ml_gsl_blas_dtrsm"
  external syrk :
    Blas_flat.uplo ->
    Blas_flat.transpose ->
    alpha:float ->
    a:Matrix_flat.matrix -> beta:float -> c:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dsyrk_bc" "ml_gsl_blas_dsyrk"
  external syr2k :
    Blas_flat.uplo ->
    Blas_flat.transpose ->
    alpha:float ->
    a:Matrix_flat.matrix ->
    b:Matrix_flat.matrix -> beta:float -> c:Matrix_flat.matrix -> unit
    = "ml_gsl_blas_dsyr2k_bc" "ml_gsl_blas_dsyr2k"
  module Complex :
    sig
      external dotu :
        Vector_complex_flat.vector ->
        Vector_complex_flat.vector -> Gsl_complex.complex
        = "ml_gsl_blas_zdotu"
      external dotc :
        Vector_complex_flat.vector ->
        Vector_complex_flat.vector -> Gsl_complex.complex
        = "ml_gsl_blas_zdotc"
      external nrm2 : Vector_complex_flat.vector -> float
        = "ml_gsl_blas_znrm2"
      external asum : Vector_complex_flat.vector -> float
        = "ml_gsl_blas_zasum"
      external iamax : Vector_complex_flat.vector -> int
        = "ml_gsl_blas_izamax"
      external swap :
        Vector_complex_flat.vector -> Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zswap"
      external copy :
        Vector_complex_flat.vector -> Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zcopy"
      external axpy :
        Gsl_complex.complex ->
        Vector_complex_flat.vector -> Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zaxpy"
      external scal :
        Gsl_complex.complex -> Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zscal"
      external zdscal : float -> Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zdscal"
      external gemv :
        Blas_flat.transpose ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        x:Vector_complex_flat.vector ->
        beta:Gsl_complex.complex -> y:Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zgemv_bc" "ml_gsl_blas_zgemv"
      external trmv :
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        Blas_flat.diag ->
        a:Matrix_complex_flat.matrix -> x:Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_ztrmv"
      external trsv :
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        Blas_flat.diag ->
        a:Matrix_complex_flat.matrix -> x:Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_ztrsv"
      external hemv :
        Blas_flat.uplo ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        x:Vector_complex_flat.vector ->
        beta:Gsl_complex.complex -> y:Vector_complex_flat.vector -> unit
        = "ml_gsl_blas_zhemv_bc" "ml_gsl_blas_zhemv"
      external geru :
        alpha:Gsl_complex.complex ->
        x:Vector_complex_flat.vector ->
        y:Vector_complex_flat.vector -> a:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zgeru"
      external gerc :
        alpha:Gsl_complex.complex ->
        x:Vector_complex_flat.vector ->
        y:Vector_complex_flat.vector -> a:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zgerc"
      external her :
        Blas_flat.uplo ->
        alpha:float ->
        x:Vector_complex_flat.vector -> a:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zher"
      external her2 :
        Blas_flat.uplo ->
        alpha:Gsl_complex.complex ->
        x:Vector_complex_flat.vector ->
        y:Vector_complex_flat.vector -> a:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zher2"
      external gemm :
        ta:Blas_flat.transpose ->
        tb:Blas_flat.transpose ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        b:Matrix_complex_flat.matrix ->
        beta:Gsl_complex.complex -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zgemm_bc" "ml_gsl_blas_zgemm"
      external symm :
        Blas_flat.side ->
        Blas_flat.uplo ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        b:Matrix_complex_flat.matrix ->
        beta:Gsl_complex.complex -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zsymm_bc" "ml_gsl_blas_zsymm"
      external syrk :
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        beta:Gsl_complex.complex -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zsyrk_bc" "ml_gsl_blas_zsyrk"
      external syr2k :
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        b:Matrix_complex_flat.matrix ->
        beta:Gsl_complex.complex -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zsyr2k_bc" "ml_gsl_blas_zsyr2k"
      external trmm :
        Blas_flat.side ->
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        Blas_flat.diag ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix -> b:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_ztrmm_bc" "ml_gsl_blas_ztrmm"
      external trsm :
        Blas_flat.side ->
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        Blas_flat.diag ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix -> b:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_ztrsm_bc" "ml_gsl_blas_ztrsm"
      external hemm :
        Blas_flat.side ->
        Blas_flat.uplo ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        b:Matrix_complex_flat.matrix ->
        beta:Gsl_complex.complex -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zhemm_bc" "ml_gsl_blas_zhemm"
      external herk :
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        alpha:float ->
        a:Matrix_complex_flat.matrix ->
        beta:float -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zherk_bc" "ml_gsl_blas_zherk"
      external her2k :
        Blas_flat.uplo ->
        Blas_flat.transpose ->
        alpha:Gsl_complex.complex ->
        a:Matrix_complex_flat.matrix ->
        b:Matrix_complex_flat.matrix ->
        beta:float -> c:Matrix_complex_flat.matrix -> unit
        = "ml_gsl_blas_zher2k_bc" "ml_gsl_blas_zher2k"
    end
end