1 /* { dg-options "-O1 -floop-parallelize-all" } */ 2 lagrange(const double x[],const double y[],long n,double xval)3double 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