OpenMP метод исключения(фон Неймана)

#include <iostream>
#include <random>
#include <omp.h>
#include <cmath>
#include <math.h>
using std::cout;
using std::endl;
 
 
 
double f(double x){
    return (x * x * x * sin(x) * sin(x)) / sqrt(x * x + 8); 
}
 
int main() {
    double sum = 0.0, a = 0.0, b = 3.14, eps = 0.007, ans, h=1.92117;
    int n = 10000000, n_in=0;
 
    
    do{
        #pragma omp parallel     
        {
            std::default_random_engine randomize;
            std::uniform_real_distribution<double> unirnd(0.0,1.0);
            randomize.seed(100);     
    
            double x,y;
            #pragma omp for private(x) reduction(+:sum) 
 
                for(int i = 0;i <= n;i++){
                    x=a + (b - a)*unirnd(randomize);    
                    y=h*unirnd(randomize);  
                        if (y<=f(x)){
                                    i_in++; 
                            #pragma omp atomic 
                            sum+=f(x)/(n_in/n);
                        }
            }
            ans=sum / n;     
        }
    }
    while(eps >= ans);  
    cout << ans;     
return 0;
}