#include <migraphx/rewrite_batchnorm.hpp>
#include <migraphx/program.hpp>
#include <migraphx/ref/target.hpp>
#include <migraphx/op/convolution.hpp>
#include <migraphx/op/reshape.hpp>
#include <migraphx/op/batch_norm_inference.hpp>
#include <migraphx/instruction.hpp>
#include <migraphx/generate.hpp>
#include <migraphx/ranges.hpp>
#include <test.hpp>
#include <migraphx/make_op.hpp>

#include <migraphx/serialize.hpp>

#include <migraphx/verify.hpp>

bool is_batch_norm(migraphx::instruction& ins) { return ins.name() == "batch_norm_inference"; }

TEST_CASE(fwd_conv_batchnorm_rewrite_test)
{
    std::vector<float> xdata = {
        0.26485917, 0.61703885, 0.32762103, 0.2503367,  0.6552712,  0.07947932, 0.95442678,
        0.70892651, 0.890563,   0.80808088, 0.89540492, 0.52657048, 0.94614791, 0.64371508,
        0.0971229,  0.2475562,  0.47405955, 0.85538928, 0.05428386, 0.993078,   0.72771973,
        0.18312255, 0.3091522,  0.51396558, 0.35158192, 0.2419852,  0.83691474, 0.36355352,
        0.04769134, 0.08312604, 0.61804092, 0.0508887,  0.30987137, 0.81307629, 0.16398955,
        0.69886166, 0.02415926, 0.60608918, 0.81907569, 0.13208211, 0.48303735, 0.87533734,
        0.92998813, 0.65553674, 0.73223327, 0.99401001, 0.09850688, 0.76972609, 0.11118327,
        0.04392097, 0.39252306, 0.91129653, 0.89078693, 0.60571206, 0.98410397, 0.15290698,
        0.86992609, 0.7575111,  0.80583525, 0.23649562, 0.7478029,  0.62888878, 0.39886601,
        0.37066793, 0.72627947, 0.8745595,  0.13568234, 0.7413787,  0.5039495,  0.18945697,
        0.87046838, 0.63970494, 0.01124038, 0.27459063, 0.65745586, 0.69182619, 0.80470603,
        0.58039348, 0.36950583, 0.43634225, 0.01694425, 0.14099377, 0.77015849, 0.35809292,
        0.40547674, 0.46538817, 0.65835358, 0.2266954,  0.39057646, 0.64642207, 0.84491134,
        0.20998067, 0.41074121, 0.73055221, 0.26424874, 0.10612507, 0.24478521, 0.24091282,
        0.52536754, 0.57292341, 0.82190903, 0.51858515, 0.17162996, 0.52048114, 0.96624787,
        0.17527163, 0.56384485, 0.91991603};
    std::vector<float> wdata = {
        -1.12125056, 0.50228441,  1.12719446,  -2.61705068, -0.2027315,  -0.82199441, 0.05337102,
        -0.62146691, -2.40572931, -1.47175612, 1.49654601,  -1.07070376, -0.65908074, -0.28457694,
        1.60046717,  0.20677642,  -1.51844486, 0.41203847,  -0.01285751, 0.07948031,  -0.91507006,
        -1.59481079, -0.12856238, 0.39970482,  -1.89015158, 0.66969754,  0.10312618};
    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 6, 6}};
    migraphx::shape ws{migraphx::shape::float_type, {1, 3, 3, 3}};
    migraphx::shape vars{migraphx::shape::float_type, {1}};

    auto create_program = [&]() {
        migraphx::program p;

        auto* mm  = p.get_main_module();
        auto x    = mm->add_literal(xs, xdata);
        auto w    = mm->add_literal(ws, wdata);
        auto conv = mm->add_instruction(
            migraphx::make_op("convolution",
                              {{"padding", {0, 0}}, {"stride", {1, 1}}, {"dilation", {1, 1}}}),
            x,
            w);
        auto scale    = mm->add_literal(migraphx::literal{vars, {3.0f}});
        auto bias     = mm->add_literal(migraphx::literal{vars, {8.1f}});
        auto mean     = mm->add_literal(migraphx::literal{vars, {4.0f}});
        auto variance = mm->add_literal(migraphx::literal{vars, {37.11f}});
        mm->add_instruction(
            migraphx::make_op("batch_norm_inference"), conv, scale, bias, mean, variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();

    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    p1.compile(migraphx::ref::target{});
    p2.compile(migraphx::ref::target{});

    auto result1 = p1.eval({}).back();
    auto result2 = p2.eval({}).back();

    std::vector<float> results_vector1;
    std::vector<float> results_vector2;
    result1.visit([&](auto output) { results_vector1.assign(output.begin(), output.end()); });
    result2.visit([&](auto output) { results_vector2.assign(output.begin(), output.end()); });
    EXPECT(migraphx::verify_range(results_vector1, results_vector2));
}

TEST_CASE(non_literal)
{

    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 8, 8}};
    migraphx::shape ws{migraphx::shape::float_type, {4, 3, 1, 1}};
    migraphx::shape vars{migraphx::shape::float_type, {4}};
    auto create_program = [&]() {
        migraphx::program p;
        auto* mm      = p.get_main_module();
        auto x        = mm->add_parameter("x", xs);
        auto w        = mm->add_parameter("w", ws);
        auto conv     = mm->add_instruction(migraphx::make_op("convolution"), x, w);
        auto scale    = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 1)));
        auto bias     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 2)));
        auto mean     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 3)));
        auto variance = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 4)));
        mm->add_instruction(
            migraphx::make_op("batch_norm_inference"), conv, scale, bias, mean, variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();

    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    EXPECT(any_of(p1, &is_batch_norm));
    EXPECT(none_of(p2, &is_batch_norm));
}

TEST_CASE(as_literal)
{

    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 8, 8}};
    migraphx::shape ws{migraphx::shape::float_type, {4, 3, 1, 1}};
    migraphx::shape vars{migraphx::shape::float_type, {4}};
    auto create_program = [&]() {
        migraphx::program p;
        auto* mm      = p.get_main_module();
        auto x        = mm->add_literal(migraphx::generate_literal(xs, 1));
        auto w        = mm->add_literal(migraphx::generate_literal(ws, 1));
        auto conv     = mm->add_instruction(migraphx::make_op("convolution"), x, w);
        auto scale    = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 1)));
        auto bias     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 2)));
        auto mean     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 3)));
        auto variance = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 4)));
        mm->add_instruction(
            migraphx::make_op("batch_norm_inference"), conv, scale, bias, mean, variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();
    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    EXPECT(any_of(p1, &is_batch_norm));
    EXPECT(none_of(p2, &is_batch_norm));

    p1.compile(migraphx::ref::target{});
    p2.compile(migraphx::ref::target{});

    auto result1 = p1.eval({}).back();
    auto result2 = p2.eval({}).back();
    visit_all(result1, result2)([&](auto r1, auto r2) { EXPECT(migraphx::verify_range(r1, r2)); });
}

TEST_CASE(as_literal_1d)
{
    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 8}};
    migraphx::shape ws{migraphx::shape::float_type, {4, 3, 1}};
    migraphx::shape vars{migraphx::shape::float_type, {4}};
    auto create_program = [&]() {
        migraphx::program p;
        auto* mm  = p.get_main_module();
        auto x    = mm->add_literal(migraphx::generate_literal(xs, 1));
        auto w    = mm->add_literal(migraphx::generate_literal(ws, 1));
        auto conv = mm->add_instruction(
            migraphx::make_op("convolution",
                              {{"padding", {0}}, {"stride", {1}}, {"dilation", {1}}}),
            x,
            w);
        auto scale    = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 1)));
        auto bias     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 2)));
        auto mean     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 3)));
        auto variance = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 4)));
        mm->add_instruction(
            migraphx::make_op("batch_norm_inference"), conv, scale, bias, mean, variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();
    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    EXPECT(any_of(p1, &is_batch_norm));
    EXPECT(none_of(p2, &is_batch_norm));

    p1.compile(migraphx::ref::target{});
    p2.compile(migraphx::ref::target{});

    auto result1 = p1.eval({}).back();
    auto result2 = p2.eval({}).back();
    visit_all(result1, result2)([&](auto r1, auto r2) { EXPECT(migraphx::verify_range(r1, r2)); });
}

TEST_CASE(as_literal_3d)
{
    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 2, 4, 8}};
    migraphx::shape ws{migraphx::shape::float_type, {4, 3, 1, 1, 1}};
    migraphx::shape vars{migraphx::shape::float_type, {4}};
    auto create_program = [&]() {
        migraphx::program p;
        auto* mm = p.get_main_module();
        migraphx::op::convolution conv_op;
        conv_op.padding  = {0, 0, 0};
        conv_op.stride   = {1, 1, 1};
        conv_op.dilation = {1, 1, 1};

        auto x        = mm->add_literal(migraphx::generate_literal(xs, 1));
        auto w        = mm->add_literal(migraphx::generate_literal(ws, 1));
        auto conv     = mm->add_instruction(conv_op, x, w);
        auto scale    = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 1)));
        auto bias     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 2)));
        auto mean     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 3)));
        auto variance = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 4)));
        mm->add_instruction(
            migraphx::make_op("batch_norm_inference"), conv, scale, bias, mean, variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();
    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    EXPECT(any_of(p1, &is_batch_norm));
    EXPECT(none_of(p2, &is_batch_norm));

    p1.compile(migraphx::ref::target{});
    p2.compile(migraphx::ref::target{});

    auto result1 = p1.eval({}).back();
    auto result2 = p2.eval({}).back();
    visit_all(result1, result2)([&](auto r1, auto r2) { EXPECT(migraphx::verify_range(r1, r2)); });
}

TEST_CASE(literal_reshape)
{
    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 8, 8}};
    migraphx::shape ws{migraphx::shape::float_type, {4, 3, 1, 1}};
    migraphx::shape vars{migraphx::shape::float_type, {4}};

    auto create_program = [&]() {
        migraphx::program p;
        auto* mm      = p.get_main_module();
        auto x        = mm->add_literal(migraphx::generate_literal(xs, 1));
        auto w        = mm->add_literal(migraphx::generate_literal(ws, 1));
        auto conv     = mm->add_instruction(migraphx::make_op("convolution"), x, w);
        auto scale    = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 1)));
        auto bias     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 2)));
        auto mean     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 3)));
        auto variance = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 4)));
        mm->add_instruction(
            migraphx::make_op("batch_norm_inference"), conv, scale, bias, mean, variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();
    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    EXPECT(any_of(p1, &is_batch_norm));
    EXPECT(none_of(p2, &is_batch_norm));

    p1.compile(migraphx::ref::target{});
    p2.compile(migraphx::ref::target{});

    auto result1 = p1.eval({}).back();
    auto result2 = p2.eval({}).back();
    visit_all(result1, result2)([&](auto r1, auto r2) { EXPECT(migraphx::verify_range(r1, r2)); });
}

TEST_CASE(literal_reshape_per_actv)
{
    migraphx::shape xs{migraphx::shape::float_type, {1, 3, 8, 7, 4}};
    migraphx::shape ws{migraphx::shape::float_type, {4, 3, 1, 1, 1}};
    migraphx::shape vars{migraphx::shape::float_type, {4, 8, 7, 4}};

    auto create_program = [&]() {
        migraphx::program p;
        auto* mm  = p.get_main_module();
        auto x    = mm->add_literal(migraphx::generate_literal(xs, 1));
        auto w    = mm->add_literal(migraphx::generate_literal(ws, 1));
        auto conv = mm->add_instruction(
            migraphx::make_op(
                "convolution",
                {{"padding", {0, 0, 0}}, {"stride", {1, 1, 1}}, {"dilation", {1, 1, 1}}}),
            x,
            w);
        auto scale    = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 1)));
        auto bias     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 2)));
        auto mean     = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 3)));
        auto variance = mm->add_literal(migraphx::abs(migraphx::generate_literal(vars, 4)));
        mm->add_instruction(
            migraphx::make_op(
                "batch_norm_inference",
                {{"epsilon", 1.0e-5},
                 {"momentum", 0.88},
                 {"bn_mode",
                  migraphx::to_value(migraphx::op::batch_norm_inference::per_activation)}}),
            conv,
            scale,
            bias,
            mean,
            variance);
        return p;
    };

    migraphx::program p1 = create_program();
    migraphx::program p2 = create_program();
    migraphx::rewrite_batchnorm opt;
    opt.apply(*p2.get_main_module());
    EXPECT(any_of(p1, &is_batch_norm));
    EXPECT(none_of(p2, &is_batch_norm));

    p1.compile(migraphx::ref::target{});
    p2.compile(migraphx::ref::target{});

    auto result1 = p1.eval({}).back();
    auto result2 = p2.eval({}).back();
    visit_all(result1, result2)([&](auto r1, auto r2) { EXPECT(migraphx::verify_range(r1, r2)); });
}

int main(int argc, const char* argv[]) { test::run(argc, argv); }
