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.