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