2
我正在处理一个分配,我必须将C代码片段转换为内联汇编。该代码是呈现Julia分形的程序的一部分。C++和内联汇编代码在一个程序中具有相同的输出,但在另一个程序中具有相同的输出
我已经测试了两个代码片段的输出,并且它们完全匹配,但是我的程序仍然输出不同的图像(适用于C代码的Julia分形,内联汇编代码的扁平粉色屏幕)。
这是函数的开始部分以及返回
COLORREF render_point(const double &a,
const double &b, int N) {
double cRe = -0.5;
double cIm = -0.05;
double x(a), y(b);
double norm = x*x+y*y;
int n;
double three = 3.0;
(loop goes here)
return HSVtoRGB(n % 256, 255 , 255 *(n<N));
}
这里的C代码
for (n = 0; norm < 4.0 && n < N; ++n)
{
double old_x = x;
double old_y = y;
x = (old_x * old_x * old_x) - (3 * old_y * old_y * old_x) + cRe;
y = (3 * old_y * old_x * old_x) - (old_y * old_y * old_y) + cIm;
norm = x*x+y*y;
}
和内联汇编代码:
for (n = 0; norm < 4.0 && n < N; ++n)
{
__asm {
// Create (old_x * old_x * old_x)
fld x;
fmul x;
fmul x;
// Create (3 * old_y * old_y * old_x)
fld three;
fmul y;
fmul y;
fmul x;
// Create the full equation for x
fsubp st(1), st(0);
fadd cRe;
// Create (3 * old_y * old_x * old_x) + cIm
fld three;
fmul y;
fmul x;
fmul x;
fadd cIm;
// Create (old_y * old_y * old_y)
fld y;
fmul y;
fmul y;
fsubp st(1), st(0); // Create the full equation for y
fst y; // Store in y to use for next loop
fmul st(0), st(0); // Get y*y
fxch st(1); // Swap places of y*y with newly calculated x
fst x; // Store in x to use for next loop
fmul st(0), st(0); // Get x*x
faddp st(1), st(0); // Get x*x + y*y
fst norm; // Set loop variable
}
}
有两个循环之间的差异可能会导致程序中的输出不同?
你是否已经开始使用调试器并观察? –
在汇编循环结束时,FPU堆栈上仍有一个值。最后的指令应该是'fstp norm'。 – 1201ProgramAlarm
@ 1201ProgramAlarm是的,修复它!我不知道我不能把东西留在FPU堆栈上。我想,案件已经结案。 – ozma