← Back to branch summary

~ubuntu-branches/ubuntu/raring/simgrid/raring

« back to all changes in this revision

Viewing changes to examples/gras/chrono/chrono2.c

  • Committer: Package Import Robot
  • Author(s): Martin Quinson
  • Date: 2013-01-31 00:24:51 UTC
  • mfrom: (10.1.6 sid)
  • Revision ID: package-import@ubuntu.com-20130131002451-krejhf7w7h24lpsc
Tags: 3.9~rc1-1
* New upstream release: the "Grasgory" release. Major changes:
  - Gras was completely removed from this version.
  - Documentation reorganization to ease browsing it.
  - New default value for the TCP_gamma parameter: 4MiB

Show diffs side-by-side

added added

removed removed

Lines of Context:
examples/gras/chrono/chrono2.c
1
 
/* chrono - demo of GRAS benchmarking features                              */
2
 
 
3
 
/* Copyright (c) 2005, 2010. The SimGrid Team.
4
 
 * All rights reserved.                                                     */
5
 
 
6
 
/* This program is free software; you can redistribute it and/or modify it
7
 
 * under the terms of the license (GNU LGPL) which comes with this package. */
8
 
 
9
 
#include "gras.h"
10
 
#include "xbt/log.h"
11
 
 
12
 
XBT_LOG_NEW_DEFAULT_CATEGORY(Chrono, "Messages specific to this example");
13
 
 
14
 
#include <cblas.h>
15
 
void cblas_dgemm(const enum CBLAS_ORDER Order,
16
 
                 const enum CBLAS_TRANSPOSE TransA,
17
 
                 const enum CBLAS_TRANSPOSE TransB, const int M,
18
 
                 const int N, const int K, const double alpha,
19
 
                 const double *A, const int lda, const double *B,
20
 
                 const int ldb, const double beta, double *C,
21
 
                 const int ldc);
22
 
 
23
 
 
24
 
/* Function prototypes */
25
 
static int mat_mult(int n)
26
 
{
27
 
  int i, j, k, l;
28
 
  double *A, *B, *C;
29
 
  double start = 0.0;
30
 
  double now = 0.0;
31
 
 
32
 
  A = malloc(n * n * sizeof(double));
33
 
  B = malloc(n * n * sizeof(double));
34
 
  C = malloc(n * n * sizeof(double));
35
 
 
36
 
  start = now = gras_os_time();
37
 
 
38
 
  XBT_INFO("Matrix size: %d", n);
39
 
/*   XBT_INFO("Before computation: %lg", start); */
40
 
 
41
 
  for (l = 0; l < 40; l++) {
42
 
    now = gras_os_time();
43
 
    GRAS_BENCH_ONCE_RUN_ONCE_BEGIN();
44
 
    for (i = 0; i < n; i++)
45
 
      for (j = 0; j < n; j++) {
46
 
        A[i * n + j] = 2 / n;
47
 
        B[i * n + j] = 1 / n;
48
 
        C[i * n + j] = 0.0;
49
 
      }
50
 
 
51
 
    cblas_dgemm(CblasRowMajor,
52
 
                CblasNoTrans, CblasNoTrans, n, n, n,
53
 
                1.0, A, n, B, n, 0.0, C, n);
54
 
    GRAS_BENCH_ONCE_RUN_ONCE_END();
55
 
    now = gras_os_time() - now;
56
 
/*     XBT_INFO("Iteration %d : %lg ", l, now); */
57
 
  }
58
 
 
59
 
  now = gras_os_time() - start;
60
 
  XBT_INFO("Duration: %lg ", now);
61
 
 
62
 
  free(A);
63
 
  free(B);
64
 
  free(C);
65
 
}
66
 
 
67
 
int multiplier(int argc, char *argv[])
68
 
{
69
 
  gras_init(&argc, argv);
70
 
 
71
 
  mat_mult(10);
72
 
  mat_mult(20);
73
 
  mat_mult(30);
74
 
  mat_mult(40);
75
 
  mat_mult(50);
76
 
  mat_mult(75);
77
 
  mat_mult(100);
78
 
  mat_mult(300);
79
 
  mat_mult(500);
80
 
  mat_mult(750);
81
 
  mat_mult(1000);
82
 
 
83
 
  return 0;
84
 
}