%%% -*- Mode: Prolog; -*- %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Date: 2010-09-28 21:04:43 +0200 (Tue, 28 Sep 2010) $ % $Revision: 4838 $ % % This file is part of ProbLog % http://dtai.cs.kuleuven.be/problog % % ProbLog was developed at Katholieke Universiteit Leuven % % Copyright 2008, 2009, 2010 % Katholieke Universiteit Leuven % % Main authors of this file: % Theofrastos Mantadelis % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Artistic License 2.0 % % Copyright (c) 2000-2006, The Perl Foundation. % % Everyone is permitted to copy and distribute verbatim copies of this % license document, but changing it is not allowed. Preamble % % This license establishes the terms under which a given free software % Package may be copied, modified, distributed, and/or % redistributed. The intent is that the Copyright Holder maintains some % artistic control over the development of that Package while still % keeping the Package available as open source and free software. % % You are always permitted to make arrangements wholly outside of this % license directly with the Copyright Holder of a given Package. If the % terms of this license do not permit the full use that you propose to % make of the Package, you should contact the Copyright Holder and seek % a different licensing arrangement. Definitions % % "Copyright Holder" means the individual(s) or organization(s) named in % the copyright notice for the entire Package. % % "Contributor" means any party that has contributed code or other % material to the Package, in accordance with the Copyright Holder's % procedures. % % "You" and "your" means any person who would like to copy, distribute, % or modify the Package. % % "Package" means the collection of files distributed by the Copyright % Holder, and derivatives of that collection and/or of those files. A % given Package may consist of either the Standard Version, or a % Modified Version. % % "Distribute" means providing a copy of the Package or making it % accessible to anyone else, or in the case of a company or % organization, to others outside of your company or organization. % % "Distributor Fee" means any fee that you charge for Distributing this % Package or providing support for this Package to another party. It % does not mean licensing fees. % % "Standard Version" refers to the Package if it has not been modified, % or has been modified only in ways explicitly requested by the % Copyright Holder. % % "Modified Version" means the Package, if it has been changed, and such % changes were not explicitly requested by the Copyright Holder. % % "Original License" means this Artistic License as Distributed with the % Standard Version of the Package, in its current version or as it may % be modified by The Perl Foundation in the future. % % "Source" form means the source code, documentation source, and % configuration files for the Package. % % "Compiled" form means the compiled bytecode, object code, binary, or % any other form resulting from mechanical transformation or translation % of the Source form. % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Permission for Use and Modification Without Distribution % % (1) You are permitted to use the Standard Version and create and use % Modified Versions for any purpose without restriction, provided that % you do not Distribute the Modified Version. % % Permissions for Redistribution of the Standard Version % % (2) You may Distribute verbatim copies of the Source form of the % Standard Version of this Package in any medium without restriction, % either gratis or for a Distributor Fee, provided that you duplicate % all of the original copyright notices and associated disclaimers. At % your discretion, such verbatim copies may or may not include a % Compiled form of the Package. % % (3) You may apply any bug fixes, portability changes, and other % modifications made available from the Copyright Holder. The resulting % Package will still be considered the Standard Version, and as such % will be subject to the Original License. % % Distribution of Modified Versions of the Package as Source % % (4) You may Distribute your Modified Version as Source (either gratis % or for a Distributor Fee, and with or without a Compiled form of the % Modified Version) provided that you clearly document how it differs % from the Standard Version, including, but not limited to, documenting % any non-standard features, executables, or modules, and provided that % you do at least ONE of the following: % % (a) make the Modified Version available to the Copyright Holder of the % Standard Version, under the Original License, so that the Copyright % Holder may include your modifications in the Standard Version. (b) % ensure that installation of your Modified Version does not prevent the % user installing or running the Standard Version. In addition, the % modified Version must bear a name that is different from the name of % the Standard Version. (c) allow anyone who receives a copy of the % Modified Version to make the Source form of the Modified Version % available to others under (i) the Original License or (ii) a license % that permits the licensee to freely copy, modify and redistribute the % Modified Version using the same licensing terms that apply to the copy % that the licensee received, and requires that the Source form of the % Modified Version, and of any works derived from it, be made freely % available in that license fees are prohibited but Distributor Fees are % allowed. % % Distribution of Compiled Forms of the Standard Version or % Modified Versions without the Source % % (5) You may Distribute Compiled forms of the Standard Version without % the Source, provided that you include complete instructions on how to % get the Source of the Standard Version. Such instructions must be % valid at the time of your distribution. If these instructions, at any % time while you are carrying out such distribution, become invalid, you % must provide new instructions on demand or cease further % distribution. If you provide valid instructions or cease distribution % within thirty days after you become aware that the instructions are % invalid, then you do not forfeit any of your rights under this % license. % % (6) You may Distribute a Modified Version in Compiled form without the % Source, provided that you comply with Section 4 with respect to the % Source of the Modified Version. % % Aggregating or Linking the Package % % (7) You may aggregate the Package (either the Standard Version or % Modified Version) with other packages and Distribute the resulting % aggregation provided that you do not charge a licensing fee for the % Package. Distributor Fees are permitted, and licensing fees for other % components in the aggregation are permitted. The terms of this license % apply to the use and Distribution of the Standard or Modified Versions % as included in the aggregation. % % (8) You are permitted to link Modified and Standard Versions with % other works, to embed the Package in a larger work of your own, or to % build stand-alone binary or bytecode versions of applications that % include the Package, and Distribute the result without restriction, % provided the result does not expose a direct interface to the Package. % % Items That are Not Considered Part of a Modified Version % % (9) Works (including, but not limited to, modules and scripts) that % merely extend or make use of the Package, do not, by themselves, cause % the Package to be a Modified Version. In addition, such works are not % considered parts of the Package itself, and are not subject to the % terms of this license. % % General Provisions % % (10) Any use, modification, and distribution of the Standard or % Modified Versions is governed by this Artistic License. By using, % modifying or distributing the Package, you accept this license. Do not % use, modify, or distribute the Package, if you do not accept this % license. % % (11) If your Modified Version has been derived from a Modified Version % made by someone other than you, you are nevertheless required to % ensure that your Modified Version complies with the requirements of % this license. % % (12) This license does not grant you the right to use any trademark, % service mark, tradename, or logo of the Copyright Holder. % % (13) This license includes the non-exclusive, worldwide, % free-of-charge patent license to make, have made, use, offer to sell, % sell, import and otherwise transfer the Package with respect to any % patent claims licensable by the Copyright Holder that are necessarily % infringed by the Package. If you institute patent litigation % (including a cross-claim or counterclaim) against any party alleging % that the Package constitutes direct or contributory patent % infringement, then this Artistic License to you shall terminate on the % date that such litigation is filed. % % (14) Disclaimer of Warranty: THE PACKAGE IS PROVIDED BY THE COPYRIGHT % HOLDER AND CONTRIBUTORS "AS IS' AND WITHOUT ANY EXPRESS OR IMPLIED % WARRANTIES. THE IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A % PARTICULAR PURPOSE, OR NON-INFRINGEMENT ARE DISCLAIMED TO THE EXTENT % PERMITTED BY YOUR LOCAL LAW. UNLESS REQUIRED BY LAW, NO COPYRIGHT % HOLDER OR CONTRIBUTOR WILL BE LIABLE FOR ANY DIRECT, INDIRECT, % INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING IN ANY WAY OUT OF THE USE % OF THE PACKAGE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% :- module(variable_elimination, [trie_check_for_and_cluster/1, trie_replace_and_cluster/2, clean_up/0, variable_elimination_stats/3]). :- use_module(library(lists), [append/3, delete/3, memberchk/2, reverse/2]). :- use_module(library(tries)). :- use_module('flags', _, [problog_define_flag/5]). :- initialization(( nb_setval(prob_fact_count, 0), problog_define_flag(variable_elimination, problog_flag_validate_boolean, 'enable variable elimination', false, variable_elimination) )). bit_encode(L, ON):- bit_encode(L, ON, 0). bit_encode([], ON, ON):-!. bit_encode([PF|T], ON, Acc):- (recorded(variable_elimination, prob_fact(PF, ID), _) -> true ; nb_getval(prob_fact_count, ID), NID is ID + 1, nb_setval(prob_fact_count, NID), recordz(variable_elimination, prob_fact(PF, ID), _) ), NAcc is Acc \/ (1 << ID), bit_encode(T, ON, NAcc). bit_decode(ON, L):- bit_decode(ON, 0, L). bit_decode(_, ID, []):- nb_getval(prob_fact_count, ID), !. bit_decode(ON, ID, [PF|L]):- 0 < ON /\ (1 << ID), recorded(variable_elimination, prob_fact(PF, ID), _), NID is ID + 1, bit_decode(ON, NID, L). bit_decode(ON, ID, L):- NID is ID + 1, bit_decode(ON, NID, L). update_table(T, ON, NT):- update_table(T, ON, NT, 0). update_table([], _ON, [], _). update_table([H|T], ON, [NH|NT], Row):- 0 is ON /\ (1 << Row), !, NH is H /\ \ ON, % this is an optional improvement NRow is Row + 1, update_table(T, ON, NT, NRow). update_table([H|T], ON, [NH|NT], Row):- NH is H /\ ON, NRow is Row + 1, update_table(T, ON, NT, NRow). make_mask(FromBit, ToBit, Mask):- Mask is (1 << (ToBit + 1) - 1) - (1 << FromBit - 1). make_table(_, 0, []):-!. make_table(ON, T, [ON|L]):- NT is T - 1, make_table(ON, NT, L). modify_table(L, OT, NT):- nb_getval(prob_fact_count, OLS), bit_encode(L, ON), nb_getval(prob_fact_count, NLS), update_table(OT, ON, L1), D is NLS - OLS, make_mask(OLS, NLS, M), NON is M /\ ON, make_table(NON, D, L2), append(L1, L2, NT). examin(T):- examin(T, 0). examin([], _Row). examin([H|T], Row):- N is 1 << Row, 0 is H /\ (N - 1), 0 < H - N, !, bit_decode(H, L), calc_prob_AND_cluster(L, P), make_prob_fact(L, P, ID), recordz(variable_elimination, and_cluster(L, ID), _), NRow is Row + 1, examin(T, NRow). examin([_H|T], Row):- NRow is Row + 1, examin(T, NRow). trie_check_for_and_cluster(T):- tries:trie_traverse_first(T, E), !, trie_check_for_and_cluster(E, []). trie_check_for_and_cluster(_T). trie_check_for_and_cluster(E, T):- tries:trie_traverse_next(E, N), !, tries:trie_get_entry(E, L), modify_table(L, T, NT), trie_check_for_and_cluster(N, NT). trie_check_for_and_cluster(E, T):- tries:trie_get_entry(E, L), modify_table(L, T, NT), examin(NT), !. trie_replace_and_cluster(To, Tn):- tries:trie_open(Tn), trie_replace_and_cluster_do(To, Tn). trie_replace_and_cluster_do(To, Tn):- trie_traverse(To, E), trie_get_entry(E, L), findall(Cluster/VarName, recorded(variable_elimination, and_cluster(Cluster, VarName), _), Clusters), foreach(Clusters, NewL, L), trie_put_entry(Tn, NewL, _), fail. trie_replace_and_cluster_do(_To, _Tn). foreach([], L, L). foreach([Cluster/VarName|Rest], L, Acc):- check_replace_cluster(Cluster, VarName, Acc, NL), foreach(Rest, L, NL). check_replace_cluster(Cluster, _VarName, L, L):- nocluster(Cluster, L), !. check_replace_cluster(Cluster, VarName, L, NewL):- replace_cluster(Cluster, VarName, L, NewL). replace_cluster(Cluster, VarName, L, Res):- first_cluster_element(L, Cluster, First), replace(L, First, VarName, NL), delete(Cluster, First, RestCluster), eliminate_list(RestCluster, NL, Res), !. replace_cluster(Cluster, VarName, _L, _Res):- throw(error(Cluster, VarName)). replace_cluster2(Cluster, VarName, L, Res):- eliminate_list(Cluster, L, NL), append(NL, [VarName], Res), !. replace_cluster2(Cluster, VarName, _L, _Res):- throw(error(Cluster, VarName)). replace_cluster3(Cluster, VarName, L, Res):- last_cluster_element(L, Cluster, Last), replace(L, Last, VarName, NL), delete(Cluster, Last, RestCluster), eliminate_list(RestCluster, NL, Res), !. replace_cluster3(Cluster, VarName, _L, _Res):- throw(error(Cluster, VarName)). first_cluster_element([], _, _). first_cluster_element([H|_T], Cluster, H):- memberchk(H, Cluster), !. first_cluster_element([_H|T], Cluster, R):- first_cluster_element(T, Cluster, R). last_cluster_element(L, Cluster, R):- reverse(L, RL), first_cluster_element(RL, Cluster, R). nocluster([], _). nocluster([H|T], L):- \+ memberchk(H, L), nocluster(T, L). eliminate_list([], L, L). eliminate_list([H|T], L, Res):- memberchk(H, L), delete(L, H, NL), eliminate_list(T, NL, Res). replace([], _, _, []). replace([H|T], H, NH, [NH|NT]):- replace(T, H, NH, NT). replace([H|T], R, NR, [H|NT]):- H \== R, replace(T, R, NR, NT). clean_up:- eraseall(variable_elimination), nb_setval(prob_fact_count, 0). variable_elimination_stats(Clusters, OrigPF, CompPF):- nb_getval(prob_fact_count, OrigPF), findall(L, (recorded(variable_elimination, and_cluster(Cluster, _), _), length(Cluster, L)), LL), sum_list(LL, EliminatedPF), length(LL, Clusters), CompPF is OrigPF - EliminatedPF + Clusters. calc_prob_AND_cluster(L, P):- multiply_list(L, P, 1.0). multiply_list([], P, P). multiply_list([H|T], Pr, A):- problog:get_fact_probability(H, P), number(P), NA is A * P, multiply_list(T, Pr, NA). make_prob_fact(L, P, ID):- (clause(problog:problog_var_elimination(ID, L, _), true) -> true ; problog:probclause_id(ID), assert_static(problog:prob_for_id(ID, P, _)), (clause(problog:problog_predicate(var_elimination, 1), true) -> true ; assertz(problog:problog_predicate(var_elimination, 1)) ), assertz(problog:problog_var_elimination(ID, L, P)) ).