C 带互斥同步的Pi计算器

在这里完成一项作业。让代码运行并计算饼图，但随机值除外，我收到以下错误：

/10 3 用3个线程中的10个项计算的pi为3.14183961892940200045 *在“./piesync”中出错：free（）：下一个大小无效（快速）：0x0000000001ca3010*

/piesync 100 5 *在“./piesync”中出错：双重释放或损坏（输出）：0x0000000000ee5040*

我知道它可能与数组或互斥有关，但我不知道是什么

代码：

//皮尼·瓦克宁
#包括
#包括
#包括
//全局变量
int N，T；
双gpie=3.0；
pthread_mutex_t mutex；
//饼函数
void*pie_runner（void*arg）
{
长j=（长）arg；
下长=（N/T）*（j-1）+1；
长上=（（N/T）*（j））；
双myPartialSum=0；
//printf（“下部=%lu上部=%lu\n”，下部，上部）；
对于（long i=lower；i您正在索引超出数组范围的内容。您在此处为
T
元素分配了一个数组

pthread_t *tids = (pthread_t *) calloc(T, sizeof(pthread_t));

但是你在这里索引不正确

for(int k = 1; k<=T; k++)
    {
        pthread_join(tids[k], NULL);
    }

你有一个除以零的分数

if(i % 2 == 0){
  myPartialSum -= 4.0/((2*i)*(2*i+1)*(2*i+2));

对于j=0

修理这个

for(long i = lower; i <= upper; i++)
{        
  if(i % 2 == 0){
    if ( ((2*i)*(2*i+1)*(2*i+2)) == 0)
      myPartialSum = 0.0;
    else 
      myPartialSum -= 4.0/((2*i)*(2*i+1)*(2*i+2)); 
  }
  else{
    myPartialSum += 4.0/((2*i)*(2*i+1)*(2*i+2));
  }
}

for（long i=lower；ic使用零索引循环和数组/指针。此外，当
i==0
in
for时，您还有一个除以零的除法（long i=lower；我用一种丑陋的方式修复了它，但现在我得到了24个线程中55个项的pi计算值是3.14159027400620292947oh抱歉，当j=0时是这样，但它不是。
j
作为最小值
1
传递，然后在
long lower=（N/t）*（j-1）中使用）+1、 
因此，在所有情况下，表达式中的
lower>=1
和
i>=1
都会导致div 0错误。啊，是的！一旦调用循环被纠正，
j
可以
0
，但是
lower
只有在
N/T==1
时才会
0
。
if(i % 2 == 0){
  myPartialSum -= 4.0/((2*i)*(2*i+1)*(2*i+2));

for(long i = lower; i <= upper; i++)
{        
  if(i % 2 == 0){
    if ( ((2*i)*(2*i+1)*(2*i+2)) == 0)
      myPartialSum = 0.0;
    else 
      myPartialSum -= 4.0/((2*i)*(2*i+1)*(2*i+2)); 
  }
  else{
    myPartialSum += 4.0/((2*i)*(2*i+1)*(2*i+2));
  }
}