examples/gras/chrono/chrono.c

   1 /* $Id$ */
   2
   3 /* chrono - demo of GRAS benchmarking features                              */
   4
   5 /* Copyright (c) 2005 Martin Quinson, Arnaud Legrand. All rights reserved.  */
   6
   7 /* This program is free software; you can redistribute it and/or modify it
   8  * under the terms of the license (GNU LGPL) which comes with this package. */
   9
  10 #include "gras.h"
  11 #include "xbt/log.h"
  12
  13 XBT_LOG_NEW_DEFAULT_CATEGORY(Chrono,"Messages specific to this example");
  14
  15
  16 /* Function prototypes */
  17 int multiplier (int argc,char *argv[]);
  18
  19 int multiplier (int argc,char *argv[])
  20 {
  21   int i,j,k,l;
  22   double *A,*B,*C;
  23   int n = 100;
  24   double start = 0.0;
  25   double now = 0.0;
  26
  27   gras_init(&argc, argv, NULL);
  28
  29   A = malloc(n*n*sizeof(double));
  30   B = malloc(n*n*sizeof(double));
  31   C = malloc(n*n*sizeof(double));
  32
  33   start=now=gras_os_time();
  34
  35   INFO1("Before computation: %lg", start);
  36
  37   for(l=0; l<4; l++) {
  38     now=gras_os_time();
  39     GRAS_BENCH_ONCE_RUN_ONCE_BEGIN();
  40     for(i=0; i<n; i++)
  41       for(j=0; j<n; j++) {
  42         A[i*n+j]=2/n;
  43         B[i*n+j]=1/n;
  44         C[i*n+j]=0.0;
  45       }
  46
  47     for(i=0; i<n; i++)
  48       for(j=0; j<n; j++)
  49         for(k=0; k<n; k++)
  50           C[i*n+j] += A[i*n+k]*B[k*n+j];
  51
  52     GRAS_BENCH_ONCE_RUN_ONCE_END();
  53     now=gras_os_time()-now;
  54     INFO2("Iteration %d : %lg ", l, now);
  55   }
  56
  57   now=gras_os_time()-start;
  58   INFO2("After computation: %lg; Duration: %lg ", gras_os_time(), now);
  59
  60   return 0;
  61 }