-1

我的程序中有一个非常简单的函数，使用割线方法来查找函数的根。它适用于下面所示的 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;
4

1 回答 1

1

printf根本不会更改您的数据。printf函数族对您的数据产生任何影响的唯一方法是，如果您使用%n格式说明符（它写入您传递的地址的变量），或者您正在做一些调用未定义行为的事情，例如传递错误的格式参数. 您在这里没有做任何这些事情（您的格式字符串是正确的），所以您的错误在其他地方。

检查您的所有数组访问是否在界限内。尝试在 Valgrind 或其他验证器中运行您的代码以尝试查找内存错误。

于 2012-07-09T19:45:23.357 回答