Могу ли я рассчитать среднее число многопоточных номеров?

Question

Я использую многопоточную функцию в C ++. Используя n многопоточность, у меня n случайный вывод. Мне нужно рассчитать среднее значение многопотокового вывода в коде. допустим, что для n=4 потоков, код

#include <omp.h>
#include <unistd.h>
#include <stdio.h>
#include <random>
#include <iostream>
#include <cmath>
#include <iomanip>
#include <array>
#include <eigen3/Eigen/Dense>

#define         W               1.0
#define         avg_disorder    10
#define         numThd          4

int main()
{
  #pragma omp parallel num_threads(numThd)
  {
    // define random numbers
    std::mt19937 rng;                  
    std::uniform_real_distribution <> dist;
    std::random_device r;
    std::array<int,624> seed_data;
    std::generate(seed_data.begin(), seed_data.end(), std::ref(r));
    std::seed_seq seq(std::begin(seed_data), std::end(seed_data));
    rng.seed(seq);

    Eigen::Array<double, -1, 1> rp; // rp= random potential
    rp  = Eigen::Array<double, -1, 1>::Zero(avg_disorder, 1);

    //List of 10 random numbers
    for(int avr = 0; avr < avg_disorder; avr++ )
    {
      rp(avr,0) = 0.5* W* (-1 + 2* dist(rng) );
    }
    // output: the mean of the rp
    #pragma omp critical
    { 
      FILE *output;
      char name[50];
      sprintf(name, "avgW%06.2f.dat", W);    
      output = fopen(name, "a");
      fprintf(output, "%e \n", rp.mean() );
      fclose(output);     
    }

  }

  return 1;
}

// run the file avgThreads
//g++ -std=c++0x -fopenmp -o avgThreads avgThreads.cpp

Вывод: 4 потока: 4 результата (4 различных средних значения rp).

1.015983e-01
4.097469e-02
-1.275186e-01
-1.243190e-01

Мой вопрос:

Можно ли рассчитать среднее значение многопотокового вывода в коде?

Answer 1

Вы можете попробовать что-то вроде этого

#include <omp.h>
#include <stdio.h>
#include <random>
#include <iostream>
#include <cmath>
#include <iomanip>
#include <array>
#include <Eigen/Dense>

#define         W          1.0
#define    avg_disorder    10
#define     numThd          4

int main()
{
    double average = 0.0;
#pragma omp parallel num_threads(numThd) reduction(+ : average)  
{
    // define random numbers
    std::mt19937 rng;
    std::uniform_real_distribution <> dist;
    std::random_device r;
    std::array<int, 624> seed_data;
    std::generate(seed_data.begin(), seed_data.end(), std::ref(r));
    std::seed_seq seq(std::begin(seed_data), std::end(seed_data));
    rng.seed(seq);

    Eigen::Array<double, -1, 1> rp; // rp= random potential
    rp = Eigen::Array<double, -1, 1>::Zero(avg_disorder, 1);

    //List of 10 random numbers
    for (int avr = 0; avr < avg_disorder; avr++)
    {
        rp(avr, 0) = 0.5* W* (-1 + 2 * dist(rng));
    }
    average += rp.mean();
    // output: the mean of the rp
#pragma omp critical
 {
     cout << rp.mean() << endl;
     //FILE *output;
     //char name[50];
     //sprintf(name, "avgW%06.2f.dat", W);
     //output = fopen(name, "a");
     //fprintf(output, "%e \n", rp.mean());
     //fclose(output);
 }

}
cout << average / (double) numThd << endl;
return 1;
}

Используется уменьшение omp.

ДЛЯ ОБЩЕГО НАЗНАЧЕНИЯ (предположим, вы не знаете, сколько использованных потоков) вы можете сделать

double operation();
double averaged_operation()
{
    double average = 0.0, count_threads=0.0;
    #pragma omp parallel reduction(+ : average, count_threads)  
    {
        average += operation();
        count_threads += 1.0;
    }
    return  average / count_threads;
}