全选主元高斯消去算法

/*
    全选主元高斯消去算法,求解线性方程组 A * x = b
    By Kingwei  2005.2.27
*/

#include <stdio.h>
#include <math.h>

#define EPS        1e-11

#define MAX_DIM        20

int GaussLEquation(double x[], int dim, double A[][MAX_DIM], double b[])
{
    int maxColPos[MAX_DIM], maxRowPos;
    double maxElem, absElem, temp;
    int k, i, j, status = 1;
    
    for (k=0; k<dim; k++)
    {
        maxElem = 0.0;
        for (i=k; i<dim; i++)
            for (j=k; j<dim; j++)
            {
                absElem = fabs(A[i][j]);
                if (maxElem < absElem)
                {
                    maxElem = absElem;
                    maxRowPos = i;
                    maxColPos[k] = j;