高斯消元函数无法将结果赋值给新结构体变量的问题

为什么当前代码无法实现需求 #

你的代码存在三个核心问题：

1. 原地修改原矩阵：gauss函数直接对传入的MATRIX[n]的数组进行修改，没有生成新的矩阵副本，原矩阵会被直接覆盖，根本没有往MATRIX[n+1]写入任何有效数据。
2. 错误的返回值与赋值逻辑：函数最后返回的MATRIX[n+1].array[k][l]是数组越界访问——C语言数组下标从0开始，array[k][l]超出了float array[k][l]的有效范围（合法下标是0~k-1和0~l-1），而且这只是单个浮点数，不可能赋值整个矩阵。
3. 调用方式完全错误：MATRIX[n+1].array[k][l]=gauss(MATRIX,rem);试图把单个浮点数赋值给数组的越界位置，本质上和复制整个矩阵没有任何关系。

修正方案：复制原矩阵后在副本上消元 #

我们需要修改gauss函数，先将源矩阵完整复制到目标结构体，再在副本上执行消元操作，这样原矩阵保持不变，结果会存在新的结构体变量中。

修改后的gauss函数 #

#include <math.h> // 用于fabs判断浮点数是否为0

// src_idx: 源矩阵的下标，dest_idx: 存储结果的目标矩阵下标
void gauss(struct matrix MATRIX[], int src_idx, int dest_idx) {
    // 1. 完整复制源矩阵到目标矩阵（静态数组支持直接结构体赋值）
    MATRIX[dest_idx] = MATRIX[src_idx];
    
    int row = 0, flag = 0;
    while(row < k){
        int j=0, cid, rid;
        flag=0;
        // 寻找最左侧非零列
        while((j < l) && (flag != 1)){
            int i = row;    
            while((i < k) && (flag != 1)){    
                // 用epsilon判断浮点数是否为0，避免精度误差
                if(fabs(MATRIX[dest_idx].array[i][j]) > 1e-6){    
                    flag = 1;    
                    cid = j; 
                    rid = i;    
                }    
                i++;    
            } 
            j++;
        }
        
        // 交换行（仅当找到非零列时）
        if(flag != 0){
            for(int j=0; j < l; j++){
                float temp = MATRIX[dest_idx].array[rid][j];
                MATRIX[dest_idx].array[rid][j] = MATRIX[dest_idx].array[row][j];
                MATRIX[dest_idx].array[row][j] = temp; 
            }
            
            // 主元归一化为1
            float temp = MATRIX[dest_idx].array[row][cid];
            if(fabs(temp) > 1e-6){
                for(int j=0; j < l; j++){
                    MATRIX[dest_idx].array[row][j] /= temp;
                }
            }
            
            // 消去下方行的该列元素
            for(int i = row+1; i < k; i++){
                float temp = MATRIX[dest_idx].array[i][cid];
                for(int j=0; j < l; j++){
                    MATRIX[dest_idx].array[i][j] += (-1 * temp * MATRIX[dest_idx].array[row][j]);
                }
            }
        }
        row++;
    }
    
    // 修正负零为正零
    for(int i=0; i < k; i++){
        for(int j=0; j < l; j++){
            if(fabs(MATRIX[dest_idx].array[i][j]) < 1e-6){
                MATRIX[dest_idx].array[i][j] = 0.0f;
            }
        }
    }
    
    printf("Done final operations:\n"); 
    printmatrix(MATRIX, dest_idx);
}

调用方式 #

// 将MATRIX[rem]的行阶梯形结果保存到MATRIX[rem+1]
gauss(MATRIX, rem, rem+1);

额外注意事项 #

• 浮点数精度问题：原代码中直接用==或!=判断浮点数是否为0会有精度误差，改用fabs(val) < 1e-6的方式更可靠。
• 未初始化变量处理：如果矩阵全为0，原代码中cid会未初始化，修改后的代码通过flag判断跳过后续操作，避免未定义行为。

内容的提问来源于stack exchange，提问作者Yannis Ravasopoulos