diff options
author | Gael Guennebaud <g.gael@free.fr> | 2008-07-09 14:04:48 +0000 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2008-07-09 14:04:48 +0000 |
commit | 28539e7597c643dbc2b8d4f49dd16bd86fb7251f (patch) | |
tree | 3bfb96ae898f0afc943a388edcc83c27d28d5b3a /bench/btl/libs/C | |
parent | 5f55ab524cf0aad34dd6884bcae09eaa6c43c247 (diff) |
imported a reworked version of BTL (Benchmark for Templated Libraries).
the modifications to initial code follow:
* changed build system from plain makefiles to cmake
* added eigen2 (4 versions: vec/novec and fixed/dynamic), GMM++, MTL4 interfaces
* added "transposed matrix * vector" product action
* updated blitz interface to use condensed products instead of hand coded loops
* removed some deprecated interfaces
* changed default storage order to column major for all libraries
* new generic bench timer strategy which is supposed to be more accurate
* various code clean-up
Diffstat (limited to 'bench/btl/libs/C')
-rw-r--r-- | bench/btl/libs/C/CMakeLists.txt | 2 | ||||
-rwxr-xr-x | bench/btl/libs/C/C_interface.hh | 96 | ||||
-rw-r--r-- | bench/btl/libs/C/main.cpp | 44 |
3 files changed, 142 insertions, 0 deletions
diff --git a/bench/btl/libs/C/CMakeLists.txt b/bench/btl/libs/C/CMakeLists.txt new file mode 100644 index 000000000..d3d2312d8 --- /dev/null +++ b/bench/btl/libs/C/CMakeLists.txt @@ -0,0 +1,2 @@ +include_directories(${PROJECT_SOURCE_DIR}/libs/f77) +btl_add_bench(btl_C main.cpp) diff --git a/bench/btl/libs/C/C_interface.hh b/bench/btl/libs/C/C_interface.hh new file mode 100755 index 000000000..b688f18fa --- /dev/null +++ b/bench/btl/libs/C/C_interface.hh @@ -0,0 +1,96 @@ +//===================================================== +// 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; + } + } + + 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 new file mode 100644 index 000000000..b15b7d9d5 --- /dev/null +++ b/bench/btl/libs/C/main.cpp @@ -0,0 +1,44 @@ +//===================================================== +// 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_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" + +int main() +{ + + bench<Action_matrix_vector_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; +} + + |