Making pymc wrapper for finite element

You are passing a function object to mu in your likelihood term. This is what is producing the error. You need to actually call the function on inputs a0, a1, a2