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
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