The rbf function takes lenghtscale parameter - which is what you are modeling in the beginning. You should not use it this way, my suggestion is to port the rbf function to theano first (again, it is equivalent as using the pm.gp module).
1 Like