朱莉娅 - Interpolations.jl抛出“ERROR：InexactError（）”为载体，而不是点

Question

我插了一个二维矩阵J：

J_itp = interpolate(J, (BSpline(Linear()), NoInterp()), OnGrid()) 

正如你看到的，我只是沿着插第一个维度。第二个维度需要在网格上。我有一个我想插入点的向量。

julia> Np_args[1,67:73] 
1x7 Array{Float64,2}: 
1.0 1.0 1.0 1.0 1.01394 1.03275 1.0517 

我所说的插补功能，像这样：

J_itp[ Np_args[1,67:73], 1] 

，并出现以下错误

ERROR: InexactError() 
in trunc at /Applications/Julia- 0.4.3.app/Contents/Resources/julia/lib/julia/sys.dylib 
[inlined code] from /Users/btengels/.julia/v0.4/Interpolations/src/b-splines/linear.jl:6 
in getindex at /Users/btengels/.julia/v0.4/Interpolations/src/b-splines/indexing.jl:39 

我可以，但是，插值使用我的输入下半部哪里等于一。所以它不可能是一个类型问题（据我所知）。

julia> J_itp[ Np_args[1,67:70], 1] 
4-element Array{Float64,1}: 
0.75366 
0.75366 
0.75366 
0.75366 

当我包含任何大于1.0的点时（如第一个例子），事情会崩溃。然而，我可以自己插入其中一个问题点。

julia> J_itp[ Np_args[1,71], 1] 
0.753702904230203 

但是，如果在像J_itp[ Np_args[1,71]*ones(3,1), 1]这样的数组中，事情会像以前一样崩溃。

这是一个错误还是我做错了什么？最近安装了Pkg.add("Interpolations")，OSX - 优胜美地。 Julia version 0.4.3

Answer 1

有点难以确切说出，因为我不知道J有什么元素类型，因此它没有足够的信息来重建您的问题。但是，我能够触发InexactError从

julia> J = rand(1:10, 100,100); 

julia> J_itp = interpolate(J, (BSpline(Linear()), NoInterp()), OnGrid()); 

julia> J_itp[rand(5)+10, 1] 
ERROR: InexactError() 
in to_index at deprecated.jl:454 
in _unsafe_getindex at multidimensional.jl:192 
in getindex at abstractarray.jl:488 

但这个工程：

Float64[J_itp[x, 1] for x in rand(5)+10] 

（您需要Float64如果你在全球范围内运行这个，如果你把它放在一个函数，你不应该需要它。）

至少在我测试的情况下，错误是由茱莉亚生成的，而不是插值;问题是在julia中支持向量索引，它希望你使用整数值索引。这可能（或可能不）在茱莉亚-0.5中改变。这通过使用理解绕过向量索引。