Matrix subtraction correctness tests results

Sample 5x5 matrices:
Matrix 1:
     1000000               4000               6000            80000              10000
           4000         1100000             16000            18000              20000
           6000             16000         1200000            28000              30000
           8000             18000             28000         1300000             40000
         10000             20000             30000             40000         1400000
Matrix 2:
      1000                  4               6                8              10
            4            1100             16              18              20
            6               16          1200              28              30
            8               18              28          1300              40
           10              20              30              40          1400
 

1. Matrix - Matrix, 1 time

LinearAlgebra

Matrix1 - Matrix2 =
  9.99e+05            4e+03        5.99e+03        7.99e+03        9.99e+03
       4e+03         1.1e+06          1.6e+04          1.8e+04             2e+04
  5.99e+03         1.6e+04          1.2e+06          2.8e+04             3e+04
  7.99e+03         1.8e+04          2.8e+04          1.3e+06             4e+04
  9.99e+03            2e+04             3e+04             4e+04          1.4e+06

LinearAlgebra (Symmetric matrices)

Matrix1 - Matrix2 =
  9.99e+05            4e+03        5.99e+03        7.99e+03        9.99e+03
       4e+03         1.1e+06          1.6e+04          1.8e+04             2e+04
  5.99e+03         1.6e+04          1.2e+06          2.8e+04             3e+04
  7.99e+03         1.8e+04          2.8e+04          1.3e+06             4e+04
  9.99e+03            2e+04             3e+04             4e+04          1.4e+06

CLHEP

Matrix1 - Matrix2 =
      999000              3996                5994               7992                9990
         3996    1.0989e+06              15984             17982              19980
         5994             15984     1.1988e+06             27972              29970
         7992             17982              27972    1.2987e+06              39960
         9990             19980              29970              39960    1.3986e+06

CLHEP (Symmetric matrices)

Matrix1 - Matrix2 =
      999000              3996                5994               7992                9990
         3996    1.0989e+06              15984             17982              19980
         5994             15984     1.1988e+06             27972              29970
         7992             17982              27972    1.2987e+06              39960
         9990             19980              29970             39960     1.3986e+06

CovMatrices (inherently symmetric)

Matrix1 - Matrix2 =
  +9.9900000000e+05
  +3.9960000000e+03   +1.0989000000e+06
  +5.9940000000e+03   +1.5984000000e+04   +1.1988000000e+06
  +7.9920000000e+03   +1.7982000000e+04   +2.7972000000e+04   +1.2987000000e+06
  +9.9900000000e+03   +1.9980000000e+04   +2.9970000000e+04   +3.9960000000e+04   +1.3986000000e+06
 

2. Matrix - Matrix, 100 times (the difference of the two matrices is stored on the first matrix)

LinearAlgebra

Matrix1 - Matrix2 =
      9e+05          3.6e+03          5.4e+03          7.2e+03             9e+03
   3.6e+03          9.9e+05        1.44e+04        1.62e+04          1.8e+04
   5.4e+03        1.44e+04        1.08e+06        2.52e+04          2.7e+04
   7.2e+03        1.62e+04        2.52e+04        1.17e+06          3.6e+04
      9e+03          1.8e+04          2.7e+04          3.6e+04        1.26e+06

LinearAlgebra (Symmetric matrices)

Matrix1 - Matrix2 =
      9e+05          3.6e+03          5.4e+03          7.2e+03             9e+03
   3.6e+03          9.9e+05        1.44e+04        1.62e+04          1.8e+04
   5.4e+03        1.44e+04        1.08e+06        2.52e+04          2.7e+04
   7.2e+03        1.62e+04        2.52e+04        1.17e+06          3.6e+04
      9e+03          1.8e+04          2.7e+04          3.6e+04        1.26e+06

CLHEP

Matrix1 - Matrix2 =
       900000            3600             5400             7200             9000
           3600        990000           14400           16200           18000
           5400          14400      1.08e+06           25200           27000
           7200          16200           25200      1.17e+06           36000
           9000          18000           27000           36000      1.26e+06

CLHEP (Symmetric Matrices)

Matrix1 - Matrix2 =
       900000            3600             5400             7200             9000
           3600        990000           14400           16200           18000
           5400          14400      1.08e+06           25200           27000
           7200          16200           25200      1.17e+06           36000
           9000          18000           27000           36000      1.26e+06

CovMatrices (Inherently Symmetric)

Matrix1 - Matrix2 =
  +9.0000000000e+05
  +3.6000000000e+03   +9.9000000000e+05
  +5.4000000000e+03   +1.4400000000e+04   +1.0800000000e+06
  +7.2000000000e+03   +1.6200000000e+04   +2.5200000000e+04   +1.1700000000e+06
  +9.0000000000e+03   +1.8000000000e+04   +2.7000000000e+04   +3.6000000000e+04   +1.2600000000e+06
 

Tests were compiled and build using GCC 3.0.1

Nick Macks