|
| 1 | +#ifndef CPPID_GM3D_DETAIL_BASIC_POINT_ND_HPP |
| 2 | +#define CPPID_GM3D_DETAIL_BASIC_POINT_ND_HPP |
| 3 | + |
| 4 | +#include <algorithm> |
| 5 | +#include <cmath> |
| 6 | +#include <cstddef> |
| 7 | +#include <functional> |
| 8 | +#include <numeric> |
| 9 | +#include <type_traits> |
| 10 | +#include <utility> |
| 11 | + |
| 12 | +#include <boost/qvm/all.hpp> |
| 13 | + |
| 14 | +#include "basic_point_vector.hpp" |
| 15 | +#include "type_traits.hpp" |
| 16 | + |
| 17 | +namespace cppid::gm3d::detail { |
| 18 | + |
| 19 | +template<typename T, std::size_t N> |
| 20 | +struct basic_point_nd : public basic_point_vector<basic_point_nd<T, N>, T, N> { |
| 21 | + using typename basic_point_vector<basic_point_nd<T, N>, T, N>::scalar_type; |
| 22 | + |
| 23 | + using basic_point_vector<basic_point_nd<T, N>, T, N>::dimensions; |
| 24 | + |
| 25 | + using basic_point_vector<basic_point_nd<T, N>, T, N>::basic_point_vector; |
| 26 | + |
| 27 | + constexpr basic_point_nd() noexcept = default; |
| 28 | + |
| 29 | + template<typename U, |
| 30 | + std::enable_if_t< |
| 31 | + detail::is_explicit_constructible_v<scalar_type, U>>* = nullptr> |
| 32 | + explicit constexpr basic_point_nd(basic_point_nd<U, dimensions> v) noexcept |
| 33 | + : basic_point_nd{v, std::make_index_sequence<dimensions>()} |
| 34 | + { |
| 35 | + } |
| 36 | + |
| 37 | + template<typename U, |
| 38 | + std::size_t... Ns, |
| 39 | + std::enable_if_t< |
| 40 | + (sizeof...(Ns) == dimensions) && |
| 41 | + detail::is_explicit_constructible_v<scalar_type, U>>* = nullptr> |
| 42 | + constexpr basic_point_nd(basic_point_nd<U, dimensions> v, |
| 43 | + std::index_sequence<Ns...>) noexcept |
| 44 | + : basic_point_nd{v.elems[Ns]...} |
| 45 | + { |
| 46 | + } |
| 47 | + |
| 48 | + template<typename U, |
| 49 | + std::enable_if_t< |
| 50 | + detail::is_implicit_constructible_v<scalar_type, U>>* = nullptr> |
| 51 | + constexpr basic_point_nd(basic_point_nd<U, dimensions> v) noexcept |
| 52 | + : basic_point_nd{v, std::make_index_sequence<dimensions>()} |
| 53 | + { |
| 54 | + } |
| 55 | + |
| 56 | + template<typename U, |
| 57 | + std::size_t... Ns, |
| 58 | + std::enable_if_t< |
| 59 | + (sizeof...(Ns) == dimensions) && |
| 60 | + detail::is_implicit_constructible_v<scalar_type, U>>* = nullptr> |
| 61 | + constexpr basic_point_nd(basic_point_nd<U, dimensions> v, |
| 62 | + std::index_sequence<Ns...>) noexcept |
| 63 | + : basic_point_nd{scalar_type(v.elems[Ns])...} |
| 64 | + { |
| 65 | + } |
| 66 | + |
| 67 | + template<typename U, |
| 68 | + typename = std::enable_if_t<std::is_convertible_v<U, scalar_type>>> |
| 69 | + constexpr auto operator=(const basic_point_nd<U, dimensions> rhs) noexcept |
| 70 | + -> basic_point_nd& |
| 71 | + { |
| 72 | + op_assign(rhs.elems, std::make_index_sequence<dimensions>()); |
| 73 | + return *this; |
| 74 | + } |
| 75 | + |
| 76 | + constexpr auto operator==(const basic_point_nd& rhs) const noexcept -> bool |
| 77 | + { |
| 78 | + return is_equal(rhs.elems, std::make_index_sequence<dimensions>()); |
| 79 | + } |
| 80 | + |
| 81 | + constexpr auto operator!=(const basic_point_nd& rhs) const noexcept -> bool |
| 82 | + { |
| 83 | + return !(*this == rhs); |
| 84 | + } |
| 85 | + |
| 86 | + friend constexpr auto operator+(const basic_point_nd& val) noexcept |
| 87 | + -> basic_point_nd |
| 88 | + { |
| 89 | + return val; |
| 90 | + } |
| 91 | + |
| 92 | + friend constexpr auto operator-(const basic_point_nd& lhs) noexcept |
| 93 | + -> basic_point_nd |
| 94 | + { |
| 95 | + return invoke_op{}( |
| 96 | + std::negate{}, lhs, std::make_index_sequence<dimensions>()); |
| 97 | + } |
| 98 | + |
| 99 | +private: |
| 100 | + struct invoke_op { |
| 101 | + template<typename P, typename U, std::size_t... Ns> |
| 102 | + constexpr auto operator()(const P op, |
| 103 | + const basic_point_nd<U, dimensions>& lhs, |
| 104 | + std::index_sequence<Ns...>) const noexcept |
| 105 | + -> basic_point_nd<std::invoke_result_t<P, U>, dimensions> |
| 106 | + { |
| 107 | + static_assert(sizeof...(Ns) == dimensions); |
| 108 | + return {op(lhs.elems[Ns])...}; |
| 109 | + } |
| 110 | + |
| 111 | + template<typename P, typename U, typename V, std::size_t... Ns> |
| 112 | + constexpr auto operator()(const P op, |
| 113 | + const basic_point_nd<U, dimensions>& lhs, |
| 114 | + const basic_point_nd<V, dimensions>& rhs, |
| 115 | + std::index_sequence<Ns...>) const noexcept |
| 116 | + -> basic_point_nd<std::invoke_result_t<P, U, V>, dimensions> |
| 117 | + { |
| 118 | + static_assert(sizeof...(Ns) == dimensions); |
| 119 | + return {op(lhs.elems[Ns], rhs.elems[Ns])...}; |
| 120 | + } |
| 121 | + }; |
| 122 | + |
| 123 | + template<typename U, typename P, std::size_t... Ns> |
| 124 | + inline constexpr void op_assign(const U (&rhs)[dimensions], |
| 125 | + P op, |
| 126 | + std::index_sequence<Ns...>) noexcept |
| 127 | + { |
| 128 | + ((this->elems[Ns] = op(this->elems[Ns], rhs[Ns])), ...); |
| 129 | + } |
| 130 | + |
| 131 | + template<typename U, typename P, std::size_t... Ns> |
| 132 | + inline constexpr void |
| 133 | + op_assign(U rhs, P op, std::index_sequence<Ns...>) noexcept |
| 134 | + { |
| 135 | + ((this->elems[Ns] = op(this->elems[Ns], rhs)), ...); |
| 136 | + } |
| 137 | + |
| 138 | + template<typename U, std::size_t... Ns> |
| 139 | + inline constexpr void op_assign(const U (&rhs)[dimensions], |
| 140 | + std::index_sequence<Ns...>) noexcept |
| 141 | + { |
| 142 | + ((this->elems[Ns] = rhs[Ns]), ...); |
| 143 | + } |
| 144 | + |
| 145 | + template<typename U, std::size_t... Ns> |
| 146 | + constexpr auto is_equal(const U (&rhs)[dimensions], |
| 147 | + std::index_sequence<Ns...>) const noexcept -> bool |
| 148 | + { |
| 149 | + return ((this->elems[Ns] == rhs[Ns]) && ... && true); |
| 150 | + } |
| 151 | + |
| 152 | + template<std::size_t... Ns> |
| 153 | + constexpr auto negate(std::index_sequence<Ns...>) const noexcept |
| 154 | + -> basic_point_nd |
| 155 | + { |
| 156 | + return basic_point_nd{(-this->elems[Ns])...}; |
| 157 | + } |
| 158 | +}; |
| 159 | + |
| 160 | +} // namespace cppid::gm3d::detail |
| 161 | + |
| 162 | +namespace boost::qvm { |
| 163 | + |
| 164 | +template<typename T, std::size_t N> |
| 165 | +struct vec_traits<cppid::gm3d::detail::basic_point_nd<T, N>> |
| 166 | +: public vec_traits< |
| 167 | + cppid::gm3d::detail:: |
| 168 | + basic_point_vector<cppid::gm3d::detail::basic_point_nd<T, N>, T, N>> { |
| 169 | +}; |
| 170 | + |
| 171 | +} // namespace boost::qvm |
| 172 | + |
| 173 | +#endif |
0 commit comments