blob: cdcfb8c5caf77dfcb4c84850dc54c2762e781c36 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
/*--------------------------------------------------------------------------*
Name __poly - generates a polynomial from arguments
Usage double __poly ( double x, int n, const double* c );
Prototype in math.h
Description __poly generates a polynomial in x, of degree n, with
coefficients c[0], c[1], ..., c[n]. For example, if n=4,
the generated polynomial is
c[4]*x^4 + c[3]*x^3 + c[2]*x^2 + c[1]*x + c[0]
The polynomial is calculated using Horner's method:
polynom = (..((x*c[n] + c[n-1])*x + c[n-2])..)*x + c[0]
Return value __poly returns the value of the polynomial as evaluated for
the given x.
A range error occurs if the result exceeds double range.
*---------------------------------------------------------------------------*/
#include <stdio.h>
#include "dietlibm.h"
double __poly ( double x, size_t n, const double* c)
{
long double ret;
size_t i;
i = n;
c += n;
ret = 0;
do
ret = ret * x + *c--;
while ( i-- );
return ret;
}
|