src/Tools/isac/MathEngBasic/rewrite.sml
author wenzelm
Mon, 21 Jun 2021 15:36:09 +0200
changeset 60309 70a1d102660d
parent 60269 74965ce81297
child 60331 40eb8aa2b0d6
permissions -rw-r--r--
more antiquotations for Isabelle/HOL consts/types, without change of semantics;
     1 (* isac's rewriter
     2    (c) Walther Neuper 2000
     3 *)
     4 
     5 signature REWRITE =
     6 sig
     7   exception NO_REWRITE
     8   val calculate_: theory -> string * Eval_Def.eval_fn -> term -> (term * (string * thm)) option
     9   val eval__true: theory -> int -> term list -> (term * term) list -> Rule_Set.T -> term list * bool
    10   val eval_prog_expr: theory -> Rule_Set.T -> term -> term
    11   val eval_true_: theory -> Rule_Set.T -> term -> bool
    12   val eval_true: theory -> term list -> Rule_Set.T -> bool
    13   val rew_sub: theory -> int -> (term * term) list -> Rule_Def.rew_ord_
    14     -> Rule_Set.T -> bool -> TermC.path -> term -> term -> term * term list * TermC.path * bool
    15   val rewrite_: theory -> Rule_Def.rew_ord_ -> Rule_Set.T -> bool -> thm ->
    16     term -> (term * term list) option
    17   val rewrite_inst_: theory -> Rule_Def.rew_ord_ -> Rule_Set.T -> bool
    18     -> (term * term) list -> thm -> term -> (term * term list) option
    19   val rewrite_set_: theory -> bool -> Rule_Set.T -> term -> (term * term list) option
    20   val rewrite_set_inst_: theory -> bool -> (term * term) list -> Rule_Set.T -> term -> (term * term list) option
    21   val rewrite_terms_: theory -> Rule_Def.rew_ord_ -> Rule_Set.T -> term list
    22     -> term -> (term * term list) option
    23 
    24   val trace_on: bool Unsynchronized.ref
    25   val depth: int Unsynchronized.ref
    26   val lim_deriv: int Unsynchronized.ref
    27 
    28 \<^isac_test>\<open>
    29   val rewrite__: theory -> int -> (term * term) list -> Rule_Def.rew_ord_ ->
    30     Rule_Set.T -> bool -> thm -> term -> (term * term list) option
    31   val rewrite__set_: theory -> int -> bool -> (term * term) list -> Rule_Set.T -> term -> (term * term list) option
    32   val app_rev: theory -> int -> Rule_Set.T -> term -> term * term list * bool
    33   val app_sub: theory -> int -> Rule_Set.T -> term -> term * term list * bool
    34   val trace1: int -> string -> unit
    35   val trace_eq1 : int -> string -> Rule_Def.rule_set -> theory -> term -> unit;
    36   val trace_eq2 : int -> string -> theory -> term -> term -> unit;
    37   val trace_in1 : int -> string -> string -> unit;
    38   val trace_in2 : int -> string -> theory -> term -> unit;
    39   val trace_in3 : int -> string -> theory -> (term * 'a) option -> unit;
    40   val trace_in4 : int -> string -> theory -> term list -> term list -> unit;
    41   val trace_in5 : int -> string -> theory -> term list -> unit;
    42 \<close>
    43 end
    44 
    45 (**)
    46 structure Rewrite(**): REWRITE(**) =
    47 struct
    48 (**)
    49 
    50 exception NO_REWRITE;
    51 
    52 val trace_on = Unsynchronized.ref false;
    53 (* depth of recursion in traces of the rewriter, if trace_on:=true *)
    54 val depth = Unsynchronized.ref 99999;
    55 (* no of rewrites exceeding this int -> NO rewrite *)
    56 val lim_deriv = Unsynchronized.ref 100;
    57 
    58 fun trace i str = 
    59   if ! trace_on andalso i < ! depth then tracing (idt "#" i ^ str) else ()
    60 fun trace_eq1 i str rrls thy t =
    61   trace i (" " ^ str ^ ": " ^ Rule_Set.id rrls ^ " on: " ^ UnparseC.term_in_thy thy t)
    62 fun trace_eq2 i str thy t t' =
    63   trace i (" " ^ str ^ ": \"" ^
    64     UnparseC.term_in_thy thy t ^ "\" > \"" ^ UnparseC.term_in_thy thy t' ^ "\"");
    65 
    66 fun trace1 i str =
    67   if ! trace_on andalso i < ! depth then tracing (idt "#" (i + 1) ^ str) else ()
    68 fun trace_in1 i str thmid =
    69   trace1 i (" " ^ str ^ ": \"" ^ thmid ^ "\"")
    70 fun trace_in2 i str thy t =
    71   trace1 i (" " ^ str ^ ": \"" ^ UnparseC.term_in_thy thy t ^ "\"");
    72 fun trace_in3 i str thy pairopt =
    73   trace1 i (" " ^ str ^ ": " ^ UnparseC.term_in_thy thy ((fst o the) pairopt));
    74 fun trace_in4 i str thy ts ts' =
    75   if ! trace_on andalso i < ! depth andalso ts <> []
    76   then tracing (idt "#" (i + 1) ^ " " ^ str ^ ": " ^ UnparseC.terms_in_thy thy ts ^
    77   	"   stored: " ^ UnparseC.terms_in_thy thy ts')
    78   else ();
    79 fun trace_in5 i str thy p' =
    80   if ! trace_on andalso i < ! depth 
    81   then tracing (idt "#" (i + 1) ^ " " ^ str ^ ": " ^ UnparseC.terms_in_thy thy p')
    82   else();
    83 
    84 fun rewrite__ thy i bdv tless rls put_asm thm ct =
    85   let
    86     val (t', asms, _(*lrd*), rew) = rew_sub thy i bdv tless rls put_asm ([(*root of the term*)]: TermC.path)
    87 		  (TermC.inst_bdv bdv (Eval.norm (Thm.prop_of thm))) ct
    88   in if rew then SOME (t', distinct op = asms) else NONE end
    89   (* one rewrite (possibly conditional, ordered) EXOR exn EXOR go into subterms *)
    90 and rew_sub thy i bdv tless rls put_asm lrd r t = 
    91   (let
    92     val (lhs, rhs) = (HOLogic.dest_eq o HOLogic.dest_Trueprop o Logic.strip_imp_concl) r
    93     val r' = (Envir.subst_term (Pattern.match thy (lhs, t) (Vartab.empty, Vartab.empty)) r)
    94       handle Pattern.MATCH => raise NO_REWRITE
    95     val p' = map HOLogic.dest_Trueprop ((fst o Logic.strip_prems) (Logic.count_prems r', [], r'))
    96     val t' = (snd o HOLogic.dest_eq o HOLogic.dest_Trueprop o Logic.strip_imp_concl) r'
    97     val _ = trace_in2 i "eval asms" thy r';
    98     val (t'', p'') =                                                      (*conditional rewriting*)
    99       let val (simpl_p', nofalse) = eval__true thy (i + 1) p' bdv rls 	     
   100 	    in
   101 	      if nofalse
   102         then (trace_in4 i "asms accepted" thy p' simpl_p'; (t',simpl_p'))(* uncond.rew.from above*)
   103         else (trace_in5 i "asms false" thy p'; raise NO_REWRITE)   (* don't go into subtm.of cond*)
   104 	    end
   105   in
   106     if TermC.perm lhs rhs andalso not (tless bdv (t', t))                     (*ordered rewriting*)
   107     then (trace_eq2 i "not >" thy t t'; raise NO_REWRITE)
   108     else (t'', p'', [], true)
   109   end
   110   ) handle NO_REWRITE =>
   111     (case t of
   112       Const(s, T) => (Const(s, T), [], lrd, false)
   113     | Free(s, T) => (Free(s, T), [], lrd, false)
   114     | Var(n, T) => (Var(n, T), [], lrd, false)
   115     | Bound i => (Bound i, [], lrd, false)
   116     | Abs(s, T, body) => 
   117       let val (t', asms, _ (*lrd*), rew) =  rew_sub thy i bdv tless rls put_asm (lrd @ [TermC.D]) r body
   118        in (Abs(s, T, t'), asms, [], rew) end
   119     | t1 $ t2 => 
   120        let val (t2', asm2, lrd, rew2) = rew_sub thy i bdv tless rls put_asm (lrd @ [TermC.R]) r t2
   121        in
   122         if rew2 then (t1 $ t2', asm2, lrd, true)
   123         else
   124           let val (t1', asm1, lrd, rew1) = rew_sub thy i bdv tless rls put_asm (lrd @ [TermC.L]) r t1
   125           in if rew1 then (t1' $ t2, asm1, lrd, true) else (t1 $ t2,[], lrd, false) end
   126     end)
   127 and eval__true thy i asms bdv rls =           (* simplify asumptions until one evaluates to false*)
   128   if asms = [@{term True}] orelse asms = [] then ([], true)
   129   else (* this allows to check Rrls with prepat = ([@{term True}], pat) *)
   130     if asms = [@{term False}] then ([], false)
   131     else
   132       let                            
   133         fun chk indets [] = (indets, true) (*return asms<>True until false*)
   134           | chk indets (a :: asms) =
   135             (case rewrite__set_ thy (i + 1) false bdv rls a of
   136               NONE => (chk (indets @ [a]) asms)
   137             | SOME (t, a') =>
   138               if t = @{term True} then (chk (indets @ a') asms) 
   139               else if t = @{term False} then ([], false)
   140             (*asm false .. thm not applied ^^^; continue until False vvv*)
   141             else chk (indets @ [t] @ a') asms);
   142       in chk [] asms end
   143 and rewrite__set_ thy _ _ _ Rule_Set.Empty t =                          (* rewrite with a rule set*)
   144     raise ERROR ("rewrite__set_ called with 'Erls' for '" ^ UnparseC.term_in_thy thy t ^ "'")
   145   | rewrite__set_ thy i _ _ (rrls as Rule_Set.Rrls _) t =    (* rewrite with a 'reverse rule set'*)
   146     let
   147       val _= trace_eq1 i "rls" rrls thy t;
   148 	    val (t', asm, rew) = app_rev thy (i + 1) rrls t                   
   149     in if rew then SOME (t', distinct op = asm) else NONE end
   150   | rewrite__set_ thy i put_asm bdv rls ct =           (* Rls, Seq containing Thms or Eval, Cal1 *)
   151     let
   152       (* attention with cp to test/..: unbound thy, i, bdv, rls; TODO1803? pull out to rewrite__*)
   153       datatype switch = Appl | Noap;
   154       fun rew_once _ asm ct Noap [] = (ct, asm) (* ?TODO unify with Prog_Expr.rew_once? *)
   155         | rew_once ruls asm ct Appl [] = 
   156           (case rls of Rule_Def.Repeat _ => rew_once ruls asm ct Noap ruls
   157           | Rule_Set.Sequence _ => (ct, asm)
   158           | rls => raise ERROR ("rew_once not appl. to \"" ^ Rule_Set.id rls ^ "\""))
   159         | rew_once ruls asm ct apno (rul :: thms) =
   160           case rul of
   161             Rule.Thm (thmid, thm) =>
   162               (trace_in1 i "try thm" thmid;
   163               case rewrite__ thy (i + 1) bdv ((snd o #rew_ord o Rule_Set.rep) rls)
   164                   ((#erls o Rule_Set.rep) rls) put_asm thm ct of
   165                 NONE => rew_once ruls asm ct apno thms
   166               | SOME (ct', asm') => 
   167                 (trace_in2 i "rewrites to" thy ct';
   168                 rew_once ruls (union (op =) asm asm') ct' Appl (rul :: thms)))
   169                 (* once again try the same rule, e.g. associativity against "()"*)
   170           | Rule.Eval (cc as (op_, _)) => 
   171             let val _= trace_in1 i "try calc" op_;
   172               val ct = TermC.uminus_to_string ct (*WN190312: superfluous?*)
   173             in case Eval.adhoc_thm thy cc ct of
   174                 NONE => rew_once ruls asm ct apno thms
   175               | SOME (_, thm') => 
   176                 let 
   177                   val pairopt = rewrite__ thy (i + 1) bdv ((snd o #rew_ord o Rule_Set.rep) rls)
   178                     ((#erls o Rule_Set.rep) rls) put_asm thm' ct;
   179                   val _ = if pairopt <> NONE then () else raise ERROR ("rewrite_set_, rewrite_ \"" ^ 
   180                     ThmC.string_of_thm thm' ^ "\" " ^ UnparseC.term_in_thy thy ct ^ " = NONE")
   181                   val _ = trace_in3 i "calc. to" thy pairopt;
   182                 in rew_once ruls asm ((fst o the) pairopt) Appl (rul :: thms) end
   183             end
   184           | Rule.Cal1 (cc as (op_, _)) => 
   185             let val _= trace_in1 i "try cal1" op_;
   186               val ct = TermC.uminus_to_string ct
   187             in case Eval.adhoc_thm1_ thy cc ct of
   188                 NONE => (ct, asm)
   189               | SOME (_, thm') =>
   190                 let 
   191                   val pairopt = rewrite__ thy (i + 1) bdv ((snd o #rew_ord o Rule_Set.rep) rls)
   192                     ((#erls o Rule_Set.rep) rls) put_asm thm' ct;
   193                   val _ = if pairopt <> NONE then () else raise ERROR ("rewrite_set_, rewrite_ \"" ^
   194                      ThmC.string_of_thm thm' ^ "\" " ^ UnparseC.term_in_thy thy ct ^ " = NONE")
   195                   val _ = trace_in3 i "cal1. to" thy pairopt;
   196                 in the pairopt end
   197             end
   198           | Rule.Rls_ rls' => 
   199             (case rewrite__set_ thy (i + 1) put_asm bdv rls' ct of
   200               SOME (t', asm') => rew_once ruls (union (op =) asm asm') t' Appl thms
   201             | NONE => rew_once ruls asm ct apno thms)
   202           | r => raise ERROR ("rew_once not appl. to \"" ^ Rule.to_string r ^ "\"");
   203       val ruls = (#rules o Rule_Set.rep) rls;
   204       val _ = trace_eq1 i "rls" rls thy ct
   205       val (ct', asm') = rew_once ruls [] ct Noap ruls;
   206 	  in if ct = ct' then NONE else SOME (ct', distinct op =  asm') end
   207 (*-------------------------------------------------------------*)
   208 and app_rev thy i rrls t =             (* apply an Rrls; if not applicable proceed with subterms*)
   209   let (* check a (precond, pattern) of a rev-set; stops with 1st true *)
   210     fun chk_prepat _ _ [] _ = true
   211       | chk_prepat thy erls prepat t =
   212         let
   213           fun chk (pres, pat) =
   214             (let 
   215               val subst: Type.tyenv * Envir.tenv =
   216                 Pattern.match thy (pat, t) (Vartab.empty, Vartab.empty)
   217              in
   218               snd (eval__true thy (i + 1) (map (Envir.subst_term subst) pres) [] erls)
   219              end) handle Pattern.MATCH => false
   220            fun scan_ _ [] = false
   221              | scan_ f (pp :: pps) =
   222                if f pp then true else scan_ f pps;
   223         in scan_ chk prepat end;
   224     (* apply the normal_form of a rev-set *)
   225     fun app_rev' thy (Rule_Set.Rrls {erls, prepat, scr = Rule.Rfuns {normal_form, ...}, ...}) t =
   226       if chk_prepat thy erls prepat t then normal_form t else NONE
   227       | app_rev' _ r _ = raise ERROR ("app_rev' not appl. to \"" ^ Rule_Set.id r ^ "\"");
   228     val opt = app_rev' thy rrls t
   229   in
   230     case opt of
   231       SOME (t', asm) => (t', asm, true)
   232     | NONE => app_sub thy i rrls t
   233   end
   234 and app_sub thy i rrls t =                                          (* apply an Rrls to subterms*)
   235   case t of
   236     Const (s, T) => (Const(s, T), [], false)
   237   | Free (s, T) => (Free(s, T), [], false)
   238   | Var (n, T) => (Var(n, T), [], false)
   239   | Bound i => (Bound i, [], false)
   240   | Abs (s, T, body) => 
   241 	  let val (t', asm, rew) = app_rev thy i rrls body
   242 	  in (Abs(s, T, t'), asm, rew) end
   243   | t1 $ t2 => 
   244     let val (t2', asm2, rew2) = app_rev thy i rrls t2
   245     in
   246       if rew2 then (t1 $ t2', asm2, true)
   247       else
   248         let val (t1', asm1, rew1) = app_rev thy i rrls t1
   249         in if rew1 then (t1' $ t2, asm1, true)
   250            else (t1 $ t2, [], false)
   251         end
   252     end;
   253 
   254 (* rewriting without argument [] for rew_ord;  WN110603: shouldnt asm<>[] lead to false? *)
   255 fun eval_true thy terms rls = (snd o (eval__true thy 1 terms [])) rls;
   256 
   257 (* rewriting without internal argument [] *)
   258 fun rewrite_ thy rew_ord erls bool thm term = rewrite__ thy 1 [] rew_ord erls bool thm term;
   259 fun rewrite_set_ thy bool rls term = rewrite__set_ thy 1 bool [] rls term;
   260 
   261 (* variants of rewrite; TODO del. put_asm *)
   262 fun rewrite_inst_  thy rew_ord rls put_asm subst thm ct =
   263   rewrite__ thy 1 subst rew_ord rls put_asm thm ct;
   264 fun rewrite_set_inst_ thy put_asm subst rls ct = rewrite__set_ thy 1 put_asm subst rls ct;
   265 
   266 (* given a list of equalities (lhs = rhs) and a term, 
   267    replace all occurrences of lhs in the term with rhs;
   268    thus the order or equalities matters: put variables in lhs first. *)
   269 fun rewrite_terms_ thy ord erls equs t =
   270   let
   271 	  fun rew_ (t', asm') [] _ = (t', asm')
   272 	    | rew_ (t', asm') (rules as r::rs) t =
   273 	        let
   274 	          val (t'', asm'', _(*lrd*), rew) = rew_sub thy 1 [] ord erls false [] (HOLogic.Trueprop $ r) t
   275 	        in 
   276 	          if rew 
   277 	          then rew_ (t'', asm' @ asm'') rules t''
   278 	          else rew_ (t', asm') rs t'
   279 	        end
   280 	  val (t'', asm'') = rew_ (TermC.empty, []) equs t
   281     in if t'' = TermC.empty then NONE else SOME (t'', asm'')
   282     end;
   283 
   284 (* search ct for adjacent numerals and calculate them by operator isa_fn *)
   285 fun calculate_ thy isa_fn ct =
   286   let val ct = TermC.uminus_to_string ct
   287     in case Eval.adhoc_thm thy isa_fn ct of
   288 	   NONE => NONE
   289 	 | SOME (thmID, thm) =>
   290 	   (let val rew = case rewrite_ thy Rewrite_Ord.dummy_ord Rule_Set.empty false thm ct of
   291          SOME (rew, _) => rew
   292        | NONE => raise ERROR ""
   293      in SOME (rew, (thmID, thm)) end)
   294 	   handle NO_REWRITE => raise ERROR ("calculate_: " ^ thmID ^ " does not rewrite")
   295   end;
   296 
   297 fun eval_prog_expr thy srls t =
   298   let val rew = rewrite_set_ thy false srls t;
   299   in case rew of SOME (res,_) => res | NONE => t end;
   300 
   301 fun eval_true_ _ _ (Const (\<^const_name>\<open>True\<close>,_)) = true
   302   | eval_true_ thy rls t =
   303     case rewrite_set_ thy false rls t of
   304 	   SOME (Const (\<^const_name>\<open>True\<close>,_),_) => true
   305 	 | _ => false;
   306 
   307 end