File:
[DVB] /
dietlibc /
libm /
poly.c
Revision
1.2:
download - view:
text,
annotated -
select for diffs
Mon Nov 18 01:16:51 2002 UTC (21 years, 6 months ago) by
fefe
Branches:
MAIN
CVS tags:
finnland_test_200301,
branch_rc14_fieldtest_finnland,
branch_rc13_fieldtest_finnland,
branch_rc12_fieldtest_finnland,
branch_rc10_fieldtest_finnland,
RELEASE_finnland_200301_1,
RC12_FIELDTEST_FINNLAND,
RC10_FIELDTEST_FINNLAND,
HEAD
remove a few more warnings
/*--------------------------------------------------------------------------*
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;
}
LinuxTV legacy CVS <linuxtv.org/cvs>