examples/trisolv/trisolv.c

   1 #include <stdio.h>
   2 #include <stdlib.h>
   3 #include <math.h>
   4 #include <sys/time.h>
   5
   6 #include "decls.h"
   7 #include "util.h"
   8
   9 void trisolv(long N)
  10 {
  11     long i,j,k;
  12
  13 #pragma scop
  14     for (i=0;i<=N-1;i++) {
  15         for (j=0;j<=N-1;j++) {
  16             for (k=0;k<=j-1;k++) {
  17                 B[j][i]=B[j][i]-L[j][k]*B[k][i];  //S1 ;
  18
  19             }
  20             B[j][i]=B[j][i]/L[j][j]; // S2 ;
  21         } // for j
  22     } // for i
  23 #pragma endscop
  24 }
  25
  26
  27 int main()
  28 {
  29     long N=NMAX;
  30     int i,j;
  31     double t_start, t_end;
  32
  33     IF_TIME(t_start = rtclock());
  34     trisolv(N);
  35     IF_TIME(t_end = rtclock());
  36     IF_TIME(fprintf(stderr, "%0.6lfs\n", t_end - t_start));
  37
  38     if (fopen(".test", "r")) {
  39         for (i = 0; i < N; i++) {
  40             for (j = 0; j < N; j++) {
  41                 fprintf(stdout, "%lf ", B[i][j]);
  42             }
  43             fprintf(stdout, "\n");
  44         }
  45     }
  46     return 0;
  47 }