lapacke
pub unsafe fn dlartgp_work( f: f64, g: f64, cs: &mut [f64], sn: &mut [f64], r: &mut [f64]) -> i32