From 32ad2ae1a0e96703553e3cf1fc265c8ba3757562 Mon Sep 17 00:00:00 2001
From: sjunges <sebastian.junges@rwth-aachen.de>
Date: Mon, 24 Mar 2014 13:34:11 +0100
Subject: [PATCH] minimal changes to make sparse matrices with polynomials work

Former-commit-id: 6688ebc9fdfae1434eb1a00fa4cf02b798d22a9d
---
 src/storage/SparseMatrix.cpp | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)

diff --git a/src/storage/SparseMatrix.cpp b/src/storage/SparseMatrix.cpp
index 0f8d0acfa..21814b643 100644
--- a/src/storage/SparseMatrix.cpp
+++ b/src/storage/SparseMatrix.cpp
@@ -8,6 +8,8 @@
 #endif
 
 #include "src/storage/SparseMatrix.h"
+#include "src/storage/parameters.h"
+
 #include "src/exceptions/InvalidStateException.h"
 
 #include "log4cplus/logger.h"
@@ -692,6 +694,14 @@ namespace storm {
             return std::make_pair(std::move(resultLU), dInvBuilder.build());
         }
         
+        template<>
+        typename std::pair<storm::storage::SparseMatrix<Polynomial>, storm::storage::SparseMatrix<Polynomial>> SparseMatrix<Polynomial>::getJacobiDecomposition() const {
+            // NOT SUPPORTED
+            // TODO do whatever storm does in such cases.
+            assert(false);
+        }
+        
+        
         template<typename T>
         std::vector<T> SparseMatrix<T>::getPointwiseProductRowSumVector(storm::storage::SparseMatrix<T> const& otherMatrix) const {
             std::vector<T> result(rowCount, storm::utility::constantZero<T>());
@@ -923,6 +933,12 @@ namespace storm {
         template class SparseMatrixBuilder<int>;
         template class SparseMatrix<int>;
         template std::ostream& operator<<(std::ostream& out, SparseMatrix<int> const& matrix);
+#ifdef STORM_HAVE_CARL
+        template class SparseMatrixBuilder<Polynomial>;
+        template class SparseMatrix<Polynomial>;
+        template std::ostream& operator<<(std::ostream& out, SparseMatrix<Polynomial> const& matrix);
+#endif
+        
         
     } // namespace storage
 } // namespace storm