public static void axpy(double a, BlockVector <double> X, BlockVector <double> Y, ref BlockVector <double> Z) { Y.CopyTo(ref Z); CBLAS.cblas_daxpy(length(X, Y), a, ref head(X), 1, ref head(Z), 1); }