0

Cholesky decomposition with OpenCV

c++
cholesky

asked 2013-10-10 09:57:38 -0600

imikbox
1 ●1

updated 2013-10-10 11:31:39 -0600

Moster
1716 ●5 ●24

Hi there,

I have a prolem using the Cholesky decomposition in OpenCV. I have the positive definite matrix A with

A =

 8312069  34511.344  -306.38861

 34511.313  8231020.5  -112.6131

 -306.38861  -112.6131  1

The function call Mat L = A.inv(DECOMP_CHOLESKY) should give me the Matrix L, which satisfies A = LL'. However, I'm getting a matrix with awkwardly small values and for which A = LL' does not hold. What could there possibly be wrong, the usage is dead simple ... ?? (The decomposition of A works fine in MATLAB).

edit retag flag offensive close merge delete

1 answer

Sort by » oldest newest most voted

2

answered 2013-10-11 01:41:05 -0600

flag of Australia

1000 ●5 ●16

updated 2013-10-11 01:41:58 -0600

L is the inverse of A (using Cholesky Decomposition), not the actual Cholesky Decomposition. To confirm this do:

cout << L * A << endl; // this is inv(A) * A.

You'll get the identity matrix (or very close to it).

edit flag offensive delete link

Comments

this makes sense. I totaly misunderstood the description .. thanks for pointing this out.

imikbox ( 2013-10-14 07:19:32 -0600 )edit

Links

Question Tools

1 follower

Stats

Asked: 2013-10-10 09:57:38 -0600

Seen: 1,338 times

Last updated: Oct 11 '13

Related questions

Copyright OpenCV foundation, 2012-2018. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.

Powered by Askbot version 0.10.2

here is how