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

//
// A rust binding for the GSL library by Guillaume Gomez (guillaume1.gomez@gmail.com)
//

use ffi;

/// Compute the covariance matrix cov = inv (J^T J) by QRP^T decomposition of J
pub fn covar(J: &::MatrixF64, epsrel: f64, covar: &mut ::MatrixF64) -> ::Value {
    ::Value::from(unsafe {
        ffi::gsl_multifit_covar(ffi::FFI::unwrap_shared(J), epsrel, ffi::FFI::unwrap_unique(covar))
    })
}

pub fn test_delta(dx: &::VectorF64, x: &::VectorF64, epsabs: f64, epsrel: f64) -> ::Value {
    ::Value::from(unsafe {
        ffi::gsl_multifit_test_delta(ffi::FFI::unwrap_shared(dx),
                                     ffi::FFI::unwrap_shared(x),
                                     epsabs,
                                     epsrel)
    })
}

pub fn gradient(J: &::MatrixF64, f: &::VectorF64, g: &mut ::VectorF64) -> ::Value {
    ::Value::from(unsafe {
        ffi::gsl_multifit_gradient(ffi::FFI::unwrap_shared(J),
                                   ffi::FFI::unwrap_shared(f),
                                   ffi::FFI::unwrap_unique(g))
    })
}