Connecting PyMC3 to external code - help with understanding Theano custom Ops

Any insight on why it is behaving that way? My guess is that the underlying algorithm can only handle scalars, and leave out information about other variables. Probably it works well enough to make a guess about what to try next in the parameter realm, in most practical cases. If that is correct, then better make sure the first output variable has a significant weight in the likelihood. Just thinking loud…