File: histo.ml

package info (click to toggle)
ocamlgsl 1.24.3-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm, bullseye
  • size: 1,244 kB
  • sloc: ml: 8,748; ansic: 7,395; makefile: 37
file content (119 lines) | stat: -rw-r--r-- 2,671 bytes parent folder | download | duplicates (2)
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"