src/Tools/isac/Knowledge/PolyMinus.ML
author Walther Neuper <neuper@ist.tugraz.at>
Wed, 25 Aug 2010 16:20:07 +0200
branchisac-update-Isa09-2
changeset 37947 22235e4dbe5f
parent 37926 src/Tools/isac/IsacKnowledge/PolyMinus.ML@e6fc98fbcb85
permissions -rw-r--r--
renamed isac's directories and Build_Isac.thy

Scripts --> ProgLang
ME --> Interpret
IsacKnowledge --> Knowledge
neuper@37906
     1
(* questionable attempts to perserve binary minus as wanted by teachers
neuper@37906
     2
   WN071207
neuper@37906
     3
   (c) due to copyright terms
neuper@37906
     4
remove_thy"PolyMinus";
neuper@37947
     5
use_thy"Knowledge/PolyMinus";
neuper@37906
     6
neuper@37947
     7
use_thy"Knowledge/Isac";
neuper@37947
     8
use"Knowledge/PolyMinus.ML";
neuper@37906
     9
*)
neuper@37906
    10
neuper@37906
    11
(** interface isabelle -- isac **)
neuper@37906
    12
theory' := overwritel (!theory', [("PolyMinus.thy",PolyMinus.thy)]);
neuper@37906
    13
neuper@37906
    14
(** eval functions **)
neuper@37906
    15
neuper@37906
    16
(*. get the identifier from specific monomials; see fun ist_monom .*)
neuper@37906
    17
(*HACK.WN080107*)
neuper@37906
    18
fun increase str = 
neuper@37906
    19
    let val s::ss = explode str
neuper@37906
    20
    in implode ((chr (ord s + 1))::ss) end;
neuper@37906
    21
fun identifier (Free (id,_)) = id                            (* 2,   a   *)
neuper@37906
    22
  | identifier (Const ("op *", _) $ Free (num, _) $ Free (id, _)) = 
neuper@37906
    23
    id                                                       (* 2*a, a*b *)
neuper@37906
    24
  | identifier (Const ("op *", _) $                          (* 3*a*b    *)
neuper@37906
    25
		     (Const ("op *", _) $
neuper@37906
    26
			    Free (num, _) $ Free _) $ Free (id, _)) = 
neuper@37906
    27
    if is_numeral num then id
neuper@37906
    28
    else "|||||||||||||"
neuper@37906
    29
  | identifier (Const ("Atools.pow", _) $ Free (base, _) $ Free (exp, _)) =
neuper@37906
    30
    if is_numeral base then "|||||||||||||"                  (* a^2      *)
neuper@37906
    31
    else (*increase*) base
neuper@37906
    32
  | identifier (Const ("op *", _) $ Free (num, _) $          (* 3*a^2    *)
neuper@37906
    33
		     (Const ("Atools.pow", _) $
neuper@37906
    34
			    Free (base, _) $ Free (exp, _))) = 
neuper@37906
    35
    if is_numeral num andalso not (is_numeral base) then (*increase*) base
neuper@37906
    36
    else "|||||||||||||"
neuper@37906
    37
  | identifier _ = "|||||||||||||"(*the "largest" string*);
neuper@37906
    38
neuper@37906
    39
(*("kleiner", ("PolyMinus.kleiner", eval_kleiner ""))*)
neuper@37906
    40
(* order "by alphabet" w.r.t. var: num < (var | num*var) > (var*var | ..) *)
neuper@37906
    41
fun eval_kleiner _ _ (p as (Const ("PolyMinus.kleiner",_) $ a $ b)) _  =
neuper@37906
    42
     if is_num b then
neuper@37906
    43
	 if is_num a then (*123 kleiner 32 = True !!!*)
neuper@37906
    44
	     if int_of_Free a < int_of_Free b then 
neuper@37926
    45
		 SOME ((term2str p) ^ " = True",
neuper@37906
    46
		       Trueprop $ (mk_equality (p, HOLogic.true_const)))
neuper@37926
    47
	     else SOME ((term2str p) ^ " = False",
neuper@37906
    48
			Trueprop $ (mk_equality (p, HOLogic.false_const)))
neuper@37906
    49
	 else (* -1 * -2 kleiner 0 *)
neuper@37926
    50
	     SOME ((term2str p) ^ " = False",
neuper@37906
    51
		   Trueprop $ (mk_equality (p, HOLogic.false_const)))
neuper@37906
    52
    else
neuper@37906
    53
	if identifier a < identifier b then 
neuper@37926
    54
	     SOME ((term2str p) ^ " = True",
neuper@37906
    55
		  Trueprop $ (mk_equality (p, HOLogic.true_const)))
neuper@37926
    56
	else SOME ((term2str p) ^ " = False",
neuper@37906
    57
		   Trueprop $ (mk_equality (p, HOLogic.false_const)))
neuper@37926
    58
  | eval_kleiner _ _ _ _ =  NONE;
neuper@37906
    59
neuper@37906
    60
fun ist_monom (Free (id,_)) = true
neuper@37906
    61
  | ist_monom (Const ("op *", _) $ Free (num, _) $ Free (id, _)) = 
neuper@37906
    62
    if is_numeral num then true else false
neuper@37906
    63
  | ist_monom _ = false;
neuper@37906
    64
(*. this function only accepts the most simple monoms       vvvvvvvvvv .*)
neuper@37906
    65
fun ist_monom (Free (id,_)) = true                          (* 2,   a   *)
neuper@37906
    66
  | ist_monom (Const ("op *", _) $ Free _ $ Free (id, _)) = (* 2*a, a*b *)
neuper@37906
    67
    if is_numeral id then false else true
neuper@37906
    68
  | ist_monom (Const ("op *", _) $                          (* 3*a*b    *)
neuper@37906
    69
		     (Const ("op *", _) $
neuper@37906
    70
			    Free (num, _) $ Free _) $ Free (id, _)) =
neuper@37906
    71
    if is_numeral num andalso not (is_numeral id) then true else false
neuper@37906
    72
  | ist_monom (Const ("Atools.pow", _) $ Free (base, _) $ Free (exp, _)) = 
neuper@37906
    73
    true                                                    (* a^2      *)
neuper@37906
    74
  | ist_monom (Const ("op *", _) $ Free (num, _) $          (* 3*a^2    *)
neuper@37906
    75
		     (Const ("Atools.pow", _) $
neuper@37906
    76
			    Free (base, _) $ Free (exp, _))) = 
neuper@37906
    77
    if is_numeral num then true else false
neuper@37906
    78
  | ist_monom _ = false;
neuper@37906
    79
neuper@37906
    80
(* is this a univariate monomial ? *)
neuper@37906
    81
(*("ist_monom", ("PolyMinus.ist'_monom", eval_ist_monom ""))*)
neuper@37906
    82
fun eval_ist_monom _ _ (p as (Const ("PolyMinus.ist'_monom",_) $ a)) _  =
neuper@37906
    83
    if ist_monom a  then 
neuper@37926
    84
	SOME ((term2str p) ^ " = True",
neuper@37906
    85
	      Trueprop $ (mk_equality (p, HOLogic.true_const)))
neuper@37926
    86
    else SOME ((term2str p) ^ " = False",
neuper@37906
    87
	       Trueprop $ (mk_equality (p, HOLogic.false_const)))
neuper@37926
    88
  | eval_ist_monom _ _ _ _ =  NONE;
neuper@37906
    89
neuper@37906
    90
neuper@37906
    91
(** rewrite order **)
neuper@37906
    92
neuper@37906
    93
(** rulesets **)
neuper@37906
    94
neuper@37906
    95
val erls_ordne_alphabetisch =
neuper@37906
    96
    append_rls "erls_ordne_alphabetisch" e_rls
neuper@37906
    97
	       [Calc ("PolyMinus.kleiner", eval_kleiner ""),
neuper@37906
    98
		Calc ("PolyMinus.ist'_monom", eval_ist_monom "")
neuper@37906
    99
		];
neuper@37906
   100
neuper@37906
   101
val ordne_alphabetisch = 
neuper@37906
   102
  Rls{id = "ordne_alphabetisch", preconds = [], 
neuper@37906
   103
      rew_ord = ("dummy_ord", dummy_ord), srls = Erls, calc = [],
neuper@37906
   104
      erls = erls_ordne_alphabetisch, 
neuper@37906
   105
      rules = [Thm ("tausche_plus",num_str tausche_plus),
neuper@37906
   106
	       (*"b kleiner a ==> (b + a) = (a + b)"*)
neuper@37906
   107
	       Thm ("tausche_minus",num_str tausche_minus),
neuper@37906
   108
	       (*"b kleiner a ==> (b - a) = (-a + b)"*)
neuper@37906
   109
	       Thm ("tausche_vor_plus",num_str tausche_vor_plus),
neuper@37906
   110
	       (*"[| b ist_monom; a kleiner b  |] ==> (- b + a) = (a - b)"*)
neuper@37906
   111
	       Thm ("tausche_vor_minus",num_str tausche_vor_minus),
neuper@37906
   112
	       (*"[| b ist_monom; a kleiner b  |] ==> (- b - a) = (-a - b)"*)
neuper@37906
   113
	       Thm ("tausche_plus_plus",num_str tausche_plus_plus),
neuper@37906
   114
	       (*"c kleiner b ==> (a + c + b) = (a + b + c)"*)
neuper@37906
   115
	       Thm ("tausche_plus_minus",num_str tausche_plus_minus),
neuper@37906
   116
	       (*"c kleiner b ==> (a + c - b) = (a - b + c)"*)
neuper@37906
   117
	       Thm ("tausche_minus_plus",num_str tausche_minus_plus),
neuper@37906
   118
	       (*"c kleiner b ==> (a - c + b) = (a + b - c)"*)
neuper@37906
   119
	       Thm ("tausche_minus_minus",num_str tausche_minus_minus)
neuper@37906
   120
	       (*"c kleiner b ==> (a - c - b) = (a - b - c)"*)
neuper@37906
   121
	       ], scr = EmptyScr}:rls;
neuper@37906
   122
neuper@37906
   123
val fasse_zusammen = 
neuper@37906
   124
    Rls{id = "fasse_zusammen", preconds = [], 
neuper@37906
   125
	rew_ord = ("dummy_ord", dummy_ord),
neuper@37906
   126
	erls = append_rls "erls_fasse_zusammen" e_rls 
neuper@37906
   127
			  [Calc ("Atools.is'_const",eval_const "#is_const_")], 
neuper@37906
   128
	srls = Erls, calc = [],
neuper@37906
   129
	rules = 
neuper@37906
   130
	[Thm ("real_num_collect",num_str real_num_collect), 
neuper@37906
   131
	 (*"[| l is_const; m is_const |]==>l * n + m * n = (l + m) * n"*)
neuper@37906
   132
	 Thm ("real_num_collect_assoc_r",num_str real_num_collect_assoc_r),
neuper@37906
   133
	 (*"[| l is_const; m..|] ==>  (k + m * n) + l * n = k + (l + m)*n"*)
neuper@37906
   134
	 Thm ("real_one_collect",num_str real_one_collect),	
neuper@37906
   135
	 (*"m is_const ==> n + m * n = (1 + m) * n"*)
neuper@37906
   136
	 Thm ("real_one_collect_assoc_r",num_str real_one_collect_assoc_r), 
neuper@37906
   137
	 (*"m is_const ==> (k + n) + m * n = k + (m + 1) * n"*)
neuper@37906
   138
neuper@37906
   139
neuper@37906
   140
	 Thm ("subtrahiere",num_str subtrahiere),
neuper@37906
   141
	 (*"[| l is_const; m is_const |] ==> m * v - l * v = (m - l) * v"*)
neuper@37906
   142
	 Thm ("subtrahiere_von_1",num_str subtrahiere_von_1),
neuper@37906
   143
	 (*"[| l is_const |] ==> v - l * v = (1 - l) * v"*)
neuper@37906
   144
	 Thm ("subtrahiere_1",num_str subtrahiere_1),
neuper@37906
   145
	 (*"[| l is_const; m is_const |] ==> m * v - v = (m - 1) * v"*)
neuper@37906
   146
neuper@37906
   147
	 Thm ("subtrahiere_x_plus_minus",num_str subtrahiere_x_plus_minus), 
neuper@37906
   148
	 (*"[| l is_const; m..|] ==> (k + m * n) - l * n = k + ( m - l) * n"*)
neuper@37906
   149
	 Thm ("subtrahiere_x_plus1_minus",num_str subtrahiere_x_plus1_minus),
neuper@37906
   150
	 (*"[| l is_const |] ==> (x + v) - l * v = x + (1 - l) * v"*)
neuper@37906
   151
	 Thm ("subtrahiere_x_plus_minus1",num_str subtrahiere_x_plus_minus1),
neuper@37906
   152
	 (*"[| m is_const |] ==> (x + m * v) - v = x + (m - 1) * v"*)
neuper@37906
   153
neuper@37906
   154
	 Thm ("subtrahiere_x_minus_plus",num_str subtrahiere_x_minus_plus), 
neuper@37906
   155
	 (*"[| l is_const; m..|] ==> (k - m * n) + l * n = k + (-m + l) * n"*)
neuper@37906
   156
	 Thm ("subtrahiere_x_minus1_plus",num_str subtrahiere_x_minus1_plus),
neuper@37906
   157
	 (*"[| l is_const |] ==> (x - v) + l * v = x + (-1 + l) * v"*)
neuper@37906
   158
	 Thm ("subtrahiere_x_minus_plus1",num_str subtrahiere_x_minus_plus1),
neuper@37906
   159
	 (*"[| m is_const |] ==> (x - m * v) + v = x + (-m + 1) * v"*)
neuper@37906
   160
neuper@37906
   161
	 Thm ("subtrahiere_x_minus_minus",num_str subtrahiere_x_minus_minus), 
neuper@37906
   162
	 (*"[| l is_const; m..|] ==> (k - m * n) - l * n = k + (-m - l) * n"*)
neuper@37906
   163
	 Thm ("subtrahiere_x_minus1_minus",num_str subtrahiere_x_minus1_minus),
neuper@37906
   164
	 (*"[| l is_const |] ==> (x - v) - l * v = x + (-1 - l) * v"*)
neuper@37906
   165
	 Thm ("subtrahiere_x_minus_minus1",num_str subtrahiere_x_minus_minus1),
neuper@37906
   166
	 (*"[| m is_const |] ==> (x - m * v) - v = x + (-m - 1) * v"*)
neuper@37906
   167
	 
neuper@37906
   168
	 Calc ("op +", eval_binop "#add_"),
neuper@37906
   169
	 Calc ("op -", eval_binop "#subtr_"),
neuper@37906
   170
	 
neuper@37906
   171
	 (*MG: Reihenfolge der folgenden 2 Thm muss so bleiben, wegen
neuper@37906
   172
           (a+a)+a --> a + 2*a --> 3*a and not (a+a)+a --> 2*a + a *)
neuper@37906
   173
	 Thm ("real_mult_2_assoc_r",num_str real_mult_2_assoc_r),
neuper@37906
   174
	 (*"(k + z1) + z1 = k + 2 * z1"*)
neuper@37906
   175
	 Thm ("sym_real_mult_2",num_str (real_mult_2 RS sym)),
neuper@37906
   176
	 (*"z1 + z1 = 2 * z1"*)
neuper@37906
   177
neuper@37906
   178
	 Thm ("addiere_vor_minus",num_str addiere_vor_minus),
neuper@37906
   179
	 (*"[| l is_const; m is_const |] ==> -(l * v) +  m * v = (-l + m) *v"*)
neuper@37906
   180
	 Thm ("addiere_eins_vor_minus",num_str addiere_eins_vor_minus),
neuper@37906
   181
	 (*"[| m is_const |] ==> -  v +  m * v = (-1 + m) * v"*)
neuper@37906
   182
	 Thm ("subtrahiere_vor_minus",num_str subtrahiere_vor_minus),
neuper@37906
   183
	 (*"[| l is_const; m is_const |] ==> -(l * v) -  m * v = (-l - m) *v"*)
neuper@37906
   184
	 Thm ("subtrahiere_eins_vor_minus",num_str subtrahiere_eins_vor_minus)
neuper@37906
   185
	 (*"[| m is_const |] ==> -  v -  m * v = (-1 - m) * v"*)
neuper@37906
   186
	 
neuper@37906
   187
	 ], scr = EmptyScr}:rls;
neuper@37906
   188
    
neuper@37906
   189
val verschoenere = 
neuper@37906
   190
  Rls{id = "verschoenere", preconds = [], 
neuper@37906
   191
      rew_ord = ("dummy_ord", dummy_ord), srls = Erls, calc = [],
neuper@37906
   192
      erls = append_rls "erls_verschoenere" e_rls 
neuper@37906
   193
			[Calc ("PolyMinus.kleiner", eval_kleiner "")], 
neuper@37906
   194
      rules = [Thm ("vorzeichen_minus_weg1",num_str vorzeichen_minus_weg1),
neuper@37906
   195
	       (*"l kleiner 0 ==> a + l * b = a - -l * b"*)
neuper@37906
   196
	       Thm ("vorzeichen_minus_weg2",num_str vorzeichen_minus_weg2),
neuper@37906
   197
	       (*"l kleiner 0 ==> a - l * b = a + -l * b"*)
neuper@37906
   198
	       Thm ("vorzeichen_minus_weg3",num_str vorzeichen_minus_weg3),
neuper@37906
   199
	       (*"l kleiner 0 ==> k + a - l * b = k + a + -l * b"*)
neuper@37906
   200
	       Thm ("vorzeichen_minus_weg4",num_str vorzeichen_minus_weg4),
neuper@37906
   201
	       (*"l kleiner 0 ==> k - a - l * b = k - a + -l * b"*)
neuper@37906
   202
neuper@37906
   203
	       Calc ("op *", eval_binop "#mult_"),
neuper@37906
   204
neuper@37906
   205
	       Thm ("real_mult_0",num_str real_mult_0),    
neuper@37906
   206
	       (*"0 * z = 0"*)
neuper@37906
   207
	       Thm ("real_mult_1",num_str real_mult_1),     
neuper@37906
   208
	       (*"1 * z = z"*)
neuper@37906
   209
	       Thm ("real_add_zero_left",num_str real_add_zero_left),
neuper@37906
   210
	       (*"0 + z = z"*)
neuper@37906
   211
	       Thm ("null_minus",num_str null_minus),
neuper@37906
   212
	       (*"0 - a = -a"*)
neuper@37906
   213
	       Thm ("vor_minus_mal",num_str vor_minus_mal)
neuper@37906
   214
	       (*"- a * b = (-a) * b"*)
neuper@37906
   215
neuper@37906
   216
	       (*Thm ("",num_str ),*)
neuper@37906
   217
	       (**)
neuper@37906
   218
	       ], scr = EmptyScr}:rls (*end verschoenere*);
neuper@37906
   219
neuper@37906
   220
val klammern_aufloesen = 
neuper@37906
   221
  Rls{id = "klammern_aufloesen", preconds = [], 
neuper@37906
   222
      rew_ord = ("dummy_ord", dummy_ord), srls = Erls, calc = [], erls = Erls, 
neuper@37906
   223
      rules = [Thm ("sym_real_add_assoc",num_str (real_add_assoc RS sym)),
neuper@37906
   224
	       (*"a + (b + c) = (a + b) + c"*)
neuper@37906
   225
	       Thm ("klammer_plus_minus",num_str klammer_plus_minus),
neuper@37906
   226
	       (*"a + (b - c) = (a + b) - c"*)
neuper@37906
   227
	       Thm ("klammer_minus_plus",num_str klammer_minus_plus),
neuper@37906
   228
	       (*"a - (b + c) = (a - b) - c"*)
neuper@37906
   229
	       Thm ("klammer_minus_minus",num_str klammer_minus_minus)
neuper@37906
   230
	       (*"a - (b - c) = (a - b) + c"*)
neuper@37906
   231
	       ], scr = EmptyScr}:rls;
neuper@37906
   232
neuper@37906
   233
val klammern_ausmultiplizieren = 
neuper@37906
   234
  Rls{id = "klammern_ausmultiplizieren", preconds = [], 
neuper@37906
   235
      rew_ord = ("dummy_ord", dummy_ord), srls = Erls, calc = [], erls = Erls, 
neuper@37906
   236
      rules = [Thm ("real_add_mult_distrib" ,num_str real_add_mult_distrib),
neuper@37906
   237
	       (*"(z1.0 + z2.0) * w = z1.0 * w + z2.0 * w"*)
neuper@37906
   238
	       Thm ("real_add_mult_distrib2",num_str real_add_mult_distrib2),
neuper@37906
   239
	       (*"w * (z1.0 + z2.0) = w * z1.0 + w * z2.0"*)
neuper@37906
   240
	       
neuper@37906
   241
	       Thm ("klammer_mult_minus",num_str klammer_mult_minus),
neuper@37906
   242
	       (*"a * (b - c) = a * b - a * c"*)
neuper@37906
   243
	       Thm ("klammer_minus_mult",num_str klammer_minus_mult)
neuper@37906
   244
	       (*"(b - c) * a = b * a - c * a"*)
neuper@37906
   245
neuper@37906
   246
	       (*Thm ("",num_str ),
neuper@37906
   247
	       (*""*)*)
neuper@37906
   248
	       ], scr = EmptyScr}:rls;
neuper@37906
   249
neuper@37906
   250
val ordne_monome = 
neuper@37906
   251
  Rls{id = "ordne_monome", preconds = [], 
neuper@37906
   252
      rew_ord = ("dummy_ord", dummy_ord), srls = Erls, calc = [], 
neuper@37906
   253
      erls = append_rls "erls_ordne_monome" e_rls
neuper@37906
   254
	       [Calc ("PolyMinus.kleiner", eval_kleiner ""),
neuper@37906
   255
		Calc ("Atools.is'_atom", eval_is_atom "")
neuper@37906
   256
		], 
neuper@37906
   257
      rules = [Thm ("tausche_mal",num_str tausche_mal),
neuper@37906
   258
	       (*"[| b is_atom; a kleiner b  |] ==> (b * a) = (a * b)"*)
neuper@37906
   259
	       Thm ("tausche_vor_mal",num_str tausche_vor_mal),
neuper@37906
   260
	       (*"[| b is_atom; a kleiner b  |] ==> (-b * a) = (-a * b)"*)
neuper@37906
   261
	       Thm ("tausche_mal_mal",num_str tausche_mal_mal),
neuper@37906
   262
	       (*"[| c is_atom; b kleiner c  |] ==> (a * c * b) = (a * b *c)"*)
neuper@37906
   263
	       Thm ("x_quadrat",num_str x_quadrat)
neuper@37906
   264
	       (*"(x * a) * a = x * a ^^^ 2"*)
neuper@37906
   265
neuper@37906
   266
	       (*Thm ("",num_str ),
neuper@37906
   267
	       (*""*)*)
neuper@37906
   268
	       ], scr = EmptyScr}:rls;
neuper@37906
   269
neuper@37906
   270
neuper@37906
   271
val rls_p_33 = 
neuper@37906
   272
    append_rls "rls_p_33" e_rls
neuper@37906
   273
	       [Rls_ ordne_alphabetisch,
neuper@37906
   274
		Rls_ fasse_zusammen,
neuper@37906
   275
		Rls_ verschoenere
neuper@37906
   276
		];
neuper@37906
   277
val rls_p_34 = 
neuper@37906
   278
    append_rls "rls_p_34" e_rls
neuper@37906
   279
	       [Rls_ klammern_aufloesen,
neuper@37906
   280
		Rls_ ordne_alphabetisch,
neuper@37906
   281
		Rls_ fasse_zusammen,
neuper@37906
   282
		Rls_ verschoenere
neuper@37906
   283
		];
neuper@37906
   284
val rechnen = 
neuper@37906
   285
    append_rls "rechnen" e_rls
neuper@37906
   286
	       [Calc ("op *", eval_binop "#mult_"),
neuper@37906
   287
		Calc ("op +", eval_binop "#add_"),
neuper@37906
   288
		Calc ("op -", eval_binop "#subtr_")
neuper@37906
   289
		];
neuper@37906
   290
neuper@37906
   291
ruleset' := 
neuper@37906
   292
overwritelthy thy (!ruleset',
neuper@37906
   293
		   [("ordne_alphabetisch", prep_rls ordne_alphabetisch),
neuper@37906
   294
		    ("fasse_zusammen", prep_rls fasse_zusammen),
neuper@37906
   295
		    ("verschoenere", prep_rls verschoenere),
neuper@37906
   296
		    ("ordne_monome", prep_rls ordne_monome),
neuper@37906
   297
		    ("klammern_aufloesen", prep_rls klammern_aufloesen),
neuper@37906
   298
		    ("klammern_ausmultiplizieren", 
neuper@37906
   299
		     prep_rls klammern_ausmultiplizieren)
neuper@37906
   300
		    ]);
neuper@37906
   301
neuper@37906
   302
(** problems **)
neuper@37906
   303
neuper@37906
   304
store_pbt
neuper@37906
   305
 (prep_pbt PolyMinus.thy "pbl_vereinf_poly" [] e_pblID
neuper@37906
   306
 (["polynom","vereinfachen"],
neuper@37926
   307
  [], Erls, NONE, []));
neuper@37906
   308
neuper@37906
   309
store_pbt
neuper@37906
   310
 (prep_pbt PolyMinus.thy "pbl_vereinf_poly_minus" [] e_pblID
neuper@37906
   311
 (["plus_minus","polynom","vereinfachen"],
neuper@37906
   312
  [("#Given" ,["term t_"]),
neuper@37906
   313
   ("#Where" ,["t_ is_polyexp",
neuper@37906
   314
	       "Not (matchsub (?a + (?b + ?c)) t_ | \
neuper@37906
   315
	       \     matchsub (?a + (?b - ?c)) t_ | \
neuper@37906
   316
	       \     matchsub (?a - (?b + ?c)) t_ | \
neuper@37906
   317
	       \     matchsub (?a + (?b - ?c)) t_ )",
neuper@37906
   318
	       "Not (matchsub (?a * (?b + ?c)) t_ | \
neuper@37906
   319
	       \     matchsub (?a * (?b - ?c)) t_ | \
neuper@37906
   320
	       \     matchsub ((?b + ?c) * ?a) t_ | \
neuper@37906
   321
	       \     matchsub ((?b - ?c) * ?a) t_ )"]),
neuper@37906
   322
   ("#Find"  ,["normalform n_"])
neuper@37906
   323
  ],
neuper@37906
   324
  append_rls "prls_pbl_vereinf_poly" e_rls 
neuper@37906
   325
	     [Calc ("Poly.is'_polyexp", eval_is_polyexp ""),
neuper@37906
   326
	      Calc ("Tools.matchsub", eval_matchsub ""),
neuper@37906
   327
	      Thm ("or_true",or_true),
neuper@37906
   328
	      (*"(?a | True) = True"*)
neuper@37906
   329
	      Thm ("or_false",or_false),
neuper@37906
   330
	      (*"(?a | False) = ?a"*)
neuper@37906
   331
	      Thm ("not_true",num_str not_true),
neuper@37906
   332
	      (*"(~ True) = False"*)
neuper@37906
   333
	      Thm ("not_false",num_str not_false)
neuper@37906
   334
	      (*"(~ False) = True"*)], 
neuper@37926
   335
  SOME "Vereinfache t_", 
neuper@37906
   336
  [["simplification","for_polynomials","with_minus"]]));
neuper@37906
   337
neuper@37906
   338
store_pbt
neuper@37906
   339
 (prep_pbt PolyMinus.thy "pbl_vereinf_poly_klammer" [] e_pblID
neuper@37906
   340
 (["klammer","polynom","vereinfachen"],
neuper@37906
   341
  [("#Given" ,["term t_"]),
neuper@37906
   342
   ("#Where" ,["t_ is_polyexp",
neuper@37906
   343
	       "Not (matchsub (?a * (?b + ?c)) t_ | \
neuper@37906
   344
	       \     matchsub (?a * (?b - ?c)) t_ | \
neuper@37906
   345
	       \     matchsub ((?b + ?c) * ?a) t_ | \
neuper@37906
   346
	       \     matchsub ((?b - ?c) * ?a) t_ )"]),
neuper@37906
   347
   ("#Find"  ,["normalform n_"])
neuper@37906
   348
  ],
neuper@37906
   349
  append_rls "prls_pbl_vereinf_poly_klammer" e_rls [Calc ("Poly.is'_polyexp", eval_is_polyexp ""),
neuper@37906
   350
	      Calc ("Tools.matchsub", eval_matchsub ""),
neuper@37906
   351
	      Thm ("or_true",or_true),
neuper@37906
   352
	      (*"(?a | True) = True"*)
neuper@37906
   353
	      Thm ("or_false",or_false),
neuper@37906
   354
	      (*"(?a | False) = ?a"*)
neuper@37906
   355
	      Thm ("not_true",num_str not_true),
neuper@37906
   356
	      (*"(~ True) = False"*)
neuper@37906
   357
	      Thm ("not_false",num_str not_false)
neuper@37906
   358
	      (*"(~ False) = True"*)], 
neuper@37926
   359
  SOME "Vereinfache t_", 
neuper@37906
   360
  [["simplification","for_polynomials","with_parentheses"]]));
neuper@37906
   361
neuper@37906
   362
store_pbt
neuper@37906
   363
 (prep_pbt PolyMinus.thy "pbl_vereinf_poly_klammer_mal" [] e_pblID
neuper@37906
   364
 (["binom_klammer","polynom","vereinfachen"],
neuper@37906
   365
  [("#Given" ,["term t_"]),
neuper@37906
   366
   ("#Where" ,["t_ is_polyexp"]),
neuper@37906
   367
   ("#Find"  ,["normalform n_"])
neuper@37906
   368
  ],
neuper@37906
   369
  append_rls "e_rls" e_rls [(*for preds in where_*)
neuper@37906
   370
			    Calc ("Poly.is'_polyexp", eval_is_polyexp "")], 
neuper@37926
   371
  SOME "Vereinfache t_", 
neuper@37906
   372
  [["simplification","for_polynomials","with_parentheses_mult"]]));
neuper@37906
   373
neuper@37906
   374
store_pbt
neuper@37906
   375
 (prep_pbt PolyMinus.thy "pbl_probe" [] e_pblID
neuper@37906
   376
 (["probe"],
neuper@37926
   377
  [], Erls, NONE, []));
neuper@37906
   378
neuper@37906
   379
store_pbt
neuper@37906
   380
 (prep_pbt PolyMinus.thy "pbl_probe_poly" [] e_pblID
neuper@37906
   381
 (["polynom","probe"],
neuper@37906
   382
  [("#Given" ,["Pruefe e_", "mitWert ws_"]),
neuper@37906
   383
   ("#Where" ,["e_ is_polyexp"]),
neuper@37906
   384
   ("#Find"  ,["Geprueft p_"])
neuper@37906
   385
  ],
neuper@37906
   386
  append_rls "prls_pbl_probe_poly" 
neuper@37906
   387
	     e_rls [(*for preds in where_*)
neuper@37906
   388
		    Calc ("Poly.is'_polyexp", eval_is_polyexp "")], 
neuper@37926
   389
  SOME "Probe e_ ws_", 
neuper@37906
   390
  [["probe","fuer_polynom"]]));
neuper@37906
   391
neuper@37906
   392
store_pbt
neuper@37906
   393
 (prep_pbt PolyMinus.thy "pbl_probe_bruch" [] e_pblID
neuper@37906
   394
 (["bruch","probe"],
neuper@37906
   395
  [("#Given" ,["Pruefe e_", "mitWert ws_"]),
neuper@37906
   396
   ("#Where" ,["e_ is_ratpolyexp"]),
neuper@37906
   397
   ("#Find"  ,["Geprueft p_"])
neuper@37906
   398
  ],
neuper@37906
   399
  append_rls "prls_pbl_probe_bruch"
neuper@37906
   400
	     e_rls [(*for preds in where_*)
neuper@37906
   401
		    Calc ("Rational.is'_ratpolyexp", eval_is_ratpolyexp "")], 
neuper@37926
   402
  SOME "Probe e_ ws_", 
neuper@37906
   403
  [["probe","fuer_bruch"]]));
neuper@37906
   404
neuper@37906
   405
neuper@37906
   406
(** methods **)
neuper@37906
   407
neuper@37906
   408
store_met
neuper@37906
   409
    (prep_met PolyMinus.thy "met_simp_poly_minus" [] e_metID
neuper@37906
   410
	      (["simplification","for_polynomials","with_minus"],
neuper@37906
   411
	       [("#Given" ,["term t_"]),
neuper@37906
   412
		("#Where" ,["t_ is_polyexp",
neuper@37906
   413
	       "Not (matchsub (?a + (?b + ?c)) t_ | \
neuper@37906
   414
	       \     matchsub (?a + (?b - ?c)) t_ | \
neuper@37906
   415
	       \     matchsub (?a - (?b + ?c)) t_ | \
neuper@37906
   416
	       \     matchsub (?a + (?b - ?c)) t_ )"]),
neuper@37906
   417
		("#Find"  ,["normalform n_"])
neuper@37906
   418
		],
neuper@37906
   419
	       {rew_ord'="tless_true", rls' = e_rls, calc = [], srls = e_rls, 
neuper@37906
   420
		prls = append_rls "prls_met_simp_poly_minus" e_rls 
neuper@37906
   421
				  [Calc ("Poly.is'_polyexp", eval_is_polyexp ""),
neuper@37906
   422
	      Calc ("Tools.matchsub", eval_matchsub ""),
neuper@37906
   423
	      Thm ("and_true",and_true),
neuper@37906
   424
	      (*"(?a & True) = ?a"*)
neuper@37906
   425
	      Thm ("and_false",and_false),
neuper@37906
   426
	      (*"(?a & False) = False"*)
neuper@37906
   427
	      Thm ("not_true",num_str not_true),
neuper@37906
   428
	      (*"(~ True) = False"*)
neuper@37906
   429
	      Thm ("not_false",num_str not_false)
neuper@37906
   430
	      (*"(~ False) = True"*)],
neuper@37906
   431
		crls = e_rls, nrls = rls_p_33},
neuper@37906
   432
"Script SimplifyScript (t_::real) =                   \
neuper@37906
   433
\  ((Repeat((Try (Rewrite_Set ordne_alphabetisch False)) @@  \
neuper@37906
   434
\           (Try (Rewrite_Set fasse_zusammen     False)) @@  \
neuper@37906
   435
\           (Try (Rewrite_Set verschoenere       False)))) t_)"
neuper@37906
   436
	       ));
neuper@37906
   437
neuper@37906
   438
store_met
neuper@37906
   439
    (prep_met PolyMinus.thy "met_simp_poly_parenth" [] e_metID
neuper@37906
   440
	      (["simplification","for_polynomials","with_parentheses"],
neuper@37906
   441
	       [("#Given" ,["term t_"]),
neuper@37906
   442
		("#Where" ,["t_ is_polyexp"]),
neuper@37906
   443
		("#Find"  ,["normalform n_"])
neuper@37906
   444
		],
neuper@37906
   445
	       {rew_ord'="tless_true", rls' = e_rls, calc = [], srls = e_rls, 
neuper@37906
   446
		prls = append_rls "simplification_for_polynomials_prls" e_rls 
neuper@37906
   447
				  [(*for preds in where_*)
neuper@37906
   448
				   Calc("Poly.is'_polyexp",eval_is_polyexp"")],
neuper@37906
   449
		crls = e_rls, nrls = rls_p_34},
neuper@37906
   450
"Script SimplifyScript (t_::real) =                          \
neuper@37906
   451
\  ((Repeat((Try (Rewrite_Set klammern_aufloesen False)) @@  \
neuper@37906
   452
\           (Try (Rewrite_Set ordne_alphabetisch False)) @@  \
neuper@37906
   453
\           (Try (Rewrite_Set fasse_zusammen     False)) @@  \
neuper@37906
   454
\           (Try (Rewrite_Set verschoenere       False)))) t_)"
neuper@37906
   455
	       ));
neuper@37906
   456
neuper@37906
   457
store_met
neuper@37906
   458
    (prep_met PolyMinus.thy "met_simp_poly_parenth_mult" [] e_metID
neuper@37906
   459
	      (["simplification","for_polynomials","with_parentheses_mult"],
neuper@37906
   460
	       [("#Given" ,["term t_"]),
neuper@37906
   461
		("#Where" ,["t_ is_polyexp"]),
neuper@37906
   462
		("#Find"  ,["normalform n_"])
neuper@37906
   463
		],
neuper@37906
   464
	       {rew_ord'="tless_true", rls' = e_rls, calc = [], srls = e_rls, 
neuper@37906
   465
		prls = append_rls "simplification_for_polynomials_prls" e_rls 
neuper@37906
   466
				  [(*for preds in where_*)
neuper@37906
   467
				   Calc("Poly.is'_polyexp",eval_is_polyexp"")],
neuper@37906
   468
		crls = e_rls, nrls = rls_p_34},
neuper@37906
   469
"Script SimplifyScript (t_::real) =                          \
neuper@37906
   470
\  ((Repeat((Try (Rewrite_Set klammern_ausmultiplizieren False)) @@ \
neuper@37906
   471
\           (Try (Rewrite_Set discard_parentheses        False)) @@ \
neuper@37906
   472
\           (Try (Rewrite_Set ordne_monome               False)) @@ \
neuper@37906
   473
\           (Try (Rewrite_Set klammern_aufloesen         False)) @@ \
neuper@37906
   474
\           (Try (Rewrite_Set ordne_alphabetisch         False)) @@ \
neuper@37906
   475
\           (Try (Rewrite_Set fasse_zusammen             False)) @@ \
neuper@37906
   476
\           (Try (Rewrite_Set verschoenere               False)))) t_)"
neuper@37906
   477
	       ));
neuper@37906
   478
neuper@37906
   479
store_met
neuper@37906
   480
    (prep_met PolyMinus.thy "met_probe" [] e_metID
neuper@37906
   481
	      (["probe"],
neuper@37906
   482
	       [],
neuper@37906
   483
	       {rew_ord'="tless_true", rls' = e_rls, calc = [], srls = e_rls, 
neuper@37906
   484
		prls = Erls, crls = e_rls, nrls = Erls}, 
neuper@37906
   485
	       "empty_script"));
neuper@37906
   486
neuper@37906
   487
store_met
neuper@37906
   488
    (prep_met PolyMinus.thy "met_probe_poly" [] e_metID
neuper@37906
   489
	      (["probe","fuer_polynom"],
neuper@37906
   490
	       [("#Given" ,["Pruefe e_", "mitWert ws_"]),
neuper@37906
   491
		("#Where" ,["e_ is_polyexp"]),
neuper@37906
   492
		("#Find"  ,["Geprueft p_"])
neuper@37906
   493
		],
neuper@37906
   494
	       {rew_ord'="tless_true", rls' = e_rls, calc = [], srls = e_rls, 
neuper@37906
   495
		prls = append_rls "prls_met_probe_bruch"
neuper@37906
   496
				  e_rls [(*for preds in where_*)
neuper@37906
   497
					 Calc ("Rational.is'_ratpolyexp", 
neuper@37906
   498
					       eval_is_ratpolyexp "")], 
neuper@37906
   499
		crls = e_rls, nrls = rechnen}, 
neuper@37906
   500
"Script ProbeScript (e_::bool) (ws_::bool list) = \
neuper@37906
   501
\ (let e_ = Take e_;                              \
neuper@37906
   502
\      e_ = Substitute ws_ e_                     \
neuper@37906
   503
\ in (Repeat((Try (Repeat (Calculate times))) @@  \
neuper@37906
   504
\            (Try (Repeat (Calculate plus ))) @@  \
neuper@37906
   505
\            (Try (Repeat (Calculate minus))))) e_)"
neuper@37906
   506
));
neuper@37906
   507
neuper@37906
   508
store_met
neuper@37906
   509
    (prep_met PolyMinus.thy "met_probe_bruch" [] e_metID
neuper@37906
   510
	      (["probe","fuer_bruch"],
neuper@37906
   511
	       [("#Given" ,["Pruefe e_", "mitWert ws_"]),
neuper@37906
   512
		("#Where" ,["e_ is_ratpolyexp"]),
neuper@37906
   513
		("#Find"  ,["Geprueft p_"])
neuper@37906
   514
		],
neuper@37906
   515
	       {rew_ord'="tless_true", rls' = e_rls, calc = [], srls = e_rls, 
neuper@37906
   516
		prls = append_rls "prls_met_probe_bruch"
neuper@37906
   517
				  e_rls [(*for preds in where_*)
neuper@37906
   518
					 Calc ("Rational.is'_ratpolyexp", 
neuper@37906
   519
					       eval_is_ratpolyexp "")], 
neuper@37906
   520
		crls = e_rls, nrls = Erls}, 
neuper@37906
   521
	       "empty_script"));