21 #include "utilities.h"
30 template<
class Interface>
41 init_vector<pseudo_random>(X_stl,_size);
42 init_vector<pseudo_random>(Y_stl,_size);
43 init_vector<null_function>(resu_stl,_size);
46 Interface::vector_from_stl(X_ref,X_stl);
47 Interface::vector_from_stl(Y_ref,Y_stl);
49 Interface::vector_from_stl(
X,X_stl);
50 Interface::vector_from_stl(
Y,Y_stl);
56 INFOS(
"illegal call to Action_axpby Copy Ctor");
65 Interface::free_vector(X_ref);
66 Interface::free_vector(Y_ref);
68 Interface::free_vector(
X);
69 Interface::free_vector(
Y);
73 static inline std::string
name(
void )
83 Interface::copy_vector(X_ref,
X,_size);
84 Interface::copy_vector(Y_ref,
Y,_size);
89 Interface::axpby(_alpha,
X,_beta,
Y,_size);
94 if (_size>128)
return;
96 Interface::vector_to_stl(
Y,resu_stl);
100 typename Interface::real_type
error=
104 INFOS(
"WRONG CALCULATION...residual=" <<
error);
111 typename Interface::stl_vector
X_stl;
112 typename Interface::stl_vector
Y_stl;
115 typename Interface::gene_vector
X_ref;
116 typename Interface::gene_vector
Y_ref;
118 typename Interface::gene_vector
X;
119 typename Interface::gene_vector
Y;
122 typename Interface::real_type
_beta;