1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119
|
(* gsl-ocaml - OCaml interface to GSL *)
(* Copyright (©) 2002-2012 - Olivier Andrieu *)
(* Distributed under the terms of the GPL version 3 *)
let () = Error.init ()
(** Histograms *)
(** The histogram type *)
type t = {
n : int ;
range : float array ;
bin : float array ;
}
let check h =
h.n > 0 && Array.length h.range = succ h.n && Array.length h.bin = h.n
(** {3 Allocating histograms} *)
let make n =
{ n = n ; range = Array.make (succ n) 0. ; bin = Array.make n 0. ; }
let copy h =
{ n = h.n ;
range = Array.copy h.range ;
bin = Array.copy h.bin ;
}
external set_ranges : t -> float array -> unit
= "ml_gsl_histogram_set_ranges"
external set_ranges_uniform : t -> xmin:float -> xmax:float -> unit
= "ml_gsl_histogram_set_ranges_uniform"
(** {3 Updating and accessing histogram elements} *)
external accumulate : t -> ?w:float -> float -> unit
= "ml_gsl_histogram_accumulate"
let get h i =
h.bin.(i)
let get_range h i =
(h.range.(i), h.range.(succ i))
let h_max h =
h.range.(h.n)
let h_min h =
h.range.(0)
let bins h =
h.n
let reset h =
Array.fill h.bin 0 h.n 0.
(** {3 Searching histogram ranges} *)
external find : t -> float -> int
= "ml_gsl_histogram_find"
(** {3 Histograms statistics } *)
external max_val : t -> float
= "ml_gsl_histogram_max_val"
external max_bin : t -> int
= "ml_gsl_histogram_max_bin"
external min_val : t -> float
= "ml_gsl_histogram_min_val"
external min_bin : t -> int
= "ml_gsl_histogram_min_bin"
external mean : t -> float
= "ml_gsl_histogram_mean"
external sigma : t -> float
= "ml_gsl_histogram_sigma"
external sum : t -> float
= "ml_gsl_histogram_sum"
(** {3 Histogram operations} *)
external equal_bins_p : t -> t -> bool
= "ml_gsl_histogram_equal_bins_p"
external add : t -> t -> unit
= "ml_gsl_histogram_add"
external sub : t -> t -> unit
= "ml_gsl_histogram_sub"
external mul : t -> t -> unit
= "ml_gsl_histogram_mul"
external div : t -> t -> unit
= "ml_gsl_histogram_div"
external scale : t -> float -> unit
= "ml_gsl_histogram_scale"
external shift : t -> float -> unit
= "ml_gsl_histogram_shift"
(** {3 Resampling} *)
type histo_pdf = {
pdf_n : int ;
pdf_range : float array ;
pdf_sum : float array ;
}
external _init : histo_pdf -> t -> unit
= "ml_gsl_histogram_pdf_init"
let init h =
let p = { pdf_n = h.n ;
pdf_range = Array.copy h.range ;
pdf_sum = Array.copy h.bin ; } in
_init p h ;
p
external sample : histo_pdf -> float -> float
= "ml_gsl_histogram_pdf_sample"
|