static void Main(string[] args)
{
IMatrix A = new SparseMatrix(2, 3);
IMatrix B = new SparseMatrix(3, 2);
IMatrix C = new SparseMatrix(2, 2);
A[0, 2] = 1;
A[1, 0] = 1;
B[1, 1] = 4;
B[2, 0] = 1;
C[0, 1] = 1;
C[1, 0] = 1;
IMatrix D = (SparseMatrix)A * (SparseMatrix)B + (SparseMatrix)C;
IMatrix E = A.Multiply(B).Add(C);
Console.WriteLine("D:\n{0}", D);
Console.WriteLine("E:\n{0}", E);
}
/// <summary>
/// 矩阵转置
/// </summary>
/// <returns>转置后的新矩阵</returns>
public IMatrix Transpose()
{
SparseMatrix temp = new SparseMatrix(cols, rows);
for (int i = 0; i < temp.Rows; i++)
for (int j = 0; j < temp.Cols; j++)
temp[i, j] = this[j, i];
return temp;
}
/// <summary>
/// 矩阵乘法
/// </summary>
/// <param name="B">与之右乘的另一个矩阵</param>
/// <returns>相乘后的新矩阵</returns>
public IMatrix Multiply(IMatrix B)
{
if (B == null)
throw new Exception("参数矩阵为null.");
if (cols != B.Rows)
throw new Exception("两矩阵不能相乘.");
SparseMatrix temp = new SparseMatrix(rows, B.Cols);
for (int i = 0; i < rows; i++)
{
for (int j = 0; j < B.Cols; j++)
{
for (int k = 0; k < cols; k++)
temp[i, j] += this[i, k] * B[k, j];
}
}
return temp;
}
/// <summary>
/// 矩阵加法
/// </summary>
/// <param name="B">与之相加的另一个矩阵</param>
/// <returns>相加后的新矩阵</returns>
public IMatrix Add(IMatrix B)
{
if (B == null)
throw new Exception("参数矩阵为null.");
if (B.Rows != this.rows || B.Cols != this.cols)
throw new Exception("两矩阵不能相加.");
SparseMatrix temp = new SparseMatrix(rows, cols);
for (int i = 0; i < rows; i++)
for (int j = 0; j < cols; j++)
temp[i, j] = this[i, j] + B[i, j];
return temp;
}