test/Tools/isac/Knowledge/polyminus.sml
author Walther Neuper <neuper@ist.tugraz.at>
Fri, 01 Oct 2010 16:29:33 +0200
branchisac-update-Isa09-2
changeset 38038 d164e328db74
parent 38037 a51a70334191
child 38050 4c52ad406c20
permissions -rw-r--r--
tuned: tests go through
neuper@37906
     1
(* tests on PolyMinus
neuper@37906
     2
   author: Walther Neuper
neuper@37906
     3
   WN071207,
neuper@37906
     4
   (c) due to copyright terms
neuper@37906
     5
neuper@37906
     6
use"../smltest/IsacKnowledge/polyminus.sml";
neuper@37906
     7
use"polyminus.sml";
neuper@37906
     8
*)
neuper@38037
     9
(*========== inhibit exn ?======================================================
neuper@38037
    10
"--------------------------------------------------------";
neuper@38037
    11
"--------------------------------------------------------";
neuper@38037
    12
"table of contents --------------------------------------";
neuper@38037
    13
"--------------------------------------------------------";
neuper@38037
    14
"----------- fun eval_ist_monom -------------------------";
neuper@38037
    15
"----------- watch order_add_mult  ----------------------";
neuper@38037
    16
"----------- build predicate for +- ordering ------------";
neuper@38037
    17
"----------- build fasse_zusammen -----------------------";
neuper@38037
    18
"----------- build verschoenere -------------------------";
neuper@38037
    19
"----------- met simplification for_polynomials with_minu";
neuper@38037
    20
"----------- pbl polynom vereinfachen p.33 --------------";
neuper@38037
    21
"----------- met probe fuer_polynom ---------------------";
neuper@38037
    22
"----------- pbl polynom probe --------------------------";
neuper@38037
    23
"----------- pbl klammer polynom vereinfachen p.34 ------";
neuper@38037
    24
"----------- try fun applyTactics -----------------------";
neuper@38037
    25
"----------- pbl binom polynom vereinfachen p.39 --------";
neuper@38037
    26
"----------- pbl binom polynom vereinfachen: cube -------";
neuper@38037
    27
"----------- refine Vereinfache -------------------------";
neuper@38037
    28
"----------- *** prep_pbt: syntax error in '#Where' of [v";
neuper@38037
    29
"--------------------------------------------------------";
neuper@38037
    30
"--------------------------------------------------------";
neuper@38037
    31
"--------------------------------------------------------";
neuper@38037
    32
neuper@37906
    33
val thy = PolyMinus.thy;
neuper@37906
    34
neuper@37906
    35
"----------- fun eval_ist_monom ----------------------------------";
neuper@37906
    36
"----------- fun eval_ist_monom ----------------------------------";
neuper@37906
    37
"----------- fun eval_ist_monom ----------------------------------";
neuper@37906
    38
ist_monom (str2term "12");
neuper@37906
    39
case eval_ist_monom 0 0 (str2term "12 ist_monom") 0 of
neuper@37926
    40
    SOME ("12 ist_monom = True", _) => ()
neuper@38031
    41
  | _ => error "polyminus.sml: 12 ist_monom = True";
neuper@37906
    42
neuper@37906
    43
case eval_ist_monom 0 0 (str2term "a ist_monom") 0 of
neuper@37926
    44
    SOME ("a ist_monom = True", _) => ()
neuper@38031
    45
  | _ => error "polyminus.sml: a ist_monom = True";
neuper@37906
    46
neuper@37906
    47
case eval_ist_monom 0 0 (str2term "(3*a) ist_monom") 0 of
neuper@37926
    48
    SOME ("3 * a ist_monom = True", _) => ()
neuper@38031
    49
  | _ => error "polyminus.sml: 3 * a ist_monom = True";
neuper@37906
    50
neuper@37906
    51
case eval_ist_monom 0 0 (str2term "(a^^^2) ist_monom") 0 of 
neuper@37926
    52
   SOME ("a ^^^ 2 ist_monom = True", _) => ()
neuper@38031
    53
  | _ => error "polyminus.sml: a^^^2 ist_monom = True";
neuper@37906
    54
neuper@37906
    55
case eval_ist_monom 0 0 (str2term "(3*a^^^2) ist_monom") 0 of
neuper@37926
    56
    SOME ("3 * a ^^^ 2 ist_monom = True", _) => ()
neuper@38031
    57
  | _ => error "polyminus.sml: 3*a^^^2 ist_monom = True";
neuper@37906
    58
neuper@37906
    59
case eval_ist_monom 0 0 (str2term "(a*b) ist_monom") 0 of
neuper@37926
    60
    SOME ("a * b ist_monom = True", _) => ()
neuper@38031
    61
  | _ => error "polyminus.sml: a*b ist_monom = True";
neuper@37906
    62
neuper@37906
    63
case eval_ist_monom 0 0 (str2term "(3*a*b) ist_monom") 0 of
neuper@37926
    64
    SOME ("3 * a * b ist_monom = True", _) => ()
neuper@38031
    65
  | _ => error "polyminus.sml: 3*a*b ist_monom = True";
neuper@37906
    66
neuper@37906
    67
neuper@37906
    68
"----------- watch order_add_mult  -------------------------------";
neuper@37906
    69
"----------- watch order_add_mult  -------------------------------";
neuper@37906
    70
"----------- watch order_add_mult  -------------------------------";
neuper@37906
    71
"----- with these simple variables it works...";
neuper@37906
    72
trace_rewrite:=true;
neuper@37906
    73
trace_rewrite:=false;
neuper@37906
    74
val t = str2term "((a + d) + c) + b";
neuper@37926
    75
val SOME (t,_) = rewrite_set_ thy false order_add_mult t; term2str t;
neuper@37906
    76
if term2str t = "a + (b + (c + d))" then ()
neuper@38031
    77
else error "polyminus.sml 1 watch order_add_mult";
neuper@37906
    78
trace_rewrite:=false;
neuper@37906
    79
neuper@37906
    80
"----- the same stepwise...";
neuper@37906
    81
val od = ord_make_polynomial true Poly.thy;
neuper@37906
    82
val t = str2term "((a + d) + c) + b";
neuper@37906
    83
"((a + d) + c) + b"; 
neuper@37971
    84
val SOME (t,_) = rewrite_ thy od e_rls true add_commute t; term2str t;
neuper@37906
    85
"b + ((a + d) + c)";
neuper@37971
    86
val SOME (t,_) = rewrite_ thy od e_rls true add_commute t; term2str t;
neuper@37906
    87
"b + (c + (a + d))";
neuper@37971
    88
val SOME (t,_) = rewrite_ thy od e_rls true add_left_commute t;term2str t;
neuper@37906
    89
"b + (a + (c + d))";
neuper@37971
    90
val SOME (t,_) = rewrite_ thy od e_rls true add_left_commute t;term2str t;
neuper@37906
    91
"a + (b + (c + d))";
neuper@37906
    92
if term2str t = "a + (b + (c + d))" then ()
neuper@38031
    93
else error "polyminus.sml 2 watch order_add_mult";
neuper@37906
    94
neuper@37906
    95
"----- if parentheses are right, left_commute is (almost) sufficient...";
neuper@37906
    96
val t = str2term "a + (d + (c + b))";
neuper@37906
    97
"a + (d + (c + b))";
neuper@37971
    98
val SOME (t,_) = rewrite_ thy od e_rls true add_left_commute t;term2str t;
neuper@37906
    99
"a + (c + (d + b))";
neuper@37971
   100
val SOME (t,_) = rewrite_ thy od e_rls true add_commute t;term2str t;
neuper@37906
   101
"a + (c + (b + d))";
neuper@37971
   102
val SOME (t,_) = rewrite_ thy od e_rls true add_left_commute t;term2str t;
neuper@37906
   103
"a + (b + (c + d))";
neuper@37906
   104
neuper@37906
   105
"----- but we do not want the parentheses at right; thus: cond.rew.";
neuper@37906
   106
"WN0712707 complicated monomials do not yet work ...";
neuper@37906
   107
val t = str2term "((5*a + 4*d) + 3*c) + 2*b";
neuper@37926
   108
val SOME (t,_) = rewrite_set_ thy false order_add_mult t; term2str t;
neuper@37906
   109
if term2str t = "2 * b + (3 * c + (4 * d + 5 * a))" then ()
neuper@38031
   110
else error "polyminus.sml: order_add_mult changed";
neuper@37906
   111
neuper@37906
   112
"----- here we see rew_sub going into subterm with ord.rew....";
neuper@37906
   113
val od = ord_make_polynomial false Poly.thy;
neuper@37906
   114
val t = str2term "b + a + c + d";
neuper@37971
   115
val SOME (t,_) = rewrite_ thy od e_rls false add_commute t; term2str t;
neuper@37971
   116
val SOME (t,_) = rewrite_ thy od e_rls false add_commute t; term2str t;
neuper@37906
   117
(*@@@ rew_sub gosub: t = d + (b + a + c)
neuper@37906
   118
  @@@ rew_sub begin: t = b + a + c*)
neuper@37906
   119
neuper@37906
   120
neuper@37906
   121
"----------- build predicate for +- ordering ---------------------";
neuper@37906
   122
"----------- build predicate for +- ordering ---------------------";
neuper@37906
   123
"----------- build predicate for +- ordering ---------------------";
neuper@37906
   124
"a" < "b";
neuper@37906
   125
"ba" < "ab";
neuper@37906
   126
"123" < "a"; (*unused due to ---vvv*)
neuper@37906
   127
"12" < "3"; (*true !!!*)
neuper@37906
   128
neuper@37906
   129
" a kleiner b ==> (b + a) = (a + b)";
neuper@37906
   130
str2term "aaa";
neuper@37906
   131
str2term "222 * aaa";
neuper@37906
   132
(*
neuper@37906
   133
case eval_kleiner 0 0 (str2term "123 kleiner 32") 0 of
neuper@37926
   134
    SOME ("12 kleiner 9 = False", _) => ()
neuper@38031
   135
  | _ => error "polyminus.sml: 12 kleiner 9 = False";
neuper@37906
   136
*)
neuper@37906
   137
case eval_kleiner 0 0 (str2term "a kleiner b") 0 of
neuper@37926
   138
    SOME ("a kleiner b = True", _) => ()
neuper@38031
   139
  | _ => error "polyminus.sml: a kleiner b = True";
neuper@37906
   140
neuper@37906
   141
case eval_kleiner 0 0 (str2term "(10*g) kleiner f") 0 of
neuper@37926
   142
    SOME ("10 * g kleiner f = False", _) => ()
neuper@38031
   143
  | _ => error "polyminus.sml: 10 * g kleiner f = False";
neuper@37906
   144
neuper@37906
   145
case eval_kleiner 0 0 (str2term "(a^^^2) kleiner b") 0 of
neuper@37926
   146
    SOME ("a ^^^ 2 kleiner b = True", _) => ()
neuper@38031
   147
  | _ => error "polyminus.sml: a ^^^ 2 kleiner b = True";
neuper@37906
   148
neuper@37906
   149
case eval_kleiner 0 0 (str2term "(3*a^^^2) kleiner b") 0 of
neuper@37926
   150
    SOME ("3 * a ^^^ 2 kleiner b = True", _) => ()
neuper@38031
   151
  | _ => error "polyminus.sml: 3 * a ^^^ 2 kleiner b = True";
neuper@37906
   152
neuper@37906
   153
case eval_kleiner 0 0 (str2term "(a*b) kleiner c") 0 of
neuper@37926
   154
    SOME ("a * b kleiner c = True", _) => ()
neuper@38031
   155
  | _ => error "polyminus.sml: a * b kleiner b = True";
neuper@37906
   156
neuper@37906
   157
case eval_kleiner 0 0 (str2term "(3*a*b) kleiner c") 0 of
neuper@37926
   158
    SOME ("3 * a * b kleiner c = True", _) => ()
neuper@38031
   159
  | _ => error "polyminus.sml: 3 * a * b kleiner b = True";
neuper@37906
   160
neuper@37906
   161
neuper@37906
   162
neuper@37906
   163
"----- compare tausche_plus with real_num_collect";
neuper@37906
   164
val od = dummy_ord;
neuper@37906
   165
neuper@37906
   166
val erls = erls_ordne_alphabetisch;
neuper@37906
   167
val t = str2term "b + a";
neuper@37926
   168
val SOME (t,_) = rewrite_ thy od erls false tausche_plus t; term2str t;
neuper@37906
   169
if term2str t = "a + b" then ()
neuper@38031
   170
else error "polyminus.sml: ordne_alphabetisch1 b + a";
neuper@37906
   171
neuper@37906
   172
val erls = Atools_erls;
neuper@37906
   173
val t = str2term "2*a + 3*a";
neuper@37926
   174
val SOME (t,_) = rewrite_ thy od erls false real_num_collect t; term2str t;
neuper@37906
   175
neuper@37906
   176
"----- test rewrite_, rewrite_set_";
neuper@37906
   177
trace_rewrite:=true;
neuper@37906
   178
val erls = erls_ordne_alphabetisch;
neuper@37906
   179
val t = str2term "b + a";
neuper@37926
   180
val SOME (t,_) = rewrite_set_ thy false ordne_alphabetisch t; term2str t;
neuper@37906
   181
if term2str t = "a + b" then ()
neuper@38031
   182
else error "polyminus.sml: ordne_alphabetisch a + b";
neuper@37906
   183
neuper@37906
   184
val t = str2term "2*b + a";
neuper@37926
   185
val SOME (t,_) = rewrite_set_ thy false ordne_alphabetisch t; term2str t;
neuper@37906
   186
if term2str t = "a + 2 * b" then ()
neuper@38031
   187
else error "polyminus.sml: ordne_alphabetisch a + 2 * b";
neuper@37906
   188
neuper@37906
   189
val t = str2term "a + c + b";
neuper@37926
   190
val SOME (t,_) = rewrite_set_ thy false ordne_alphabetisch t; term2str t;
neuper@37906
   191
if term2str t = "a + b + c" then ()
neuper@38031
   192
else error "polyminus.sml: ordne_alphabetisch a + b + c";
neuper@37906
   193
neuper@37906
   194
"----- rewrite goes into subterms";
neuper@37906
   195
val t = str2term "a + c + b + d";
neuper@37926
   196
val SOME (t,_) = rewrite_ thy od erls false tausche_plus_plus t; term2str t;
neuper@37906
   197
if term2str t = "a + b + c + d" then ()
neuper@38031
   198
else error "polyminus.sml: ordne_alphabetisch1 a + b + c + d";
neuper@37906
   199
neuper@37906
   200
val t = str2term "a + c + d + b";
neuper@37926
   201
val SOME (t,_) = rewrite_set_ thy false ordne_alphabetisch t; term2str t;
neuper@37906
   202
if term2str t = "a + b + c + d" then ()
neuper@38031
   203
else error "polyminus.sml: ordne_alphabetisch2 a + b + c + d";
neuper@37906
   204
neuper@37906
   205
"----- here we see rew_sub going into subterm with cond.rew....";
neuper@37906
   206
val t = str2term "b + a + c + d";
neuper@37926
   207
val SOME (t,_) = rewrite_ thy od erls false tausche_plus t; term2str t;
neuper@37906
   208
if term2str t = "a + b + c + d" then ()
neuper@38031
   209
else error "polyminus.sml: ordne_alphabetisch3 a + b + c + d";
neuper@37906
   210
neuper@37906
   211
"----- compile rls for the most complicated terms";
neuper@37906
   212
val t = str2term "5*e + 6*f - 8*g - 9 - 7*e - 4*f + 10*g + 12";
neuper@37906
   213
"5 * e + 6 * f - 8 * g - 9 - 7 * e - 4 * f + 10 * g + 12";
neuper@37926
   214
val SOME (t,_) = rewrite_set_ thy false ordne_alphabetisch t; 
neuper@37906
   215
if term2str t = "- 9 + 12 + 5 * e - 7 * e + 6 * f - 4 * f - 8 * g + 10 * g"
neuper@38031
   216
then () else error "polyminus.sml: ordne_alphabetisch finished";
neuper@37906
   217
neuper@37906
   218
neuper@37906
   219
"----------- build fasse_zusammen --------------------------------";
neuper@37906
   220
"----------- build fasse_zusammen --------------------------------";
neuper@37906
   221
"----------- build fasse_zusammen --------------------------------";
neuper@37906
   222
val t = str2term "- 9 + 12 + 5 * e - 7 * e + 6 * f - 4 * f - 8 * g + 10 * g";
neuper@37926
   223
val SOME (t,_) = rewrite_set_ thy false fasse_zusammen t;
neuper@37906
   224
if term2str t = "3 + -2 * e + 2 * f + 2 * g" then ()
neuper@38031
   225
else error "polyminus.sml: fasse_zusammen finished";
neuper@37906
   226
neuper@37906
   227
"----------- build verschoenere ----------------------------------";
neuper@37906
   228
"----------- build verschoenere ----------------------------------";
neuper@37906
   229
"----------- build verschoenere ----------------------------------";
neuper@37906
   230
val t = str2term "3 + -2 * e + 2 * f + 2 * g";
neuper@37926
   231
val SOME (t,_) = rewrite_set_ thy false verschoenere t;
neuper@37906
   232
if term2str t = "3 - 2 * e + 2 * f + 2 * g" then ()
neuper@38031
   233
else error "polyminus.sml: verschoenere 3 + -2 * e ...";
neuper@37906
   234
neuper@37906
   235
trace_rewrite:=true;
neuper@37906
   236
trace_rewrite:=false;
neuper@37906
   237
neuper@37906
   238
"----------- met simplification for_polynomials with_minus -------";
neuper@37906
   239
"----------- met simplification for_polynomials with_minus -------";
neuper@37906
   240
"----------- met simplification for_polynomials with_minus -------";
neuper@37906
   241
val str = 
neuper@37906
   242
"Script SimplifyScript (t_::real) =                \
neuper@37906
   243
\  (((Try (Rewrite_Set ordne_alphabetisch False)) @@     \
neuper@37906
   244
\    (Try (Rewrite_Set fasse_zusammen False)) @@     \
neuper@37906
   245
\    (Try (Rewrite_Set verschoenere False))) t_)"
neuper@37906
   246
val sc = ((inst_abs thy) o term_of o the o (parse thy)) str;
neuper@37906
   247
atomty sc;
neuper@37906
   248
neuper@37906
   249
neuper@37906
   250
"----------- pbl polynom vereinfachen p.33 -----------------------";
neuper@37906
   251
"----------- pbl polynom vereinfachen p.33 -----------------------";
neuper@37906
   252
"----------- pbl polynom vereinfachen p.33 -----------------------";
neuper@37906
   253
"----------- 140 c ---";
neuper@37906
   254
states:=[];
neuper@37967
   255
CalcTree [(["TERM (5*e + 6*f - 8*g - 9 - 7*e - 4*f + 10*g + 12)",
neuper@37906
   256
	    "normalform N"],
neuper@37991
   257
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   258
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   259
moveActiveRoot 1;
neuper@37906
   260
autoCalculate 1 CompleteCalc;
neuper@37906
   261
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   262
if p = ([], Res) andalso 
neuper@37906
   263
   term2str (get_obj g_res pt (fst p)) = "3 - 2 * e + 2 * f + 2 * g"
neuper@38031
   264
then () else error "polyminus.sml: Vereinfache (3 - 2 * e + 2 * f...";
neuper@37906
   265
neuper@37906
   266
"----------- 140 d ---";
neuper@37906
   267
states:=[];
neuper@37967
   268
CalcTree [(["TERM (-r - 2*s - 3*t + 5 + 4*r + 8*s - 5*t - 2)",
neuper@37906
   269
	    "normalform N"],
neuper@37991
   270
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   271
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   272
moveActiveRoot 1;
neuper@37906
   273
autoCalculate 1 CompleteCalc;
neuper@37906
   274
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   275
if p = ([], Res) andalso 
neuper@37906
   276
   term2str (get_obj g_res pt (fst p)) = "3 + 3 * r + 6 * s - 8 * t"
neuper@38031
   277
then () else error "polyminus.sml: Vereinfache 140 d)";
neuper@37906
   278
neuper@37906
   279
neuper@37906
   280
"----------- 139 c ---";
neuper@37906
   281
states:=[];
neuper@37967
   282
CalcTree [(["TERM (3*e - 6*f - 8*e - 4*f + 5*e + 7*f)",
neuper@37906
   283
	    "normalform N"],
neuper@37991
   284
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   285
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   286
moveActiveRoot 1;
neuper@37906
   287
autoCalculate 1 CompleteCalc;
neuper@37906
   288
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   289
if p = ([], Res) andalso 
neuper@37906
   290
   term2str (get_obj g_res pt (fst p)) = "- (3 * f)"
neuper@38031
   291
then () else error "polyminus.sml: Vereinfache 139 c)";
neuper@37906
   292
neuper@37906
   293
"----------- 139 b ---";
neuper@37906
   294
states:=[];
neuper@37967
   295
CalcTree [(["TERM (8*u - 5*v - 5*u + 7*v - 6*u - 3*v)",
neuper@37906
   296
	    "normalform N"],
neuper@37991
   297
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   298
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   299
moveActiveRoot 1;
neuper@37906
   300
autoCalculate 1 CompleteCalc;
neuper@37906
   301
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   302
if p = ([], Res) andalso 
neuper@37906
   303
   term2str (get_obj g_res pt (fst p)) = "-3 * u - v"
neuper@38031
   304
then () else error "polyminus.sml: Vereinfache 139 b)";
neuper@37906
   305
neuper@37906
   306
"----------- 138 a ---";
neuper@37906
   307
states:=[];
neuper@37967
   308
CalcTree [(["TERM (2*u - 3*v - 6*u + 5*v)",
neuper@37906
   309
	    "normalform N"],
neuper@37991
   310
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   311
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   312
moveActiveRoot 1;
neuper@37906
   313
autoCalculate 1 CompleteCalc;
neuper@37906
   314
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   315
if p = ([], Res) andalso 
neuper@37906
   316
   term2str (get_obj g_res pt (fst p)) = "-4 * u + 2 * v"
neuper@38031
   317
then () else error "polyminus.sml: Vereinfache 138 a)";
neuper@37906
   318
neuper@37906
   319
neuper@37906
   320
"----------- met probe fuer_polynom ------------------------------";
neuper@37906
   321
"----------- met probe fuer_polynom ------------------------------";
neuper@37906
   322
"----------- met probe fuer_polynom ------------------------------";
neuper@37906
   323
val str = 
neuper@37982
   324
"Script ProbeScript (e_e::bool) (ws_::bool list) =\
neuper@37991
   325
\ (let e_e = Take e_e;                             \
neuper@37981
   326
\      e_e = Substitute ws_ e_e                    \
neuper@37975
   327
\ in (Repeat((Try (Repeat (Calculate TIMES))) @@  \
neuper@37975
   328
\            (Try (Repeat (Calculate PLUS ))) @@  \
neuper@37981
   329
\            (Try (Repeat (Calculate MINUS))))) e_e)"
neuper@37906
   330
val sc = ((inst_abs thy) o term_of o the o (parse thy)) str;
neuper@37906
   331
atomty sc;
neuper@37906
   332
neuper@37906
   333
neuper@37906
   334
"----------- pbl polynom probe -----------------------------------";
neuper@37906
   335
"----------- pbl polynom probe -----------------------------------";
neuper@37906
   336
"----------- pbl polynom probe -----------------------------------";
neuper@37906
   337
states:=[];
neuper@37906
   338
CalcTree [(["Pruefe (5*e + 6*f - 8*g - 9 - 7*e - 4*f + 10*g + 12 =\
neuper@37906
   339
	    \3 - 2 * e + 2 * f + 2 * g)",
neuper@37906
   340
	    "mitWert [e = 1, f = 2, g = 3]",
neuper@37906
   341
	    "Geprueft b"],
neuper@37991
   342
	   ("PolyMinus",["polynom","probe"],
neuper@37906
   343
	    ["probe","fuer_polynom"]))];
neuper@37906
   344
moveActiveRoot 1;
neuper@37906
   345
autoCalculate 1 CompleteCalc;
neuper@37906
   346
(* autoCalculate 1 CompleteCalcHead;
neuper@37906
   347
   autoCalculate 1 (Step 1);
neuper@37906
   348
   autoCalculate 1 (Step 1);
neuper@37906
   349
   val ((pt,p),_) = get_calc 1; term2str (get_obj g_res pt (fst p));
neuper@37906
   350
@@@@@WN081114 gives "??.empty", all "Pruefe" are the same,
neuper@37906
   351
although analogies work in interface.sml: FIXME.WN081114 in "Pruefe"*)
neuper@37906
   352
val ((pt,p),_) = get_calc 1;
neuper@37906
   353
if p = ([], Res) andalso term2str (get_obj g_res pt (fst p)) = "11 = 11"
neuper@38031
   354
then () else error "polyminus.sml: Probe 11 = 11";
neuper@37906
   355
show_pt pt;
neuper@37906
   356
neuper@37906
   357
neuper@37906
   358
"----------- pbl klammer polynom vereinfachen p.34 ---------------";
neuper@37906
   359
"----------- pbl klammer polynom vereinfachen p.34 ---------------";
neuper@37906
   360
"----------- pbl klammer polynom vereinfachen p.34 ---------------";
neuper@37906
   361
states:=[];
neuper@37967
   362
CalcTree [(["TERM (2*u - 5 - (3 - 4*u) + (8*u + 9))",
neuper@37906
   363
	    "normalform N"],
neuper@37991
   364
	   ("PolyMinus",["klammer","polynom","vereinfachen"],
neuper@37906
   365
	    ["simplification","for_polynomials","with_parentheses"]))];
neuper@37906
   366
moveActiveRoot 1;
neuper@37906
   367
autoCalculate 1 CompleteCalc;
neuper@37906
   368
val ((pt,p),_) = get_calc 1;
neuper@37906
   369
if p = ([], Res) andalso 
neuper@37906
   370
   term2str (get_obj g_res pt (fst p)) = "1 + 14 * u"
neuper@38031
   371
then () else error "polyminus.sml: Vereinfache (2*u - 5 - (3 - ...";
neuper@37906
   372
show_pt pt;
neuper@37906
   373
neuper@37906
   374
"----- probe p.34 -----";
neuper@37906
   375
states:=[];
neuper@37906
   376
CalcTree [(["Pruefe (2*u - 5 - (3 - 4*u) + (8*u + 9) = 1 + 14 * u)",
neuper@37906
   377
	    "mitWert [u = 2]",
neuper@37906
   378
	    "Geprueft b"],
neuper@37991
   379
	   ("PolyMinus",["polynom","probe"],
neuper@37906
   380
	    ["probe","fuer_polynom"]))];
neuper@37906
   381
moveActiveRoot 1;
neuper@37906
   382
autoCalculate 1 CompleteCalc;
neuper@37906
   383
val ((pt,p),_) = get_calc 1;
neuper@37906
   384
if p = ([], Res) andalso term2str (get_obj g_res pt (fst p)) = "29 = 29"
neuper@38031
   385
then () else error "polyminus.sml: Probe 29 = 29";
neuper@37906
   386
show_pt pt;
neuper@37906
   387
neuper@37906
   388
neuper@37906
   389
"----------- try fun applyTactics --------------------------------";
neuper@37906
   390
"----------- try fun applyTactics --------------------------------";
neuper@37906
   391
"----------- try fun applyTactics --------------------------------";
neuper@37906
   392
states:=[];
neuper@37967
   393
CalcTree [(["TERM (5*e + 6*f - 8*g - 9 - 7*e - 4*f + 10*g + 12)",
neuper@37906
   394
	    "normalform N"],
neuper@37991
   395
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   396
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   397
moveActiveRoot 1;
neuper@37906
   398
autoCalculate 1 CompleteCalcHead;
neuper@37906
   399
autoCalculate 1 (Step 1);
neuper@37906
   400
autoCalculate 1 (Step 1);
neuper@37906
   401
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   402
"----- 1 ^^^";
neuper@37906
   403
fetchApplicableTactics 1 0 p;
neuper@37906
   404
val appltacs = sel_appl_atomic_tacs pt p;
neuper@37906
   405
applyTactic 1 p (hd appltacs) (*addiere_x_plus_minus*);
neuper@37906
   406
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   407
"----- 2 ^^^";
neuper@37906
   408
trace_rewrite := true;
neuper@37906
   409
val erls = erls_ordne_alphabetisch;
neuper@37906
   410
val t = str2term "- 9 + 12 + 5 * e - 7 * e + (6 - 4) * f - 8 * g + 10 * g";
neuper@37926
   411
val SOME (t',_) = 
neuper@37906
   412
    rewrite_ Isac.thy e_rew_ord erls false tausche_minus t;
neuper@37906
   413
term2str t';     "- 9 + 12 + 5 * e - 7 * e + (- 4 + 6) * f - 8 * g + 10 * g";
neuper@37906
   414
neuper@37906
   415
val t = str2term "- 9 + 12 + 5 * e - 7 * e + (6 - 4) * f - 8 * g + 10 * g";
neuper@37926
   416
val NONE = 
neuper@37906
   417
    rewrite_ Isac.thy e_rew_ord erls false tausche_minus_plus t;
neuper@37906
   418
neuper@37906
   419
val t = str2term "- 9 + 12 + 5 * e - 7 * e + (6 - 4) * f - 8 * g + 10 * g";
neuper@37926
   420
val SOME (t',_) = 
neuper@37906
   421
    rewrite_set_ Isac.thy false ordne_alphabetisch t;
neuper@37906
   422
term2str t';     "- 9 + 12 + 5 * e - 7 * e - 8 * g + 10 * g + (- 4 + 6) * f";
neuper@37906
   423
trace_rewrite := false;
neuper@37906
   424
neuper@37906
   425
neuper@37906
   426
applyTactic 1 p (hd (sel_appl_atomic_tacs pt p)) (*tausche_minus*);
neuper@37906
   427
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   428
"----- 3 ^^^";
neuper@37906
   429
applyTactic 1 p (hd (sel_appl_atomic_tacs pt p)) (**);
neuper@37906
   430
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   431
"----- 4 ^^^";
neuper@37906
   432
applyTactic 1 p (hd (sel_appl_atomic_tacs pt p)) (**);
neuper@37906
   433
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   434
"----- 5 ^^^";
neuper@37906
   435
applyTactic 1 p (hd (sel_appl_atomic_tacs pt p)) (**);
neuper@37906
   436
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   437
"----- 6 ^^^";
neuper@37906
   438
neuper@37906
   439
(*<CALCMESSAGE> failure </CALCMESSAGE>
neuper@37906
   440
applyTactic 1 p (hd (sel_appl_atomic_tacs pt p)) (**);
neuper@37906
   441
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   442
"----- 7 ^^^";
neuper@37906
   443
*)
neuper@37906
   444
neuper@37906
   445
autoCalculate 1 CompleteCalc;
neuper@37906
   446
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   447
(*independent from failure above: met_simp_poly_minus not confluent:
neuper@37906
   448
(([9], Res), - (8 * g) + 10 * g + (3 - 2 * e + 2 * f)),
neuper@37906
   449
(([], Res), - (8 * g) + 10 * g + (3 - 2 * e + 2 * f))]
neuper@37906
   450
~~~~~~~~~~~###~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*)
neuper@37906
   451
neuper@37906
   452
states:=[];
neuper@37967
   453
CalcTree [(["TERM (- (8 * g) + 10 * g + h)",
neuper@37906
   454
	    "normalform N"],
neuper@37991
   455
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   456
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   457
moveActiveRoot 1;
neuper@37906
   458
autoCalculate 1 CompleteCalc;
neuper@37906
   459
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   460
if p = ([], Res) andalso term2str (get_obj g_res pt (fst p)) = "2 * g + h"
neuper@38031
   461
then () else error "polyminus.sml: addiere_vor_minus";
neuper@37906
   462
neuper@37906
   463
neuper@37906
   464
states:=[];
neuper@37967
   465
CalcTree [(["TERM (- (8 * g) + 10 * g + f)",
neuper@37906
   466
	    "normalform N"],
neuper@37991
   467
	   ("PolyMinus",["plus_minus","polynom","vereinfachen"],
neuper@37906
   468
	    ["simplification","for_polynomials","with_minus"]))];
neuper@37906
   469
moveActiveRoot 1;
neuper@37906
   470
autoCalculate 1 CompleteCalc;
neuper@37906
   471
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   472
if p = ([], Res) andalso term2str (get_obj g_res pt (fst p)) = "f + 2 * g"
neuper@38031
   473
then () else error "polyminus.sml: tausche_vor_plus";
neuper@37906
   474
neuper@37906
   475
neuper@37906
   476
"----------- pbl binom polynom vereinfachen p.39 -----------------";
neuper@37906
   477
"----------- pbl binom polynom vereinfachen p.39 -----------------";
neuper@37906
   478
"----------- pbl binom polynom vereinfachen p.39 -----------------";
neuper@37906
   479
val rls = klammern_ausmultiplizieren;
neuper@37906
   480
val t = str2term "(3 * a + 2) * (4 * a - 1)";
neuper@37926
   481
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   482
"3 * a * (4 * a) - 3 * a * 1 + (2 * (4 * a) - 2 * 1)";
neuper@37906
   483
val rls = discard_parentheses;
neuper@37926
   484
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   485
"3 * a * 4 * a - 3 * a * 1 + (2 * 4 * a - 2 * 1)";
neuper@37906
   486
val rls = ordne_monome;
neuper@37926
   487
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   488
"3 * 4 * a * a - 1 * 3 * a + (2 * 4 * a - 1 * 2)";
neuper@37906
   489
(*
neuper@37906
   490
val t = str2term "3 * a * 4 * a";
neuper@37906
   491
val rls = ordne_monome;
neuper@37926
   492
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   493
*)
neuper@37906
   494
val rls = klammern_aufloesen;
neuper@37926
   495
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   496
"3 * 4 * a * a - 1 * 3 * a + 2 * 4 * a - 1 * 2";
neuper@37906
   497
val rls = ordne_alphabetisch;
neuper@37906
   498
(*TODO: make is_monom more general, a*a=a^2, ...*)
neuper@37926
   499
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   500
"3 * 4 * a * a - 1 * 2 - 1 * 3 * a + 2 * 4 * a";
neuper@37906
   501
(*GOON.WN080104
neuper@37906
   502
val rls = fasse_zusammen;
neuper@37926
   503
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   504
val rls = verschoenere;
neuper@37926
   505
val SOME (t,_) = rewrite_set_ thy false rls t; term2str t;
neuper@37906
   506
*)
neuper@37906
   507
neuper@37906
   508
neuper@37906
   509
trace_rewrite := true;
neuper@37906
   510
trace_rewrite := false;
neuper@37906
   511
neuper@37906
   512
(*@@@@@@@*)
neuper@37906
   513
states:=[];
neuper@37967
   514
CalcTree [(["TERM ((3*a + 2) * (4*a - 1))",
neuper@37906
   515
	    "normalform N"],
neuper@37991
   516
	   ("PolyMinus",["binom_klammer","polynom","vereinfachen"],
neuper@37906
   517
	    ["simplification","for_polynomials","with_parentheses_mult"]))];
neuper@37906
   518
moveActiveRoot 1;
neuper@37906
   519
autoCalculate 1 CompleteCalc;
neuper@37906
   520
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   521
neuper@37906
   522
(*
neuper@37906
   523
if p = ([], Res) andalso 
neuper@37906
   524
   term2str (get_obj g_res pt (fst p)) = "1 + 14 * u"
neuper@38031
   525
then () else error "polyminus.sml: Vereinfache (2*u - 5 - (3 - ...";
neuper@37906
   526
*)
neuper@37906
   527
neuper@37906
   528
neuper@37906
   529
"----------- pbl binom polynom vereinfachen: cube ----------------";
neuper@37906
   530
"----------- pbl binom polynom vereinfachen: cube ----------------";
neuper@37906
   531
"----------- pbl binom polynom vereinfachen: cube ----------------";
neuper@37906
   532
states:=[];
neuper@37967
   533
CalcTree [(["TERM (8*(a - q) + a - 2*q + 3*(a - 2*q))",
neuper@37906
   534
	    "normalform N"],
neuper@37991
   535
	   ("PolyMinus",["binom_klammer","polynom","vereinfachen"],
neuper@37906
   536
	    ["simplification","for_polynomials","with_parentheses_mult"]))];
neuper@37906
   537
moveActiveRoot 1;
neuper@37906
   538
autoCalculate 1 CompleteCalc;
neuper@37906
   539
val ((pt,p),_) = get_calc 1; show_pt pt;
neuper@37906
   540
neuper@37906
   541
neuper@37906
   542
"----------- refine Vereinfache ----------------------------------";
neuper@37906
   543
"----------- refine Vereinfache ----------------------------------";
neuper@37906
   544
"----------- refine Vereinfache ----------------------------------";
neuper@37967
   545
val fmz = ["TERM (8*(a - q) + a - 2*q + 3*(a - 2*q))",
neuper@37906
   546
	    "normalform N"];
neuper@37906
   547
print_depth 11;
neuper@37906
   548
val matches = refine fmz ["vereinfachen"];
neuper@37906
   549
print_depth 3;
neuper@37906
   550
neuper@37906
   551
"----- go into details, if it seems not to work -----";
neuper@37906
   552
"--- does the predicate evaluate correctly ?";
neuper@37906
   553
val t = str2term 
neuper@37906
   554
	    "matchsub (?a * (?b - ?c)) (8 * (a - q) + a - 2 * q + \
neuper@37906
   555
	    \3 * (a - 2 * q))";
neuper@37906
   556
val ma = eval_matchsub "" "Tools.matchsub" t thy;
neuper@37906
   557
case ma of
neuper@37926
   558
    SOME ("matchsub (?a * (?b - ?c)) (8 * (a - q) + \
neuper@37906
   559
	  \a - 2 * q + 3 * (a - 2 * q)) = True", _) => ()
neuper@38031
   560
  | _ => error "polyminus.sml matchsub (?a * (?b - ?c)...A";
neuper@37906
   561
neuper@37906
   562
"--- does the respective prls rewrite ?";
neuper@37906
   563
val prls = append_rls "prls_pbl_vereinf_poly" e_rls 
neuper@37906
   564
	     [Calc ("Poly.is'_polyexp", eval_is_polyexp ""),
neuper@37906
   565
	      Calc ("Tools.matchsub", eval_matchsub ""),
neuper@37906
   566
	      Thm ("or_true",or_true),
neuper@37906
   567
	      (*"(?a | True) = True"*)
neuper@37906
   568
	      Thm ("or_false",or_false),
neuper@37906
   569
	      (*"(?a | False) = ?a"*)
neuper@37970
   570
	      Thm ("not_true",num_str @{thm not_true}),
neuper@37906
   571
	      (*"(~ True) = False"*)
neuper@37970
   572
	      Thm ("not_false",num_str @{thm not_false})
neuper@37906
   573
	      (*"(~ False) = True"*)];
neuper@37906
   574
trace_rewrite := true;
neuper@37926
   575
val SOME (t', _) = rewrite_set_ thy false prls t;
neuper@37906
   576
trace_rewrite := false;
neuper@37906
   577
neuper@37906
   578
"--- does the respective prls rewrite the whole predicate ?";
neuper@37906
   579
val t = str2term 
neuper@37906
   580
	    "Not (matchsub (?a * (?b + ?c)) (8 * (a - q) + a - 2 * q) | \
neuper@37906
   581
	    \     matchsub (?a * (?b - ?c)) (8 * (a - q) + a - 2 * q) | \
neuper@37906
   582
	    \     matchsub ((?b + ?c) * ?a) (8 * (a - q) + a - 2 * q) | \
neuper@37906
   583
	    \     matchsub ((?b - ?c) * ?a) (8 * (a - q) + a - 2 * q) )";
neuper@37906
   584
trace_rewrite := true;
neuper@37926
   585
val SOME (t', _) = rewrite_set_ thy false prls t;
neuper@37906
   586
trace_rewrite := false;
neuper@37906
   587
if term2str t' = "False" then ()
neuper@38031
   588
else error "polyminus.sml Not (matchsub (?a * (?b + ?c)) (8 ...";
neuper@37906
   589
neuper@38037
   590
============ inhibit exn ?====================================================*)
neuper@37906
   591
neuper@37906
   592
neuper@38037
   593
"----------- *** prep_pbt: syntax error in '#Where' of [v";
neuper@38037
   594
"----------- *** prep_pbt: syntax error in '#Where' of [v";
neuper@38037
   595
"----------- *** prep_pbt: syntax error in '#Where' of [v";
neuper@38037
   596
(*see test/../termC.sml for details*)
neuper@38037
   597
val t = parse_patt thy "t_t is_polyexp";
neuper@38037
   598
val t = parse_patt thy ("Not (matchsub (?a + (?b + ?c)) t_t | " ^
neuper@38037
   599
	                "     matchsub (?a + (?b - ?c)) t_t | " ^
neuper@38037
   600
	                "     matchsub (?a - (?b + ?c)) t_t | " ^
neuper@38037
   601
	                "     matchsub (?a + (?b - ?c)) t_t )");
neuper@38038
   602
show_types := true;
neuper@38037
   603
if term2str t = "~ (matchsub ((?a::real) + ((?b::real) + (?c::real))) (t_t::real) |\n   matchsub (?a + (?b - ?c)) t_t |\n   matchsub (?a - (?b + ?c)) t_t | matchsub (?a + (?b - ?c)) t_t)"
neuper@38038
   604
then () else error "polyminus.sml type-structure of \"?a :: real\" changed 1";
neuper@38037
   605
show_types := false;
neuper@37906
   606
neuper@38037
   607
(*========== inhibit exn =======================================================
neuper@38037
   608
============ inhibit exn =====================================================*)
neuper@37906
   609
neuper@38037
   610
(*-.-.-.-.-.-isolate response.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
neuper@38037
   611
-.-.-.-.-.-.-isolate response.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.*)