sddsinterp returns error when interpolating on the grid

Moderators: michael_borland, soliday

Post Reply
wguo
Posts: 41
Joined: 13 Jun 2008, 12:54

sddsinterp returns error when interpolating on the grid

Post by wguo » 13 Sep 2011, 10:10

Hi Michael,

When I use sddsinterp and set -order=4, the program returns an error message: zero denominator in LagrangeInterp when
the interpolation happens to be on the grid. What should I do?

Many thanks,

Weiming

michael_borland
Posts: 1927
Joined: 19 May 2008, 09:33
Location: Argonne National Laboratory
Contact:

Re: sddsinterp returns error when interpolating on the grid

Post by michael_borland » 13 Sep 2011, 10:23

Weiming,

Can you send me (or post) the data and the full command used?

--Michael

wguo
Posts: 41
Joined: 13 Jun 2008, 12:54

Re: sddsinterp returns error when interpolating on the grid

Post by wguo » 13 Sep 2011, 15:04

Hi Michael,

The error circumstance is not what I thought.
Use the attched file: test0.sdds and the following command:

sddsinterp test0.sdds test1.sdds -col=x,Byn10 -atValues=0.00013942073576 -order=4

The reason it quits is because there are many repeated values as following:

x: -2,-1,0,0,0,0,0,0,1,2
Byn10: 4002,4001,4000,4000,4000,4000,4000,4000,3999,3998

Can I claim it as a bug, or, should I fix the data myself?

Many thanks,

Weiming
Attachments
test0.sdds
(27.12 KiB) Downloaded 3281 times

michael_borland
Posts: 1927
Joined: 19 May 2008, 09:33
Location: Argonne National Laboratory
Contact:

Re: sddsinterp returns error when interpolating on the grid

Post by michael_borland » 13 Sep 2011, 15:37

Weiming,

The program doesn't give a clear enough error message, so we'll fix it. Thanks for reporting the problem.

By the way, you can easily fix the problem with
% sddssort test0.sdds -column=x -unique

--Michael

wguo
Posts: 41
Joined: 13 Jun 2008, 12:54

Re: sddsinterp returns error when interpolating on the grid

Post by wguo » 13 Sep 2011, 15:39

Hi Michael,

Thanks for showing me the trick!

Weiming

Post Reply