From 352e86db9858192c6c56709d58c2e3612daa5ec8 Mon Sep 17 00:00:00 2001
From: Neil Gershenfeld <gersh@cba.mit.edu>
Date: Fri, 21 Jun 2019 23:40:05 -0400
Subject: [PATCH] wip

---
 OpenMP/mppi.c | 6 ++----
 1 file changed, 2 insertions(+), 4 deletions(-)

diff --git a/OpenMP/mppi.c b/OpenMP/mppi.c
index 430f23c..da10b94 100644
--- a/OpenMP/mppi.c
+++ b/OpenMP/mppi.c
@@ -15,7 +15,6 @@ double pi;
 
 void main() {
    int i;
-   printf("%d\n",omp_get_max_threads());
    double a,b,c,dt,mflops;
    struct timespec tstart,tend;
    clock_gettime(CLOCK_REALTIME,&tstart);
@@ -24,13 +23,12 @@ void main() {
    c = 0.25;
    pi = 0;
    pi = 0;
-   #pragma omp parallel
-   #pragma omp for reduction(+:pi)
+   #pragma omp parallel for reduction(+:pi)
    for (i = 1; i <= NPTS; ++i)
       pi += a/((i-b)*(i-c));
    clock_gettime(CLOCK_REALTIME,&tend);
    dt = (tend.tv_sec+tend.tv_nsec/1e9)-(tstart.tv_sec+tstart.tv_nsec/1e9);
    mflops = NPTS*5.0/(dt*1e6);
-   printf("NPTS = %d, pi = %f\n",NPTS,pi);
+   printf("NPTS = %d, pi = %f, threads = %d\n",NPTS,pi,omp_get_max_threads());
    printf("time = %f, estimated MFlops = %f\n",dt,mflops);
    }
-- 
GitLab