blob: f0a214919d35832e5568bcd3a7cf4a8ffc81f0d5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
/********************* */
/*! \file constraints.cpp
** \verbatim
** Top contributors (to current version):
** Gereon Kremer
** This file is part of the CVC4 project.
** Copyright (c) 2009-2020 by the authors listed in the file AUTHORS
** in the top-level source directory and their institutional affiliations.
** All rights reserved. See the file COPYING in the top-level source
** directory for licensing information.\endverbatim
**
** \brief Implements a container for CAD constraints.
**
** Implements a container for CAD constraints.
**/
#include "theory/arith/nl/cad/constraints.h"
#ifdef CVC4_POLY_IMP
#include <algorithm>
#include "theory/arith/nl/poly_conversion.h"
#include "util/poly_util.h"
namespace CVC4 {
namespace theory {
namespace arith {
namespace nl {
namespace cad {
void Constraints::addConstraint(const poly::Polynomial& lhs,
poly::SignCondition sc,
Node n)
{
d_constraints.emplace_back(lhs, sc, n);
sortConstraints();
}
void Constraints::addConstraint(Node n)
{
auto c = as_poly_constraint(n, d_varMapper);
addConstraint(c.first, c.second, n);
sortConstraints();
}
const Constraints::ConstraintVector& Constraints::getConstraints() const
{
return d_constraints;
}
void Constraints::reset() { d_constraints.clear(); }
void Constraints::sortConstraints()
{
using Tpl = std::tuple<poly::Polynomial, poly::SignCondition, Node>;
std::sort(d_constraints.begin(),
d_constraints.end(),
[](const Tpl& at, const Tpl& bt) {
// Check if a is smaller than b
const poly::Polynomial& a = std::get<0>(at);
const poly::Polynomial& b = std::get<0>(bt);
bool ua = is_univariate(a);
bool ub = is_univariate(b);
if (ua != ub) return ua;
std::size_t tda = poly_utils::totalDegree(a);
std::size_t tdb = poly_utils::totalDegree(b);
if (tda != tdb) return tda < tdb;
return degree(a) < degree(b);
});
for (auto& c : d_constraints)
{
auto* p = std::get<0>(c).get_internal();
lp_polynomial_set_external(p);
}
}
} // namespace cad
} // namespace nl
} // namespace arith
} // namespace theory
} // namespace CVC4
#endif
|