类型转换繁重程序中省略小数点

Question

我已经为高斯消除方法编写了这段代码：

#include <stdio.h>
int main()
{
    float arr[3][4]={{0}}, i, j, k, p;
    printf("Enter the coefficients:-\n");
    for(i=0; i<3;i++)
    {
        printf("Row #(%d)\n", (int)(i+1));
        for(j=0; j<4; j++)
        {
            scanf("%f", &arr[(int)i][(int)j]);
        }
    }
    for(k=1, p=(k-1); (k<=2); k++, p++)
    {
        for(i=k; i<=2; i++)
        {
            for(j=0; j<=3; j++)
            {
                arr[(int)i][(int)j]-=((arr[(int)i][(int)p]/arr[(int)p][(int)p])*arr[(int)p][(int)j]);
            }
        }
    }
    for(i=0; i<3;i++)
    {
        for(j=0; j<4; j++)
        {
            printf("%0.3f  ", arr[(int)i][(int)j]);
        }
        printf("\n");
    }
    return 0;
}

但是，它不会打印正确的输出值。对输出的分析表明这些值已被四舍五入或其他东西。我猜在某个地方已经考虑了二维数组的整数值，即使我将该数组作为双精度数组并且仅对语法有效的元素地址使用类型转换。你可以自己试试这个矩阵：

R1-   10  -1  2  4
R2-   1   10 -1  3
R3-   2   3  20  7

正确的输出应该是：

R1-   10    -1        2         4
R2-   0    10.1     -1.2       2.6
R3-   0     0       19.980    5.376

我正在使用 GNU GCC 编译器。

Answer 1

1. 用整数类型声明用于数组索引的所有变量。
2. 去掉所有的演员表。您不需要其中任何一个。
3. 用您只定义一次的Rows 和Columns 替换行号和列号的所有常量。
4. 在循环范围内保持一致；不要在某些地方使用j < 4，在其他地方不要使用j <= 3。 （两者都将变为j < Columns。）
5. 摆脱p。
6. k 上的循环迭代源行，即包含将用于取消其他元素的数据的行。它运行for (k = 0; k < Rows-1; ++k)。
7. i 中的循环迭代您为每个 k 取消的目标行。它运行for (i = k+1; i < Rows; ++i)。
8. 取消时要使用的因子是列k 的行i 和k 中的元素之间的比率。计算一次并调用它factor。
9. 对于i 行中的每一列，从那里的元素中减去factor 乘以行k 中同一列中的元素。

我如上所述编写了代码，它会打印出您所说的预期结果。

Answer 2

问题可能与舍入无关（您传递给数组的索引不应影响数组本身的任何舍入）。

我在主循环中添加了一个打印语句，只是为了查看在执行每行减法时使用了哪些系数，并且它们在每列之间似乎有所不同，因此您的数学很可能是错误的。我怀疑这与您的p 变量有关，该变量始终固定为k 上方的行，我认为这对您没有帮助。

您想要获取每对行的行系数，然后减去您用于消除的行乘以该系数。系数为 M(1,0)/M(0,0) 用于第二行逐第一行消除，M(2,1)/M(1,1) 用于第三行逐第二行消除，一般而言，M(I,K)/M(K,K)（其中 K 是正在使用的行要消除，而我是您要消除的那一行）。

最终你想做的事：

For each row K excluding the last row (0 to 1)...
    For each row I after it (K+1 to 2)...
        Get the row factor = M(I, K) / M(K, K).
        For each column J (0 to 3)...
            M(I, J) -= M(K, J) * factor

将其（几乎直接）翻译成代码，结果应该是：

10.000000  -1.000000  2.000000   4.000000
-0.000000  10.100000  -1.200000  2.600000
-0.000000  -0.000000  19.980198  5.376237