30 #ifndef _SPXDEFINES_H_ 31 #define _SPXDEFINES_H_ 42 #define SOPLEX_VERSION 301 43 #define SOPLEX_SUBVERSION 0 44 #define SOPLEX_APIVERSION 0 45 #define SOPLEX_COPYRIGHT "Copyright (c) 1996-2017 Konrad-Zuse-Zentrum fuer Informationstechnik Berlin (ZIB)" 62 #if defined (NDEBUG) && defined (WITH_WARNINGS) 63 #define ASSERT_WARN( prefix, expr ) \ 68 << " failed assertion on line " << __LINE__ \ 69 << " in file " << __FILE__ << ": " \ 73 #else // just a normal assert 74 #define ASSERT_WARN( prefix, expr ) ( assert( expr ) ) 90 #ifdef DISABLE_VERBOSITY 91 #define DO_WITH_TMP_VERBOSITY( verbosity, spxout, do_something ) {} 92 #define DO_WITH_ERR_VERBOSITY( do_something ) {} 94 #define DO_WITH_TMP_VERBOSITY( verbosity, spxout, do_something ) \ 96 if( &spxout != NULL ) \ 98 if( verbosity <= spxout.getVerbosity() ) \ 100 const SPxOut::Verbosity old_verbosity = spxout.getVerbosity(); \ 101 spxout.setVerbosity( verbosity ); \ 103 spxout.setVerbosity( old_verbosity ); \ 107 #define DO_WITH_ERR_VERBOSITY( do_something ) { do_something; } 111 #define MSG_ERROR(x) { DO_WITH_ERR_VERBOSITY( x ) } 113 #define MSG_WARNING(spxout, x) { DO_WITH_TMP_VERBOSITY( SPxOut::WARNING, spxout, x ) } 115 #define MSG_INFO1(spxout, x) { DO_WITH_TMP_VERBOSITY( SPxOut::INFO1, spxout, x ) } 117 #define MSG_INFO2(spxout, x) { DO_WITH_TMP_VERBOSITY( SPxOut::INFO2, spxout, x ) } 119 #define MSG_INFO3(spxout, x) { DO_WITH_TMP_VERBOSITY( SPxOut::INFO3, spxout, x ) } 121 extern bool msginconsistent(
const char* name,
const char* file,
int line);
123 #define MSGinconsistent(name) msginconsistent(name, __FILE__, __LINE__) 125 #if defined(SOPLEX_DEBUG) 127 #define MSG_DEBUG(x) { x; } 138 #if !defined(THREADLOCAL) 139 #if defined(_MSC_VER) && _MSC_VER < 1900 142 #define THREADLOCAL thread_local 151 #ifdef WITH_LONG_DOUBLE 154 typedef long double Real;
158 #define REAL_FORMAT "Lf" 161 #ifndef DEFAULT_BND_VIOL 162 #define DEFAULT_BND_VIOL 1e-12L 165 #ifndef DEFAULT_EPS_ZERO 166 #define DEFAULT_EPS_ZERO 1e-28L 169 #ifndef DEFAULT_EPS_FACTOR 170 #define DEFAULT_EPS_FACTOR 1e-30L 173 #ifndef DEFAULT_EPS_UPDATE 174 #define DEFAULT_EPS_UPDATE 1e-26L 176 #ifndef DEFAULT_EPS_PIVOT 177 #define DEFAULT_EPS_PIVOT 1e-20L 180 #define DEFAULT_INFINITY 1e100L 191 #define REAL_FORMAT "f" 194 #ifndef DEFAULT_BND_VIOL 195 #define DEFAULT_BND_VIOL 1e-1 198 #ifndef DEFAULT_EPS_ZERO 199 #define DEFAULT_EPS_ZERO 1e-7 201 #ifndef DEFAULT_EPS_FACTOR 202 #define DEFAULT_EPS_FACTOR 1e-7 204 #ifndef DEFAULT_EPS_UPDATE 205 #define DEFAULT_EPS_UPDATE 1e-6 207 #ifndef DEFAULT_EPS_PIVOT 208 #define DEFAULT_EPS_PIVOT 1e-6 210 #define DEFAULT_INFINITY 1e100 219 #define REAL_FORMAT "lf" 222 #ifndef DEFAULT_BND_VIOL 223 #define DEFAULT_BND_VIOL 1e-6 226 #ifndef DEFAULT_EPS_ZERO 227 #define DEFAULT_EPS_ZERO 1e-16 229 #ifndef DEFAULT_EPS_FACTOR 230 #define DEFAULT_EPS_FACTOR 1e-20 232 #ifndef DEFAULT_EPS_UPDATE 233 #define DEFAULT_EPS_UPDATE 1e-16 235 #ifndef DEFAULT_EPS_PIVOT 236 #define DEFAULT_EPS_PIVOT 1e-10 238 #define DEFAULT_INFINITY 1e100 240 #endif // !WITH_FLOAT 241 #endif // !WITH_LONG_DOUBLE 243 #define MAXIMUM(x,y) ((x)>(y) ? (x) : (y)) 244 #define MINIMUM(x,y) ((x)<(y) ? (x) : (y)) 289 #ifdef WITH_LONG_DOUBLE 291 inline Real
spxAbs(Real a)
303 #ifndef SOPLEX_LEGACY 306 return nextafterl(x,y);
311 inline Real
spxLdexp(Real x,
int exp)
313 return ldexpl(x,exp);
317 inline Real
spxFrexp(Real y,
int* exp)
319 return frexpl(y, exp);
335 #ifndef SOPLEX_LEGACY 339 return nextafter(x,y);
341 return _nextafter(x,y);
355 return frexp(y, exp);
362 const Real absa =
spxAbs(a);
363 const Real absb =
spxAbs(b);
365 return absa > absb ? absa : absb;
371 return (a - b) / (
maxAbs(a, b) > 1.0 ?
maxAbs(a, b) : 1.0);
377 return spxAbs(a - b) <= eps;
383 return spxAbs(a - b) > eps;
389 return (a - b) < -eps;
395 return (a - b) < eps;
401 return (a - b) > eps;
407 return (a - b) > -eps;
459 #endif // _SPXDEFINES_H_ Rational spxAbs(const Rational &r)
Absolute.
bool isNotZero(Real a, Real eps=Param::epsilon())
returns true iff |a| > eps
bool LTrel(Real a, Real b, Real eps=Param::epsilon())
returns true iff relDiff(a,b) <= -eps
bool GE(Real a, Real b, Real eps=Param::epsilon())
returns true iff a >= b + eps
static void setEpsilon(Real eps)
THREADLOCAL const Real infinity
Real maxAbs(Real a, Real b)
returns max(|a|,|b|)
bool msginconsistent(const char *name, const char *file, int line)
bool LE(Real a, Real b, Real eps=Param::epsilon())
returns true iff a <= b + eps
bool LT(Real a, Real b, Real eps=Param::epsilon())
returns true iff a < b + eps
bool NE(Real a, Real b, Real eps=Param::epsilon())
returns true iff |a-b| > eps
Real spxFrexp(Real y, int *exp)
Real spxLdexp(Real x, int exp)
returns x * 2^exp
bool GT(Real a, Real b, Real eps=Param::epsilon())
returns true iff a > b + eps
static void setEpsilonPivot(Real eps)
static THREADLOCAL Real s_epsilon_pivot
epsilon for pivot zero tolerance in factorization
static Real epsilonUpdate()
Real spxSqrt(Real a)
returns square root
bool GTrel(Real a, Real b, Real eps=Param::epsilon())
returns true iff relDiff(a,b) > eps
bool GErel(Real a, Real b, Real eps=Param::epsilon())
returns true iff relDiff(a,b) > -eps
bool EQ(Real a, Real b, Real eps=Param::epsilon())
returns true iff |a-b| <= eps
static THREADLOCAL Real s_epsilon
default allowed additive zero: 1.0 + EPS_ZERO == 1.0
Everything should be within this namespace.
Real relDiff(Real a, Real b)
returns (a-b) / max(|a|,|b|,1.0)
bool EQrel(Real a, Real b, Real eps=Param::epsilon())
returns true iff |relDiff(a,b)| <= eps
Real spxNextafter(Real x, Real y)
static void setEpsilonUpdate(Real eps)
bool NErel(Real a, Real b, Real eps=Param::epsilon())
returns true iff |relDiff(a,b)| > eps
static THREADLOCAL Real s_epsilon_update
epsilon for factorization update
static Real epsilonPivot()
bool LErel(Real a, Real b, Real eps=Param::epsilon())
returns true iff relDiff(a,b) <= eps
static void setEpsilonFactorization(Real eps)
bool isZero(Real a, Real eps=Param::epsilon())
returns true iff |a| <= eps
static Real epsilonFactorization()
static THREADLOCAL Real s_epsilon_factorization
epsilon for factorization
#define THREADLOCAL
SOPLEX_DEBUG.