An experimental demo project of mine in numerical linear algebra and the "Implicit QR Method" in particular.
./qr
ndim = 5
ntimes = 1
*************************************************************************************************************************
-----------------------------------------------H-----------------------------------------------
4.668401669892632 1.102789907490572 0.000000000000000 -0.000000000000000 -0.000000000000000
1.102789907490572 1.331996569198935 0.030895825080788 -0.000000000000000 -0.000000000000000
0.000000000000000 0.030895825080787 2.999601760908436 -0.000000000000000 -0.000000000000000
0.000000000000000 0.000000000000000 0.000000000000000 2.002541127257135 0.071244629176677
0.000000000000000 0.000000000000000 0.000000000000000 0.071244629176677 3.997458872742871
perf: 3 3 7.2000000000E-05
---------------------------H---------------------------
4.999999998394087 0.000080147732139 0.000000000000000
0.000080147732139 1.000000001605916 0.000000000000001
0.000000000000000 0.000000000000000 3.000000000000003
Error tolerance = 1.000000000000000E-13
Maximum relative error = 2.664535259100376E-15
Maximum absolute error = 2.664535259100376E-15
--------Λ------- ------Λref------ ---Rel. error--- ---Abs. error---
1.000000000000002 1.000000000000000 1.5543122345E-15 1.5543122345E-15
2.000000000000005 2.000000000000000 2.6645352591E-15 5.3290705182E-15
3.000000000000003 3.000000000000000 8.8817841970E-16 2.6645352591E-15
4.000000000000002 4.000000000000000 4.4408920985E-16 1.7763568394E-15
5.000000000000002 5.000000000000000 4.4408920985E-16 1.7763568394E-15