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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
#include <Eigen/Core>
#include <iostream>
#include <iomanip>
#include <string>
using Eigen::MatrixXd;
using Eigen::VectorXd;
void eigs_sym_F77(MatrixXd &M, VectorXd &init_resid, int k, int m, double &time_used, double &prec_err, int &nops);
void eigs_gen_F77(MatrixXd &M, VectorXd &init_resid, int k, int m, double &time_used, double &prec_err, int &nops);
void eigs_sym_Cpp(MatrixXd &M, VectorXd &init_resid, int k, int m, double &time_used, double &prec_err, int &nops);
void eigs_gen_Cpp(MatrixXd &M, VectorXd &init_resid, int k, int m, double &time_used, double &prec_err, int &nops);
void print_header(std::string title)
{
const int width = 80;
const char sep = ' ';
std::cout << std::endl
<< std::string(width, '=') << std::endl;
std::cout << std::string((width - title.length()) / 2, ' ') << title << std::endl;
std::cout << std::string(width, '-') << std::endl;
std::cout << std::left << std::setw(7) << std::setfill(sep) << "size";
std::cout << std::left << std::setw(10) << std::setfill(sep) << "dataset";
std::cout << std::left << std::setw(11) << std::setfill(sep) << "F77/time";
std::cout << std::left << std::setw(13) << std::setfill(sep) << "error";
std::cout << std::left << std::setw(7) << std::setfill(sep) << "nops";
std::cout << std::left << std::setw(11) << std::setfill(sep) << "C++/time";
std::cout << std::left << std::setw(13) << std::setfill(sep) << "error";
std::cout << std::left << std::setw(7) << std::setfill(sep) << "nops";
std::cout << std::endl;
std::cout << std::string(width, '-') << std::endl;
}
void print_row(int n, int dataset,
double time_f77, double err_f77, int nops_f77,
double time_cpp, double err_cpp, int nops_cpp)
{
const char sep = ' ';
std::cout.precision(5);
std::cout << std::left << std::setw(7) << std::setfill(sep) << n;
std::cout << std::left << std::setw(10) << std::setfill(sep) << dataset;
std::cout << std::left << std::setw(11) << std::setfill(sep) << time_f77;
std::cout << std::left << std::setw(13) << std::setfill(sep) << err_f77;
std::cout << std::left << std::setw(7) << std::setfill(sep) << nops_f77;
std::cout << std::left << std::setw(11) << std::setfill(sep) << time_cpp;
std::cout << std::left << std::setw(13) << std::setfill(sep) << err_cpp;
std::cout << std::left << std::setw(7) << std::setfill(sep) << nops_cpp;
std::cout << std::endl;
}
void print_footer()
{
const int width = 80;
std::cout << std::string(width, '=') << std::endl
<< std::endl;
}
void run_eigs_sym(int n_experiment, int n_replicate, int n, int k, int m)
{
double time_f77, time_cpp;
double err_f77, err_cpp;
int nops_f77, nops_cpp;
for (int i = 0; i < n_experiment; i++)
{
MatrixXd A = MatrixXd::Random(n, n);
MatrixXd M = A.transpose() + A;
VectorXd init_resid = VectorXd::Random(M.cols());
init_resid.array() -= 0.5;
init_resid = M * init_resid;
for (int j = 0; j < n_replicate; j++)
{
eigs_sym_F77(M, init_resid, k, m, time_f77, err_f77, nops_f77);
eigs_sym_Cpp(M, init_resid, k, m, time_cpp, err_cpp, nops_cpp);
print_row(n, i + 1, time_f77, err_f77, nops_f77, time_cpp, err_cpp, nops_cpp);
}
}
}
void run_eigs_gen(int n_experiment, int n_replicate, int n, int k, int m)
{
double time_f77, time_cpp;
double err_f77, err_cpp;
int nops_f77, nops_cpp;
for (int i = 0; i < n_experiment; i++)
{
MatrixXd A = MatrixXd::Random(n, n);
VectorXd init_resid = VectorXd::Random(A.cols());
init_resid.array() -= 0.5;
init_resid = A * init_resid;
for (int j = 0; j < n_replicate; j++)
{
eigs_gen_F77(A, init_resid, k, m, time_f77, err_f77, nops_f77);
eigs_gen_Cpp(A, init_resid, k, m, time_cpp, err_cpp, nops_cpp);
print_row(n, i + 1, time_f77, err_f77, nops_f77, time_cpp, err_cpp, nops_cpp);
}
}
}
int main()
{
std::srand(123);
int n_experiment = 10;
int n_replicate = 10;
print_header("eigs_sym");
run_eigs_sym(n_experiment, n_replicate, 100, 10, 20);
run_eigs_sym(n_experiment, n_replicate, 1000, 10, 30);
print_footer();
print_header("eigs_gen");
run_eigs_gen(n_experiment, n_replicate, 100, 10, 20);
run_eigs_gen(n_experiment, n_replicate, 1000, 10, 30);
print_footer();
return 0;
}
|
