pub fn gemm(alpha: f64, a: &Matrix, b: &Matrix, beta: f64, c: &mut Matrix)
Expand description
GEMM wrapper for Matrixmultiply
§Examples
#[macro_use]
extern crate peroxide;
use peroxide::prelude::*;
fn main() {
let a = ml_matrix("1 2 3;4 5 6");
let b = ml_matrix("1 2;3 4;5 6");
let mut c1 = zeros(2, 2);
let mut c2 = matrix(vec![1f64; 9], 3, 3, Col);
gemm(1f64, &a, &b, 0f64, &mut c1);
gemm(1f64, &b, &a, 2f64, &mut c2);
assert_eq!(c1, ml_matrix("22 28; 49 64"));
assert_eq!(c2, ml_matrix("11 14 17;21 28 35;31 42 53"));
}