pthread实现梯形积分法

下面是使用 pthread 实现梯形积分法的示例代码：

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>

#define MAX_THREADS 8

double a = 0, b = 1;    // 积分区间 [a, b]
int n = 100000000;      // 分割数
double h;               // 步长
double result = 0;      // 积分结果

void *integrate(void *tid) {
    int mytid = *((int*)tid);
    double local_result = 0;

    int left = mytid * (n / MAX_THREADS);
    int right = (mytid + 1) * (n / MAX_THREADS);
    if (mytid == MAX_THREADS - 1) {
        right = n;
    }

    for (int i = left; i < right; i++) {
        double x = a + i * h;
        local_result += (x * x + (x + h) * (x + h));
    }

    local_result *= h / 2.0;
    result += local_result;

    pthread_exit(NULL);
}

int main(int argc, char *argv[]) {
    h = (b - a) / n;

    pthread_t threads[MAX_THREADS];
    int tids[MAX_THREADS];

    for (int i = 0; i < MAX_THREADS; i++) {
        tids[i] = i;
        pthread_create(&threads[i], NULL, integrate, (void*)&tids[i]);
    }

    for (int i = 0; i < MAX_THREADS; i++) {
        pthread_join(threads[i], NULL);
    }

    printf("Result: %f\n", result);

    return 0;
}

该程序将积分区间 [a, b] 分成 n 个小区间，然后使用梯形积分法对每个小区间进行积分，最终将所有小区间的积分结果累加得到总积分结果。程序使用了 pthread 来实现多线程，在这里我们将积分区间划分为 MAX_THREADS 个部分，每个线程计算其中一个部分的积分，最后把所有线程的计算结果累加得到总积分结果。

需要注意的是，线程数量不能太多，否则可能会导致线程切换频繁而降低效率。在这里我们将线程数量设置为 8。

另外，需要注意的是，在每个线程中计算局部积分结果时，计算的小区间数量可能会有所不同。为了保证精确性，我们需要对最后一个线程计算的小区间数量进行特殊处理。在这里，我们将最后一个线程负责计算剩余的所有小区间。