public void SetIdentity()
{
MatrixD m = new MatrixD(3, 3, 12);
m.SetIdentity();
for(int i=0; i<m.NumberOfRows; i++)
{
for (int j = 0; j < m.NumberOfColumns; j++)
{
if (i == j)
Assert.AreEqual(1, m[i, j]);
else
Assert.AreEqual(0, m[i, j]);
}
}
m = new MatrixD(10, 4, 12);
m.SetIdentity();
for (int i = 0; i < m.NumberOfRows; i++)
{
for (int j = 0; j < m.NumberOfColumns; j++)
{
if (i == j)
Assert.AreEqual(1, m[i, j]);
else
Assert.AreEqual(0, m[i, j]);
}
}
m = new MatrixD(2, 5, 12);
m.SetIdentity();
for (int i = 0; i < m.NumberOfRows; i++)
{
for (int j = 0; j < m.NumberOfColumns; j++)
{
if (i == j)
Assert.AreEqual(1, m[i, j]);
else
Assert.AreEqual(0, m[i, j]);
}
}
}