comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: SVDFIT Problems
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: SVDFIT Problems [message #31936 is a reply to message #31931] Fri, 30 August 2002 10:08 Go to previous message
Paul Van Delst[1] is currently offline  Paul Van Delst[1]
Messages: 1157
Registered: April 2002
Senior Member
Bill wrote:
>
> In my occaisional attempts to use SVDFIT I have had problems in that it
> reports singular solutions when other fitting routines, e.g., REGRESS,
> have had not problems.

Deja vu! I was using the SVDFIT routine just yesterday to fit a second order polynomial. I kept
getting an error report for a singular value (even with /DOUBLE). Used POLY_FIT with no problems
(SVDFIT _was_ overkill for what I wanted to do.)

<snip>

> Unfortunately it appeas that the code currently relies on a simplified
> version of the Numerical Recipes routine SVDCMP, that does not return
> the singular value vector that Numerical Recipes uses to compute wmax. A
> rewrite to give SVDFIT the proper properties is more work than I am
> willing to put in at this time.

Thanks for taking the time to nut the problem out and report it here. It really pains me (just
having finished the yearly license renewal tooth-pulling process) that these sorts of things are
*still* found in incarnations of such well-worn linear algebra problems like SVD. I spent about half
a day trying to figure out what I was doing wrong until I tossed the lot and surfed over to
netlib.org for some fortran routines. I thought IDL was supposed to *save* time for this sort of
thing?

paulv

--
Paul van Delst
CIMSS @ NOAA/NCEP/EMC
Ph: (301)763-8000 x7274
Fax:(301)763-8545
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Previous Topic: CHECK_MATH and exp()
Next Topic: mesh clipping

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 23:30:51 PDT 2025

Total time taken to generate the page: 1.44070 seconds