nrfct.h

Go to the documentation of this file.

 
#ifndef nr_fct_h
#define nr_fct_h
 
int ludcmp(double **a, int n, int *indx, double *d);
int ludcmp_0(double **a, int n, int *indx, double *d);
 
void lubksb(double **a, int n, int *indx, double b[]);
void lubksb_0(double **a, int n, int *indx, double b[]);
 
int choldc(double **a, int n, double p[]);
int choldc_0(double **a, int n, double *p);
 
void cholsl(double **a, int n, double p[], double b[], double x[]);
void cholsl_0(double **a, int n, double *p, double *b, double *x);
 
void svdcmp(double **a, int m, int n, double w[], double **v);
 
void svdcmp_0(double **a, int m, int n, double w[], double **v);
 
void svbksb(double **u, double w[], double **v, int m, int n, double b[], double x[]);
void svbksb_0(double **u, double w[], double **v, int m, int n, double b[], double x[]);
 
void gaussj(double **a, int n, double **b, int m);
 
#define SIGN(a,b) ((b) >= 0.0 ? fabs(a) : -fabs(a))
//static int iminarg1,iminarg2;
#define IMIN(a,b) (iminarg1=(a),iminarg2=(b),(iminarg1) < (iminarg2) ?\
        (iminarg1) : (iminarg2))
//static double dmaxarg1,dmaxarg2;
#define DMAX(a,b) (dmaxarg1=(a),dmaxarg2=(b),(dmaxarg1) > (dmaxarg2) ?\
        (dmaxarg1) : (dmaxarg2))
//static double dsqrarg;
#define DSQR(a) ((dsqrarg=(a)) == 0.0 ? 0.0 : dsqrarg*dsqrarg)
 
 
#endif