19 #ifndef ACTION_PARTIAL_LU
20 #define ACTION_PARTIAL_LU
21 #include "utilities.h"
30 template<
class Interface>
39 MESSAGE(
"Action_partial_lu Ctor");
42 init_matrix<pseudo_random>(X_stl,_size);
43 init_matrix<null_function>(C_stl,_size);
46 for (
int i=0;
i<_size; ++
i)
47 X_stl[
i][
i] = X_stl[
i][
i] * 1e2 + 1;
50 Interface::matrix_from_stl(X_ref,X_stl);
51 Interface::matrix_from_stl(
X,X_stl);
52 Interface::matrix_from_stl(
C,C_stl);
61 INFOS(
"illegal call to Action_partial_lu Copy Ctor");
69 MESSAGE(
"Action_partial_lu Dtor");
72 Interface::free_matrix(X_ref,_size);
73 Interface::free_matrix(
X,_size);
74 Interface::free_matrix(
C,_size);
79 static inline std::string
name(
void )
93 Interface::partial_lu_decomp(
X,
C,_size);
114 typename Interface::stl_matrix
X_stl;
115 typename Interface::stl_matrix
C_stl;
117 typename Interface::gene_matrix
X_ref;
118 typename Interface::gene_matrix
X;
119 typename Interface::gene_matrix
C;