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
46
47
48
49
50
51

/// Scalar vector multiplication
#[no_mangle]
#[inline(always)]
pub fn dscal(dx: &mut [f64], n: usize, incx: usize, da: &f64) {
    /*
    Purpose : scalar vector multiplication

    dx = da * dx


    --- Input ---

    n    : number of elements in input vector
    da   : double scale factor
    dx   : double vector with n+1 elements, dx[0] is not used
    incx : storage spacing between elements of dx


    --- Output ---

    dx = da * dx, unchanged if n <= 0


    For i = 0 to n-1, replace dx[1+i*incx] with
    da * dx[1+i*incx].

    */
    /* Code for increments not equal to 1.  */
    if incx != 0 {
        for i in (1..n * incx + 1).step_by(incx) {
            dx[i] *= da;
        }
        return;
    }
    /* Code for increments not equal to 1.  */
    let m = n % 4;
    if m != 0 {
        for i in 1..m + 1 {
            dx[i] *= da;
        }
        if n < 4 {
            return;
        }
    }
    for i in (m + 1..n + 1).step_by(4) {
        dx[i] *= da;
        dx[i + 1] *= da;
        dx[i + 2] *= da;
        dx[i + 3] *= da;
    }
}