Insert more code.

MathieuDutSik · Feb 28, 2024 · 8593a87 · 8593a87
1 parent fc45fc9
commit 8593a87
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src_matrix/MAT_MatrixInt.h b/src_matrix/MAT_MatrixInt.h
@@ -234,6 +234,15 @@ GCD_dot<T> ComputeGcdDot(MyVector<T> const& x) {
   return {std::move(V), gcd};
 }
 
+template<typename T>
+GCD_dot<T> PositivityNormalizeGcdDot(GCD_dot<T> const& x) {
+  if (x.gcd > 0) {
+    return x;
+  }
+  return {-x.V, -x.gcd};
+}
+
+
 // This is the return type for the second extended GCD computations
 // In input a list of entries x=(x_0, ...., x_m)
 // In return we have: