Can you provide more details about the type of functions you are trying to test?
You can always test the output of a "PyMC"function by compiling your expressions within a theano.function([inputs], output) but maybe you can do something simpler by checking the nodes returned by your function…