int thrd_create(thrd_t *thr, thrd_start_t func, void *arg);

int thrd_join(thrd_t thr, int *result);

int thrd_detach(thrd_t thr);

#include <stdbool.h>
#include <threads.h>
#define MAX_THREADS 8      // 1、2、4、8……所創(chuàng)建線程數(shù)量的最大值
#define MIN_BLOCK_SIZE 100   // 一個數(shù)組塊的最小值
typedef struct         // 函數(shù)parallel_sum()的參數(shù)
{
  float *start;    // 傳遞給parallel_sum()的數(shù)組塊的起始地址
  int len;      // 數(shù)組塊長度
  int block_size;   // 最小數(shù)組塊的大小
  double sum;         // 求和結果
} Sum_arg;
int parallel_sum(void *arg);  // 線程函數(shù)的原型
// ---------------------------------------------------------------
// 計算數(shù)組元素的和，并寫入*sumPtr
// sum()調(diào)用函數(shù)parallel_sum()進行并行處理
// 返回值：如果沒有發(fā)生錯誤，則返回true；否則，返回false
bool sum(float arr[], int len, double* sumPtr)
{
  int block_size = len / MAX_THREADS;
  if (block_size < MIN_BLOCK_SIZE) block_size = len;
  Sum_arg args = { arr, len, block_size, 0.0 };
  if (parallel_sum(&args))
  { *sumPtr = args.sum; return true; }
  else
    return false;
}
// ---------------------------------------------------------------
// 遞歸輔助函數(shù)，用以將工作分解到幾個線程中處理
int parallel_sum(void *arg)
{
  Sum_arg *argp = (Sum_arg*)arg;       // 指向參數(shù)的指針
  if (argp->len <= argp->block_size)         // 如果length <= block_size，
                            // 對所有元素求和
  {                                 
    for (int i = 0; i < argp->len; ++i)
      argp->sum += argp->start[i];
    return 1;
  }
  else                        // 如果length > block_size，
                            // 分解數(shù)組
  {                                      
    int mid = argp->len / 2;
    Sum_arg arg2 = { argp->start+mid, argp->len-mid,
             argp->block_size, 0};    // 指定后一半數(shù)組
    argp->len = mid;           // 前一半數(shù)組的長度
    thrd_t th;               // 在新線程中處理前一半數(shù)組
    int res = 0;
    if (thrd_create(&th, parallel_sum, arg) != thrd_success)
      return 0;              // 沒能成功創(chuàng)建新線程
    if (!parallel_sum(&arg2))  // 在當前線程下，以遞歸方式處理后一半數(shù)組
    {
      thrd_detach(th); return 0;     // 遞歸調(diào)用失敗
    }
    thrd_join(th, &res);
    if (!res)
      return 0;      // 同級線程報告執(zhí)行失敗
    argp->sum += arg2.sum;
    return 1; 
  }
}