]> matita.cs.unibo.it Git - helm.git/blobdiff - helm/software/components/tactics/closeCoercionGraph.ml
1. composition of coercions with saturations > 0 is now implemented
[helm.git] / helm / software / components / tactics / closeCoercionGraph.ml
index a009944f8c26fe4b14caf41059507e0022e72ddf..68e4e7cb2c2b2c176e274e76b6257ea7e09cf2a0 100644 (file)
 
 (* $Id: cicCoercion.ml 7077 2006-12-05 15:44:54Z fguidi $ *)
 
-let debug = false
+let debug = false 
 let debug_print s = if debug then prerr_endline (Lazy.force s) else ()
 
 (* given the new coercion uri from src to tgt returns the list 
- * of new coercions to create. hte list elements are
+ * of new coercions to create. the list elements are
  * (source, list of coercions to follow, target)
  *)
 let get_closure_coercions src tgt uri coercions =
-  let eq_carr s t = 
+  let eq_carr ?exact s t = 
     try
-      CoercDb.eq_carr s t
+      CoercDb.eq_carr ?exact s t
     with
     | CoercDb.EqCarrNotImplemented _ | CoercDb.EqCarrOnNonMetaClosed -> false
   in
@@ -52,15 +52,23 @@ let get_closure_coercions src tgt uri coercions =
           coercions 
       in
         (HExtlib.flatten_map 
-          (fun (_,t,ul) -> List.map (fun u -> src,[uri; u],t) ul) c_from_tgt) @
+          (fun (_,t,ul) -> 
+             if eq_carr ~exact:true src t then [] else
+             List.map (fun u -> src,[uri; u],t) ul) c_from_tgt) @
         (HExtlib.flatten_map 
-          (fun (s,_,ul) -> List.map (fun u -> s,[u; uri],tgt) ul) c_to_src) @
+          (fun (s,_,ul) -> 
+             if eq_carr ~exact:true s tgt then [] else
+             List.map (fun u -> s,[u; uri],tgt) ul) c_to_src) @
         (HExtlib.flatten_map 
           (fun (s,_,u1l) ->
             HExtlib.flatten_map 
               (fun (_,t,u2l) ->
                 HExtlib.flatten_map
                   (fun u1 ->
+                    if  eq_carr ~exact:true s t
+                     || eq_carr ~exact:true s tgt
+                     || eq_carr ~exact:true src t
+                    then [] else
                     List.map 
                       (fun u2 -> (s,[u1;uri;u2],t)) 
                       u2l)
@@ -76,25 +84,24 @@ exception UnableToCompose
 
 (* generate_composite (c2 (c1 s)) in the universe graph univ
  * both living in the same context and metasenv *)
-let generate_composite c1 c2 context metasenv univ arity last_lam_with_inn_arg =
-  let module RT = RefinementTool in
+let generate_composite' (c1,sat1) (c2,sat2) context metasenv univ arity =
+  let original_metasenv = metasenv in 
   let c1_ty,univ = CicTypeChecker.type_of_aux' metasenv context c1 univ in
   let c2_ty,univ = CicTypeChecker.type_of_aux' metasenv context c2 univ in
   let rec mk_implicits = function
     | 0 -> [] | n -> (Cic.Implicit None) :: mk_implicits (n-1)
   in
-  let rec mk_lambda_spline c namer = function
+  let rec mk_lambda_spine c namer = function
     | 0 -> c
     | n -> 
         Cic.Lambda 
           (namer n,
            (Cic.Implicit None), 
-           mk_lambda_spline (CicSubstitution.lift 1 c) namer (n-1))
+           mk_lambda_spine (CicSubstitution.lift 1 c) namer (n-1))
   in 
-  let count_saturations_needed t arity = 
+  let count_pis t arity = 
     let rec aux acc n = function
-      | Cic.Prod (name,src, ((Cic.Prod _) as t)) -> 
-          aux (acc@[name]) (n+1) t
+      | Cic.Prod (name,src,tgt) -> aux (acc@[name]) (n+1) tgt
       | _ -> n,acc
     in
     let len,names = aux [] 0 t in
@@ -104,36 +111,18 @@ let generate_composite c1 c2 context metasenv univ arity last_lam_with_inn_arg =
       names
   in
   let compose c1 nc1 c2 nc2 =
-    Cic.Lambda 
-      (Cic.Name "x", (Cic.Implicit None),
-          (Cic.Appl (  CicSubstitution.lift 1 c2 :: mk_implicits nc2 @ 
-            [ Cic.Appl (  CicSubstitution.lift 1 c1 :: mk_implicits nc1 @ 
-             [if last_lam_with_inn_arg then Cic.Rel 1 else Cic.Implicit None])
-            ])))
+   Cic.Appl (CicSubstitution.lift 1 c2 :: mk_implicits (nc2 - sat2 - 1) @
+     Cic.Appl (CicSubstitution.lift 1 c1 :: mk_implicits nc1 ) ::
+     mk_implicits sat2)
   in
-(*
-  let order_metasenv metasenv = 
-    let module OT = struct type t = int let compare = Pervasives.compare end in
-    let module S = HTopoSort.Make (OT) in
-    let dep i = 
-      let _,_,ty = List.find (fun (j,_,_) -> j=i) metasenv in
-      let metas = List.map fst (CicUtil.metas_of_term ty) in
-      HExtlib.list_uniq (List.sort Pervasives.compare metas)
-    in
-    let om =  
-      S.topological_sort (List.map (fun (i,_,_) -> i) metasenv) dep 
-    in
-    List.map (fun i -> List.find (fun (j,_,_) -> i=j) metasenv) om
-  in 
-*)
   let rec create_subst_from_metas_to_rels n = function 
     | [] -> []
     | (metano, ctx, ty)::tl -> 
-        (metano,(ctx,Cic.Rel (n+1),ty)) ::
+        (metano,(ctx,Cic.Rel n,ty)) ::
           create_subst_from_metas_to_rels (n-1) tl
   in
   let split_metasenv metasenv n =
-    List.partition (fun (_,ctx,_) -> List.length ctx > n) metasenv
+    List.partition (fun (_,ctx,_) -> List.length ctx >= n) metasenv
   in
   let purge_unused_lambdas metasenv t =
     let rec aux = function
@@ -146,7 +135,7 @@ let generate_composite c1 c2 context metasenv univ arity last_lam_with_inn_arg =
     in
     aux t
   in
-  let order_body_menv term body_metasenv =
+  let order_body_menv term body_metasenv c1_pis c2_pis =
     let rec purge_lambdas = function
       | Cic.Lambda (_,_,t) -> purge_lambdas t
       | t -> t
@@ -158,28 +147,52 @@ let generate_composite c1 c2 context metasenv univ arity last_lam_with_inn_arg =
           let metas = CicUtil.metas_of_term t in
           let metas = List.map fst metas in
           let metas = 
-            List.filter 
-              (fun i -> List.for_all (fun (j,_) -> j<>i) acc) 
-              metas 
-          in
+            List.filter (fun i -> List.for_all (fun (j,_) -> j<>i) acc) metas in
           let metas = List.map (fun i -> i,n) metas in
           metas @ acc, n+1)
         ([],0) l
     in
     let l_c2 = skip_appl (purge_lambdas term) in
-    let l_c1 = 
-      match HExtlib.list_last l_c2 with
-      | Cic.Appl l -> List.tl l
-      | _ -> assert false
+    let l_c2_b,l_c2_a =
+     try
+      HExtlib.split_nth (c2_pis - sat2 - 1) l_c2
+     with
+      Failure _ -> assert false in
+    let l_c1,l_c2_a =
+     match l_c2_a with
+        Cic.Appl (_::l_c1)::tl -> l_c1,tl
+      | _ -> assert false in
+    let meta_to_be_coerced =
+     try
+      match List.nth l_c1 (c1_pis - sat1 - 1) with
+       | Cic.Meta (i,_) -> i
+       | _ -> assert false
+     with
+      Failure _ -> assert false
     in
-    (* i should cut off the laet elem of l_c2 *)
-    let meta2no = fst (metas_that_saturate (l_c1 @ l_c2)) in
-    List.sort 
+    let meta2no = fst (metas_that_saturate (l_c2_b @ l_c1 @ l_c2_a)) in
+    let sorted =
+     List.sort 
       (fun (i,ctx1,ty1) (j,ctx1,ty1) -> 
           try List.assoc i meta2no -  List.assoc j meta2no 
-          with Not_found -> 
-            assert false) 
+          with Not_found -> assert false) 
       body_metasenv
+    in
+    let rec position_of n acc =
+     function
+        [] -> assert false
+      | (i,_,_)::_ when i = n -> acc
+      | _::tl -> position_of n (acc + 1) tl
+    in
+    debug_print
+     (lazy ("META_TO_BE_COERCED: " ^ string_of_int meta_to_be_coerced));
+    let position_of_meta_to_be_coerced =
+     position_of meta_to_be_coerced 0 sorted in
+     debug_print (lazy ("POSITION_OF_META_TO_BE_COERCED: " ^
+      string_of_int position_of_meta_to_be_coerced));
+     debug_print (lazy ("SATURATIONS: " ^
+      string_of_int (List.length sorted - position_of_meta_to_be_coerced - 1)));
+     sorted, List.length sorted - position_of_meta_to_be_coerced - 1
   in
   let namer l n = 
     let l = List.map (function Cic.Name s -> s | _ -> "A") l in
@@ -197,14 +210,16 @@ let generate_composite c1 c2 context metasenv univ arity last_lam_with_inn_arg =
   debug_print (lazy ("\nCOMPOSING"));
   debug_print (lazy (" c1= "^CicPp.ppterm c1 ^"  :  "^ CicPp.ppterm c1_ty));
   debug_print (lazy (" c2= "^CicPp.ppterm c2 ^"  :  "^ CicPp.ppterm c2_ty));
-  let saturations_for_c1, names_c1 = count_saturations_needed c1_ty 0 in 
-  let saturations_for_c2, names_c2 = count_saturations_needed c2_ty arity in
-  let c = compose c1 saturations_for_c1 c2 saturations_for_c2 in
-  let spline_len = saturations_for_c1 + saturations_for_c2 in
-  let c = mk_lambda_spline c (namer (names_c1 @ names_c2)) spline_len in
+  let c1_pis, names_c1 = count_pis c1_ty 0 in 
+  let c2_pis, names_c2 = count_pis c2_ty arity in
+  let c = compose c1 c1_pis c2 c2_pis in
+  let spine_len = c1_pis + c2_pis in
+  let c = mk_lambda_spine c (namer (names_c1 @ names_c2)) spine_len in
   debug_print (lazy ("COMPOSTA: " ^ CicPp.ppterm c));
-  let c, metasenv, univ = 
+  let old_insert_coercions = !CicRefine.insert_coercions in
+  let c, metasenv, univ, saturationsres =
     try
+      CicRefine.insert_coercions := false;
       let term, ty, metasenv, ugraph = 
         CicRefine.type_of_aux' metasenv context c univ
       in
@@ -212,32 +227,58 @@ let generate_composite c1 c2 context metasenv univ arity last_lam_with_inn_arg =
 (*       let metasenv = order_metasenv metasenv in *)
 (*       debug_print(lazy("ORDERED MENV: "^CicMetaSubst.ppmetasenv [] metasenv)); *)
       let body_metasenv, lambdas_metasenv = 
-        split_metasenv metasenv (spline_len + List.length context)
+        split_metasenv metasenv (spine_len + List.length context)
       in
       debug_print(lazy("B_MENV: "^CicMetaSubst.ppmetasenv [] body_metasenv));
       debug_print(lazy("L_MENV: "^CicMetaSubst.ppmetasenv [] lambdas_metasenv));
-      let body_metasenv = order_body_menv term body_metasenv in
+      let body_metasenv, saturationsres =
+       order_body_menv term body_metasenv c1_pis c2_pis
+      in
       debug_print(lazy("ORDERED_B_MENV: "^CicMetaSubst.ppmetasenv [] body_metasenv));
-      let subst = create_subst_from_metas_to_rels spline_len body_metasenv in
+      let subst = create_subst_from_metas_to_rels spine_len body_metasenv in
       debug_print (lazy("SUBST: "^CicMetaSubst.ppsubst body_metasenv subst));
       let term = CicMetaSubst.apply_subst subst term in
+      let metasenv = CicMetaSubst.apply_subst_metasenv subst metasenv in
       debug_print (lazy ("COMPOSED SUBSTITUTED: " ^ CicPp.ppterm term));
       let term, ty, metasenv, ugraph = 
         CicRefine.type_of_aux' metasenv context term ugraph
       in
       let body_metasenv, lambdas_metasenv = 
-        split_metasenv metasenv (spline_len + List.length context)
+        split_metasenv metasenv (spine_len + List.length context)
+      in
+      let lambdas_metasenv = 
+        List.filter 
+          (fun (i,_,_) -> 
+            List.for_all (fun (j,_,_) -> i <> j) original_metasenv)
+          lambdas_metasenv
       in
       let term = purge_unused_lambdas lambdas_metasenv term in
+      let metasenv = 
+        List.filter 
+          (fun (i,_,_) -> 
+            List.for_all 
+              (fun (j,_,_) ->
+                i <> j || List.exists (fun (j,_,_) -> j=i) original_metasenv) 
+              lambdas_metasenv) 
+          metasenv 
+      in
+      debug_print (lazy ("####################"));
       debug_print (lazy ("COMPOSED: " ^ CicPp.ppterm term));
-      debug_print(lazy("MENV: "^CicMetaSubst.ppmetasenv [] metasenv));
-      term, metasenv, ugraph
+      debug_print (lazy ("SATURATIONS: " ^ string_of_int saturationsres));
+      debug_print (lazy ("MENV: "^CicMetaSubst.ppmetasenv [] metasenv));
+      debug_print (lazy ("####################"));
+      CicRefine.insert_coercions := old_insert_coercions;
+      term, metasenv, ugraph, saturationsres
     with
     | CicRefine.RefineFailure s 
     | CicRefine.Uncertain s -> debug_print s; 
+        CicRefine.insert_coercions := old_insert_coercions;
         raise UnableToCompose
-  in  
-  c, metasenv, univ 
+    | exn ->
+        CicRefine.insert_coercions := old_insert_coercions;
+        raise exn
+  in
+  c, metasenv, univ, saturationsres
 ;;
 
 let build_obj c univ arity =
@@ -264,7 +305,7 @@ let filter_duplicates l coercions =
           CoercDb.eq_carr s src && 
           CoercDb.eq_carr t tgt &&
           try 
-            List.for_all2 (fun u1 u2 -> UriManager.eq u1 u2) l1 l2
+            List.for_all2 (fun (u1,_) (u2,_) -> UriManager.eq u1 u2) l1 l2
           with
           | Invalid_argument "List.for_all2" -> false)
         coercions))
@@ -301,8 +342,8 @@ let number_if_already_defined buri name l =
     if List.exists (UriManager.eq uri) l then retry ()
     else
       try
-        let _  = Http_getter.resolve' ~writable:true uri in
-        if Http_getter.exists' uri then retry () else uri
+        let _  = Http_getter.resolve' ~local:true ~writable:true uri in
+        if Http_getter.exists' ~local:true uri then retry () else uri
       with 
       | Http_getter_types.Key_not_found _ -> uri
       | Http_getter_types.Unresolvable_URI _ -> assert false
@@ -313,10 +354,10 @@ let number_if_already_defined buri name l =
 (* given a new coercion uri from src to tgt returns 
  * a list of (new coercion uri, coercion obj, universe graph) 
  *)
-let close_coercion_graph _ src tgt uri baseuri =
+let close_coercion_graph src tgt uri saturations baseuri =
   (* check if the coercion already exists *)
   let coercions = CoercDb.to_list () in
-  let todo_list = get_closure_coercions src tgt uri coercions in
+  let todo_list = get_closure_coercions src tgt (uri,saturations) coercions in
   let todo_list = filter_duplicates todo_list coercions in
   try
     let new_coercions = 
@@ -325,42 +366,44 @@ let close_coercion_graph _ src tgt uri baseuri =
           try 
             (match l with
             | [] -> assert false 
-            | he :: tl ->
+            | (he,saturations1) :: tl ->
                 let arity = match tgt with CoercDb.Fun n -> n | _ -> 0 in
                 let first_step = 
                   Cic.Constant ("", 
                     Some (CoercDb.term_of_carr (CoercDb.Uri he)),
-                    Cic.Sort Cic.Prop, [], obj_attrs arity)
+                    Cic.Sort Cic.Prop, [], obj_attrs arity), saturations1
                 in
                 let o,_ = 
-                  List.fold_left (fun (o,univ) coer ->
+                  List.fold_left (fun (o,univ) (coer,saturations) ->
                     match o with 
-                    | Cic.Constant (_,Some c,_,[],_) ->
-                        let t, menv, univ = 
-                          generate_composite c 
-                            (CoercDb.term_of_carr (CoercDb.Uri coer)
-                            [] [] univ arity true
+                    | Cic.Constant (_,Some u,_,[],_),saturations1 ->
+                        let t, menv, univ, saturationsres = 
+                          generate_composite' (u,saturations1) 
+                            (CoercDb.term_of_carr (CoercDb.Uri coer),
+                             saturations) [] [] univ arity
                         in
-                        assert (menv = []);
-                        build_obj t univ arity
+                        if (menv = []) then
+                          HLog.warn "MENV non empty after composing coercions";
+                        let o,univ = build_obj t univ arity in
+                         (o,saturationsres),univ
                     | _ -> assert false 
                   ) (first_step, CicUniv.empty_ugraph) tl
                 in
                 let name_src = CoercDb.name_of_carr src in
                 let name_tgt = CoercDb.name_of_carr tgt in
-                let by = List.map UriManager.name_of_uri l in
+                let by = List.map (fun u,_ -> UriManager.name_of_uri u) l in
                 let name = mangle name_tgt name_src by in
                 let c_uri = 
                   number_if_already_defined baseuri name 
-                    (List.map (fun (_,_,u,_) -> u) acc) 
+                    (List.map (fun (_,_,u,_,_) -> u) acc) 
                 in
-                let named_obj = 
+                let named_obj,saturations = 
                   match o with
-                  | Cic.Constant (_,bo,ty,vl,attrs) ->
-                      Cic.Constant (name,bo,ty,vl,attrs)
+                  | Cic.Constant (_,bo,ty,vl,attrs),saturations ->
+                      Cic.Constant (name,bo,ty,vl,attrs),saturations
                   | _ -> assert false 
                 in
-                  (src,tgt,c_uri,named_obj))::acc
+                  (src,tgt,c_uri,saturations,named_obj))::acc
           with UnableToCompose -> acc
       ) [] todo_list
     in
@@ -369,3 +412,12 @@ let close_coercion_graph _ src tgt uri baseuri =
 ;;
 
 CicCoercion.set_close_coercion_graph close_coercion_graph;;
+
+(* generate_composite (c2 (c1 s)) in the universe graph univ
+ * both living in the same context and metasenv *)
+let generate_composite c1 c2 context metasenv univ arity =
+ let a,b,c,_ =
+  generate_composite' (c1,0) (c2,0) context metasenv univ arity
+ in
+  a,b,c
+;;