- added an optional algorithm-switch to lsolve().

  Ayeee, this breaks the interface again!

diff --git a/ginac/inifcns.cpp b/ginac/inifcns.cpp
index 3258182e3a65aeb163e0feb18fbdd785b263bfb7..5ebd3d02139b153efa647263a232b956fb7bf5cc 100644 (file)
--- a/ginac/inifcns.cpp
+++ b/ginac/inifcns.cpp
@@ -433,7 +433,7 @@ REGISTER_FUNCTION(Order, eval_func(Order_eval).
 // Solve linear system
 //////////
 
 // 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)) {
 {
        // 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 {
        
        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
        } 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 0defba9fbca56f41354f7c0abec922d4f175dea7..e0ddc8bbe6ef0d885cdf614a49f757b3a1771ab9 100644 (file)
--- 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)
 
 /** 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)
 
 /** Check whether a function is the Order (O(n)) function. */
 inline bool is_order_function(const ex & e)