src/ZF/arith_data.ML
author wenzelm
Sat Nov 10 16:25:17 2001 +0100 (2001-11-10)
changeset 12134 7049eead7a50
parent 9874 0aa0874ab66b
child 12206 60d52181840c
permissions -rw-r--r--
use Tactic.prove;
     1 (*  Title:      ZF/arith_data.ML
     2     ID:         $Id$
     3     Author:     Lawrence C Paulson, Cambridge University Computer Laboratory
     4     Copyright   2000  University of Cambridge
     5 
     6 Arithmetic simplification: cancellation of common terms
     7 *)
     8 
     9 signature ARITH_DATA =
    10 sig
    11   (*the main outcome*)
    12   val nat_cancel: simproc list
    13   (*tools for use in similar applications*)
    14   val gen_trans_tac: thm -> thm option -> tactic
    15   val prove_conv: string -> tactic list -> Sign.sg -> 
    16                   thm list -> term * term -> thm option
    17   val simplify_meta_eq: thm list -> thm -> thm
    18   (*debugging*)
    19   structure EqCancelNumeralsData   : CANCEL_NUMERALS_DATA
    20   structure LessCancelNumeralsData : CANCEL_NUMERALS_DATA
    21   structure DiffCancelNumeralsData : CANCEL_NUMERALS_DATA
    22 end;
    23 
    24 
    25 structure ArithData: ARITH_DATA =
    26 struct
    27 
    28 val iT = Ind_Syntax.iT;
    29 
    30 val zero = Const("0", iT);
    31 val succ = Const("succ", iT --> iT);
    32 fun mk_succ t = succ $ t;
    33 val one = mk_succ zero;
    34 
    35 val mk_plus = FOLogic.mk_binop "Arith.add";
    36 
    37 (*Thus mk_sum[t] yields t+#0; longer sums don't have a trailing zero*)
    38 fun mk_sum []        = zero
    39   | mk_sum [t,u]     = mk_plus (t, u)
    40   | mk_sum (t :: ts) = mk_plus (t, mk_sum ts);
    41 
    42 (*this version ALWAYS includes a trailing zero*)
    43 fun long_mk_sum []        = zero
    44   | long_mk_sum (t :: ts) = mk_plus (t, mk_sum ts);
    45 
    46 val dest_plus = FOLogic.dest_bin "Arith.add" iT;
    47 
    48 (* dest_sum *)
    49 
    50 fun dest_sum (Const("0",_)) = []
    51   | dest_sum (Const("succ",_) $ t) = one :: dest_sum t
    52   | dest_sum (Const("Arith.add",_) $ t $ u) = dest_sum t @ dest_sum u
    53   | dest_sum tm = [tm];
    54 
    55 (*Apply the given rewrite (if present) just once*)
    56 fun gen_trans_tac th2 None      = all_tac
    57   | gen_trans_tac th2 (Some th) = ALLGOALS (rtac (th RS th2));
    58 
    59 (*Use <-> or = depending on the type of t*)
    60 fun mk_eq_iff(t,u) =
    61   if fastype_of t = iT then FOLogic.mk_eq(t,u)
    62                        else FOLogic.mk_iff(t,u);
    63 
    64 (*We remove equality assumptions because they confuse the simplifier and
    65   because only type-checking assumptions are necessary.*)
    66 fun is_eq_thm th = 
    67     can FOLogic.dest_eq (FOLogic.dest_Trueprop (#prop (rep_thm th)));
    68 
    69 fun add_chyps chyps ct = Drule.list_implies (map cprop_of chyps, ct);
    70 
    71 fun prove_conv name tacs sg hyps (t,u) =
    72   if t aconv u then None
    73   else
    74   let val hyps' = filter (not o is_eq_thm) hyps
    75       val goal = Logic.list_implies (map (#prop o Thm.rep_thm) hyps',
    76         FOLogic.mk_Trueprop (mk_eq_iff (t, u)));
    77   in Some (hyps' MRS Tactic.prove sg [] [] goal (K (EVERY tacs)))
    78       handle ERROR => 
    79 	(warning ("Cancellation failed: no typing information? (" ^ name ^ ")"); None)
    80   end;
    81 
    82 fun prep_simproc (name, pats, proc) = Simplifier.mk_simproc name pats proc;
    83 fun prep_pat s = Thm.read_cterm (Theory.sign_of (the_context ()))
    84                       (s, TypeInfer.anyT ["logic"]);
    85 val prep_pats = map prep_pat;
    86 
    87 
    88 (*** Use CancelNumerals simproc without binary numerals, 
    89      just for cancellation ***)
    90 
    91 val mk_times = FOLogic.mk_binop "Arith.mult";
    92 
    93 fun mk_prod [] = one
    94   | mk_prod [t] = t
    95   | mk_prod (t :: ts) = if t = one then mk_prod ts
    96                         else mk_times (t, mk_prod ts);
    97 
    98 val dest_times = FOLogic.dest_bin "Arith.mult" iT;
    99 
   100 fun dest_prod t =
   101       let val (t,u) = dest_times t
   102       in  dest_prod t @ dest_prod u  end
   103       handle TERM _ => [t];
   104 
   105 (*Dummy version: the only arguments are 0 and 1*)
   106 fun mk_coeff (0, t) = zero
   107   | mk_coeff (1, t) = t
   108   | mk_coeff _       = raise TERM("mk_coeff", []);
   109 
   110 (*Dummy version: the "coefficient" is always 1.
   111   In the result, the factors are sorted terms*)
   112 fun dest_coeff t = (1, mk_prod (sort Term.term_ord (dest_prod t)));
   113 
   114 (*Find first coefficient-term THAT MATCHES u*)
   115 fun find_first_coeff past u [] = raise TERM("find_first_coeff", [])
   116   | find_first_coeff past u (t::terms) =
   117         let val (n,u') = dest_coeff t
   118         in  if u aconv u' then (n, rev past @ terms)
   119                           else find_first_coeff (t::past) u terms
   120         end
   121         handle TERM _ => find_first_coeff (t::past) u terms;
   122 
   123 
   124 (*Simplify #1*n and n*#1 to n*)
   125 val add_0s = [add_0_natify, add_0_right_natify];
   126 val add_succs = [add_succ, add_succ_right];
   127 val mult_1s = [mult_1_natify, mult_1_right_natify];
   128 val tc_rules = [natify_in_nat, add_type, diff_type, mult_type];
   129 val natifys = [natify_0, natify_ident, add_natify1, add_natify2,
   130                diff_natify1, diff_natify2];
   131 
   132 (*Final simplification: cancel + and **)
   133 fun simplify_meta_eq rules =
   134     mk_meta_eq o
   135     simplify (FOL_ss addeqcongs[eq_cong2,iff_cong2] 
   136                      delsimps iff_simps (*these could erase the whole rule!*)
   137 		     addsimps rules);
   138 
   139 val final_rules = add_0s @ mult_1s @ [mult_0, mult_0_right];
   140 
   141 structure CancelNumeralsCommon =
   142   struct
   143   val mk_sum            = mk_sum
   144   val dest_sum          = dest_sum
   145   val mk_coeff          = mk_coeff
   146   val dest_coeff        = dest_coeff
   147   val find_first_coeff  = find_first_coeff []
   148   val norm_tac_ss1 = ZF_ss addsimps add_0s@add_succs@mult_1s@add_ac
   149   val norm_tac_ss2 = ZF_ss addsimps add_ac@mult_ac@tc_rules@natifys
   150   val norm_tac = ALLGOALS (asm_simp_tac norm_tac_ss1)
   151                  THEN ALLGOALS (asm_simp_tac norm_tac_ss2)
   152   val numeral_simp_tac_ss = ZF_ss addsimps add_0s@tc_rules@natifys
   153   val numeral_simp_tac  = ALLGOALS (asm_simp_tac numeral_simp_tac_ss)
   154   val simplify_meta_eq  = simplify_meta_eq final_rules
   155   end;
   156 
   157 (** The functor argumnets are declared as separate structures
   158     so that they can be exported to ease debugging. **)
   159 
   160 structure EqCancelNumeralsData = 
   161   struct
   162   open CancelNumeralsCommon
   163   val prove_conv = prove_conv "nateq_cancel_numerals"
   164   val mk_bal   = FOLogic.mk_eq
   165   val dest_bal = FOLogic.dest_eq
   166   val bal_add1 = eq_add_iff RS iff_trans
   167   val bal_add2 = eq_add_iff RS iff_trans
   168   val trans_tac = gen_trans_tac iff_trans
   169   end;
   170 
   171 structure EqCancelNumerals = CancelNumeralsFun(EqCancelNumeralsData);
   172 
   173 structure LessCancelNumeralsData = 
   174   struct
   175   open CancelNumeralsCommon
   176   val prove_conv = prove_conv "natless_cancel_numerals"
   177   val mk_bal   = FOLogic.mk_binrel "Ordinal.op <"
   178   val dest_bal = FOLogic.dest_bin "Ordinal.op <" iT
   179   val bal_add1 = less_add_iff RS iff_trans
   180   val bal_add2 = less_add_iff RS iff_trans
   181   val trans_tac = gen_trans_tac iff_trans
   182   end;
   183 
   184 structure LessCancelNumerals = CancelNumeralsFun(LessCancelNumeralsData);
   185 
   186 structure DiffCancelNumeralsData = 
   187   struct
   188   open CancelNumeralsCommon
   189   val prove_conv = prove_conv "natdiff_cancel_numerals"
   190   val mk_bal   = FOLogic.mk_binop "Arith.diff"
   191   val dest_bal = FOLogic.dest_bin "Arith.diff" iT
   192   val bal_add1 = diff_add_eq RS trans
   193   val bal_add2 = diff_add_eq RS trans
   194   val trans_tac = gen_trans_tac trans
   195   end;
   196 
   197 structure DiffCancelNumerals = CancelNumeralsFun(DiffCancelNumeralsData);
   198 
   199 
   200 val nat_cancel =
   201       map prep_simproc
   202        [("nateq_cancel_numerals",
   203 	 prep_pats ["l #+ m = n", "l = m #+ n",
   204 		    "l #* m = n", "l = m #* n",
   205 		    "succ(m) = n", "m = succ(n)"],
   206 	 EqCancelNumerals.proc),
   207 	("natless_cancel_numerals",
   208 	 prep_pats ["l #+ m < n", "l < m #+ n",
   209 		    "l #* m < n", "l < m #* n",
   210 		    "succ(m) < n", "m < succ(n)"],
   211 	 LessCancelNumerals.proc),
   212 	("natdiff_cancel_numerals",
   213 	 prep_pats ["(l #+ m) #- n", "l #- (m #+ n)",
   214 		    "(l #* m) #- n", "l #- (m #* n)",
   215 		    "succ(m) #- n", "m #- succ(n)"],
   216 	 DiffCancelNumerals.proc)];
   217 
   218 end;
   219 
   220 (*examples:
   221 print_depth 22;
   222 set timing;
   223 set trace_simp;
   224 fun test s = (Goal s; by (Asm_simp_tac 1));
   225 
   226 test "x #+ y = x #+ z";
   227 test "y #+ x = x #+ z";
   228 test "x #+ y #+ z = x #+ z";
   229 test "y #+ (z #+ x) = z #+ x";
   230 test "x #+ y #+ z = (z #+ y) #+ (x #+ w)";
   231 test "x#*y #+ z = (z #+ y) #+ (y#*x #+ w)";
   232 
   233 test "x #+ succ(y) = x #+ z";
   234 test "x #+ succ(y) = succ(z #+ x)";
   235 test "succ(x) #+ succ(y) #+ z = succ(z #+ y) #+ succ(x #+ w)";
   236 
   237 test "(x #+ y) #- (x #+ z) = w";
   238 test "(y #+ x) #- (x #+ z) = dd";
   239 test "(x #+ y #+ z) #- (x #+ z) = dd";
   240 test "(y #+ (z #+ x)) #- (z #+ x) = dd";
   241 test "(x #+ y #+ z) #- ((z #+ y) #+ (x #+ w)) = dd";
   242 test "(x#*y #+ z) #- ((z #+ y) #+ (y#*x #+ w)) = dd";
   243 
   244 (*BAD occurrence of natify*)
   245 test "(x #+ succ(y)) #- (x #+ z) = dd";
   246 
   247 test "x #* y2 #+ y #* x2 = y #* x2 #+ x #* y2";
   248 
   249 test "(x #+ succ(y)) #- (succ(z #+ x)) = dd";
   250 test "(succ(x) #+ succ(y) #+ z) #- (succ(z #+ y) #+ succ(x #+ w)) = dd";
   251 
   252 (*use of typing information*)
   253 test "x : nat ==> x #+ y = x";
   254 test "x : nat --> x #+ y = x";
   255 test "x : nat ==> x #+ y < x";
   256 test "x : nat ==> x < y#+x";
   257 
   258 (*fails: no typing information isn't visible*)
   259 test "x #+ y = x";
   260 
   261 test "x #+ y < x #+ z";
   262 test "y #+ x < x #+ z";
   263 test "x #+ y #+ z < x #+ z";
   264 test "y #+ z #+ x < x #+ z";
   265 test "y #+ (z #+ x) < z #+ x";
   266 test "x #+ y #+ z < (z #+ y) #+ (x #+ w)";
   267 test "x#*y #+ z < (z #+ y) #+ (y#*x #+ w)";
   268 
   269 test "x #+ succ(y) < x #+ z";
   270 test "x #+ succ(y) < succ(z #+ x)";
   271 test "succ(x) #+ succ(y) #+ z < succ(z #+ y) #+ succ(x #+ w)";
   272 
   273 test "x #+ succ(y) le succ(z #+ x)";
   274 *)