verify_grad (tensor/basic.py) was implemented to test gradient computation in theano by the Finite Difference Method. Multiple_outputs_numeric_grad (scan/tests/test_scan.py) serves the same scope, but it also supports ops with multiple outputs ( by using a random projection over the outputs and summing them up).
Somebody should double check the multiple_outputs_numeric_grad implementation and replace verify_grad with it.