pub unsafe fn scalar_mul_add(a: f64, b: f64, c: f64) -> f64