added a normal() to compare the output from lsolve with what is expected

diff --git a/check/linear_solve.cpp b/check/linear_solve.cpp
index 0faabb0b6ad7b0a750576d3f4dbacebf3142e8d9..1a5a0bd39aa7fec33f497c28f2dce1711a579bab 100644 (file)
--- a/check/linear_solve.cpp
+++ b/check/linear_solve.cpp
@@ -61,8 +61,8 @@ static unsigned lsolve2a(void)
     ex sol_y = sol.op(1).rhs();  // rhs of solution for second variable (y)
     
     // It should have returned [x==(3+b^2)/(a+b),y==(3-a*b)/(a+b)]
-    if (!(sol_x - (3+pow(b,2))/(a+b)).is_zero() ||
-        !(sol_y - (3-a*b)/(a+b)).is_zero()) {
+    if (!normal(sol_x - (3+pow(b,2))/(a+b)).is_zero() ||
+        !normal(sol_y - (3-a*b)/(a+b)).is_zero()) {
         result++;
         clog << "solution of the system " << eqns << " for " << vars
              << " erroneously returned " << sol << endl;