rusoda 0.1.0

Pure rust and just need log dependencies implementation of odepack dlsoda.
Documentation 
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

/// Operate muladd of 2 vectors and a scalar(the first element of vector will not used).
/// # Examples
///
/// '''
/// let n = 2;
/// let da = 2.0;
/// let dx = vec![0.0, 1.0, 2.0];
/// let dy = vec![0.0, 3.0, 4.0];
/// let axpy = ddot(n, &dx, &mut dy, 0, 0, &da);
/// assert_eq!(vec![5.0, 8.0], axpy);
#[no_mangle]
pub fn daxpy(n: usize, dx: &[f64], dy: &mut [f64], incx: usize, incy: usize, da: &f64) {
    let m: usize;
    let mut ix: usize;
    let mut iy: usize;

    if *da == 0.0 {
        return;
    }
    if incx == 0 && incy == 0 {
        m = n % 4;
        if m != 0 {
            for i in 1..m + 1 {
                dy[i] += da * dx[i];
            }
        }
        if n < 4 {
            return;
        }
        for i in (m + 1..n + 1).step_by(4) {
            dy[i] += da * dx[i];
            dy[i + 1] += da * dx[i + 1];
            dy[i + 2] += da * dx[i + 2];
            dy[i + 3] += da * dx[i + 3];
        }
    } else {
        ix = 1;
        iy = 1;
        for _ in 1..n + 1 {
            dy[iy] += da * dx[ix];
            ix += incx;
            iy += incy;
        }
    }
}