【发布时间】:2023-03-16 12:15:01
【问题描述】:
我的程序中有一个非常简单的函数,它使用割线方法来查找函数的根。它适用于下面所示的 printf()。但是如果我把它注释掉,循环就会无休止地重复。我不知道为什么...
我已经阅读了有关 printf 更改变量的信息,但我没有看到任何内容更改了变量的存储。我错过了什么吗?打印它不是一个很好的解决方案,因为迭代并不重要,而且函数被调用了数百万次。
double guess1=500.;
double y1=estimater(r,volume,guess1,adm,tm,rr[r]);
double guess2=adm/30.;
double y2=estimater(r,volume,guess2,adm,tm,rr[r]);
int i;
double guess3=0.;
double y3;
double tol =heightTOL;
int secmax=SECANTMAX;
for(i=1;i<=secmax;i++){
guess3=guess2-y2*(guess2-guess1)/(y2-y1);
if(guess3>adm/2.){
guess3=adm/2.;
}
if(guess3<=0.){
guess3=0.;
}
y3=estimater(r,volume,guess3,adm,tm,rr[r]);
y1=y2;
y2=y3;
guess1=guess2;
guess2=guess3;
if(fabs(guess2-guess1)<tol){
break;
}
if(i==secmax){
printf("\nRan out of iterations in height finder\n");
}
printf("%d %f",i,guess3); //THIS ONE HERE!!!!!!!!
}
return guess3;
【问题讨论】:
-
您可能在某处调用了未定义的行为。你能构造一个我们可以编译和运行的minimal test-case吗?
-
我最好的猜测:你在数组边界之外写,
rr[r]是我看到的唯一候选对象,它会覆盖i,因此i == secmax永远不会变为真,除非调用printf会导致布局更改,将i移动到它还没有(尚未)踩到的地方。 -
你在调试器中运行过这个吗?或者如果你在 linux 上,通过 valgrind 来检查内存错误?