diff options
| author | Tavian Barnes <tavianator@gmail.com> | 2010-10-31 23:28:14 -0400 |
|---|---|---|
| committer | Tavian Barnes <tavianator@gmail.com> | 2010-10-31 23:28:14 -0400 |
| commit | 80388b952b5844261e711e8da5a59c34b4919da7 (patch) | |
| tree | 22d7b414b392cea5003029e5360efb3b4b30c78f /bench | |
| parent | 70ddf1e89aca531b5a2b4d71d2873c6883a61953 (diff) | |
| download | dimension-80388b952b5844261e711e8da5a59c34b4919da7.tar.xz | |
Numerical fixes for polynomial.c.
Diffstat (limited to 'bench')
| -rw-r--r-- | bench/libdimension/polynomial.c | 10 |
1 files changed, 5 insertions, 5 deletions
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; |