// Copyright (c) 2021 CINN Authors. All Rights Reserved.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

#include <gtest/gtest.h>

#include "paddle/cinn/frontend/decomposer/test_helper.h"

namespace cinn::frontend {

TEST(Decomposer, elementwise_add_bcast0) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {4, 1, 20, 10});
  auto y = builder.CreateInput(Float(32), {10, 20});
  auto out = builder.Add(x, y, 1);

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{4, 10, 20, 10}};
  RunAndCheckShape<float>(&builder, input_names, output_names, output_shapes);
}

TEST(Decomposer, elementwise_add_bcase1) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {10, 20});
  auto y = builder.CreateInput(Float(32), {4, 1, 20, 10});
  auto out = builder.Add(x, y, 1);

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{4, 10, 20, 10}};
  RunAndCheckShape<float>(&builder, input_names, output_names, output_shapes);
}

TEST(Decomposer, elementwise_add_grad_bcast0) {
  NetBuilder builder("elementwise_add_grad");
  auto dout = builder.CreateInput(Float(32), {4, 10, 20, 10});
  auto x = builder.CreateInput(Float(32), {4, 1, 20, 10});
  auto y = builder.CreateInput(Float(32), {10, 20});
  auto out_grads = builder.ElementwiseAddGrad(dout, x, y, 1);

  std::vector<std::string> input_names = {dout.id().data()};
  std::vector<std::string> output_names = {out_grads[0]->id, out_grads[1]->id};
  std::vector<std::vector<int>> output_shapes = {{4, 1, 20, 10}, {10, 20}};
  RunAndCheckShape<float>(&builder, input_names, output_names, output_shapes);
}

TEST(Decomposer, elementwise_add_bcast1) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {32, 64, 32, 32});
  auto y = builder.CreateInput(Float(32), {64});
  auto out = builder.Add(x, y, 1);

  auto add_cpu = [](const std::vector<size_t>& lengths,
                    const std::vector<void*>& ptrs) {
    float* x = static_cast<float*>(ptrs[0]);
    float* y = static_cast<float*>(ptrs[1]);
    float* out = static_cast<float*>(ptrs[2]);
    for (size_t i = 0; i < 32; ++i) {
      for (size_t j = 0; j < 64; ++j) {
        for (size_t k = 0; k < 32 * 32; ++k) {
          out[(i * 64 + j) * 32 * 32 + k] =
              x[(i * 64 + j) * 32 * 32 + k] + y[j];
        }
      }
    }
  };

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{32, 64, 32, 32}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_cpu);
}

TEST(Decomposer, elementwise_add_bcast1_2) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {64});
  auto y = builder.CreateInput(Float(32), {32, 64, 32, 32});
  auto out = builder.Add(x, y, 1);

  auto add_cpu = [](const std::vector<size_t>& lengths,
                    const std::vector<void*>& ptrs) {
    float* x = static_cast<float*>(ptrs[0]);
    float* y = static_cast<float*>(ptrs[1]);
    float* out = static_cast<float*>(ptrs[2]);
    for (size_t i = 0; i < 32; ++i) {
      for (size_t j = 0; j < 64; ++j) {
        for (size_t k = 0; k < 32 * 32; ++k) {
          out[(i * 64 + j) * 32 * 32 + k] =
              y[(i * 64 + j) * 32 * 32 + k] + x[j];
        }
      }
    }
  };

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{32, 64, 32, 32}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_cpu);
}

TEST(Decomposer, elementwise_add_grad_bcast1) {
  NetBuilder builder("elementwise_add_grad");
  auto dout = builder.CreateInput(Float(32), {32, 64, 32, 32});
  auto x = builder.CreateInput(Float(32), {32, 64, 32, 32});
  auto y = builder.CreateInput(Float(32), {64});
  auto out_grads = builder.ElementwiseAddGrad(dout, x, y, 1);

  auto add_grad_cpu = [](const std::vector<size_t>& lengths,
                         const std::vector<void*>& ptrs) {
    float* dout = static_cast<float*>(ptrs[0]);
    float* dx = static_cast<float*>(ptrs[1]);
    float* dy = static_cast<float*>(ptrs[2]);
    for (size_t j = 0; j < 64; ++j) {
      dy[j] = 0;
    }
    for (size_t i = 0; i < 32; ++i) {
      for (size_t j = 0; j < 64; ++j) {
        for (size_t k = 0; k < 32 * 32; ++k) {
          dx[(i * 64 + j) * 32 * 32 + k] = dout[(i * 64 + j) * 32 * 32 + k];
          dy[j] = dy[j] + dout[(i * 64 + j) * 32 * 32 + k];
        }
      }
    }
  };

  std::vector<std::string> input_names = {dout.id().data()};
  std::vector<std::string> output_names = {out_grads[0]->id, out_grads[1]->id};
  std::vector<std::vector<int>> output_shapes = {{32, 64, 32, 32}, {64}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_grad_cpu);
}

TEST(Decomposer, elementwise_add_bcast2) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {32, 16});
  auto y = builder.CreateInput(Float(32), {1});
  auto out = builder.Add(x, y);

  auto add_cpu = [](const std::vector<size_t>& lengths,
                    const std::vector<void*>& ptrs) {
    size_t n = lengths[0];
    float* x = static_cast<float*>(ptrs[0]);
    float* y = static_cast<float*>(ptrs[1]);
    float* out = static_cast<float*>(ptrs[2]);
    float y_data = y[0];
    for (size_t i = 0; i < n; ++i) {
      out[i] = x[i] + y_data;
    }
  };

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{32, 16}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_cpu);
}

TEST(Decomposer, elementwise_add_bcast2_2) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {1});
  auto y = builder.CreateInput(Float(32), {32, 16});
  auto out = builder.Add(x, y);

  auto add_cpu = [](const std::vector<size_t>& lengths,
                    const std::vector<void*>& ptrs) {
    size_t n = 32 * 16;
    float* x = static_cast<float*>(ptrs[0]);
    float* y = static_cast<float*>(ptrs[1]);
    float* out = static_cast<float*>(ptrs[2]);
    float x_data = x[0];
    for (size_t i = 0; i < n; ++i) {
      out[i] = y[i] + x_data;
    }
  };

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{32, 16}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_cpu);
}

TEST(Decomposer, elementwise_add_bcast2_3) {
  constexpr int kLength = 64;
  using int_ty = int64_t;
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Int(kLength), {32, 16});
  auto y = builder.CreateInput(Int(kLength), {1});
  auto out = builder.Add(x, y);

  auto add_cpu = [](const std::vector<size_t>& lengths,
                    const std::vector<void*>& ptrs) {
    size_t n = lengths[0];
    int_ty* x = static_cast<int_ty*>(ptrs[0]);
    int_ty* y = static_cast<int_ty*>(ptrs[1]);
    int_ty* out = static_cast<int_ty*>(ptrs[2]);
    int_ty y_data = y[0];
    for (size_t i = 0; i < n; ++i) {
      out[i] = x[i] + y_data;
    }
  };

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{32, 16}};
  RunAndCheck<int_ty>(
      &builder, input_names, output_names, output_shapes, add_cpu);
}

TEST(Decomposer, elementwise_add_grad_bcast2) {
  NetBuilder builder("elementwise_add_grad");
  auto dout = builder.CreateInput(Float(32), {32, 16});
  auto x = builder.CreateInput(Float(32), {32, 16});
  auto y = builder.CreateInput(Float(32), {1});
  auto out_grads = builder.ElementwiseAddGrad(dout, x, y);

  auto add_grad_cpu = [](const std::vector<size_t>& lengths,
                         const std::vector<void*>& ptrs) {
    size_t n = lengths[0];
    float* dout = static_cast<float*>(ptrs[0]);
    float* dx = static_cast<float*>(ptrs[1]);
    float* dy = static_cast<float*>(ptrs[2]);
    for (size_t i = 0; i < n; ++i) {
      float tmp = dout[i];
      dx[i] = tmp;
      dy[0] += tmp;
    }
  };

  std::vector<std::string> input_names = {dout.id().data()};
  std::vector<std::string> output_names = {out_grads[0]->id, out_grads[1]->id};
  std::vector<std::vector<int>> output_shapes = {{32, 16}, {1}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_grad_cpu);
}

TEST(Decomposer, elementwise_add_same_dims) {
  NetBuilder builder("elementwise_add");
  auto x = builder.CreateInput(Float(32), {32, 16});
  auto y = builder.CreateInput(Float(32), {32, 16});
  auto out = builder.Add(x, y);

  auto add_cpu = [](const std::vector<size_t>& lengths,
                    const std::vector<void*>& ptrs) {
    size_t n = lengths[0];
    float* x = static_cast<float*>(ptrs[0]);
    float* y = static_cast<float*>(ptrs[1]);
    float* out = static_cast<float*>(ptrs[2]);
    for (size_t i = 0; i < n; ++i) {
      out[i] = x[i] + y[i];
    }
  };

  std::vector<std::string> input_names = {x.id().data(), y.id().data()};
  std::vector<std::string> output_names = {out->id};
  std::vector<std::vector<int>> output_shapes = {{32, 16}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_cpu);
}

TEST(Decomposer, elementwise_add_grad_same_dims) {
  NetBuilder builder("elementwise_add_grad");
  auto dout = builder.CreateInput(Float(32), {32, 16});
  auto x = builder.CreateInput(Float(32), {32, 16});
  auto y = builder.CreateInput(Float(32), {32, 16});
  auto out_grads = builder.ElementwiseAddGrad(dout, x, y);

  auto add_grad_cpu = [](const std::vector<size_t>& lengths,
                         const std::vector<void*>& ptrs) {
    size_t n = lengths[0];
    float* dout = static_cast<float*>(ptrs[0]);
    float* dx = static_cast<float*>(ptrs[1]);
    float* dy = static_cast<float*>(ptrs[2]);
    for (size_t i = 0; i < n; ++i) {
      float tmp = dout[i];
      dx[i] = tmp;
      dy[i] = tmp;
    }
  };

  std::vector<std::string> input_names = {dout.id().data()};
  std::vector<std::string> output_names = {out_grads[0]->id, out_grads[1]->id};
  std::vector<std::vector<int>> output_shapes = {{32, 16}, {32, 16}};
  RunAndCheck<float>(
      &builder, input_names, output_names, output_shapes, add_grad_cpu);
}

}  // namespace cinn::frontend