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.