| 1 | % (c) 2019-2024 Lehrstuhl fuer Softwaretechnik und Programmiersprachen, | |
| 2 | % Heinrich Heine Universitaet Duesseldorf | |
| 3 | % This software is licenced under EPL 1.0 (http://www.eclipse.org/org/documents/epl-v10.html) | |
| 4 | ||
| 5 | :- module(cdclt_solver, [cdcl_sat_solve_predicate/2, | |
| 6 | cdcl_sat_solve_predicate_in_state/3, | |
| 7 | cdclt_solve_predicate/3, | |
| 8 | cdclt_solve_predicate/4, | |
| 9 | cdclt_solve_predicate_in_state/4, | |
| 10 | cdclt_solve_predicate_in_state/5, | |
| 11 | get_amount_of_sat_variables/2]). | |
| 12 | ||
| 13 | :- use_module(library(lists)). | |
| 14 | :- use_module(library(timeout)). | |
| 15 | :- use_module(library(sets), [subtract/3,add_element/3]). | |
| 16 | :- use_module(library(clpfd), [fd_var/1, fd_min/2, fd_max/2]). | |
| 17 | ||
| 18 | :- use_module(smt_solvers_interface(ast_optimizer_for_smt)). | |
| 19 | ||
| 20 | :- use_module(cdclt_solver('symmetry_check/smt_symmetry_breaking')). | |
| 21 | :- use_module(cdclt_solver('cdclt_stats')). | |
| 22 | :- use_module(cdclt_solver('cdclt_preprocessing')). | |
| 23 | :- use_module(cdclt_solver('cdclt_pred_to_sat')). | |
| 24 | :- use_module(cdclt_solver('cdclt_sat_solver')). | |
| 25 | :- use_module(cdclt_solver('difference_logic/difference_logic_solver')). | |
| 26 | :- use_module(cdclt_solver('difference_logic/ast_to_difference_logic')). | |
| 27 | :- use_module(cdclt_solver('cdclt_settings')). | |
| 28 | ||
| 29 | :- use_module(probsrc(debug)). | |
| 30 | :- use_module(probsrc(tools), [ajoin/2]). | |
| 31 | :- use_module(probsrc(b_ast_cleanup), [clean_up_pred/3]). | |
| 32 | :- use_module(probsrc(preferences), [get_preference/2, | |
| 33 | temporary_set_preference/2, | |
| 34 | reset_temporary_preference/1]). | |
| 35 | :- use_module(probsrc(bmachine), [b_get_machine_variables/1,b_get_machine_constants/1]). | |
| 36 | :- use_module(probsrc(b_global_sets), [list_contains_unfixed_deferred_set_id/1]). | |
| 37 | :- use_module(probsrc(b_compiler), [b_compile/6]). | |
| 38 | :- use_module(probsrc(b_enumerate), [b_tighter_enumerate_all_values/2]). | |
| 39 | :- use_module(probsrc(tools_meta), [safe_time_out/3]). | |
| 40 | :- use_module(probsrc(translate), [print_bexpr/1]). | |
| 41 | :- use_module(probsrc(tools_platform), [host_platform/1]). | |
| 42 | :- use_module(probsrc(tools), [start_ms_timer/1,stop_ms_timer_with_debug_msg/2]). | |
| 43 | :- use_module(probsrc(kernel_waitflags), [init_wait_flags/2,ground_wait_flags/1]). | |
| 44 | :- use_module(probsrc(kernel_objects), [infer_value_type/2,contains_any/1]). | |
| 45 | :- use_module(probsrc(b_interpreter), [set_up_typed_localstate/6, | |
| 46 | b_test_boolean_expression/4, | |
| 47 | b_convert_bool_timeout/7]). | |
| 48 | :- use_module(probsrc(error_manager), [error_occurred_in_error_scope/0, | |
| 49 | add_error/3, | |
| 50 | add_error_and_fail/3, | |
| 51 | add_internal_error/2, | |
| 52 | add_message/3, | |
| 53 | add_message/4, | |
| 54 | check_error_occured/2, | |
| 55 | enter_new_error_scope/2, | |
| 56 | exit_error_scope/3, | |
| 57 | clear_enumeration_warnings/0, | |
| 58 | clear_wd_errors/0, | |
| 59 | critical_enumeration_warning_occured_in_error_scope/4]). | |
| 60 | :- use_module(probsrc(bsyntaxtree), [find_typed_identifier_uses/2, | |
| 61 | find_typed_identifier_uses/3, | |
| 62 | find_identifier_uses/3, | |
| 63 | check_ast/1, | |
| 64 | get_texpr_info/2, | |
| 65 | conjunct_predicates/2, | |
| 66 | safe_create_texpr/4]). | |
| 67 | %:- use_module(probsrc('well_def/well_def_analyser'), [analyse_wd_for_expr/3]). | |
| 68 | %:- use_module(probsrc('debug'), [set_silent_mode/1]). | |
| 69 | :- use_module(probsrc(module_information), [module_info/2]). | |
| 70 | :- use_module(extension('banditfuzz/welldef')). | |
| 71 | ||
| 72 | :- module_info(group, cdclt). | |
| 73 | :- module_info(description,'This module provides a CDCL(T) based solver for B.'). | |
| 74 | ||
| 75 | :- dynamic initial_solution/0, unfixed_deferred_set_error_after_grounding/0, unfixed_deferred_set_has_been_propagated/0, solve_in_state/0, idl_candidate_constants/2, grounding_timeout_occurred/0, grounding_fallback_z3/0, additional_z3_reification/0. | |
| 76 | ||
| 77 | %% Use a small timeout for grounding with ProB and fall back to Z3 if time limit exceeded. | |
| 78 | grounding_fallback_z3. | |
| 79 | % Additionally reify constraints with Z3 using the incremental solver. | |
| 80 | %additional_z3_reification. | |
| 81 | ||
| 82 | debug_format_cdclt(_, _) :- | |
| 83 | print_logs(false), | |
| 84 | !. | |
| 85 | debug_format_cdclt(Msg, Vars) :- | |
| 86 | append(Msg, " (CDCL(T) Solver)~n", NCodes), | |
| 87 | format(NCodes, Vars), !. | |
| 88 | ||
| 89 | debug_format_cdclt(_, _, _) :- | |
| 90 | print_logs(false), | |
| 91 | !. | |
| 92 | debug_format_cdclt(Msg, Vars, Pred) :- | |
| 93 | format(Msg, Vars), | |
| 94 | translate:print_bexpr(Pred), nl, !. | |
| 95 | ||
| 96 | init :- | |
| 97 | % enumeration has to be linear for symmetry breaking but also set to false by default | |
| 98 | temporary_set_preference(randomise_enumeration_order, false), | |
| 99 | cdclt_pred_to_sat:reset_sat_var_id, | |
| 100 | cdclt_preprocessing:reset_artificial_id_counter, | |
| 101 | reset_stats, | |
| 102 | ( additional_z3_reification | |
| 103 | -> smt_solvers_interface:reset_smt_supported_interpreter | |
| 104 | ; true | |
| 105 | ), | |
| 106 | welldef:filter_typing_pos(true), | |
| 107 | welldef:preprocess_pos_for_cdclt(true), | |
| 108 | retractall(initial_solution), | |
| 109 | retractall(grounding_timeout_occurred), | |
| 110 | retractall(unfixed_deferred_set_error_after_grounding), | |
| 111 | retractall(unfixed_deferred_set_has_been_propagated), | |
| 112 | retractall(solve_in_state), | |
| 113 | retractall(idl_candidate_constants(_,_)), | |
| 114 | %temporary_set_preference(use_chr_solver, true), | |
| 115 | temporary_set_preference(unsat_core_algorithm, divide_and_conquer), % default is linear_greedy | |
| 116 | temporary_set_preference(allow_improving_wd_mode, true), | |
| 117 | %temporary_set_preference(solver_strength, 20000), | |
| 118 | temporary_set_preference(use_smt_mode, false), | |
| 119 | temporary_set_preference(use_clpfd_solver, true), | |
| 120 | temporary_set_preference(optimize_ast, true), | |
| 121 | temporary_set_preference(use_common_subexpression_elimination, false), | |
| 122 | temporary_set_preference(add_wd_pos_for_z3, false), | |
| 123 | temporary_set_preference(normalize_ast_sort_commutative, false), | |
| 124 | temporary_set_preference(normalize_ast, false). | |
| 125 | ||
| 126 | reset_preferences :- | |
| 127 | reset_temporary_preference(randomise_enumeration_order), | |
| 128 | %reset_temporary_preference(solver_strength), | |
| 129 | reset_temporary_preference(unsat_core_algorithm), | |
| 130 | reset_temporary_preference(allow_improving_wd_mode), | |
| 131 | reset_temporary_preference(use_smt_mode), | |
| 132 | reset_temporary_preference(use_clpfd_solver), | |
| 133 | reset_temporary_preference(optimize_ast), | |
| 134 | reset_temporary_preference(use_common_subexpression_elimination), | |
| 135 | reset_temporary_preference(add_wd_pos_for_z3), | |
| 136 | reset_temporary_preference(normalize_ast_sort_commutative), | |
| 137 | reset_temporary_preference(normalize_ast). | |
| 138 | ||
| 139 | set_current_state(Pred, TypedIds, State, NPred) :- | |
| 140 | reset_optimizer_state, | |
| 141 | assert_state_id_values(TypedIds, State), | |
| 142 | replace_ids_with_ground_values(Pred, 0, [], NPred), | |
| 143 | reset_optimizer_state, !. | |
| 144 | ||
| 145 | %% cdcl_sat_solve_predicate_in_state(+BoolFormula, +State, -Result). | |
| 146 | cdcl_sat_solve_predicate_in_state(BoolFormula, State, Result) :- | |
| 147 | include(ground_state_binding, State, GState), | |
| 148 | get_typed_ids_in_scope_of_pred(BoolFormula, TypedIds), | |
| 149 | set_current_state(BoolFormula, TypedIds, GState, NPred), | |
| 150 | get_equalities_from_bindings(GState, TypedIds, EqConj), | |
| 151 | safe_create_texpr(conjunct(EqConj,NPred), pred, [], Conj), | |
| 152 | cdcl_sat_solve_predicate(Conj, Result). | |
| 153 | ||
| 154 | %% cdcl_sat_solve_predicate(+BoolFormula, -Result). | |
| 155 | % Pure SAT solving without theories. | |
| 156 | cdcl_sat_solve_predicate(BoolFormula, Result) :- | |
| 157 | init, | |
| 158 | get_preference(time_out, Timeout), | |
| 159 | safe_time_out(cdcl_sat_solve_predicate_no_timeout(BoolFormula, TResult), | |
| 160 | Timeout, | |
| 161 | TimeoutResult), | |
| 162 | ( TimeoutResult == time_out | |
| 163 | -> Result = time_out | |
| 164 | ; Result = TResult | |
| 165 | ). | |
| 166 | ||
| 167 | analyze_binary_clauses(Clauses, VarsFromBinaryClauses) :- | |
| 168 | analyze_binary_clauses(Clauses, 0, 0, [], VarsFromBinaryClauses). | |
| 169 | ||
| 170 | analyze_binary_clauses([], Len, Binary, Acc, Acc) :- | |
| 171 | debug_format_cdclt("Found ~w clauses ~w of which are binary clauses~n", [Len,Binary]). | |
| 172 | analyze_binary_clauses([[_-_-_-N1,_-_-_-N2]|T], Len, Binary, Acc, VarsFromBinaryClauses) :- | |
| 173 | !, | |
| 174 | add_element(N1, Acc, NAcc1), | |
| 175 | add_element(N2, NAcc1, NAcc), | |
| 176 | Len1 is Len + 1, | |
| 177 | Binary1 is Binary + 1, | |
| 178 | analyze_binary_clauses(T, Len1, Binary1, NAcc, VarsFromBinaryClauses). | |
| 179 | analyze_binary_clauses([_|T], Len, Binary, Acc, VarsFromBinaryClauses) :- | |
| 180 | Len1 is Len + 1, | |
| 181 | analyze_binary_clauses(T, Len1, Binary, Acc, VarsFromBinaryClauses). | |
| 182 | ||
| 183 | cdcl_sat_solve_predicate_no_timeout(BoolFormula, Result) :- | |
| 184 | find_typed_identifier_uses(BoolFormula, [], SatVars), | |
| 185 | % rewriting similar to Tseitin | |
| 186 | optimize_clause_size_by_rewriting(BoolFormula, SatVars, TOptBoolFormula, NewSatVars, NewVarConjList), | |
| 187 | conjunct_predicates(NewVarConjList, NewVarConj), | |
| 188 | % add new top-level Boolean variables introduced by Tseitin rewriting | |
| 189 | conjunct_predicates([NewVarConj,TOptBoolFormula], FullBoolFormula), | |
| 190 | get_bindings_from_ids(NewSatVars, StackBindings, SatBindings, SatVarNames, _IdPrologVarTuples, PrologSatVarTriple), | |
| 191 | b_to_cnf_safe(FullBoolFormula, StackBindings, CnfBoolFormula), | |
| 192 | analyze_binary_clauses(CnfBoolFormula, VarsFromBinaryClauses), | |
| 193 | cdclt_sat_solver:create_solver_state(SatStateMutable), | |
| 194 | SolverEnv = [sat_vars:NewSatVars,sat_varnames:SatVarNames,sat_bindings:SatBindings,t2b_env:_,pl_vars:PrologSatVarTriple,binary_clause_vars:VarsFromBinaryClauses], | |
| 195 | ( ( cdclt_sat_solver:unit_propagate_cnf(SolverEnv, SatStateMutable, CnfBoolFormula, TCnfBoolFormula2), | |
| 196 | remove_subsumed_clauses(TCnfBoolFormula2, TCnfBoolFormula3), | |
| 197 | sort(TCnfBoolFormula3, CnfBoolFormula2), | |
| 198 | cdclt_sat_solver:sat(default_sat, SolverEnv, SatStateMutable, CnfBoolFormula2) | |
| 199 | ) | |
| 200 | -> Result = solution(SatBindings) | |
| 201 | ; Result = contradiction_found | |
| 202 | ), | |
| 203 | print_stats. | |
| 204 | ||
| 205 | %% cdclt_solve_predicate(+Pred, -SolvedPred, -Result). | |
| 206 | % Main interface predicate for CDCL(T) based solver which times out | |
| 207 | % with respect to ProB's time_out preference. | |
| 208 | cdclt_solve_predicate(Pred, SolvedPred, Result) :- | |
| 209 | cdclt_solve_predicate(default, Pred, SolvedPred, Result). | |
| 210 | ||
| 211 | %% cdclt_solve_predicate(+SolverName, +Pred, -SolvedPred, -Result). | |
| 212 | cdclt_solve_predicate(SolverName, Pred, SolvedPred, Result) :- | |
| 213 | init, | |
| 214 | cdclt_solve_predicate_timeout(SolverName, Pred, SolvedPred, Result). | |
| 215 | ||
| 216 | %% cdclt_solve_predicate_in_state(+Pred, +Pred, -Result). | |
| 217 | cdclt_solve_predicate_in_state(Pred, State, SolvedPred, Result) :- | |
| 218 | cdclt_solve_predicate_in_state(default, Pred, State, SolvedPred, Result). | |
| 219 | ||
| 220 | %% cdclt_solve_predicate_in_state(+SolverName, +Pred, +Pred, -Result). | |
| 221 | cdclt_solve_predicate_in_state(SolverName, Pred, State, SolvedPred, Result) :- | |
| 222 | init, | |
| 223 | include(ground_state_binding, State, GState), | |
| 224 | get_typed_ids_in_scope_of_pred(Pred, TypedIds), | |
| 225 | set_current_state(Pred, TypedIds, GState, NPred), | |
| 226 | get_equalities_from_bindings(GState, TypedIds, EqConj), | |
| 227 | safe_create_texpr(conjunct(EqConj,NPred), pred, [], Conj), | |
| 228 | asserta(solve_in_state), | |
| 229 | cdclt_solve_predicate_timeout(SolverName, Conj, SolvedPred, Result). | |
| 230 | ||
| 231 | cdclt_solve_predicate_timeout(SolverName, Pred, SolvedPred, Result) :- | |
| 232 | get_preference(time_out, Timeout), | |
| 233 | safe_time_out(cdclt_solve_predicate_initialized(SolverName, Pred, TSolvedPred, SolverResult), | |
| 234 | Timeout, | |
| 235 | TimeoutResult), | |
| 236 | ( TimeoutResult == time_out | |
| 237 | -> Result = time_out, | |
| 238 | SolvedPred = Pred | |
| 239 | ; ground(SolverResult), | |
| 240 | Result = SolverResult, | |
| 241 | SolvedPred = TSolvedPred | |
| 242 | ). | |
| 243 | ||
| 244 | cdclt_solve_predicate_initialized(SolverName, Pred, SolvedPred, Result) :- | |
| 245 | ( check_ast(Pred) | |
| 246 | -> ( (error_occurred_in_error_scope, check_error_occured(check_ast_typing, _)) | |
| 247 | -> % input AST is not correctly typed | |
| 248 | SolvedPred = Pred, | |
| 249 | Result = error | |
| 250 | ; cdclt_solve_predicate_initialized_checked(SolverName, Pred, SolvedPred, Result) | |
| 251 | ) | |
| 252 | ; add_message(cdclt_solve_predicate_no_timeout, 'Input predicate is not well-defined~nProB\'s SMT solver transforms the input to a well-defined predicate before constraint solving~nFor: ',SolverName), | |
| 253 | cdclt_solve_predicate_initialized_checked(SolverName, Pred, SolvedPred, Result) | |
| 254 | ). | |
| 255 | ||
| 256 | cdclt_solve_predicate_initialized_checked(SolverName, Pred, SolvedPred, Result) :- | |
| 257 | call_cleanup(catch(cdclt_solve_predicate_no_timeout(SolverName, Pred, SolvedPred, Result), | |
| 258 | Exception, % TO DO: catch specific exception instead of all | |
| 259 | handle_clpfd_overflow(SolverName, Exception, Pred, SolvedPred, Result) | |
| 260 | ), | |
| 261 | reset_preferences | |
| 262 | ). | |
| 263 | ||
| 264 | %% handle_clpfd_overflow(+SolverName, +Exception, +Pred, -SolvedPred, -Result). | |
| 265 | % Disable CLP(FD) interface and restart CDCL(T) if an overflow occurred. | |
| 266 | handle_clpfd_overflow(SolverName, Exception, Pred, SolvedPred, Result) :- | |
| 267 | Exception = error(representation_error(Err),_), | |
| 268 | memberchk(Err, ['CLPFD integer overflow','max_clpfd_integer','min_clpfd_integer']), | |
| 269 | !, | |
| 270 | debug_format_cdclt("Handle CLP(FD) overflow error and restart CDCL(T) without CLP(FD).", []), | |
| 271 | temporary_set_preference(use_clpfd_solver, false), | |
| 272 | cdclt_solve_predicate_initialized(SolverName, Pred, SolvedPred, Result). | |
| 273 | handle_clpfd_overflow(SolverName, Exception, Pred, SolvedPred, unknown) :- | |
| 274 | ajoin(['Exception occured during solving with ',SolverName,':'],Msg), | |
| 275 | add_error(cdclt,Msg,Exception), | |
| 276 | SolvedPred = Pred. | |
| 277 | ||
| 278 | %% cdclt_solve_predicate(+SolverName, +Pred, -SolvedPred, -SolverResult). | |
| 279 | % Interface predicate for CDCL(T) based solver without a timeout. | |
| 280 | cdclt_solve_predicate_no_timeout(SolverName, Pred, SolvedPred, SolverResult) :- | |
| 281 | debug_format_cdclt("Parsed predicate: ", [], Pred), | |
| 282 | debug_format_cdclt("Start preprocessing", []), | |
| 283 | cdclt_solve_predicate_no_timeout_cleanup(SolverName, Pred, NPred), | |
| 284 | debug_format_cdclt("After preprocessing: ", [], NPred), | |
| 285 | cdclt_solve_predicate_no_timeout_clean(SolverName, NPred, SolvedPred, SolverResult). | |
| 286 | ||
| 287 | cdclt_solve_predicate_no_timeout_clean(_, b(truth,pred,Info), SolvedPred, SolverResult) :- | |
| 288 | !, | |
| 289 | SolvedPred = b(truth,pred,Info), | |
| 290 | SolverResult = solution([]). | |
| 291 | cdclt_solve_predicate_no_timeout_clean(_, b(falsity,pred,Info), SolvedPred, SolverResult) :- | |
| 292 | !, | |
| 293 | SolvedPred = b(falsity,pred,Info), | |
| 294 | SolverResult = contradiction_found. | |
| 295 | cdclt_solve_predicate_no_timeout_clean(SolverName, CleanPred, SolvedPred, SolverResult) :- | |
| 296 | cdclt_solve_predicate_no_timeout_setup(SolverName, CleanPred, Env, SmtVars, SatVars, WDPosBoolFormula, AnalysisBoolFormula, SmtBoolFormula), | |
| 297 | cdclt_sat_solver:cnf_to_smt(SmtBoolFormula, SolvedPred), | |
| 298 | start_ms_timer(Timer), | |
| 299 | enter_new_error_scope(ScopeID, cdclt_solve_predicate), | |
| 300 | call_cleanup(conflict_driven_clause_learning_online(SolverName, Env, SmtVars, SatVars, WDPosBoolFormula, AnalysisBoolFormula, SmtBoolFormula, SolverResult), | |
| 301 | exit_error_scope(ScopeID, _, cdclt_solve_predicate)), | |
| 302 | stop_ms_timer_with_debug_msg(Timer, cdclt_solving_success), | |
| 303 | print_stats. | |
| 304 | ||
| 305 | %% cdclt_solve_predicate_no_timeout_cleanup(+SolverName, +Pred, -NPred). | |
| 306 | % Some rewriting and cleanup for CDCL(T). | |
| 307 | cdclt_solve_predicate_no_timeout_cleanup(SolverName, Pred, NPred) :- | |
| 308 | start_ms_timer(Timer0), | |
| 309 | simplify_negation(Pred, SPred), | |
| 310 | debug_format_cdclt("After simplifying negations: ", [], SPred), | |
| 311 | reset_optimizer_state, | |
| 312 | assert_ground_id_values(0, SPred), | |
| 313 | replace_ids_with_ground_values(SPred, 0, [], AstOpt), | |
| 314 | precompute_values_non_recursive([instantiate_quantifier_limit(0),instantiate_sets_limit(1000)], AstOpt, AstPrecomputed), | |
| 315 | debug_format_cdclt("After precomputing values: ", [], AstOpt), | |
| 316 | stop_ms_timer_with_debug_msg(Timer0, precompute_values_for_cdclt), | |
| 317 | find_identifier_uses(AstPrecomputed, [], UsedIds), | |
| 318 | temporary_set_preference(data_validation_mode, true), | |
| 319 | catch(b_compile(AstPrecomputed, UsedIds, [], [], AstCompiled, no_wf_available), | |
| 320 | enumeration_warning(_,_,_,_,_), % cancel if enumeration warning has occurred | |
| 321 | AstCompiled = AstPrecomputed), | |
| 322 | clear_wd_errors, % b_compile might throw a wd error | |
| 323 | reset_temporary_preference(data_validation_mode), | |
| 324 | debug_format_cdclt("After compiling values (used ids ~w): ", [UsedIds], AstCompiled), | |
| 325 | start_ms_timer(Timer2), | |
| 326 | ( clean_up_pred(AstCompiled, _, NPred) | |
| 327 | -> true | |
| 328 | ; add_internal_error('Clean up failed ',SolverName), | |
| 329 | NPred = AstCompiled | |
| 330 | ), | |
| 331 | stop_ms_timer_with_debug_msg(Timer2, clean_up_pred_for_cdclt). | |
| 332 | ||
| 333 | %% cdclt_solve_predicate_no_timeout_setup(+SolverName, +CleanPred, -Env, -SmtVars, -SatVars, -WDPosBoolFormula, -AnalysisBoolFormula, -SmtBoolFormula). | |
| 334 | cdclt_solve_predicate_no_timeout_setup(SolverName, CleanPred, Env, SmtVars, SatVars, WDPosBoolFormula, AnalysisBoolFormula, SmtBoolFormula) :- | |
| 335 | start_ms_timer(Timer), | |
| 336 | cdclt_settings:static_syntax_analysis(PerformStaticAnalysis), | |
| 337 | is_rewrite_to_idl(RewriteToIdl), | |
| 338 | add_symmetry_breaking_predicates_cond(CleanPred, CleanPredSym), | |
| 339 | preprocess_predicate(PerformStaticAnalysis, RewriteToIdl, CleanPredSym, LiftedPred, InferredImplsConj, CandidateImpls), | |
| 340 | log_inferred_implications(InferredImplsConj), | |
| 341 | debug_format_cdclt("Lifted predicate: ", [], LiftedPred), | |
| 342 | predicate_to_sat(normal_make_wd, LiftedPred, Env1, WDPOs1, SmtBoolFormula, SatVars1), !, | |
| 343 | ( PerformStaticAnalysis == true | |
| 344 | -> predicate_to_sat(only_reuse, WDPOs1, SatVars1, Env1, InferredImplsConj, Env3, WDPOsList, ImplBoolFormula, SatVars3), | |
| 345 | debug_format_cdclt("Static analysis inferred: ", [], ImplBoolFormula) | |
| 346 | ; ImplBoolFormula = b(truth,pred,[]), | |
| 347 | SatVars3 = SatVars1, Env3 = Env1, WDPOsList = WDPOs1 | |
| 348 | ), | |
| 349 | % well-definedness implication on the top-level to encode well-definedness in the SAT solver | |
| 350 | conjunct_predicates(WDPOsList, WDPOs), | |
| 351 | preprocess_predicate(false, RewriteToIdl, WDPOs, LiftedWDPOs, _, WDCandidateImpls), | |
| 352 | % theory deduction for wd | |
| 353 | get_wd_theory_implications(CandidateImpls, WDCandidateImpls, WDTheoryImpls), | |
| 354 | log_inferred_wd_theory_implications(InferredImplsConj), | |
| 355 | predicate_to_sat(normal_make_wd, [], SatVars3, Env3, WDTheoryImpls, Env4, _, WDTheoryImplsBoolFormula, SatVars4), | |
| 356 | debug_format_cdclt("WD Theory Deduction Implications: ", [], WDTheoryImpls), | |
| 357 | predicate_to_sat(normal_make_wd, [], SatVars4, Env4, LiftedWDPOs, Env5, WDPOsList2, WDPosBoolFormula1, SatVars5), | |
| 358 | % there can be WD POs that also have WD POs | |
| 359 | conjunct_predicates(WDPOsList2, WDPOs2), | |
| 360 | preprocess_predicate(false, RewriteToIdl, WDPOs2, LiftedWDPOs2, _, _), | |
| 361 | predicate_to_sat(normal_make_wd, [], SatVars5, Env5, LiftedWDPOs2, Env, _, WDPosBoolFormula2, SatVars), | |
| 362 | conjunct_predicates([WDPosBoolFormula1,WDPosBoolFormula2], WDPosBoolFormula), | |
| 363 | ( SatVars3 \== SatVars | |
| 364 | -> %set_silent_mode(on), | |
| 365 | %analyse_wd_for_expr(CleanPred, _ResStr, _IsWd), % prints missing WD conditions | |
| 366 | %set_silent_mode(off), | |
| 367 | add_message(cdclt_solve_predicate_no_timeout_setup, 'Input not well-defined: automatically added WD POs for solver ',SolverName) | |
| 368 | ; true | |
| 369 | ), | |
| 370 | debug_format_cdclt("WD Implication: ", [], LiftedWDPOs), | |
| 371 | debug_format_cdclt("WD Implication as Boolean formula: ", [], WDPosBoolFormula), | |
| 372 | !, | |
| 373 | conjunct_predicates([WDTheoryImplsBoolFormula,ImplBoolFormula], AnalysisBoolFormula), | |
| 374 | find_typed_identifier_uses(LiftedPred, SmtVars), % no global sets or constants are considered here | |
| 375 | debug_format_cdclt("Check LiftedPred", []), | |
| 376 | ( check_ast(LiftedPred) | |
| 377 | -> true | |
| 378 | ; add_internal_error('AST is missing well-definedness information for solver ', SolverName) | |
| 379 | ), | |
| 380 | %translate:nested_print_bexpr(LiftedPred),nl, | |
| 381 | stop_ms_timer_with_debug_msg(Timer, preprocessing_for_cdclt), | |
| 382 | debug_format_cdclt("End preprocessing", []), | |
| 383 | !. | |
| 384 | cdclt_solve_predicate_no_timeout_setup(_, Pred, _, _, _, _, _, _) :- | |
| 385 | add_error_and_fail(cdclt_solve_predicate, 'Cannot create SAT formula from B predicate:', [Pred]). | |
| 386 | ||
| 387 | %% add_symmetry_breaking_predicates_cond(+Pred, -SymPred). | |
| 388 | add_symmetry_breaking_predicates_cond(Pred, SymPred) :- | |
| 389 | cdclt_settings:static_smt_symmetry_breaking(true), | |
| 390 | !, | |
| 391 | debug_format_cdclt("Start symmetry breaking..", []), | |
| 392 | ( add_symmetry_breaking_predicates(Pred, SymPred), | |
| 393 | get_amount_of_found_sbps(FoundSBPs), | |
| 394 | log_symmetry_breaking_stats(FoundSBPs), | |
| 395 | debug_format_cdclt("done.", []) | |
| 396 | -> true | |
| 397 | ; add_message(smt_symmetry_breaking, 'Symmetry breaking failed for: ',Pred), | |
| 398 | SymPred = Pred | |
| 399 | ). | |
| 400 | add_symmetry_breaking_predicates_cond(Pred, Pred). | |
| 401 | ||
| 402 | is_rewrite_to_idl(Res) :- | |
| 403 | get_preference(cdclt_use_idl_theory_solver,true), | |
| 404 | !, | |
| 405 | Res = true. | |
| 406 | is_rewrite_to_idl(false). | |
| 407 | ||
| 408 | %% conflict_driven_clause_learning_online(+SolverName, +Env, +SmtVars, +SatVars, +WDPosBoolFormula, +AnalysisBoolFormula, +BoolFormula, -SolverResult). | |
| 409 | % Conflict-driven clause learning from incomplete assignments of the boolean formula by setting up reification constraints (b_interpreter:b_convert_bool_timeout/7) | |
| 410 | % connecting the SAT and theory solver. IDL constraint solver uses custom coroutines for the reification (cdclt_solver:idl_solver_interface/3). | |
| 411 | conflict_driven_clause_learning_online(_, _, _, _, _, _, b(truth,pred,_), SolverResult) :- | |
| 412 | !, | |
| 413 | SolverResult = solution([]). | |
| 414 | conflict_driven_clause_learning_online(_, _, _, _, _, _, b(falsity,pred,_), SolverResult) :- | |
| 415 | !, | |
| 416 | SolverResult = contradiction_found. | |
| 417 | conflict_driven_clause_learning_online(SolverName, Env, SmtVars, SatVars, WDPosBoolFormula, AnalysisBoolFormula, BoolFormula, SolverResult) :- | |
| 418 | % rewriting similar to Tseitin | |
| 419 | optimize_clause_size_by_rewriting(BoolFormula, SatVars, TOptBoolFormula, NewSatVars, NewVarConjList), | |
| 420 | conjunct_predicates(NewVarConjList, NewVarConj), | |
| 421 | conjunct_predicates([AnalysisBoolFormula,TOptBoolFormula], TFullBoolFormula), | |
| 422 | % add new top-level Boolean variables introduced by Tseitin rewriting | |
| 423 | ( NewVarConj = b(truth,pred,_) | |
| 424 | -> FullBoolFormula = TFullBoolFormula | |
| 425 | ; safe_create_texpr(conjunct(NewVarConj,TFullBoolFormula), pred, [], FullBoolFormula) | |
| 426 | ), | |
| 427 | get_bindings_from_ids(NewSatVars, StackBindings, SatBindings, SatVarNames, IdPrologVarTuples, PlVars), | |
| 428 | b_to_cnf_safe(FullBoolFormula, StackBindings, CnfBoolFormula), | |
| 429 | setup_theory_wf_store(SmtVars, SmtBindings, WfStoreSmt), | |
| 430 | ( get_preference(cdclt_use_idl_theory_solver,true) | |
| 431 | -> difference_logic_solver:init_idl_solver(IdlGraphMut) | |
| 432 | ; true | |
| 433 | ), | |
| 434 | ( WDPosBoolFormula \= b(truth,pred,_) | |
| 435 | -> % add WD PO implications | |
| 436 | b_to_cnf_safe(WDPosBoolFormula, StackBindings, CnfWDPosBoolFormula), | |
| 437 | append(CnfWDPosBoolFormula, CnfBoolFormula, TNCnfBoolFormula) | |
| 438 | ; TNCnfBoolFormula = CnfBoolFormula | |
| 439 | ), | |
| 440 | cdclt_sat_solver:create_solver_state(SatStateMutable), | |
| 441 | SolverEnv = [sat_vars:NewSatVars,sat_varnames:SatVarNames,sat_bindings:SatBindings,t2b_env:Env,pl_vars:PlVars], | |
| 442 | log_unique_predicates(SatVarNames), | |
| 443 | cdclt_sat_solver:unit_propagate_cnf(SolverEnv, SatStateMutable, TNCnfBoolFormula, TNCnfBoolFormula2), | |
| 444 | remove_subsumed_clauses(TNCnfBoolFormula2, TNCnfBoolFormula3), | |
| 445 | sort(TNCnfBoolFormula3, NCnfBoolFormula), | |
| 446 | ( debug_mode(off) | |
| 447 | -> true | |
| 448 | ; cdclt_sat_solver:portray_cnf(NCnfBoolFormula) | |
| 449 | ), | |
| 450 | analyze_binary_clauses(CnfBoolFormula, VarsFromBinaryClauses), | |
| 451 | NSolverEnv = [sat_vars:NewSatVars,sat_varnames:SatVarNames,sat_bindings:SatBindings,t2b_env:Env,pl_vars:PlVars,binary_clause_vars:VarsFromBinaryClauses], | |
| 452 | setup_reification(SatStateMutable, IdlGraphMut, IdPrologVarTuples, SmtBindings, WfStoreSmt), | |
| 453 | cdclt_sat_solver:sat(SolverName, NSolverEnv, SatStateMutable, NCnfBoolFormula, [allow_partial_model]), | |
| 454 | retractall(unfixed_deferred_set_error_after_grounding), | |
| 455 | cdclt_sat_solver:clear_pending_theory_propagations, | |
| 456 | ( severe_error_occurred | |
| 457 | -> % error occurred during SAT solving | |
| 458 | !, | |
| 459 | SolverResult = error | |
| 460 | ; debug_format_cdclt("Ground waitflags", []), | |
| 461 | remove_pending_theory_propagations_on_bt, | |
| 462 | ( conflict_driven_clause_learning_online_grounding(IdlGraphMut, SmtVars, PlVars, SmtBindings, WfStoreSmt, SolverResult) | |
| 463 | -> true | |
| 464 | ; ( unfixed_deferred_set_has_been_propagated | |
| 465 | -> % a contradiction has been found but at least one unfixed deferred set has been propagated | |
| 466 | % we thus do not know if it's a genuine contradiction | |
| 467 | asserta(unfixed_deferred_set_error_after_grounding), | |
| 468 | fail | |
| 469 | ; % this is a genuine contradiction after grounding | |
| 470 | fail | |
| 471 | ) | |
| 472 | ) | |
| 473 | ). | |
| 474 | conflict_driven_clause_learning_online(_, _, _, _, _, _, _, Res) :- | |
| 475 | \+ wd_error_occurred_in_error_scope(_), | |
| 476 | critical_enumeration_warning_occured_in_error_scope(A, B, C, D), | |
| 477 | !, | |
| 478 | clear_enumeration_warnings, | |
| 479 | Res = no_solution_found(enumeration_warning(A,B,C,D,critical)). | |
| 480 | conflict_driven_clause_learning_online(_, _, _, _, _, _, _, Res) :- | |
| 481 | severe_error_occurred, | |
| 482 | !, | |
| 483 | Res = error. | |
| 484 | conflict_driven_clause_learning_online(_, _, _, _, _, _, _, error) :- | |
| 485 | \+ initial_solution, | |
| 486 | wd_error_occurred_in_error_scope(Reason), | |
| 487 | add_error(well_definedness_error, Reason, 'conflict_driven_clause_learning_online'). | |
| 488 | conflict_driven_clause_learning_online(_, _, _, _, _, _, _, Res) :- | |
| 489 | ( grounding_timeout_occurred | |
| 490 | ; cdclt_sat_solver:assignment_timeout_occurred | |
| 491 | ), | |
| 492 | !, | |
| 493 | Res = no_solution_found(solver_answered_unknown). | |
| 494 | conflict_driven_clause_learning_online(_, _, _, _, _, _, _, no_solution_found(unfixed_deferred_sets)) :- | |
| 495 | \+ solve_in_state, | |
| 496 | ( cdclt_sat_solver:unfixed_deferred_set_error_occurred | |
| 497 | ; unfixed_deferred_set_error_after_grounding | |
| 498 | ; unfixed_deferred_set_has_been_propagated | |
| 499 | ). | |
| 500 | conflict_driven_clause_learning_online(_, _, _, _, _, _, _, contradiction_found) :- | |
| 501 | \+ initial_solution, | |
| 502 | ( solve_in_state | |
| 503 | -> true | |
| 504 | ; \+ unfixed_deferred_set_error_after_grounding, | |
| 505 | \+ cdclt_sat_solver:unfixed_deferred_set_error_occurred, | |
| 506 | \+ unfixed_deferred_set_has_been_propagated | |
| 507 | ), | |
| 508 | \+ wd_error_occurred_in_error_scope(_). | |
| 509 | ||
| 510 | % Backtracking after grounding: We have either found a solution or a conflict. | |
| 511 | % If a conflict has been found, it can be a theory conflict in the SAT solver | |
| 512 | % (pending props have already been used for CDCL) or we just backtrack to the last decision. | |
| 513 | % In all cases, remove all pending theory propagations now. | |
| 514 | remove_pending_theory_propagations_on_bt :- | |
| 515 | ( true | |
| 516 | ; cdclt_sat_solver:clear_pending_theory_propagations, | |
| 517 | debug_format_cdclt("backtrack after grounding waitflags", []), | |
| 518 | fail | |
| 519 | ). | |
| 520 | ||
| 521 | :- use_module(smt_solvers_interface(smt_solvers_interface)). | |
| 522 | ||
| 523 | %% conflict_driven_clause_learning_online_grounding(+IdlGraphMut, +SmtVars, +PlVars, +SmtBindings, +WfStoreSmt, -SolverResult). | |
| 524 | % Use IDL theory solver if option is set and enter ProB's grounding phase. | |
| 525 | conflict_driven_clause_learning_online_grounding(IdlGraphMut, SmtVars, _, SmtBindings, WfStoreSmt, SolverResult) :- | |
| 526 | get_preference(cdclt_use_idl_theory_solver, true), | |
| 527 | !, | |
| 528 | add_typed_ids_in_scope(SmtVars, StateTypedIds), | |
| 529 | propagate_fd_bounds_to_idl_solver(IdlGraphMut, SmtBindings), | |
| 530 | get_idl_solution_bindings(IdlGraphMut, IdlBindings), | |
| 531 | ( % try one idl solution first since ProB could timeout, but do not enumerate idl solutions here | |
| 532 | set_bindings(IdlBindings, SmtBindings), | |
| 533 | ground_wait_flags(WfStoreSmt) | |
| 534 | ; exclude_idl_solution(IdlBindings, SmtBindings, StateTypedIds, WfStoreSmt), | |
| 535 | ( ground_wait_flags(WfStoreSmt) | |
| 536 | ; % fallback to idl solver and possibly enumerate all solutions | |
| 537 | critical_enumeration_warning_occured_in_error_scope(_, _, _, _), | |
| 538 | clear_enumeration_warnings, | |
| 539 | propagate_idl_solution_to_bindings_bt(IdlGraphMut, SmtBindings), | |
| 540 | ground_wait_flags(WfStoreSmt) | |
| 541 | ) | |
| 542 | ), | |
| 543 | finalize_smt_solution(SmtBindings, SolverResult). | |
| 544 | conflict_driven_clause_learning_online_grounding(_, _, PlVars, SmtBindings, WfStoreSmt, SolverResult) :- | |
| 545 | grounding_fallback_z3, | |
| 546 | \+ host_platform(windows), | |
| 547 | !, | |
| 548 | ( grounding_timeout_occurred | |
| 549 | -> % use a smaller time limit if a timeout already occurred | |
| 550 | GroundTimeout = 5000 | |
| 551 | ; GroundTimeout = 30000 | |
| 552 | ), | |
| 553 | safe_time_out(ground_wait_flags(WfStoreSmt), | |
| 554 | GroundTimeout, | |
| 555 | TimeOutRes), | |
| 556 | ( TimeOutRes == time_out | |
| 557 | -> debug_format_cdclt("Timeout when grounding with ProB. Fallback to Z3.", []), | |
| 558 | sat_bindings_to_smt_formula(PlVars, SmtFormula), | |
| 559 | safe_time_out(smt_solve_predicate(z3, [check_sat_skeleton(0)], SmtFormula, _, Result), | |
| 560 | 10000, | |
| 561 | TimeOutResZ3), | |
| 562 | ( ( TimeOutResZ3 == time_out | |
| 563 | ; Result = no_solution_found(_) | |
| 564 | ) | |
| 565 | -> asserta(grounding_timeout_occurred), | |
| 566 | fail | |
| 567 | ; Result = solution(NSmtBindings), | |
| 568 | debug_format_cdclt("Validating Z3 solution with ProB.", []), | |
| 569 | copy_smt_solution(SmtBindings,NSmtBindings), | |
| 570 | ground_wait_flags(WfStoreSmt) % also avoid Call Residue warnings | |
| 571 | ) | |
| 572 | ; NSmtBindings = SmtBindings | |
| 573 | ), | |
| 574 | \+ severe_error_occurred, | |
| 575 | finalize_smt_solution(NSmtBindings, SolverResult). | |
| 576 | conflict_driven_clause_learning_online_grounding(_, _, _, SmtBindings, WfStoreSmt, SolverResult) :- | |
| 577 | ground_wait_flags(WfStoreSmt), | |
| 578 | \+ severe_error_occurred, | |
| 579 | finalize_smt_solution(SmtBindings, SolverResult). | |
| 580 | ||
| 581 | :- use_module(probsrc(kernel_objects),[equal_object/3]). | |
| 582 | copy_smt_solution([],_). | |
| 583 | copy_smt_solution([bind(ID,Val)|T],NSmtBindings) :- member(binding(ID,Z3Val,_),NSmtBindings),!, | |
| 584 | equal_object(Val,Z3Val,copy_smt_solution), | |
| 585 | copy_smt_solution(T,NSmtBindings). | |
| 586 | copy_smt_solution([H|T],NSmtBindings) :- add_internal_error('No SMT solution for:',H), | |
| 587 | copy_smt_solution(T,NSmtBindings). | |
| 588 | ||
| 589 | finalize_smt_solution(SmtBindings, SolverResult) :- | |
| 590 | SolverResult = solution(SmtBindings), | |
| 591 | log_solution, | |
| 592 | announce_bt_from_smt_solution. | |
| 593 | ||
| 594 | wd_error_occurred_in_error_scope(Reason) :- | |
| 595 | error_occurred_in_error_scope, | |
| 596 | check_error_occured(well_definedness_error, Reason). | |
| 597 | ||
| 598 | %% setup_reification(+SatStateMutable, +IdlGraphMut, +IdPrologVarTuples, +SmtBindings, +WfStoreSmt). | |
| 599 | % Set up reification constraints to connect SAT and SMT solver. | |
| 600 | % Propagates in both directions (b_convert_bool_timeout/7 implements theory propagation). | |
| 601 | setup_reification(_, _, [], _, _). | |
| 602 | setup_reification(SatStateMutable, IdlGraphMut, [(SatId,SatPrologVar,StackInfo)|T], SmtBindings, WfStoreSmt) :- | |
| 603 | get_texpr_info(SatId, Info), | |
| 604 | memberchk(smt_formula(SmtFormula), Info), | |
| 605 | !, | |
| 606 | SatId = b(identifier(SatVarName),_,_), | |
| 607 | setup_reification_for_solver(SatStateMutable, StackInfo, IdlGraphMut, SmtFormula, SmtBindings, WfStoreSmt, SatVarName, SatPrologVar), | |
| 608 | %sat_smt_reification(SatId, SatPrologVar, SmtBindings, WfStoreSmt) | |
| 609 | setup_reification(SatStateMutable, IdlGraphMut, T, SmtBindings, WfStoreSmt). | |
| 610 | setup_reification(SatStateMutable, IdlGraphMut, [_|T], SmtBindings, WfStoreSmt) :- | |
| 611 | % plain SAT variable which is not reified with a theory solver, e.g., introduced by CNF optimization rewriting | |
| 612 | setup_reification(SatStateMutable, IdlGraphMut, T, SmtBindings, WfStoreSmt). | |
| 613 | ||
| 614 | /*:- block sat_smt_reification(?, -, ?, ?). | |
| 615 | % without theory propagation | |
| 616 | sat_smt_reification(SatId, SatPrologVar, SmtBindings, WfStoreSmt) :- | |
| 617 | get_texpr_info(SatId, Info), | |
| 618 | memberchk(smt_formula(TSmtFormula), Info), | |
| 619 | ( SatPrologVar == pred_true | |
| 620 | -> SmtFormula = TSmtFormula | |
| 621 | ; SmtFormula = b(negation(TSmtFormula),pred,[]) | |
| 622 | ), | |
| 623 | b_test_boolean_expression(SmtFormula, SmtBindings, [], WfStoreSmt).*/ | |
| 624 | ||
| 625 | setup_reification_for_solver(SatStateMutable, StackInfo, IdlGraphMut, SmtFormula, SmtBindings, WfStoreSmt, SatVarName, SatPrologVar) :- | |
| 626 | find_typed_identifier_uses(SmtFormula, [], UsedIds), | |
| 627 | ( list_contains_unfixed_deferred_set_id(UsedIds) | |
| 628 | -> ContainsUnfixed = true, | |
| 629 | ( ground(SatPrologVar) | |
| 630 | -> DetUnfixedProp = true | |
| 631 | ; DetUnfixedProp = false | |
| 632 | ) | |
| 633 | ; ContainsUnfixed = false, | |
| 634 | DetUnfixedProp = false | |
| 635 | ), | |
| 636 | log_theory_propagation_and_unfixed_deferred_sets(SatStateMutable, SatVarName, StackInfo, ContainsUnfixed, DetUnfixedProp, SatPrologVar), | |
| 637 | sat_var_assignment_timeout(ATO), | |
| 638 | ( get_preference(cdclt_use_idl_theory_solver,true) | |
| 639 | -> ( ast_to_difference_logic:rewrite_to_idl(SmtFormula, [DLConstraint]) | |
| 640 | -> debug_format_cdclt("Constraint for IDL: ", [], SmtFormula), | |
| 641 | idl_solver_interface(IdlGraphMut, SatPrologVar, [DLConstraint]) | |
| 642 | ; log_idl_candidates_from_constraint(SmtFormula) | |
| 643 | -> propagate_non_idl_to_idl(IdlGraphMut, SatPrologVar, SmtFormula) | |
| 644 | ; true | |
| 645 | ), | |
| 646 | sat_debug_msg(SatVarName,SmtFormula,SatPrologVar), | |
| 647 | b_convert_bool_timeout(SmtFormula, SmtBindings, [], WfStoreSmt, SatPrologVar, ATO, TORes), | |
| 648 | cdclt_sat_solver:log_det_theory_timeout(TORes) | |
| 649 | ; sat_debug_msg(SatVarName,SmtFormula,SatPrologVar), | |
| 650 | b_convert_bool_timeout(SmtFormula, SmtBindings, [], WfStoreSmt, SatPrologVar, ATO, TORes), | |
| 651 | cdclt_sat_solver:log_det_theory_timeout(TORes), | |
| 652 | ( additional_z3_reification | |
| 653 | -> kernel_waitflags:get_wait_flag1(smt_call, WfStoreSmt, BeforeEnumWF), | |
| 654 | gensym:gensym(smt_assertion_name, Symbol), | |
| 655 | z3_coroutine(SatPrologVar, BeforeEnumWF, Symbol, SmtBindings, SmtFormula) | |
| 656 | ; true | |
| 657 | ) | |
| 658 | ). | |
| 659 | ||
| 660 | :- block z3_coroutine(-, ?, ?, ?, ?). | |
| 661 | z3_coroutine(SatPrologVar, BeforeEnumWF, Symbol, SmtBindings, SmtFormula) :- | |
| 662 | SatPrologVar == pred_true, | |
| 663 | smt_add_predicate(BeforeEnumWF, SmtFormula, SmtBindings, SmtBindings, Symbol). | |
| 664 | z3_coroutine(SatPrologVar, BeforeEnumWF, Symbol, SmtBindings, SmtFormula) :- | |
| 665 | SatPrologVar == pred_false, | |
| 666 | safe_create_texpr(negation(SmtFormula), pred, [], NSmtFormula), | |
| 667 | smt_add_predicate(BeforeEnumWF, NSmtFormula, SmtBindings, SmtBindings, Symbol). | |
| 668 | ||
| 669 | sat_debug_msg(SatVarName, SmtFormula, SatPrologVar) :- | |
| 670 | ( debug_mode(off) | |
| 671 | -> true | |
| 672 | ; format('% SAT variable ~w (~w) for ProB pred: ',[SatVarName,SatPrologVar]), | |
| 673 | print_bexpr(SmtFormula), nl | |
| 674 | ). | |
| 675 | ||
| 676 | :- block log_theory_propagation_and_unfixed_deferred_sets(?, ?, ?, ?, ?, -). | |
| 677 | log_theory_propagation_and_unfixed_deferred_sets(SatStateMutable, SatVarName, StackInfo, ContainsUnfixed, DetUnfixedProp, SatPrologVar) :- | |
| 678 | ground(SatPrologVar), | |
| 679 | ( ContainsUnfixed | |
| 680 | -> % log unfixed deferred set propagation but undo when backtracking if not deterministic unit propagation | |
| 681 | asserta(unfixed_deferred_set_has_been_propagated), | |
| 682 | ( DetUnfixedProp | |
| 683 | -> true | |
| 684 | ; (true; (retract(unfixed_deferred_set_has_been_propagated), !, fail)) | |
| 685 | ) | |
| 686 | ; true | |
| 687 | ), | |
| 688 | StackInfo = (PropType,_,_,_), | |
| 689 | ( ( PropType == unit; PropType == branch) % has been propagated by the SAT solver | |
| 690 | -> true | |
| 691 | ; cdclt_sat_solver:log_theory_propagation_sat_stack(SatStateMutable, SatVarName, SatPrologVar, StackInfo) | |
| 692 | ). | |
| 693 | ||
| 694 | log_idl_candidates_from_constraint(Constraint) :- | |
| 695 | get_ids_and_int_constants(Constraint, Ids, IntConstants), | |
| 696 | Ids \== [], | |
| 697 | IntConstants \== [], | |
| 698 | !, | |
| 699 | assert_idl_candidates(Ids, IntConstants). | |
| 700 | log_idl_candidates_from_constraint(_). | |
| 701 | ||
| 702 | assert_idl_candidates([], _). | |
| 703 | assert_idl_candidates([Id|T], IntConstants) :- | |
| 704 | retract(idl_candidate_constants(Id, Candidates)), | |
| 705 | !, | |
| 706 | subtract(IntConstants, Candidates, New), | |
| 707 | append(Candidates, New, NewCandidates), | |
| 708 | asserta(idl_candidate_constants(Id, NewCandidates)), | |
| 709 | assert_idl_candidates(T, IntConstants). | |
| 710 | assert_idl_candidates([Id|T], IntConstants) :- | |
| 711 | asserta(idl_candidate_constants(Id, IntConstants)), | |
| 712 | assert_idl_candidates(T, IntConstants). | |
| 713 | ||
| 714 | get_ids_and_int_constants(b(Node,pred,_), Ids, IntConstants) :- | |
| 715 | comparison_op(Node, Lhs, Rhs), | |
| 716 | get_ids_and_int_constants_expr(Lhs, [], [], IdsAcc, IntConstantsAcc), | |
| 717 | get_ids_and_int_constants_expr(Rhs, IdsAcc, IntConstantsAcc, Ids, IntConstants). | |
| 718 | ||
| 719 | get_ids_and_int_constants_expr(Id, IdsAcc, IntConstantsAcc, [Name|IdsAcc], IntConstantsAcc) :- | |
| 720 | Id = b(identifier(Name),integer,_). | |
| 721 | get_ids_and_int_constants_expr(Int, IdsAcc, IntConstantsAcc, IdsAcc, [Value,Value1|IntConstantsAcc]) :- | |
| 722 | ( Int = b(integer(Value),integer,_) | |
| 723 | ; Int = b(unary_minus(b(integer(Value),integer,_)),integer,_) | |
| 724 | ), | |
| 725 | % select positive and negative value | |
| 726 | Value1 is Value * -1. | |
| 727 | get_ids_and_int_constants_expr(b(Node,integer,_), IdsAcc, IntConstantsAcc, Ids, IntConstants) :- | |
| 728 | arithmetic_expr(Node, Lhs, Rhs), | |
| 729 | get_ids_and_int_constants_expr(Lhs, IdsAcc, IntConstantsAcc, NIdsAcc, NIntConstantsAcc), | |
| 730 | get_ids_and_int_constants_expr(Rhs, NIdsAcc, NIntConstantsAcc, Ids, IntConstants). | |
| 731 | ||
| 732 | comparison_op(less(Lhs,Rhs), Lhs, Rhs). | |
| 733 | comparison_op(less_equal(Lhs,Rhs), Lhs, Rhs). | |
| 734 | comparison_op(equal(Lhs,Rhs), Lhs, Rhs). | |
| 735 | comparison_op(not_equal(Lhs,Rhs), Lhs, Rhs). | |
| 736 | ||
| 737 | arithmetic_expr(minus(Lhs,Rhs), Lhs, Rhs). | |
| 738 | arithmetic_expr(add(Lhs,Rhs), Lhs, Rhs). | |
| 739 | arithmetic_expr(div(Lhs,Rhs), Lhs, Rhs). | |
| 740 | arithmetic_expr(floored_div(Lhs,Rhs), Lhs, Rhs). | |
| 741 | arithmetic_expr(power_of(Lhs,Rhs), Lhs, Rhs). | |
| 742 | arithmetic_expr(multiplication(Lhs,Rhs), Lhs, Rhs). | |
| 743 | arithmetic_expr(modulo(Lhs,Rhs), Lhs, Rhs). | |
| 744 | ||
| 745 | :- block propagate_non_idl_to_idl(?, -, ?). | |
| 746 | propagate_non_idl_to_idl(IdlGraphMut, SatPrologVar, SmtFormula) :- | |
| 747 | ( SatPrologVar == pred_true | |
| 748 | % e.g., x:NAT -> x>=0 | |
| 749 | -> ( infer_constraints_for_idl_solver(SmtFormula, ConjList) | |
| 750 | -> register_constraints(IdlGraphMut, ConjList) | |
| 751 | ; true | |
| 752 | ) | |
| 753 | ; true % TO DO: also propagate pred_false? probably not beneficial | |
| 754 | ). | |
| 755 | ||
| 756 | %% idl_solver_interface(+IdlGraphMut, +SatPrologVar, +Constraint). | |
| 757 | :- block idl_solver_interface(?, -, ?). | |
| 758 | idl_solver_interface(IdlGraphMut, SatPrologVar, [DLConstraint]) :- | |
| 759 | SatPrologVar == pred_true, | |
| 760 | difference_logic_solver:remove_unsat_core, | |
| 761 | cdclt_sat_solver:remove_idl_unsat_core, | |
| 762 | ( difference_logic_solver:register_constraint(IdlGraphMut, DLConstraint) | |
| 763 | ; \+ cdclt_sat_solver:is_backjumping, | |
| 764 | propagate_idl_unsat_core_to_sat_solver(IdlGraphMut), | |
| 765 | fail | |
| 766 | ). | |
| 767 | idl_solver_interface(IdlGraphMut, SatPrologVar, [DLConstraint]) :- | |
| 768 | SatPrologVar == pred_false, | |
| 769 | difference_logic_solver:remove_unsat_core, | |
| 770 | cdclt_sat_solver:remove_idl_unsat_core, | |
| 771 | ( difference_logic_solver:register_constraint(IdlGraphMut, b(negation(DLConstraint),pred,[])) | |
| 772 | ; \+ cdclt_sat_solver:is_backjumping, | |
| 773 | propagate_idl_unsat_core_to_sat_solver(IdlGraphMut), | |
| 774 | fail | |
| 775 | ). | |
| 776 | ||
| 777 | propagate_idl_unsat_core_to_sat_solver(IdlGraphMut) :- | |
| 778 | difference_logic_solver:get_solver_result(IdlGraphMut, Res), | |
| 779 | Res == contradiction_found, | |
| 780 | difference_logic_solver:get_unsat_core(IdlCore), | |
| 781 | difference_logic_solver:remove_unsat_core, | |
| 782 | cdclt_sat_solver:store_idl_unsat_core(IdlCore), !. | |
| 783 | ||
| 784 | %% get_idl_solution_bindings(+IdlGraphMut, -Bindings). | |
| 785 | get_idl_solution_bindings(IdlGraphMut, Bindings) :- | |
| 786 | get_solver_result(IdlGraphMut, SolverResult), | |
| 787 | SolverResult = solution(Bindings). | |
| 788 | ||
| 789 | %% propagate_idl_solution_to_bindings_bt(+IdlGraphMut, +SmtBindings). | |
| 790 | propagate_idl_solution_to_bindings_bt(IdlGraphMut, SmtBindings) :- | |
| 791 | get_preference(cdclt_use_idl_theory_solver,true), | |
| 792 | !, | |
| 793 | % TO DO: try the current solution first | |
| 794 | ( get_candidate_bounds_from_non_idl_constraints(CandidateTuples), | |
| 795 | difference_logic_solver:try_candidate_bounds(IdlGraphMut, CandidateTuples, Result) | |
| 796 | ; difference_logic_solver:get_all_solutions_on_bt(IdlGraphMut, Result) | |
| 797 | ), | |
| 798 | Result = solution(Bindings), | |
| 799 | set_bindings(Bindings, SmtBindings). | |
| 800 | propagate_idl_solution_to_bindings_bt(_, _). | |
| 801 | ||
| 802 | %% propagate_fd_bounds_to_idl_solver(+IdlGraphMut, +SmtBindings). | |
| 803 | propagate_fd_bounds_to_idl_solver(IdlGraphMut, SmtBindings) :- | |
| 804 | get_preference(cdclt_use_idl_theory_solver,true), | |
| 805 | !, | |
| 806 | difference_logic_solver:get_registered_vars(IdlGraphMut, Vars), | |
| 807 | get_constraints_from_fd_bounds(SmtBindings, Vars, [], ConjList), | |
| 808 | difference_logic_solver:register_constraints(IdlGraphMut, ConjList). | |
| 809 | propagate_fd_bounds_to_idl_solver(_, _). | |
| 810 | ||
| 811 | get_constraints_from_fd_bounds(_, [], Acc, Acc). | |
| 812 | get_constraints_from_fd_bounds(SmtBindings, [Var|T], Acc, ConjList) :- | |
| 813 | memberchk(bind(Var, int(Int)), SmtBindings), | |
| 814 | integer(Int), | |
| 815 | !, | |
| 816 | Zero = b(identifier('_zero'),integer,[]), | |
| 817 | % Var = Int but represented as Var - _zero <= Int & _zero - Var <= -Int | |
| 818 | Conj1 = b(less_equal(b(minus(b(identifier(Var),integer,[]),Zero),integer,[]),b(integer(Int),integer,[])),pred,[]), | |
| 819 | Conj2 = b(less_equal(b(minus(Zero,b(identifier(Var),integer,[])),integer,[]),b(unary_minus(b(integer(Int),integer,[])),integer,[])),pred,[]), | |
| 820 | get_constraints_from_fd_bounds(SmtBindings, T, [Conj1,Conj2|Acc], ConjList). | |
| 821 | get_constraints_from_fd_bounds(SmtBindings, [Var|T], Acc, ConjList) :- | |
| 822 | memberchk(bind(Var, int(FDVar)), SmtBindings), | |
| 823 | fd_var(FDVar), | |
| 824 | fd_min(FDVar, Min), | |
| 825 | fd_max(FDVar, Max), | |
| 826 | !, | |
| 827 | ( var_geq_min(Var, Min, MinBound) | |
| 828 | -> ( var_leq_max(Var, Max, MaxBound) | |
| 829 | -> NAcc = [MinBound,MaxBound|Acc] | |
| 830 | ; NAcc = [MinBound|Acc] | |
| 831 | ) | |
| 832 | ; ( var_leq_max(Var, Max, MaxBound) | |
| 833 | -> NAcc = [MaxBound|Acc] | |
| 834 | ; NAcc = Acc | |
| 835 | ) | |
| 836 | ), | |
| 837 | get_constraints_from_fd_bounds(SmtBindings, T, NAcc, ConjList). | |
| 838 | get_constraints_from_fd_bounds(SmtBindings, [_|T], Acc, ConjList) :- | |
| 839 | get_constraints_from_fd_bounds(SmtBindings, T, Acc, ConjList). | |
| 840 | ||
| 841 | var_geq_min(Var, Min, MinBound) :- | |
| 842 | % -v <= -min | |
| 843 | integer(Min), | |
| 844 | Min1 is Min * -1, | |
| 845 | Minus = b(minus(b(identifier('_zero'),integer,[]),b(identifier(Var),integer,[])),integer,[]), | |
| 846 | MinBound = b(less_equal(Minus,b(integer(Min1),integer,[])),pred,[]). | |
| 847 | ||
| 848 | var_leq_max(Var, Max, MaxBound) :- | |
| 849 | % v <= Max | |
| 850 | integer(Max), | |
| 851 | MaxBound = b(less_equal(b(identifier(Var),integer,[]),b(integer(Max),integer,[])),pred,[]). | |
| 852 | ||
| 853 | get_candidate_bounds_from_non_idl_constraints(CandidateTuples) :- | |
| 854 | findall((VarName,Candidates), idl_candidate_constants(VarName, Candidates), CandidateTuples). | |
| 855 | ||
| 856 | set_bindings([], _). | |
| 857 | set_bindings([binding(VarName,Val,_)|T], SmtBindings) :- | |
| 858 | member(bind(VarName,Val), SmtBindings), | |
| 859 | set_bindings(T, SmtBindings). | |
| 860 | ||
| 861 | %% exclude_idl_solution(+IdlBindings, +SmtBindings, +StateTypedIds, +WfStoreSmt). | |
| 862 | exclude_idl_solution(IdlBindings, SmtBindings, StateTypedIds, WfStoreSmt) :- | |
| 863 | initial_solution, | |
| 864 | exclude_solution(IdlBindings, StateTypedIds, Exclusion), | |
| 865 | b_test_boolean_expression(Exclusion, SmtBindings, [], WfStoreSmt). | |
| 866 | exclude_idl_solution(_, _, _, _) :- | |
| 867 | \+ initial_solution. | |
| 868 | ||
| 869 | %% exclude_solution(+Bindings, +StateTypedIds, -Exclusion). | |
| 870 | exclude_solution([], _, b(truth,pred,[])). | |
| 871 | exclude_solution([Binding|T], StateTypedIds, Exclusion) :- | |
| 872 | get_equality_from_binding(Binding, StateTypedIds, EQ), | |
| 873 | exclude_solution(T, EQ, StateTypedIds, Exclusion). | |
| 874 | ||
| 875 | exclude_solution([], Acc, _, Acc). | |
| 876 | exclude_solution([Binding|T], Acc, StateTypedIds, Exclusion) :- | |
| 877 | get_equality_from_binding(Binding, StateTypedIds, EQ), | |
| 878 | safe_create_texpr(negation(EQ), pred, [], Neg), | |
| 879 | safe_create_texpr(disjunct(Neg,Acc), pred, [], NAcc), | |
| 880 | exclude_solution(T, NAcc, StateTypedIds, Exclusion). | |
| 881 | ||
| 882 | %% get_equality_from_binding(+Binding, +TypedIds, -EQ). | |
| 883 | get_equality_from_binding(Binding, TypedIds, EQ) :- | |
| 884 | ( Binding = bind(VarName,Val) | |
| 885 | ; Binding = binding(VarName,Val,_) | |
| 886 | ), | |
| 887 | ( member(b(identifier(VarName),TType,_), TypedIds) | |
| 888 | -> Type = TType | |
| 889 | ; infer_value_type(Val, Type) | |
| 890 | ), | |
| 891 | ( contains_any(Type) | |
| 892 | -> add_error_and_fail(get_equality_from_binding, 'Cannot infer type of identifier to setup state for ', VarName) | |
| 893 | ; EQ = b(equal(b(identifier(VarName),Type,[]),b(value(Val),Type,[])),pred,[]) | |
| 894 | ). | |
| 895 | ||
| 896 | %% setup_theory_wf_store(+SmtVars, -SmtBindings, -WFStoreSMT). | |
| 897 | setup_theory_wf_store(SmtVars, SmtBindings, WFStoreSMT) :- | |
| 898 | set_up_typed_localstate(SmtVars, _, SmtTypedVals, [], SmtBindings, positive), | |
| 899 | init_wait_flags(WFStoreSMT, [wf_smt_cdcl]), | |
| 900 | b_tighter_enumerate_all_values(SmtTypedVals, WFStoreSMT),!. | |
| 901 | ||
| 902 | log_solution :- | |
| 903 | \+ initial_solution, | |
| 904 | !, | |
| 905 | asserta(initial_solution). | |
| 906 | log_solution. | |
| 907 | ||
| 908 | ground_state_binding(Binding) :- | |
| 909 | Binding = bind(_,_), | |
| 910 | ground(Binding). | |
| 911 | ||
| 912 | get_typed_ids_in_scope_of_pred(Pred, TypedIds) :- | |
| 913 | find_typed_identifier_uses(Pred, PTypedIds), | |
| 914 | add_typed_ids_in_scope(PTypedIds, TypedIds). | |
| 915 | ||
| 916 | add_typed_ids_in_scope(PTypedIds, TypedIds) :- | |
| 917 | b_get_machine_variables(MachineVars), | |
| 918 | b_get_machine_constants(MachineConstants), | |
| 919 | append([PTypedIds,MachineVars,MachineConstants], TypedIds). | |
| 920 | ||
| 921 | get_equalities_from_bindings([], _, b(truth,pred,[])). | |
| 922 | get_equalities_from_bindings([bind(Id,Val)|T], TypedIds, EqConj) :- | |
| 923 | get_equality_from_binding(bind(Id,Val), TypedIds, Eq), | |
| 924 | get_equalities_from_bindings(T, TypedIds, Eq, EqConj). | |
| 925 | ||
| 926 | get_equalities_from_bindings([], _, EqConj, EqConj). | |
| 927 | get_equalities_from_bindings([bind(Id,Val)|T], TypedIds, EqAcc, EqConj) :- | |
| 928 | get_equality_from_binding(bind(Id,Val), TypedIds, Eq), | |
| 929 | safe_create_texpr(conjunct(Eq,EqAcc), pred, [], Conj), | |
| 930 | get_equalities_from_bindings(T, TypedIds, Conj, EqConj). | |
| 931 | ||
| 932 | % only for computing the amount of sat variables of benchmarks | |
| 933 | get_amount_of_sat_variables(Pred, AmountOfSatVars) :- | |
| 934 | simplify_negation(Pred, SPred), | |
| 935 | reset_optimizer_state, | |
| 936 | assert_ground_id_values(0, SPred), | |
| 937 | replace_ids_with_ground_values(SPred, 0, [], AstCardOpt), | |
| 938 | precompute_values(AstCardOpt, [instantiate_quantifier_limit(0)], AstPrecomputed), | |
| 939 | (clean_up_pred(AstPrecomputed, _, CleanPred) -> true | |
| 940 | ; add_internal_error('Clean up failed ', cdclt), CleanPred=AstPrecomputed), | |
| 941 | ( CleanPred = b(truth,pred,_) | |
| 942 | -> AmountOfSatVars = 1 | |
| 943 | ; CleanPred = b(falsity,pred,_) | |
| 944 | -> AmountOfSatVars = 1 | |
| 945 | ; preprocess_predicate(false, false, CleanPred, LiftedPred, _, _), | |
| 946 | predicate_to_sat(normal_make_wd, LiftedPred, _, _, _, SatVars1) | |
| 947 | ), | |
| 948 | length(SatVars1, AmountOfSatVars). | |
| 949 | ||
| 950 | /* currently not used: | |
| 951 | unfold_let(Pred, Unfolded) :- | |
| 952 | unfold_let([], Pred, Unfolded). | |
| 953 | ||
| 954 | unfold_let(IdTuples, b(identifier(Name),_,_), Replacement) :- | |
| 955 | member((b(identifier(Name),_,_),Ast), IdTuples), | |
| 956 | !, | |
| 957 | Replacement = Ast. | |
| 958 | unfold_let(IdTuples, b(Node,_,_), Replacement) :- | |
| 959 | ( Node = let_predicate(Ids,Asts,Body) | |
| 960 | ; Node = let_expression(Ids,Asts,Body) | |
| 961 | ), | |
| 962 | !, | |
| 963 | maplist(unfold_let(IdTuples), Asts, NAsts), | |
| 964 | zip_acc(Ids, NAsts, IdTuples, NewIdTuples), % TO DO: ideally remove IdTuples which are clashing with Ids | |
| 965 | unfold_let(NewIdTuples, Body, Replacement). | |
| 966 | unfold_let(IdTuples, b(Node,Type,Info), Replacement) :- | |
| 967 | syntaxtransformation(Node,Subs,Names,NSubs,NewNode),!, | |
| 968 | Replacement = b(NewNode,Type,Info), | |
| 969 | exclude(hidden_by_local_var(Names),IdTuples,NIdTuples), | |
| 970 | l_unfold_let(Subs,NIdTuples,NSubs). | |
| 971 | unfold_let(I, Ast, _) :- add_internal_error('Not a typed expression:', unfold_let(I,Ast,_)),fail. | |
| 972 | ||
| 973 | hidden_by_local_var(Names,(TID,_)) :- def_get_texpr_id(TID,ID), member(ID,Names). | |
| 974 | ||
| 975 | l_unfold_let([],_,[]). | |
| 976 | l_unfold_let([H|T],IdTuples,[NH|NT]) :- | |
| 977 | unfold_let(H,IdTuples,NH), | |
| 978 | l_unfold_let(T,IdTuples,NT). | |
| 979 | ||
| 980 | zip_acc([], [], Acc, Acc). | |
| 981 | zip_acc([A|T1], [B|T2], Acc, Zipped) :- | |
| 982 | zip_acc(T1, T2, [(A,B)|Acc], Zipped). % TO DO: should we just store ID and not typed ID in A? | |
| 983 | ||
| 984 | */ | |
| 985 | ||
| 986 | %% For debugging: transform bool formula back to SMT | |
| 987 | /*bool_formula_to_smt(b(truth,pred,I), b(truth,pred,I)). | |
| 988 | bool_formula_to_smt(b(falsity,pred,I), b(falsity,pred,I)). | |
| 989 | bool_formula_to_smt(b(equal(Id,Bool),pred,_), Smt) :- | |
| 990 | Id = b(identifier(_),boolean,IdInfo), | |
| 991 | member(smt_formula(TSmt), IdInfo), | |
| 992 | !, | |
| 993 | ( Bool = b(boolean_true,boolean,_) | |
| 994 | -> Smt = TSmt | |
| 995 | ; Smt = b(negation(TSmt),pred,[]) | |
| 996 | ). | |
| 997 | bool_formula_to_smt(b(equal(Id,Bool),pred,I), Out) :- % bool ids from Tseitin optimization | |
| 998 | !, | |
| 999 | Out = b(equal(Id,Bool),pred,I). | |
| 1000 | bool_formula_to_smt(b(conjunct(A,B),pred,I), Out) :- | |
| 1001 | !, | |
| 1002 | bool_formula_to_smt(A, SmtA), | |
| 1003 | bool_formula_to_smt(B, SmtB), | |
| 1004 | Out = b(conjunct(SmtA,SmtB),pred,I). | |
| 1005 | bool_formula_to_smt(b(disjunct(A,B),pred,I), Out) :- | |
| 1006 | !, | |
| 1007 | bool_formula_to_smt(A, SmtA), | |
| 1008 | bool_formula_to_smt(B, SmtB), | |
| 1009 | Out = b(disjunct(SmtA,SmtB),pred,I).*/ | |
| 1010 | %% |