1  /* { dg-options "-O1 -floop-parallelize-all" } */
       2  
       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  }