Skip to main content

Interpolation function precision in multi dimentional case



I have this interpolation example


zzz = Interpolation[{
    {{1}, 2`100*Pi},
    {{2}, 3`100*Pi},
    {{3}, 4`100*Pi}
    }, InterpolationOrder -> 1];
zzz[2]
zzz[2.5]
zzz[2.5`100]


and I am getting these resilts


9.424777960769379715387930149838508652591508198125317462924833776923449218858626995884104476026351204

10.9956

10.99557428756427633461925184147826009469009289781287037341230607307735742200173149519812188869740974

first number as expected with 100 digits, second with only 4 digits, but if I add `100 to the argument I can get it with 100 digits as well.


Now I use multidimensional example


testData = {

   {{-2.*10^-6, 0, -18}, -1.00853429340956742965680093076966025423734590973983188199800080460508043197286933569747093680631648*10^6},
   {{-2.*10^-6, 0, -17}, 
    952504.61044236930112440429386078472672512393090882097358420180063200222870428099604804394164585046},
   {{-2.*10^-6, 1, -18}, -1.008544941816832295575440188033483155333800668013690157106923425838109768669333223372232072904210136*10^6},
   {{-2.*10^-6, 1, -17}, -952514.6672734985329986920711230513627204920640531655413765067694156739368659196153982113188871999717},

   {{-1.78947*10^-6, 0, -18}, -1.00854109353837536345550608172429712197251598047769246281925070613484574642246277456744911673832572*10^6},
   {{-1.78947*10^-6, 0, -17}, -952511.41057117723492310944481542159446029400164668155440545170216176754315387443491802212157785970},

   {{-1.78947*10^-6, 1, -18}, -1.008551741385092646935524440345434727928230066236270646432464156954808417875408736480254883876917917*10^6},

   {{-1.78947*10^-6, 1, -17}, -952521.4668417583449142805833940440259658101867440498419481296004930320450422321944211465035742782570}
   };
ff[x1_, x2_, x3_] = Interpolation[testData, InterpolationOrder -> 1][x1, x2, x3]

and for this example, I cannot get the results with 100 digits, even if I use exact existing values


ff[-2.0`100*10`100^-6, 0.0`100, -18.0`100]
ff[-2.0`100*10`100^-6, 1.0`100, -18.0`100]

Out[350]= -1.00853*10^6
Out[351]= -1.00854*10^6


How to force interpolation to use required precision? There is no WorkingPrecision argument for interpolation.


UPDATE 1


even if I put `100 everywhere like this, it is still the same result


testData = {
   {{-2`100*10^-6, 0.0`100, -18`100}, -1.00853429340956742965680093076966025423734590973983188199800080460508043197286933569747093680631648`100*10^6},
   {{-2`100*10^-6, 0.0`100, -17`100}, 952504.61044236930112440429386078472672512393090882097358420180063200222870428099604804394164585046`100},

   {{-2`100*10^-6, 1`100, -18`100}, -1.008544941816832295575440188033483155333800668013690157106923425838109768669333223372232072904210136`100*10^6},
   {{-2`100*10^-6, 1`100, -17`100}, -952514.6672734985329986920711230513627204920640531655413765067694156739368659196153982113188871999717`100},


   {{-1.78947`100*10^-6, 0.0`100, -18`100}, -1.00854109353837536345550608172429712197251598047769246281925070613484574642246277456744911673832572`100*10^6},
   {{-1.78947`100*10^-6, 0.`100, -17`100}, -952511.41057117723492310944481542159446029400164668155440545170216176754315387443491802212157785970`100},

   {{-1.78947`100*10^-6, 1`100, -18`100}, -1.008551741385092646935524440345434727928230066236270646432464156954808417875408736480254883876917917`100*10^6},
   {{-1.78947`100*10^-6, 1`100, -17`100}, -952521.4668417583449142805833940440259658101867440498419481296004930320450422321944211465035742782570`100}
};
ff[x1_, x2_, x3_] = 
 Interpolation[testData, InterpolationOrder -> 1][x1, x2, x3]

Answer




Evaluate 0.0`100 and you'll get a machine precision zero. Use the exact 0 instead:


fff = Interpolation[testData /. x_ /; x == 0 -> 0, InterpolationOrder -> 1];

fff[-2.`100.*^-6, 0, -18.`100.]
(*
-1.0085342934095674296568009307696602542373459097398318819980008046050\
804319728693356974709368063165*10^6
*)

Comments

Post a Comment

Popular posts from this blog

plotting - How to draw lines between specified dots on ListPlot?

I would like to create a plot where I have unconnected dots and some connected. So far, I have figured out how to draw the dots. My code is the following: ListPlot[{{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4,13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full] I have thought using ListLinePlot command, but I don't know how to specify to the command to draw only selected lines between the dots. Do have any suggestions/hints on how to do that? Thank you. Answer One possibility would be to use Epilog with Line : ListPlot[ {{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4, 13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full, Epilog -> { Line[ ...
Post a Comment
Read more

list manipulation - Selecting multiple columns from a matrix?

Sample data: data = { {{2013, 1, 1}, 24.13, 167.67, 231.82}, {{2013, 1, 2}, 32.15, 170.92, 225.99}, {{2013, 1, 3}, 35.43, 172.68, 221.67}, {{2013, 1, 4}, 36.73, 173.05, 218.32}, {{2013, 1, 5}, 58.19, 165.96, 197.05}, {{2013, 1, 6}, 69.99, 163.50, 187.52}, {{2013, 1, 7}, 71.37, 154.21, 175.58}, {{2013, 1, 8}, 72.51, 149.66, 163.25}}; I want a DateListPlot with three graphs, so for a matrix formed by columns 1 and 2, one for columns 1 and 3, and 1 for columns 1 and 4. At the moment I'm using this code: data2 = Transpose[{data[[All, 1]], data[[All, 2]]}]; data3 = Transpose[{data[[All, 1]], data[[All, 3]]}]; data4 = Transpose[{data[[All, 1]], data[[All, 4]]}]; DateListPlot[{data2, data3, data4}, Joined -> True, Filling -> {3 -> {1}}] but I have a hunch that this can be done more efficiently. I don't like the Transpose s in particular. Any ideas? edit (for extra credit) What if I need to multiply the second column by 2, which in my solution is simp...
Post a Comment
Read more

dynamic - How can I make a clickable ArrayPlot that returns input?

I would like to create a dynamic ArrayPlot so that the rectangles, when clicked, provide the input. Can I use ArrayPlot for this? Or is there something else I should have to use? Answer ArrayPlot is much more than just a simple array like Grid : it represents a ranged 2D dataset, and its visualization can be finetuned by options like DataReversed and DataRange . These features make it quite complicated to reproduce the same layout and order with Grid . Here I offer AnnotatedArrayPlot which comes in handy when your dataset is more than just a flat 2D array. The dynamic interface allows highlighting individual cells and possibly interacting with them. AnnotatedArrayPlot works the same way as ArrayPlot and accepts the same options plus Enabled , HighlightCoordinates , HighlightStyle and HighlightElementFunction . data = {{Missing["HasSomeMoreData"], GrayLevel[ 1], {RGBColor[0, 1, 1], RGBColor[0, 0, 1], GrayLevel[1]}, RGBColor[0, 1, 0]}, {GrayLevel[0], GrayLevel...
Post a Comment
Read more