我正在計算使用Cholesky分解使用不同庫的方陣的倒數。然而我的結果並不像我期待的那樣。我不是數學專家,但我期待得到更接近的結果。矩陣求逆Cholesky分解 - >結果不準確
我使用MLK,岩漿和CULA庫來計算CPU和GPU中矩陣的逆矩陣。 經過計算後,我發現這些庫的結果總是不同的一個元素。說我要來計算A= [0.237306,0.000458;0.000458,0.238497]
倒數:
A[0] = 0.237306
A[1] = 0.000458
A[2] = 0.000458
A[3] = 0.238497
我得到的結果是:
inv(A)[0] = 4.213983
inv(A)[1] = -0.008092
inv(A)[2] = 0.000458
inv(A)[3] = 4.192946
但是,正確的結果應該是
4.2139841 -0.0080924
-0.0080924 4.1929404
正如你所看到的,inv(A)[3]
是不同的,儘管其餘的都很好。那麼Cholesky Inversion是如何工作的?這是一個正確的/近似的結果,還是我在這裏做錯了什麼?
要麼提供使用庫的信息,要麼嘗試http://math.stackexchange.com/ – ipavlic 2012-07-20 12:02:14
您確定這些值都應該可用嗎?典型的對稱矩陣庫只用作矩陣的一邊輸出,依靠你通過對稱完成另一邊。請張貼導致結果的代碼大綱。 – 2012-07-20 12:07:44