Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019 #ifndef ACTION_CHOLESKY
00020 #define ACTION_CHOLESKY
00021 #include "utilities.h"
00022 #include "STL_interface.hh"
00023 #include <string>
00024 #include "init/init_function.hh"
00025 #include "init/init_vector.hh"
00026 #include "init/init_matrix.hh"
00027
00028 using namespace std;
00029
00030 template<class Interface>
00031 class Action_cholesky {
00032
00033 public :
00034
00035
00036
00037 Action_cholesky( int size ):_size(size)
00038 {
00039 MESSAGE("Action_cholesky Ctor");
00040
00041
00042 init_matrix_symm<pseudo_random>(X_stl,_size);
00043 init_matrix<null_function>(C_stl,_size);
00044
00045
00046 for (int i=0; i<_size; ++i)
00047 X_stl[i][i] = std::abs(X_stl[i][i]) * 1e2 + 100;
00048
00049
00050 Interface::matrix_from_stl(X_ref,X_stl);
00051 Interface::matrix_from_stl(X,X_stl);
00052 Interface::matrix_from_stl(C,C_stl);
00053
00054 _cost = 0;
00055 for (int j=0; j<_size; ++j)
00056 {
00057 double r = std::max(_size - j -1,0);
00058 _cost += 2*(r*j+r+j);
00059 }
00060 }
00061
00062
00063
00064 Action_cholesky( const Action_cholesky & )
00065 {
00066 INFOS("illegal call to Action_cholesky Copy Ctor");
00067 exit(1);
00068 }
00069
00070
00071
00072 ~Action_cholesky( void ){
00073
00074 MESSAGE("Action_cholesky Dtor");
00075
00076
00077 Interface::free_matrix(X_ref,_size);
00078 Interface::free_matrix(X,_size);
00079 Interface::free_matrix(C,_size);
00080 }
00081
00082
00083
00084 static inline std::string name( void )
00085 {
00086 return "cholesky_"+Interface::name();
00087 }
00088
00089 double nb_op_base( void ){
00090 return _cost;
00091 }
00092
00093 inline void initialize( void ){
00094 Interface::copy_matrix(X_ref,X,_size);
00095 }
00096
00097 inline void calculate( void ) {
00098 Interface::cholesky(X,C,_size);
00099 }
00100
00101 void check_result( void ){
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113 }
00114
00115 private :
00116
00117 typename Interface::stl_matrix X_stl;
00118 typename Interface::stl_matrix C_stl;
00119
00120 typename Interface::gene_matrix X_ref;
00121 typename Interface::gene_matrix X;
00122 typename Interface::gene_matrix C;
00123
00124 int _size;
00125 double _cost;
00126 };
00127
00128 #endif