FWMath.h

/* -*- mode:c++ -*- ********************************************************
 * file:        FWMath.h
 *
 * author:      Christian Frank
 *
 * copyright:   (C) 2004 Telecommunication Networks Group (TKN) at
 *              Technische Universitaet Berlin, Germany.
 *
 *              This program is free software; you can redistribute it
 *              and/or modify it under the terms of the GNU General Public
 *              License as published by the Free Software Foundation; either
 *              version 2 of the License, or (at your option) any later
 *              version.
 *              For further information see file COPYING
 *              in the top level directory
 ***************************************************************************
 * part of:     framework implementation developed by tkn
 **************************************************************************/

#ifndef FWMATH_H
#define FWMATH_H

//
// Support functions for mathematical operations
//

#include <math.h>


/* windows math.h doesn't define the PI variable so we have to do it
   by hand*/
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif

/* Constant for comparing doubles. Two doubles at most epsilon apart
   are declared equal.*/
#ifndef EPSILON
#define EPSILON 0.001
#endif


class FWMath {

public:

    static double mod(double dividend, double divisor) {
        double i;
        return modf(dividend/divisor, &i)*divisor;
    }

    static double div(double dividend, double divisor) {
        double i;
        double f;
        f=modf(dividend/divisor, &i);
        return i;
    }

    static bool close(double one, double two) {
        return fabs(one-two)<EPSILON;
    }

    static int stepfunction(double i) { return (i>EPSILON) ? 1 : close(i,0) ? 0 :-1; };


    static int sign(double i) { return (i>=0)? 1 : -1; };

    static int round(double d) { return (int)(ceil(d-0.5)); }

    static double max(double a, double b) { return (a<b)? b : a; }

    static double dBm2mW(double dBm){
        return pow(10.0, dBm/10.0);
    }

    static inline double torDist (double x1, double x2, double y1, double y2) {
        return (x1-x2) * (x1-x2) + (y1-y2) * (y1-y2);
    };
};

#endif

---