tests/c_files/gauss.c

   1 int n;
   2
   3 void main(void)
   4 {
   5   int i,j,k,l;
   6   double **a,s;
   7
   8 #pragma scop
   9   for(k = 1; k <= n; k++)
  10     {
  11       s = 1/a[k][k];
  12
  13       for(l = k+1; l <= n; l++)
  14         a[l][k] = a[l][k]*s;
  15
  16       for(j = k+1; j <= n+1; j++)
  17         for(i = k+1; i <= n; i++)
  18           a[i][j] = a[i][j] - a[k][j] * a[i][k];
  19     }
  20 #pragma endscop
  21 }