test/Tools/isac/ProgLang/rewrite.sml
author Walther Neuper <wneuper@ist.tugraz.at>
Sat, 24 Feb 2018 11:14:56 +0100
changeset 59381 cb7e75507c05
parent 59358 8509ca4e79ec
child 59382 364ce4699452
permissions -rw-r--r--
Rewrite: cleanup source file
neuper@38036
     1
(* Title: tests for ProgLang/rewrite.sml
neuper@37906
     2
   TODO.WN0509 collect typical tests from systest here !!!!!
neuper@38036
     3
   Author: Walther Neuper 050908
neuper@37906
     4
   (c) copyright due to lincense terms.
neuper@37906
     5
neuper@38022
     6
12345678901234567890123456789012345678901234567890123456789012345678901234567890
neuper@38022
     7
        10        20        30        40        50        60        70        80
neuper@37906
     8
*)
neuper@37906
     9
neuper@38022
    10
"--------------------------------------------------------";
neuper@38022
    11
"table of contents --------------------------------------";
neuper@38022
    12
"--------------------------------------------------------";
neuper@38022
    13
"----------- assemble rewrite ---------------------------";
neuper@38022
    14
"----------- test rewriting without Isac's thys ---------";
neuper@38022
    15
"----------- conditional rewriting without Isac's thys --";
neuper@38022
    16
"----------- step through 'and rew_sub': ----------------";
neuper@38025
    17
"----------- step through 'fun rewrite_terms_'  ---------";
neuper@38022
    18
"----------- rewrite_inst_ bdvs -------------------------";
neuper@38022
    19
"----------- check diff 2002--2009-3 --------------------";
neuper@38039
    20
"----------- compare all prepat's existing 2010 ---------";
neuper@38039
    21
"----------- fun app_rev, Rrls, -------------------------";
neuper@42223
    22
"----------- 2011 thms with axclasses -------------------";
neuper@42394
    23
"----------- repair NO asms from rls RatEq_eliminate ----";
wneuper@59110
    24
"----------- fun assoc_thm' -----------------------------";
wneuper@59252
    25
"----------- fun rewrite_ down to Pattern.match ------------------------------------------------";
wneuper@59381
    26
"----------- fun mk_thm ------------------------------------------------------------------------";
neuper@38022
    27
"--------------------------------------------------------";
neuper@38022
    28
"--------------------------------------------------------";
neuper@38022
    29
"--------------------------------------------------------";
neuper@37906
    30
neuper@38036
    31
neuper@38022
    32
"----------- assemble rewrite ---------------------------";
neuper@38022
    33
"----------- assemble rewrite ---------------------------";
neuper@38022
    34
"----------- assemble rewrite ---------------------------";
neuper@37906
    35
"===== rewriting by thm with 'a";
neuper@41924
    36
(*show_types := true;*)
akargl@42188
    37
neuper@37906
    38
val thy = @{theory Complex_Main};
neuper@37906
    39
val ctxt = @{context};
wneuper@59112
    40
val thm = @{thm add.commute};
wneuper@59188
    41
val t = (Thm.term_of o the) (parse thy "((r + u) + t) + s");
neuper@37906
    42
"----- from old: fun rewrite__";
neuper@37906
    43
val bdv = [];
wneuper@59188
    44
val r = (((inst_bdv bdv) o norm o #prop o Thm.rep_thm) thm);
neuper@37906
    45
"----- from old: and rew_sub";
neuper@41942
    46
val (LHS,RHS) = (dest_equals' o strip_trueprop 
neuper@37906
    47
   	      o Logic.strip_imp_concl) r;
neuper@37906
    48
(* old
neuper@41942
    49
val insts = Pattern.match thy (LHS,t) (Vartab.empty, Vartab.empty);*)
neuper@37906
    50
"----- fun match_rew in Pure/pattern.ML";
neuper@41942
    51
val rtm = the_default RHS (Term.rename_abs LHS t RHS);
neuper@37906
    52
neuper@38022
    53
writeln(Syntax.string_of_term ctxt rtm);
neuper@41942
    54
writeln(Syntax.string_of_term ctxt LHS);
neuper@38022
    55
writeln(Syntax.string_of_term ctxt t);
neuper@37906
    56
neuper@41942
    57
(Pattern.match thy (LHS, t) (Vartab.empty, Vartab.empty));
neuper@41942
    58
val (rew, RHS) = (Envir.subst_term 
neuper@41942
    59
  (Pattern.match thy (LHS, t) (Vartab.empty, Vartab.empty)) rtm, rtm);
neuper@37906
    60
(*lookup in isabelle?trace?response...*)
neuper@37906
    61
writeln(Syntax.string_of_term ctxt rew);
neuper@41942
    62
writeln(Syntax.string_of_term ctxt RHS);
neuper@37906
    63
"===== rewriting: prep insertion into rew_sub";
neuper@37906
    64
val thy = @{theory Complex_Main};
neuper@37906
    65
val ctxt = @{context};
wneuper@59358
    66
val thm =  @{thm nonzero_divide_mult_cancel_right};
neuper@37906
    67
val r = Thm.prop_of thm;
wneuper@59358
    68
val tm = @{term "x / (2 * x)::real"};
neuper@37906
    69
"----- and rew_sub";
neuper@41942
    70
val (LHS, RHS) = (HOLogic.dest_eq o HOLogic.dest_Trueprop
neuper@37906
    71
                  o Logic.strip_imp_concl) r;
neuper@41942
    72
val r' = Envir.subst_term (Pattern.match thy (LHS, tm) 
neuper@37906
    73
                                (Vartab.empty, Vartab.empty)) r;
neuper@37906
    74
val p' = (fst o Logic.strip_prems) (Logic.count_prems r', [], r');
neuper@37906
    75
val t' = (snd o HOLogic.dest_eq o HOLogic.dest_Trueprop 
neuper@37906
    76
            o Logic.strip_imp_concl) r';
neuper@37906
    77
neuper@37906
    78
(*is displayed on top of <response> buffer...*)
neuper@48761
    79
Pretty.writeln (Proof_Context.pretty_term_abbrev @{context} r');
neuper@48761
    80
Pretty.writeln (Proof_Context.pretty_term_abbrev @{context} t');
akargl@42188
    81
(**)
neuper@37906
    82
neuper@38022
    83
"----------- test rewriting without Isac's thys ---------";
neuper@38022
    84
"----------- test rewriting without Isac's thys ---------";
neuper@38022
    85
"----------- test rewriting without Isac's thys ---------";
neuper@38022
    86
neuper@37906
    87
"===== rewriting with Isabelle2009-1 only, i.e without isac-hacks";
neuper@37906
    88
val thy = @{theory Complex_Main};
neuper@37906
    89
val ctxt = @{context};
wneuper@59112
    90
val thm =  @{thm add.commute};
neuper@37906
    91
val tm = @{term "x + y*z::real"};
neuper@37906
    92
neuper@37906
    93
val SOME (r,_) = (rewrite_ thy dummy_ord e_rls false thm tm)
neuper@38022
    94
  handle _ => error "rewrite.sml diff.behav. in rewriting";
neuper@37906
    95
(*is displayed on _TOP_ of <response> buffer...*)
neuper@48761
    96
Pretty.writeln (Proof_Context.pretty_term_abbrev @{context} r);
neuper@38022
    97
if r = @{term "y*z + x::real"}
neuper@38022
    98
then () else error "rewrite.sml diff.result in rewriting";
neuper@37906
    99
neuper@37906
   100
"----- rewriting a subterm";
neuper@37906
   101
val tm = @{term "w*(x + y*z)::real"};
neuper@37906
   102
neuper@37906
   103
val SOME (r,_) = (rewrite_ thy dummy_ord e_rls false thm tm)
neuper@38022
   104
  handle _ => error "rewrite.sml diff.behav. in rew_sub";
neuper@37906
   105
neuper@37906
   106
"----- ordered rewriting";
neuper@38022
   107
fun tord (_:subst) pp = Term_Ord.termless pp;
neuper@37906
   108
if tord [] (@{term "x + y*z::real"}, @{term "y*z + x::real"}) then ()
neuper@38022
   109
else error "rewrite.sml diff.behav. in ord.rewr.";
neuper@37906
   110
neuper@37906
   111
val NONE = (rewrite_ thy tord e_rls false thm tm)
neuper@38022
   112
  handle _ => error "rewrite.sml diff.behav. in rewriting";
neuper@37906
   113
(*is displayed on _TOP_ of <response> buffer...*)
neuper@48761
   114
Pretty.writeln (Proof_Context.pretty_term_abbrev @{context} r);
neuper@37906
   115
neuper@37906
   116
val tm = @{term "x*y + z::real"};
neuper@37906
   117
val SOME (r,_) = (rewrite_ thy tord e_rls false thm tm)
neuper@38022
   118
  handle _ => error "rewrite.sml diff.behav. in rewriting";
neuper@37906
   119
neuper@37906
   120
neuper@38022
   121
"----------- conditional rewriting without Isac's thys --";
neuper@38022
   122
"----------- conditional rewriting without Isac's thys --";
neuper@38022
   123
"----------- conditional rewriting without Isac's thys --";
akargl@42224
   124
neuper@37906
   125
"===== prepr cond.rew. with Pattern.match";
neuper@37906
   126
val thy = @{theory Complex_Main};
neuper@37906
   127
val ctxt = @{context};
wneuper@59358
   128
val thm =  @{thm nonzero_divide_mult_cancel_right};
neuper@37906
   129
val rule = Thm.prop_of thm;
wneuper@59358
   130
val tm = @{term "x / (2 * x)::real"};
neuper@37906
   131
neuper@37906
   132
val prem = Logic.strip_imp_prems rule;
neuper@37906
   133
val nps = Logic.count_prems rule;
neuper@37906
   134
val prems = Logic.strip_prems (nps, [], rule);
neuper@37906
   135
neuper@37906
   136
val eq = Logic.strip_imp_concl rule;
neuper@41942
   137
val (LHS, RHS) = (HOLogic.dest_eq o HOLogic.dest_Trueprop) eq;
neuper@37906
   138
neuper@41942
   139
val mtcs = Pattern.match thy (LHS, tm) (Vartab.empty, Vartab.empty);
neuper@37906
   140
val rule' = Envir.subst_term mtcs rule;
neuper@37906
   141
neuper@37906
   142
val prems' = (fst o Logic.strip_prems) 
neuper@37906
   143
              (Logic.count_prems rule', [], rule');
neuper@41942
   144
val RHS' = (snd o HOLogic.dest_eq o HOLogic.dest_Trueprop 
neuper@37906
   145
            o Logic.strip_imp_concl) rule';
neuper@37906
   146
neuper@37906
   147
"----- conditional rewriting creating an assumption";
neuper@37906
   148
"----- conditional rewriting creating an assumption";
wneuper@59358
   149
val thm =  @{thm nonzero_divide_mult_cancel_right};
wneuper@59358
   150
val tm = @{term "x / (2 * x)::real"};
neuper@38022
   151
val SOME (rew, asm) = (rewrite_ thy dummy_ord e_rls false thm tm)
neuper@38022
   152
  handle _ => error "rewrite.sml diff.behav. in cond.rew.";
neuper@37906
   153
wneuper@59358
   154
if rew = @{term "1 / 2::real"} then () (* the rewrite is _NO_ Trueprop *)
neuper@38022
   155
else error "rewrite.sml diff.result in cond.rew.";
neuper@37906
   156
wneuper@59358
   157
if hd asm = @{term "x \<noteq> (0::real)"} (* dropped Trueprop $ ...*)
neuper@38022
   158
then () else error "rewrite.sml diff.asm in cond.rew.";
neuper@38022
   159
"----- conditional rewriting immediately: can only be done with ";
neuper@38022
   160
"------Isabelle numerals, because erls cannot handle them yet.";
neuper@37906
   161
neuper@37906
   162
neuper@38022
   163
"----------- step through 'and rew_sub': ----------------";
neuper@38022
   164
"----------- step through 'and rew_sub': ----------------";
neuper@38022
   165
"----------- step through 'and rew_sub': ----------------";
neuper@38022
   166
(*and make asms without Trueprop, beginning with the result:*)
wneuper@59358
   167
val tm = @{term "x / (2 * x)::real"};
wneuper@59188
   168
val (t', asm, _, _) = rew_sub thy 0 [] dummy_ord e_rls true [] (Thm.prop_of thm) tm;
neuper@41924
   169
(*show_types := false;*)
neuper@38022
   170
"----- evaluate arguments";
neuper@38022
   171
val (thy, i, bdv, tless, rls, put_asm, lrd, r, t) = 
wneuper@59188
   172
    (thy, 0, [], dummy_ord, e_rls, true, [], (Thm.prop_of thm), tm);
neuper@41942
   173
"----- step 1: LHS, RHS of rule";
neuper@41942
   174
     val (LHS, RHS) = (HOLogic.dest_eq o HOLogic.dest_Trueprop
neuper@38022
   175
                       o Logic.strip_imp_concl) r;
wneuper@59358
   176
term2str r = "?b \<noteq> (0::?'a) \<Longrightarrow> ?b / (?a * ?b) = (1::?'a) / ?a";
wneuper@59358
   177
term2str LHS = "?b / (?a * ?b)"; term2str RHS = "(1::?'a) / ?a";
neuper@38022
   178
"----- step 2: the rule instantiated";
neuper@38025
   179
     val r' = Envir.subst_term 
neuper@41942
   180
                  (Pattern.match thy (LHS, t) (Vartab.empty, Vartab.empty)) r;
wneuper@59358
   181
term2str r' = "x \<noteq> 0 \<Longrightarrow> x / (2 * x) = 1 / 2";
neuper@38022
   182
"----- step 3: get the (instantiated) assumption(s)";
neuper@38022
   183
     val p' = (fst o Logic.strip_prems) (Logic.count_prems r', [], r');
wneuper@59358
   184
term2str (hd p') = "x \<noteq> 0";
neuper@38022
   185
"=====vvv make asms without Trueprop ---vvv";
neuper@38022
   186
     val p' = map HOLogic.dest_Trueprop ((fst o Logic.strip_prems) 
neuper@38022
   187
                                             (Logic.count_prems r', [], r'));
neuper@38022
   188
case p' of
wneuper@59358
   189
    [Const ("HOL.Not", _) $ (Const ("HOL.eq", _) 
wneuper@59358
   190
      $ Free ("x", _) $ Const ("Groups.zero_class.zero", _))] => ()
neuper@38022
   191
  | _ => error "rewrite.sml assumption changed";
neuper@38022
   192
"=====^^^ make asms without Trueprop ---^^^";
neuper@41942
   193
"----- step 4: get the (instantiated) RHS";
neuper@38022
   194
     val t' = (snd o HOLogic.dest_eq o HOLogic.dest_Trueprop 
neuper@38022
   195
               o Logic.strip_imp_concl) r';
wneuper@59358
   196
term2str t' = "1 / 2";
neuper@37906
   197
neuper@38025
   198
"----------- step through 'fun rewrite_terms_'  ---------";
neuper@38025
   199
"----------- step through 'fun rewrite_terms_'  ---------";
neuper@38025
   200
"----------- step through 'fun rewrite_terms_'  ---------";
neuper@38025
   201
"----- step 0: args for rewrite_terms_, local fun";
neuper@38025
   202
val (thy, ord, erls, equs, t) =
neuper@41928
   203
    (@{theory "Biegelinie"}, dummy_ord, Erls, [str2term "x = 0"],
neuper@38025
   204
     str2term "M_b x = -1 * q_0 * x ^^^ 2 / 2 + x * c + c_2");
neuper@38025
   205
"----- step 1: args for rew_";
neuper@38025
   206
val ((t', asm'), (rules as r::rs), t) = ((e_term, []), equs, t);
neuper@38025
   207
"----- step 2: rew_sub";
neuper@38025
   208
rew_sub thy 1 [] ord erls false [] (Trueprop $ r) t;
neuper@38025
   209
"----- step 3: step through rew_sub -- inefficient: goes into subterms";
neuper@38025
   210
neuper@38025
   211
neuper@38025
   212
val SOME (t', _) = rewrite_terms_ thy dummy_ord Erls equs t;
neuper@38025
   213
writeln "----------- rewrite_terms_  1---------------------------";
neuper@37906
   214
if term2str t' = "M_b 0 = -1 * q_0 * 0 ^^^ 2 / 2 + 0 * c + c_2" then ()
neuper@38031
   215
else error "rewrite.sml rewrite_terms_ [x = 0]";
neuper@37906
   216
neuper@38025
   217
val equs = [str2term "M_b 0 = 0"];
neuper@38025
   218
val t = str2term "M_b 0 = -1 * q_0 * 0 ^^^ 2 / 2 + 0 * c + c_2";
neuper@38025
   219
val SOME (t', _) = rewrite_terms_ thy dummy_ord Erls equs t;
neuper@38025
   220
writeln "----------- rewrite_terms_  2---------------------------";
neuper@37906
   221
if term2str t' = "0 = -1 * q_0 * 0 ^^^ 2 / 2 + 0 * c + c_2" then ()
neuper@38031
   222
else error "rewrite.sml rewrite_terms_ [M_b 0 = 0]";
neuper@37906
   223
neuper@38025
   224
val equs = [str2term "x = 0", str2term"M_b 0 = 0"];
neuper@38025
   225
val t = str2term "M_b x = -1 * q_0 * x ^^^ 2 / 2 + x * c + c_2";
neuper@38025
   226
val SOME (t', _) = rewrite_terms_ thy dummy_ord Erls equs t;
neuper@38025
   227
writeln "----------- rewrite_terms_  3---------------------------";
neuper@37906
   228
if term2str t' = "0 = -1 * q_0 * 0 ^^^ 2 / 2 + 0 * c + c_2" then ()
neuper@38031
   229
else error "rewrite.sml rewrite_terms_ [x = 0, M_b 0 = 0]";
neuper@37906
   230
neuper@37906
   231
neuper@38022
   232
"----------- rewrite_inst_ bdvs -------------------------";
neuper@38022
   233
"----------- rewrite_inst_ bdvs -------------------------";
neuper@38022
   234
"----------- rewrite_inst_ bdvs -------------------------";
neuper@37906
   235
(*see smltest/Scripts/term_G.sml: inst_bdv 2*)
neuper@37906
   236
val t = str2term"-1 * (q_0 * L ^^^ 2) / 2 + (L * c_3 + c_4) = 0";
neuper@37906
   237
val bdvs = [(str2term"bdv_1",str2term"c"),
neuper@37906
   238
	    (str2term"bdv_2",str2term"c_2"),
neuper@37906
   239
	    (str2term"bdv_3",str2term"c_3"),
neuper@37906
   240
	    (str2term"bdv_4",str2term"c_4")];
neuper@37906
   241
(*------------ outcommented WN071210, after inclusion into ROOT.ML 
neuper@37926
   242
val SOME (t,_) = 
neuper@37906
   243
    rewrite_inst_ thy e_rew_ord 
neuper@37906
   244
		  (append_rls "erls_isolate_bdvs" e_rls 
neuper@37906
   245
			      [(Calc ("EqSystem.occur'_exactly'_in", 
neuper@37906
   246
				      eval_occur_exactly_in 
neuper@37906
   247
					  "#eval_occur_exactly_in_"))
neuper@37906
   248
			       ]) 
neuper@37967
   249
		  false bdvs (num_str @{separate_bdvs_add) t;
neuper@37906
   250
(writeln o term2str) t;
neuper@37906
   251
if term2str t = "L * c_3 + c_4 = 0 + -1 * (-1 * (q_0 * L ^^^ 2) / 2)"
neuper@38031
   252
then () else error "rewrite.sml rewrite_inst_ bdvs";
neuper@37906
   253
trace_rewrite:=true;
neuper@37906
   254
trace_rewrite:=false;--------------------------------------------*)
neuper@37906
   255
neuper@38025
   256
neuper@38022
   257
"----------- check diff 2002--2009-3 --------------------";
neuper@38022
   258
"----------- check diff 2002--2009-3 --------------------";
neuper@38022
   259
"----------- check diff 2002--2009-3 --------------------";
neuper@38022
   260
(*----- 2002 -------------------------------------------------------------------
neuper@38022
   261
#  rls: norm_Rational_noadd_fractions on: 1 / EI * (L * q_0 * x / 2 + -1 *
neuper@38022
   262
q_0 * x ^^^ 2 / 2)
neuper@38022
   263
##  rls: discard_minus_ on: 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2
neuper@38022
   264
/ 2)
neuper@38022
   265
###  try thm: real_diff_minus
neuper@38022
   266
###  try thm: sym_real_mult_minus1
neuper@38022
   267
##  rls: rat_mult_poly on: 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2
neuper@38022
   268
/ 2)
neuper@38022
   269
###  try thm: rat_mult_poly_l
neuper@38022
   270
###  try thm: rat_mult_poly_r
neuper@38022
   271
####  eval asms: L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2 is_polyexp
neuper@38022
   272
==> 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2) =
neuper@38022
   273
    1 * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2) / EI
neuper@38022
   274
#####  rls: e_rls-is_polyexp on: L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2
neuper@38022
   275
is_polyexp
neuper@38022
   276
######  try calc: Poly.is'_polyexp'
neuper@38022
   277
======  calc. to: False
neuper@38022
   278
######  try calc: Poly.is'_polyexp'
neuper@38022
   279
######  try calc: Poly.is'_polyexp'
neuper@38022
   280
####  asms false: ["L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2 is_polyexp"]
neuper@38022
   281
----- 2002 NONE rewrite --------------------------------------------------------
neuper@38022
   282
----- 2009 should maintain this behaviour, but: --------------------------------
neuper@38022
   283
#  rls: norm_Rational_noadd_fractions on: 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2)
neuper@38022
   284
##  rls: discard_minus on: 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2)
neuper@38022
   285
###  try thm: real_diff_minus
neuper@38022
   286
###  try thm: sym_real_mult_minus1
neuper@38022
   287
##  rls: rat_mult_poly on: 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2)
neuper@38022
   288
###  try thm: rat_mult_poly_l
neuper@38022
   289
###  try thm: rat_mult_poly_r
neuper@38022
   290
####  eval asms: L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2 is_polyexp
neuper@38022
   291
==> 1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2) =
neuper@38022
   292
    1 * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2) / EI
neuper@38022
   293
#####  rls: e_rls-is_polyexp on: L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2 is_polyexp
neuper@38022
   294
######  try calc: Poly.is'_polyexp'
neuper@38022
   295
======  calc. to: False
neuper@38022
   296
######  try calc: Poly.is'_polyexp'
neuper@38022
   297
######  try calc: Poly.is'_polyexp'
neuper@38022
   298
####  asms accepted: ["L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2 is_polyexp"]   stored: ["False"]
neuper@38022
   299
===  rewrites to: 1 * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2) / EI
neuper@38022
   300
----- 2009 -------------------------------------------------------------------*)
neuper@38022
   301
neuper@38022
   302
(*the situation as was before repair (asm without Trueprop) is outcommented*)
neuper@41928
   303
val thy = @{theory "Isac"};
neuper@38022
   304
"===== example which raised the problem =================";
neuper@38022
   305
val t = @{term "1 / EI * (L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2)"};
neuper@38022
   306
"----- rewrite_set_inst_ norm_Rational_noadd_fractions--";
neuper@38022
   307
val subs = [(str2term "bdv", str2term "x")];
neuper@38022
   308
val rls = norm_Rational_noadd_fractions;
wneuper@59112
   309
val SOME (t', asms) = rewrite_set_inst_ thy true subs rls t;
wneuper@59112
   310
if term2str t' = "1 / EI * (L * q_0 * x / 2 + -1 * 1 * q_0 * x ^^^ 2 / 2)" andalso
wneuper@59112
   311
  terms2str asms = "[]" then {}
wneuper@59112
   312
else error "this was NONE with Isabelle2013-2 ?!?"
neuper@38022
   313
"----- rewrite_ rat_mult_poly_r--------------------------";
neuper@38022
   314
val thm = @{thm rat_mult_poly_r};
neuper@38022
   315
         "?c::real is_polyexp ==> (?a::real) / (?b::real) * ?c = ?a * ?c / ?b";
neuper@38022
   316
val erls = append_rls "e_rls-is_polyexp" e_rls 
neuper@38022
   317
                      [Calc ("Poly.is'_polyexp", eval_is_polyexp "")];
neuper@38022
   318
val NONE (*SOME (t'', _)*) = rewrite_ thy dummy_ord erls true thm t;
neuper@38022
   319
(*t' = t''; (*false because of further rewrites in t'*)*)
neuper@38022
   320
"----- rew_sub  --------------------------------";
wneuper@59188
   321
val (t''', _, _, _) = rew_sub thy 0 [] dummy_ord erls true [] (Thm.prop_of thm) t;
neuper@38022
   322
(*t'' = t'''; (*true*)*)
neuper@38022
   323
"----- rewrite_set_ erls --------------------------------";
neuper@38022
   324
val cond = @{term "(L * q_0 * x / 2 + -1 * q_0 * x ^^^ 2 / 2)"};
neuper@38022
   325
val NONE = rewrite_set_ thy true erls cond; 
neuper@38022
   326
(* ^^^^^ goes with '======  calc. to: False' above from beginning*)
neuper@38022
   327
neuper@38022
   328
writeln "===== maximally simplified example =====================";
neuper@38022
   329
val t = @{term "a / b * (x / x ^^^ 2)"};
neuper@38022
   330
         "?c::real is_polyexp ==> (?a::real) / (?b::real) * ?c = ?a * ?c / ?b";
neuper@38022
   331
writeln "----- rewrite_set_inst_ norm_Rational_noadd_fractions--";
neuper@38022
   332
val (*NONE*) SOME (t', asm) = rewrite_set_inst_ thy true subs rls t;
neuper@38022
   333
term2str t' = "a * x / (b * x ^^^ 2)"; (*rew. by thm outside test case*)
neuper@38022
   334
(*checked visually: trace_rewrite looks like above for 2009*)
neuper@38022
   335
neuper@38022
   336
writeln "----- rewrite_ rat_mult_poly_r--------------------------";
neuper@38022
   337
val NONE (*SOME (t'', _)*) = rewrite_ thy dummy_ord erls true thm t;
neuper@38022
   338
(*t' = t''; (*false because of further rewrites in t'*)*)
neuper@38022
   339
writeln "----- rew_sub  --------------------------------";
wneuper@59188
   340
val (t''', _, _, _) = rew_sub thy 0 [] dummy_ord erls true [] (Thm.prop_of thm) t;
neuper@38022
   341
(*t'' = t'''; (*true*)*)
neuper@38022
   342
writeln "----- rewrite_set_ erls --------------------------------";
neuper@38022
   343
val cond = @{term "(x / x ^^^ 2)"};
neuper@38022
   344
val NONE = rewrite_set_ thy true erls cond; 
neuper@38022
   345
(* ^^^^^ goes with '======  calc. to: False' above from beginning*)
neuper@42201
   346
neuper@38025
   347
neuper@38039
   348
"----------- compare all prepat's existing 2010 ---------";
neuper@38039
   349
"----------- compare all prepat's existing 2010 ---------";
neuper@38039
   350
"----------- compare all prepat's existing 2010 ---------";
neuper@41928
   351
val thy = @{theory "Isac"};
neuper@38036
   352
val t = @{term "a + b * x = (0 ::real)"};
neuper@38036
   353
val pat = parse_patt thy "?l = (?r ::real)";
neuper@38036
   354
val precond = parse_patt thy "is_polynomial (?l::real)";(*no infix def*)
neuper@38039
   355
val precond = parse_patt thy "(?l::real) is_expanded"; 
neuper@38036
   356
neuper@38036
   357
val inst = Pattern.match thy (pat, t) (Vartab.empty, Vartab.empty);
neuper@38036
   358
val preinst = Envir.subst_term inst precond;
neuper@38036
   359
term2str preinst;
neuper@38036
   360
neuper@38036
   361
"===== Rational.thy: cancel ===";
neuper@38036
   362
(* pat matched with the current term gives an environment 
neuper@38036
   363
   (or not: hen the Rrls not applied);
neuper@48760
   364
   if pre1 and pre2 evaluate to @{term True} in this environment,
neuper@38036
   365
   then the Rrls is applied. *)
neuper@38036
   366
val t = str2term "(a + b) / c ::real";
neuper@38036
   367
val pat = parse_patt thy "?r / ?s ::real";
neuper@38036
   368
val pres = [parse_patt thy "?r is_expanded", parse_patt thy "?s is_expanded"];
neuper@38036
   369
val prepat = [(pres, pat)];
neuper@38036
   370
val erls = rational_erls;
neuper@38036
   371
(* erls got from Rrls {erls, prepat, scr = Rfuns {normal_form, ...}, ...} *)
neuper@38036
   372
neuper@38036
   373
val subst = Pattern.match thy (pat, t) (Vartab.empty, Vartab.empty);
neuper@38036
   374
val asms = map (Envir.subst_term subst) pres;
neuper@38036
   375
if terms2str asms = "[\"a + b is_expanded\",\"c is_expanded\"]"
neuper@38036
   376
then () else error "rewrite.sml: prepat cancel subst";
neuper@38036
   377
if ([], true) = eval__true thy 0 asms [] erls
neuper@38036
   378
then () else error "rewrite.sml: prepat cancel eval__true";
neuper@38036
   379
neuper@52105
   380
"===== Rational.thy: add_fractions_p ===";
neuper@38036
   381
(* if each pat* matches with the current term, the Rrls is applied
neuper@48760
   382
   (there are no preconditions to be checked, they are @{term True}) *)
neuper@38036
   383
val t = str2term "a / b + 1 / 2";
neuper@38036
   384
val pat = parse_patt thy "?r / ?s + ?u / ?v";
neuper@48760
   385
val pres = [@{term True}];
neuper@38036
   386
val prepat = [(pres, pat)];
neuper@38036
   387
val erls = rational_erls;
neuper@38036
   388
(* erls got from Rrls {erls, prepat, scr = Rfuns {normal_form, ...}, ...} *)
neuper@38036
   389
neuper@38036
   390
val subst = Pattern.match thy (pat, t) (Vartab.empty, Vartab.empty);
neuper@38036
   391
if ([], true) = eval__true thy 0 (map (Envir.subst_term subst) pres) [] erls
neuper@52105
   392
then () else error "rewrite.sml: prepat add_fractions_p";
neuper@38036
   393
neuper@38036
   394
"===== Poly.thy: order_mult_ ===";
neuper@38036
   395
          (* ?p matched with the current term gives an environment,
neuper@38036
   396
             which evaluates (the instantiated) "p is_multUnordered" to true*)
neuper@38036
   397
val t = str2term "x^^^2 * x";
neuper@38036
   398
val pat = parse_patt thy "?p :: real"
neuper@38036
   399
val pres = [parse_patt thy "?p is_multUnordered"];
neuper@38036
   400
val prepat = [(pres, pat)];
neuper@38036
   401
val erls = append_rls "e_rls-is_multUnordered" e_rls
neuper@38036
   402
		      [Calc ("Poly.is'_multUnordered", 
neuper@38036
   403
                             eval_is_multUnordered "")];
neuper@38036
   404
neuper@38036
   405
val subst = Pattern.match thy (pat, t) (Vartab.empty, Vartab.empty);
neuper@38036
   406
val asms = map (Envir.subst_term subst) pres;
neuper@38036
   407
if terms2str asms = "[\"x ^^^ 2 * x is_multUnordered\"]"
neuper@38036
   408
then () else error "rewrite.sml: prepat order_mult_ subst";
neuper@38036
   409
if ([], true) = eval__true thy 0 asms [] erls
neuper@38036
   410
then () else error "rewrite.sml: prepat order_mult_ eval__true";
neuper@38036
   411
neuper@38036
   412
neuper@38039
   413
"----------- fun app_rev, Rrls, -------------------------";
neuper@38039
   414
"----------- fun app_rev, Rrls, -------------------------";
neuper@38039
   415
"----------- fun app_rev, Rrls, -------------------------";
neuper@38039
   416
val t = str2term "x^^^2 * x";
neuper@41928
   417
neuper@38039
   418
if is_multUnordered t then () else error "rewrite.sml diff. is_multUnordered 2";
neuper@38039
   419
val tm = str2term "(x^^^2 * x) is_multUnordered";
neuper@41928
   420
eval_is_multUnordered "testid" "" tm thy;
neuper@41928
   421
neuper@38039
   422
case eval_is_multUnordered "testid" "" tm thy of
neuper@41924
   423
    SOME (_, Const ("HOL.Trueprop", _) $ 
neuper@41922
   424
                   (Const ("HOL.eq", _) $
neuper@38039
   425
                          (Const ("Poly.is'_multUnordered", _) $ _) $ 
neuper@41928
   426
                          Const ("HOL.True", _))) => ()
neuper@41928
   427
  | _ => error "rewrite.sml diff. eval_is_multUnordered 2b";
neuper@38039
   428
neuper@38039
   429
tracing "----- begin rewrite x^^^2 * x ---"; trace_rewrite := true;
neuper@38039
   430
val SOME (t', _) = rewrite_set_ thy true order_mult_ t;
neuper@38039
   431
tracing "----- end rewrite x^^^2 * x ---"; trace_rewrite := false;
neuper@38039
   432
if term2str t' = "x * x ^^^ 2" then ()
neuper@38039
   433
else error "rewrite.sml Poly.is'_multUnordered doesn't work";
neuper@38039
   434
neuper@38039
   435
(* for achieving the previous result, the following code was taken apart *)
neuper@38039
   436
"----- rewrite__set_ ---";
neuper@38039
   437
val (thy, i, _, _, (rrls as Rrls _), t) = (thy, 0, true, [], order_mult_, tm);
neuper@38039
   438
	val (t', asm, rew) = app_rev thy (i+1) rrls t;
neuper@38039
   439
"----- app_rev ---";
neuper@38039
   440
val (thy, i, rrls, t) = (thy, (i+1), rrls, t);
neuper@38039
   441
	fun chk_prepat thy erls [] t = true
neuper@38039
   442
	  | chk_prepat thy erls prepat t =
neuper@38039
   443
	    let fun chk (pres, pat) =
neuper@38039
   444
		    (let val subst: Type.tyenv * Envir.tenv = 
neuper@38039
   445
			     Pattern.match thy (pat, t)
neuper@38039
   446
					    (Vartab.empty, Vartab.empty)
neuper@38039
   447
		     in snd (eval__true thy (i+1) 
neuper@38039
   448
					(map (Envir.subst_term subst) pres)
neuper@38039
   449
					[] erls)
neuper@38039
   450
		     end)
neuper@38039
   451
		    handle _ => false
neuper@38039
   452
		fun scan_ f [] = false (*scan_ NEVER called by []*)
neuper@38039
   453
		  | scan_ f (pp::pps) = if f pp then true
neuper@38039
   454
					else scan_ f pps;
neuper@38039
   455
	    in scan_ chk prepat end;
neuper@38039
   456
neuper@38039
   457
	(*.apply the normal_form of a rev-set.*)
neuper@38039
   458
	fun app_rev' thy (Rrls{erls,prepat,scr=Rfuns{normal_form,...},...}) t =
neuper@38039
   459
	    if chk_prepat thy erls prepat t
neuper@38039
   460
	    then ((*tracing("### app_rev': t = "^(term2str t));*)
neuper@38039
   461
                  normal_form t)
neuper@38039
   462
	    else NONE;
neuper@38039
   463
(*fixme val NONE = app_rev' thy rrls t;*)
neuper@38039
   464
"----- app_rev' ---";
neuper@38039
   465
val (thy, Rrls{erls,prepat,scr=Rfuns{normal_form,...},...}, t) = (thy, rrls, t);
neuper@38039
   466
(*fixme false*)   chk_prepat thy erls prepat t;
neuper@38039
   467
"----- chk_prepat ---";
neuper@38039
   468
val (thy, erls, prepat, t) = (thy, erls, prepat, t);
neuper@38039
   469
                fun chk (pres, pat) =
neuper@38039
   470
		    (let val subst: Type.tyenv * Envir.tenv = 
neuper@38039
   471
			     Pattern.match thy (pat, t)
neuper@38039
   472
					    (Vartab.empty, Vartab.empty)
neuper@38039
   473
		     in snd (eval__true thy (i+1) 
neuper@38039
   474
					(map (Envir.subst_term subst) pres)
neuper@38039
   475
					[] erls)
neuper@38039
   476
		     end)
neuper@38039
   477
		    handle _ => false;
neuper@38039
   478
		fun scan_ f [] = false (*scan_ NEVER called by []*)
neuper@38039
   479
		  | scan_ f (pp::pps) = if f pp then true
neuper@38039
   480
					else scan_ f pps;
neuper@38039
   481
tracing "=== poly.sml: scan_ chk prepat begin";
neuper@38039
   482
scan_ chk prepat;
neuper@38039
   483
tracing "=== poly.sml: scan_ chk prepat end";
neuper@38039
   484
neuper@38039
   485
"----- chk ---";
neuper@38039
   486
(*reestablish...*) val t = str2term "x ^^^ 2 * x";
neuper@38039
   487
val [(pres, pat)] = prepat;
neuper@38039
   488
                         val subst: Type.tyenv * Envir.tenv = 
neuper@38039
   489
			     Pattern.match thy (pat, t)
neuper@38039
   490
					    (Vartab.empty, Vartab.empty);
neuper@38039
   491
neuper@38039
   492
(*fixme: asms = ["p is_multUnordered"]...instantiate*)
neuper@38039
   493
"----- eval__true ---";
neuper@38039
   494
val asms = (map (Envir.subst_term subst) pres);
neuper@38039
   495
if terms2str asms = "[\"x ^^^ 2 * x is_multUnordered\"]" then ()
neuper@38039
   496
else error "rewrite.sml: diff. is_multUnordered, asms";
neuper@38039
   497
val (thy, i, asms, bdv, rls) = 
neuper@38039
   498
    (thy, (i+1), [str2term "(x^^^2 * x) is_multUnordered"], 
neuper@38039
   499
     [] : (term * term) list, erls);
neuper@38039
   500
case eval__true thy i asms bdv rls of 
neuper@38039
   501
    ([], true) => ()
neuper@38039
   502
  | _ => error "rewrite.sml: diff. is_multUnordered, eval__true";
neuper@42223
   503
neuper@42223
   504
"----------- 2011 thms with axclasses -------------------";
neuper@42223
   505
"----------- 2011 thms with axclasses -------------------";
neuper@42223
   506
"----------- 2011 thms with axclasses -------------------";
wneuper@59348
   507
val thm = num_str @{thm div_by_1};
wneuper@59188
   508
val prop = Thm.prop_of thm;
neuper@42223
   509
atomty prop;                                     (*Type 'a*)
neuper@42223
   510
val t = str2term "(2*x)/1";                      (*Type real*)
neuper@42223
   511
neuper@42223
   512
val (thy, ro, er, inst) = 
neuper@42223
   513
    (@{theory "Isac"}, tless_true, eval_rls, [(@{term "bdv::real"}, @{term "x::real"})]);
neuper@42223
   514
val SOME (t, _) = rewrite_ thy ro er true thm t; (*real matches 'a ?via ring? etc*)
neuper@42223
   515
neuper@42394
   516
"----------- repair NO asms from rls RatEq_eliminate ----";
neuper@42394
   517
"----------- repair NO asms from rls RatEq_eliminate ----";
neuper@42394
   518
"----------- repair NO asms from rls RatEq_eliminate ----";
neuper@42394
   519
val t = str2term "1 / x = 5";
neuper@42394
   520
trace_rewrite := true;
neuper@42394
   521
val SOME (t', asm) = rewrite_ thy e_rew_ord e_rls true @{thm rat_mult_denominator_right} t;
neuper@42394
   522
term2str t' = "1 = 5 * x";
neuper@42394
   523
terms2str asm = "[\"x ~= 0\"]";
neuper@42394
   524
(*
neuper@42394
   525
 ##  eval asms: x ~= 0 ==> (1 / x = 5) = (1 = 5 * x) 
neuper@42394
   526
 ###  rls: e_rls on: x ~= 0 
neuper@42394
   527
 ##  asms accepted: ["x ~= 0"]   stored: ["x ~= 0"] 
neuper@42394
   528
*)
neuper@42394
   529
trace_rewrite := false;
neuper@42394
   530
neuper@52101
   531
trace_rewrite := false;
neuper@42394
   532
val SOME (t', []) = rewrite_set_ thy true RatEq_eliminate t; (*= [] must be = "x ~= 0"*)
neuper@42394
   533
term2str t' = "1 = 5 * x";
neuper@42394
   534
(*
neuper@42394
   535
 :
neuper@42394
   536
 ####  rls: rateq_erls on: x ~= 0 
neuper@42394
   537
 :
neuper@42394
   538
 #####  try calc: HOL.eq'    <<<------------------------------- here the error comes from
neuper@42394
   539
 =====  calc. to: ~ False 
neuper@42394
   540
 #####  try calc: HOL.eq' 
neuper@42394
   541
 #####  try thm: not_true 
neuper@42394
   542
 #####  try thm: not_false 
neuper@42394
   543
 =====  rewrites to: True 
neuper@42394
   544
 :
neuper@42394
   545
 ###  asms accepted: ["x ~= 0"]   stored: []
neuper@42394
   546
 :
neuper@42394
   547
*)
neuper@42394
   548
trace_rewrite := false;
neuper@42394
   549
(* WN120317.TODO dropped rateq: the above error is the same in 2002 *)
neuper@42394
   550
wneuper@59110
   551
"----------- fun assoc_thm' -----------------------------";
wneuper@59110
   552
"----------- fun assoc_thm' -----------------------------";
wneuper@59110
   553
"----------- fun assoc_thm' -----------------------------";
wneuper@59110
   554
val thy = @{theory ProgLang}
wneuper@59110
   555
wneuper@59110
   556
val tth = assoc_thm' thy ("sym_#mult_2_3","6 = 2 * 3");
wneuper@59110
   557
if string_of_thm' thy tth = "6 = 2 * 3" then ()
wneuper@59110
   558
else error "assoc_thm' (sym_#mult_2_3, 6 = 2 * 3) changed";
wneuper@59110
   559
wneuper@59110
   560
val tth = assoc_thm' thy ("add_0_left","");
wneuper@59110
   561
if string_of_thm' thy tth = "0 + ?a = ?a" then ()
wneuper@59110
   562
else error "assoc_thm' (add_0_left,\"\") changed";
wneuper@59110
   563
wneuper@59110
   564
val tth = assoc_thm' thy ("sym_add_0_left","");
wneuper@59110
   565
if string_of_thm' thy tth = "?t = 0 + ?t" then ()
wneuper@59110
   566
else error "assoc_thm' (sym_add_0_left,\"\") changed";
wneuper@59110
   567
wneuper@59110
   568
val tth = assoc_thm' thy ("add_commute","");
wneuper@59110
   569
if string_of_thm' thy tth = "?a + ?b = ?b + ?a" then ()
wneuper@59110
   570
else error "assoc_thm' (add_commute,\"\") changed"
wneuper@59110
   571
wneuper@59252
   572
"----------- fun rewrite_ down to Pattern.match ------------------------------------------------";
wneuper@59252
   573
"----------- fun rewrite_ down to Pattern.match ------------------------------------------------";
wneuper@59252
   574
"----------- fun rewrite_ down to Pattern.match ------------------------------------------------";
wneuper@59252
   575
"~~~~~ fun rewrite_ , args:"; val (thy, rew_ord, erls, bool, thm, term) =
wneuper@59252
   576
  (@{theory}, dummy_ord, e_rls, false, @{thm distrib_left}, @{term "x * (y + z) :: int"});
wneuper@59252
   577
"~~~~~ fun rewrite__ , args:"; val (thy, i, bdv, tless, rls, put_asm, thm, ct) =
wneuper@59252
   578
  (thy, 1, [], rew_ord, erls, bool, thm, term);
wneuper@59252
   579
"~~~~~ and rew_sub , args:"; val (thy, i, bdv, tless, rls, put_asm, lrd, r, t) =
wneuper@59252
   580
  (thy, i, bdv, tless, rls, put_asm, [], (((inst_bdv bdv) o norm o #prop o Thm.rep_thm) thm), ct)
wneuper@59252
   581
     val (lhss, rhss) = (HOLogic.dest_eq o HOLogic.dest_Trueprop o Logic.strip_imp_concl) r
wneuper@59252
   582
     val r' = Envir.subst_term (Pattern.match thy (lhss, t) (Vartab.empty, Vartab.empty)) r
wneuper@59252
   583
     val p' = map HOLogic.dest_Trueprop ((fst o Logic.strip_prems) (Logic.count_prems r', [], r'))
wneuper@59252
   584
     val t' = (snd o HOLogic.dest_eq o HOLogic.dest_Trueprop o Logic.strip_imp_concl) r'
wneuper@59252
   585
;
wneuper@59252
   586
if term2str lhss = "?a * (?b + ?c)" andalso term2str t = "x * (y + z)" then ()
wneuper@59252
   587
else error "ARGS FOR Pattern.match CHANGED";
wneuper@59252
   588
val match = Pattern.match thy (lhss, t) (Vartab.empty, Vartab.empty);
wneuper@59252
   589
if (Envir.subst_term match r |> term2str) = "x * (y + z) = x * y + x * z" then ()
wneuper@59252
   590
  else error "Pattern.match CHANGED";
wneuper@59381
   591
wneuper@59381
   592
"----------- fun mk_thm ------------------------------------------------------------------------";
wneuper@59381
   593
"----------- fun mk_thm ------------------------------------------------------------------------";
wneuper@59381
   594
"----------- fun mk_thm ------------------------------------------------------------------------";
wneuper@59381
   595
(*
wneuper@59381
   596
  val str = "?r ^^^ 2 = ?r * ?r";
wneuper@59381
   597
  val thm = realpow_twoI;
wneuper@59381
   598
wneuper@59381
   599
  val t1 = (#prop o rep_thm) (num_str thm);
wneuper@59381
   600
  val t2 = Trueprop $ ((Thm.term_of o the o (parse thy)) str);
wneuper@59381
   601
  t1 = t2;
wneuper@59381
   602
val it = true : bool      ... !!!
wneuper@59381
   603
  val th1 = (num_str thm);
wneuper@59381
   604
  val th2 = ((*num_str*) (mk_thm thy str)) handle e => print_exn e;
wneuper@59381
   605
  th1 = th2;
wneuper@59381
   606
ML> val it = false : bool ... HIDDEN DIFFERENCES IRRELEVANT FOR ISAC ?!
wneuper@59381
   607
wneuper@59381
   608
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
wneuper@59381
   609
  val str = "k ~= 0 ==> m * k / (n * k) = m / n";
wneuper@59381
   610
  val thm = real_mult_div_cancel2;
wneuper@59381
   611
wneuper@59381
   612
  val t1 = (#prop o rep_thm) (num_str thm);
wneuper@59381
   613
  val t2 = ((Thm.term_of o the o (parse thy)) str);
wneuper@59381
   614
  t1 = t2;
wneuper@59381
   615
val it = false : bool     ... Var .. Free
wneuper@59381
   616
  val th1 = (num_str thm);
wneuper@59381
   617
  val th2 = ((*num_str*) (mk_thm thy str)) handle e => print_exn e;
wneuper@59381
   618
  th1 = th2;
wneuper@59381
   619
ML> val it = false : bool ... PLUS HIDDEN DIFFERENCES IRRELEVANT FOR ISAC ?!
wneuper@59381
   620
*)