summaryrefslogtreecommitdiff
path: root/libs/numeric/odeint/performance/mpi
diff options
context:
space:
mode:
Diffstat (limited to 'libs/numeric/odeint/performance/mpi')
-rw-r--r--libs/numeric/odeint/performance/mpi/Jamfile.v220
-rw-r--r--libs/numeric/odeint/performance/mpi/osc_chain_1d.cpp112
-rw-r--r--libs/numeric/odeint/performance/mpi/osc_chain_1d_system.hpp88
-rwxr-xr-xlibs/numeric/odeint/performance/mpi/osc_chain_speedup.gnu38
-rwxr-xr-xlibs/numeric/odeint/performance/mpi/osc_chain_speedup.sh24
5 files changed, 0 insertions, 282 deletions
diff --git a/libs/numeric/odeint/performance/mpi/Jamfile.v2 b/libs/numeric/odeint/performance/mpi/Jamfile.v2
deleted file mode 100644
index d80e860d8..000000000
--- a/libs/numeric/odeint/performance/mpi/Jamfile.v2
+++ /dev/null
@@ -1,20 +0,0 @@
-# Copyright 2011-2013 Mario Mulansky
-# Copyright 2012 Karsten Ahnert
-# Copyright 2013 Pascal Germroth
-# Distributed under the Boost Software License, Version 1.0. (See
-# accompanying file LICENSE_1_0.txt or copy at
-# http://www.boost.org/LICENSE_1_0.txt)
-
-import mpi ;
-
-project
- : requirements
- <include>../../../../..
- <include>..
- <define>BOOST_ALL_NO_LIB=1
- <library>/boost//timer
- <library>/boost//mpi
- <library>/boost//program_options
- ;
-
-exe osc_chain_1d : osc_chain_1d.cpp ;
diff --git a/libs/numeric/odeint/performance/mpi/osc_chain_1d.cpp b/libs/numeric/odeint/performance/mpi/osc_chain_1d.cpp
deleted file mode 100644
index ef2ce7dbd..000000000
--- a/libs/numeric/odeint/performance/mpi/osc_chain_1d.cpp
+++ /dev/null
@@ -1,112 +0,0 @@
-/* Boost libs/numeric/odeint/performance/openmp/osc_chain_1d.cpp
-
- Copyright 2013 Karsten Ahnert
- Copyright 2013 Mario Mulansky
- Copyright 2013 Pascal Germroth
-
- stronlgy nonlinear hamiltonian lattice in 2d
-
- Distributed under the Boost Software License, Version 1.0.
-(See accompanying file LICENSE_1_0.txt or
- copy at http://www.boost.org/LICENSE_1_0.txt)
- */
-
-#include <iostream>
-#include <vector>
-
-#include <boost/numeric/odeint.hpp>
-#include <boost/numeric/odeint/external/mpi/mpi.hpp>
-
-#include <boost/program_options.hpp>
-#include <boost/random.hpp>
-#include <boost/timer/timer.hpp>
-#include <boost/foreach.hpp>
-#include <boost/accumulators/accumulators.hpp>
-#include <boost/accumulators/statistics/stats.hpp>
-#include <boost/accumulators/statistics/mean.hpp>
-#include <boost/accumulators/statistics/median.hpp>
-#include "osc_chain_1d_system.hpp"
-
-using namespace std;
-using namespace boost::numeric::odeint;
-using namespace boost::accumulators;
-using namespace boost::program_options;
-
-using boost::timer::cpu_timer;
-
-const double p_kappa = 3.3;
-const double p_lambda = 4.7;
-
-int main( int argc , char* argv[] )
-{
- boost::mpi::environment env(argc, argv);
- boost::mpi::communicator world;
-
- size_t N, steps, repeat;
- bool dump;
- options_description desc("Options");
- desc.add_options()
- ("help,h", "show this help")
- ("length", value(&N)->default_value(1024), "length of chain")
- ("steps", value(&steps)->default_value(100), "simulation steps")
- ("repeat", value(&repeat)->default_value(25), "repeat runs")
- ("dump", bool_switch(&dump), "dump final state to stderr (on node 0)")
- ;
- variables_map vm;
- store(command_line_parser(argc, argv).options(desc).run(), vm);
- notify(vm);
- if(vm.count("help"))
- {
- if(world.rank() == 0)
- cerr << desc << endl;
- return EXIT_FAILURE;
- }
- cout << "length\tsteps\tthreads\ttime" << endl;
-
- accumulator_set< double, stats<tag::mean, tag::median> > acc_time;
-
- vector<double> p( N ), q( N, 0 );
- if(world.rank() == 0) {
- boost::random::uniform_real_distribution<double> distribution;
- boost::random::mt19937 engine( 0 );
- generate( p.begin() , p.end() , boost::bind( distribution , engine ) );
- }
-
- typedef vector<double> inner_state_type;
- typedef mpi_state< inner_state_type > state_type;
- typedef symplectic_rkn_sb3a_mclachlan<
- state_type , state_type , double
- > stepper_type;
- state_type p_split( world ), q_split( world );
- split(p, p_split);
- split(q, q_split);
-
- for(size_t n_run = 0 ; n_run != repeat ; n_run++) {
- cpu_timer timer;
- world.barrier();
- integrate_n_steps( stepper_type() , osc_chain( p_kappa , p_lambda ) ,
- make_pair( boost::ref(q_split) , boost::ref(p_split) ) ,
- 0.0 , 0.01 , steps );
- world.barrier();
- if(world.rank() == 0) {
- double run_time = static_cast<double>(timer.elapsed().wall) * 1.0e-9;
- acc_time(run_time);
- cout << N << '\t' << steps << '\t' << world.size() << '\t' << run_time << endl;
- }
- }
-
- if(dump) {
- unsplit(p_split, p);
- if(world.rank() == 0) {
- copy(p.begin(), p.end(), ostream_iterator<double>(cerr, "\t"));
- cerr << endl;
- }
- }
-
- if(world.rank() == 0)
- cout << "# mean=" << mean(acc_time)
- << " median=" << median(acc_time) << endl;
-
- return 0;
-}
-
diff --git a/libs/numeric/odeint/performance/mpi/osc_chain_1d_system.hpp b/libs/numeric/odeint/performance/mpi/osc_chain_1d_system.hpp
deleted file mode 100644
index 697208458..000000000
--- a/libs/numeric/odeint/performance/mpi/osc_chain_1d_system.hpp
+++ /dev/null
@@ -1,88 +0,0 @@
-/* Boost libs/numeric/odeint/performance/openmp/osc_chain_1d_system.hpp
-
- Copyright 2013 Karsten Ahnert
- Copyright 2013 Mario Mulansky
- Copyright 2013 Pascal Germroth
-
- stronlgy nonlinear hamiltonian lattice
-
- Distributed under the Boost Software License, Version 1.0.
-(See accompanying file LICENSE_1_0.txt or
- copy at http://www.boost.org/LICENSE_1_0.txt)
- */
-
-#ifndef SYSTEM_HPP
-#define SYSTEM_HPP
-
-#include <vector>
-#include <cmath>
-#include <iostream>
-
-#include <boost/math/special_functions/sign.hpp>
-#include <boost/numeric/odeint/external/mpi/mpi.hpp>
-
-namespace checked_math {
- inline double pow( double x , double y )
- {
- if( x==0.0 )
- // 0**y = 0, don't care for y = 0 or NaN
- return 0.0;
- using std::pow;
- using std::abs;
- return pow( abs(x) , y );
- }
-}
-
-double signed_pow( double x , double k )
-{
- using boost::math::sign;
- return checked_math::pow( x , k ) * sign(x);
-}
-
-struct osc_chain {
- const double m_kap, m_lam;
- osc_chain( const double kap , const double lam )
- : m_kap( kap ) , m_lam( lam ) { }
-
- void operator()( const boost::numeric::odeint::mpi_state< std::vector<double> > &q ,
- boost::numeric::odeint::mpi_state< std::vector<double> > &dpdt ) const
- {
- const bool have_left = q.world.rank() > 0;
- const bool have_right = q.world.rank() + 1 < q.world.size();
- double q_left = 0, q_right = 0;
- boost::mpi::request r_left, r_right;
- if(have_left)
- {
- q.world.isend(q.world.rank() - 1, 0, q().front());
- r_left = q.world.irecv(q.world.rank() - 1, 0, q_left);
- }
- if(have_right)
- {
- q.world.isend(q.world.rank() + 1, 0, q().back());
- r_right = q.world.irecv(q.world.rank() + 1, 0, q_right);
- }
-
- double coupling_lr = 0;
- if(have_left)
- {
- r_left.wait();
- coupling_lr = signed_pow( q_left - q()[0] , m_lam-1 );
- }
- const size_t N = q().size();
- for(size_t i = 0 ; i < N-1 ; ++i)
- {
- dpdt()[i] = -signed_pow( q()[i] , m_kap-1 ) + coupling_lr;
- coupling_lr = signed_pow( q()[i] - q()[i+1] , m_lam-1 );
- dpdt()[i] -= coupling_lr;
- }
- dpdt()[N-1] = -signed_pow( q()[N-1] , m_kap-1 ) + coupling_lr;
- if(have_right)
- {
- r_right.wait();
- dpdt()[N-1] -= signed_pow( q()[N-1] - q_right , m_lam-1 );
- }
- }
- //]
-};
-
-#endif
diff --git a/libs/numeric/odeint/performance/mpi/osc_chain_speedup.gnu b/libs/numeric/odeint/performance/mpi/osc_chain_speedup.gnu
deleted file mode 100755
index 7f3877257..000000000
--- a/libs/numeric/odeint/performance/mpi/osc_chain_speedup.gnu
+++ /dev/null
@@ -1,38 +0,0 @@
-#!/usr/bin/env gnuplot
-
-set terminal pngcairo size 1000,1000
-set output "osc_chain_speedup.png"
-
-set multiplot layout 2,2
-
-set key left
-
-set xrange [1:16]
-set x2range [1:16]
-set x2tics 8 format ""
-set grid x2tics
-set yrange [0:8]
-
-set title "short: speedup"
-plot \
- "osc_chain_speedup-short.dat" i 0 u "block":"mul" w lp t "MPI" , \
- (x < 4 ? x : 4) lc 0 lt 0 t "target"
-
-unset key
-
-set title "long: speedup"
-plot \
- "osc_chain_speedup-long.dat" i 0 u "block":"mul" w lp, \
- (x < 4 ? x : 4) lc 0 lt 0
-
-set yrange [0:*]
-
-set title "short: time[s]"
-plot \
- "osc_chain_speedup-short.dat" i 0 u "block":"med" w lp
-
-set title "long: time[s]"
-plot \
- "osc_chain_speedup-long.dat" i 0 u "block":"med" w lp
-
-unset multiplot
diff --git a/libs/numeric/odeint/performance/mpi/osc_chain_speedup.sh b/libs/numeric/odeint/performance/mpi/osc_chain_speedup.sh
deleted file mode 100755
index 7a0ff7891..000000000
--- a/libs/numeric/odeint/performance/mpi/osc_chain_speedup.sh
+++ /dev/null
@@ -1,24 +0,0 @@
-#!/bin/bash
-
-bench="bin/clang-linux-3.2/release/threading-multi/osc_chain_1d"
-
-repeat=5
-maxnodes=16
-
-function run {
- n=$1
- steps=$2
- for ((nodes=1 ; nodes < $maxnodes ; nodes++)) ; do
- # swap stderr & stdout
- mpirun -np $nodes $bench $n $steps $repeat 3>&1 1>&2 2>&3
- done
-}
-
-function run_all {
- printf "n\tsteps\tnodes\ttime\n"
- run 256 1024
- run 4096 1024
- run 4194304 1
-}
-
-run_all | tee osc_chain_speedup.dat