#ifndef TH_GENERIC_FILE
#define TH_GENERIC_FILE "generic/cpu.c"
#else

void scatter_(add)(THTensor *output, THLongTensor *index, THTensor *input, int dim) {
  int64_t idx;
  TH_TENSOR_DIM_APPLY3(real, output, real, input, int64_t, index, dim,
    TH_TENSOR_DIM_APPLY3_SIZE_EQ_EXCEPT_DIM,
    for (int64_t i = 0; i < THLongTensor_size(index, dim); i++) {
      idx = *(index_data + i * index_stride);
      assertIndexInBoundaries(idx, output_size, TH_TENSOR_DIM_APPLY_counter);
      output_data[idx] += *(input_data + i * input_stride);
    })
}

void scatter_(sub)(THTensor *output, THLongTensor *index, THTensor *input, int dim) {
  int64_t idx;
  TH_TENSOR_DIM_APPLY3(real, output, real, input, int64_t, index, dim,
    TH_TENSOR_DIM_APPLY3_SIZE_EQ_EXCEPT_DIM,
    for (int64_t i = 0; i < THLongTensor_size(index, dim); i++) {
      idx = *(index_data + i * index_stride);
      assertIndexInBoundaries(idx, output_size, TH_TENSOR_DIM_APPLY_counter);
      output_data[idx] += *(input_data + i * input_stride);
    })
}

void scatter_(mul)(THTensor *output, THLongTensor *index, THTensor *input, int dim) {
  int64_t idx;
  TH_TENSOR_DIM_APPLY3(real, output, real, input, int64_t, index, dim,
    TH_TENSOR_DIM_APPLY3_SIZE_EQ_EXCEPT_DIM,
    for (int64_t i = 0; i < THLongTensor_size(index, dim); i++) {
      idx = *(index_data + i * index_stride);
      assertIndexInBoundaries(idx, output_size, TH_TENSOR_DIM_APPLY_counter);
      output_data[idx] += *(input_data + i * input_stride);
    })
}

void scatter_(div)(THTensor *output, THLongTensor *index, THTensor *input, int dim) {
  int64_t idx;
  TH_TENSOR_DIM_APPLY3(real, output, real, input, int64_t, index, dim,
    TH_TENSOR_DIM_APPLY3_SIZE_EQ_EXCEPT_DIM,
    for (int64_t i = 0; i < THLongTensor_size(index, dim); i++) {
      idx = *(index_data + i * index_stride);
      assertIndexInBoundaries(idx, output_size, TH_TENSOR_DIM_APPLY_counter);
      output_data[idx] += *(input_data + i * input_stride);
    })
}

#endif