linearalgebra.h

Go to the documentation of this file.

#define NRANSI
#define THRESH 0
#define ITOL 3
#define TOL 0.00001
#define ITMAX 1000
#define SQR(a) ((sqrarg=(a)) == 0.0 ? 0.0 : sqrarg*sqrarg)
#define SWAP(a,b) tempr=(a);(a)=(b);(b)=tempr
#define EPS 1.0e-14
#define TINY 1.0E-20
static float sqrarg; /* this variables is used in file linearalgebra.h, it must not be removed, comment by Emanuele Cordano on 1 September 2009 */

void convlv(FLOATVECTOR *data,unsigned long n,FLOATVECTOR *respns,unsigned long m,int isign,
            double delta_time,FLOATVECTOR *ans);


void four1(FLOATVECTOR *data,unsigned long nn,int isign);


void ludcmp(SHORTVECTOR *indx, DOUBLEMATRIX *var);

void lubksb(DOUBLEMATRIX *var, SHORTVECTOR *indx,DOUBLEVECTOR *gam);

void realft(FLOATVECTOR *data,unsigned long n,int isign);

void twofft(FLOATVECTOR *data1,FLOATVECTOR *data2,FLOATVECTOR *fft1,
                 FLOATVECTOR *fft2,unsigned long n);

/* commentata per l'errore (dovuto ad una errata definizione di asolve):
error: conflicting types for `asolve'
void ris_sistema (double d[], double ds[], double di[], double b[], double x[], int n);

*/
void sprsin(double **,int,float,long,double *,long *);

/* commentata per l'errore (dovuto ad una errata definizione di asolve):
error: conflicting types for `asolve'
void linbcg(long, double *,double *, int, double, int, int *,double *,
             double *, long *);

*/
DOUBLEMATRIX *vett_mat (double *d,double *ds,double *di,int n);


double snrm(long n, double sx[], int itol);


void atimes(long n, double x[], double r[], int itrnsp,double sa[], long ija[]);

void asolve(long n, double b[], double x[],double sa[]);

void dsprsax(double sa[], long ija[], double x[], double b[], long n);


void dsprstx(double sa[], long ija[], double x[], double b[], long n);

float integration ( float ( *fun )( float ) , float , float, int);