]> www.ginac.de Git - ginac.git/blobdiff - ginac/matrix.h
- expairseq.cpp: moved expairseq::to_rational to...
[ginac.git] / ginac / matrix.h
index c0e35dcdaf5e7c59e648f694f6d3756cd5fb4f74..a36276f4d332b06f57afdc6747e7a5f54697cf72 100644 (file)
@@ -89,7 +89,7 @@ public:
     matrix transpose(void) const;
     ex determinant(void) const;
     ex trace(void) const;
-    ex charpoly(const ex & lambda) const;
+    ex charpoly(const symbol & lambda) const;
     matrix inverse(void) const;
     matrix fraction_free_elim(const matrix & vars, const matrix & v) const;
     matrix solve(const matrix & vars, const matrix & rhs) const;
@@ -97,9 +97,10 @@ public:
 protected:
     ex determinant_numeric(void) const;
     ex determinant_minor(void) const;
-    ex determinant_perm(void) const;
+    ex determinant_bareiss(bool normalize=true) const;
     int gauss_elimination(void);
     int fraction_free_elimination(void);
+    int division_free_elimination(void);
     int pivot(unsigned ro, bool symbolic=true);
 private:  // FIXME: these should be obsoleted
     void ffe_swap(unsigned r1, unsigned c1, unsigned r2 ,unsigned c2);
@@ -150,7 +151,7 @@ inline ex determinant(const matrix & m)
 inline ex trace(const matrix & m)
 { return m.trace(); }
 
-inline ex charpoly(const matrix & m, const ex & lambda)
+inline ex charpoly(const matrix & m, const symbol & lambda)
 { return m.charpoly(lambda); }
 
 inline matrix inverse(const matrix & m)
@@ -162,6 +163,8 @@ inline const matrix &ex_to_matrix(const ex &e)
        return static_cast<const matrix &>(*e.bp);
 }
 
+extern ex lst_to_matrix(const ex &l);
+
 #ifndef NO_NAMESPACE_GINAC
 } // namespace GiNaC
 #endif // ndef NO_NAMESPACE_GINAC