bugGNU Scientific Library - Bugs: bug #46576, polynomial interpolation failure

 
 

You are not allowed to post comments on this tracker with your current authentication level.

bug #46576: polynomial interpolation failure

Submitter:  Patrick Alken <psa>
Submitted:  Mon 30 Nov 2015 11:37:10 PM UTC
   
 
Category:  Runtime error Severity:  3 - Normal
Operating System:  Status:  Wont Fix
Assigned to:  None Open/Closed:  Closed
Release: 

Tue 01 Dec 2015 04:06:27 PM UTC, comment #1: 

The problem is the program is using a degree 1000 polynomial which cannot be computed with numerical stability. Suggest to use cubic splines instead.

Patrick Alken <psa>
Group administrator
Mon 30 Nov 2015 11:37:10 PM UTC, original submission:  

from adam =dot= petrus =at= baesystems =dot= com

I am attempting to use the polynomial interpolation method with GSL. However I am getting very strange results.

If I run the attached program I get the below result. If I change the method to linear or cspline the interpolation works fine.

Not in previous gsl versions the polynomial interpolation has worked!

[cid:image001.png@01D12B93.B5E118F0]

Is there something I am missing?

Yours aye

Adam

Patrick Alken <psa>
Group administrator

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attached Files
file #35604:  main.c added by psa (785B - text/x-c)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by psa (Submitted the item)
  •  

    Follow 3 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2015-12-01 psa StatusConfirmed Wont Fix
        Open/ClosedOpen Closed
    2015-11-30 psa Attached File- Added main.c, #35604

    Back to the top

    Powered by Savane 3.13-4448.
    Corresponding source code