1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This program is free software; you can redistribute it and/or 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// modify it under the terms of the GNU General Public License 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// as published by the Free Software Foundation; either version 2 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// of the License, or (at your option) any later version. 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This program is distributed in the hope that it will be useful, 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// but WITHOUT ANY WARRANTY; without even the implied warranty of 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// GNU General Public License for more details. 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// You should have received a copy of the GNU General Public License 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// along with this program; if not, write to the Free Software 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifndef ACTION_ROT 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define ACTION_ROT 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "utilities.h" 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "STL_interface.hh" 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <string> 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "init/init_function.hh" 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "init/init_vector.hh" 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "init/init_matrix.hh" 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathusing namespace std; 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<class Interface> 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathclass Action_rot { 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathpublic : 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Ctor 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_DONT_INLINE Action_rot( int size ):_size(size) 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MESSAGE("Action_rot Ctor"); 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // STL matrix and vector initialization 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_matrix tmp; 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_vector<pseudo_random>(A_stl,_size); 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_vector<pseudo_random>(B_stl,_size); 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // generic matrix and vector initialization 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::vector_from_stl(A_ref,A_stl); 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::vector_from_stl(A,A_stl); 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::vector_from_stl(B_ref,B_stl); 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::vector_from_stl(B,B_stl); 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // invalidate copy ctor 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Action_rot( const Action_rot & ) 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INFOS("illegal call to Action_rot Copy Ctor"); 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath exit(1); 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Dtor 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_DONT_INLINE ~Action_rot( void ){ 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MESSAGE("Action_rot Dtor"); 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_vector(A); 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_vector(B); 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_vector(A_ref); 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_vector(B_ref); 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // action name 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline std::string name( void ) 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return "rot_" + Interface::name(); 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double nb_op_base( void ){ 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return 6.0*_size; 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_DONT_INLINE void initialize( void ){ 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::copy_vector(A_ref,A,_size); 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::copy_vector(B_ref,B,_size); 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_DONT_INLINE void calculate( void ) { 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_ASM_COMMENT("#begin rot"); 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::rot(A,B,0.5,0.6,_size); 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_ASM_COMMENT("end rot"); 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath BTL_DONT_INLINE void check_result( void ){ 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // calculation check 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Interface::vector_to_stl(X,resu_stl); 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// STL_interface<typename Interface::real_type>::rot(A_stl,B_stl,X_stl,_size); 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// typename Interface::real_type error= 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// STL_interface<typename Interface::real_type>::norm_diff(X_stl,resu_stl); 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// if (error>1.e-3){ 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// INFOS("WRONG CALCULATION...residual=" << error); 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// exit(0); 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// } 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathprivate : 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_vector A_stl; 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_vector B_stl; 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_vector A_ref; 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_vector B_ref; 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_vector A; 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_vector B; 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int _size; 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif 117