1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//===================================================== 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// File : action_matrix_matrix_product.hh 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Author : L. Plagne <laurent.plagne@edf.fr)> 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) EDF R&D, lun sep 30 14:23:19 CEST 2002 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//===================================================== 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This program is free software; you can redistribute it and/or 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// modify it under the terms of the GNU General Public License 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// as published by the Free Software Foundation; either version 2 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// of the License, or (at your option) any later version. 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This program is distributed in the hope that it will be useful, 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// but WITHOUT ANY WARRANTY; without even the implied warranty of 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// GNU General Public License for more details. 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// You should have received a copy of the GNU General Public License 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// along with this program; if not, write to the Free Software 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifndef ACTION_TRMM 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define ACTION_TRMM 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "utilities.h" 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "STL_interface.hh" 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <string> 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "init/init_function.hh" 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "init/init_vector.hh" 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "init/init_matrix.hh" 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathusing namespace std; 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<class Interface> 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathclass Action_trmm { 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathpublic : 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Ctor 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Action_trmm( int size ):_size(size) 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MESSAGE("Action_trmm Ctor"); 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // STL matrix and vector initialization 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_matrix<pseudo_random>(A_stl,_size); 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_matrix<pseudo_random>(B_stl,_size); 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_matrix<null_function>(X_stl,_size); 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_matrix<null_function>(resu_stl,_size); 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0; j<_size; ++j) 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0; i<j; ++i) 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A_stl[j][i] = 0; 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A_stl[j][j] += 3; 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // generic matrix and vector initialization 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::matrix_from_stl(A_ref,A_stl); 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::matrix_from_stl(B_ref,B_stl); 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::matrix_from_stl(X_ref,X_stl); 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::matrix_from_stl(A,A_stl); 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::matrix_from_stl(B,B_stl); 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::matrix_from_stl(X,X_stl); 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath _cost = 0; 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0; j<_size; ++j) 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath _cost += 2*j + 1; 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath _cost *= _size; 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // invalidate copy ctor 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Action_trmm( const Action_trmm & ) 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INFOS("illegal call to Action_trmm Copy Ctor"); 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath exit(0); 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Dtor 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ~Action_trmm( void ){ 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MESSAGE("Action_trmm Dtor"); 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // deallocation 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_matrix(A,_size); 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_matrix(B,_size); 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_matrix(X,_size); 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_matrix(A_ref,_size); 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_matrix(B_ref,_size); 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::free_matrix(X_ref,_size); 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // action name 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline std::string name( void ) 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return "trmm_"+Interface::name(); 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double nb_op_base( void ){ 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return _cost; 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline void initialize( void ){ 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::copy_matrix(A_ref,A,_size); 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::copy_matrix(B_ref,B,_size); 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::copy_matrix(X_ref,X,_size); 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline void calculate( void ) { 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Interface::trmm(A,B,X,_size); 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath void check_result( void ){ 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // calculation check 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Interface::matrix_to_stl(X,resu_stl); 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// STL_interface<typename Interface::real_type>::matrix_matrix_product(A_stl,B_stl,X_stl,_size); 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// typename Interface::real_type error= 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// STL_interface<typename Interface::real_type>::norm_diff(X_stl,resu_stl); 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// if (error>1.e-6){ 135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// INFOS("WRONG CALCULATION...residual=" << error); 136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// // exit(1); 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// } 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathprivate : 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_matrix A_stl; 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_matrix B_stl; 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_matrix X_stl; 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::stl_matrix resu_stl; 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_matrix A_ref; 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_matrix B_ref; 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_matrix X_ref; 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_matrix A; 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_matrix B; 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typename Interface::gene_matrix X; 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int _size; 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double _cost; 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif 163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 166