forum.alglib.net

ALGLIB forum
It is currently Thu Mar 28, 2024 8:15 pm

All times are UTC


Forum rules


1. This forum can be used for discussion of both ALGLIB-related and general numerical analysis questions
2. This forum is English-only - postings in other languages will be removed.



Post new topic Reply to topic  [ 1 post ] 
Author Message
 Post subject: Inverse of matrix square root
PostPosted: Sun Mar 09, 2014 2:36 am 
Offline

Joined: Sun Mar 09, 2014 2:21 am
Posts: 1
Hi,
I was wondering which routine to use to get the fastest results with the least amount of computation...

I have a symmetric 4x4 matrix, A, but there are two cases.
(1) The matrix, A, is positive semi-definite i.e. may have a eigenvalue of 0.
(2) The matrix, A, is positive definite.

In either case I know ahead of time which type it is, so I can use the optimal routine when needed.

I'm looking to solve for the matrix B, where B is defined as follows:
B = P*D^(-1/2)*P^T

where A = P*D*P^T, D is the diagonal matrix of eigenvalues of A.

I will not need to know P or D after I obtain B.

I would also like to know if there are refinements to this process in the case that A is a symmetric banded matrix, though this is less important.

I'm considering using the VBA version of AlgLib, but if you were to point me to any of the routines in any language I can find the equivalent.


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 44 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group