* Interface to symbolic matrices */
/*
- * GiNaC Copyright (C) 1999-2015 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2019 Johannes Gutenberg University Mainz, Germany
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
matrix(unsigned r, unsigned c, const lst & l);
matrix(std::initializer_list<std::initializer_list<ex>> l);
- matrix_init<ex, exvector::iterator> operator=(const ex & x) deprecated;
+ matrix_init<ex, exvector::iterator> operator=(const ex & x) attribute_deprecated;
protected:
matrix(unsigned r, unsigned c, const exvector & m2);
matrix(unsigned r, unsigned c, exvector && m2);
size_t nops() const override;
ex op(size_t i) const override;
ex & let_op(size_t i) override;
- ex eval(int level=0) const override;
ex evalm() const override {return *this;}
ex subs(const exmap & m, unsigned options = 0) const override;
ex eval_indexed(const basic & i) const override;
ex trace() const;
ex charpoly(const ex & lambda) const;
matrix inverse() const;
+ matrix inverse(unsigned algo) const;
matrix solve(const matrix & vars, const matrix & rhs,
unsigned algo = solve_algo::automatic) const;
unsigned rank() const;
+ unsigned rank(unsigned solve_algo) const;
bool is_zero_matrix() const;
protected:
ex determinant_minor() const;
+ std::vector<unsigned> echelon_form(unsigned algo, int n);
int gauss_elimination(const bool det = false);
int division_free_elimination(const bool det = false);
int fraction_free_elimination(const bool det = false);
+ std::vector<unsigned> markowitz_elimination(unsigned n);
int pivot(unsigned ro, unsigned co, bool symbolic = true);
void print_elements(const print_context & c, const char *row_start, const char *row_end, const char *row_sep, const char *col_sep) const;
inline ex expand(const matrix & m, unsigned options = 0)
{ return m.expand(options); }
-inline ex eval(const matrix & m, int level = 0)
-{ return m.eval(level); }
-
-inline ex evalf(const matrix & m, int level = 0)
-{ return m.evalf(level); }
+inline ex evalf(const matrix & m)
+{ return m.evalf(); }
inline unsigned rows(const matrix & m)
{ return m.rows(); }
{ return m.charpoly(lambda); }
inline matrix inverse(const matrix & m)
-{ return m.inverse(); }
+{ return m.inverse(solve_algo::automatic); }
+inline matrix inverse(const matrix & m, unsigned algo)
+{ return m.inverse(algo); }
inline unsigned rank(const matrix & m)
{ return m.rank(); }
+inline unsigned rank(const matrix & m, unsigned solve_algo)
+{ return m.rank(solve_algo); }
// utility functions