/// <summary>
/// Computes the matrix-vector product y = A*x;
/// </summary>
/// <param name="A">A source matrix of dimension MxN</param>
/// <param name="x">A source vector of length N</param>
/// <param name="y">A target vector of length M</param>
/// <typeparam name="M">The row dimension type of A</typeparam>
/// <typeparam name="N">The column dimension type of A</typeparam>
public static ref Block256 <M, double> gemv <M, N>(Matrix256 <M, N, double> A, Block256 <N, double> x, ref Block256 <M, double> y)
where M : unmanaged, ITypeNat
where N : unmanaged, ITypeNat
{
var m = nat32i <M>();
var n = nat32i <N>();
var lda = n;
CBLAS.cblas_dgemv(RowMajor, NoTranspose, m, n, alpha: 1.0, ref head(A), lda, ref head(x), incX: 1, beta: 0, ref head(y), incY: 1);
return(ref y);
}
/// <summary>
/// Computes the matrix-vector product y = A*x;
/// </summary>
/// <param name="A">A source matrix of dimension MxN</param>
/// <param name="x">A source vector of length N</param>
/// <param name="y">A target vector of length M</param>
/// <typeparam name="M">The row dimension type of A</typeparam>
/// <typeparam name="N">The column dimension type of A</typeparam>
public static ref BlockVector <M, double> gemv <M, N>(BlockMatrix <M, N, double> A, BlockVector <N, double> x, ref BlockVector <M, double> y)
where M : ITypeNat, new()
where N : ITypeNat, new()
{
var m = nati <M>();
var n = nati <N>();
var lda = n;
CBLAS.cblas_dgemv(RowMajor, NoTranspose, m, n, alpha: 1.0, ref head(A), lda, ref head(x), incX: 1, beta: 0, ref head(y), incY: 1);
return(ref y);
}
