Forrest-Tomlin

src/basis_factorization/.gitignore

@@ -0,0 +1 @@
+*.cxx

src/basis_factorization/AlmostIdentityMatrix.h

@@ -18,19 +18,29 @@
 class AlmostIdentityMatrix
 {
 public:
+    unsigned _row;
+    unsigned _column;
+    double _value;
+
     AlmostIdentityMatrix()
-        : _identity( true )
     {
     }
 
-    /*
-      True iff this is the identity matrix.
-    */
-    bool _identity;
+    AlmostIdentityMatrix( const AlmostIdentityMatrix &other )
+        : _row( other._row )
+        , _column( other._column )
+        , _value( other._value )
+    {
+    }
+
+    AlmostIdentityMatrix &operator=( const AlmostIdentityMatrix &other )
+    {
+        _row = other._row;
+        _column = other._column;
+        _value = other._value;
 
-    unsigned _row;
-    unsigned _column;
-    double _value;
+        return *this;
+    }
 };
 
 #endif // __AlmostIdentityMatrix_h__

src/basis_factorization/BasisFactorizationError.h

@@ -22,6 +22,8 @@ class BasisFactorizationError : public Error
         ALLOCATION_FAILED = 0,
         CANT_INVERT_BASIS_BECAUSE_OF_ETAS = 1,
         UNKNOWN_BASIS_FACTORIZATION_TYPE = 2,
+        CORRUPT_PERMUATION_MATRIX = 3,
+        CANT_INVERT_BASIS_BECAUSE_BASIS_ISNT_AVAILABLE = 4,
     };
 
     BasisFactorizationError( BasisFactorizationError::Code code ) :

src/basis_factorization/EtaMatrix.cpp

@@ -85,7 +85,7 @@ EtaMatrix::~EtaMatrix()
     }
 }
 
-void EtaMatrix::dump()
+void EtaMatrix::dump() const
 {
     printf( "Dumping eta matrix\n" );
     printf( "\tm = %u. column index = %u\n", _m, _columnIndex );

src/basis_factorization/EtaMatrix.h

@@ -27,7 +27,7 @@ class EtaMatrix
     EtaMatrix &operator=( const EtaMatrix &other );
 
     ~EtaMatrix();
-    void dump();
+    void dump() const;
 	void toMatrix( double *A );
 
     void resetToIdentity();