} // private constructor for singleton pattern
#region BLAS Level 1
/// <summary>
/// See http://www.dotnumerics.com/NumericalLibraries/LinearAlgebra/CSharpCodeFiles/daxpy.aspx
/// </summary>
public void Daxpy(int n, double alpha, double[] x, int offsetX, int incX, double[] y, int offsetY, int incY)
=> daxpy.Run(n, alpha, x, offsetX, incX, ref y, offsetY, incY);
} // private constructor for singleton pattern
public void Daxpby(int n, double alpha, double[] x, int offsetX, int incX, double beta, double[] y, int offsetY, int incY)
{
dscal.Run(n, beta, ref y, offsetY, incY);
daxpy.Run(n, alpha, x, offsetX, incX, ref y, offsetY, incY);
}