utils.h

Go to the documentation of this file.

 
#ifndef UTILS_H
#define UTILS_H
 
#include "mycomplex.h"
#include "def.h"
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <cmath>
#include <stdio.h>
#include <string>
#include <exception>
#include <csignal>
#include <unistd.h>
#include <stdlib.h>
#include <stdio.h>
#include <iostream>
#include <signal.h>
#include <dirent.h>
 
using namespace std;
using namespace softsusy;
 
inline double zeroSqrt(double f){ if (f > 0.) return sqrt(f); 
  else return EPSTOL; 
}
double frexp(const Complex & c, int * i);
void FPE_ExceptionHandler(int nErrType);
int theta(double a);
void outputCharacteristics(int);
inline double sqr(double a) { return a * a; }
inline double maximum(double a, double b) { return ((a > b) ? a : b); }
inline double minimum(double a, double b) { return ((a < b) ? a : b); }
inline int minimum(int a, int b) { return ((a < b) ? a : b); }
double toleranceCheck(double sTin, double sTout);
 
inline double minimum(double a, double b, double c) {
   if(a <= b && a <= c) return a;
   else if(b <= a && b <= c) return b;
   else return c;
}
 
inline double maximum(double a, double b, double c) {
   if(a >= b && a >= c) return a;
   else if(b >= a && b >= c) return b;
   else return c;
}
 
inline double middle(double a, double b, double c) {
   if((a >= b && a <= c) ||(a >= c && a <= b)   ) return a;
   else if((b >= a && b <= c) ||(b >= c && b <= a) ) return b;
   else  return c;
}
 
inline int massorder(double & a, double & b, double & c) {
   double anew, bnew, cnew;
   anew = minimum(a,b,c);
   cnew= maximum(a,b,c);
   bnew = middle(a,b,c);
 
   a = anew;
   b = bnew;
   c = cnew;
   return 0;
}
 
bool close(double m1, double m2, double tol);
 
bool close(Complex m1, Complex m2, double tol);
 
inline double sign(double a, double b) { return ((b) >= 0.0 ? 
                                                 fabs(a) : -fabs(a)); }
 
inline int sgn(double x)
{ return (x >= 0.0 ? 1 : -1); }
 
void printRow(ostream & o, double x);
 
bool testNan(double f);
 
// For DEBUGING
#ifdef DEBUG
#define dout cout
#else
#define dout if(0)cout
#endif
 
 
#endif