| /* { 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; | |
| } |