Skip to content

keksipurkki/keksipurkki-demos-qr

Repository files navigation

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

About

Implicit QR method (Fortran 2018)

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages