diff --git a/lib/proof_utils.cpp b/lib/proof_utils.cpp new file mode 100644 index 000000000..35e49c4e1 --- /dev/null +++ b/lib/proof_utils.cpp @@ -0,0 +1,389 @@ +#include "dl_util.h" +#include "proof_utils.h" +#include "ast_smt2_pp.h" + +class reduce_hypotheses { + typedef obj_hashtable expr_set; + ast_manager& m; + expr_ref_vector m_refs; + obj_map m_cache; + obj_map m_units; + ptr_vector m_units_trail; + unsigned_vector m_limits; + obj_map m_hypmap; + ptr_vector m_hyprefs; + ptr_vector m_literals; + + void reset() { + m_refs.reset(); + m_cache.reset(); + m_units.reset(); + m_units_trail.reset(); + m_limits.reset(); + std::for_each(m_hyprefs.begin(), m_hyprefs.end(), delete_proc()); + m_hypmap.reset(); + m_hyprefs.reset(); + m_literals.reset(); + } + + void push() { + m_limits.push_back(m_units_trail.size()); + } + + void pop() { + unsigned sz = m_limits.back(); + while (m_units_trail.size() > sz) { + m_units.remove(m_units_trail.back()); + m_units_trail.pop_back(); + } + m_limits.pop_back(); + } + + void get_literals(expr* clause) { + m_literals.reset(); + if (m.is_or(clause)) { + m_literals.append(to_app(clause)->get_num_args(), to_app(clause)->get_args()); + } + else { + m_literals.push_back(clause); + } + } + + void add_hypotheses(proof* p) { + expr_set* hyps = 0; + bool inherited = false; + if (p->get_decl_kind() == PR_HYPOTHESIS) { + hyps = alloc(expr_set); + hyps->insert(m.get_fact(p)); + m_hyprefs.push_back(hyps); + } + else { + for (unsigned i = 0; i < m.get_num_parents(p); ++i) { + expr_set* hyps1 = m_hypmap.find(m.get_parent(p, i)); + if (hyps1) { + if (!hyps) { + hyps = hyps1; + inherited = true; + continue; + } + if (inherited) { + hyps = alloc(expr_set,*hyps); + m_hyprefs.push_back(hyps); + inherited = false; + } + datalog::set_union(*hyps, *hyps1); + } + } + } + m_hypmap.insert(p, hyps); + } + + expr_ref complement_lit(expr* e) { + expr* e1; + if (m.is_not(e, e1)) { + return expr_ref(e1, m); + } + else { + return expr_ref(m.mk_not(e), m); + } + } + + bool in_hypotheses(expr* e, expr_set* hyps) { + if (!hyps) { + return false; + } + expr_ref not_e = complement_lit(e); + return hyps->contains(not_e); + } + + bool contains_hypothesis(proof* p) { + ptr_vector todo; + ast_mark visit; + todo.push_back(p); + while (!todo.empty()) { + p = todo.back(); + todo.pop_back(); + if (visit.is_marked(p)) { + continue; + } + visit.mark(p, true); + if (PR_HYPOTHESIS == p->get_decl_kind()) { + return true; + } + for (unsigned i = 0; i < m.get_num_parents(p); ++i) { + todo.push_back(m.get_parent(p, i)); + } + } + return false; + } + + bool is_closed(proof* p) { + expr_set* hyps = m_hypmap.find(p); + return !hyps || hyps->empty(); + } + +public: + reduce_hypotheses(ast_manager& m): m(m), m_refs(m) {} + + void operator()(proof_ref& pr) { + proof_ref tmp(m); + tmp = pr; + elim(pr); + reset(); + CTRACE("proof_utils", contains_hypothesis(pr), + tout << "Contains hypothesis:\n"; + tout << mk_ismt2_pp(tmp, m) << "\n====>\n"; + tout << mk_ismt2_pp(pr, m) << "\n";); + + } + + void elim(proof_ref& p) { + proof_ref tmp(m); + proof* result = p.get(); + if (m_cache.find(p, result)) { + p = result; + return; + } + switch(p->get_decl_kind()) { + case PR_HYPOTHESIS: + if (!m_units.find(m.get_fact(p), result)) { + result = p.get(); + } + add_hypotheses(result); + break; + case PR_LEMMA: { + SASSERT(m.get_num_parents(p) == 1); + tmp = m.get_parent(p, 0); + elim(tmp); + get_literals(m.get_fact(p)); + expr_set* hyps = m_hypmap.find(tmp); + expr_set* new_hyps = 0; + if (hyps) { + new_hyps = alloc(expr_set, *hyps); + } + for (unsigned i = 0; i < m_literals.size(); ++i) { + expr* e = m_literals[i]; + if (!in_hypotheses(e, hyps)) { + m_literals[i] = m_literals.back(); + m_literals.pop_back(); + --i; + } + else { + SASSERT(new_hyps); + expr_ref not_e = complement_lit(e); + SASSERT(new_hyps->contains(not_e)); + new_hyps->remove(not_e); + } + } + if (m_literals.empty()) { + result = tmp; + } + else { + expr_ref clause(m); + if (m_literals.size() == 1) { + clause = m_literals[0]; + } + else { + clause = m.mk_or(m_literals.size(), m_literals.c_ptr()); + } + tmp = m.mk_lemma(tmp, clause); + m_refs.push_back(tmp); + result = tmp; + } + if (new_hyps && new_hyps->empty()) { + dealloc(new_hyps); + new_hyps = 0; + } + m_hypmap.insert(result, new_hyps); + m_hyprefs.push_back(new_hyps); + TRACE("proof_utils", + { + tout << "New lemma: " << mk_pp(m.get_fact(p), m) + << "\n==>\n" + << mk_pp(m.get_fact(result), m) << "\n"; + if (hyps) { + expr_set::iterator it = hyps->begin(), end = hyps->end(); + for (; it != end; ++it) { + tout << "Hypothesis: " << mk_pp(*it, m) << "\n"; + } + } + }); + } + break; + } + case PR_UNIT_RESOLUTION: { + proof_ref_vector parents(m); + parents.push_back(m.get_parent(p, 0)); + push(); + bool found_false = false; + for (unsigned i = 1; i < m.get_num_parents(p); ++i) { + tmp = m.get_parent(p, i); + elim(tmp); + if (m.is_false(m.get_fact(tmp))) { + result = tmp; + found_false = true; + break; + } + SASSERT(m.get_fact(tmp) == m.get_fact(m.get_parent(p, i))); + parents.push_back(tmp); + if (is_closed(tmp) && !m_units.contains(m.get_fact(tmp))) { + m_units.insert(m.get_fact(tmp), tmp); + m_units_trail.push_back(m.get_fact(tmp)); + } + } + if (found_false) { + pop(); + break; + } + tmp = m.get_parent(p, 0); + elim(tmp); + parents[0] = tmp; + expr* clause = m.get_fact(tmp); + if (m.is_false(clause)) { + m_refs.push_back(tmp); + result = tmp; + pop(); + break; + } + get_literals(clause); + for (unsigned i = 1; i < parents.size(); ++i) { + bool found = false; + for (unsigned j = 0; j < m_literals.size(); ++j) { + if (m.is_complement(m_literals[j], m.get_fact(parents[i].get()))) { + found = true; + break; + } + } + if (!found) { + // literal was removed as hypothesis. + parents[i] = parents.back(); + parents.pop_back(); + --i; + } + } + if (parents.size() == 1) { + result = parents[0].get(); + } + else { + result = m.mk_unit_resolution(parents.size(), parents.c_ptr()); + m_refs.push_back(result); + add_hypotheses(result); + } + pop(); + break; + } + default: { + ptr_buffer args; + bool change = false; + bool found_false = false; + for (unsigned i = 0; i < m.get_num_parents(p); ++i) { + tmp = m.get_parent(p, i); + elim(tmp); + if (m.is_false(m.get_fact(tmp))) { + result = tmp; + found_false = true; + break; + } + SASSERT(m.get_fact(tmp) == m.get_fact(m.get_parent(p, i))); + change = change || (tmp != m.get_parent(p, i)); + args.push_back(tmp); + } + if (found_false) { + break; + } + if (m.has_fact(p)) { + args.push_back(m.get_fact(p)); + } + if (change) { + tmp = m.mk_app(p->get_decl(), args.size(), args.c_ptr()); + m_refs.push_back(tmp); + } + else { + tmp = p; + } + result = tmp; + add_hypotheses(result); + break; + } + } + SASSERT(m_hypmap.contains(result)); + m_cache.insert(p, result); + p = result; + } +}; + +void proof_utils::reduce_hypotheses(proof_ref& pr) { + class reduce_hypotheses reduce(pr.get_manager()); + reduce(pr); + SASSERT(is_closed(pr.get_manager(), pr)); +} + +class proof_is_closed { + ast_manager& m; + ptr_vector m_literals; + ast_mark m_visit; + + void reset() { + m_literals.reset(); + m_visit.reset(); + } + + bool check(proof* p) { + // really just a partial check because nodes may be visited + // already under a different lemma scope. + if (m_visit.is_marked(p)) { + return true; + } + bool result = false; + m_visit.mark(p, true); + switch(p->get_decl_kind()) { + case PR_LEMMA: { + unsigned sz = m_literals.size(); + expr* cls = m.get_fact(p); + m_literals.push_back(cls); + if (m.is_or(cls)) { + m_literals.append(to_app(cls)->get_num_args(), to_app(cls)->get_args()); + } + SASSERT(m.get_num_parents(p) == 1); + result = check(m.get_parent(p, 0)); + m_literals.resize(sz); + break; + } + case PR_HYPOTHESIS: { + expr* fact = m.get_fact(p); + for (unsigned i = 0; i < m_literals.size(); ++i) { + if (m.is_complement(m_literals[i], fact)) { + result = true; + break; + } + } + break; + } + default: + result = true; + for (unsigned i = 0; i < m.get_num_parents(p); ++i) { + if (!check(m.get_parent(p, i))) { + result = false; + break; + } + } + break; + } + + return result; + } + +public: + proof_is_closed(ast_manager& m): m(m) {} + + bool operator()(proof *p) { + bool ok = check(p); + reset(); + return ok; + } +}; + +bool proof_utils::is_closed(ast_manager& m, proof* p) { + proof_is_closed checker(m); + return checker(p); +} diff --git a/lib/proof_utils.h b/lib/proof_utils.h new file mode 100644 index 000000000..10a73f388 --- /dev/null +++ b/lib/proof_utils.h @@ -0,0 +1,36 @@ +/*++ +Copyright (c) 2012 Microsoft Corporation + +Module Name: + + proof_utils.h + +Abstract: + + Utilities for proofs. + +Author: + + Nikolaj Bjorner (nbjorner) 2012-10-12. + +Revision History: + +--*/ +#ifndef _PROOF_UTILS_H_ +#define _PROOF_UTILS_H_ + +class proof_utils { +public: + /** + \brief reduce the set of hypotheses used in the proof. + */ + static void reduce_hypotheses(proof_ref& pr); + + /** + \brief Check that a proof does not contain open hypotheses. + */ + static bool is_closed(ast_manager& m, proof* p); + +}; + +#endif _PROOF_UTILS_H_