Add support for SclarOp.c_support_code in GpuElemwise.

theano/sandbox/cuda/elemwise.py

@@ -864,7 +864,7 @@ nd_collapse_[i]=0;
         return defines + kernels
 
     def c_support_code(self):
-        raise gof.utils.MethodNotDefined()
+        return self.scalar_op.c_support_code()
 
     def c_code(self, node, nodename, inputs, outputs, sub):
         d = dict(sub)

theano/scalar/basic_scipy.py

@@ -220,9 +220,16 @@ class Psi(UnaryScalarOp):
     def c_support_code(self):
         return (
 """
+// For GPU support
+#ifdef __CUDACC__
+#define DEVICE __device__
+#else
+#define DEVICE
+#endif
+
 #ifndef _PSIFUNCDEFINED
 #define _PSIFUNCDEFINED
-double _psi(double x){
+DEVICE double _psi(double x){
 
     /*taken from
     Bernardo, J. M. (1976). Algorithm AS 103: