4.7 求一般行列式的值 sdet.c

许士良常用算法程序集C语言,包括c++一些常用算法代码

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