朴素高斯消除

时间:2014-10-12 20:13:39

标签: c++ double gaussian

在C ++中执行朴素高斯消元时,我可以得到一个上三角矩阵,但当方程组的系数大于5 * 5时,它不起作用。知道为什么会这样吗?我使用double作为系数的数据类型。

double factor =0;
    for(int k=0;k<(count-1);k++){
        for(int i=(k+1);i<count;i++){

            factor = (doubleCoefficient[i][k])/(doubleCoefficient[k][k]);

            for(int j=0;j<count;j++){
                doubleCoefficient[i][j] -= (doubleCoefficient[k][j]*factor);
            }
        }
    }

0 个答案:

没有答案