20 #ifndef ACTION_AAT_PRODUCT
21 #define ACTION_AAT_PRODUCT
22 #include "utilities.h"
31 template<
class Interface>
40 MESSAGE(
"Action_aat_product Ctor");
44 init_matrix<pseudo_random>(A_stl,_size);
45 init_matrix<null_function>(X_stl,_size);
46 init_matrix<null_function>(resu_stl,_size);
50 Interface::matrix_from_stl(A_ref,A_stl);
51 Interface::matrix_from_stl(X_ref,X_stl);
53 Interface::matrix_from_stl(
A,A_stl);
54 Interface::matrix_from_stl(
X,X_stl);
62 INFOS(
"illegal call to Action_aat_product Copy Ctor");
70 MESSAGE(
"Action_aat_product Dtor");
74 Interface::free_matrix(
A,_size);
75 Interface::free_matrix(
X,_size);
77 Interface::free_matrix(A_ref,_size);
78 Interface::free_matrix(X_ref,_size);
84 static inline std::string
name(
void )
90 return double(_size)*double(_size)*double(_size);
102 Interface::aat_product(
A,
X,_size);
107 if (_size>128)
return;
110 Interface::matrix_to_stl(
X,resu_stl);
114 typename Interface::real_type
error=
118 INFOS(
"WRONG CALCULATION...residual=" <<
error);
126 typename Interface::stl_matrix
A_stl;
127 typename Interface::stl_matrix
X_stl;
130 typename Interface::gene_matrix
A_ref;
131 typename Interface::gene_matrix
X_ref;
133 typename Interface::gene_matrix
A;
134 typename Interface::gene_matrix
X;