1
      2
      3
      4
      5
      6
      7
      8
      9
     10
     11
     12
     13
     14
     15
     16
     17
     18
     19
     20
     21
     22
     23
     24
     25
     26
     27
     28
     29
     30
     31
#include <complex>
#include <cmath>
#include <cstdio>
#include "bench.h"

int
main()
{
    bench::WallTime<double>  wallTime;

//
//  Initialize a matrix of doubles
//
    wallTime.tic();

    double *A = new double[1000*1000];
    bench::initMatrix(1000, 1000, A, 1, 1000);
    delete [] A;

    std::printf("time elapsed: %lf\n", wallTime.toc());

//
//  Initialize a matrix of complex doubles
//
    wallTime.tic();
    std::complex<double> *B = new std::complex<double>[1000*1000];
    bench::initMatrix(1000, 1000, B, 1, 1000);
    delete [] B;
    std::printf("time elapsed: %lf\n", wallTime.toc());

}