1 /* { dg-options "-O1 -floop-parallelize-all" } */
2 
lagrange(const double x[],const double y[],long n,double xval)3 double lagrange(const double x[],
4                 const double y[],
5                 long n,
6                 double xval)
7 {
8   long i, j;
9   double yval = 0.;
10 
11   for( i=0; i < n; i++ )
12     {
13       double l = 1.;
14       for( j=0; j < n; j++ )
15 	if( i != j )
16 	  l *= (xval-x[j])/(x[i]-x[j]);
17       yval += y[i]*l;
18     }
19   return yval;
20 }
21