- added an optional algorithm-switch to lsolve().
authorRichard Kreckel <Richard.Kreckel@uni-mainz.de>
Tue, 14 Aug 2001 18:13:16 +0000 (18:13 +0000)
committerRichard Kreckel <Richard.Kreckel@uni-mainz.de>
Tue, 14 Aug 2001 18:13:16 +0000 (18:13 +0000)
  Ayeee, this breaks the interface again!

ginac/inifcns.cpp
ginac/inifcns.h

index 3258182e3a65aeb163e0feb18fbdd785b263bfb7..5ebd3d02139b153efa647263a232b956fb7bf5cc 100644 (file)
@@ -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
index 0defba9fbca56f41354f7c0abec922d4f175dea7..e0ddc8bbe6ef0d885cdf614a49f757b3a1771ab9 100644 (file)
@@ -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)