define multivector negation

Change-Id: Ia277a4e73b42baac8295cd1dfc74061a6a150b97
oliverlee · Jan 12, 2024 · e1c9127 · e1c9127
1 parent 8f806d6
commit e1c9127
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 1 deletion.
diff --git a/geometry/src/algebra.hpp b/geometry/src/algebra.hpp
@@ -284,7 +284,7 @@ struct algebra
 
     /// construct a multivector from blades
     ///
-    constexpr multivector(Bs... bs) : Bs{bs}... {}
+    constexpr explicit multivector(Bs... bs) : Bs{bs}... {}
 
     /// equality comparison
     ///
@@ -302,6 +302,15 @@ struct algebra
       return not(x == y);
     }
     /// @}
+
+    /// unary negation
+    ///
+    [[nodiscard]]
+    friend constexpr auto
+    operator-(const multivector& x) -> multivector
+    {
+      return multivector{-static_cast<Bs>(x)...};
+    }
   };
 };
 

diff --git a/test/algebra_multivector_test.cpp b/test/algebra_multivector_test.cpp
@@ -51,4 +51,11 @@ auto main() -> int
             eq(multivector(2 * e<0>, 3 * e<1>), 3 * e<1> + 2 * e<0>) and
             eq(multivector(-2 * e<0>, 3 * e<1>), 3 * e<1> - 2 * e<0>));
       };
+
+  "multivector negation"_ctest =  //
+      [] {
+        return expect(
+            eq(e<0> + e<1>, -(-e<0> - e<1>)) and
+            eq(e<0, 1> + e<1, 2>, -(-e<0, 1> - e<1, 2>)));
+      };
 }