File: pr42326.c

package info (click to toggle)
gcc-arm-none-eabi 15%3A7-2018-q2-5
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 539,188 kB
  • sloc: ansic: 2,739,935; cpp: 848,238; ada: 602,637; makefile: 62,919; asm: 55,635; xml: 46,238; exp: 23,020; sh: 19,616; python: 6,371; pascal: 3,889; awk: 3,278; perl: 2,691; yacc: 316; ml: 285; f90: 234; lex: 198; objc: 194; haskell: 119
file content (20 lines) | stat: -rw-r--r-- 381 bytes parent folder | download | duplicates (10) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
/* { 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;
}