File: permut.ml

package info (click to toggle)
ocamlgsl 0.3.5-3
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 3,444 kB
  • ctags: 2,901
  • sloc: ml: 7,956; ansic: 6,796; makefile: 303; sh: 87; awk: 13
file content (31 lines) | stat: -rw-r--r-- 740 bytes parent folder | download | duplicates (4)
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
open Bigarray

let _ = 
  Gsl_error.init () ;
  Random.self_init ()

let print_arr arr = 
  Array.iter (fun i -> Printf.printf "% 4d " i) arr ;
  print_newline ()

let print_barr arr = 
  for i=0 to pred (Array1.dim arr) do
    Printf.printf "% 4d " arr.{i}
  done ;
  print_newline ()

let _ = 
  let p = Gsl_permut.make 5 in
  Gsl_permut.next p ;
  print_string "permut :" ; print_arr (Gsl_permut.to_array p) ;
  let a = Array.init 5 (fun _ -> Random.int 10) in
  print_string "arr :" ; print_arr a ;
  Gsl_permut.permute p a ;
  print_string "arr :" ; print_arr a ;
  
  let a1 = Array1.of_array int c_layout a in
  print_string "arr1 :" ; print_barr a1 ;
  Gsl_permut.permute_barr p a1 ;
  print_string "arr1 :" ; print_barr a1