1D interpolation with eval_linear about interpolation.py

[albop]

Here is a working version:

grid = np.linspace(0., np.pi, 10),
value = np.sin(grid[0])                   # a 1d array
xx = np.linspace(0.,np.pi,2)          # a 1d array
eval_linear(grid, value, xx[:,None]) # xx[:,None] is a 2x1 array

[albop]

In the example you give, grid is 1d so values should be représented by a 1d Vector. In your example it id a 2d Vector. Also, grid should be a tupe of 1d components.

[grburgess]

@albop here I will try a better example.

1. the comma after the grid definition packs the grid in a tuple

grid = np.linspace(0., np.pi, 10),

value = np.sin(grid).reshape((1,10))

xx = np.linspace(0.,np.pi,2)

eval_linear(grid, value, np.atleast_1d(xx))

output:

array([0.00000000e+00, 3.42020143e-01, 6.42787610e-01, 8.66025404e-01,
       9.84807753e-01, 9.84807753e-01, 8.66025404e-01, 6.42787610e-01,
       3.42020143e-01, 1.22464680e-16])

It seems it only evals at the grid point and not the new points.

[albop]

Second instruction still looks looks wrong to me: - how Can the argument to np.sin be a tupe ? It shouldn't work. - values should be a 1d array with 10 values not a 1 times 10 Matrix. As for the last one since you interpolate on 10 1d points, you should pass a 10 Times 1 Matrix. Le mar. 22 juin 2021 à 09:48, J. Michael Burgess ***@***.***> a écrit :

…

@albop <https://github.com/albop> here I will try a better example. 1. the comma after the grid definition packs the grid in a tuple grid = np.linspace(0., np.pi, 10), value = np.sin(grid).reshape((1,10)) xx = np.linspace(0.,np.pi,2) eval_linear(grid, value, np.atleast_1d(xx)) output: array([0.00000000e+00, 3.42020143e-01, 6.42787610e-01, 8.66025404e-01, 9.84807753e-01, 9.84807753e-01, 8.66025404e-01, 6.42787610e-01, 3.42020143e-01, 1.22464680e-16]) It seems it only evals at the grid point and not the new points. — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#89 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AACDSKJICNGNZSLD3NBJNN3TUA56HANCNFSM47B3OCPQ> .

[grburgess]

Thanks! Now I get the logic