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
|
/* ocamlgsl - OCaml interface to GSL */
/* Copyright (©) 2002-2005 - Olivier Andrieu */
/* distributed under the terms of the GPL version 2 */
#include <gsl/gsl_diff.h>
#include <caml/mlvalues.h>
#include <caml/memory.h>
#include "wrappers.h"
#include "mlgsl_fun.h"
value ml_gsl_diff_central(value f, value x)
{
CAMLparam1(f);
double result,abserr;
GSLFUN_CLOSURE(gf, f);
gsl_diff_central(&gf, Double_val(x),
&result, &abserr);
CAMLreturn(copy_two_double_arr(result, abserr));
}
value ml_gsl_diff_forward(value f, value x)
{
CAMLparam1(f);
double result,abserr;
GSLFUN_CLOSURE(gf, f);
gsl_diff_forward(&gf, Double_val(x),
&result, &abserr);
CAMLreturn(copy_two_double_arr(result, abserr));
}
value ml_gsl_diff_backward(value f, value x)
{
CAMLparam1(f);
double result,abserr;
GSLFUN_CLOSURE(gf, f);
gsl_diff_backward(&gf, Double_val(x),
&result, &abserr);
CAMLreturn(copy_two_double_arr(result, abserr));
}
|
