Insert one more function.

MathieuDutSik · Jul 22, 2024 · 76f3d57 · 76f3d57
1 parent 75eee0f
commit 76f3d57
Showing 1 changed file with 19 additions and 0 deletions.
diff --git a/src_matrix/MAT_Matrix.h b/src_matrix/MAT_Matrix.h
@@ -2909,6 +2909,25 @@ MyMatrix<T> MatrixFromVectorFamily(std::vector<MyVector<T>> const &ListVect) {
   return M;
 }
 
+template<typename T>
+MyMatrix<T> MatrixFromPairVector(std::pair<MyVector<T>, MyVector<T>> const& pair) {
+  int n1 = pair.first.size();
+  int n2 = pair.second.size();
+  if (n1 != n2) {
+    std::cerr << "n1=" << n1 << " n2=" << n2 << "\n";
+    std::cerr << "Vector lengths are not homogeneous\n";
+    throw TerminalException{1};
+  }
+  int n = n1;
+  MyMatrix<T> M(2,n);
+  for (int i=0; i<n; i++) {
+    M(0, i) = pair.first(i);
+    M(1, i) = pair.second(i);
+  }
+  return M;
+}
+
+
 template <typename T>
 MyMatrix<T>
 MatrixFromVectorFamilyDim(int const &dim,