Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Add generated files to the CVS so that the use of autotools is not necessary to build...
[simgrid.git] / examples / gras / chrono / chrono2.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 #include <cblas.h>
16 void cblas_dgemm(const enum CBLAS_ORDER Order,
17                  const enum CBLAS_TRANSPOSE TransA,
18                  const enum CBLAS_TRANSPOSE TransB, const int M,
19                  const int N, const int K, const double alpha,
20                  const double *A, const int lda, const double *B,
21                  const int ldb, const double beta, double *C,
22                  const int ldc);
23
24
25 /* Function prototypes */
26 static int mat_mult(int n)
27 {
28   int i,j,k,l;
29   double *A,*B,*C;
30   double start = 0.0;
31   double now = 0.0;
32
33   A = malloc(n*n*sizeof(double));
34   B = malloc(n*n*sizeof(double));
35   C = malloc(n*n*sizeof(double));
36
37   start=now=gras_os_time();
38
39   INFO1("Matrix size: %d",n);
40 /*   INFO1("Before computation: %lg", start); */
41
42   for(l=0; l<40; l++) {
43     now=gras_os_time();
44     GRAS_BENCH_ONCE_RUN_ONCE_BEGIN();
45     for(i=0; i<n; i++)
46       for(j=0; j<n; j++) {
47         A[i*n+j]=2/n;
48         B[i*n+j]=1/n;
49         C[i*n+j]=0.0;
50       }
51
52     cblas_dgemm(CblasRowMajor,
53                 CblasNoTrans, CblasNoTrans, n, n, n,
54                 1.0, A, n, B, n, 0.0, C, n);    
55     GRAS_BENCH_ONCE_RUN_ONCE_END();
56     now=gras_os_time()-now;
57 /*     INFO2("Iteration %d : %lg ", l, now); */
58   }
59
60   now=gras_os_time()-start;
61   INFO1("Duration: %lg ", now);
62
63   free(A);
64   free(B);
65   free(C);
66 }
67
68 int multiplier (int argc,char *argv[])
69 {
70   gras_init(&argc, argv);
71
72   mat_mult(10);
73   mat_mult(20);
74   mat_mult(30);
75   mat_mult(40);
76   mat_mult(50);
77   mat_mult(75);
78   mat_mult(100);
79   mat_mult(300);
80   mat_mult(500);
81   mat_mult(750);
82   mat_mult(1000);
83
84   return 0;
85 }