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
open Fun
open Vector
module NoDeriv =
struct
type kind =
| HYBRIDS
| HYBRID
| DNEWTON
| BROYDEN
type t
external _alloc : kind -> int -> t
= "ml_gsl_multiroot_fsolver_alloc"
external _free : t -> unit
= "ml_gsl_multiroot_fsolver_free"
external _set : t -> multi_fun -> vector -> unit
= "ml_gsl_multiroot_fsolver_set"
let make kind dim f x =
let s = _alloc kind dim in
Gc.finalise _free s ;
_set s f x ;
s
external name : t -> string
= "ml_gsl_multiroot_fsolver_name"
external iterate : t -> unit
= "ml_gsl_multiroot_fsolver_iterate"
external root : t -> vector -> unit
= "ml_gsl_multiroot_fsolver_root"
external get_state : t ->
?x:vector -> ?f:vector ->
?dx:vector -> unit -> unit
= "ml_gsl_multiroot_fsolver_get_state"
external test_delta : t -> epsabs:float -> epsrel:float -> bool
= "ml_gsl_multiroot_test_delta_f"
external test_residual : t -> epsabs:float -> bool
= "ml_gsl_multiroot_test_residual_f"
end
module Deriv =
struct
type kind =
| HYBRIDSJ
| HYBRIDJ
| NEWTON
| GNEWTON
type t
external _alloc : kind -> int -> t
= "ml_gsl_multiroot_fdfsolver_alloc"
external _free : t -> unit
= "ml_gsl_multiroot_fdfsolver_free"
external _set : t -> multi_fun_fdf -> vector -> unit
= "ml_gsl_multiroot_fdfsolver_set"
let make kind dim f x =
let s = _alloc kind dim in
Gc.finalise _free s ;
_set s f x ;
s
external name : t -> string
= "ml_gsl_multiroot_fdfsolver_name"
external root : t -> vector -> unit
= "ml_gsl_multiroot_fdfsolver_root"
external iterate : t -> unit
= "ml_gsl_multiroot_fdfsolver_iterate"
external get_state : t ->
?x:vector -> ?f:vector ->
?j:Matrix.matrix -> ?dx:vector -> unit -> unit
= "ml_gsl_multiroot_fdfsolver_get_state_bc" "ml_gsl_multiroot_fdfsolver_get_state"
external test_delta : t -> epsabs:float -> epsrel:float -> bool
= "ml_gsl_multiroot_test_delta_fdf"
external test_residual : t -> epsabs:float -> bool
= "ml_gsl_multiroot_test_residual_fdf"
end