quadric inversion

55f2deea · Philip Trettner · e208eb52 · 55f2deea
Commit 55f2deea authored Mar 16, 2021 by Philip Trettner
Hide whitespace changes
Inline Side-by-side

Showing with 16 additions and 0 deletions

src/typed-geometry/types/quadric.hh src/typed-geometry/types/quadric.hh +16 -0

No files found.
--- a/src/typed-geometry/types/quadric.hh
+++ b/src/typed-geometry/types/quadric.hh
@@ -314,6 +314,22 @@ public:
               + c;                                   // + c
    }

+    [[nodiscard]] constexpr quadric operator-() const
+    {
+        quadric q;
+        q.A00 = -A00;
+        q.A11 = -A11;
+        q.A22 = -A22;
+        q.A01 = -A01;
+        q.A02 = -A02;
+        q.A12 = -A12;
+        q.b0 = -b0;
+        q.b1 = -b1;
+        q.b2 = -b2;
+        q.c = -c;
+        return q;
+    }
+
    // TODO: fix me
    template <class T>
    friend constexpr T distance_sqr(pos<3, T> const& p, quadric<3, T> const& q);