C 中浮点乘法的意外输出

Question

我正在用 C 语言执行乘法和加法运算，但得到了错误的结果。

我正在执行以下操作：

#include <stdio.h>
#include <strlib.h>
#include <math.h>

int main() 
{
    float** fpimag;
    float** fpreal;

    // Some code to assign fpimag and fpreal
    // fpimag and fpreal are two dimensional arrays that are
    // passed by reference to some functions to assign data to them
    allocatememory(&fpimag,&fpreal);
    storedata(&fpimag,&fpreal);


    autofunction(&fpimag,&fpreal);
} // main()

void autofunction(float*** fpimag,float***fpreal) {
    float expreal;
    float expimag;

    expreal = cos((*afph_correct)[row][column]);
    expimag = sin((*afph_correct)[row][column]);

    (*fpimag)[row][column] = (*fpimag)[row][column] * expreal + (*fpreal)[row][column] * expimag;

 printf("Operation looks like this\n");
 printf("%f*%f + %f*%f\n",(*fpimag)[row][column],expreal,(*fpreal)[row][column],expimag);
 printf("The value is %f\n",(*fpimag)[row][column]);
} // autofunction

这是我得到的输出：

Operation looks like this
-0.003095*-0.431162 + 0.000027*-0.902275
The value is 0.003865

但是正确的答案应该是：

0.0013101

谢谢！

Answer 1

尝试 long double，比 float 更精确：
http://en.wikipedia.org/wiki/Long_double

请注意，您应该更改 printf （%Lf 表示长双精度）

Answer 2

我很惊讶它没有崩溃。事实上，我很惊讶它竟然能编译。如果它在没有警告的情况下编译，我会大吃一惊。

• 没有为 fpreal 和 fpimag 分配内存。
• fpreal 中的值在初始化之前使用。
• 使用二维数组访问但不给出任何维度。
• 将指针 fpreal 和 fpimag 作为浮点数打印出来。
• 尝试打印出操作前的值，但 fpimag 中的值已被覆盖。

（前两个可能被未显示的代码覆盖 - 其他三个看起来是明确的问题）

Answer 3

看下面两行：

(*fpimag)[row][column] = (*fpimag)[row][column] * expreal +
                         (*fpreal)[row][column] * expimag;
printf("%f*%f + %f*%f\n",fpimag,expreal,fpreal,expimag);

你看不出fpimag 和(*fpimag)[行][列]不一样吗？ （与 fpreal 相同）您正在编写指向浮点指针值的指针，而不是指针指向的值。试试这个：

float imag = (*fpimag)[row][column];
float real = (*fpreal)[row][column];
(*fpimag)[row][column] = imag * expreal + real * expimag;
printf("%f*%f + %f*%f\n",imag,expreal,real,expimag);

Answer 4

嘿，我发现问题非常明显：

我覆盖了 fpreal 的值，然后尝试在另一个计算中使用它。

因此我做了：

(*fpreal)[row][column] = (*fpreal)[row][column]*expreal - (*fpimag)[row][column]*expimag;
(*fpimag)[row][column] = (*fpimag)[row][column]*expreal + (*fpreal)[row][column]*expimag;

我已经更新了代码，以便将 fpreal 和 fpimag 的值存储在一些临时变量中，并使用它们进行计算，如下所示：

oldreal = (*fpreal)[row][column];
oldimag = (*fpimag)[row][column];

(*fpreal)[row][column] = oldreal*expreal - oldimag*expimag;
(*fpimag)[row][column] = oldimag*expreal + oldreal*expimag;

谢谢大家的帮助！