r493:81627fa1b007
Wed, 14 Jan 2009 16:55:29
[Not Reviewed]
0 comments (0 inline)
0
2
0
| ... | ... |
@@ -588,29 +588,29 @@ |
| 588 | 588 |
const Value &lowerBound() const { return _lb; }
|
| 589 | 589 |
///Reference to the upper bound. |
| 590 | 590 |
|
| 591 | 591 |
///\return |
| 592 | 592 |
///- \ref INF "INF": the constraint is upper unbounded. |
| 593 | 593 |
///- \ref NaN "NaN": upper bound has not been set. |
| 594 | 594 |
///- finite number: the upper bound |
| 595 | 595 |
Value &upperBound() { return _ub; }
|
| 596 | 596 |
///The const version of \ref upperBound() |
| 597 | 597 |
const Value &upperBound() const { return _ub; }
|
| 598 | 598 |
///Is the constraint lower bounded? |
| 599 | 599 |
bool lowerBounded() const {
|
| 600 |
return _lb != -INF && ! |
|
| 600 |
return _lb != -INF && !isnan(_lb); |
|
| 601 | 601 |
} |
| 602 | 602 |
///Is the constraint upper bounded? |
| 603 | 603 |
bool upperBounded() const {
|
| 604 |
return _ub != INF && ! |
|
| 604 |
return _ub != INF && !isnan(_ub); |
|
| 605 | 605 |
} |
| 606 | 606 |
|
| 607 | 607 |
}; |
| 608 | 608 |
|
| 609 | 609 |
///Linear expression of rows |
| 610 | 610 |
|
| 611 | 611 |
///This data structure represents a column of the matrix, |
| 612 | 612 |
///thas is it strores a linear expression of the dual variables |
| 613 | 613 |
///(\ref Row "Row"s). |
| 614 | 614 |
/// |
| 615 | 615 |
///There are several ways to access and modify the contents of this |
| 616 | 616 |
///container. |
| ... | ... |
@@ -1657,61 +1657,61 @@ |
| 1657 | 1657 |
inline LpBase::Constr operator==(const LpBase::Expr &e, |
| 1658 | 1658 |
const LpBase::Expr &f) {
|
| 1659 | 1659 |
return LpBase::Constr(0, f - e, 0); |
| 1660 | 1660 |
} |
| 1661 | 1661 |
|
| 1662 | 1662 |
///Create constraint |
| 1663 | 1663 |
|
| 1664 | 1664 |
///\relates LpBase::Constr |
| 1665 | 1665 |
/// |
| 1666 | 1666 |
inline LpBase::Constr operator<=(const LpBase::Value &n, |
| 1667 | 1667 |
const LpBase::Constr &c) {
|
| 1668 | 1668 |
LpBase::Constr tmp(c); |
| 1669 |
LEMON_ASSERT( |
|
| 1669 |
LEMON_ASSERT(isnan(tmp.lowerBound()), "Wrong LP constraint"); |
|
| 1670 | 1670 |
tmp.lowerBound()=n; |
| 1671 | 1671 |
return tmp; |
| 1672 | 1672 |
} |
| 1673 | 1673 |
///Create constraint |
| 1674 | 1674 |
|
| 1675 | 1675 |
///\relates LpBase::Constr |
| 1676 | 1676 |
/// |
| 1677 | 1677 |
inline LpBase::Constr operator<=(const LpBase::Constr &c, |
| 1678 | 1678 |
const LpBase::Value &n) |
| 1679 | 1679 |
{
|
| 1680 | 1680 |
LpBase::Constr tmp(c); |
| 1681 |
LEMON_ASSERT( |
|
| 1681 |
LEMON_ASSERT(isnan(tmp.upperBound()), "Wrong LP constraint"); |
|
| 1682 | 1682 |
tmp.upperBound()=n; |
| 1683 | 1683 |
return tmp; |
| 1684 | 1684 |
} |
| 1685 | 1685 |
|
| 1686 | 1686 |
///Create constraint |
| 1687 | 1687 |
|
| 1688 | 1688 |
///\relates LpBase::Constr |
| 1689 | 1689 |
/// |
| 1690 | 1690 |
inline LpBase::Constr operator>=(const LpBase::Value &n, |
| 1691 | 1691 |
const LpBase::Constr &c) {
|
| 1692 | 1692 |
LpBase::Constr tmp(c); |
| 1693 |
LEMON_ASSERT( |
|
| 1693 |
LEMON_ASSERT(isnan(tmp.upperBound()), "Wrong LP constraint"); |
|
| 1694 | 1694 |
tmp.upperBound()=n; |
| 1695 | 1695 |
return tmp; |
| 1696 | 1696 |
} |
| 1697 | 1697 |
///Create constraint |
| 1698 | 1698 |
|
| 1699 | 1699 |
///\relates LpBase::Constr |
| 1700 | 1700 |
/// |
| 1701 | 1701 |
inline LpBase::Constr operator>=(const LpBase::Constr &c, |
| 1702 | 1702 |
const LpBase::Value &n) |
| 1703 | 1703 |
{
|
| 1704 | 1704 |
LpBase::Constr tmp(c); |
| 1705 |
LEMON_ASSERT( |
|
| 1705 |
LEMON_ASSERT(isnan(tmp.lowerBound()), "Wrong LP constraint"); |
|
| 1706 | 1706 |
tmp.lowerBound()=n; |
| 1707 | 1707 |
return tmp; |
| 1708 | 1708 |
} |
| 1709 | 1709 |
|
| 1710 | 1710 |
///Addition |
| 1711 | 1711 |
|
| 1712 | 1712 |
///\relates LpBase::DualExpr |
| 1713 | 1713 |
/// |
| 1714 | 1714 |
inline LpBase::DualExpr operator+(const LpBase::DualExpr &a, |
| 1715 | 1715 |
const LpBase::DualExpr &b) {
|
| 1716 | 1716 |
LpBase::DualExpr tmp(a); |
| 1717 | 1717 |
tmp+=b; |
| ... | ... |
@@ -46,18 +46,27 @@ |
| 46 | 46 |
const long double LN10 = 2.3025850929940456840179914546843642L; |
| 47 | 47 |
/// pi |
| 48 | 48 |
const long double PI = 3.1415926535897932384626433832795029L; |
| 49 | 49 |
/// pi/2 |
| 50 | 50 |
const long double PI_2 = 1.5707963267948966192313216916397514L; |
| 51 | 51 |
/// pi/4 |
| 52 | 52 |
const long double PI_4 = 0.7853981633974483096156608458198757L; |
| 53 | 53 |
/// sqrt(2) |
| 54 | 54 |
const long double SQRT2 = 1.4142135623730950488016887242096981L; |
| 55 | 55 |
/// 1/sqrt(2) |
| 56 | 56 |
const long double SQRT1_2 = 0.7071067811865475244008443621048490L; |
| 57 | 57 |
|
| 58 |
///Check whether the parameter is NaN or not |
|
| 59 |
|
|
| 60 |
///This function checks whether the parameter is NaN or not. |
|
| 61 |
///Is should be equivalent with std::isnan(), but it is not |
|
| 62 |
///provided by all compilers. |
|
| 63 |
inline bool isnan(double v) |
|
| 64 |
{
|
|
| 65 |
return v!=v; |
|
| 66 |
} |
|
| 58 | 67 |
|
| 59 | 68 |
/// @} |
| 60 | 69 |
|
| 61 | 70 |
} //namespace lemon |
| 62 | 71 |
|
| 63 | 72 |
#endif //LEMON_TOLERANCE_H |
0 comments (0 inline)