/****************************************************************************
 * Copyright (c) 2017-2022 by the ArborX authors                            *
 * All rights reserved.                                                     *
 *                                                                          *
 * This file is part of the ArborX library. ArborX is                       *
 * distributed under a BSD 3-clause license. For the licensing terms see    *
 * the LICENSE file in the top-level directory.                             *
 *                                                                          *
 * SPDX-License-Identifier: BSD-3-Clause                                    *
 ****************************************************************************/
#include <misc/ArborX_Vector.hpp>

#include <boost/test/unit_test.hpp>

BOOST_AUTO_TEST_CASE(make_euclidean_vector)
{
  using Point = ArborX::Point<3, float>;
  using Vector = ArborX::Details::Vector<3, float>;
  static_assert(Point{1, 2, 3} - Point{0, 0, 0} == Vector{1, 2, 3});
  static_assert(Point{4, 5, 6} - Point{1, 2, 3} == Vector{3, 3, 3});
}

BOOST_AUTO_TEST_CASE(vector_dot_product)
{
  using Vector = ArborX::Details::Vector<3, float>;
  static_assert(Vector{1, 0, 0}.dot(Vector{1, 0, 0}) == 1);
  static_assert(Vector{1, 0, 0}.dot(Vector{0, 1, 0}) == 0);
  static_assert(Vector{1, 0, 0}.dot(Vector{0, 0, 1}) == 0);
  static_assert(Vector{1, 1, 1}.dot(Vector{1, 1, 1}) == 3);

  static_assert(Vector{1, 1e-7, 0}.dot(Vector{1, 1e-7, 0}) == 1);
  static_assert(Vector{1, 1e-7, 0}.dot<double>(Vector{1, 1e-7, 0}) > 1);
}

BOOST_AUTO_TEST_CASE(vector_norm)
{
  using Vector = ArborX::Details::Vector<3, float>;
  BOOST_TEST((Vector{3, 4}.norm()) == 5);
  BOOST_TEST((Vector{6, 13, 18}.norm()) == 23);

  BOOST_TEST((Vector{1, 1e-7, 0}.norm()) == 1);
  BOOST_TEST((Vector{1, 1e-7, 0}.norm<double>()) > 1);
}

BOOST_AUTO_TEST_CASE(vector_normalize)
{
  using Vector = ArborX::Details::Vector<2, float>;
  using ArborX::Details::normalize;

  Vector v{3, 0};
  BOOST_TEST((normalize(v) == Vector{1, 0}));

  Vector w{3, 4};
  BOOST_TEST((normalize<double>(w) == Vector{0.6f, 0.8f}));
}

BOOST_AUTO_TEST_CASE(vector_cross_product)
{
  using Vector = ArborX::Details::Vector<3, float>;

  // clang-format off
  static_assert(Vector{1, 0, 0}.cross(Vector{1, 0, 0}) == Vector{0, 0, 0});
  static_assert(Vector{1, 0, 0}.cross(Vector{0, 1, 0}) == Vector{0, 0, 1});
  static_assert(Vector{1, 0, 0}.cross(Vector{0, 0, 1}) == Vector{0, -1, 0});
  static_assert(Vector{0, 1, 0}.cross(Vector{1, 0, 0}) == Vector{0, 0, -1});
  static_assert(Vector{0, 1, 0}.cross(Vector{0, 1, 0}) == Vector{0, 0, 0});
  static_assert(Vector{0, 1, 0}.cross(Vector{0, 0, 1}) == Vector{1, 0, 0});
  static_assert(Vector{1, 1, 1}.cross(Vector{1, 1, 1}) == Vector{0, 0, 0});
  // clang-format on
}
#undef static_assert