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
sddsinterp returns error when interpolating on the grid
Moderators: michael_borland, soliday
-
- Posts: 1951
- Joined: 19 May 2008, 09:33
- Location: Argonne National Laboratory
- Contact:
Re: sddsinterp returns error when interpolating on the grid
Weiming,
Can you send me (or post) the data and the full command used?
--Michael
Can you send me (or post) the data and the full command used?
--Michael
Re: sddsinterp returns error when interpolating on the grid
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
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 3311 times
-
- Posts: 1951
- Joined: 19 May 2008, 09:33
- Location: Argonne National Laboratory
- Contact:
Re: sddsinterp returns error when interpolating on the grid
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
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
Re: sddsinterp returns error when interpolating on the grid
Hi Michael,
Thanks for showing me the trick!
Weiming
Thanks for showing me the trick!
Weiming