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.