sig
  type double_mat_flat = private {
    data : float array;
    off : int;
    dim1 : int;
    dim2 : int;
    tda : int;
  }
  type matrix = Gsl_matrix_flat.double_mat_flat
  val create : ?init:float -> int -> int -> Gsl_matrix_flat.matrix
  val dims : Gsl_matrix_flat.matrix -> int * int
  val of_array : float array -> int -> int -> Gsl_matrix_flat.matrix
  val of_arrays : float array array -> Gsl_matrix_flat.matrix
  val to_array : matrix -> float array
  val to_arrays : Gsl_matrix_flat.matrix -> float array array
  val to_array : Gsl_matrix_flat.matrix -> float array
  val get : Gsl_matrix_flat.matrix -> int -> int -> float
  val set : Gsl_matrix_flat.matrix -> int -> int -> float -> unit
  val set_all : Gsl_matrix_flat.matrix -> float -> unit
  val set_zero : Gsl_matrix_flat.matrix -> unit
  val set_id : Gsl_matrix_flat.matrix -> unit
  val memcpy :
    src:Gsl_matrix_flat.matrix -> dst:Gsl_matrix_flat.matrix -> unit
  val copy : Gsl_matrix_flat.matrix -> Gsl_matrix_flat.matrix
  external add : Gsl_matrix_flat.matrix -> Gsl_matrix_flat.matrix -> unit
    = "ml_gsl_matrix_add"
  external sub : Gsl_matrix_flat.matrix -> Gsl_matrix_flat.matrix -> unit
    = "ml_gsl_matrix_sub"
  external mul_elements :
    Gsl_matrix_flat.matrix -> Gsl_matrix_flat.matrix -> unit
    = "ml_gsl_matrix_mul"
  external div_elements :
    Gsl_matrix_flat.matrix -> Gsl_matrix_flat.matrix -> unit
    = "ml_gsl_matrix_div"
  external scale : Gsl_matrix_flat.matrix -> float -> unit
    = "ml_gsl_matrix_scale"
  external add_constant : Gsl_matrix_flat.matrix -> float -> unit
    = "ml_gsl_matrix_add_constant"
  external add_diagonal : Gsl_matrix_flat.matrix -> float -> unit
    = "ml_gsl_matrix_add_diagonal"
  external is_null : Gsl_matrix_flat.matrix -> bool = "ml_gsl_matrix_isnull"
  external swap_rows : Gsl_matrix_flat.matrix -> int -> int -> unit
    = "ml_gsl_matrix_swap_rows"
  external swap_columns : Gsl_matrix_flat.matrix -> int -> int -> unit
    = "ml_gsl_matrix_swap_columns"
  external swap_rowcol : Gsl_matrix_flat.matrix -> int -> int -> unit
    = "ml_gsl_matrix_swap_rowcol"
  external transpose :
    Gsl_matrix_flat.matrix -> Gsl_matrix_flat.matrix -> unit
    = "ml_gsl_matrix_transpose_memcpy"
  external transpose_in_place : Gsl_matrix_flat.matrix -> unit
    = "ml_gsl_matrix_transpose"
  val submatrix :
    Gsl_matrix_flat.matrix ->
    k1:int -> k2:int -> n1:int -> n2:int -> Gsl_matrix_flat.matrix
  val row : Gsl_matrix_flat.matrix -> int -> Gsl_vector_flat.vector
  val column : Gsl_matrix_flat.matrix -> int -> Gsl_vector_flat.vector
  val diagonal : Gsl_matrix_flat.matrix -> Gsl_vector_flat.vector
  val subdiagonal : Gsl_matrix_flat.matrix -> int -> Gsl_vector_flat.vector
  val superdiagonal : Gsl_matrix_flat.matrix -> int -> Gsl_vector_flat.vector
  val view_array :
    float array ->
    ?off:int -> int -> ?tda:int -> int -> Gsl_matrix_flat.matrix
  val view_vector :
    Gsl_vector_flat.vector ->
    ?off:int -> int -> ?tda:int -> int -> Gsl_matrix_flat.matrix
end