sig
  type double_mat_bigarr =
      (float, Bigarray.float64_elt, Bigarray.c_layout) Bigarray.Array2.t
  type matrix = Matrix.double_mat_bigarr
  val create : ?init:float -> int -> int -> Matrix.matrix
  val dims : Matrix.matrix -> int * int
  val of_array : float array -> int -> int -> Matrix.matrix
  val of_arrays : float array array -> Matrix.matrix
  val to_array : Matrix.matrix -> float array
  val to_arrays : Matrix.matrix -> float array array
  val get : Matrix.matrix -> int -> int -> float
  val set : Matrix.matrix -> int -> int -> float -> unit
  val set_all : Matrix.matrix -> float -> unit
  val set_zero : Matrix.matrix -> unit
  val set_id : Matrix.matrix -> unit
  val memcpy : src:Matrix.matrix -> dst:Matrix.matrix -> unit
  val copy : Matrix.matrix -> Matrix.matrix
  val row : Matrix.matrix -> int -> Vector.vector
  external add : Matrix.matrix -> Matrix.matrix -> unit = "ml_gsl_matrix_add"
  external sub : Matrix.matrix -> Matrix.matrix -> unit = "ml_gsl_matrix_sub"
  external mul_elements : Matrix.matrix -> Matrix.matrix -> unit
    = "ml_gsl_matrix_mul"
  external div_elements : Matrix.matrix -> Matrix.matrix -> unit
    = "ml_gsl_matrix_div"
  external scale : Matrix.matrix -> float -> unit = "ml_gsl_matrix_scale"
  external add_constant : Matrix.matrix -> float -> unit
    = "ml_gsl_matrix_add_constant"
  external add_diagonal : Matrix.matrix -> float -> unit
    = "ml_gsl_matrix_add_diagonal"
  external is_null : Matrix.matrix -> bool = "ml_gsl_matrix_isnull"
  external swap_rows : Matrix.matrix -> int -> int -> unit
    = "ml_gsl_matrix_swap_rows"
  external swap_columns : Matrix.matrix -> int -> int -> unit
    = "ml_gsl_matrix_swap_columns"
  external swap_rowcol : Matrix.matrix -> int -> int -> unit
    = "ml_gsl_matrix_swap_rowcol"
  external transpose : Matrix.matrix -> Matrix.matrix -> unit
    = "ml_gsl_matrix_transpose_memcpy"
  external transpose_in_place : Matrix.matrix -> unit
    = "ml_gsl_matrix_transpose"
  module Single :
    sig
      type float_mat_bigarr =
          (float, Bigarray.float32_elt, Bigarray.c_layout) Bigarray.Array2.t
      type matrix = Matrix.Single.float_mat_bigarr
      val create : ?init:float -> int -> int -> Matrix.Single.matrix
      val dims : Matrix.Single.matrix -> int * int
      val of_array : float array -> int -> int -> Matrix.Single.matrix
      val of_arrays : float array array -> Matrix.Single.matrix
      val to_array : Matrix.Single.matrix -> float array
      val to_arrays : Matrix.Single.matrix -> float array array
      val get : Matrix.Single.matrix -> int -> int -> float
      val set : Matrix.Single.matrix -> int -> int -> float -> unit
      val set_all : Matrix.Single.matrix -> float -> unit
      val set_zero : Matrix.Single.matrix -> unit
      val set_id : Matrix.Single.matrix -> unit
      val memcpy :
        src:Matrix.Single.matrix -> dst:Matrix.Single.matrix -> unit
      val copy : Matrix.Single.matrix -> Matrix.Single.matrix
      val row : Matrix.Single.matrix -> int -> Vector.Single.vector
      external add : Matrix.Single.matrix -> Matrix.Single.matrix -> unit
        = "ml_gsl_matrix_float_add"
      external sub : Matrix.Single.matrix -> Matrix.Single.matrix -> unit
        = "ml_gsl_matrix_float_sub"
      external mul_elements :
        Matrix.Single.matrix -> Matrix.Single.matrix -> unit
        = "ml_gsl_matrix_float_mul"
      external div_elements :
        Matrix.Single.matrix -> Matrix.Single.matrix -> unit
        = "ml_gsl_matrix_float_div"
      external scale : Matrix.Single.matrix -> float -> unit
        = "ml_gsl_matrix_float_scale"
      external add_constant : Matrix.Single.matrix -> float -> unit
        = "ml_gsl_matrix_float_add_constant"
      external add_diagonal : Matrix.Single.matrix -> float -> unit
        = "ml_gsl_matrix_float_add_diagonal"
      external is_null : Matrix.Single.matrix -> bool
        = "ml_gsl_matrix_float_isnull"
      external swap_rows : Matrix.Single.matrix -> int -> int -> unit
        = "ml_gsl_matrix_float_swap_rows"
      external swap_columns : Matrix.Single.matrix -> int -> int -> unit
        = "ml_gsl_matrix_float_swap_columns"
      external swap_rowcol : Matrix.Single.matrix -> int -> int -> unit
        = "ml_gsl_matrix_float_swap_rowcol"
      external transpose :
        Matrix.Single.matrix -> Matrix.Single.matrix -> unit
        = "ml_gsl_matrix_float_transpose_memcpy"
      external transpose_in_place : Matrix.Single.matrix -> unit
        = "ml_gsl_matrix_float_transpose"
    end
end