/// <summary>
/// Multiply for a symmetric matrix S. c = S * b
/// </summary>
/// <param name="A"></param>
/// <param name="b"></param>
/// <param name="c"></param>
public static void Symmetric_Multiply_Update(matrix S, vector b, vector c)
{
Blas.symm(Layout.ColMajor, Side.Left, UpLo.Lower, b.Length, 1, 1, S.Array, S.Rows, b.Array, b.Length, 0, c.Array, c.Length);
}