indexing issue with ArrayParition arguments

[mkarikom]

I'm not sure if this is supported, but it seemed similar to one of the multivariate tests:

MWE:

Y=rand(10)
X=rand(10,3)
J=ones(10)
n=3
ar1 = ArrayPartition(rand(),rand(),rand(3),rand(3,3),rand(),rand())
ar2 = ArrayPartition(rand(),rand(),rand(3),rand(3,3),rand(),rand())
function Lik(X,Y,J,n,ap)
    lik = (ap.x[5]*ap.x[6]+(Y-ap.x[1]*J-X*ap.x[3])'*inv(I+ap.x[2]*J*J'+X*ap.x[4]*X')*(Y-ap.x[1]*J-X*ap.x[3]))^(-(n+ap.x[5])/2)
end

Lik(X,Y,J,n,ar1)

0.33995995508984683

res = optimize(ap -> Lik(X,Y,J,n,ap), ar2):

ERROR: MethodError: no method matching setindex!(::Float64, ::Float64, ::Int64)
Stacktrace:
 [1] setindex!(::ArrayPartition{Float64,Tuple{Float64,Float64,Array{Float64,1},Array{Float64,2},Float64,Float64}}, ::Float64, ::Int64) at /home/me/.julia/packages/RecursiveArrayTools/OAIEc/src/array_partition.jl:179
 [2] simplexer(::Optim.AffineSimplexer, ::ArrayPartition{Float64,Tuple{Float64,Float64,Array{Float64,1},Array{Float64,2},Float64,Float64}}) at /home/me/.julia/packages/Optim/Agd3B/src/multivariate/solvers/zeroth_order/nelder_mead.jl:13
 [3] initial_state(::NelderMead{Optim.AffineSimplexer,Optim.AdaptiveParameters}, ::Optim.Options{Float64,Nothing}, ::NonDifferentiable{Float64,ArrayPartition{Float64,Tuple{Array{Float64,0},Array{Float64,0},Array{Float64,1},Array{Float64,2},Array{Float64,0},Array{Float64,0}}}}, ::ArrayPartition{Float64,Tuple{Float64,Float64,Array{Float64,1},Array{Float64,2},Float64,Float64}}) at /home/me/.julia/packages/Optim/Agd3B/src/multivariate/solvers/zeroth_order/nelder_mead.jl:155
 [4] optimize(::NonDifferentiable{Float64,ArrayPartition{Float64,Tuple{Array{Float64,0},Array{Float64,0},Array{Float64,1},Array{Float64,2},Array{Float64,0},Array{Float64,0}}}}, ::ArrayPartition{Float64,Tuple{Float64,Float64,Array{Float64,1},Array{Float64,2},Float64,Float64}}, ::NelderMead{Optim.AffineSimplexer,Optim.AdaptiveParameters}, ::Optim.Options{Float64,Nothing}) at /home/me/.julia/packages/Optim/Agd3B/src/multivariate/optimize/optimize.jl:33
 [5] optimize(::Function, ::ArrayPartition{Float64,Tuple{Float64,Float64,Array{Float64,1},Array{Float64,2},Float64,Float64}}; inplace::Bool, autodiff::Symbol, kwargs::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}) at /home/me/.julia/packages/Optim/Agd3B/src/multivariate/optimize/interface.jl:64
 [6] optimize(::Function, ::ArrayPartition{Float64,Tuple{Float64,Float64,Array{Float64,1},Array{Float64,2},Float64,Float64}}) at /home/me/.julia/packages/Optim/Agd3B/src/multivariate/optimize/interface.jl:58
 [7] top-level scope at /mnt/evo/Software_a1/testCUDA/testoptimT.jl:122

[pkofod]

This is because the partitions are scalars not vectors here, so it's trying to setindex! into a scalar. You should use rand(1) and adjust your code accordingly. Thanks!