I think the problem with this version here is VectorDomain(8,lower=[0,0,0,0,-Inf,-Inf,-Inf,-Inf]). Could you please try this version with something like VectorDomain(8,lower=[0,0,0,0,-1e6,-1e6,-1e6,-1e6]) or alternatively VectorDomain(8,lower=[0,0,0,0,-Pumas.∞,-Pumas.∞,-Pumas.∞,-Pumas.∞])? Using Inf in the bounds causes Inf to propagate through the computations causing the error you are getting.