pub fn cdotu_sub<T>(
    N: i32,
    x: &[T],
    incx: i32,
    y: &[T],
    incy: i32,
    dotu: &mut [T]
)