copy variants for sdf3

Author: MarcG-TUE

include/maxplus/algebra/mpmatrix.h

@@ -171,7 +171,8 @@ class Matrix {
 
     void pasteRowVector(unsigned int top_row, unsigned int left_column, const Vector *pastedVector);
 
-    [[nodiscard]] virtual std::shared_ptr<Matrix> createCopy() const;
+    [[nodiscard]] virtual std::shared_ptr<Matrix> createCopyPtr() const;
+    [[nodiscard]] virtual Matrix createCopy() const;
 
     [[nodiscard]] virtual std::shared_ptr<Matrix> getTransposedCopy() const;
 

src/algebra/mpmatrix.cc

@@ -573,9 +573,9 @@ Matrix Matrix::mp_power(const unsigned int p) const {
 }
 
 /**
-    * Matrix copy.
-    */
-std::shared_ptr<Matrix> Matrix::createCopy() const
+ * Matrix copy.
+ */
+std::shared_ptr<Matrix> Matrix::createCopyPtr() const
 {
     std::shared_ptr<Matrix> newMatrix = std::make_shared<Matrix>(this->getRows(), this->getCols());
     unsigned int nEls = this->getRows() * this->getCols();
@@ -584,6 +584,18 @@ std::shared_ptr<Matrix> Matrix::createCopy() const
     }
     return newMatrix;
 }
+/**
+ * Matrix copy.
+ */
+Matrix Matrix::createCopy() const
+{
+    Matrix newMatrix(this->getRows(), this->getCols());
+    unsigned int nEls = this->getRows() * this->getCols();
+    for (unsigned int pos = 0; pos < nEls; pos++) {
+        newMatrix.table[pos] = this->table[pos];
+    }
+    return newMatrix;
+}
 
 /**
     * Matrix transposed copy.

src/graph/smpls.cc

@@ -473,7 +473,7 @@ namespace MaxPlus::SMPLS {
 					throw MPException("mode " + output + " not found in dissected matrices!");
 				}
 
-				sMatrix = disSm->core.begin()->second->createCopy();
+				sMatrix = disSm->core.begin()->second->createCopyPtr();
 
 				// look if it has any events to emit
 				if (!disSm->eventRows.empty())