From dc2510946d9ce577aab2bd3e5d2f62c50d3faa30 Mon Sep 17 00:00:00 2001 From: Richard Kreckel Date: Tue, 14 Aug 2001 18:13:16 +0000 Subject: [PATCH] - added an optional algorithm-switch to lsolve(). Ayeee, this breaks the interface again! --- ginac/inifcns.cpp | 4 ++-- ginac/inifcns.h | 2 +- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/ginac/inifcns.cpp b/ginac/inifcns.cpp index 3258182e..5ebd3d02 100644 --- a/ginac/inifcns.cpp +++ b/ginac/inifcns.cpp @@ -433,7 +433,7 @@ REGISTER_FUNCTION(Order, eval_func(Order_eval). // Solve linear system ////////// -ex lsolve(const ex &eqns, const ex &symbols) +ex lsolve(const ex &eqns, const ex &symbols, unsigned options) { // solve a system of linear equations if (eqns.info(info_flags::relation_equal)) { @@ -493,7 +493,7 @@ ex lsolve(const ex &eqns, const ex &symbols) matrix solution; try { - solution = sys.solve(vars,rhs); + solution = sys.solve(vars,rhs,options); } catch (const std::runtime_error & e) { // Probably singular matrix or otherwise overdetermined system: // It is consistent to return an empty list diff --git a/ginac/inifcns.h b/ginac/inifcns.h index 0defba9f..e0ddc8bb 100644 --- a/ginac/inifcns.h +++ b/ginac/inifcns.h @@ -132,7 +132,7 @@ DECLARE_FUNCTION_2P(binomial) /** Order term function (for truncated power series). */ DECLARE_FUNCTION_1P(Order) -ex lsolve(const ex &eqns, const ex &symbols); +ex lsolve(const ex &eqns, const ex &symbols, unsigned options = determinant_algo::automatic); /** Check whether a function is the Order (O(n)) function. */ inline bool is_order_function(const ex & e) -- 2.44.0