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
let () = Error.init ()
type ws
external _alloc : k:int -> nbreak:int -> ws = "ml_gsl_bspline_alloc"
external _free : ws -> unit = "ml_gsl_bspline_free"
let make ~k ~nbreak =
let ws = _alloc ~k ~nbreak in
Gc.finalise _free ws;
ws
external ncoeffs : ws -> int = "ml_gsl_bspline_ncoeffs" [@@noalloc]
open Vectmat
external knots : [< vec ] -> ws -> unit = "ml_gsl_bspline_knots"
external knots_uniform : a:float -> b:float -> ws -> unit
= "ml_gsl_bspline_knots_uniform"
external _eval : float -> [< vec ] -> ws -> unit = "ml_gsl_bspline_eval"
let eval ws x =
let n = ncoeffs ws in
let v = `V (Vector.create n) in
_eval x v ws;
v