diff options
Diffstat (limited to 'bench/btl/libs/C')
| -rw-r--r-- | bench/btl/libs/C/CMakeLists.txt | 3 | ||||
| -rwxr-xr-x | bench/btl/libs/C/C_interface.hh | 117 | ||||
| -rw-r--r-- | bench/btl/libs/C/main.cpp | 48 |
3 files changed, 0 insertions, 168 deletions
diff --git a/bench/btl/libs/C/CMakeLists.txt b/bench/btl/libs/C/CMakeLists.txt deleted file mode 100644 index 3d4d24cee..000000000 --- a/bench/btl/libs/C/CMakeLists.txt +++ /dev/null @@ -1,3 +0,0 @@ -include_directories(${PROJECT_SOURCE_DIR}/libs/f77) -btl_add_bench(btl_C main.cpp OFF) -# set_target_properties(btl_C PROPERTIES COMPILE_FLAGS "-fpeel-loops")
\ No newline at end of file diff --git a/bench/btl/libs/C/C_interface.hh b/bench/btl/libs/C/C_interface.hh deleted file mode 100755 index d6092517d..000000000 --- a/bench/btl/libs/C/C_interface.hh +++ /dev/null @@ -1,117 +0,0 @@ -//===================================================== -// File : C_interface.hh -// Author : L. Plagne <laurent.plagne@edf.fr)> -// Copyright (C) EDF R&D, lun sep 30 14:23:23 CEST 2002 -//===================================================== -// -// This program is free software; you can redistribute it and/or -// modify it under the terms of the GNU General Public License -// as published by the Free Software Foundation; either version 2 -// of the License, or (at your option) any later version. -// -// This program is distributed in the hope that it will be useful, -// but WITHOUT ANY WARRANTY; without even the implied warranty of -// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the -// GNU General Public License for more details. -// You should have received a copy of the GNU General Public License -// along with this program; if not, write to the Free Software -// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. -// -#ifndef C_INTERFACE_HH -#define C_INTERFACE_HH - -#include "f77_interface.hh" - -template<class real> -class C_interface : public f77_interface_base<real> { - -public : - - typedef typename f77_interface_base<real>::gene_matrix gene_matrix; - typedef typename f77_interface_base<real>::gene_vector gene_vector; - - static inline std::string name() { return "C"; } - - static inline void matrix_vector_product(const gene_matrix & A, const gene_vector & B, gene_vector & X, int N) - { -// for (int i=0;i<N;i++) -// { -// real somme = 0.0; -// for (int j=0;j<N;j++) -// somme += A[j*N+i] * B[j]; -// X[i] = somme; -// } - for (int i=0;i<N;i++) - X[i] = 0; - for (int i=0;i<N;i++) - { - real tmp = B[i]; - int iN = i*N; - for (int j=0;j<N;j++) - X[j] += tmp * A[j+iN]; - } - } - - static inline void atv_product(const gene_matrix & A, const gene_vector & B, gene_vector & X, int N) - { - for (int i=0;i<N;i++) - { - int iN = i*N; - real somme = 0.0; - for (int j=0;j<N;j++) - somme += A[iN+j] * B[j]; - X[i] = somme; - } - } - - static inline void matrix_matrix_product(const gene_matrix & A, const gene_matrix & B, gene_matrix & X, int N) - { - real somme; - for (int i=0;i<N;i++){ - for (int j=0;j<N;j++){ - somme=0.0; - for (int k=0;k<N;k++){ - somme += A[i+k*N] * B[k+j*N]; - } - X[i+j*N] = somme; - } - } - } - - static inline void ata_product(const gene_matrix & A, gene_matrix & X, int N) - { - - real somme; - for (int i=0;i<N;i++){ - for (int j=0;j<N;j++){ - somme=0.0; - for (int k=0;k<N;k++){ - somme+=A[k+i*N]*A[k+j*N]; - } - X[i+j*N]=somme; - } - } - } - - static inline void aat_product(const gene_matrix & A, gene_matrix & X, int N){ - real somme; - for (int i=0;i<N;i++){ - for (int j=0;j<N;j++){ - somme=0.0; - for (int k=0;k<N;k++){ - somme+=A[i+k*N]*A[j+k*N]; - } - X[i+j*N] = somme; - } - } - } - - static inline void axpy(real coef, const gene_vector & X, gene_vector & Y, int N){ - for (int i=0;i<N;i++) - Y[i]+=coef*X[i]; - } - - -}; - -#endif diff --git a/bench/btl/libs/C/main.cpp b/bench/btl/libs/C/main.cpp deleted file mode 100644 index f0a0e5c65..000000000 --- a/bench/btl/libs/C/main.cpp +++ /dev/null @@ -1,48 +0,0 @@ -//===================================================== -// File : main.cpp -// Author : L. Plagne <laurent.plagne@edf.fr)> -// Copyright (C) EDF R&D, lun sep 30 14:23:23 CEST 2002 -//===================================================== -// -// This program is free software; you can redistribute it and/or -// modify it under the terms of the GNU General Public License -// as published by the Free Software Foundation; either version 2 -// of the License, or (at your option) any later version. -// -// This program is distributed in the hope that it will be useful, -// but WITHOUT ANY WARRANTY; without even the implied warranty of -// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the -// GNU General Public License for more details. -// You should have received a copy of the GNU General Public License -// along with this program; if not, write to the Free Software -// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. -// -#include "utilities.h" -#include "bench.hh" -#include "C_interface.hh" -#include "action_matrix_vector_product.hh" -#include "action_atv_product.hh" -#include "action_matrix_matrix_product.hh" -#include "action_axpy.hh" -#include "action_ata_product.hh" -#include "action_aat_product.hh" -//#include "action_lu_solve.hh" -#include "timers/mixed_perf_analyzer.hh" - -BTL_MAIN; - -int main() -{ - - bench<Action_matrix_vector_product<C_interface<REAL_TYPE> > >(MIN_MV,MAX_MV,NB_POINT); - bench<Action_atv_product<C_interface<REAL_TYPE> > >(MIN_MV,MAX_MV,NB_POINT); - bench<Action_matrix_matrix_product<C_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT); - bench<Action_aat_product<C_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT); - bench<Action_ata_product<C_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT); - bench<Action_axpy<C_interface<REAL_TYPE> > >(MIN_AXPY,MAX_AXPY,NB_POINT); - - - return 0; -} - - |