Revision 69e4dffe991aafe28a7457fbe7cc3978a8074020 authored by Andrew Adams on 21 February 2021, 17:31:49 UTC, committed by Andrew Adams on 21 February 2021, 17:31:49 UTC
Simplify_EQ.cpp
#include "Simplify_Internal.h"
namespace Halide {
namespace Internal {
Expr Simplify::visit(const EQ *op, ExprInfo *bounds) {
if (truths.count(op)) {
return const_true(op->type.lanes());
} else if (falsehoods.count(op)) {
return const_false(op->type.lanes());
}
if (!may_simplify(op->a.type())) {
Expr a = mutate(op->a, nullptr);
Expr b = mutate(op->b, nullptr);
if (a.same_as(op->a) && b.same_as(op->b)) {
return op;
} else {
return EQ::make(a, b);
}
}
if (op->a.type().is_bool()) {
Expr a = mutate(op->a, nullptr);
Expr b = mutate(op->b, nullptr);
if (should_commute(a, b)) {
std::swap(a, b);
}
const int lanes = op->type.lanes();
auto rewrite = IRMatcher::rewriter(IRMatcher::eq(a, b), op->type);
if (rewrite(x == 1, x)) {
return rewrite.result;
} else if (rewrite(x == 0, !x)) {
return mutate(rewrite.result, bounds);
} else if (rewrite(x == x, const_true(lanes))) {
return rewrite.result;
} else if (a.same_as(op->a) && b.same_as(op->b)) {
return op;
} else {
return EQ::make(a, b);
}
}
ExprInfo delta_bounds;
Expr delta = mutate(op->a - op->b, &delta_bounds);
const int lanes = op->type.lanes();
// If the delta is 0, then it's just x == x
if (is_const_zero(delta)) {
return const_true(lanes);
}
// Attempt to disprove using bounds analysis
if (delta_bounds.min_defined && delta_bounds.min > 0) {
return const_false(lanes);
}
if (delta_bounds.max_defined && delta_bounds.max < 0) {
return const_false(lanes);
}
// Attempt to disprove using modulus remainder analysis
if (delta_bounds.alignment.remainder != 0) {
return const_false(lanes);
}
auto rewrite = IRMatcher::rewriter(IRMatcher::eq(delta, 0), op->type, delta.type());
if (rewrite(broadcast(x, c0) == 0, broadcast(x == 0, c0)) ||
(no_overflow(delta.type()) && rewrite(x * y == 0, (x == 0) || (y == 0))) ||
rewrite(select(x, 0, y) == 0, x || (y == 0)) ||
rewrite(select(x, c0, y) == 0, !x && (y == 0), c0 != 0) ||
rewrite(select(x, y, 0) == 0, !x || (y == 0)) ||
rewrite(select(x, y, c0) == 0, x && (y == 0), c0 != 0) ||
rewrite(select(x, c0, y) + c1 == 0, x || (y == fold(-c1)), c0 + c1 == 0) ||
rewrite(select(x, y, c0) + c1 == 0, !x || (y == fold(-c1)), c0 + c1 == 0) ||
rewrite(select(x, c0, y) + c1 == 0, !x && (y == fold(-c1)), c0 + c1 != 0) ||
rewrite(select(x, y, c0) + c1 == 0, x && (y == fold(-c1)), c0 + c1 != 0) ||
rewrite(max(x, y) - y == 0, x <= y) ||
rewrite(min(x, y) - y == 0, y <= x) ||
rewrite(max(y, x) - y == 0, x <= y) ||
rewrite(min(y, x) - y == 0, y <= x) ||
rewrite(y - max(x, y) == 0, x <= y) ||
rewrite(y - min(x, y) == 0, y <= x) ||
rewrite(y - max(y, x) == 0, x <= y) ||
rewrite(y - min(y, x) == 0, y <= x) ||
rewrite(max(x, c0) + c1 == 0, x == fold(-c1), c0 + c1 < 0) ||
rewrite(min(x, c0) + c1 == 0, x == fold(-c1), c0 + c1 > 0) ||
rewrite(max(x, c0) + c1 == 0, false, c0 + c1 > 0) ||
rewrite(min(x, c0) + c1 == 0, false, c0 + c1 < 0) ||
rewrite(max(x, c0) + c1 == 0, x <= c0, c0 + c1 == 0) ||
rewrite(min(x, c0) + c1 == 0, c0 <= x, c0 + c1 == 0) ||
// Special case the above where c1 == 0
rewrite(max(x, c0) == 0, x == 0, c0 < 0) ||
rewrite(min(x, c0) == 0, x == 0, c0 > 0) ||
rewrite(max(x, c0) == 0, false, c0 > 0) ||
rewrite(min(x, c0) == 0, false, c0 < 0) ||
rewrite(max(x, 0) == 0, x <= 0) ||
rewrite(min(x, 0) == 0, 0 <= x) ||
false) {
return mutate(rewrite.result, bounds);
}
if (rewrite(c0 == 0, fold(c0 == 0)) ||
rewrite((x - y) + c0 == 0, x == y + fold(-c0)) ||
rewrite(x + c0 == 0, x == fold(-c0)) ||
rewrite(c0 - x == 0, x == c0) ||
rewrite(x - y == 0, x == y) ||
rewrite(x == 0, x == 0)) {
const EQ *eq = rewrite.result.as<EQ>();
if (eq &&
eq->a.same_as(op->a) &&
eq->b.same_as(op->b)) {
return op;
} else {
return rewrite.result;
}
}
// Unreachable. That last rewrite catches everything and
// constructs delta == 0 for us.
return Expr();
}
// ne redirects to not eq
Expr Simplify::visit(const NE *op, ExprInfo *bounds) {
if (!may_simplify(op->a.type())) {
Expr a = mutate(op->a, nullptr);
Expr b = mutate(op->b, nullptr);
if (a.same_as(op->a) && b.same_as(op->b)) {
return op;
} else {
return NE::make(a, b);
}
}
Expr mutated = mutate(Not::make(EQ::make(op->a, op->b)), bounds);
if (const NE *ne = mutated.as<NE>()) {
if (ne->a.same_as(op->a) && ne->b.same_as(op->b)) {
return op;
}
}
return mutated;
}
} // namespace Internal
} // namespace Halide
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...