10 (* compound operator *)
11 type compound_operator = op list
13 let string_of_cop op =
14 if op = [] then "id" else String.concat "" (List.map string_of_op op)
16 let dot_of_cop op = "\"" ^ string_of_cop op ^ "\""
18 let rec matita_of_cop v =
21 | I::tl -> "i (" ^ matita_of_cop v tl ^ ")"
22 | C::tl -> "c (" ^ matita_of_cop v tl ^ ")"
23 | M::tl -> "m (" ^ matita_of_cop v tl ^ ")"
25 (* representative, other elements in the equivalence class,
26 leq classes, geq classes *)
27 type equivalence_class =
28 compound_operator * compound_operator list *
29 equivalence_class list ref * equivalence_class list ref
31 let string_of_equivalence_class (repr,others,leq,_) =
32 String.concat " = " (List.map string_of_cop (repr::others)) ^
37 (function (repr',_,_,_) ->
38 string_of_cop repr ^ " <= " ^ string_of_cop repr') !leq)
42 let dot_of_equivalence_class (repr,others,leq,_) =
44 let eq = String.concat " = " (List.map string_of_cop (repr::others)) in
45 dot_of_cop repr ^ "[label=\"" ^ eq ^ "\"];" ^
46 if !leq = [] then "" else "\n"
47 else if !leq = [] then
53 (function (repr',_,_,_) ->
54 dot_of_cop repr' ^ " -> " ^ dot_of_cop repr ^ ";") !leq)
56 (* set of equivalence classes *)
57 type set = equivalence_class list
60 String.concat "\n" (List.map string_of_equivalence_class s)
62 let ps_of_set (to_be_considered,under_consideration,news) ?processing s =
63 let ch = open_out "xxx.dot" in
64 output_string ch "digraph G {\n";
65 (match under_consideration with
68 output_string ch (dot_of_cop repr ^ " [color=yellow];"));
70 (function repr -> output_string ch (dot_of_cop repr ^ " [color=green];")
73 (function repr -> output_string ch (dot_of_cop repr ^ " [color=navy];")
75 output_string ch (String.concat "\n" (List.map dot_of_equivalence_class s));
76 output_string ch "\n";
77 (match processing with
79 | Some (repr,rel,repr') ->
80 output_string ch (dot_of_cop repr ^ " [color=red];");
82 (dot_of_cop repr' ^ " -> " ^ dot_of_cop repr ^
84 (if rel="=" then "arrowhead=none " else "") ^
86 output_string ch "}\n";
88 ignore (Unix.system "dot -Tps xxx.dot > xxx.ps")
90 let test to_be_considered_and_now set rel candidate repr =
91 ps_of_set to_be_considered_and_now ~processing:(candidate,rel,repr) set;
93 (string_of_cop candidate ^ " " ^ rel ^ " " ^ string_of_cop repr ^ "? ");
95 assert (Unix.system "cp formal_topology.ma xxx.ma" = Unix.WEXITED 0);
96 let ch = open_out_gen [Open_append] 0 "xxx.ma" in
99 (function (repr,others,leq,_) ->
104 ("axiom ax" ^ string_of_int !i ^
106 matita_of_cop "A" repr ^ " = " ^ matita_of_cop "A" repr' ^ ".\n");
109 (function (repr',_,_,_) ->
112 ("axiom ax" ^ string_of_int !i ^
114 matita_of_cop "A" repr ^ " ⊆ " ^ matita_of_cop "A" repr' ^ ".\n");
118 ("theorem foo: \\forall A." ^ matita_of_cop "A" candidate ^ " " ^ rel ^ " " ^
119 matita_of_cop "A" repr ^ ". intros; auto size=6 depth=4. qed.\n");
122 Unix.system "../../../matitac.opt xxx.ma >> log 2>&1" = Unix.WEXITED 0
124 print_endline (if res then "y" else "n");
127 let normalize to_be_considered_and_now candidate set =
130 [] -> raise Not_found
131 | (repr,others,leq,geq) as eqclass :: tl ->
132 if test to_be_considered_and_now set "=" candidate repr then
133 (repr,others@[candidate],leq,geq)::tl
140 let locate to_be_considered_and_now ((repr,_,leq,geq) as node) set =
144 | (repr',_,leq',geq') as node' :: tl ->
145 if repr = repr' then ()
146 else if test to_be_considered_and_now set "⊆" repr repr' then
148 leq := node' :: !leq;
149 geq' := node :: !geq'
151 else if test to_be_considered_and_now set "⊆" repr' repr then
153 geq := node' :: !geq;
154 leq' := node :: !leq'
161 let analyze_one to_be_considered repr hecandidate (news,set) =
162 let candidate = hecandidate::repr in
163 if List.length (List.filter ((=) M) candidate) > 1 then
167 let set = normalize (to_be_considered,Some repr,news) candidate set in
173 let node = candidate,[],leq,geq in
174 let set = node::set in
175 locate (to_be_considered,Some repr,news) node set;
179 let rec explore i set news =
180 let rec aux news set =
185 List.fold_right (analyze_one tl repr) [I;C;M] (news,set)
189 let news,set = aux [] set news in
192 print_endline ("PUNTO FISSO RAGGIUNTO! i=" ^ string_of_int i);
193 print_endline (string_of_set set ^ "\n----------------");
194 ps_of_set ([],None,[]) set
198 print_endline ("NUOVA ITERAZIONE, i=" ^ string_of_int i);
199 print_endline (string_of_set set ^ "\n----------------");
200 explore (i+1) set news
204 let set = [id,[],ref [], ref []] in
205 print_endline ("PRIMA ITERAZIONE, i=0, j=0");
206 print_endline (string_of_set set ^ "\n----------------");
207 ignore (Unix.system "rm -f log");
208 ps_of_set ([id],None,[]) set;