OpenMP: Was passiert, wenn die Schleifenvariable keine Ganzzahl ist, wie in diesem Beispiel?

Question

#pargma omp parallel for erfordert, dass die Schleifenvariable eine ganze Zahl ist. Was passiert also, wenn die Schleifenvariable keine Ganzzahl ist, wie in diesem Beispiel?

#include <stdio.h> 
    #include <math.h> 

    #define START (M_PI/2) 
    #define END (M_PI*2) 

    double f(double x) 
    { return sin(x/2)+1; 
    } 

    int main(int argc, char *argv[]) 
    { double total = 0, x; 
    int partitions; 
    double slice; 

    printf("How many partitions? "); fflush(stdout); 
    scanf("%d", &partitions); 
    slice = (END-START)/partitions; 
    for (x = START + (slice/2); x < END; x = x + slice) 
    total = total + f(x); 
    total = total * slice; 

    printf("The integration is %1.20f\n", total); 
    } 

Wie konvertiert man dieses Programm in ein OpenMP? Danke

Answer 1

Zuerst sollten Sie ein wenig vorsichtiger mit dem Stil Ihres Codes sein, vor allem mit der Einrückung. Code wie

for (x = START + (slice/2); x < END; x = x + slice) 
total = total + f(x); 
total = total * slice; 

kann sehr irreführend sein, und Sie Ihr Leben viel einfacher, wenn es als

for (x = START + (slice/2); x < END; x = x + slice) { 
    total = total + f(x); 
} 
total = total * slice; 

statt zu schreiben.

In Bezug auf Ihr Problem, wie von Andrew erwähnt, können Sie es ganz durch Schleifen über die Partitionen mit einer Integer-Variablen vermeiden kann, das heißt,

#pragma omp parallel for private(x) reduction(+:total) 
for(int i = 0; i < partitions; i++) { 
    x = START + slice/2 + i*slice; 
    total = total + f(x); 
} 
total = total * slice;