浮点数高斯消元的模板

//接口：Gauss();解存在x数组中。equ，var分别表示方程个数和变量的个数，（因为是变量 的个数不算常数的个数）要赋值。a就是行列式，其中常数是在方程右侧的符号。
int const maxn = 2000;
const double eps=1e-10;
double a[maxn][maxn],x[maxn];
int equ,var;
int Gauss()
{
    int i, j, k, max_r, col;
    double tmp;
    col = 0;

    for(k = 0; k<equ && col<var; k++, col++)
    {
        max_r = k;
        for(i = k+1; i < equ; i++)
            if(fabs(a[i][col])-fabs(a[max_r][col]) > eps)
            max_r = i;

        if(max_r != k)
            for(j = k; j < var+1; j++)
            swap(a[k][j], a[max_r][j]);

        if(fabs(a[k][col]) < eps)
        {
            k--;
            continue;
        }
        for(i = k+1; i < equ; i++)
        {
            if(fabs(a[i][col]) > eps)
            {
                double t = a[i][col]/a[k][col];  //这里和整型的不同。
                a[i][col] = 0.0;

                for(j = col; j < var+1; j++)
                a[i][j] -= a[k][j]*t;
            }
        }
    }
    for(i = var-1; i >= 0; i--)
    {
        if(fabs(a[i][i]) < eps) continue;
        tmp = a[i][var];
        for(j = i+1; j < var; j++)
        if(a[i][j] != 0)
        tmp -= a[i][j]*x[j];

        //if(tmp%a[i][i] != 0) return -2;
        x[i] = tmp/a[i][i];
    }
    return 0;
}