gaussj.h

eC++编译器源码

#pragma GaussJ
/*NOT TESTED YET*/
  const unsigned int
    MAX = 100,        /* Maximum number of data points in array */
    MAxNoTerms = 15;  /* Max number of terms in fitted polynomial expression */

    typedef struct {
                  double coeff[MAxNoTerms+1];
                    /* coefficients for fitted polynomial expression */
                  int noTerms;  /* actual number of terms used in array */
                  int listA[MAxNoTerms+1];
                    /* array to keep order and track of terms to be fitted */
                  int minFit;   /* minimum number of terms to fit */
                    /* also number of terms of listA actually used */
                  double covarMatrix[MAX+1][MAxNoTerms+1];
                    /* max size of covariance matrix from function fitting */
                  int subI, subJ;
                    /* ordinates for actual size of covariance matrix */
                  double chi2;        /* chi-square */
     } PolyReturn;
     typedef double PolyArray[MAxNoTerms+1];

  void GaussJordan(PolyReturn &returnRec,
                        PolyArray &betaFunc,
                        int m, int mp);