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

Submitted by:  Patrick Alken <psa>
Submitted on:  Mon 30 Nov 2015 11:37:10 PM UTC  
 
Category: Runtime errorSeverity: 3 - Normal
Operating System: Status: Wont Fix
Assigned to: NoneOpen/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>
Project 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>
Project Administrator

 

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
  • -unavailable- added by psa (Submitted the item)
  •  

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 3 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Tue 01 Dec 2015 04:06:27 PM UTCpsaStatusConfirmed=>Wont Fix
      Open/ClosedOpen=>Closed
    Mon 30 Nov 2015 11:37:10 PM UTCpsaAttached File-=>Added main.c, #35604

    Back to the top


    Powered by Savane 3.1-cleanup1