Function matrixmultiply::dgemm

pub unsafe fn dgemm(m: usize,
                    k: usize,
                    n: usize,
                    alpha: f64,
                    a: *const f64,
                    rsa: isize,
                    csa: isize,
                    b: *const f64,
                    rsb: isize,
                    csb: isize,
                    beta: f64,
                    c: *mut f64,
                    rsc: isize,
                    csc: isize)

General matrix multiplication (f64)

C ← α A B + β C

• m, k, n: dimensions
• a, b, c: pointer to the first element in the matrix
• A: m by k matrix
• B: k by n matrix
• C: m by n matrix
• rsx: row stride of x
• csx: col stride of x

Strides for A and B may be arbitrary. Strides for C must not result in elements that alias each other, for example they can not be zero.

If β is zero, then C does not need to be initialized.