src/Pure/pattern.ML
author nipkow
Thu Oct 10 19:02:23 2002 +0200 (2002-10-10 ago)
changeset 13642 a3d97348ceb6
parent 13195 98975cc13d28
child 13645 430310b12c5b
permissions -rw-r--r--
added failure trace information to pattern unification
wenzelm@12784
     1
(*  Title:      Pure/pattern.ML
clasohm@0
     2
    ID:         $Id$
wenzelm@12784
     3
    Author:     Tobias Nipkow, Christine Heinzelmann, and Stefan Berghofer
wenzelm@12784
     4
    License:    GPL (GNU GENERAL PUBLIC LICENSE)
clasohm@0
     5
clasohm@0
     6
Unification of Higher-Order Patterns.
clasohm@0
     7
clasohm@0
     8
See also:
clasohm@0
     9
Tobias Nipkow. Functional Unification of Higher-Order Patterns.
clasohm@0
    10
In Proceedings of the 8th IEEE Symposium Logic in Computer Science, 1993.
nipkow@678
    11
nipkow@678
    12
TODO: optimize red by special-casing it
clasohm@0
    13
*)
clasohm@0
    14
paulson@2751
    15
infix aeconv;
paulson@2751
    16
clasohm@0
    17
signature PATTERN =
paulson@1501
    18
  sig
clasohm@0
    19
  type type_sig
clasohm@0
    20
  type sg
clasohm@0
    21
  type env
nipkow@13642
    22
  val trace_unify_fail  : bool ref
paulson@2751
    23
  val aeconv            : term * term -> bool
paulson@2725
    24
  val eta_contract      : term -> term
wenzelm@12781
    25
  val beta_eta_contract : term -> term
paulson@2725
    26
  val eta_contract_atom : term -> term
paulson@2725
    27
  val match             : type_sig -> term * term
paulson@2725
    28
                          -> (indexname*typ)list * (indexname*term)list
nipkow@4820
    29
  val first_order_match : type_sig -> term * term
nipkow@4820
    30
                          -> (indexname*typ)list * (indexname*term)list
paulson@2725
    31
  val matches           : type_sig -> term * term -> bool
nipkow@4667
    32
  val matches_subterm   : type_sig -> term * term -> bool
paulson@2725
    33
  val unify             : sg * env * (term * term)list -> env
nipkow@4820
    34
  val first_order       : term -> bool
nipkow@4820
    35
  val pattern           : term -> bool
berghofe@13195
    36
  val rewrite_term      : type_sig -> (term * term) list -> (term -> term option) list
berghofe@13195
    37
                          -> term -> term
clasohm@0
    38
  exception Unif
clasohm@0
    39
  exception MATCH
clasohm@0
    40
  exception Pattern
paulson@1501
    41
  end;
clasohm@0
    42
paulson@1501
    43
structure Pattern : PATTERN =
clasohm@0
    44
struct
clasohm@0
    45
clasohm@0
    46
type type_sig = Type.type_sig
clasohm@0
    47
type sg = Sign.sg
clasohm@0
    48
type env = Envir.env
clasohm@0
    49
clasohm@0
    50
exception Unif;
clasohm@0
    51
exception Pattern;
clasohm@0
    52
nipkow@13642
    53
val trace_unify_fail = ref false;
nipkow@13642
    54
nipkow@13642
    55
(* Only for communication between unification and error message printing *)
nipkow@13642
    56
val sgr = ref Sign.pre_pure
nipkow@13642
    57
nipkow@13642
    58
fun string_of_term env binders t = Sign.string_of_term (!sgr)
nipkow@13642
    59
  (Envir.norm_term env (subst_bounds(map Free binders,t)));
nipkow@13642
    60
nipkow@13642
    61
fun bname binders i = fst(nth_elem(i,binders));
nipkow@13642
    62
fun bnames binders is = space_implode " " (map (bname binders) is);
nipkow@13642
    63
nipkow@13642
    64
fun norm_typ tye =
nipkow@13642
    65
  let fun chase(v,s) =
nipkow@13642
    66
      (case  Vartab.lookup (tye, v) of
nipkow@13642
    67
        Some U => norm_typ tye U
nipkow@13642
    68
      | None => TVar(v,s))
nipkow@13642
    69
  in map_type_tvar chase end
nipkow@13642
    70
nipkow@13642
    71
fun typ_clash(tye,T,U) =
nipkow@13642
    72
  if !trace_unify_fail
nipkow@13642
    73
  then let val t = Sign.string_of_typ (!sgr) (norm_typ tye T)
nipkow@13642
    74
           and u = Sign.string_of_typ (!sgr) (norm_typ tye U)
nipkow@13642
    75
       in tracing("The following types do not unify:\n" ^ t ^ "\n" ^ u) end
nipkow@13642
    76
  else ()
nipkow@13642
    77
nipkow@13642
    78
fun clash a b =
nipkow@13642
    79
  if !trace_unify_fail then tracing("Clash: " ^ a ^ " =/= " ^ b) else ()
nipkow@13642
    80
nipkow@13642
    81
fun boundVar binders i =
nipkow@13642
    82
  "bound variable " ^ bname binders i ^ " (depth " ^ string_of_int i ^ ")";
nipkow@13642
    83
nipkow@13642
    84
fun clashBB binders i j =
nipkow@13642
    85
  if !trace_unify_fail then clash (boundVar binders i) (boundVar binders j)
nipkow@13642
    86
  else ()
nipkow@13642
    87
nipkow@13642
    88
fun clashB binders i s =
nipkow@13642
    89
  if !trace_unify_fail then clash (boundVar binders i) s
nipkow@13642
    90
  else ()
nipkow@13642
    91
nipkow@13642
    92
fun proj_fail(env,binders,F,is,t) =
nipkow@13642
    93
  if !trace_unify_fail
nipkow@13642
    94
  then let val f = Syntax.string_of_vname F
nipkow@13642
    95
           val xs = bnames binders is
nipkow@13642
    96
           val u = string_of_term env binders t
nipkow@13642
    97
           val ys = bnames binders (loose_bnos t \\ is)
nipkow@13642
    98
       in tracing("Cannot unify variable " ^ f ^
nipkow@13642
    99
               " (depending on bound variables " ^ xs ^ ")\nwith term " ^ u ^
nipkow@13642
   100
               "\nTerm contains additional bound variable(s) " ^ ys)
nipkow@13642
   101
       end
nipkow@13642
   102
  else ()
nipkow@13642
   103
nipkow@13642
   104
fun ocheck_fail(F,t,binders,env) =
nipkow@13642
   105
  if !trace_unify_fail
nipkow@13642
   106
  then let val f = Syntax.string_of_vname F
nipkow@13642
   107
           val u = string_of_term env binders t
nipkow@13642
   108
       in tracing("Variable " ^ f ^ " occurs in term\n" ^ u ^
nipkow@13642
   109
                  "\nCannot unify!\n")
nipkow@13642
   110
       end
nipkow@13642
   111
  else ()
nipkow@13642
   112
wenzelm@12784
   113
fun occurs(F,t,env) =
clasohm@0
   114
    let fun occ(Var(G,_))   = (case Envir.lookup(env,G) of
clasohm@0
   115
                                 Some(t) => occ t
clasohm@0
   116
                               | None    => F=G)
clasohm@0
   117
          | occ(t1$t2)      = occ t1 orelse occ t2
clasohm@0
   118
          | occ(Abs(_,_,t)) = occ t
clasohm@0
   119
          | occ _           = false
clasohm@0
   120
    in occ t end;
clasohm@0
   121
clasohm@0
   122
clasohm@0
   123
fun mapbnd f =
clasohm@0
   124
    let fun mpb d (Bound(i))     = if i < d then Bound(i) else Bound(f(i-d)+d)
clasohm@0
   125
          | mpb d (Abs(s,T,t))   = Abs(s,T,mpb(d+1) t)
nipkow@678
   126
          | mpb d ((u1 $ u2))    = (mpb d u1)$(mpb d u2)
nipkow@678
   127
          | mpb _ atom           = atom
clasohm@0
   128
    in mpb 0 end;
clasohm@0
   129
clasohm@0
   130
fun idx [] j     = ~10000
clasohm@0
   131
  | idx(i::is) j = if i=j then length is else idx is j;
clasohm@0
   132
clasohm@0
   133
fun at xs i = nth_elem (i,xs);
clasohm@0
   134
clasohm@0
   135
fun mkabs (binders,is,t)  =
clasohm@0
   136
    let fun mk(i::is) = let val (x,T) = nth_elem(i,binders)
wenzelm@12784
   137
                        in Abs(x,T,mk is) end
clasohm@0
   138
          | mk []     = t
clasohm@0
   139
    in mk is end;
clasohm@0
   140
clasohm@0
   141
val incr = mapbnd (fn i => i+1);
clasohm@0
   142
clasohm@0
   143
fun ints_of []             = []
wenzelm@12784
   144
  | ints_of (Bound i ::bs) =
clasohm@0
   145
      let val is = ints_of bs
berghofe@1576
   146
      in if i mem_int is then raise Pattern else i::is end
clasohm@0
   147
  | ints_of _              = raise Pattern;
clasohm@0
   148
berghofe@12232
   149
fun ints_of' env ts = ints_of (map (Envir.head_norm env) ts);
berghofe@12232
   150
clasohm@0
   151
clasohm@0
   152
fun app (s,(i::is)) = app (s$Bound(i),is)
clasohm@0
   153
  | app (s,[])      = s;
clasohm@0
   154
clasohm@0
   155
fun red (Abs(_,_,s)) (i::is) js = red s is (i::js)
nipkow@678
   156
  | red t            []      [] = t
nipkow@678
   157
  | red t            is      jn = app (mapbnd (at jn) t,is);
nipkow@678
   158
clasohm@0
   159
clasohm@0
   160
(* split_type ([T1,....,Tn]---> T,n,[]) = ([Tn,...,T1],T) *)
clasohm@0
   161
fun split_type (T,0,Ts)                    = (Ts,T)
clasohm@0
   162
  | split_type (Type ("fun",[T1,T2]),n,Ts) = split_type (T2,n-1,T1::Ts)
nipkow@678
   163
  | split_type _                           = error("split_type");
clasohm@0
   164
clasohm@0
   165
fun type_of_G (T,n,is) =
clasohm@0
   166
  let val (Ts,U) = split_type(T,n,[]) in map(at Ts)is ---> U end;
clasohm@0
   167
clasohm@0
   168
fun mkhnf (binders,is,G,js) = mkabs (binders, is, app(G,js));
clasohm@0
   169
clasohm@0
   170
fun mknewhnf(env,binders,is,F as (a,_),T,js) =
clasohm@0
   171
  let val (env',G) = Envir.genvar a (env,type_of_G(T,length is,js))
clasohm@0
   172
  in Envir.update((F,mkhnf(binders,is,G,js)),env') end;
clasohm@0
   173
clasohm@0
   174
clasohm@0
   175
(* mk_proj_list(is) = [ |is| - k | 1 <= k <= |is| and is[k] >= 0 ] *)
clasohm@0
   176
fun mk_proj_list is =
clasohm@0
   177
    let fun mk(i::is,j) = if i >= 0 then j :: mk(is,j-1) else mk(is,j-1)
clasohm@0
   178
          | mk([],_)    = []
clasohm@0
   179
    in mk(is,length is - 1) end;
clasohm@0
   180
clasohm@0
   181
fun proj(s,env,binders,is) =
clasohm@0
   182
    let fun trans d i = if i<d then i else (idx is (i-d))+d;
berghofe@12232
   183
        fun pr(s,env,d,binders) = (case Envir.head_norm env s of
clasohm@0
   184
              Abs(a,T,t) => let val (t',env') = pr(t,env,d+1,((a,T)::binders))
clasohm@0
   185
                            in (Abs(a,T,t'),env') end
clasohm@0
   186
            | t => (case strip_comb t of
clasohm@0
   187
                (c as Const _,ts) =>
clasohm@0
   188
                         let val (ts',env') = prs(ts,env,d,binders)
clasohm@0
   189
                         in (list_comb(c,ts'),env') end
clasohm@0
   190
                 | (f as Free _,ts) =>
clasohm@0
   191
                         let val (ts',env') = prs(ts,env,d,binders)
clasohm@0
   192
                         in (list_comb(f,ts'),env') end
clasohm@0
   193
                 | (Bound(i),ts) =>
clasohm@0
   194
                         let val j = trans d i
clasohm@0
   195
                         in if j < 0 then raise Unif
clasohm@0
   196
                            else let val (ts',env') = prs(ts,env,d,binders)
clasohm@0
   197
                                 in (list_comb(Bound j,ts'),env') end
clasohm@0
   198
                         end
clasohm@0
   199
                 | (Var(F as (a,_),Fty),ts) =>
berghofe@12232
   200
                      let val js = ints_of' env ts;
clasohm@0
   201
                          val js' = map (trans d) js;
clasohm@0
   202
                          val ks = mk_proj_list js';
clasohm@0
   203
                          val ls = filter (fn i => i >= 0) js'
clasohm@0
   204
                          val Hty = type_of_G(Fty,length js,ks)
clasohm@0
   205
                          val (env',H) = Envir.genvar a (env,Hty)
clasohm@0
   206
                          val env'' =
clasohm@0
   207
                                Envir.update((F,mkhnf(binders,js,H,ks)),env')
clasohm@0
   208
                      in (app(H,ls),env'') end
clasohm@0
   209
                 | _  => raise Pattern))
clasohm@0
   210
        and prs(s::ss,env,d,binders) =
clasohm@0
   211
              let val (s',env1) = pr(s,env,d,binders)
clasohm@0
   212
                  val (ss',env2) = prs(ss,env1,d,binders)
clasohm@0
   213
              in (s'::ss',env2) end
clasohm@0
   214
          | prs([],env,_,_) = ([],env)
clasohm@0
   215
   in if downto0(is,length binders - 1) then (s,env)
clasohm@0
   216
      else pr(s,env,0,binders)
clasohm@0
   217
   end;
clasohm@0
   218
clasohm@0
   219
clasohm@0
   220
(* mk_ff_list(is,js) = [ length(is) - k | 1 <= k <= |is| and is[k] = js[k] ] *)
wenzelm@12784
   221
fun mk_ff_list(is,js) =
wenzelm@12784
   222
    let fun mk([],[],_)        = []
clasohm@0
   223
          | mk(i::is,j::js, k) = if i=j then k :: mk(is,js,k-1)
clasohm@0
   224
                                        else mk(is,js,k-1)
nipkow@678
   225
          | mk _               = error"mk_ff_list"
clasohm@0
   226
    in mk(is,js,length is-1) end;
clasohm@0
   227
clasohm@0
   228
fun flexflex1(env,binders,F,Fty,is,js) =
clasohm@0
   229
  if is=js then env
clasohm@0
   230
  else let val ks = mk_ff_list(is,js)
clasohm@0
   231
       in mknewhnf(env,binders,is,F,Fty,ks) end;
clasohm@0
   232
clasohm@0
   233
fun flexflex2(env,binders,F,Fty,is,G,Gty,js) =
clasohm@0
   234
  let fun ff(F,Fty,is,G as (a,_),Gty,js) =
berghofe@1576
   235
            if js subset_int is
clasohm@0
   236
            then let val t= mkabs(binders,is,app(Var(G,Gty),map (idx is) js))
clasohm@0
   237
                 in Envir.update((F,t),env) end
berghofe@1576
   238
            else let val ks = is inter_int js
clasohm@0
   239
                     val Hty = type_of_G(Fty,length is,map (idx is) ks)
clasohm@0
   240
                     val (env',H) = Envir.genvar a (env,Hty)
clasohm@0
   241
                     fun lam(is) = mkabs(binders,is,app(H,map (idx is) ks));
clasohm@0
   242
                 in Envir.update((G,lam js), Envir.update((F,lam is),env'))
clasohm@0
   243
                 end;
clasohm@0
   244
  in if xless(G,F) then ff(F,Fty,is,G,Gty,js) else ff(G,Gty,js,F,Fty,is) end
clasohm@0
   245
clasohm@0
   246
val tsgr = ref(Type.tsig0);
clasohm@0
   247
clasohm@0
   248
fun unify_types(T,U, env as Envir.Envir{asol,iTs,maxidx}) =
clasohm@0
   249
  if T=U then env
wenzelm@12527
   250
  else let val (iTs',maxidx') = Type.unify (!tsgr) (iTs, maxidx) (U, T)
nipkow@1435
   251
       in Envir.Envir{asol=asol,maxidx=maxidx',iTs=iTs'} end
nipkow@13642
   252
       handle Type.TUNIFY => (typ_clash(iTs,T,U); raise Unif);
clasohm@0
   253
berghofe@12232
   254
fun unif binders (env,(s,t)) = case (Envir.head_norm env s, Envir.head_norm env t) of
clasohm@0
   255
      (Abs(ns,Ts,ts),Abs(nt,Tt,tt)) =>
clasohm@0
   256
         let val name = if ns = "" then nt else ns
clasohm@0
   257
         in unif ((name,Ts)::binders) (env,(ts,tt)) end
clasohm@0
   258
    | (Abs(ns,Ts,ts),t) => unif ((ns,Ts)::binders) (env,(ts,(incr t)$Bound(0)))
clasohm@0
   259
    | (t,Abs(nt,Tt,tt)) => unif ((nt,Tt)::binders) (env,((incr t)$Bound(0),tt))
clasohm@0
   260
    | p => cases(binders,env,p)
clasohm@0
   261
clasohm@0
   262
and cases(binders,env,(s,t)) = case (strip_comb s,strip_comb t) of
wenzelm@12784
   263
       ((Var(F,Fty),ss),(Var(G,Gty),ts)) =>
berghofe@12232
   264
         if F = G then flexflex1(env,binders,F,Fty,ints_of' env ss,ints_of' env ts)
berghofe@12232
   265
                  else flexflex2(env,binders,F,Fty,ints_of' env ss,G,Gty,ints_of' env ts)
berghofe@12232
   266
      | ((Var(F,_),ss),_)             => flexrigid(env,binders,F,ints_of' env ss,t)
berghofe@12232
   267
      | (_,(Var(F,_),ts))             => flexrigid(env,binders,F,ints_of' env ts,s)
clasohm@0
   268
      | ((Const c,ss),(Const d,ts))   => rigidrigid(env,binders,c,d,ss,ts)
clasohm@0
   269
      | ((Free(f),ss),(Free(g),ts))   => rigidrigid(env,binders,f,g,ss,ts)
wenzelm@12784
   270
      | ((Bound(i),ss),(Bound(j),ts)) => rigidrigidB (env,binders,i,j,ss,ts)
clasohm@0
   271
      | ((Abs(_),_),_)                => raise Pattern
clasohm@0
   272
      | (_,(Abs(_),_))                => raise Pattern
nipkow@13642
   273
      | ((Const(c,_),_),(Free(f,_),_)) => (clash c f; raise Unif)
nipkow@13642
   274
      | ((Const(c,_),_),(Bound i,_))   => (clashB binders i c; raise Unif)
nipkow@13642
   275
      | ((Free(f,_),_),(Const(c,_),_)) => (clash f c; raise Unif)
nipkow@13642
   276
      | ((Free(f,_),_),(Bound i,_))    => (clashB binders i f; raise Unif)
nipkow@13642
   277
      | ((Bound i,_),(Const(c,_),_))   => (clashB binders i c; raise Unif)
nipkow@13642
   278
      | ((Bound i,_),(Free(f,_),_))    => (clashB binders i f; raise Unif)
nipkow@13642
   279
clasohm@0
   280
clasohm@0
   281
and rigidrigid (env,binders,(a,Ta),(b,Tb),ss,ts) =
nipkow@13642
   282
      if a<>b then (clash a b; raise Unif)
clasohm@0
   283
      else foldl (unif binders) (unify_types(Ta,Tb,env), ss~~ts)
clasohm@0
   284
clasohm@0
   285
and rigidrigidB (env,binders,i,j,ss,ts) =
nipkow@13642
   286
     if i <> j then (clashBB binders i j; raise Unif)
nipkow@13642
   287
     else foldl (unif binders) (env ,ss~~ts)
clasohm@0
   288
nipkow@13642
   289
and flexrigid (params as (env,binders,F,is,t)) =
nipkow@13642
   290
      if occurs(F,t,env) then (ocheck_fail(F,t,binders,env); raise Unif)
nipkow@13642
   291
      else (let val (u,env') = proj(t,env,binders,is)
nipkow@13642
   292
            in Envir.update((F,mkabs(binders,is,u)),env') end
nipkow@13642
   293
            handle Unif => (proj_fail params; raise Unif));
clasohm@0
   294
nipkow@13642
   295
fun unify(sg,env,tus) = (sgr := sg; tsgr := #tsig(Sign.rep_sg sg);
clasohm@0
   296
                         foldl (unif []) (env,tus));
clasohm@0
   297
clasohm@0
   298
paulson@2725
   299
(*Eta-contract a term (fully)*)
nipkow@2792
   300
wenzelm@12797
   301
fun eta_contract t =
wenzelm@12797
   302
  let
wenzelm@12797
   303
    exception SAME;
wenzelm@12797
   304
    fun eta (Abs (a, T, body)) = 
wenzelm@12797
   305
      ((case eta body of
wenzelm@12797
   306
          body' as (f $ Bound 0) => 
wenzelm@12797
   307
	    if loose_bvar1 (f, 0) then Abs(a, T, body')
wenzelm@12797
   308
	    else incr_boundvars ~1 f
wenzelm@12797
   309
        | body' => Abs (a, T, body')) handle SAME =>
wenzelm@12797
   310
       (case body of
wenzelm@12797
   311
          (f $ Bound 0) => 
wenzelm@12797
   312
	    if loose_bvar1 (f, 0) then raise SAME
wenzelm@12797
   313
	    else incr_boundvars ~1 f
wenzelm@12797
   314
        | _ => raise SAME))
wenzelm@12797
   315
      | eta (f $ t) =
wenzelm@12797
   316
          (let val f' = eta f
wenzelm@12797
   317
           in f' $ etah t end handle SAME => f $ eta t)
wenzelm@12797
   318
      | eta _ = raise SAME
wenzelm@12797
   319
    and etah t = (eta t handle SAME => t)
wenzelm@12797
   320
  in etah t end;
clasohm@0
   321
wenzelm@12781
   322
val beta_eta_contract = eta_contract o Envir.beta_norm;
wenzelm@12781
   323
nipkow@6539
   324
(*Eta-contract a term from outside: just enough to reduce it to an atom
nipkow@6539
   325
DOESN'T QUITE WORK!
nipkow@6539
   326
*)
wenzelm@12784
   327
fun eta_contract_atom (t0 as Abs(a, T, body)) =
paulson@2725
   328
      (case  eta_contract2 body  of
wenzelm@12784
   329
        body' as (f $ Bound 0) =>
wenzelm@12784
   330
            if loose_bvar1(f,0) then Abs(a,T,body')
wenzelm@12784
   331
            else eta_contract_atom (incr_boundvars ~1 f)
paulson@2725
   332
      | _ => t0)
paulson@2725
   333
  | eta_contract_atom t = t
paulson@2725
   334
and eta_contract2 (f$t) = f $ eta_contract_atom t
paulson@2725
   335
  | eta_contract2 t     = eta_contract_atom t;
paulson@2725
   336
paulson@2725
   337
nipkow@678
   338
(*Tests whether 2 terms are alpha/eta-convertible and have same type.
nipkow@678
   339
  Note that Consts and Vars may have more than one type.*)
paulson@2751
   340
fun t aeconv u = aconv_aux (eta_contract_atom t, eta_contract_atom u)
paulson@2751
   341
and aconv_aux (Const(a,T), Const(b,U)) = a=b  andalso  T=U
paulson@2751
   342
  | aconv_aux (Free(a,T),  Free(b,U))  = a=b  andalso  T=U
paulson@2751
   343
  | aconv_aux (Var(v,T),   Var(w,U))   = eq_ix(v,w) andalso  T=U
wenzelm@12784
   344
  | aconv_aux (Bound i,    Bound j)    = i=j
paulson@2751
   345
  | aconv_aux (Abs(_,T,t), Abs(_,U,u)) = (t aeconv u)  andalso  T=U
wenzelm@12784
   346
  | aconv_aux (f$t,        g$u)        = (f aeconv g)  andalso (t aeconv u)
paulson@2751
   347
  | aconv_aux _ =  false;
nipkow@678
   348
nipkow@678
   349
nipkow@4820
   350
(*** Matching ***)
nipkow@4820
   351
nipkow@4820
   352
exception MATCH;
nipkow@4820
   353
nipkow@4820
   354
fun typ_match tsig args = (Type.typ_match tsig args)
nipkow@4820
   355
                          handle Type.TYPE_MATCH => raise MATCH;
nipkow@4820
   356
nipkow@4820
   357
(*First-order matching;
nipkow@4820
   358
  fomatch tsig (pattern, object) returns a (tyvar,typ)list and (var,term)list.
nipkow@4820
   359
  The pattern and object may have variables in common.
nipkow@4820
   360
  Instantiation does not affect the object, so matching ?a with ?a+1 works.
nipkow@4820
   361
  Object is eta-contracted on the fly (by eta-expanding the pattern).
nipkow@4820
   362
  Precondition: the pattern is already eta-contracted!
nipkow@4820
   363
  Note: types are matched on the fly *)
nipkow@4820
   364
fun fomatch tsig =
nipkow@4820
   365
  let
nipkow@4820
   366
    fun mtch (instsp as (tyinsts,insts)) = fn
nipkow@4820
   367
        (Var(ixn,T), t)  =>
wenzelm@12784
   368
          if loose_bvar(t,0) then raise MATCH
wenzelm@12784
   369
          else (case assoc_string_int(insts,ixn) of
wenzelm@12784
   370
                  None => (typ_match tsig (tyinsts, (T, fastype_of t)),
wenzelm@12784
   371
                           (ixn,t)::insts)
wenzelm@12784
   372
                | Some u => if t aeconv u then instsp else raise MATCH)
nipkow@4820
   373
      | (Free (a,T), Free (b,U)) =>
wenzelm@12784
   374
          if a=b then (typ_match tsig (tyinsts,(T,U)), insts) else raise MATCH
nipkow@4820
   375
      | (Const (a,T), Const (b,U))  =>
wenzelm@12784
   376
          if a=b then (typ_match tsig (tyinsts,(T,U)), insts) else raise MATCH
nipkow@4820
   377
      | (Bound i, Bound j)  =>  if  i=j  then  instsp  else raise MATCH
nipkow@4820
   378
      | (Abs(_,T,t), Abs(_,U,u))  =>
wenzelm@12784
   379
          mtch (typ_match tsig (tyinsts,(T,U)),insts) (t,u)
nipkow@4820
   380
      | (f$t, g$u) => mtch (mtch instsp (f,g)) (t, u)
nipkow@4820
   381
      | (t, Abs(_,U,u))  =>  mtch instsp ((incr t)$(Bound 0), u)
nipkow@4820
   382
      | _ => raise MATCH
nipkow@4820
   383
  in mtch end;
nipkow@4820
   384
berghofe@8406
   385
fun first_order_match tsig = apfst Vartab.dest o fomatch tsig (Vartab.empty, []);
berghofe@8406
   386
nipkow@4820
   387
(* Matching of higher-order patterns *)
nipkow@4820
   388
nipkow@4820
   389
fun match_bind(itms,binders,ixn,is,t) =
nipkow@4820
   390
  let val js = loose_bnos t
nipkow@4820
   391
  in if null is
nipkow@4820
   392
     then if null js then (ixn,t)::itms else raise MATCH
nipkow@4820
   393
     else if js subset_int is
nipkow@4820
   394
          then let val t' = if downto0(is,length binders - 1) then t
nipkow@4820
   395
                            else mapbnd (idx is) t
nipkow@4820
   396
               in (ixn, mkabs(binders,is,t')) :: itms end
nipkow@4820
   397
          else raise MATCH
nipkow@4820
   398
  end;
nipkow@4820
   399
nipkow@678
   400
fun match tsg (po as (pat,obj)) =
nipkow@678
   401
let
nipkow@4820
   402
  (* Pre: pat and obj have same type *)
nipkow@4820
   403
  fun mtch binders (env as (iTs,itms),(pat,obj)) =
nipkow@4820
   404
    case pat of
nipkow@4820
   405
      Abs(ns,Ts,ts) =>
nipkow@4820
   406
        (case obj of
nipkow@4820
   407
           Abs(nt,Tt,tt) => mtch ((nt,Tt)::binders) (env,(ts,tt))
berghofe@8406
   408
         | _ => let val Tt = typ_subst_TVars_Vartab iTs Ts
nipkow@4820
   409
                in mtch((ns,Tt)::binders)(env,(ts,(incr obj)$Bound(0))) end)
nipkow@4820
   410
    | _ => (case obj of
nipkow@4820
   411
              Abs(nt,Tt,tt) =>
nipkow@4820
   412
                mtch((nt,Tt)::binders)(env,((incr pat)$Bound(0),tt))
nipkow@4820
   413
            | _ => cases(binders,env,pat,obj))
nipkow@678
   414
nipkow@4820
   415
  and cases(binders,env as (iTs,itms),pat,obj) =
nipkow@4820
   416
    let val (ph,pargs) = strip_comb pat
nipkow@4820
   417
        fun rigrig1(iTs,oargs) =
nipkow@4820
   418
              foldl (mtch binders) ((iTs,itms), pargs~~oargs)
nipkow@4820
   419
        fun rigrig2((a,Ta),(b,Tb),oargs) =
nipkow@4820
   420
              if a<> b then raise MATCH
nipkow@4820
   421
              else rigrig1(typ_match tsg (iTs,(Ta,Tb)), oargs)
nipkow@4820
   422
    in case ph of
nipkow@4820
   423
         Var(ixn,_) =>
nipkow@4820
   424
           let val is = ints_of pargs
nipkow@4820
   425
           in case assoc_string_int(itms,ixn) of
nipkow@4820
   426
                None => (iTs,match_bind(itms,binders,ixn,is,obj))
nipkow@4820
   427
              | Some u => if obj aeconv (red u is []) then env
nipkow@4820
   428
                          else raise MATCH
nipkow@4820
   429
           end
nipkow@4820
   430
       | _ =>
nipkow@4820
   431
           let val (oh,oargs) = strip_comb obj
nipkow@4820
   432
           in case (ph,oh) of
nipkow@4820
   433
                (Const c,Const d) => rigrig2(c,d,oargs)
nipkow@4820
   434
              | (Free f,Free g)   => rigrig2(f,g,oargs)
nipkow@4820
   435
              | (Bound i,Bound j) => if i<>j then raise MATCH
nipkow@4820
   436
                                     else rigrig1(iTs,oargs)
nipkow@4820
   437
              | (Abs _, _)        => raise Pattern
nipkow@4820
   438
              | (_, Abs _)        => raise Pattern
nipkow@4820
   439
              | _                 => raise MATCH
nipkow@4820
   440
           end
nipkow@4820
   441
    end;
nipkow@678
   442
nipkow@4820
   443
  val pT = fastype_of pat
nipkow@4820
   444
  and oT = fastype_of obj
berghofe@8406
   445
  val iTs = typ_match tsg (Vartab.empty, (pT,oT))
nipkow@4820
   446
  val insts2 = (iTs,[])
nipkow@678
   447
berghofe@8406
   448
in apfst Vartab.dest (mtch [] (insts2, po)
berghofe@8406
   449
   handle Pattern => fomatch tsg insts2 po)
nipkow@678
   450
end;
clasohm@0
   451
clasohm@0
   452
(*Predicate: does the pattern match the object?*)
nipkow@678
   453
fun matches tsig po = (match tsig po; true) handle MATCH => false;
clasohm@0
   454
nipkow@4667
   455
(* Does pat match a subterm of obj? *)
nipkow@4667
   456
fun matches_subterm tsig (pat,obj) =
nipkow@4667
   457
  let fun msub(bounds,obj) = matches tsig (pat,obj) orelse
nipkow@4667
   458
            case obj of
nipkow@4667
   459
              Abs(x,T,t) => let val y = variant bounds x
nipkow@4667
   460
                                val f = Free(":" ^ y,T)
nipkow@4667
   461
                            in msub(x::bounds,subst_bound(f,t)) end
nipkow@4667
   462
            | s$t => msub(bounds,s) orelse msub(bounds,t)
nipkow@4667
   463
            | _ => false
nipkow@4667
   464
  in msub([],obj) end;
nipkow@4667
   465
nipkow@4820
   466
fun first_order(Abs(_,_,t)) = first_order t
nipkow@4820
   467
  | first_order(t $ u) = first_order t andalso first_order u andalso
nipkow@4820
   468
                         not(is_Var t)
nipkow@4820
   469
  | first_order _ = true;
nipkow@4820
   470
nipkow@4820
   471
fun pattern(Abs(_,_,t)) = pattern t
nipkow@4820
   472
  | pattern(t) = let val (head,args) = strip_comb t
nipkow@4820
   473
                 in if is_Var head
nipkow@4820
   474
                    then let val _ = ints_of args in true end
nipkow@4820
   475
                         handle Pattern => false
nipkow@4820
   476
                    else forall pattern args
nipkow@4820
   477
                 end;
nipkow@4820
   478
wenzelm@12784
   479
wenzelm@12784
   480
(* rewriting -- simple but fast *)
wenzelm@12784
   481
berghofe@13195
   482
fun rewrite_term tsig rules procs tm =
wenzelm@12784
   483
  let
berghofe@13195
   484
    val skel0 = Bound 0;
berghofe@13195
   485
wenzelm@12797
   486
    val rhs_names =
wenzelm@12797
   487
      foldr (fn ((_, rhs), names) => add_term_free_names (rhs, names)) (rules, []);
berghofe@13195
   488
wenzelm@12797
   489
    fun variant_absfree (x, T, t) =
wenzelm@12797
   490
      let
wenzelm@12797
   491
        val x' = variant (add_term_free_names (t, rhs_names)) x;
wenzelm@12797
   492
        val t' = subst_bound (Free (x', T), t);
wenzelm@12797
   493
      in (fn u => Abs (x, T, abstract_over (Free (x', T), u)), t') end;
wenzelm@12797
   494
wenzelm@12784
   495
    fun match_rew tm (tm1, tm2) =
berghofe@13195
   496
      let val rtm = if_none (Term.rename_abs tm1 tm tm2) tm2
berghofe@13195
   497
      in Some (subst_vars (match tsig (tm1, tm)) rtm, rtm)
berghofe@13195
   498
        handle MATCH => None
berghofe@13195
   499
      end;
wenzelm@12784
   500
berghofe@13195
   501
    fun rew (Abs (_, _, body) $ t) = Some (subst_bound (t, body), skel0)
berghofe@13195
   502
      | rew tm = (case get_first (match_rew tm) rules of
berghofe@13195
   503
          None => apsome (rpair skel0) (get_first (fn p => p tm) procs)
berghofe@13195
   504
        | x => x);
berghofe@13195
   505
berghofe@13195
   506
    fun rew1 (Var _) _ = None
berghofe@13195
   507
      | rew1 skel tm = (case rew2 skel tm of
wenzelm@12817
   508
          Some tm1 => (case rew tm1 of
berghofe@13195
   509
              Some (tm2, skel') => Some (if_none (rew1 skel' tm2) tm2)
wenzelm@12784
   510
            | None => Some tm1)
wenzelm@12817
   511
        | None => (case rew tm of
berghofe@13195
   512
              Some (tm1, skel') => Some (if_none (rew1 skel' tm1) tm1)
wenzelm@12784
   513
            | None => None))
wenzelm@12784
   514
berghofe@13195
   515
    and rew2 skel (tm1 $ tm2) = (case tm1 of
wenzelm@12784
   516
            Abs (_, _, body) =>
wenzelm@12784
   517
              let val tm' = subst_bound (tm2, body)
berghofe@13195
   518
              in Some (if_none (rew2 skel0 tm') tm') end
berghofe@13195
   519
          | _ => 
berghofe@13195
   520
            let val (skel1, skel2) = (case skel of
berghofe@13195
   521
                skel1 $ skel2 => (skel1, skel2)
berghofe@13195
   522
              | _ => (skel0, skel0))
berghofe@13195
   523
            in case rew1 skel1 tm1 of
berghofe@13195
   524
                Some tm1' => (case rew1 skel2 tm2 of
berghofe@13195
   525
                    Some tm2' => Some (tm1' $ tm2')
berghofe@13195
   526
                  | None => Some (tm1' $ tm2))
berghofe@13195
   527
              | None => (case rew1 skel2 tm2 of
berghofe@13195
   528
                    Some tm2' => Some (tm1 $ tm2')
berghofe@13195
   529
                  | None => None)
berghofe@13195
   530
            end)
berghofe@13195
   531
      | rew2 skel (Abs (x, T, tm)) =
berghofe@13195
   532
          let
berghofe@13195
   533
            val (abs, tm') = variant_absfree (x, T, tm);
berghofe@13195
   534
            val skel' = (case skel of Abs (_, _, skel') => skel' | _ => skel0)
berghofe@13195
   535
          in case rew1 skel' tm' of
wenzelm@12797
   536
              Some tm'' => Some (abs tm'')
berghofe@13195
   537
            | None => None
wenzelm@12797
   538
          end
berghofe@13195
   539
      | rew2 _ _ = None
wenzelm@12784
   540
berghofe@13195
   541
  in if_none (rew1 skel0 tm) tm end;
wenzelm@12784
   542
clasohm@0
   543
end;
nipkow@13642
   544
nipkow@13642
   545
val trace_unify_fail = Pattern.trace_unify_fail;