Task 1 (Stochastic PI)

A possible solution can be found at:

Please tell me about any errors or problems you may encounter.

Scaffold Head: #include <iostream> #include <iomanip> #include <cstdint> #include <random> #include <mpi.h> constexpr double pi = 3.1415926535897932384626433832795028841971693993751058; double residue (const double& result) { return (result > pi ? result-pi : pi-result); } template <typename engine_t> double uniform_random(engine_t& engine) { std::uniform_real_distribution<double> dist(0.0, 1.0); return dist(engine); } int main (int argc, char * argv[]) { MPI::Init(argc , argv ); uint64_t num_samples=0; const uint64_t num_ranks = MPI::COMM_WORLD.Get_size(); const uint64_t rank = MPI::COMM_WORLD.Get_rank(); // actual value of num_samples only known to the root if (rank == 0) num_samples= 1UL<<24;
Scaffold Foot: if (rank == 0) { const double error = residue(result); std::cout << std::setprecision(16) << "# pi: " << result << " residue: " << error << "\nParallel programming is " << (error < 1E-3 ? "fun!" : "error-prone!") << std::endl; } MPI::Finalize(); }

Tests