A= sprand(10000, 10000, 0.1);
II= rand(10000, 2); # typeof(II) is Matrix{Float64} (alias for Array{Float64, 2})
V= zeros(10000, 2);
ps = MKLPardisoSolver()
solve!(ps, V, A, II); # it works
II= vec(II); # typeof(V) is Vector{Float64} (alias for Array{Float64, 1})
V = vec(V);
solve!(ps, V, A, II); # it does not work
julia> ERROR: DimensionMismatch: matrix has 10000 rows, RHS has 20000 rows.