21 lines
381 B
C
21 lines
381 B
C
/* { dg-options "-O1 -floop-parallelize-all" } */
|
|
|
|
double lagrange(const double x[],
|
|
const double y[],
|
|
long n,
|
|
double xval)
|
|
{
|
|
long i, j;
|
|
double yval = 0.;
|
|
|
|
for( i=0; i < n; i++ )
|
|
{
|
|
double l = 1.;
|
|
for( j=0; j < n; j++ )
|
|
if( i != j )
|
|
l *= (xval-x[j])/(x[i]-x[j]);
|
|
yval += y[i]*l;
|
|
}
|
|
return yval;
|
|
}
|