Question

我需要在c中使用pthread制作Leibniz算法，现在我有了这个代码，但是目前线程实现需要顺序实现的同一时间，我认为它不是并发运行的。有人可以看到错误。

谢谢！

#include<stdio.h>
#include<math.h>
#include<pthread.h>
#include<stdlib.h>
#define NUM_THREADS 2
#define ITERATIONS 100000000
double result = 0.0;
void *leibniz(void *threadid){
  int size = ITERATIONS/NUM_THREADS;
  int start = (long)threadid * size;
  int end = ((long)threadid+1) * size;
  int i;
  for(i = start; i<end; i++){
    int denom = 2*i+1;
    result += pow(-1.0, i) * (1.0/denom);
  }
}

int main(){
  pthread_t threads[NUM_THREADS];
  long t;
  int rc;

  // CREATE
  for(t=0;t<NUM_THREADS;t++){
    rc = pthread_create(&threads[t], NULL, leibniz, (void *)t);
    if(rc){
      printf("ERROR: return code %d\n", rc);
    }
  }
  // JOIN
  for(t=0;t<NUM_THREADS;t++){
    rc = pthread_join(threads[t], NULL);
    if(rc){
      printf("ERROR: return code %d\n", rc);
      exit(-1);
    }
  }
  printf("Pi %f\n", result*4);
  exit(0);

}

感谢Jean-FrançoisFabre，我做了这些改变，现在它有效了！

double result=0.0;

void *leibniz(void *threadid){
  double local = 0.0;
  int size = ITERATIONS/NUM_THREADS;
  int start = (long)threadid * size;
  int end = ((long)threadid+1) * size;
  int i;
  for(i = start; i<end; i++){
    local += (i%2==0 ? 1 : -1) * (1.0/(2*i+1));
  }
  result += local*4;
}

Answer 1

我会尝试回答。

即使您的应用程序是多线程的，也不能保证每个核心有1个FPU。我对此知之甚少，但我认为某些AMD处理器实际上共享核心之间的FPU。

由于您的循环基本上是添加pow，因此它的99％FPU计算，所以如果FPU在您的计算机上共享，它就解释了瓶颈。

您可以通过不调用-1来计算1或-1来减少FPU的使用，这可能是一个标量操作，然后可能会产生影响。如果i是奇数，则使用1，否则使用double result = 0.0; void *leibniz(void *threadid){ double local = 0.0; int size = ITERATIONS/NUM_THREADS; int start = (long)threadid * size; int end = ((long)threadid+1) * size; int i; for(i = start; i<end; i++){ int denom = 2*i+1; // using a ternary/scalar speeds up the "pow" computation, multithread or not local += (i%2 ? -1 : 1) * (1.0/denom); } // you may want to protect that addition with a pthread_mutex // start of critical section result += local; // end of critical section }，或者在每次迭代时否定外部1 / -1变量。

另外，为了避免竞争条件，将结果累加到局部结果中，最后添加它（最后通过互斥体保护添加会更好）

{{1}}

http://wccftech.com/amd-one-fpu-per-core-design-zen-processors/

Answer 2

我在Windows上运行Visual Studio，并且我没有安装pthread，因此我使用Windows线程创建了一个测试程序。我将计算分为一个计算所有正项的函数和一个计算所有负项的函数。双精度不是问题，因为正和是＆lt; 22，负和是> -19。

处理器是英特尔3770K 3.5ghz（每个核心都有自己的FPU）。我测试连续调用两个函数而不是为第二个函数使用单独的线程，两个线程的情况是单线程情况的两倍，单线程~0.360秒，双线程〜= 0.180秒。

#include <stdio.h>
#include <time.h>
#include <windows.h>

static HANDLE ht1;                      /* thread handle */

static DWORD WINAPI Thread0(LPVOID);    /* thread functions */
static DWORD WINAPI Thread1(LPVOID);

static clock_t ctTimeStart;             /* clock values */
static clock_t ctTimeStop;
static double  dTime;

static double pip;              /* sum of positive terms */
static double pim;              /* sum of negative terms */
static double pi;               /* pi */

int main()
{
    ctTimeStart = clock();
    Thread0(NULL);
    Thread1(NULL);
    ctTimeStop = clock();
    dTime = (double)(ctTimeStop - ctTimeStart) / (double)(CLOCKS_PER_SEC);  
    pip *= 4.;          /* pip <  22 after *= 4. */
    pim *= 4.;          /* pim > -19 after *= 4. */
    pi = pip + pim;
    printf("%.16lf %.16lf %.16lf %2.5lf secs\n", pi, pip, pim, dTime);

    ctTimeStart = clock();
    ht1 = CreateThread(NULL, 0, Thread1, 0, 0, 0);
    Thread0(NULL);
    WaitForSingleObject(ht1, INFINITE); // wait for thead 1
    ctTimeStop = clock();
    dTime = (double)(ctTimeStop - ctTimeStart) / (double)(CLOCKS_PER_SEC);  
    pip *= 4.;          /* pip <  22 after *= 4. */
    pim *= 4.;          /* pim > -19 after *= 4. */
    pi = pip + pim;
    printf("%.16lf %.16lf %.16lf %2.5lf secs\n", pi, pip, pim, dTime);

    CloseHandle(ht1);

    return 0;
}

DWORD WINAPI Thread0(LPVOID lpVoid)
{
double pp = 0.;                         /* local sum */
int j;
    for(j = 200000001; j >= 0; j -= 4)
        pp += 1. / (double)(j);
    pip = pp;                           /* store sum */
    return 0;
}

DWORD WINAPI Thread1(LPVOID lpVoid)
{
double pm = 0.;                         /* local sum */
int j;
    for(j = 200000003; j >= 0; j -= 4)
        pm -= 1. / (double)(j);
    pim = pm;                           /* store sum */
    return 0;
}

pthread并发运行线程c

2 个答案: