/* COPYRIGHT (2011-2012) by: Kevin Marco Erler (author), http://www.kevinerler.de AIU-FSU Jena (co-owner), http://www.astro.uni-jena.de SBSZ Jena-Göschwitz (co-owner), http://www.sbsz-jena.de BSZ-Hermsdorf (co-owner), http://www.bszh.de Advanced Licensing (dual license: COPYRIGHT and following licenses): License (international): CC-BY v3.0-unported or later - link: http://creativecommons.org/licenses/by/3.0/deed.en License (Germany): CC-BY v3.0-DE or later - link: http://creativecommons.org/licenses/by/3.0/de/ ------------------ Compilation requirements: Packages (x86-64): GCC >v4.2, compat. libstdc++ and GOMP v3.0 Normal-Compile with g++-Compiler (Red Hat GCC 4.4.5-6 x86-64 tested) + OpenMP v3.0 ([lib]GOMP v3.0 x86-64 tested) g++ -std=c++0x -m64 -fopenmp -Wall -Wextra -pedantic -pedantic-errors -lgomp -lm -s <source.cpp> -o <dest> Release-Compile with g++-Compiler (Red Hat GCC 4.4.5-6 x86-64 tested) + OpenMP v3.0 ([lib]GOMP v3.0 x86-64 tested) g++ -std=c++0x -m64 -fopenmp -Wall -Wextra -pedantic -pedantic-errors -lgomp -lm -O3 -s <source.cpp> -o <dest> Debug-Compile with g++-Compiler (Red Hat GCC 4.4.5-6 x86-64 tested) + OpenMP v3.0 ([lib]GOMP v3.0 x86-64 tested) g++ -std=c++0x -m64 -fopenmp -Wall -Wextra -pedantic -pedantic-errors -lgomp -lm -g -ggdb3 <source.cpp> -o <dest> */ // Includes of C/C++-Librarys for INTs, REAL/FLOATs, STRINGS, Math-Calc and I/O #include <climits> #include <cstdint> #include <cinttypes> #include <cfloat> #include <cwchar> #include <string> //std:string #include <string.h> #include <cstring> #include <cstdlib> #include <cstdio> #include <iostream> #include <sstream> #include <cmath> // Conditional compilation (conditional include) of the OpenMP-Mainlib for OpenMP-Support #ifdef _OPENMP #include <omp.h> #endif using namespace std; #define free(x) free(x); *x=NULL #define PRId128 "s" #define PRIi128 "s" #define PRIu128 "s" const uint64_t UINT64_MIN = 0; const __int128_t INT128_MIN = (__int128_t)((-170141183460469231731.687303715884105728) * pow(10,18)); const __int128_t INT128_MAX = (__int128_t)(( 170141183460469231731.687303715884105727) * pow(10,18)); const __uint128_t UINT128_MAX = (__uint128_t)((340282366920938463463.374607431768211455) * pow(10,18)); const __uint128_t UINT128_MIN = 0/* * pow(10,18)*/; std::ostream &operator<<(std::ostream &out, __uint128_t x) { if(x >= 10) { out << x / 10; } return out << static_cast<unsigned>(x % 10); } std::ostream &operator<<(std::ostream &out, __int128_t x) { if(x < 0) { out << '-'; x = -x; } return out << static_cast<__uint128_t>(x); } string INT128ToSTR(__int128_t x) { std::stringstream sstr; sstr<<x; return sstr.str(); } #define INT128ToCSTR(x) (INT128ToSTR(x)).c_str() string UINT128ToSTR(__uint128_t x) { std::stringstream sstr; sstr<<x; return sstr.str(); } #define UINT128ToCSTR(x) (UINT128ToSTR(x)).c_str() const uint64_t NF = 1000000000ULL; // Anzahl Brüche (num fractions) int main(int argc, char *argv[]) { // Runtime manipulation of OpenMP-state variables //omp_set_num_threads(4); omp_set_dynamic(0); // data declarations and implementations double starttime = 0.00, sdelay = 0.00, pdelay = 0.00; long double PiSingle = 0.000000000000000000, PiParallel = 0.000000000000000000; std::cout.precision(20); std::cout << "PiCalc (Leibniz-Reihe) (64-Bit)\n" << "=======================================================================\n" << "SERIELLE AUSFÜHRUNG:"; //--------------------------Begin: CPU-serial execution of algorithm------------------------------- starttime = omp_get_wtime(); //CPU-serial algorithm: for(uint64_t i = 0ULL;i<NF;++i) { PiSingle+=(pow(-1.000000000000000000,i)/((2.000000000000000000*i)+1.000000000000000000)); } PiSingle*=4.000000000000000000; sdelay = omp_get_wtime()-starttime; std::cout << " done\n"; //serial //--------------------------End: CPU-serial execution of algorithm--------------------------------- //--------------------------Begin: CPU-parallel OpenMP-execution of algorithm---------------------- std::cout << "PARALLELE AUSFÜHRUNG mit "; //long double *Threads = new long double[omp_get_max_threads()]; // create parallel region: shared(Threads,...) if using "Threads"-variable #pragma omp parallel default(none) shared(std::cout, starttime, pdelay, PiParallel) { #pragma omp master { std::cout << omp_get_num_threads() << " Threads:"; starttime = omp_get_wtime(); } //OpenMP-CPU-parallel algorithm with the Reduction-clause: #pragma omp flush #pragma omp for schedule(static) reduction(+: PiParallel) for(uint64_t k = 0ULL;k<NF;++k) { PiParallel+=(pow(-1.000000000000000000,k)/((2.000000000000000000*k)+1.000000000000000000)); //Threads[omp_get_thread_num()]=PiParallel; //Partial results } #pragma omp master { PiParallel*=4.000000000000000000; pdelay = omp_get_wtime()-starttime; if(omp_get_num_threads() >= 10) { std::cout << " done\n"; //parallel } else { std::cout << " done\n"; //parallel } /* //Partial results of each CPU-Thread for(int j=0;j<omp_get_num_threads();++j) { if(j<9) { if((Threads[j]*4.000000000000000000)>=0) { std::cout << "Teilergebnis von Thread " << (j+1) << ": " << (Threads[j]*4.000000000000000000) << '\n'; } else { std::cout << "Teilergebnis von Thread " << (j+1) << ":" << (Threads[j]*4.000000000000000000) << '\n'; } } else { if((Threads[j]*4.000000000000000000)>=0) { std::cout << "Teilergebnis von Thread " << (j+1) << ": " << (Threads[j]*4.000000000000000000) << '\n'; } else { std::cout << "Teilergebnis von Thread " << (j+1) << ":" << (Threads[j]*4.000000000000000000) << '\n'; } } }*/ } } //--------------------------End: CPU-parallel OpenMP-execution of algorithm------------------------ //--------------------------Analysis of results---------------------------------------------------- std::cout << "\nAuswertung:\n" << "***********************************************************************\n" << "Anzahl Brüche: " << NF << '\n' << "PI - SERIELL: " << PiSingle << "\n" << "PI - PARALLEL: " << PiParallel << "\n" << "Referenzwert: 3.141592653589793238462643383279502884197169399375105\n" << "Dauer - SERIELL: " << sdelay << " sec\n" << "Dauer - PARALLEL: " << pdelay << " sec\n"; //delete []Threads; getchar(); return 0; }