From 80388b952b5844261e711e8da5a59c34b4919da7 Mon Sep 17 00:00:00 2001 From: Tavian Barnes Date: Sun, 31 Oct 2010 23:28:14 -0400 Subject: Numerical fixes for polynomial.c. --- bench/libdimension/polynomial.c | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'bench') diff --git a/bench/libdimension/polynomial.c b/bench/libdimension/polynomial.c index bcf9281..eaf7efe 100644 --- a/bench/libdimension/polynomial.c +++ b/bench/libdimension/polynomial.c @@ -27,29 +27,29 @@ main(void) #define NPOLY 5 double p[NPOLY][NPOLY + 1], x[NPOLY]; - // p[0][] = x - 0.5; + /* p[0][] = x - 0.5; */ p[0][1] = 1.0; p[0][0] = -0.5; - // p[1][] = (x + 0.5)*(x - 0.5) + /* p[1][] = (x + 0.5)*(x - 0.5) */ p[1][2] = 1.0; p[1][1] = 0.0; p[1][0] = -0.25; - // p[2][] = (x + 1)*(x - 1.2345)*(x - 100) + /* p[2][] = (x + 1)*(x - 1.2345)*(x - 100) */ p[2][3] = 1.0; p[2][2] = -100.2345; p[2][1] = 22.2155; p[2][0] = 123.45; - // p[3][] = (x + 1)*(x - 1.2345)*(x - 5)*(x - 100) + /* p[3][] = (x + 1)*(x - 1.2345)*(x - 5)*(x - 100) */ p[3][4] = 1.0; p[3][3] = -105.2345; p[3][2] = 523.388; p[3][1] = 12.3725; p[3][0] = -617.25; - // p[4][] = (x + 1)*(x - 1.2345)*(x - 2.3456)*(x - 5)*(x - 100) + /* p[4][] = (x + 1)*(x - 1.2345)*(x - 2.3456)*(x - 5)*(x - 100) */ p[4][5] = 1.0; p[4][4] = -107.5801; p[4][3] = 770.2260432; -- cgit v1.2.3