From 5226e19e3bf2e082b02dec9e856e2613744c8571 Mon Sep 17 00:00:00 2001 From: Tavian Barnes Date: Mon, 1 Nov 2010 02:49:52 -0400 Subject: More numerical tweaks. --- libdimension/polynomial.c | 11 +++++++++++ 1 file changed, 11 insertions(+) (limited to 'libdimension/polynomial.c') diff --git a/libdimension/polynomial.c b/libdimension/polynomial.c index 6e2aed2..ab80f69 100644 --- a/libdimension/polynomial.c +++ b/libdimension/polynomial.c @@ -57,6 +57,15 @@ dmnsn_real_degree(const double poly[], size_t degree) return 0; } +/* Divide each coefficient by the leading coefficient */ +static inline void +dmnsn_normalize_polynomial(double poly[], size_t degree) +{ + for (size_t i = 0; i <= degree; ++i) { + poly[i] /= poly[degree]; + } +} + /* Eliminate trivial zero roots from poly[] */ static inline void dmnsn_eliminate_zero_roots(double poly[], size_t *degree) @@ -285,6 +294,8 @@ dmnsn_solve_polynomial(const double poly[], size_t degree, double x[]) /* Account for leading zero coefficients */ degree = dmnsn_real_degree(p, degree); + /* Normalize the leading coefficient to 1.0 */ + dmnsn_normalize_polynomial(p, degree); /* Eliminate simple zero roots */ dmnsn_eliminate_zero_roots(p, °ree); -- cgit v1.2.3