public void ComplexMatrixAddition()
{
/*
* MATLAB:
* sum_cc = ca3x2 + cb3x2
* diff_cc = ca3x2 - cb3x2
* sum_cm = ca3x2 + mb3x2
* diff_cm = ca3x2 - mb3x2
* sum_cs = ca3x2 + s
* diff_cs = ca3x2 - s
* neg_c = -ca3x2
*/
// ComplexMatrix + ComplexMatrix
ComplexMatrix sumCmCm = new ComplexMatrix(new Complex[][] {
new Complex[] { 15 + (72 * j), 1.5 + (16.5 * j) },
new Complex[] { -2 - (37 * j), 9.5 - (43.5 * j) },
new Complex[] { 32 + (137 * j), 11 - (96 * j) }
});
Assert.That(_ca3X2 + _cb3X2, NumericIs.AlmostEqualTo(sumCmCm), "sum cc 1");
ComplexMatrix sumCmCmInplace = _ca3X2.Clone();
Assert.That(sumCmCmInplace.Add(_cb3X2), NumericIs.AlmostEqualTo(sumCmCm), "sum cc 2");
sumCmCmInplace.AddInplace(_cb3X2);
Assert.That(sumCmCmInplace, NumericIs.AlmostEqualTo(sumCmCm), "sum cc 3");
// ComplexMatrix - ComplexMatrix
ComplexMatrix diffCmCm = new ComplexMatrix(new Complex[][] {
new Complex[] { -3 - (96 * j), -1.5 - (16.5 * j) },
new Complex[] { 6 + (29 * j), 14.5 - (4.5 * j) },
new Complex[] { -4 - (193 * j), 25 + (24 * j) }
});
Assert.That(_ca3X2 - _cb3X2, NumericIs.AlmostEqualTo(diffCmCm), "diff cc 1");
ComplexMatrix diffCmCmInplace = _ca3X2.Clone();
Assert.That(diffCmCmInplace.Subtract(_cb3X2), NumericIs.AlmostEqualTo(diffCmCm), "diff cc 2");
diffCmCmInplace.SubtractInplace(_cb3X2);
Assert.That(diffCmCmInplace, NumericIs.AlmostEqualTo(diffCmCm), "diff cc 3");
// ComplexMatrix + Matrix
ComplexMatrix sumCmM = new ComplexMatrix(new Complex[][] {
new Complex[] { 16 - (12 * j), 2.5 },
new Complex[] { -1 - (4 * j), 10.5 - (24 * j) },
new Complex[] { 33 - (28 * j), 12 - (36 * j) }
});
Assert.That(_ca3X2 + _mb3X2, NumericIs.AlmostEqualTo(sumCmM), "sum cm 1");
ComplexMatrix sumCmMInplace = _ca3X2.Clone();
Assert.That(sumCmMInplace.Add(_mb3X2), NumericIs.AlmostEqualTo(sumCmM), "sum cm 2");
sumCmMInplace.AddInplace(_mb3X2);
Assert.That(sumCmMInplace, NumericIs.AlmostEqualTo(sumCmM), "sum cm 3");
// ComplexMatrix - Matrix
ComplexMatrix diffCmM = new ComplexMatrix(new Complex[][] {
new Complex[] { -4 - (12 * j), -2.5 },
new Complex[] { 5 - (4 * j), 13.5 - (24 * j) },
new Complex[] { -5 - (28 * j), 24 - (36 * j) }
});
Assert.That(_ca3X2 - _mb3X2, NumericIs.AlmostEqualTo(diffCmM), "diff cm 1");
ComplexMatrix diffCmMInplace = _ca3X2.Clone();
Assert.That(diffCmMInplace.Subtract(_mb3X2), NumericIs.AlmostEqualTo(diffCmM), "diff cm 2");
diffCmMInplace.SubtractInplace(_mb3X2);
Assert.That(diffCmMInplace, NumericIs.AlmostEqualTo(diffCmM), "diff cm 3");
// ComplexMatrix + Complex
ComplexMatrix sumCmC = new ComplexMatrix(new Complex[][] {
new Complex[] { 7 - (11 * j), 1 + j },
new Complex[] { 3 - (3 * j), 13 - (23 * j) },
new Complex[] { 15 - (27 * j), 19 - (35 * j) }
});
Assert.That(_ca3X2 + _s, NumericIs.AlmostEqualTo(sumCmC), "sum cs 1");
ComplexMatrix sumCmCInplace = _ca3X2.Clone();
Assert.That(sumCmCInplace.Add(_s), NumericIs.AlmostEqualTo(sumCmC), "sum cs 2");
sumCmCInplace.AddInplace(_s);
Assert.That(sumCmCInplace, NumericIs.AlmostEqualTo(sumCmC), "sum cs 3");
// ComplexMatrix - Complex
ComplexMatrix diffCmC = new ComplexMatrix(new Complex[][] {
new Complex[] { 5 - (13 * j), -1 - j },
new Complex[] { 1 - (5 * j), 11 - (25 * j) },
new Complex[] { 13 - (29 * j), 17 - (37 * j) }
});
Assert.That(_ca3X2 - _s, NumericIs.AlmostEqualTo(diffCmC), "diff cs 1");
ComplexMatrix diffCmCInplace = _ca3X2.Clone();
Assert.That(diffCmCInplace.Subtract(_s), NumericIs.AlmostEqualTo(diffCmC), "diff cs 2");
diffCmCInplace.SubtractInplace(_s);
Assert.That(diffCmCInplace, NumericIs.AlmostEqualTo(diffCmC), "diff cs 3");
// ComplexMatrix Negate
ComplexMatrix negateCm = new ComplexMatrix(new Complex[][] {
new Complex[] { -6 + (12 * j), 0 },
new Complex[] { -2 + (4 * j), -12 + (24 * j) },
new Complex[] { -14 + (28 * j), -18 + (36 * j) }
});
Assert.That(-_ca3X2, NumericIs.AlmostEqualTo(negateCm), "neg c 1");
ComplexMatrix negCmInplace = _ca3X2.Clone();
Assert.That(negCmInplace.Negate(), NumericIs.AlmostEqualTo(negateCm), "neg c 2");
negCmInplace.NegateInplace();
Assert.That(negCmInplace, NumericIs.AlmostEqualTo(negateCm), "neg c 3");
}
public void TestComplexMatrix_AdditiveTranspose()
{
/*
* MATLAB:
* sum_cc = ca3x2 + cb3x2
* diff_cc = ca3x2 - cb3x2
* sum_cm = ca3x2 + mb3x2
* diff_cm = ca3x2 - mb3x2
* sum_cs = ca3x2 + s
* diff_cs = ca3x2 - s
* neg_c = -ca3x2
* conj_c = conj(ca3x2)
* trans_c = ca3x2.'
* htrans_c = ca3x2'
* trans_c2 = cc2x2.'
* htrans_c2 = cc2x2'
*/
// ComplexMatrix + ComplexMatrix
ComplexMatrix sum_cc = new ComplexMatrix(new Complex[][] {
new Complex[] { 15 + 72 * j, 1.5 + 16.5 * j },
new Complex[] { -2 - 37 * j, 9.5 - 43.5 * j },
new Complex[] { 32 + 137 * j, 11 - 96 * j }
});
NumericAssert.AreAlmostEqual(sum_cc, ca3x2 + cb3x2, "sum cc 1");
ComplexMatrix sum_cc_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(sum_cc, sum_cc_inplace.Add(cb3x2), "sum cc 2");
sum_cc_inplace.AddInplace(cb3x2);
NumericAssert.AreAlmostEqual(sum_cc, sum_cc_inplace, "sum cc 3");
// ComplexMatrix - ComplexMatrix
ComplexMatrix diff_cc = new ComplexMatrix(new Complex[][] {
new Complex[] { -3 - 96 * j, -1.5 - 16.5 * j },
new Complex[] { 6 + 29 * j, 14.5 - 4.5 * j },
new Complex[] { -4 - 193 * j, 25 + 24 * j }
});
NumericAssert.AreAlmostEqual(diff_cc, ca3x2 - cb3x2, "diff cc 1");
ComplexMatrix diff_cc_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(diff_cc, diff_cc_inplace.Subtract(cb3x2), "diff cc 2");
diff_cc_inplace.SubtractInplace(cb3x2);
NumericAssert.AreAlmostEqual(diff_cc, diff_cc_inplace, "diff cc 3");
// ComplexMatrix + Matrix
ComplexMatrix sum_cm = new ComplexMatrix(new Complex[][] {
new Complex[] { 16 - 12 * j, 2.5 },
new Complex[] { -1 - 4 * j, 10.5 - 24 * j },
new Complex[] { 33 - 28 * j, 12 - 36 * j }
});
NumericAssert.AreAlmostEqual(sum_cm, ca3x2 + mb3x2, "sum cm 1");
ComplexMatrix sum_cm_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(sum_cm, sum_cm_inplace.Add(mb3x2), "sum cm 2");
sum_cm_inplace.AddInplace(mb3x2);
NumericAssert.AreAlmostEqual(sum_cm, sum_cm_inplace, "sum cm 3");
// ComplexMatrix - Matrix
ComplexMatrix diff_cm = new ComplexMatrix(new Complex[][] {
new Complex[] { -4 - 12 * j, -2.5 },
new Complex[] { 5 - 4 * j, 13.5 - 24 * j },
new Complex[] { -5 - 28 * j, 24 - 36 * j }
});
NumericAssert.AreAlmostEqual(diff_cm, ca3x2 - mb3x2, "diff cm 1");
ComplexMatrix diff_cm_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(diff_cm, diff_cm_inplace.Subtract(mb3x2), "diff cm 2");
diff_cm_inplace.SubtractInplace(mb3x2);
NumericAssert.AreAlmostEqual(diff_cm, diff_cm_inplace, "diff cm 3");
// ComplexMatrix + Complex
ComplexMatrix sum_cs = new ComplexMatrix(new Complex[][] {
new Complex[] { 7 - 11 * j, 1 + j },
new Complex[] { 3 - 3 * j, 13 - 23 * j },
new Complex[] { 15 - 27 * j, 19 - 35 * j }
});
NumericAssert.AreAlmostEqual(sum_cs, ca3x2 + s, "sum cs 1");
ComplexMatrix sum_cs_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(sum_cs, sum_cs_inplace.Add(s), "sum cs 2");
sum_cs_inplace.AddInplace(s);
NumericAssert.AreAlmostEqual(sum_cs, sum_cs_inplace, "sum cs 3");
// ComplexMatrix - Complex
ComplexMatrix diff_cs = new ComplexMatrix(new Complex[][] {
new Complex[] { 5 - 13 * j, -1 - j },
new Complex[] { 1 - 5 * j, 11 - 25 * j },
new Complex[] { 13 - 29 * j, 17 - 37 * j }
});
NumericAssert.AreAlmostEqual(diff_cs, ca3x2 - s, "diff cs 1");
ComplexMatrix diff_cs_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(diff_cs, diff_cs_inplace.Subtract(s), "diff cs 2");
diff_cs_inplace.SubtractInplace(s);
NumericAssert.AreAlmostEqual(diff_cs, diff_cs_inplace, "diff cs 3");
// ComplexMatrix Negate
ComplexMatrix neg_c = new ComplexMatrix(new Complex[][] {
new Complex[] { -6 + 12 * j, 0 },
new Complex[] { -2 + 4 * j, -12 + 24 * j },
new Complex[] { -14 + 28 * j, -18 + 36 * j }
});
NumericAssert.AreAlmostEqual(neg_c, -ca3x2, "neg c 1");
ComplexMatrix neg_c_inplace = ca3x2.Clone();
NumericAssert.AreAlmostEqual(neg_c, neg_c_inplace.Negate(), "neg c 2");
neg_c_inplace.NegateInplace();
NumericAssert.AreAlmostEqual(neg_c, neg_c_inplace, "neg c 3");
// ComplexMatrix Conjugate
ComplexMatrix conj_c = new ComplexMatrix(new Complex[][] {
new Complex[] { 6 + 12 * j, 0 },
new Complex[] { 2 + 4 * j, 12 + 24 * j },
new Complex[] { 14 + 28 * j, 18 + 36 * j }
});
NumericAssert.AreAlmostEqual(conj_c, ca3x2.Conjugate(), "conj c 1");
ComplexMatrix conj_c_inplace = ca3x2.Clone();
conj_c_inplace.ConjugateInplace();
NumericAssert.AreAlmostEqual(conj_c, conj_c_inplace, "conj c 2");
// ComplexMatrix Transpose (Non-Conjugated)
ComplexMatrix trans_c = new ComplexMatrix(new Complex[][] {
new Complex[] { 6 - 12 * j, 2 - 4 * j, 14 - 28 * j },
new Complex[] { 0, 12 - 24 * j, 18 - 36 * j }
});
NumericAssert.AreAlmostEqual(trans_c, ca3x2.Transpose(), "trans c 1");
// ComplexMatrix Hermitian Transpose (Conjugated)
ComplexMatrix htrans_c = new ComplexMatrix(new Complex[][] {
new Complex[] { 6 + 12 * j, 2 + 4 * j, 14 + 28 * j },
new Complex[] { 0, 12 + 24 * j, 18 + 36 * j }
});
NumericAssert.AreAlmostEqual(htrans_c, ca3x2.HermitianTranspose(), "htrans c 1");
// ComplexMatrix Transpose (Non-Conjugated) (Square)
ComplexMatrix trans_c2 = new ComplexMatrix(new Complex[][] {
new Complex[] { 8 - 24 * j, 10 - 30 * j },
new Complex[] { 9 - 27 * j, 11 - 33 * j }
});
NumericAssert.AreAlmostEqual(trans_c2, cc2x2.Transpose(), "trans c2 1");
ComplexMatrix trans_c2_inplace = cc2x2.Clone();
trans_c2_inplace.TransposeInplace();
NumericAssert.AreAlmostEqual(trans_c2, trans_c2_inplace, "trans c2 2");
// ComplexMatrix Hermitian Transpose (Conjugated) (Square)
ComplexMatrix htrans_c2 = new ComplexMatrix(new Complex[][] {
new Complex[] { 8 + 24 * j, 10 + 30 * j },
new Complex[] { 9 + 27 * j, 11 + 33 * j }
});
NumericAssert.AreAlmostEqual(htrans_c2, cc2x2.HermitianTranspose(), "htrans c2 1");
ComplexMatrix htrans_c2_inplace = cc2x2.Clone();
htrans_c2_inplace.HermitianTransposeInplace();
NumericAssert.AreAlmostEqual(htrans_c2, htrans_c2_inplace, "htrans c2 2");
}