tools.cpp File Reference

#include <stdlib.h>
#include <iostream>
#include <time.h>
#include <string.h>
#include "tools.h"
#include <unistd.h>

Include dependency graph for tools.cpp:

Go to the source code of this file.

Macros
#define	EOL1   13
#define	EOL2   10

Functions
unsigned long	choose (unsigned n, unsigned k)
double	rand1 ()
void	initRandom (int i)
void	error (char *s)
double	euclidianNorm (int i, double *xp)
char	isEmpty (const char *line)
const char *	skipSpaces (const char *t)
char	isEmptyline (const char *line)
char *	GetRidOfTheEOL (char *tline)
char *	removeQuotes (char *t)
char **	getNameTable (const char *line, int *ncolumn)
char *	stringDuplicate (const char *line)
char *	removeAllEOL (char *t)
const char *	skipLine (const char *t)
char *	loadFile (FILE *f)
void	deleteFile (const char *line)

Variables
unsigned long	mysrand

Macro Definition Documentation

EOL1

#define EOL1   13

Definition at line 162 of file tools.cpp.

EOL2

#define EOL2   10

Definition at line 163 of file tools.cpp.

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);
}

deleteFile()

void deleteFile

(

const char *

line

)

Definition at line 280 of file tools.cpp.

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

error()

void error

(

char *

s

)

Definition at line 107 of file tools.cpp.

{
   printf("Error due to %s.", s);
   getchar(); exit(255);
}

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

)

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);
}

isEmptyline()

char isEmptyline

(

const char *

line

)

Definition at line 177 of file tools.cpp.

{
    if (*line==';') return 1;
    line=skipSpaces(line);
    if ((*line==EOL1)||(*line==EOL2)||(*line=='\0')) return 1;
    return 0;
}

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;
}

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;
}

skipLine()

const char* skipLine

(

const char *

t

)

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 *

t

)

Definition at line 171 of file tools.cpp.

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

stringDuplicate()

char* stringDuplicate

(

const char *

line

)

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;
}

Variable Documentation

mysrand

unsigned long mysrand

Definition at line 87 of file tools.cpp.