Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- --------------------TEST HILBERT---------------------------
- For Hilbert Matrix size n = 1:
- Rank = 1, cond = 1.0
- Gauss error: 0.0
- Inverse error: 0.0
- For Hilbert Matrix size n = 2:
- Rank = 2, cond = 19.281470067903967
- Gauss error: 5.661048867003676e-16
- Inverse error: 1.1240151438116956e-15
- For Hilbert Matrix size n = 3:
- Rank = 3, cond = 524.056777586062
- Gauss error: 8.118051169482656e-15
- Inverse error: 1.7907430486334138e-14
- For Hilbert Matrix size n = 4:
- Rank = 4, cond = 15513.738738929262
- Gauss error: 3.349632515431573e-13
- Inverse error: 2.268815452633455e-13
- For Hilbert Matrix size n = 5:
- Rank = 5, cond = 476607.2502421033
- Gauss error: 2.8186181571329407e-13
- Inverse error: 3.654697861370179e-12
- For Hilbert Matrix size n = 6:
- Rank = 6, cond = 1.495105864172721e7
- Gauss error: 2.344229118278564e-10
- Inverse error: 2.1526924552188031e-10
- For Hilbert Matrix size n = 7:
- Rank = 7, cond = 4.753673559839011e8
- Gauss error: 7.410208490784287e-9
- Inverse error: 1.1004084486340423e-8
- For Hilbert Matrix size n = 8:
- Rank = 8, cond = 1.525757550554701e10
- Gauss error: 3.736355234341925e-7
- Inverse error: 3.044718586687223e-7
- For Hilbert Matrix size n = 9:
- Rank = 9, cond = 4.931538348163301e11
- Gauss error: 1.0485260733621061e-5
- Inverse error: 5.958958489267314e-6
- For Hilbert Matrix size n = 10:
- Rank = 10, cond = 1.6025337742793652e13
- Gauss error: 0.0001902892061817039
- Inverse error: 0.00023929944568637062
- For Hilbert Matrix size n = 11:
- Rank = 10, cond = 5.219813567997335e14
- Gauss error: 0.004840612082131825
- Inverse error: 0.006720407723584139
- For Hilbert Matrix size n = 12:
- Rank = 11, cond = 1.6546640506383568e16
- Gauss error: 0.0515312199356675
- Inverse error: 0.06988898126789012
- For Hilbert Matrix size n = 13:
- Rank = 11, cond = 2.376786926717726e18
- Gauss error: 0.8773261344142039
- Inverse error: 1.5119571944035046
- For Hilbert Matrix size n = 14:
- Rank = 11, cond = 2.44167173619644e17
- Gauss error: 3.934776094572252
- Inverse error: 5.840670521642857
- For Hilbert Matrix size n = 15:
- Rank = 12, cond = 2.4022870188034854e17
- Gauss error: 35.372359774035836
- Inverse error: 31.13950131676681
- For Hilbert Matrix size n = 16:
- Rank = 12, cond = 7.301849062715712e17
- Gauss error: 5.345783450414308
- Inverse error: 4.99917180206163
- For Hilbert Matrix size n = 17:
- Rank = 12, cond = 4.7507771590947725e17
- Gauss error: 14.381303216697646
- Inverse error: 15.299960220214384
- For Hilbert Matrix size n = 18:
- Rank = 12, cond = 1.0621139181205069e18
- Gauss error: 5.454180864329791
- Inverse error: 6.942728611787584
- For Hilbert Matrix size n = 19:
- Rank = 13, cond = 4.215819618491786e18
- Gauss error: 7.007414463759739
- Inverse error: 16.272559569359128
- For Hilbert Matrix size n = 20:
- Rank = 13, cond = 2.365121119732083e18
- Gauss error: 11.373002753850809
- Inverse error: 16.695512931170512
- --------------------TEST RANDOM----------------------------
- For n = 5:
- +For c = 1.0: For Random Matrix size n = 5:
- Rank = 5, cond = 1.0000000000000007
- Gauss error: 1.9860273225978183e-16
- Inverse error: 1.4043333874306804e-16
- +For c = 10.0: For Random Matrix size n = 5:
- Rank = 5, cond = 9.999999999999984
- Gauss error: 8.067281120720504e-16
- Inverse error: 7.497088734768789e-16
- +For c = 1000.0: For Random Matrix size n = 5:
- Rank = 5, cond = 999.9999999999893
- Gauss error: 1.4569578196833943e-14
- Inverse error: 1.5409959571464003e-14
- +For c = 1.0e7: For Random Matrix size n = 5:
- Rank = 5, cond = 9.999999995016571e6
- Gauss error: 6.509002890638254e-11
- Inverse error: 7.16802786267915e-11
- +For c = 1.0e12: For Random Matrix size n = 5:
- Rank = 5, cond = 1.0000667836867411e12
- Gauss error: 1.5754201967782774e-5
- Inverse error: 1.1473277192184816e-5
- +For c = 1.0e16: For Random Matrix size n = 5:
- Rank = 4, cond = 6.684139511958563e15
- Gauss error: 0.2889029295986325
- Inverse error: 0.36543586339438777
- For n = 10:
- +For c = 1.0: For Random Matrix size n = 10:
- Rank = 10, cond = 1.0000000000000007
- Gauss error: 3.293453726225543e-16
- Inverse error: 1.7554167342883504e-16
- +For c = 10.0: For Random Matrix size n = 10:
- Rank = 10, cond = 9.999999999999988
- Gauss error: 2.8522145930998397e-16
- Inverse error: 3.439900227959406e-16
- +For c = 1000.0: For Random Matrix size n = 10:
- Rank = 10, cond = 1000.0000000000089
- Gauss error: 2.281882411185135e-14
- Inverse error: 1.9081016541524968e-14
- +For c = 1.0e7: For Random Matrix size n = 10:
- Rank = 10, cond = 9.999999992953647e6
- Gauss error: 1.0675497117591031e-10
- Inverse error: 9.083790051135667e-11
- +For c = 1.0e12: For Random Matrix size n = 10:
- Rank = 10, cond = 1.0000160016557174e12
- Gauss error: 9.796100712610881e-6
- Inverse error: 1.1734966723427664e-5
- +For c = 1.0e16: For Random Matrix size n = 10:
- Rank = 9, cond = 1.0678961887500974e16
- Gauss error: 0.30751924637911515
- Inverse error: 0.2721239637520368
- For n = 20:
- +For c = 1.0: For Random Matrix size n = 20:
- Rank = 20, cond = 1.0000000000000007
- Gauss error: 4.2494139782533216e-16
- Inverse error: 6.392224159095141e-16
- +For c = 10.0: For Random Matrix size n = 20:
- Rank = 20, cond = 10.000000000000004
- Gauss error: 5.874748045952207e-16
- Inverse error: 5.800848902387804e-16
- +For c = 1000.0: For Random Matrix size n = 20:
- Rank = 20, cond = 999.9999999999872
- Gauss error: 7.431365537108807e-15
- Inverse error: 3.874582920122523e-15
- +For c = 1.0e7: For Random Matrix size n = 20:
- Rank = 20, cond = 1.0000000000664406e7
- Gauss error: 1.1071229293974872e-10
- Inverse error: 1.5049099906006415e-10
- +For c = 1.0e12: For Random Matrix size n = 20:
- Rank = 20, cond = 9.999686665866554e11
- Gauss error: 1.3652012009378222e-5
- Inverse error: 1.7112407435987304e-5
- +For c = 1.0e16: For Random Matrix size n = 20:
- Rank = 19, cond = 1.9729393524098988e16
- Gauss error: 0.02766854560080551
- Inverse error: 0.034307645676184245
- Process finished with exit code 0
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement