tools.h File Reference

#include <math.h>
#include <iosfwd>
#include "Vector.h"
#include "Matrix.h"

Include dependency graph for tools.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Macros
#define	maxDWORD   4294967295
#define	INF   1.7E+308
#define	EOL   10
#define	PI   3.1415926535897932384626433832795
#define	ROOT2   1.41421356

Functions
double	condorAbs (const double t1)
double	sign (const double a)
double	sign (const double t1, const double t2)
double	isInt (const double a)
double	mmin (const double t1, const double t2)
unsigned	mmin (const unsigned t1, const unsigned t2)
int	mmin (const int t1, const int t2)
double	mmax (const double t1, const double t2)
unsigned	mmax (const unsigned t1, const unsigned t2)
int	mmax (int t1, int t2)
double	sqr (const double &t)
double	CondorRound (double a)
unsigned long	choose (unsigned n, unsigned k)
double	rand1 ()
void	initRandom (int i=0)
double	euclidianNorm (int i, double *xp)
char	emptyline (const char *line)
char *	GetRidOfTheEOL (char *tline)
char *	removeQuotes (char *t)
char **	getNameTable (const char *line, int *ncolumn)
char *	stringDuplicate (const char *l)
char	isEmpty (const char *line)
char *	removeAllEOL (char *t)
char *	loadFile (FILE *f)
const char *	skipLine (const char *l)
const char *	skipSpaces (const char *l)
void	deleteFile (const char *l)

Macro Definition Documentation

EOL

#define EOL   10

Definition at line 41 of file tools.h.

INF

#define INF   1.7E+308

Definition at line 40 of file tools.h.

maxDWORD

#define maxDWORD   4294967295

Definition at line 39 of file tools.h.

PI

#define PI   3.1415926535897932384626433832795

Definition at line 43 of file tools.h.

ROOT2

#define ROOT2   1.41421356

Definition at line 45 of file tools.h.

Function Documentation

choose()

unsigned long choose	(	unsigned	n,
		unsigned	k
	)

Definition at line 38 of file tools.cpp.

{
    const unsigned long uupSize = 100;
    static unsigned long uup[uupSize];
    unsigned long *up;
    static unsigned long Nold = 0;
    static unsigned long Kold = 0;
    unsigned long l,m;
    unsigned i,j;

    if ( (n < k) || !n ) return 0;

    if ( (n == k) || !k )   // includes n == 1
        return 1;

    if ( k > (n >> 1) ) // Only lower half
        k = n-k;

    if ( (Nold == n) && (k < Kold) )    // We did it last time ...
        return *(uup + k - 1);

    if ( k > uupSize )
    {
        printf( "choose( unsigned, unsigned) :  overflow\n");
        getchar(); exit(-1);
    }

    Nold=n; Kold=k;

    *(up=uup)=2;
    for (i=2; i<n; i++)     // Pascal's triangle
    {
        // todo: remove next line:
        *(up+1)=1;
        l=1;
        m=*(up=uup);
        for (j=0; j<mmin(i,k); j++)
        {
            *up=m+l;
            l=m;
            m=*(++up);
        }
        // todo: remove next line:
        *up=1;
    }

    return *(uup + k - 1);
}

condorAbs()

double condorAbs ( const double  t1 )

inline

Definition at line 47 of file tools.h.

{
    return t1 > 0.0 ? t1 : -t1;
}

CondorRound()

double CondorRound ( double  a )

inline

Definition at line 104 of file tools.h.

{
    return (int)(a+.5);
}

deleteFile()

void deleteFile

(

const char *

l

)

Definition at line 280 of file tools.cpp.

{
#ifdef WIN32
    remove(line);
#else
    unlink(line);
#endif
}

emptyline()

char emptyline

(

const char *

line

)

euclidianNorm()

double euclidianNorm	(	int	i,
		double *	xp
	)

Definition at line 113 of file tools.cpp.

{
// no tested
// same code for the Vector eucilidian norm and for the Matrix Froebenis norm
/*
    double sum=0;
    while (i--) sum+=sqr(*(xp++));
    return sqrt(sum);
*/
    const double SMALL=5.422e-20, BIG=1.304e19/((double)i);
    double s1=0,s2=0,s3=0, x1max=0, x3max=0, xabs;

    while (i--)
    {
        xabs=condorAbs(*(xp++));
        
        if (xabs>BIG)
        {
            if (xabs>x1max)
            {
                s1=1.0+s1*sqr(x1max/xabs);
                x1max=xabs;
                continue;
            }
            s1+=sqr(xabs/x1max);
            continue;
        }
        if (xabs<SMALL)
        {
            if (xabs>x3max)
            {
               s3=1.0+s3*sqr(x3max/xabs);
               x3max=xabs;
               continue;
            }
            if (xabs!=0) s3+=sqr(xabs/x3max);
            continue;
        }
        s2+=sqr(xabs);
    };
    if (s1!=0) return x1max*sqrt(s1+(s2/x1max)/x1max);
    if (s2!=0)
    {
        if (s2>=x3max) return sqrt(s2*(1.0+(x3max/s2)*(x3max*s3)));
        return sqrt(x3max*((s2/x3max)+(x3max*s3)));
    }
    return x3max*sqrt(s3);
}

getNameTable()

char** getNameTable	(	const char *	line,
		int *	ncolumn
	)

Definition at line 204 of file tools.cpp.

{
    char *n=(char*)line;
    int j=0,nc;
    if (ncolumn) nc=*ncolumn;
    if (nc==0)
    {
        while (*n) 
        {
            while ((*n)&&(*n!='\t')&&(*n!=' ')&&(*n!=13)&&(*n!=10)) n++;
            nc++;
            while ((*n)&&((*n=='\t')||(*n==' ')||(*n==13)||(*n==10))) n++;
        }
    }
    if (ncolumn) *ncolumn=nc;
    char **names=(char**)malloc(nc*sizeof(char**));
    n=(char*)malloc(strlen(line)+1);
    strcpy(n,line);
    
    for (j=0; j<nc-1; j++)
    {
        names[j]=n;
        while ((*n)&&(*n!='\t')&&(*n!=' ')&&(*n!=13)&&(*n!=10)) n++;
        if (*n) { *n=0; n++; }
        while ((*n)&&((*n=='\t')||(*n==' ')||(*n==13)||(*n==10))) n++;
    }
    names[j]=n;
    return names;
}

GetRidOfTheEOL()

char* GetRidOfTheEOL

(

char *

tline

)

Definition at line 185 of file tools.cpp.

{
   char *t;
   t=tline=(char*)skipSpaces(tline);
   while ((*tline!=EOL1)&&(*tline!=EOL2)&&(*tline)) tline++;
   *tline='\0';
   return t;
}

initRandom()

void initRandom

(

int

i = 0

)

Definition at line 101 of file tools.cpp.

{
    if (i) { mysrand=i; return; }
    mysrand=(unsigned long) (clock());
}

isEmpty()

char isEmpty

(

const char *

line

)

Definition at line 165 of file tools.cpp.

{
    line=skipSpaces(line);
    return (*line==0);
}

isInt()

double isInt ( const double  a )

inline

Definition at line 64 of file tools.h.

{
    return condorAbs(a-floor( a + 0.5 ))<1e-4;
}

loadFile()

char* loadFile

(

FILE *

f

)

Definition at line 265 of file tools.cpp.

{
//    FILE *f=fopen(filename,"rb");
    fseek(f,0,SEEK_END);
    int l=ftell(f);
    fseek(f,0,SEEK_SET);
    char *buf=(char*)malloc(l+1);
    if (fread(buf,l,1,f) != l) {
        std::cerr << "error while loading file (no more info available)" << std::endl; exit(255); // FIXME implement properly
    }
    fclose(f);
    buf[l]=0;
    return buf;
}

mmax() [1/3]

double mmax ( const double  t1, const double  t2  )

inline

Definition at line 84 of file tools.h.

{
    return t1 > t2 ? t1 : t2;
}

mmax() [2/3]

unsigned mmax ( const unsigned  t1, const unsigned  t2  )

inline

Definition at line 89 of file tools.h.

{
    return t1 > t2 ? t1 : t2;
}

mmax() [3/3]

int mmax ( int  t1, int  t2  )

inline

Definition at line 94 of file tools.h.

{
    return t1 > t2 ? t1 : t2;
}

mmin() [1/3]

double mmin ( const double  t1, const double  t2  )

inline

Definition at line 69 of file tools.h.

{
    return t1 < t2 ? t1 : t2;
}

mmin() [2/3]

unsigned mmin ( const unsigned  t1, const unsigned  t2  )

inline

Definition at line 74 of file tools.h.

{
    return t1 < t2 ? t1 : t2;
}

mmin() [3/3]

int mmin ( const int  t1, const int  t2  )

inline

Definition at line 79 of file tools.h.

{
    return t1 < t2 ? t1 : t2;
}

rand1()

double rand1

(

)

Definition at line 88 of file tools.cpp.

{
    mysrand=1664525*mysrand+1013904223L;
    double r=((double)mysrand)/4294967297.0;
    //if (r>.52)
    //{
    //    printf("whoups\n");
    //}
    return r;

}

removeAllEOL()

char* removeAllEOL

(

char *

t

)

Definition at line 246 of file tools.cpp.

{
    char *t2=t;
    while (*t)
    {
        if ((*t=='\r')||(*t=='\n')) *t=' '; 
        t++;
    }
    return t2;
}

removeQuotes()

char* removeQuotes

(

char *

t

)

Definition at line 194 of file tools.cpp.

{
    if ((*t=='\'')||(*t=='"'))
    {
        t[strlen(t)-1]=0;
        return t+1;
    }
    return t;
}

sign() [1/2]

double sign ( const double  a )

inline

Definition at line 52 of file tools.h.

{
    return a<0?-1:1;
}

sign() [2/2]

double sign ( const double  t1, const double  t2  )

inline

Definition at line 58 of file tools.h.

{
    if(t2>=0) return condorAbs(t1);
    return -condorAbs(t1);
}

skipLine()

const char* skipLine

(

const char *

l

)

Definition at line 257 of file tools.cpp.

{
    while ((*t!='\r')&&(*t!='\n')&&(*t)) t++; 
    if (*t=='\n') { if (*(t+1)=='\r') t++; }
    else if (*t=='\r') { if (*(t+1)=='\n') t++; }
    return t+1;
}

skipSpaces()

const char* skipSpaces

(

const char *

l

)

Definition at line 171 of file tools.cpp.

{
    while ((*t==' ')||(*t=='\t')) t++;
    return t;
}

sqr()

double sqr ( const double &  t )

inline

Definition at line 99 of file tools.h.

{
    return t*t;
}

stringDuplicate()

char* stringDuplicate

(

const char *

l

)

Definition at line 234 of file tools.cpp.

{
    int l=(int)strlen(line);
    // remove quotes:
    if ((*line=='\'')||(*line=='"')) { l-=2; line++; }
    char *t=(char*)malloc(l+1);
    memcpy(t,line,l);
    t[l]=0;
    GetRidOfTheEOL(t);
    return t;
}