| Re: IDL inverse matrix problem?? [message #90079 is a reply to message #90073] |
Fri, 23 January 2015 10:32   |
Craig Markwardt
Messages: 1869
Registered: November 1996
|
Senior Member |
|
|
On Friday, January 23, 2015 at 1:58:37 AM UTC-5, Amin Farhang wrote:
> Dear Craig,
>
> Since my original matrix, which is a positive-definite 1000x1000 matrix, before sending my post, I try to solve the linear equations using Cholesky decomposition, but again IDL could not handle the precision. On the other hand since IDL is widely used in Astronomy and you know the values in Astronomy are huge, I think it is very important to be sure that IDL is reliable for this kind of computations. Anyway thanks and I try to re-write my code in Python or matlab.
Yes, please let us know how it turns out with Python or Matlab. By the way, you were using LA_INVERT, which is from the LAPACK package. This is one of the premier linear algebra package for computers. If LAPACK is having trouble, what does that say about your problem?
Craig
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
