大家先来看两张对比的图，着一张是正确的输出结果：
再来看一张错误的输出结果：（这一个是我需要调的程序）
然后就开始要求我GF赶紧对公式，对比正确输出所用程序的代码，看看两者到底有什么地方不一样。由于两个程序所用的坐标系不一样，所以对比起来比较困难，最后费了好大的劲，把坐标系修改成完全一样的之后，发现问题依旧。这个时候另外一个现象更加困惑了我：在修改错误程序的坐标系之后，发现图形的缩小程度明显改善了，但是仍然是严重收敛的。这个现象着实打击了我一下，因为我一开始只是以为是误差量的积累造成的。难道说就是调整了一下几个公式的位置以及一些正负号就会造成误差累积量的不同？还真是前所未见的情况！
这里整个调试过程是非常痛苦的，因为总共大约有几千万次的循环，计算量之大不可想象。而输出正确和错误的程序之间，在大约每60次循环才会产生大约1x10-7细微差别。难为我GF在那里单步调试了几十万次的循环，真是可怜啊。到底是哪里出了问题呢？多亏我用VC7来运行了一遍，才发现问题所在。（之前都是我在抽空Review代码，GF在1000km+的距离以外用VC6在调试。）下面这段代码就是存在错误的地方：
void MultiMatr(int s, int n, int t, double *arr1, double *arr2, double *result)
{
double *tmp; for (int i = 0; i < s; i++)
{ for (int j = 0; j < t; j++) { tmp = result + t * i + j; *tmp = 0; for (int k = 0; k < n; k++) { *tmp += ((*(arr1 + n * i + k)) * (*(arr2 + t * k + j))); 
} } }
}/**//***************************** update ***************************/void UpdateQ(double WnbbA[3]){ int i; double deltaSita0, angSin, angCos; double deltaSita[3], qChange[3], qTrans[4][4]; deltaSita0 = 0; for (i = 0; i < 3; i++) { deltaSita[i] = WnbbA[i] * quaDel; deltaSita0 += deltaSita[i] * deltaSita[i]; } deltaSita0 = sqrt(deltaSita0); angCos = cos(deltaSita0 / 2); if (deltaSita0 == 0) angSin = 0.5; else angSin = sin(deltaSita0 / 2.0)/deltaSita0; qTrans[0][0] = qTrans[1][1] = qTrans[2][2] = qTrans[3][3] = angCos; qTrans[0][1] = qTrans[3][2] = -(angSin * deltaSita[0]); qTrans[1][0] = qTrans[2][3] = -qTrans[0][1]; qTrans[0][2] = qTrans[1][3] = -(angSin * deltaSita[1]); qTrans[2][0] = qTrans[3][1] = -qTrans[0][2]; qTrans[0][3] = qTrans[2][1] = -(angSin * deltaSita[2]); qTrans[3][0] = qTrans[1][2] = -qTrans[0][3]; MultiMatr(4, 4, 1, (double*)qTrans, (double*)q, (double*)qChange); for (i = 0; i < 4; i++) q[i] = qChange[i];}大家看出来那里错了吗？真是稀松平常的错误啊！害死人了。
打印 | 张贴于 2004-09-17 00:45:00 | Tag:其他
留言反馈
if(i=0){
.....
}
2。没有看到任何一行的校验，特别是对传入指针的较验
调试这种程序，称之为时间杀手，在公司里写出这样的代码被review到，简直不可想像。
写程序与性别无关，不能是MM就网开一面啊～呵呵
q[i] = qChange[i];
明显数组越界了。。。
声明的qChange[3]
这是纯C，表乱说~
汇编还没灭亡呢（主要是我不让），那理论得到c/c++？嗬嗬
就不要对女性的要求太高了，基本上这个程序片段你能够看得明白就不错了，不信的话下一次我把她师姐的程序拿出来给大家参考参考，你一定会郁闷死的。
基本上我查出来的错误就是数组/指针越界。MultiMatr本身就会因为越界修改了qTrans数组，然后对q的赋值也自然会有问题。至于和0比较，肯定没有问题，这一句是用来处理极限情况的。而q等变量，是全局变量，如果都贴出来太长了。这个错误要使光看代码，在茫茫C海当中还是比较难发现的，有点跟现在时下比较时髦的游戏“大家来找碴”的难度差不多。
其实我还是比较关心C++的，并不是说C++要灭亡了，标题当然要危言耸听才会有人看啦:P
不过很不幸的消息是，即使修正了这个错误，数据仍然不正确！哎，还要继续慢慢看代码……
不知道两位有没有试一下浮点和0比较会不会不准确，我是试过了，没有任何问题
浮点比较容易出问题的致因有两个，第一个是浮点的有效位，第二个是由于10进制小数不能精确的表示成2进制小数
这两个问题导致从形式的浮点数转化到存储的浮点数是可能有损失，因此，在下列情况下可以使用==操作符比较浮点数，1、立即数有效位数小于等于该浮点型的有效位数范围，2、该立即数是0或者2的整数次方，具体为什么
angSin = 0.5;
else
angSin = sin(deltaSita0 / 2.0)/deltaSita0;
总是执行else子句，if是执行不到的。
1、c/c++还是很好用的，相信这年头还有人用c++
2、想找出问题先看懂程序比较好
void MultiMatr(int s, int n, int t, double *arr1, double *arr2, double *result);
矩阵乘法，将s*n的矩阵arr1乘以n*t的矩阵arr2结果放在result中，隐含一个约束，就是resule是s*t的矩阵
void UpdateQ(double WnbbA[3]);
这里q是4维坐标系的点，4维方阵qTrans是变换矩阵，总体功能就是对q点进行变换
3、问题显而易见，1楼高人“叶君临”已经指出，访问越界了，由于没有具体给出q的定义，暂时仅能确定对qChange的访问越界
是不是这里的问题？
double qChange[3]-> double qChange[4]
我觉得你的GF基础不扎实，与C/C++无关。
中，q在MultiMatr已被引用，*(arr2 + t * k + j）=〉、q[0~3]
是不是提供的程序不全？
if (deltaSita0 == 0)
MultiMatr(4, 4, 1, (double*)qTrans, (double*)q, (double*)qChange);就出问题了qChange下标最多到2(range0,1,2),而
tmp = result + t * i + j; 可以指到result[3]。
double deltaSita[3], qChange[**3**], qTrans[4][4];
for (i = 0; i < **4**; i++)
q[i] = qChange[i];
so, q[3] = qTrans[0][0]
不知道是否因此引起错误？