Update makefile ad update basic_vec to allow abstract algebraic structures

Unfortunately this requires C++17 as `static constexpr` members
have an undefined behaviour before this version.

1 files changed, 9 insertions, 4 deletions

diff --git a/vector/vector.cpp b/vector/vector.cpp
index 9e0ca6c..d5f0c36 100644
--- a/vector/vector.cpp
+++ b/vector/vector.cpp
@@ -6,10 +6,15 @@
 #include <initializer_list>
 #include <algorithm>
 
+
 template<typename T, std::size_t d>
 struct basic_vec : public std::array<T, d> {
     static constexpr std::size_t dimensions = d;
 
+    static constexpr T null_element = static_cast<T>(0);
+    static constexpr T unit_element = static_cast<T>(1);
+    static constexpr T unit_additive_inverse_element = static_cast<T>(-1);
+
     basic_vec();
     basic_vec(const std::initializer_list<T> l);
     template<std::size_t n> basic_vec(const basic_vec<T, n>& other);
@@ -19,7 +24,7 @@ struct basic_vec : public std::array<T, d> {
 
 template<typename T, std::size_t d>
 basic_vec<T, d>::basic_vec() : std::array<T, d>() {
-    this->fill(static_cast<T>(0));
+    this->fill(basic_vec<T,d>::null_element);
 }
 
 template<typename T, std::size_t d>
@@ -44,7 +49,7 @@ basic_vec<T, d>::basic_vec(const basic_vec<T, n>& other) {
 
 template<typename T, std::size_t d>
 T basic_vec<T, d>::length() const {
-    T res = static_cast<T>(0);
+    T res = basic_vec<T, d>::null_element;
     for (const T& val : *this) {
         res += val * val;
     }
@@ -77,12 +82,12 @@ basic_vec<T, d> operator*(const T& rhs, const basic_vec<T, d>& lhs) {
 
 template<typename T, std::size_t d>
 basic_vec<T, d> operator-(const basic_vec<T, d>& rhs, const basic_vec<T, d>& lhs) {
-    return rhs + static_cast<T>(-1) * lhs;
+    return rhs + basic_vec<T, d>::unit_additive_inverse_element * lhs;
 }
 
 template<typename T, std::size_t d>
 T operator*(const basic_vec<T, d>& rhs, const basic_vec<T, d>& lhs) {
-    T res = static_cast<T>(0);
+    T res = basic_vec<T, d>::null_element;
     for (std::size_t i = 0; i < d; i++) {
         res += rhs[i] * lhs[i];
     }