test/Tools/isac/Knowledge/polyminus.sml
author Walther Neuper <neuper@ist.tugraz.at>
Wed, 29 Dec 2010 20:07:52 +0100
branchdecompose-isar
changeset 38080 53ee777684ca
parent 38050 4c52ad406c20
child 38083 a1d13f3de312
permissions -rw-r--r--
intermed. in course/T3_MathEngine.thy

two problems encountered
(1) isabelle jedit does NOT accept syntax in Test_Isac.thy,
which is accepted by isabelle emacs
Thus unnecessarily edited test/../rational.sml
(2) There is an error at Model_Problem in me poly_minus.
This will be fixed before continuing T3_MathEngine.thy

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