码迷,mamicode.com
首页 > 其他好文 > 详细

高斯消元 模板

时间:2017-09-01 13:43:20      阅读:239      评论:0      收藏:0      [点我收藏+]

标签:double   直接   open   lap   bsp   oid   ++   none   pre   

  直接上模板:

技术分享
//n是方程的个数
void gauss(Matrix A,int n)
{
    int i,j,k,r;
    for(int i=0; i<n; i++)
    {
        r=i;
        for( j=i+1; j<n; j++)
            if(fabs(A[j][i])>fabs(A[r][i]))r=j;
        if(r!=i)
            for(j=0; j<=n; j++)swap(A[r][j],A[i][j]);
        for(k=i+1; k<n; k++)
        {
            double f=A[k][i]/A[i][i];
            for(j=i; j<=n; j++)
                A[k][j]-=f*A[i][j];
        }
    }
    for(i=n-1; i>=0; i--)
    {
        for(j=i+1; j<n; j++)
            A[i][n]-=A[j][n]*A[i][j];
        A[i][n]/=A[i][i];
    }
}
template

 

高斯消元 模板

标签:double   直接   open   lap   bsp   oid   ++   none   pre   

原文地址:http://www.cnblogs.com/FriskyPuppy/p/7462694.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!