From 0cc840f49d87e01b31ad5a848a8bf5470ea43275 Mon Sep 17 00:00:00 2001 From: Alexander Frink Date: Mon, 21 Feb 2000 23:02:43 +0000 Subject: [PATCH] added a normal() to compare the output from lsolve with what is expected --- check/linear_solve.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/check/linear_solve.cpp b/check/linear_solve.cpp index 0faabb0b..1a5a0bd3 100644 --- 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; -- 2.44.0