 }
 
 inline __host__ __device__ half operator*(half a, half b) {
-  return mul(a, b);
+  return THCNumerics<half>::mul(a, b);
 }
 
 inline __host__ __device__ float operator*(half a, float b) {

}

inline __host__ __device__ half operator/(half a, half b) {
  return div(a, b);
}

inline __host__ __device__ float operator/(float a, half b) {


#ifdef CUDA_HALF_TENSOR

// arithmetic functions

inline __host__ __device__ half operator+(half a, half b) {
