]> matita.cs.unibo.it Git - helm.git/blobdiff - helm/software/components/ng_refiner/nCicRefiner.ml
arithmetics for λδ
[helm.git] / helm / software / components / ng_refiner / nCicRefiner.ml
index 656955808d5abaa231668c99ede556c9325627ac..705d80cd451611015abcd374639959408ac91e11 100644 (file)
@@ -18,16 +18,37 @@ exception AssertFailure of string Lazy.t;;
 module C = NCic
 module Ref = NReference
 
+let debug = ref false;;
 let indent = ref "";;
-let inside c = indent := !indent ^ String.make 1 c;;
-let outside () = indent := String.sub !indent 0 (String.length !indent -1);;
-
-
-let pp s = 
+let times = ref [];;
+let pp s =
+ if !debug then
   prerr_endline (Printf.sprintf "%-20s" !indent ^ " " ^ Lazy.force s)
-;;  
+;;
+let inside c =
+ if !debug then
+  begin
+   let time1 = Unix.gettimeofday () in
+   indent := !indent ^ String.make 1 c;
+   times := time1 :: !times;
+   prerr_endline ("{{{" ^ !indent ^ " ")
+  end
+;;
+let outside ok =
+ if !debug then
+  begin
+   let time2 = Unix.gettimeofday () in
+   let time1 =
+    match !times with time1::tl -> times := tl; time1 | [] -> assert false in
+   prerr_endline ("}}} " ^ string_of_float (time2 -. time1));
+   if not ok then prerr_endline "exception raised!";
+   try
+    indent := String.sub !indent 0 (String.length !indent -1)
+   with
+    Invalid_argument _ -> indent := "??"; ()
+ end
+;;
 
-let pp _ = ();;
 
 let wrap_exc msg = function
   | NCicUnification.Uncertain _ -> Uncertain msg
@@ -36,12 +57,34 @@ let wrap_exc msg = function
   | e -> raise e
 ;;
 
-let exp_implicit ~localise metasenv context expty t =
- let foo x = match expty with Some t -> `WithType t | None -> x in
+let exp_implicit rdb ~localise metasenv subst context with_type t =
  function      
-  | `Closed -> NCicMetaSubst.mk_meta metasenv [] (foo `Term)
-  | `Type -> NCicMetaSubst.mk_meta metasenv context (foo `Type)
-  | `Term -> NCicMetaSubst.mk_meta metasenv context (foo `Term)
+  | `Closed ->
+      let metasenv,subst,expty =
+       match with_type with
+          None -> metasenv,subst,None
+        | Some typ ->
+           let (metasenv,subst),typ =
+            try
+             NCicMetaSubst.delift
+              ~unify:(fun m s c t1 t2 ->
+                try Some (NCicUnification.unify rdb m s c t1 t2)
+                with NCicUnification.UnificationFailure _ | NCicUnification.Uncertain _ -> None)
+              metasenv subst context 0 (0,NCic.Irl 0) typ
+            with
+              NCicMetaSubst.MetaSubstFailure _
+            | NCicMetaSubst.Uncertain _ ->
+               raise (RefineFailure (lazy (localise t,"Trying to create a closed meta with a non closed type: " ^ NCicPp.ppterm ~metasenv ~subst ~context typ)))
+
+           in
+            metasenv,subst,Some typ
+      in
+       NCicMetaSubst.mk_meta metasenv [] ?with_type:expty `IsTerm,subst
+  | `Type -> NCicMetaSubst.mk_meta metasenv context ?with_type `IsType,subst
+  | `Term -> NCicMetaSubst.mk_meta metasenv context ?with_type `IsTerm,subst
+  | `Tagged s ->
+      NCicMetaSubst.mk_meta 
+        ~attrs:[`Name s] metasenv context ?with_type `IsTerm,subst
   | `Vector ->
       raise (RefineFailure (lazy (localise t, "A vector of implicit terms " ^
        "can only be used in argument position")))
@@ -63,7 +106,7 @@ let check_allowed_sort_elimination rdb localise r orig =
      match arity1 with
      | C.Prod (name,so1,de1) (* , t ==?== C.Prod _ *) ->
         let metasenv, _, meta, _ = 
-          NCicMetaSubst.mk_meta metasenv ((name,C.Decl so1)::context) `Type 
+          NCicMetaSubst.mk_meta metasenv ((name,C.Decl so1)::context) `IsType 
         in
         let metasenv, subst = 
           try NCicUnification.unify rdb metasenv subst context 
@@ -76,7 +119,7 @@ let check_allowed_sort_elimination rdb localise r orig =
         aux metasenv subst ((name, C.Decl so1)::context)
          (mkapp (NCicSubstitution.lift 1 ind) (C.Rel 1)) de1 meta
      | C.Sort _ (* , t ==?== C.Prod _ *) ->
-        let metasenv, _, meta, _ = NCicMetaSubst.mk_meta metasenv [] `Type in
+        let metasenv, _, meta, _ = NCicMetaSubst.mk_meta metasenv [] `IsSort in
         let metasenv, subst = 
           try NCicUnification.unify rdb metasenv subst context 
                 arity2 (C.Prod ("_", ind, meta)) 
@@ -91,7 +134,7 @@ let check_allowed_sort_elimination rdb localise r orig =
         with exc -> raise (wrap_exc (lazy (localise orig, 
           "Sort elimination not allowed ")) exc))
      | _ -> assert false
-     (*D*)in outside(); rc with exc -> outside (); raise exc
+     (*D*)in outside true; rc with exc -> outside false; raise exc
    in
     aux
 ;;
@@ -116,7 +159,7 @@ let rec typeof rdb
              let metasenv, subst =
     (*D*)inside 'U'; try let rc = 
                NCicUnification.unify rdb metasenv subst context infty expty 
-    (*D*)in outside(); rc with exc -> outside (); raise exc
+    (*D*)in outside true; rc with exc -> outside false; raise exc
              in
              metasenv, subst, t, expty
            with 
@@ -125,7 +168,7 @@ let rec typeof rdb
              try_coercions rdb ~localise 
                metasenv subst context t orig infty expty true exc)
     | None -> metasenv, subst, t, infty
-    (*D*)in outside(); rc with exc -> outside (); raise exc
+    (*D*)in outside true; rc with exc -> outside false; raise exc
   in
   let rec typeof_aux metasenv subst context expty = 
     fun t as orig -> 
@@ -152,24 +195,19 @@ let rec typeof rdb
               raise (RefineFailure (lazy (localise t, Lazy.force msg)))
          | NCicEnvironment.AssertFailure msg -> raise (AssertFailure msg))
     | C.Implicit infos -> 
-         let metasenv,_,t,ty =
-           exp_implicit ~localise metasenv context expty t infos
+         let (metasenv,_,t,ty),subst =
+           exp_implicit rdb ~localise metasenv subst context expty t infos
          in
+         if expty = None then
+          typeof_aux metasenv subst context None t
+         else
           metasenv, subst, t, ty 
     | C.Meta (n,l) as t -> 
-       let ty =
+       let metasenv, ty =
         try
-         let _,_,_,ty = NCicUtils.lookup_subst n subst in ty 
-        with NCicUtils.Subst_not_found _ -> try
-         let _,_,ty = NCicUtils.lookup_meta n metasenv in 
-         match ty with C.Implicit _ -> 
-                 prerr_endline (string_of_int n);
-                 prerr_endline (NCicPp.ppmetasenv ~subst metasenv);
-                 prerr_endline (NCicPp.ppsubst ~metasenv subst);
-                 assert false | _ -> ty 
-        with NCicUtils.Meta_not_found _ ->
-         raise (AssertFailure (lazy (Printf.sprintf
-          "%s not found" (NCicPp.ppterm ~subst ~metasenv ~context t))))
+         let _,_,_,ty = NCicUtils.lookup_subst n subst in metasenv, ty 
+        with NCicUtils.Subst_not_found _ ->
+         NCicMetaSubst.extend_meta metasenv n
        in
        metasenv, subst, t, NCicSubstitution.subst_meta l ty
     | C.Const _ -> 
@@ -198,23 +236,34 @@ let rec typeof rdb
              | C.Prod (_,s,t) -> Some s, Some t, false
              | _ -> None, None, true 
        in
-       let metasenv, subst, s, ty_s = 
-         typeof_aux metasenv subst context None s in
-       let metasenv, subst, s, _ = 
-         force_to_sort rdb 
-           metasenv subst context s orig_s localise ty_s in
-       let (metasenv,subst), exp_ty_t = 
+       let (metasenv,subst), s = 
          match exp_s with 
-         | Some exp_s -> 
-             (try 
-               pp(lazy("Force source to: "^NCicPp.ppterm ~metasenv ~subst
-                  ~context exp_s));
-               NCicUnification.unify rdb metasenv subst context s exp_s,exp_ty_t
-             with exc -> raise (wrap_exc (lazy (localise orig_s, Printf.sprintf
-               "Source type %s was expected to be %s" (NCicPp.ppterm ~metasenv
-               ~subst ~context s) (NCicPp.ppterm ~metasenv ~subst ~context
-               exp_s))) exc))
-         | None -> (metasenv, subst), None
+         | Some exp_s ->
+             (* optimized case: implicit and implicitly typed meta
+              * the optimization prevents proliferation of metas  *)
+             (match s with
+              | C.Implicit _ -> (metasenv,subst),exp_s
+              | _ ->
+                  let metasenv, subst, s = match s with 
+                    | C.Meta (n,_) 
+                        when (try (match NCicUtils.lookup_meta n metasenv with
+                              | _,_,C.Implicit (`Typeof _) -> true
+                              | _ -> false)
+                        with 
+                          | _ -> false) -> metasenv, subst, s
+                    | _ ->  check_type rdb ~localise metasenv subst context s in
+                  (try 
+                    pp(lazy("Force source to: "^NCicPp.ppterm ~metasenv ~subst
+                       ~context exp_s));
+                    NCicUnification.unify ~test_eq_only:true rdb metasenv subst context s exp_s,s
+                  with exc -> raise (wrap_exc (lazy (localise orig_s, Printf.sprintf
+                    "Source type %s was expected to be %s" (NCicPp.ppterm ~metasenv
+                    ~subst ~context s) (NCicPp.ppterm ~metasenv ~subst ~context
+                    exp_s))) exc)))
+         | None -> 
+             let metasenv, subst, s = 
+               check_type rdb ~localise metasenv subst context s in
+             (metasenv, subst), s
        in
        let context_for_t = (n,C.Decl s) :: context in
        let metasenv, subst, t, ty_t = 
@@ -225,12 +274,9 @@ let rec typeof rdb
            (C.Lambda(n,s,t)) (C.Prod (n,s,ty_t)) expty
        else 
          metasenv, subst, C.Lambda(n,s,t), C.Prod (n,s,ty_t)
-    | C.LetIn (n,(ty as orig_ty),t,bo) ->
-       let metasenv, subst, ty, ty_ty = 
-         typeof_aux metasenv subst context None ty in
-       let metasenv, subst, ty, _ = 
-         force_to_sort rdb
-           metasenv subst context ty orig_ty localise ty_ty in
+    | C.LetIn (n,ty,t,bo) ->
+       let metasenv, subst, ty = 
+         check_type rdb ~localise metasenv subst context ty in
        let metasenv, subst, t, _ = 
          typeof_aux metasenv subst context (Some ty) t in
        let context1 = (n, C.Def (t,ty)) :: context in
@@ -240,7 +286,8 @@ let rec typeof rdb
          | Some x -> 
              let m, s, x = 
                NCicUnification.delift_type_wrt_terms 
-                 rdb metasenv subst context x [t]
+                rdb metasenv subst context1 (NCicSubstitution.lift 1 x)
+                [NCicSubstitution.lift 1 t]
              in
                m, s, Some x
        in
@@ -254,7 +301,7 @@ let rec typeof rdb
        let hbr t = 
          if upto > 0 then NCicReduction.head_beta_reduce ~upto t else t 
        in
-       let refine_appl () =
+       let refine_appl metasenv subst args =
          let metasenv, subst, he, ty_he = 
             typeof_aux metasenv subst context None he in
          let metasenv, subst, t, ty = 
@@ -269,8 +316,53 @@ let rec typeof rdb
            let metasenv, subst, he, ty_he = 
               typeof_aux metasenv subst context expty he in
            metasenv, subst, hbr he, ty_he
-         with Uncertain _ | RefineFailure _ -> refine_appl ()
-       else refine_appl ()
+         with Uncertain _ | RefineFailure _ -> refine_appl metasenv subst args
+       else
+        (* CSC: whd can be useful on he too... *)
+        (match he with
+            C.Const (Ref.Ref (uri1,Ref.Con _)) -> (
+             match
+              HExtlib.map_option (NCicReduction.whd ~subst context) expty
+             with
+                Some (C.Appl(C.Const(Ref.Ref (uri2,Ref.Ind _) as ref)::expargs))
+                when NUri.eq uri1 uri2 ->
+                (try
+                 let _,leftno,_,_,_ = NCicEnvironment.get_checked_indtys ref in
+                 let leftexpargs,_ = HExtlib.split_nth leftno expargs in
+                  let rec instantiate metasenv subst revargs args =
+                   function
+                     [] ->
+                      (* some checks are re-done here, but it would be complex
+                         to avoid them (e.g. we did not check yet that the
+                         constructor is a constructor for that inductive type)*)
+                      refine_appl metasenv subst ((List.rev revargs)@args)
+                   | (exparg::expargs) as allexpargs ->
+                      match args with
+                         [] -> raise (Failure "Not enough args")
+                       | (C.Implicit `Vector::args) as allargs ->
+                          (try
+                            instantiate metasenv subst revargs args allexpargs
+                           with
+                              Sys.Break as exn -> raise exn
+                            | _ ->
+                             instantiate metasenv subst revargs
+                              (C.Implicit `Term :: allargs) allexpargs)
+                       | arg::args ->
+                          let metasenv,subst,arg,_ =
+                           typeof_aux metasenv subst context None arg in
+                          let metasenv,subst =
+                           NCicUnification.unify rdb metasenv subst context
+                            arg exparg
+                          in
+                           instantiate metasenv subst(arg::revargs) args expargs
+                  in
+                   instantiate metasenv subst [] args leftexpargs 
+                with
+                 | Sys.Break as exn -> raise exn
+                 | _ ->
+                    refine_appl metasenv subst args (* to try coercions *))
+              | _ -> refine_appl metasenv subst args)
+          | _ -> refine_appl metasenv subst args)
    | C.Appl _ -> raise (AssertFailure (lazy "Appl of length < 2"))
    | C.Match (Ref.Ref (_,Ref.Ind (_,tyno,_)) as r,
           outtype,(term as orig_term),pl) as orig ->
@@ -299,16 +391,11 @@ let rec typeof rdb
       (* next lines are to do a subst-expansion of the outtype, so
          that when it becomes a beta-abstraction, the beta-redex is
          fired during substitution *)
-      (*CSC: this is instantiate! should we move it from tactics to the
-             refiner? I think so! *)
-      let metasenv,metanoouttype,newouttype,metaoutsort =
-       NCicMetaSubst.mk_meta metasenv context `Term in
-      let metasenv,subst =
-       NCicUnification.unify rdb metasenv subst context outsort metaoutsort in
-      let metasenv =
-       List.filter (function (j,_) -> j <> metanoouttype) metasenv in
+      let _,fresh_metanoouttype,newouttype,_ =
+       NCicMetaSubst.mk_meta metasenv context `IsTerm in
       let subst =
-       (metanoouttype,(Some "outtype",context,outtype,metaoutsort))::subst in
+       (fresh_metanoouttype,([`Name "outtype"],context,outtype,outsort))
+         ::subst in
       let outtype = newouttype in
 
       (* let's control if the sort elimination is allowed: [(I q1 ... qr)|B] *)
@@ -362,19 +449,37 @@ let rec typeof rdb
     pp (lazy (NCicPp.ppterm ~metasenv ~subst ~context t ^ " ::inf:: "^
          NCicPp.ppterm ~metasenv ~subst ~context infty ));
       force_ty true true metasenv subst context orig t infty expty
-    (*D*)in outside(); rc with exc -> outside (); raise exc
+    (*D*)in outside true; rc with exc -> outside false; raise exc
   in
     typeof_aux metasenv subst context expty term
 
+and check_type rdb ~localise metasenv subst context (ty as orig_ty) =
+  let metasenv, subst, ty, sort = 
+    typeof rdb ~localise metasenv subst context ty None
+  in
+  let metasenv, subst, ty, _ = 
+    force_to_sort rdb metasenv subst context ty orig_ty localise sort
+  in
+    metasenv, subst, ty
+
 and try_coercions rdb 
   ~localise 
   metasenv subst context t orig_t infty expty perform_unification exc 
 =
+  let cs_subst = NCicSubstitution.subst ~avoid_beta_redexes:true in
+  try 
+    pp (lazy "WWW: trying coercions -- preliminary unification");
+    let metasenv, subst = 
+      NCicUnification.unify rdb metasenv subst context infty expty
+    in metasenv, subst, t, expty
+  with
+  | exn ->
   (* we try with a coercion *)
   let rec first exc = function
-  | [] ->         
+  | [] ->   
+      pp (lazy "WWW: no more coercions!");      
       raise (wrap_exc (lazy (localise orig_t, Printf.sprintf
-        "The term %s has type %s but is here used with type %s"
+        "The term\n%s\nhas type\n%s\nbut is here used with type\n%s"
         (NCicPp.ppterm ~metasenv ~subst ~context t)
         (NCicPp.ppterm ~metasenv ~subst ~context infty)
         (NCicPp.ppterm ~metasenv ~subst ~context expty))) exc)
@@ -407,7 +512,332 @@ and try_coercions rdb
       with 
       | NCicUnification.UnificationFailure _ -> first exc tl
       | NCicUnification.Uncertain _ as exc -> first exc tl
-  in 
+  in
+  
+  try 
+    pp (lazy "WWW: trying coercions -- inner check");
+    match infty, expty, t with
+    | _,_, NCic.Match (Ref.Ref (_,Ref.Ind (_,tyno,leftno)) as r,outty,m,pl) ->
+        (*{{{*) pp (lazy "CASE");
+        (* {{{ helper functions *)
+        let get_cl_and_left_p refit outty =
+          match refit with
+          | NReference.Ref (uri, NReference.Ind (_,tyno,leftno)) ->
+             let _, leftno, itl, _, _ = NCicEnvironment.get_checked_indtys r in
+             let _, _, ty, cl = List.nth itl tyno in
+             let constructorsno = List.length cl in
+              let count_pis t =
+                let rec aux ctx t = 
+                  match NCicReduction.whd ~subst ~delta:max_int ctx t with
+                  | NCic.Prod (name,src,tgt) ->
+                      let ctx = (name, (NCic.Decl src)) :: ctx in
+                      1 + aux ctx tgt
+                  | _ -> 0
+                in
+                  aux [] t
+              in
+              let rec skip_lambda_delifting t n =
+                match t,n with
+                | _,0 -> t
+                | NCic.Lambda (_,_,t),n ->
+                    pp (lazy ("WWW: failing term? "^ NCicPp.ppterm ~subst:[] ~metasenv:[] ~context:[] t)); 
+                    skip_lambda_delifting
+                      (* substitute dangling indices with a dummy *)
+                      (NCicSubstitution.subst (NCic.Sort NCic.Prop) t) (n - 1)
+                | _ -> assert false
+              in
+              let get_l_r_p n = function
+                | NCic.Lambda (_,NCic.Const (Ref.Ref (_,Ref.Ind (_,_,_))),_) -> [],[]
+                | NCic.Lambda (_,NCic.Appl 
+                    (NCic.Const (Ref.Ref (_,Ref.Ind (_,_,_))) :: args),_) ->
+                    HExtlib.split_nth n args
+                | _ -> assert false
+              in
+              let pis = count_pis ty in
+              let rno = pis - leftno in
+              let t = skip_lambda_delifting outty rno in
+              let left_p, _ = get_l_r_p leftno t in
+              let instantiate_with_left cl =
+                List.map 
+                  (fun ty -> 
+                    List.fold_left 
+                      (fun t p -> match t with
+                        | NCic.Prod (_,_,t) ->
+                            cs_subst p t
+                        | _-> assert false)
+                      ty left_p) 
+                  cl 
+              in
+              let cl = instantiate_with_left (List.map (fun (_,_,x) -> x) cl) in
+              cl, left_p, leftno, rno
+          | _ -> raise exn
+        in
+        let rec keep_lambdas_and_put_expty ctx t bo right_p matched n =
+          match t,n with
+          | _,0 ->
+            let rec mkr n = function 
+              | [] -> [] | _::tl -> NCic.Rel n :: mkr (n+1) tl
+            in
+            pp (lazy ("input replace: " ^ NCicPp.ppterm ~context:ctx ~metasenv:[] ~subst:[] bo));
+            let bo =
+            NCicRefineUtil.replace_lifting
+              ~equality:(fun _ -> NCicRefineUtil.alpha_equivalence)
+              ~context:ctx
+              ~what:(matched::right_p)
+              ~with_what:(NCic.Rel 1::List.rev (mkr 2 right_p))
+              ~where:bo
+            in
+            pp (lazy ("output replace: " ^ NCicPp.ppterm ~context:ctx ~metasenv:[] ~subst:[] bo));
+            bo
+          | NCic.Lambda (name, src, tgt),_ ->
+              NCic.Lambda (name, src,
+               keep_lambdas_and_put_expty 
+                ((name, NCic.Decl src)::ctx) tgt (NCicSubstitution.lift 1 bo)
+                (List.map (NCicSubstitution.lift 1) right_p) (NCicSubstitution.lift 1 matched) (n-1))
+          | _ -> assert false
+        in
+        (*let eq = NCic.Const (NReference.reference_of_string ("cic:/matita/ng/Plogic/equality/eq.ind(1,0,2)")) in
+        let eq_refl = NCic.Const (NReference.reference_of_string ("cic:/matita/ng/Plogic/equality/eq.con(0,1,2)")) in*)
+        let eq = NCic.Const (NReference.reference_of_string ("cic:/matita/ng/Plogic/jmeq/jmeq.ind(1,0,2)")) in
+        let eq_refl = NCic.Const (NReference.reference_of_string ("cic:/matita/ng/Plogic/jmeq/jmeq.con(0,1,2)")) in
+        let add_params 
+          metasenv subst context cty outty mty m i 
+        =
+          let rec aux context outty par j mty m = function
+            | NCic.Prod (name, src, tgt) ->
+                let t,k = 
+                  aux 
+                    ((name, NCic.Decl src) :: context)
+                    (NCicSubstitution.lift 1 outty) (NCic.Rel j::par) (j+1) 
+                    (NCicSubstitution.lift 1 mty) (NCicSubstitution.lift 1 m) tgt
+                in
+                NCic.Prod (name, src, t), k
+            | NCic.Const (Ref.Ref (_,Ref.Ind (_,_,leftno)) as r) ->
+                let k = 
+                  let k = NCic.Const(Ref.mk_constructor i r) in
+                  NCicUntrusted.mk_appl k par
+                in
+                (* the type has no parameters, so kty = mty
+                let kty = 
+                  try NCicTypechecker.typeof ~subst ~metasenv context k
+                  with NCicTypeChecker.TypeCheckerFailure _ -> assert false
+                in *)
+                NCic.Prod ("p", 
+                 NCic.Appl [eq; mty; m; mty; k],
+                 (NCicSubstitution.lift 1
+                  (NCicReduction.head_beta_reduce ~delta:max_int 
+                   (NCicUntrusted.mk_appl outty [k])))),[mty,m,mty,k]
+            | NCic.Appl (NCic.Const (Ref.Ref (_,Ref.Ind (_,_,leftno)) as r)::pl) ->
+                let left_p,right_p = HExtlib.split_nth leftno pl in
+                let has_rights = right_p <> [] in
+                let k = 
+                  let k = NCic.Const(Ref.mk_constructor i r) in
+                  NCicUntrusted.mk_appl k (left_p@par)
+                in
+                let right_p = 
+                  try match 
+                    NCicTypeChecker.typeof ~subst ~metasenv context k
+                  with
+                  | NCic.Appl (NCic.Const (Ref.Ref (_,Ref.Ind (_,_,_)))::args) ->
+                      snd (HExtlib.split_nth leftno args)
+                  | _ -> assert false
+                  with NCicTypeChecker.TypeCheckerFailure _-> assert false
+                in
+                let orig_right_p =
+                  match mty with
+                  | NCic.Appl (NCic.Const (Ref.Ref (_,Ref.Ind (_,_,_)))::args) ->
+                      snd (HExtlib.split_nth leftno args)
+                  | _ -> assert false
+                in
+                List.fold_right2 
+                  (fun x y (tacc,pacc) ->
+                    let xty = 
+                   try NCicTypeChecker.typeof ~subst ~metasenv context x
+                      with NCicTypeChecker.TypeCheckerFailure _ -> assert false
+                    in
+                    let yty = 
+                   try NCicTypeChecker.typeof ~subst ~metasenv context y
+                      with NCicTypeChecker.TypeCheckerFailure _ -> assert false
+                    in
+                    NCic.Prod ("_", NCic.Appl [eq;xty;x;yty;y],
+                     NCicSubstitution.lift 1 tacc), (xty,x,yty,y)::pacc) 
+                  (orig_right_p @ [m]) (right_p @ [k]) 
+                  (NCicReduction.head_beta_reduce ~delta:max_int
+                      (NCicUntrusted.mk_appl outty (right_p@[k])), [])  
+
+  (*              if has_rights then
+                  NCicReduction.head_beta_reduce ~delta:max_int
+                    (NCic.Appl (outty ::right_p @ [k])),k
+                else
+                  NCic.Prod ("p", 
+                   NCic.Appl [eq; mty; m; k],
+                   (NCicSubstitution.lift 1
+                    (NCicReduction.head_beta_reduce ~delta:max_int 
+                     (NCic.Appl (outty ::right_p @ [k]))))),k*)
+            | _ -> assert false
+          in
+            aux context outty [] 1 mty m cty
+        in
+        let add_lambda_for_refl_m pbo eqs cty =
+          (* k lives in the right context *)
+          (* if rno <> 0 then pbo else *)
+          let rec aux = function 
+            | NCic.Lambda (name,src,bo), NCic.Prod (_,_,ty) ->
+                NCic.Lambda (name,src,
+                (aux (bo,ty)))
+            | t,_ -> snd (List.fold_right
+                (fun (xty,x,yty,y) (n,acc) -> n+1, NCic.Lambda ("p" ^ string_of_int n,
+                  NCic.Appl [eq; xty; x; yty; y],
+                  NCicSubstitution.lift 1 acc)) eqs (1,t))
+                (*NCic.Lambda ("p",
+                  NCic.Appl [eq; mty; m; k],NCicSubstitution.lift 1 t)*)
+          in
+          aux (pbo,cty)
+        in
+        let add_pi_for_refl_m context new_outty mty m lno rno =
+          (*if rno <> 0 then new_outty else*)
+            let rec aux context mty m = function
+              | NCic.Lambda (name, src, tgt) ->
+                  let context = (name, NCic.Decl src)::context in
+                  NCic.Lambda (name, src, aux context (NCicSubstitution.lift 1 mty) (NCicSubstitution.lift 1 m) tgt)
+              | t -> 
+                  let lhs =
+                    match mty with
+                    | NCic.Appl (_::params) -> (snd (HExtlib.split_nth lno params))@[m]
+                    | _ -> [m]
+                  in
+                  let rhs = 
+                    List.map (fun x -> NCic.Rel x) 
+                      (List.rev (HExtlib.list_seq 1 (rno+2))) in
+                  List.fold_right2
+                    (fun x y acc ->
+                      let xty = 
+                   try NCicTypeChecker.typeof ~subst ~metasenv context x
+                        with NCicTypeChecker.TypeCheckerFailure _ -> assert false
+                      in
+                      let yty = 
+                   try NCicTypeChecker.typeof ~subst ~metasenv context y
+                        with NCicTypeChecker.TypeCheckerFailure _ -> assert false
+                      in
+                      NCic.Prod
+                      ("_", NCic.Appl [eq;xty;x;yty;y],
+                       (NCicSubstitution.lift 1 acc)))
+                    lhs rhs t
+  (*                NCic.Prod 
+                    ("_", NCic.Appl [eq;mty;m;NCic.Rel 1],
+                    NCicSubstitution.lift 1 t)*)
+            in
+              aux context mty m new_outty
+        in (* }}} end helper functions *)
+        (* constructors types with left params already instantiated *)
+        let outty = NCicUntrusted.apply_subst subst context outty in
+        let cl, left_p, leftno,rno = 
+          get_cl_and_left_p r outty
+        in
+        let right_p, mty = 
+          try
+            match 
+              NCicTypeChecker.typeof ~subst ~metasenv context m
+            with
+            | (NCic.Const (Ref.Ref (_,Ref.Ind (_,_,_))) | NCic.Meta _) as mty -> [], mty
+            | NCic.Appl ((NCic.Const (Ref.Ref (_,Ref.Ind (_,_,_)))|NCic.Meta _)::args) as mty ->
+                snd (HExtlib.split_nth leftno args), mty
+            | _ -> assert false
+          with NCicTypeChecker.TypeCheckerFailure _ -> 
+             raise (AssertFailure(lazy "already ill-typed matched term"))
+        in
+        let mty = NCicUntrusted.apply_subst subst context mty in
+        let outty = NCicUntrusted.apply_subst subst context outty in
+        let expty = NCicUntrusted.apply_subst subst context expty in
+        let new_outty =
+         keep_lambdas_and_put_expty context outty expty right_p m (rno+1)
+        in
+        pp 
+          (lazy ("CASE: new_outty: " ^ NCicPp.ppterm 
+           ~context ~metasenv ~subst new_outty));
+        let _,pl,subst,metasenv = 
+          List.fold_right2
+            (fun cty pbo (i, acc, s, menv) -> 
+               pp (lazy ("begin iteration"));
+              (* Pi k_par, (Pi H:m=(K_i left_par k_par)), 
+               *   (new_)outty right_par (K_i left_par k_par) *)
+               let infty_pbo, _ = 
+                 add_params menv s context cty outty mty m i in
+               pp 
+                (lazy ("CASE: infty_pbo: "^ NCicPp.ppterm
+                 ~context ~metasenv ~subst infty_pbo));
+               let expty_pbo, eqs = (* k is (K_i left_par k_par) *)
+                 add_params menv s context cty new_outty mty m i in
+               pp 
+                (lazy ("CASE: expty_pbo: "^ NCicPp.ppterm
+                 ~context ~metasenv ~subst expty_pbo));
+               let pbo = add_lambda_for_refl_m pbo eqs cty in
+               pp 
+                 (lazy ("CASE: pbo: " ^ NCicPp.ppterm
+                 ~context ~metasenv ~subst pbo));
+               let metasenv, subst, pbo, _ =
+                 try_coercions rdb ~localise menv s context pbo 
+                orig_t (*??*) infty_pbo expty_pbo perform_unification exc
+               in
+               pp 
+                 (lazy ("CASE: pbo2: " ^ NCicPp.ppterm 
+                 ~context ~metasenv ~subst pbo));
+               (i-1, pbo::acc, subst, metasenv))
+            cl pl (List.length pl, [], subst, metasenv)
+        in
+        (*let metasenv, subst = 
+          try 
+            NCicUnification.unify rdb metasenv subst context outty new_outty 
+          with _ -> raise (RefineFailure (lazy (localise orig_t, "try_coercions")))
+        in*)
+        let new_outty = add_pi_for_refl_m context new_outty mty m leftno rno in
+        pp (lazy ("CASE: new_outty: " ^ (NCicPp.ppterm 
+           ~metasenv ~subst ~context new_outty)));
+        let right_tys = 
+          List.map 
+            (fun t -> NCicTypeChecker.typeof ~subst ~metasenv context t) right_p in
+        let eqs = 
+          List.map2 (fun x y -> NCic.Appl[eq_refl;x;y]) (right_tys@[mty]) 
+            (right_p@[m]) in
+        let t = NCic.Appl (NCic.Match(r,new_outty,m,pl) :: eqs) 
+        in
+        metasenv, subst, t, expty (*}}}*)
+    | NCic.Prod (nameprod, src, ty),NCic.Prod (_, src2, ty2), _ -> 
+        let rec mk_fresh_name ctx firstch n =
+          let candidate = (String.make 1 firstch) ^ (string_of_int n) in
+          if (List.for_all (fun (s,_) -> s <> candidate) ctx) then candidate
+          else mk_fresh_name ctx firstch (n+1)
+        in
+        (*{{{*) pp (lazy "LAM");
+        pp (lazy ("LAM: t = " ^ NCicPp.ppterm ~metasenv ~subst ~context t));
+        let name_con = mk_fresh_name context 'c' 0
+          (*FreshNamesGenerator.mk_fresh_name 
+            ~subst metasenv context ~typ:src2 Cic.Anonymous*)
+        in
+        let context_src2 = ((name_con, NCic.Decl src2) :: context) in
+        (* contravariant part: the argument of f:src->ty *)
+        let metasenv, subst, rel1, _ = 
+          try_coercions rdb ~localise metasenv subst context_src2
+           (NCic.Rel 1) orig_t (NCicSubstitution.lift 1 src2) 
+           (NCicSubstitution.lift 1 src) perform_unification exc
+        in
+        (* covariant part: the result of f(c x); x:src2; (c x):src *)
+        let name_t, bo = 
+          match t with
+          | NCic.Lambda (n,_,bo) -> n, cs_subst rel1 (NCicSubstitution.lift_from 2 1 bo)
+          | _ -> name_con, NCicUntrusted.mk_appl (NCicSubstitution.lift 1 t) [rel1]
+        in
+        (* we fix the possible dependency problem in the source ty *)
+        let ty = cs_subst rel1 (NCicSubstitution.lift_from 2 1 ty) in
+        let metasenv, subst, bo, _ = 
+          try_coercions rdb ~localise metasenv subst context_src2
+            bo orig_t ty ty2 perform_unification exc
+        in
+        let coerced = NCic.Lambda (name_t,src2, bo) in
+        pp (lazy ("LAM: coerced = " ^ NCicPp.ppterm ~metasenv ~subst ~context coerced));
+        metasenv, subst, coerced, expty (*}}}*)
+    | _ -> raise exc
+  with exc2 ->    
   pp(lazy("try_coercion " ^ 
     NCicPp.ppterm ~metasenv ~subst ~context infty ^ " |---> " ^
     NCicPp.ppterm ~metasenv ~subst ~context expty));
@@ -416,26 +846,21 @@ and try_coercions rdb
         rdb metasenv subst context infty expty)
 
 and force_to_sort rdb metasenv subst context t orig_t localise ty =
-  match NCicReduction.whd ~subst context ty with
-  | C.Meta (_,(0,(C.Irl 0 | C.Ctx []))) as ty -> 
-     metasenv, subst, t, ty
-  | C.Meta (_i,(_,(C.Irl 0 | C.Ctx []))) -> assert false (*CSC: ???
-     metasenv, subst, t, C.Meta(i,(0,C.Irl 0)) *)
-  | C.Meta (i,(_,lc)) ->
-     let len = match lc with C.Irl len->len | C.Ctx l->List.length l in
-     let metasenv, subst, newmeta = 
-       if len > 0 then
-         NCicMetaSubst.restrict metasenv subst i (HExtlib.list_seq 1 (len+1)) 
-       else metasenv, subst, i
-     in
-     metasenv, subst, t, C.Meta (newmeta,(0,C.Irl 0))
-  | C.Sort _ as ty -> metasenv, subst, t, ty 
-  | ty -> 
-      try_coercions rdb ~localise metasenv subst context
-        t orig_t ty (NCic.Sort (NCic.Type NCicEnvironment.type0)) false 
-         (Uncertain (lazy (localise orig_t, 
-         "The type of " ^ NCicPp.ppterm ~metasenv ~subst ~context t
-         ^ " is not a sort: " ^ NCicPp.ppterm ~metasenv ~subst ~context ty)))
+  try 
+    let metasenv, subst, ty = 
+      NCicUnification.sortfy (Failure "sortfy") metasenv subst context ty in
+      metasenv, subst, t, ty
+  with
+      Failure _ -> 
+        let ty = NCicReduction.whd ~subst context ty in
+          try_coercions rdb ~localise metasenv subst context
+            t orig_t ty (NCic.Sort (NCic.Type 
+              (match NCicEnvironment.get_universes () with
+               | x::_ -> x 
+               | _ -> assert false))) false 
+             (Uncertain (lazy (localise orig_t, 
+             "The type of " ^ NCicPp.ppterm ~metasenv ~subst ~context t ^
+             " is not a sort: " ^ NCicPp.ppterm ~metasenv ~subst ~context ty)))
 
 and sort_of_prod 
   localise metasenv subst context orig_s orig_t (name,s) t (t1, t2) 
@@ -475,7 +900,9 @@ and guess_name subst ctx ty =
 
 and eat_prods rdb ~localise force_ty metasenv subst context expty orig_t orig_he he ty_he args =
   (*D*)inside 'E'; try let rc = 
-  let rec aux metasenv subst args_so_far he ty_he = function 
+  let rec aux metasenv subst args_so_far he ty_he xxx =
+  (*D*)inside 'V'; try let rc = 
+   match xxx with
   | [] ->
      let res = NCicUntrusted.mk_appl he (List.rev args_so_far) in
      pp(lazy("FORCE FINAL APPL: " ^ 
@@ -501,30 +928,30 @@ and eat_prods rdb ~localise force_ty metasenv subst context expty orig_t orig_he
            aux metasenv subst args_so_far he ty_he
             (NCic.Implicit `Term :: NCic.Implicit `Vector :: tl)
           with
-           Uncertain msg | RefineFailure msg -> raise (wrap_exc msg exc)))
+           Uncertain msg | RefineFailure msg -> raise (wrap_exc msg exc))
+     | RefineFailure msg when not (has_some_more_pis ty_he) ->
+        (* instantiating the head could change the has_some_more_pis flag *)
+        raise (Uncertain msg))
   | arg::tl ->
       match NCicReduction.whd ~subst context ty_he with 
       | C.Prod (_,s,t) ->
           let metasenv, subst, arg, _ = 
-            typeof rdb ~localise 
-              metasenv subst context arg (Some s) in
+            typeof rdb ~localise metasenv subst context arg (Some s) in
           let t = NCicSubstitution.subst ~avoid_beta_redexes:true arg t in
           aux metasenv subst (arg :: args_so_far) he t tl
       | C.Meta _
       | C.Appl (C.Meta _ :: _) as t ->
           let metasenv, subst, arg, ty_arg = 
-            typeof rdb ~localise 
-              metasenv subst context arg None in
+            typeof rdb ~localise metasenv subst context arg None in
           let name = guess_name subst context ty_arg in
           let metasenv, _, meta, _ = 
             NCicMetaSubst.mk_meta metasenv 
-              ((name,C.Decl ty_arg) :: context) `Type
+              ((name,C.Decl ty_arg) :: context) `IsType
           in
           let flex_prod = C.Prod (name, ty_arg, meta) in
           (* next line grants that ty_args is a type *)
           let metasenv,subst, flex_prod, _ =
-           typeof rdb ~localise metasenv subst
-            context flex_prod None in
+           typeof rdb ~localise metasenv subst context flex_prod None in
 (*
           pp (lazy ( "UNIFICATION in CTX:\n"^ 
             NCicPp.ppcontext ~metasenv ~subst context
@@ -568,6 +995,7 @@ and eat_prods rdb ~localise force_ty metasenv subst context expty orig_t orig_he
                (List.length args) (List.length args_so_far))))
           in
            aux metasenv subst [] newhead newheadty (arg :: tl)
+  (*D*)in outside true; rc with exc -> outside false; raise exc
   in
    (* We need to reverse the order of the new created metas since they
       are pushed on top of the metasenv in the wrong order *)
@@ -578,7 +1006,7 @@ and eat_prods rdb ~localise force_ty metasenv subst context expty orig_t orig_he
     List.partition (fun (i,_) -> i <= highest_meta) metasenv
    in
     (List.rev metasenv_new) @ metasenv_old, subst, newhead, newheadty
-  (*D*)in outside(); rc with exc -> outside (); raise exc
+  (*D*)in outside true; rc with exc -> outside false; raise exc
 ;;
 
 let rec first f l1 l2 =
@@ -610,13 +1038,15 @@ let relocalise old_localise dt t add =
 ;;
 
 let undebruijnate inductive ref t rev_fl =
+  let len = List.length rev_fl in 
   NCicSubstitution.psubst (fun x -> x) 
-   (List.rev (HExtlib.list_mapi 
+   (HExtlib.list_mapi 
       (fun (_,_,rno,_,_,_) i -> 
+         let i = len - i - 1 in
          NCic.Const 
            (if inductive then NReference.mk_fix i rno ref
             else NReference.mk_cofix i ref))
-      rev_fl))
+      rev_fl)
     t
 ;; 
 
@@ -624,18 +1054,10 @@ let undebruijnate inductive ref t rev_fl =
 let typeof_obj 
   rdb ?(localise=fun _ -> Stdpp.dummy_loc) (uri,height,metasenv,subst,obj)
 = 
-  let check_type metasenv subst context (ty as orig_ty) =  (* XXX fattorizza *)
-    let metasenv, subst, ty, sort = 
-      typeof rdb ~localise metasenv subst context ty None
-    in
-    let metasenv, subst, ty, sort = 
-      force_to_sort rdb metasenv subst context ty orig_ty localise sort
-    in
-      metasenv, subst, ty, sort
-  in
   match obj with 
-  | C.Constant (relevance, name, bo, ty , attr) ->
-       let metasenv, subst, ty, _ = check_type metasenv subst [] ty in
+  | C.Constant (relevance, name, bo, ty, attr) ->
+       let metasenv, subst, ty = 
+         check_type rdb ~localise metasenv subst [] ty in
        let metasenv, subst, bo, ty, height = 
          match bo with
          | Some bo ->
@@ -652,7 +1074,8 @@ let typeof_obj
       let types, metasenv, subst, rev_fl =
         List.fold_left
          (fun (types, metasenv, subst, fl) (relevance,name,k,ty,bo) ->
-           let metasenv, subst, ty, _ = check_type metasenv subst [] ty in
+           let metasenv, subst, ty = 
+             check_type rdb ~localise metasenv subst [] ty in
            let dbo = NCicTypeChecker.debruijn uri len [] ~subst bo in
            let localise = relocalise localise dbo bo in
             (name,C.Decl ty)::types,
@@ -688,7 +1111,8 @@ let typeof_obj
      let metasenv,subst,rev_itl,tys =
       List.fold_left
        (fun (metasenv,subst,res,ctx) (relevance,n,ty,cl) ->
-          let metasenv, subst, ty, _ = check_type metasenv subst [] ty in
+          let metasenv, subst, ty = 
+            check_type rdb ~localise metasenv subst [] ty in
           metasenv,subst,(relevance,n,ty,cl)::res,(n,NCic.Decl ty)::ctx
        ) (metasenv,subst,[],[]) itl in
      let metasenv,subst,itl,_ =
@@ -703,7 +1127,8 @@ let typeof_obj
               try snd (HExtlib.split_nth leftno k_relev)
               with Failure _ -> k_relev in
              let te = NCicTypeChecker.debruijn uri len [] ~subst te in
-             let metasenv, subst, te, _ = check_type metasenv subst tys te in
+             let metasenv, subst, te = 
+               check_type rdb ~localise metasenv subst tys te in
              let context,te = NCicReduction.split_prods ~subst tys leftno te in
              let _,chopped_context_rev =
               HExtlib.split_nth (List.length tys) (List.rev context) in
@@ -750,6 +1175,29 @@ let typeof_obj
                     metasenv,subst,item1::context
                 ) (metasenv,subst,tys) sx_context_ty_rev sx_context_te_rev
               with Invalid_argument "List.fold_left2" -> assert false in
+             let metasenv, subst =
+                let rec aux context (metasenv,subst) = function
+                  | NCic.Meta _ -> metasenv, subst
+                  | NCic.Implicit _ -> metasenv, subst
+                  | NCic.Appl (NCic.Rel i :: args) as t
+                      when i > List.length context - len ->
+                      let lefts, _ = HExtlib.split_nth leftno args in
+                      let ctxlen = List.length context in
+                      let (metasenv, subst), _ = 
+                        List.fold_left
+                        (fun ((metasenv, subst),l) arg ->
+                          NCicUnification.unify rdb 
+                           ~test_eq_only:true metasenv subst context arg 
+                             (NCic.Rel (ctxlen - len - l)), l+1
+                          )
+                        ((metasenv, subst), 0) lefts
+                      in
+                      metasenv, subst
+                  | t -> NCicUtils.fold (fun e c -> e::c) context aux
+                  (metasenv,subst) t
+                in
+               aux context (metasenv,subst) te
+             in
              let con_sort= NCicTypeChecker.typeof ~subst ~metasenv context te in
               (match
                 NCicReduction.whd ~subst context con_sort,