sig
  type order = RowMajor | ColMajor
  type transpose = NoTrans | Trans | ConjTrans
  type uplo = Upper | Lower
  type diag = NonUnit | Unit
  type side = Left | Right
  external dot : Gsl_vector.vector -> Gsl_vector.vector -> float
    = "ml_gsl_blas_ddot"
  external nrm2 : Gsl_vector.vector -> float = "ml_gsl_blas_dnrm2"
  external asum : Gsl_vector.vector -> float = "ml_gsl_blas_dasum"
  external iamax : Gsl_vector.vector -> int = "ml_gsl_blas_idamax"
  external swap : Gsl_vector.vector -> Gsl_vector.vector -> unit
    = "ml_gsl_blas_dswap"
  external copy : Gsl_vector.vector -> Gsl_vector.vector -> unit
    = "ml_gsl_blas_dcopy"
  external axpy : float -> Gsl_vector.vector -> Gsl_vector.vector -> unit
    = "ml_gsl_blas_daxpy"
  external rot :
    Gsl_vector.vector -> Gsl_vector.vector -> float -> float -> unit
    = "ml_gsl_blas_drot"
  external scal : float -> Gsl_vector.vector -> unit = "ml_gsl_blas_dscal"
  external gemv :
    Gsl_blas.transpose ->
    alpha:float ->
    a:Gsl_matrix.matrix ->
    x:Gsl_vector.vector -> beta:float -> y:Gsl_vector.vector -> unit
    = "ml_gsl_blas_dgemv_bc" "ml_gsl_blas_dgemv"
  external trmv :
    Gsl_blas.uplo ->
    Gsl_blas.transpose ->
    Gsl_blas.diag -> a:Gsl_matrix.matrix -> x:Gsl_vector.vector -> unit
    = "ml_gsl_blas_dtrmv"
  external trsv :
    Gsl_blas.uplo ->
    Gsl_blas.transpose ->
    Gsl_blas.diag -> a:Gsl_matrix.matrix -> x:Gsl_vector.vector -> unit
    = "ml_gsl_blas_dtrsv"
  external symv :
    Gsl_blas.uplo ->
    alpha:float ->
    a:Gsl_matrix.matrix ->
    x:Gsl_vector.vector -> beta:float -> y:Gsl_vector.vector -> unit
    = "ml_gsl_blas_dsymv_bc" "ml_gsl_blas_dsymv"
  external dger :
    alpha:float ->
    x:Gsl_vector.vector -> y:Gsl_vector.vector -> a:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dger"
  external syr :
    Gsl_blas.uplo ->
    alpha:float -> x:Gsl_vector.vector -> a:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dsyr"
  external syr2 :
    Gsl_blas.uplo ->
    alpha:float ->
    x:Gsl_vector.vector -> y:Gsl_vector.vector -> a:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dsyr2"
  external gemm :
    ta:Gsl_blas.transpose ->
    tb:Gsl_blas.transpose ->
    alpha:float ->
    a:Gsl_matrix.matrix ->
    b:Gsl_matrix.matrix -> beta:float -> c:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dgemm_bc" "ml_gsl_blas_dgemm"
  external symm :
    Gsl_blas.side ->
    Gsl_blas.uplo ->
    alpha:float ->
    a:Gsl_matrix.matrix ->
    b:Gsl_matrix.matrix -> beta:float -> c:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dsymm_bc" "ml_gsl_blas_dsymm"
  external trmm :
    Gsl_blas.side ->
    Gsl_blas.uplo ->
    Gsl_blas.transpose ->
    Gsl_blas.diag ->
    alpha:float -> a:Gsl_matrix.matrix -> b:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dtrmm_bc" "ml_gsl_blas_dtrmm"
  external trsm :
    Gsl_blas.side ->
    Gsl_blas.uplo ->
    Gsl_blas.transpose ->
    Gsl_blas.diag ->
    alpha:float -> a:Gsl_matrix.matrix -> b:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dtrsm_bc" "ml_gsl_blas_dtrsm"
  external syrk :
    Gsl_blas.uplo ->
    Gsl_blas.transpose ->
    alpha:float ->
    a:Gsl_matrix.matrix -> beta:float -> c:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dsyrk_bc" "ml_gsl_blas_dsyrk"
  external syr2k :
    Gsl_blas.uplo ->
    Gsl_blas.transpose ->
    alpha:float ->
    a:Gsl_matrix.matrix ->
    b:Gsl_matrix.matrix -> beta:float -> c:Gsl_matrix.matrix -> unit
    = "ml_gsl_blas_dsyr2k_bc" "ml_gsl_blas_dsyr2k"
  module Single :
    sig
      external sdsdot :
        alpha:float ->
        Gsl_vector.Single.vector -> Gsl_vector.Single.vector -> float
        = "ml_gsl_blas_sdsdot"
      external dsdot :
        Gsl_vector.Single.vector -> Gsl_vector.Single.vector -> float
        = "ml_gsl_blas_dsdot"
      external dot :
        Gsl_vector.Single.vector -> Gsl_vector.Single.vector -> float
        = "ml_gsl_blas_sdot"
      external nrm2 : Gsl_vector.Single.vector -> float = "ml_gsl_blas_snrm2"
      external asum : Gsl_vector.Single.vector -> float = "ml_gsl_blas_sasum"
      external iamax : Gsl_vector.Single.vector -> int = "ml_gsl_blas_isamax"
      external swap :
        Gsl_vector.Single.vector -> Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_sswap"
      external copy :
        Gsl_vector.Single.vector -> Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_scopy"
      external axpy :
        float -> Gsl_vector.Single.vector -> Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_saxpy"
      external rot :
        Gsl_vector.Single.vector ->
        Gsl_vector.Single.vector -> float -> float -> unit
        = "ml_gsl_blas_srot"
      external scal : float -> Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_sscal"
      external gemv :
        Gsl_blas.transpose ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix ->
        x:Gsl_vector.Single.vector ->
        beta:float -> y:Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_sgemv_bc" "ml_gsl_blas_sgemv"
      external trmv :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        a:Gsl_matrix.Single.matrix -> x:Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_strmv"
      external trsv :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        a:Gsl_matrix.Single.matrix -> x:Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_strsv"
      external symv :
        Gsl_blas.uplo ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix ->
        x:Gsl_vector.Single.vector ->
        beta:float -> y:Gsl_vector.Single.vector -> unit
        = "ml_gsl_blas_ssymv_bc" "ml_gsl_blas_ssymv"
      external dger :
        alpha:float ->
        x:Gsl_vector.Single.vector ->
        y:Gsl_vector.Single.vector -> a:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_sger"
      external syr :
        Gsl_blas.uplo ->
        alpha:float ->
        x:Gsl_vector.Single.vector -> a:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_ssyr"
      external syr2 :
        Gsl_blas.uplo ->
        alpha:float ->
        x:Gsl_vector.Single.vector ->
        y:Gsl_vector.Single.vector -> a:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_ssyr2"
      external gemm :
        ta:Gsl_blas.transpose ->
        tb:Gsl_blas.transpose ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix ->
        b:Gsl_matrix.Single.matrix ->
        beta:float -> c:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_sgemm_bc" "ml_gsl_blas_sgemm"
      external symm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix ->
        b:Gsl_matrix.Single.matrix ->
        beta:float -> c:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_ssymm_bc" "ml_gsl_blas_ssymm"
      external syrk :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix ->
        beta:float -> c:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_ssyrk_bc" "ml_gsl_blas_ssyrk"
      external syr2k :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix ->
        b:Gsl_matrix.Single.matrix ->
        beta:float -> c:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_ssyr2k_bc" "ml_gsl_blas_ssyr2k"
      external trmm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix -> b:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_strmm_bc" "ml_gsl_blas_strmm"
      external trsm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        alpha:float ->
        a:Gsl_matrix.Single.matrix -> b:Gsl_matrix.Single.matrix -> unit
        = "ml_gsl_blas_strsm_bc" "ml_gsl_blas_strsm"
    end
  module Complex :
    sig
      external dotu :
        Gsl_vector_complex.vector ->
        Gsl_vector_complex.vector -> Gsl_complex.complex
        = "ml_gsl_blas_zdotu"
      external dotc :
        Gsl_vector_complex.vector ->
        Gsl_vector_complex.vector -> Gsl_complex.complex
        = "ml_gsl_blas_zdotc"
      external nrm2 : Gsl_vector_complex.vector -> float
        = "ml_gsl_blas_znrm2"
      external asum : Gsl_vector_complex.vector -> float
        = "ml_gsl_blas_zasum"
      external iamax : Gsl_vector_complex.vector -> int
        = "ml_gsl_blas_izamax"
      external swap :
        Gsl_vector_complex.vector -> Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zswap"
      external copy :
        Gsl_vector_complex.vector -> Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zcopy"
      external axpy :
        Gsl_complex.complex ->
        Gsl_vector_complex.vector -> Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zaxpy"
      external scal :
        Gsl_complex.complex -> Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zscal"
      external zdscal : float -> Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zdscal"
      external gemv :
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        x:Gsl_vector_complex.vector ->
        beta:Gsl_complex.complex -> y:Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zgemv_bc" "ml_gsl_blas_zgemv"
      external trmv :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        a:Gsl_matrix_complex.matrix -> x:Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_ztrmv"
      external trsv :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        a:Gsl_matrix_complex.matrix -> x:Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_ztrsv"
      external hemv :
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        x:Gsl_vector_complex.vector ->
        beta:Gsl_complex.complex -> y:Gsl_vector_complex.vector -> unit
        = "ml_gsl_blas_zhemv_bc" "ml_gsl_blas_zhemv"
      external geru :
        alpha:Gsl_complex.complex ->
        x:Gsl_vector_complex.vector ->
        y:Gsl_vector_complex.vector -> a:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zgeru"
      external gerc :
        alpha:Gsl_complex.complex ->
        x:Gsl_vector_complex.vector ->
        y:Gsl_vector_complex.vector -> a:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zgerc"
      external her :
        Gsl_blas.uplo ->
        alpha:float ->
        x:Gsl_vector_complex.vector -> a:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zher"
      external her2 :
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        x:Gsl_vector_complex.vector ->
        y:Gsl_vector_complex.vector -> a:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zher2"
      external gemm :
        ta:Gsl_blas.transpose ->
        tb:Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        b:Gsl_matrix_complex.matrix ->
        beta:Gsl_complex.complex -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zgemm_bc" "ml_gsl_blas_zgemm"
      external symm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        b:Gsl_matrix_complex.matrix ->
        beta:Gsl_complex.complex -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zsymm_bc" "ml_gsl_blas_zsymm"
      external syrk :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        beta:Gsl_complex.complex -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zsyrk_bc" "ml_gsl_blas_zsyrk"
      external syr2k :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        b:Gsl_matrix_complex.matrix ->
        beta:Gsl_complex.complex -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zsyr2k_bc" "ml_gsl_blas_zsyr2k"
      external trmm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix -> b:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_ztrmm_bc" "ml_gsl_blas_ztrmm"
      external trsm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix -> b:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_ztrsm_bc" "ml_gsl_blas_ztrsm"
      external hemm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        b:Gsl_matrix_complex.matrix ->
        beta:Gsl_complex.complex -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zhemm_bc" "ml_gsl_blas_zhemm"
      external herk :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:float ->
        a:Gsl_matrix_complex.matrix ->
        beta:float -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zherk_bc" "ml_gsl_blas_zherk"
      external her2k :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.matrix ->
        b:Gsl_matrix_complex.matrix ->
        beta:float -> c:Gsl_matrix_complex.matrix -> unit
        = "ml_gsl_blas_zher2k_bc" "ml_gsl_blas_zher2k"
    end
  module Complex_Single :
    sig
      external dotu :
        Gsl_vector_complex.Single.vector ->
        Gsl_vector_complex.Single.vector -> Gsl_complex.complex
        = "ml_gsl_blas_cdotu"
      external dotc :
        Gsl_vector_complex.Single.vector ->
        Gsl_vector_complex.Single.vector -> Gsl_complex.complex
        = "ml_gsl_blas_cdotc"
      external nrm2 : Gsl_vector_complex.Single.vector -> float
        = "ml_gsl_blas_scnrm2"
      external asum : Gsl_vector_complex.Single.vector -> float
        = "ml_gsl_blas_scasum"
      external iamax : Gsl_vector_complex.Single.vector -> int
        = "ml_gsl_blas_icamax"
      external swap :
        Gsl_vector_complex.Single.vector ->
        Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_cswap"
      external copy :
        Gsl_vector_complex.Single.vector ->
        Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_ccopy"
      external axpy :
        Gsl_complex.complex ->
        Gsl_vector_complex.Single.vector ->
        Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_caxpy"
      external scal :
        Gsl_complex.complex -> Gsl_vector_complex.Single.vector -> unit
        = "ml_gsl_blas_cscal"
      external csscal : float -> Gsl_vector_complex.Single.vector -> unit
        = "ml_gsl_blas_csscal"
      external gemv :
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        x:Gsl_vector_complex.Single.vector ->
        beta:Gsl_complex.complex ->
        y:Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_cgemv_bc"
        "ml_gsl_blas_cgemv"
      external trmv :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        a:Gsl_matrix_complex.Single.matrix ->
        x:Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_ctrmv"
      external trsv :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        a:Gsl_matrix_complex.Single.matrix ->
        x:Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_ctrsv"
      external hemv :
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        x:Gsl_vector_complex.Single.vector ->
        beta:Gsl_complex.complex ->
        y:Gsl_vector_complex.Single.vector -> unit = "ml_gsl_blas_chemv_bc"
        "ml_gsl_blas_chemv"
      external geru :
        alpha:Gsl_complex.complex ->
        x:Gsl_vector_complex.Single.vector ->
        y:Gsl_vector_complex.Single.vector ->
        a:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_cgeru"
      external gerc :
        alpha:Gsl_complex.complex ->
        x:Gsl_vector_complex.Single.vector ->
        y:Gsl_vector_complex.Single.vector ->
        a:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_cgerc"
      external her :
        Gsl_blas.uplo ->
        alpha:float ->
        x:Gsl_vector_complex.Single.vector ->
        a:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_cher"
      external her2 :
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        x:Gsl_vector_complex.Single.vector ->
        y:Gsl_vector_complex.Single.vector ->
        a:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_cher2"
      external gemm :
        ta:Gsl_blas.transpose ->
        tb:Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix ->
        beta:Gsl_complex.complex ->
        c:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_cgemm_bc"
        "ml_gsl_blas_cgemm"
      external symm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix ->
        beta:Gsl_complex.complex ->
        c:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_csymm_bc"
        "ml_gsl_blas_csymm"
      external syrk :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        beta:Gsl_complex.complex ->
        c:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_csyrk_bc"
        "ml_gsl_blas_csyrk"
      external syr2k :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix ->
        beta:Gsl_complex.complex ->
        c:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_csyr2k_bc"
        "ml_gsl_blas_csyr2k"
      external trmm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_ctrmm_bc"
        "ml_gsl_blas_ctrmm"
      external trsm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        Gsl_blas.diag ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_ctrsm_bc"
        "ml_gsl_blas_ctrsm"
      external hemm :
        Gsl_blas.side ->
        Gsl_blas.uplo ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix ->
        beta:Gsl_complex.complex ->
        c:Gsl_matrix_complex.Single.matrix -> unit = "ml_gsl_blas_chemm_bc"
        "ml_gsl_blas_chemm"
      external herk :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:float ->
        a:Gsl_matrix_complex.Single.matrix ->
        beta:float -> c:Gsl_matrix_complex.Single.matrix -> unit
        = "ml_gsl_blas_cherk_bc" "ml_gsl_blas_cherk"
      external her2k :
        Gsl_blas.uplo ->
        Gsl_blas.transpose ->
        alpha:Gsl_complex.complex ->
        a:Gsl_matrix_complex.Single.matrix ->
        b:Gsl_matrix_complex.Single.matrix ->
        beta:float -> c:Gsl_matrix_complex.Single.matrix -> unit
        = "ml_gsl_blas_cher2k_bc" "ml_gsl_blas_cher2k"
    end
end