高斯约旦消元被卡精度怎么办

void guass(ll n)
{
    ll heng = 1;
    for(int i = 1 ; i <= n ; i++)
    {
        ll mx = heng;
        for(int j = heng + 1 ; j <= n ; j++)if(fabs(a[j][i]) - fabs(a[mx][i]) > eps)mx = j;
        if(fabs(a[mx][i]) < eps)continue;
        for(int j = 1 ; j <= n + 1 ; j++)swap(a[heng][j],a[mx][j]);
        for(int j = 1 ; j <= n ; j++)
        {
            if(heng == j)continue;
            db t = a[j][i] / a[heng][i];
            for(int k = 1 ; k <= n + 1 ; k++)a[j][k] -= a[heng][k] * t;
        }
        heng++;
    }
}