---

MathUt.hpp

#ifndef _MathUt_H
#define _MathUt_H



#include <math.h>
#include "RComplex.hpp"

class MathUt
{
public:
        
//  ordering
static int sign(double a){ if(a==0.0) return 0.0; return (a >= 0.0 ? 1 : -1); }

static int minsign(double a, double b){return (fabs(a) <= fabs(b) ? sign(a) : sign(b));}

// static double min(double a, double b){ return (a < b ? a : b); } // predefined as macro
// static double max(double a, double b){ return (a > b ? a : b); } // predefined as macro

static double minEq(double a, double b){ return (a <= b ? a : b); }
static double maxEq(double a, double b){ return (a >= b ? a : b); }

static bool equal(double a, double b, double eps=1e-12){ return (a-eps < b && a+eps > b); }

//  minmod for difference computation
static double minmod(double a, double b)
        {
        return 0.5* (sign(a) + sign(b)) * (fabs(a) < fabs(b) ? fabs(a) : fabs(b));
        }

static double maxmod(double a, double b)
        {
        return 0.5* (sign(a) + sign(b)) * (fabs(a) > fabs(b) ? fabs(a) : fabs(b));
        }

static RComplex minmod(RComplex a, RComplex b)
        {
        return RComplex(minmod(a.a,b.a), minmod(a.b,b.b));
        }

};

#endif

Restoreinpaint sourceforge project `C++/Java Image Processing, Restoration, Inpainting Project'.

Bernard De Cuyper: Open Project Leader: Concept, design and development. Bernard De Cuyper & Eddy Fraiha 2002, 2003. Bernard De Cuyper 2004. Open and free, for friendly usage only. Modifications on Belgium ground of this piece of artistic work, by governement institutions or companies, must be notified to Bernard De Cuyper.