4.8 求矩阵的秩 rank.c

来自「许士良常用算法程序集C语言,包括c++一些常用算法代码」· C语言 代码 · 共 56 行

#include "math.h"
int rank(a,m,n)
int m,n;
double a[];
{ 
	int i,j,k,nn,is,js,l,ll,u,v;
    double q,d;
    nn=m;
    if (m>=n) nn=n;
    k=0;
    for (l=0; l<=nn-1; l++)
    { 
		q=0.0;
        for (i=l; i<=m-1; i++)
			for (j=l; j<=n-1; j++)
			{ 
				ll=i*n+j; d=fabs(a[ll]);
				if (d>q) 
				{ 
					q=d; is=i; js=j;
				}
			}
        if (q+1.0==1.0) return(k);
        k=k+1;
        if (is!=l)
        { 
			for (j=l; j<=n-1; j++)
            { 
				u=l*n+j; v=is*n+j;
                d=a[u]; a[u]=a[v]; a[v]=d;
            }
        }
        if (js!=l)
        { 
			for (i=l; i<=m-1; i++)
            { 
				u=i*n+js; v=i*n+l;
                d=a[u]; a[u]=a[v]; a[v]=d;
            }
        }
        ll=l*n+l;
        for (i=l+1; i<=n-1; i++)
        { 
			d=a[i*n+l]/a[ll];
            for (j=l+1; j<=n-1; j++)
            { 
				u=i*n+j;
                a[u]=a[u]-d*a[l*n+j];
            }
        }
    }
    return(k);
}