(* WN.020812: theorems in the Reals,

   necessary for special rule sets, in addition to Isabelle2002.

   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

   !!! THIS IS THE _least_ NUMBER OF ADDITIONAL THEOREMS !!!

   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

   xxxI contain \<up> instead of ^ in the respective theorem xxx in 2002

   changed by: Richard Lang 020912

*)

theory Poly imports Simplify begin

subsection \<open>remark on term-structure of polynomials\<close>

text \<open>

WN190319:

the code below reflects missing coordination between two authors:

* ML: built the equation solver; simple rule-sets, programs; better predicates for specifications.

* MG: built simplification of polynomials with AC rewriting by ML code

WN020919:

*** there are 5 kinds of expanded normalforms ***

[1] 'complete polynomial' (Komplettes Polynom), univariate

   a_0 + a_1.x^1 +...+ a_n.x^n   not (a_n = 0)

	        not (a_n = 0), some a_i may be zero (DON'T disappear),

                variables in monomials lexicographically ordered and complete,

                x written as 1*x^1, ...

[2] 'polynomial' (Polynom), univariate and multivariate

   a_0 + a_1.x +...+ a_n.x^n   not (a_n = 0)

   a_0 + a_1.x_1.x_2^n_12...x_m^n_1m +...+  a_n.x_1^n.x_2^n_n2...x_m^n_nm

	        not (a_n = 0), some a_i may be zero (ie. monomials disappear),

                exponents and coefficients equal 1 are not (WN060904.TODO in cancel_p_)shown,

                and variables in monomials are lexicographically ordered  

   examples: [1]: "1 + (-10) * x \<up> 1 + 25 * x \<up> 2"

	     [1]: "11 + 0 * x \<up> 1 + 1 * x \<up> 2"

	     [2]: "x + (-50) * x \<up> 3"

	     [2]: "(-1) * x * y \<up> 2 + 7 * x \<up> 3"

[3] 'expanded_term' (Ausmultiplizierter Term):

   pull out unary minus to binary minus, 

   as frequently exercised in schools; other conditions for [2] hold however

   examples: "a \<up> 2 - 2 * a * b + b \<up> 2"

	     "4 * x \<up> 2 - 9 * y \<up> 2"

[4] 'polynomial_in' (Polynom in): 

   polynomial in 1 variable with arbitrary coefficients

   examples: "2 * x + (-50) * x \<up> 3"                     (poly in x)

	     "(u + v) + (2 * u \<up> 2) * a + (-u) * a \<up> 2 (poly in a)

[5] 'expanded_in' (Ausmultiplizierter Termin in): 

   analoguous to [3] with binary minus like [3]

   examples: "2 * x - 50 * x \<up> 3"                     (expanded in x)

	     "(u + v) + (2 * u \<up> 2) * a - u * a \<up> 2 (expanded in a)

\<close>

subsection \<open>consts definition for predicates in specifications\<close>

consts

  is_expanded_in :: "[real, real] => bool" ("_ is'_expanded'_in _") 

  is_poly_in     :: "[real, real] => bool" ("_ is'_poly'_in _")   (*RL DA *)

  has_degree_in  :: "[real, real] => real" ("_ has'_degree'_in _")(*RL DA *)

  is_polyrat_in  :: "[real, real] => bool" ("_ is'_polyrat'_in _")(*RL030626*)

  is_multUnordered:: "real => bool" ("_ is'_multUnordered") 

  is_addUnordered :: "real => bool" ("_ is'_addUnordered") (*WN030618*)

  is_polyexp      :: "real => bool" ("_ is'_polyexp") 

subsection \<open>theorems not yet adopted from Isabelle\<close>

axiomatization where (*.not contained in Isabelle2002,

         stated as axioms, TODO: prove as theorems;

         theorem-IDs 'xxxI' with \<up> instead of ^ in 'xxx' in Isabelle2002.*)

  realpow_pow:             "(a \<up> b) \<up> c = a \<up> (b * c)" and

  realpow_addI:            "r \<up> (n + m) = r \<up> n * r \<up> m" and

  realpow_addI_assoc_l:    "r \<up> n * (r \<up> m * s) = r \<up> (n + m) * s" and

  realpow_addI_assoc_r:    "s * r \<up> n * r \<up> m = s * r \<up> (n + m)" and

  realpow_oneI:            "r \<up> 1 = r" and

  realpow_zeroI:            "r \<up> 0 = 1" and

  realpow_eq_oneI:         "1 \<up> n = 1" and

  realpow_multI:           "(r * s) \<up> n = r \<up> n * s \<up> n"  and

  realpow_multI_poly:      "[| r is_polyexp; s is_polyexp |] ==>

			      (r * s) \<up> n = r \<up> n * s \<up> n"  and

  realpow_minus_oneI:      "(- 1) \<up> (2 * n) = 1"  and 

  real_diff_0:		         "0 - x = - (x::real)" and

  realpow_twoI:            "r \<up> 2 = r * r" and

  realpow_twoI_assoc_l:	  "r * (r * s) = r \<up> 2 * s" and

  realpow_twoI_assoc_r:	  "s * r * r = s * r \<up> 2" and

  realpow_two_atom:        "r is_atom ==> r * r = r \<up> 2" and

  realpow_plus_1:          "r * r \<up> n = r \<up> (n + 1)"   and       

  realpow_plus_1_assoc_l:  "r * (r \<up> m * s) = r \<up> (1 + m) * s"  and

  realpow_plus_1_assoc_l2: "r \<up> m * (r * s) = r \<up> (1 + m) * s"  and

  realpow_plus_1_assoc_r:  "s * r * r \<up> m = s * r \<up> (1 + m)" and

  realpow_plus_1_atom:     "r is_atom ==> r * r \<up> n = r \<up> (1 + n)" and

  realpow_def_atom:        "[| Not (r is_atom); 1 < n |]

			   ==> r \<up> n = r * r \<up> (n + -1)" and

  realpow_addI_atom:       "r is_atom ==> r \<up> n * r \<up> m = r \<up> (n + m)" and

  realpow_minus_even:	     "n is_even ==> (- r) \<up> n = r \<up> n" and

  realpow_minus_odd:       "Not (n is_even) ==> (- r) \<up> n = -1 * r \<up> n" and

(* RL 020914 *)

  real_pp_binom_times:     "(a + b)*(c + d) = a*c + a*d + b*c + b*d" and

  real_pm_binom_times:     "(a + b)*(c - d) = a*c - a*d + b*c - b*d" and

  real_mp_binom_times:     "(a - b)*(c + d) = a*c + a*d - b*c - b*d" and

  real_mm_binom_times:     "(a - b)*(c - d) = a*c - a*d - b*c + b*d" and

  real_plus_binom_pow3:    "(a + b) \<up> 3 = a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3" and

  real_plus_binom_pow3_poly: "[| a is_polyexp; b is_polyexp |] ==> 

			    (a + b) \<up> 3 = a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3" and

  real_minus_binom_pow3:   "(a - b) \<up> 3 = a \<up> 3 - 3*a \<up> 2*b + 3*a*b \<up> 2 - b \<up> 3" and

  real_minus_binom_pow3_p: "(a + -1 * b) \<up> 3 = a \<up> 3 + -3*a \<up> 2*b + 3*a*b \<up> 2 +

                           -1*b \<up> 3" and

(* real_plus_binom_pow:        "[| n is_const;  3 < n |] ==>

			       (a + b) \<up> n = (a + b) * (a + b)\<up>(n - 1)" *)

  real_plus_binom_pow4:   "(a + b) \<up> 4 = (a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3)

                           *(a + b)" and

  real_plus_binom_pow4_poly: "[| a is_polyexp; b is_polyexp |] ==> 

			   (a + b) \<up> 4 = (a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3)

                           *(a + b)" and

  real_plus_binom_pow5:    "(a + b) \<up> 5 = (a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3)

                           *(a \<up> 2 + 2*a*b + b \<up> 2)" and

  real_plus_binom_pow5_poly: "[| a is_polyexp; b is_polyexp |] ==> 

			        (a + b) \<up> 5 = (a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 

                                + b \<up> 3)*(a \<up> 2 + 2*a*b + b \<up> 2)" and

  real_diff_plus:          "a - b = a + -b" (*17.3.03: do_NOT_use*) and

  real_diff_minus:         "a - b = a + -1 * b" and

  real_plus_binom_times:   "(a + b)*(a + b) = a \<up> 2 + 2*a*b + b \<up> 2" and

  real_minus_binom_times:  "(a - b)*(a - b) = a \<up> 2 - 2*a*b + b \<up> 2" and

  (*WN071229 changed for Schaerding -----vvv*)

  (*real_plus_binom_pow2:  "(a + b) \<up> 2 = a \<up> 2 + 2*a*b + b \<up> 2"*)

  real_plus_binom_pow2:    "(a + b) \<up> 2 = (a + b) * (a + b)" and

  (*WN071229 changed for Schaerding -----\<up>*)

  real_plus_binom_pow2_poly: "[| a is_polyexp; b is_polyexp |] ==>

			       (a + b) \<up> 2 = a \<up> 2 + 2*a*b + b \<up> 2" and

  real_minus_binom_pow2:      "(a - b) \<up> 2 = a \<up> 2 - 2*a*b + b \<up> 2" and

  real_minus_binom_pow2_p:    "(a - b) \<up> 2 = a \<up> 2 + -2*a*b + b \<up> 2" and

  real_plus_minus_binom1:     "(a + b)*(a - b) = a \<up> 2 - b \<up> 2" and

  real_plus_minus_binom1_p:   "(a + b)*(a - b) = a \<up> 2 + -1*b \<up> 2" and

  real_plus_minus_binom1_p_p: "(a + b)*(a + -1 * b) = a \<up> 2 + -1*b \<up> 2" and

  real_plus_minus_binom2:     "(a - b)*(a + b) = a \<up> 2 - b \<up> 2" and

  real_plus_minus_binom2_p:   "(a - b)*(a + b) = a \<up> 2 + -1*b \<up> 2" and

  real_plus_minus_binom2_p_p: "(a + -1 * b)*(a + b) = a \<up> 2 + -1*b \<up> 2" and

  real_plus_binom_times1:     "(a +  1*b)*(a + -1*b) = a \<up> 2 + -1*b \<up> 2" and

  real_plus_binom_times2:     "(a + -1*b)*(a +  1*b) = a \<up> 2 + -1*b \<up> 2" and

  real_num_collect:           "[| l is_const; m is_const |] ==>

			      l * n + m * n = (l + m) * n" and

(* FIXME.MG.0401: replace 'real_num_collect_assoc' 

	by 'real_num_collect_assoc_l' ... are equal, introduced by MG ! *)

  real_num_collect_assoc:     "[| l is_const; m is_const |] ==> 

			      l * n + (m * n + k) = (l + m) * n + k" and

  real_num_collect_assoc_l:   "[| l is_const; m is_const |] ==>

			      l * n + (m * n + k) = (l + m)

				* n + k" and

  real_num_collect_assoc_r:   "[| l is_const; m is_const |] ==>

			      (k + m * n) + l * n = k + (l + m) * n" and

  real_one_collect:           "m is_const ==> n + m * n = (1 + m) * n" and

(* FIXME.MG.0401: replace 'real_one_collect_assoc' 

	by 'real_one_collect_assoc_l' ... are equal, introduced by MG ! *)

  real_one_collect_assoc:     "m is_const ==> n + (m * n + k) = (1 + m)* n + k" and

  real_one_collect_assoc_l:   "m is_const ==> n + (m * n + k) = (1 + m) * n + k" and

  real_one_collect_assoc_r:  "m is_const ==> (k + n) +  m * n = k + (1 + m) * n" and

(* FIXME.MG.0401: replace 'real_mult_2_assoc' 

	by 'real_mult_2_assoc_l' ... are equal, introduced by MG ! *)

  real_mult_2_assoc:          "z1 + (z1 + k) = 2 * z1 + k" and

  real_mult_2_assoc_l:        "z1 + (z1 + k) = 2 * z1 + k" and

  real_mult_2_assoc_r:        "(k + z1) + z1 = k + 2 * z1" and

  real_mult_left_commute: "z1 * (z2 * z3) = z2 * (z1 * z3)" and

  real_mult_minus1:       "-1 * z = - (z::real)" and

  real_mult_2:            "2 * z = z + (z::real)" and

  real_add_mult_distrib_poly: "w is_polyexp ==> (z1 + z2) * w = z1 * w + z2 * w" and

  real_add_mult_distrib2_poly:"w is_polyexp ==> w * (z1 + z2) = w * z1 + w * z2"

subsection \<open>auxiliary functions\<close>

ML \<open>

val thy = @{theory};

val poly_consts =

  ["Groups.plus_class.plus", "Groups.minus_class.minus",

  "Rings.divide_class.divide", "Groups.times_class.times",

  "Transcendental.powr"];

\<close>

subsubsection \<open>for predicates in specifications (ML)\<close>

ML \<open>

(*--- auxiliary for is_expanded_in, is_poly_in, has_degree_in ---*)

(*. a 'monomial t in variable v' is a term t with

  either (1) v NOT existent in t, or (2) v contained in t,

  if (1) then degree 0

  if (2) then v is a factor on the very right, ev. with exponent.*)

fun factor_right_deg (*case 2*)

	    (Const ("Groups.times_class.times", _) $ t1 $ (Const ("Transcendental.powr",_) $ vv $ Free (d, _))) v =

	   if vv = v andalso not (Prog_Expr.occurs_in v t1) then SOME (TermC.int_of_str d) else NONE

  | factor_right_deg (Const ("Transcendental.powr",_) $ vv $ Free (d,_)) v =

	   if (vv = v) then SOME (TermC.int_of_str d) else NONE

  | factor_right_deg (Const ("Groups.times_class.times",_) $ t1 $ vv) v = 

	   if vv = v andalso not (Prog_Expr.occurs_in v t1) then SOME 1 else NONE

  | factor_right_deg vv v =

	  if (vv = v) then SOME 1 else NONE;    

fun mono_deg_in m v =  (*case 1*)

	if not (Prog_Expr.occurs_in v m) then (*case 1*) SOME 0 else factor_right_deg m v;

fun expand_deg_in t v =

	let

    fun edi ~1 ~1 (Const ("Groups.plus_class.plus", _) $ t1 $ t2) =

          (case mono_deg_in t2 v of (* $ is left associative*)

            SOME d' => edi d' d' t1 | NONE => NONE)

      | edi ~1 ~1 (Const ("Groups.minus_class.minus", _) $ t1 $ t2) =

          (case mono_deg_in t2 v of

            SOME d' => edi d' d' t1 | NONE => NONE)

      | edi d dmax (Const ("Groups.minus_class.minus", _) $ t1 $ t2) =

          (case mono_deg_in t2 v of (*(d = 0 andalso d' = 0) handle 3+4-...4 +x*)

	        SOME d' => if d > d' orelse (d = 0 andalso d' = 0) then edi d' dmax t1 else NONE

          | NONE => NONE)

      | edi d dmax (Const ("Groups.plus_class.plus",_) $ t1 $ t2) =

          (case mono_deg_in t2 v of

            SOME d' =>    (*RL (d = 0 andalso d' = 0) need to handle 3+4-...4 +x*)

              if d > d' orelse (d = 0 andalso d' = 0) then edi d' dmax t1 else NONE

          | NONE => NONE)

      | edi ~1 ~1 t =

          (case mono_deg_in t v of d as SOME _ => d | NONE => NONE)

      | edi d dmax t = (*basecase last*)

    	    (case mono_deg_in t v of

    	      SOME d' => if d > d' orelse (d = 0 andalso d' = 0) then SOME dmax else NONE

		      | NONE => NONE)

	in edi ~1 ~1 t end;

fun poly_deg_in t v =

	let

    fun edi ~1 ~1 (Const ("Groups.plus_class.plus",_) $ t1 $ t2) =

		    (case mono_deg_in t2 v of (* $ is left associative *)

		      SOME d' => edi d' d' t1

        | NONE => NONE)

	    | edi d dmax (Const ("Groups.plus_class.plus",_) $ t1 $ t2) =

		    (case mono_deg_in t2 v of

	        SOME d' =>    (*RL (d = 0 andalso (d' = 0)) handle 3+4-...4 +x*)

            if d > d' orelse (d = 0 andalso d' = 0) then edi d' dmax t1 else NONE

        | NONE => NONE)

	    | edi ~1 ~1 t =

        (case mono_deg_in t v of

		      d as SOME _ => d

        | NONE => NONE)

	    | edi d dmax t = (* basecase last *)

		    (case mono_deg_in t v of

		      SOME d' =>

            if d > d' orelse (d = 0 andalso d' = 0) then SOME dmax else NONE

        | NONE => NONE)

	in edi ~1 ~1 t end;

\<close>

subsubsection \<open>for hard-coded AC rewriting (MG)\<close>

ML \<open>

(**. MG.03: make_polynomial_ ... uses SML-fun for ordering .**)

(*FIXME.0401: make SML-order local to make_polynomial(_) *)

(*FIXME.0401: replace 'make_polynomial'(old) by 'make_polynomial_'(MG) *)

(* Polynom --> List von Monomen *) 

fun poly2list (Const ("Groups.plus_class.plus",_) $ t1 $ t2) = 

    (poly2list t1) @ (poly2list t2)

  | poly2list t = [t];

(* Monom --> Liste von Variablen *)

fun monom2list (Const ("Groups.times_class.times",_) $ t1 $ t2) = 

    (monom2list t1) @ (monom2list t2)

  | monom2list t = [t];

(* liefert Variablenname (String) einer Variablen und Basis bei Potenz *)

fun get_basStr (Const ("Transcendental.powr",_) $ Free (str, _) $ _) = str

  | get_basStr (Free (str, _)) = str

  | get_basStr _ = "|||"; (* gross gewichtet; für Brüch ect. *)

(*| get_basStr t = 

    raise ERROR("get_basStr: called with t= "^(UnparseC.term t));*)

(* liefert Hochzahl (String) einer Variablen bzw Gewichtstring (zum Sortieren) *)

fun get_potStr (Const ("Transcendental.powr",_) $ Free _ $ Free (str, _)) = str

  | get_potStr (Const ("Transcendental.powr",_) $ Free _ $ _ ) = "|||" (* gross gewichtet *)

  | get_potStr (Free (_, _)) = "---" (* keine Hochzahl --> kleinst gewichtet *)

  | get_potStr _ = "||||||"; (* gross gewichtet; für Brüch ect. *)

(*| get_potStr t = 

    raise ERROR("get_potStr: called with t= "^(UnparseC.term t));*)

(* Umgekehrte string_ord *)

val string_ord_rev =  rev_order o string_ord;

 (* Ordnung zum lexikographischen Vergleich zweier Variablen (oder Potenzen) 

    innerhalb eines Monomes:

    - zuerst lexikographisch nach Variablenname 

    - wenn gleich: nach steigender Potenz *)

fun var_ord (a,b: term) = prod_ord string_ord string_ord 

    ((get_basStr a, get_potStr a), (get_basStr b, get_potStr b));

(* Ordnung zum lexikographischen Vergleich zweier Variablen (oder Potenzen); 

   verwendet zum Sortieren von Monomen mittels Gesamtgradordnung:

   - zuerst lexikographisch nach Variablenname 

   - wenn gleich: nach sinkender Potenz*)

fun var_ord_revPow (a,b: term) = prod_ord string_ord string_ord_rev 

    ((get_basStr a, get_potStr a), (get_basStr b, get_potStr b));

(* Ordnet ein Liste von Variablen (und Potenzen) lexikographisch *)

val sort_varList = sort var_ord;

(* Entfernet aeussersten Operator (Wurzel) aus einem Term und schreibt 

   Argumente in eine Liste *)

fun args u : term list =

    let fun stripc (f$t, ts) = stripc (f, t::ts)

	  | stripc (t as Free _, ts) = (t::ts)

	  | stripc (_, ts) = ts

    in stripc (u, []) end;

(* liefert True, falls der Term (Liste von Termen) nur Zahlen 

   (keine Variablen) enthaelt *)

fun filter_num [] = true

  | filter_num [Free x] = if (TermC.is_num (Free x)) then true

				else false

  | filter_num ((Free _)::_) = false

  | filter_num ts =

    (filter_num o (filter_out TermC.is_num) o flat o (map args)) ts;

(* liefert True, falls der Term nur Zahlen (keine Variablen) enthaelt 

   dh. er ist ein numerischer Wert und entspricht einem Koeffizienten *)

fun is_nums t = filter_num [t];

(* Berechnet den Gesamtgrad eines Monoms *)

local 

    fun counter (n, []) = n

      | counter (n, x :: xs) = 

	if (is_nums x) then

	    counter (n, xs) 

	else 

	    (case x of 

		 (Const ("Transcendental.powr", _) $ Free _ $ Free (str_h, T)) => 

		     if (is_nums (Free (str_h, T))) then

			 counter (n + (the (TermC.int_opt_of_string str_h)), xs)

		     else counter (n + 1000, xs) (*FIXME.MG?!*)

	       | (Const ("Transcendental.powr", _) $ Free _ $ _ ) => 

		     counter (n + 1000, xs) (*FIXME.MG?!*)

	       | (Free _) => counter (n + 1, xs)

	     (*| _ => raise ERROR("monom_degree: called with factor: "^(UnparseC.term x)))*)

	       | _ => counter (n + 10000, xs)) (*FIXME.MG?! ... Brüche ect.*)

in  

    fun monom_degree l = counter (0, l) 

end;(*local*)

(* wie Ordnung dict_ord (lexicographische Ordnung zweier Listen, mit Vergleich 

   der Listen-Elemente mit elem_ord) - Elemente die Bedingung cond erfuellen, 

   werden jedoch dabei ignoriert (uebersprungen)  *)

fun dict_cond_ord _ _ ([], []) = EQUAL

  | dict_cond_ord _ _ ([], _ :: _) = LESS

  | dict_cond_ord _ _ (_ :: _, []) = GREATER

  | dict_cond_ord elem_ord cond (x :: xs, y :: ys) =

    (case (cond x, cond y) of 

	 (false, false) => (case elem_ord (x, y) of 

				EQUAL => dict_cond_ord elem_ord cond (xs, ys) 

			      | ord => ord)

       | (false, true)  => dict_cond_ord elem_ord cond (x :: xs, ys)

       | (true, false)  => dict_cond_ord elem_ord cond (xs, y :: ys)

       | (true, true)  =>  dict_cond_ord elem_ord cond (xs, ys) );

(* Gesamtgradordnung zum Vergleich von Monomen (Liste von Variablen/Potenzen):

   zuerst nach Gesamtgrad, bei gleichem Gesamtgrad lexikographisch ordnen - 

   dabei werden Koeffizienten ignoriert (2*3*a \<up> 2*4*b gilt wie a \<up> 2*b) *)

fun degree_ord (xs, ys) =

	    prod_ord int_ord (dict_cond_ord var_ord_revPow is_nums) 

	    ((monom_degree xs, xs), (monom_degree ys, ys));

fun hd_str str = substring (str, 0, 1);

fun tl_str str = substring (str, 1, (size str) - 1);

(* liefert nummerischen Koeffizienten eines Monoms oder NONE *)

fun get_koeff_of_mon [] =  raise ERROR("get_koeff_of_mon: called with l = []")

  | get_koeff_of_mon (x::_) = if is_nums x then SOME x else NONE;

(* wandelt Koeffizient in (zum sortieren geeigneten) String um *)

fun koeff2ordStr (SOME x) = (case x of 

				 (Free (str, _)) => 

				     if (hd_str str) = "-" then (tl_str str)^"0" (* 3 < -3 *)

				     else str

			       | _ => "aaa") (* "num.Ausdruck" --> gross *)

  | koeff2ordStr NONE = "---"; (* "kein Koeff" --> kleinste *)

(* Order zum Vergleich von Koeffizienten (strings): 

   "kein Koeff" < "0" < "1" < "-1" < "2" < "-2" < ... < "num.Ausdruck" *)

fun compare_koeff_ord (xs, ys) = 

    string_ord ((koeff2ordStr o get_koeff_of_mon) xs,

		(koeff2ordStr o get_koeff_of_mon) ys);

(* Gesamtgradordnung degree_ord + Ordnen nach Koeffizienten falls EQUAL *)

fun koeff_degree_ord (xs, ys) =

	    prod_ord degree_ord compare_koeff_ord ((xs, xs), (ys, ys));

(* Ordnet ein Liste von Monomen (Monom = Liste von Variablen) mittels 

   Gesamtgradordnung *)

val sort_monList = sort koeff_degree_ord;

(* Alternativ zu degree_ord koennte auch die viel einfachere und 

   kuerzere Ordnung simple_ord verwendet werden - ist aber nicht 

   fuer unsere Zwecke geeignet!

fun simple_ord (al,bl: term list) = dict_ord string_ord 

	 (map get_basStr al, map get_basStr bl); 

val sort_monList = sort simple_ord; *)

(* aus 2 Variablen wird eine Summe bzw ein Produkt erzeugt 

   (mit gewuenschtem Typen T) *)

fun plus T = Const ("Groups.plus_class.plus", [T,T] ---> T);

fun mult T = Const ("Groups.times_class.times", [T,T] ---> T);

fun binop op_ t1 t2 = op_ $ t1 $ t2;

fun create_prod T (a,b) = binop (mult T) a b;

fun create_sum T (a,b) = binop (plus T) a b;

(* löscht letztes Element einer Liste *)

fun drop_last l = take ((length l)-1,l);

(* Liste von Variablen --> Monom *)

fun create_monom T vl = foldr (create_prod T) (drop_last vl, last_elem vl);

(* Bemerkung: 

   foldr bewirkt rechtslastige Klammerung des Monoms - ist notwendig, damit zwei 

   gleiche Monome zusammengefasst werden können (collect_numerals)! 

   zB: 2*(x*(y*z)) + 3*(x*(y*z)) --> (2+3)*(x*(y*z))*)

(* Liste von Monomen --> Polynom *)	

fun create_polynom T ml = foldl (create_sum T) (hd ml, tl ml);

(* Bemerkung: 

   foldl bewirkt linkslastige Klammerung des Polynoms (der Summanten) - 

   bessere Darstellung, da keine Klammern sichtbar! 

   (und discard_parentheses in make_polynomial hat weniger zu tun) *)

(* sorts the variables (faktors) of an expanded polynomial lexicographical *)

fun sort_variables t = 

    let

	val ll =  map monom2list (poly2list t);

	val lls = map sort_varList ll; 

	val T = type_of t;

	val ls = map (create_monom T) lls;

    in create_polynom T ls end;

(* sorts the monoms of an expanded and variable-sorted polynomial 

   by total_degree *)

fun sort_monoms t = 

    let

	val ll =  map monom2list (poly2list t);

	val lls = sort_monList ll;

	val T = type_of t;

	val ls = map (create_monom T) lls;

    in create_polynom T ls end;

\<close>

subsubsection \<open>rewrite order for hard-coded AC rewriting\<close>

ML \<open>

local (*. for make_polynomial .*)

open Term;  (* for type order = EQUAL | LESS | GREATER *)

fun pr_ord EQUAL = "EQUAL"

  | pr_ord LESS  = "LESS"

  | pr_ord GREATER = "GREATER";

fun dest_hd' (Const (a, T)) =                          (* ~ term.ML *)

  (case a of

     "Transcendental.powr" => ((("|||||||||||||", 0), T), 0)    (*WN greatest string*)

   | _ => (((a, 0), T), 0))

  | dest_hd' (Free (a, T)) = (((a, 0), T), 1)

  | dest_hd' (Var v) = (v, 2)

  | dest_hd' (Bound i) = ((("", i), dummyT), 3)

  | dest_hd' (Abs (_, T, _)) = ((("", 0), T), 4)

  | dest_hd' t = raise TERM ("dest_hd'", [t]);

fun size_of_term' (Const(str,_) $ t) =

  if "Transcendental.powr"= str then 1000 + size_of_term' t else 1+size_of_term' t(*WN*)

  | size_of_term' (Abs (_,_,body)) = 1 + size_of_term' body

  | size_of_term' (f$t) = size_of_term' f  +  size_of_term' t

  | size_of_term' _ = 1;

fun term_ord' pr thy (Abs (_, T, t), Abs(_, U, u)) =       (* ~ term.ML *)

    (case term_ord' pr thy (t, u) of EQUAL => Term_Ord.typ_ord (T, U) | ord => ord)

  | term_ord' pr thy (t, u) =

    (if pr then 

	   let

       val (f, ts) = strip_comb t and (g, us) = strip_comb u;

       val _ = tracing ("t= f@ts= \"" ^ UnparseC.term_in_thy thy f ^ "\" @ \"[" ^

         commas (map (UnparseC.term_in_thy thy) ts) ^ "]\"");

       val _ = tracing("u= g@us= \"" ^ UnparseC.term_in_thy thy g ^ "\" @ \"[" ^

         commas (map (UnparseC.term_in_thy thy) us) ^ "]\"");

       val _ = tracing ("size_of_term(t,u)= (" ^ string_of_int (size_of_term' t) ^ ", " ^

         string_of_int (size_of_term' u) ^ ")");

       val _ = tracing ("hd_ord(f,g)      = " ^ (pr_ord o hd_ord) (f,g));

       val _ = tracing ("terms_ord(ts,us) = " ^ (pr_ord o terms_ord str false) (ts, us));

       val _ = tracing ("-------");

     in () end

       else ();

	 case int_ord (size_of_term' t, size_of_term' u) of

	   EQUAL =>

	     let val (f, ts) = strip_comb t and (g, us) = strip_comb u in

	       (case hd_ord (f, g) of EQUAL => (terms_ord str pr) (ts, us) 

	     | ord => ord)

	     end

	 | ord => ord)

and hd_ord (f, g) =                                        (* ~ term.ML *)

  prod_ord (prod_ord Term_Ord.indexname_ord Term_Ord.typ_ord) int_ord (dest_hd' f, dest_hd' g)

and terms_ord _ pr (ts, us) = 

    list_ord (term_ord' pr (ThyC.get_theory "Isac_Knowledge"))(ts, us);

in

fun ord_make_polynomial (pr:bool) thy (_: subst) tu = 

    (term_ord' pr thy(***) tu = LESS );

end;(*local*)

Rewrite_Ord.rew_ord' := overwritel (! Rewrite_Ord.rew_ord', (* TODO: make analogous to KEStore_Elems.add_mets *)

[("termlessI", termlessI), ("ord_make_polynomial", ord_make_polynomial false thy)]);

\<close>

subsection \<open>predicates\<close>

subsubsection \<open>in specifications\<close>

ML \<open>

(* is_polyrat_in becomes true, if no bdv is in the denominator of a fraction*)

fun is_polyrat_in t v = 

  let

   	fun finddivide (_ $ _ $ _ $ _) _ = raise ERROR("is_polyrat_in:")

	    (* at the moment there is no term like this, but ....*)

	  | finddivide (Const ("Rings.divide_class.divide",_) $ _ $ b) v = not (Prog_Expr.occurs_in v b)

	  | finddivide (_ $ t1 $ t2) v = finddivide t1 v orelse finddivide t2 v

	  | finddivide (_ $ t1) v = finddivide t1 v

	  | finddivide _ _ = false;

  in finddivide t v end;

fun is_expanded_in t v = case expand_deg_in t v of SOME _ => true | NONE => false;

fun is_poly_in t v =     case poly_deg_in t v of SOME _ => true | NONE => false;

fun has_degree_in t v =  case expand_deg_in t v of SOME d => d | NONE => ~1;

(*.the expression contains + - * ^ only ?

   this is weaker than 'is_polynomial' !.*)

fun is_polyexp (Free _) = true

  | is_polyexp (Const _) = true (* potential danger: bdv is not considered *)

  | is_polyexp (Const ("Groups.plus_class.plus",_) $ Free _ $ Free _) = true

  | is_polyexp (Const ("Groups.minus_class.minus",_) $ Free _ $ Free _) = true

  | is_polyexp (Const ("Groups.times_class.times",_) $ Free _ $ Free _) = true

  | is_polyexp (Const ("Transcendental.powr",_) $ Free _ $ Free _) = true

  | is_polyexp (Const ("Groups.plus_class.plus",_) $ t1 $ t2) = 

               ((is_polyexp t1) andalso (is_polyexp t2))

  | is_polyexp (Const ("Groups.minus_class.minus",_) $ t1 $ t2) = 

               ((is_polyexp t1) andalso (is_polyexp t2))

  | is_polyexp (Const ("Groups.times_class.times",_) $ t1 $ t2) = 

               ((is_polyexp t1) andalso (is_polyexp t2))

  | is_polyexp (Const ("Transcendental.powr",_) $ t1 $ t2) = 

               ((is_polyexp t1) andalso (is_polyexp t2))

  | is_polyexp _ = false;

\<close>

subsubsection \<open>for hard-coded AC rewriting\<close>

ML \<open>

(* auch Klammerung muss übereinstimmen;

   sort_variables klammert Produkte rechtslastig*)

fun is_multUnordered t = ((is_polyexp t) andalso not (t = sort_variables t));

fun is_addUnordered t = ((is_polyexp t) andalso not (t = sort_monoms t));

\<close>

subsection \<open>evaluations functions\<close>

subsubsection \<open>for predicates\<close>

ML \<open>

fun eval_is_polyrat_in _ _(p as (Const ("Poly.is_polyrat_in",_) $ t $ v)) _  =

    if is_polyrat_in t v 

    then SOME ((UnparseC.term p) ^ " = True",

	        HOLogic.Trueprop $ (TermC.mk_equality (p, @{term True})))

    else SOME ((UnparseC.term p) ^ " = True",

	        HOLogic.Trueprop $ (TermC.mk_equality (p, @{term False})))

  | eval_is_polyrat_in _ _ _ _ = ((*tracing"### no matches";*) NONE);

(*("is_expanded_in", ("Poly.is_expanded_in", eval_is_expanded_in ""))*)

fun eval_is_expanded_in _ _ 

       (p as (Const ("Poly.is_expanded_in",_) $ t $ v)) _ =

    if is_expanded_in t v

    then SOME ((UnparseC.term p) ^ " = True",

	        HOLogic.Trueprop $ (TermC.mk_equality (p, @{term True})))

    else SOME ((UnparseC.term p) ^ " = True",

	        HOLogic.Trueprop $ (TermC.mk_equality (p, @{term False})))

  | eval_is_expanded_in _ _ _ _ = NONE;

(*("is_poly_in", ("Poly.is_poly_in", eval_is_poly_in ""))*)

fun eval_is_poly_in _ _ 

       (p as (Const ("Poly.is_poly_in",_) $ t $ v)) _ =

    if is_poly_in t v

    then SOME ((UnparseC.term p) ^ " = True",

	        HOLogic.Trueprop $ (TermC.mk_equality (p, @{term True})))

    else SOME ((UnparseC.term p) ^ " = True",

	        HOLogic.Trueprop $ (TermC.mk_equality (p, @{term False})))

  | eval_is_poly_in _ _ _ _ = NONE;

(*("has_degree_in", ("Poly.has_degree_in", eval_has_degree_in ""))*)

fun eval_has_degree_in _ _ 

	     (p as (Const ("Poly.has_degree_in",_) $ t $ v)) _ =

    let val d = has_degree_in t v

	val d' = TermC.term_of_num HOLogic.realT d

    in SOME ((UnparseC.term p) ^ " = " ^ (string_of_int d),

	      HOLogic.Trueprop $ (TermC.mk_equality (p, d')))

    end

  | eval_has_degree_in _ _ _ _ = NONE;

(*("is_polyexp", ("Poly.is_polyexp", eval_is_polyexp ""))*)

fun eval_is_polyexp (thmid:string) _ 

		       (t as (Const("Poly.is_polyexp", _) $ arg)) thy = 

    if is_polyexp arg

    then SOME (TermC.mk_thmid thmid (UnparseC.term_in_thy thy arg) "", 

	         HOLogic.Trueprop $ (TermC.mk_equality (t, @{term True})))

    else SOME (TermC.mk_thmid thmid (UnparseC.term_in_thy thy arg) "", 

	         HOLogic.Trueprop $ (TermC.mk_equality (t, @{term False})))

  | eval_is_polyexp _ _ _ _ = NONE; 

\<close>

subsubsection \<open>for hard-coded AC rewriting\<close>

ML \<open>

(*WN.18.6.03 *)

(*("is_addUnordered", ("Poly.is_addUnordered", eval_is_addUnordered ""))*)

fun eval_is_addUnordered (thmid:string) _ 

		       (t as (Const("Poly.is_addUnordered", _) $ arg)) thy = 

    if is_addUnordered arg

    then SOME (TermC.mk_thmid thmid (UnparseC.term_in_thy thy arg) "", 

	         HOLogic.Trueprop $ (TermC.mk_equality (t, @{term True})))

    else SOME (TermC.mk_thmid thmid (UnparseC.term_in_thy thy arg) "", 

	         HOLogic.Trueprop $ (TermC.mk_equality (t, @{term False})))

  | eval_is_addUnordered _ _ _ _ = NONE; 

fun eval_is_multUnordered (thmid:string) _ 

		       (t as (Const("Poly.is_multUnordered", _) $ arg)) thy = 

    if is_multUnordered arg

    then SOME (TermC.mk_thmid thmid (UnparseC.term_in_thy thy arg) "", 

	         HOLogic.Trueprop $ (TermC.mk_equality (t, @{term True})))

    else SOME (TermC.mk_thmid thmid (UnparseC.term_in_thy thy arg) "", 

	         HOLogic.Trueprop $ (TermC.mk_equality (t, @{term False})))

  | eval_is_multUnordered _ _ _ _ = NONE; 

\<close>

setup \<open>KEStore_Elems.add_calcs

  [("is_polyrat_in", ("Poly.is_polyrat_in",

		    eval_is_polyrat_in "#eval_is_polyrat_in")),

    ("is_expanded_in", ("Poly.is_expanded_in", eval_is_expanded_in "")),

    ("is_poly_in", ("Poly.is_poly_in", eval_is_poly_in "")),

    ("has_degree_in", ("Poly.has_degree_in", eval_has_degree_in "")),

    ("is_polyexp", ("Poly.is_polyexp", eval_is_polyexp "")),

    ("is_multUnordered", ("Poly.is_multUnordered", eval_is_multUnordered"")),

    ("is_addUnordered", ("Poly.is_addUnordered", eval_is_addUnordered ""))]\<close>

subsection \<open>rule-sets\<close>

subsubsection \<open>without specific order\<close>

ML \<open>

(* used only for merge *)

val calculate_Poly = Rule_Set.append_rules "calculate_PolyFIXXXME.not.impl." Rule_Set.empty [];

(*.for evaluation of conditions in rewrite rules.*)

val Poly_erls = Rule_Set.append_rules "Poly_erls" Atools_erls

  [Rule.Eval ("HOL.eq", Prog_Expr.eval_equal "#equal_"),

  Rule.Thm  ("real_unari_minus", ThmC.numerals_to_Free @{thm real_unari_minus}),

  Rule.Eval ("Groups.plus_class.plus", eval_binop "#add_"),

  Rule.Eval ("Groups.minus_class.minus", eval_binop "#sub_"),

  Rule.Eval ("Groups.times_class.times", eval_binop "#mult_"),

  Rule.Eval ("Transcendental.powr", eval_binop "#power_")];

val poly_crls = Rule_Set.append_rules "poly_crls" Atools_crls

  [Rule.Eval ("HOL.eq", Prog_Expr.eval_equal "#equal_"),

  Rule.Thm ("real_unari_minus", ThmC.numerals_to_Free @{thm real_unari_minus}),

  Rule.Eval ("Groups.plus_class.plus", eval_binop "#add_"),

  Rule.Eval ("Groups.minus_class.minus", eval_binop "#sub_"),

  Rule.Eval ("Groups.times_class.times", eval_binop "#mult_"),

  Rule.Eval ("Transcendental.powr" , eval_binop "#power_")];

\<close>

ML \<open>

val expand =

  Rule_Def.Repeat {id = "expand", preconds = [], rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty,srls = Rule_Set.Empty, calc = [], errpatts = [],

      rules = [Rule.Thm ("distrib_right" , ThmC.numerals_to_Free @{thm distrib_right}),

	       (*"(z1.0 + z2.0) * w = z1.0 * w + z2.0 * w"*)

	       Rule.Thm ("distrib_left", ThmC.numerals_to_Free @{thm distrib_left})

	       (*"w * (z1.0 + z2.0) = w * z1.0 + w * z2.0"*)

	       ], scr = Rule.Empty_Prog};

val discard_minus =

  Rule_Def.Repeat {id = "discard_minus", preconds = [], rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty, srls = Rule_Set.Empty, calc = [], errpatts = [],

      rules =

       [Rule.Thm ("real_diff_minus", ThmC.numerals_to_Free @{thm real_diff_minus}),

          (*"a - b = a + -1 * b"*)

	        Rule.Thm ("sym_real_mult_minus1", ThmC.numerals_to_Free (@{thm real_mult_minus1} RS @{thm sym}))

	          (*- ?z = "-1 * ?z"*)],

	      scr = Rule.Empty_Prog};

val expand_poly_ = 

  Rule_Def.Repeat{id = "expand_poly_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty,srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      rules =

        [Rule.Thm ("real_plus_binom_pow4", ThmC.numerals_to_Free @{thm real_plus_binom_pow4}),

	           (*"(a + b) \<up> 4 = ... "*)

	         Rule.Thm ("real_plus_binom_pow5",ThmC.numerals_to_Free @{thm real_plus_binom_pow5}),

	           (*"(a + b) \<up> 5 = ... "*)

	         Rule.Thm ("real_plus_binom_pow3",ThmC.numerals_to_Free @{thm real_plus_binom_pow3}),

	           (*"(a + b) \<up> 3 = a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3" *)

	         (*WN071229 changed/removed for Schaerding -----vvv*)

	         (*Rule.Thm ("real_plus_binom_pow2",ThmC.numerals_to_Free @{thm real_plus_binom_pow2}),*)

	           (*"(a + b) \<up> 2 = a \<up> 2 + 2*a*b + b \<up> 2"*)

	         Rule.Thm ("real_plus_binom_pow2",ThmC.numerals_to_Free @{thm real_plus_binom_pow2}),

	           (*"(a + b) \<up> 2 = (a + b) * (a + b)"*)

	         (*Rule.Thm ("real_plus_minus_binom1_p_p", ThmC.numerals_to_Free @{thm real_plus_minus_binom1_p_p}),*)

	           (*"(a + b)*(a + -1 * b) = a \<up> 2 + -1*b \<up> 2"*)

	         (*Rule.Thm ("real_plus_minus_binom2_p_p", ThmC.numerals_to_Free @{thm real_plus_minus_binom2_p_p}),*)

	           (*"(a + -1 * b)*(a + b) = a \<up> 2 + -1*b \<up> 2"*)

	         (*WN071229 changed/removed for Schaerding -----\<up>*)

	         Rule.Thm ("distrib_right" ,ThmC.numerals_to_Free @{thm distrib_right}),

	           (*"(z1.0 + z2.0) * w = z1.0 * w + z2.0 * w"*)

	         Rule.Thm ("distrib_left",ThmC.numerals_to_Free @{thm distrib_left}),

	           (*"w * (z1.0 + z2.0) = w * z1.0 + w * z2.0"*)

	         Rule.Thm ("realpow_multI", ThmC.numerals_to_Free @{thm realpow_multI}),

	           (*"(r * s) \<up> n = r \<up> n * s \<up> n"*)

	         Rule.Thm ("realpow_pow",ThmC.numerals_to_Free @{thm realpow_pow})

	           (*"(a \<up> b) \<up> c = a \<up> (b * c)"*)

	       ], scr = Rule.Empty_Prog};

val expand_poly_rat_ = 

  Rule_Def.Repeat{id = "expand_poly_rat_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls =  Rule_Set.append_rules "Rule_Set.empty-is_polyexp" Rule_Set.empty

	        [Rule.Eval ("Poly.is_polyexp", eval_is_polyexp "")

		 ],

      srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      rules = 

        [Rule.Thm ("real_plus_binom_pow4_poly", ThmC.numerals_to_Free @{thm real_plus_binom_pow4_poly}),

	     (*"[| a is_polyexp; b is_polyexp |] ==> (a + b) \<up> 4 = ... "*)

	 Rule.Thm ("real_plus_binom_pow5_poly", ThmC.numerals_to_Free @{thm real_plus_binom_pow5_poly}),

	     (*"[| a is_polyexp; b is_polyexp |] ==> (a + b) \<up> 5 = ... "*)

	 Rule.Thm ("real_plus_binom_pow2_poly",ThmC.numerals_to_Free @{thm real_plus_binom_pow2_poly}),

	     (*"[| a is_polyexp; b is_polyexp |] ==>

		            (a + b) \<up> 2 = a \<up> 2 + 2*a*b + b \<up> 2"*)

	 Rule.Thm ("real_plus_binom_pow3_poly",ThmC.numerals_to_Free @{thm real_plus_binom_pow3_poly}),

	     (*"[| a is_polyexp; b is_polyexp |] ==> 

			(a + b) \<up> 3 = a \<up> 3 + 3*a \<up> 2*b + 3*a*b \<up> 2 + b \<up> 3" *)

	 Rule.Thm ("real_plus_minus_binom1_p_p",ThmC.numerals_to_Free @{thm real_plus_minus_binom1_p_p}),

	     (*"(a + b)*(a + -1 * b) = a \<up> 2 + -1*b \<up> 2"*)

	 Rule.Thm ("real_plus_minus_binom2_p_p",ThmC.numerals_to_Free @{thm real_plus_minus_binom2_p_p}),

	     (*"(a + -1 * b)*(a + b) = a \<up> 2 + -1*b \<up> 2"*)

	 Rule.Thm ("real_add_mult_distrib_poly",

               ThmC.numerals_to_Free @{thm real_add_mult_distrib_poly}),

	       (*"w is_polyexp ==> (z1.0 + z2.0) * w = z1.0 * w + z2.0 * w"*)

	 Rule.Thm("real_add_mult_distrib2_poly",

              ThmC.numerals_to_Free @{thm real_add_mult_distrib2_poly}),

	     (*"w is_polyexp ==> w * (z1.0 + z2.0) = w * z1.0 + w * z2.0"*)

	 Rule.Thm ("realpow_multI_poly", ThmC.numerals_to_Free @{thm realpow_multI_poly}),

	     (*"[| r is_polyexp; s is_polyexp |] ==> 

		            (r * s) \<up> n = r \<up> n * s \<up> n"*)

	  Rule.Thm ("realpow_pow",ThmC.numerals_to_Free @{thm realpow_pow})

	      (*"(a \<up> b) \<up> c = a \<up> (b * c)"*)

	 ], scr = Rule.Empty_Prog};

val simplify_power_ = 

  Rule_Def.Repeat{id = "simplify_power_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty, srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      rules = [(*MG: Reihenfolge der folgenden 2 Rule.Thm muss so bleiben, wegen

		a*(a*a) --> a*a \<up> 2 und nicht a*(a*a) --> a \<up> 2*a *)

	       Rule.Thm ("sym_realpow_twoI",

                     ThmC.numerals_to_Free (@{thm realpow_twoI} RS @{thm sym})),	

	       (*"r * r = r \<up> 2"*)

	       Rule.Thm ("realpow_twoI_assoc_l",ThmC.numerals_to_Free @{thm realpow_twoI_assoc_l}),

	       (*"r * (r * s) = r \<up> 2 * s"*)

	       Rule.Thm ("realpow_plus_1",ThmC.numerals_to_Free @{thm realpow_plus_1}),		

	       (*"r * r \<up> n = r \<up> (n + 1)"*)

	       Rule.Thm ("realpow_plus_1_assoc_l",

                     ThmC.numerals_to_Free @{thm realpow_plus_1_assoc_l}),

	       (*"r * (r \<up> m * s) = r \<up> (1 + m) * s"*)

	       (*MG 9.7.03: neues Rule.Thm wegen a*(a*(a*b)) --> a \<up> 2*(a*b) *)

	       Rule.Thm ("realpow_plus_1_assoc_l2",

                     ThmC.numerals_to_Free @{thm realpow_plus_1_assoc_l2}),

	       (*"r \<up> m * (r * s) = r \<up> (1 + m) * s"*)

	       Rule.Thm ("sym_realpow_addI",

               ThmC.numerals_to_Free (@{thm realpow_addI} RS @{thm sym})),

	       (*"r \<up> n * r \<up> m = r \<up> (n + m)"*)

	       Rule.Thm ("realpow_addI_assoc_l",ThmC.numerals_to_Free @{thm realpow_addI_assoc_l}),

	       (*"r \<up> n * (r \<up> m * s) = r \<up> (n + m) * s"*)

	       (* ist in expand_poly - wird hier aber auch gebraucht, wegen: 

		  "r * r = r \<up> 2" wenn r=a \<up> b*)

	       Rule.Thm ("realpow_pow",ThmC.numerals_to_Free @{thm realpow_pow})

	       (*"(a \<up> b) \<up> c = a \<up> (b * c)"*)

	       ], scr = Rule.Empty_Prog};

val calc_add_mult_pow_ = 

  Rule_Def.Repeat{id = "calc_add_mult_pow_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Atools_erls(*erls3.4.03*),srls = Rule_Set.Empty,

      calc = [("PLUS"  , ("Groups.plus_class.plus", eval_binop "#add_")), 

	      ("TIMES" , ("Groups.times_class.times", eval_binop "#mult_")),

	      ("POWER", ("Transcendental.powr", eval_binop "#power_"))

	      ],

      errpatts = [],

      rules = [Rule.Eval ("Groups.plus_class.plus", eval_binop "#add_"),

	       Rule.Eval ("Groups.times_class.times", eval_binop "#mult_"),

	       Rule.Eval ("Transcendental.powr", eval_binop "#power_")

	       ], scr = Rule.Empty_Prog};

val reduce_012_mult_ = 

  Rule_Def.Repeat{id = "reduce_012_mult_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty,srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      rules = [(* MG: folgende Rule.Thm müssen hier stehen bleiben: *)

               Rule.Thm ("mult_1_right",ThmC.numerals_to_Free @{thm mult_1_right}),

	       (*"z * 1 = z"*) (*wegen "a * b * b \<up> (-1) + a"*) 

	       Rule.Thm ("realpow_zeroI",ThmC.numerals_to_Free @{thm realpow_zeroI}),

	       (*"r \<up> 0 = 1"*) (*wegen "a*a \<up> (-1)*c + b + c"*)

	       Rule.Thm ("realpow_oneI",ThmC.numerals_to_Free @{thm realpow_oneI}),

	       (*"r \<up> 1 = r"*)

	       Rule.Thm ("realpow_eq_oneI",ThmC.numerals_to_Free @{thm realpow_eq_oneI})

	       (*"1 \<up> n = 1"*)

	       ], scr = Rule.Empty_Prog};

val collect_numerals_ = 

  Rule_Def.Repeat{id = "collect_numerals_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Atools_erls, srls = Rule_Set.Empty,

      calc = [("PLUS"  , ("Groups.plus_class.plus", eval_binop "#add_"))

	      ], errpatts = [],

      rules = 

        [Rule.Thm ("real_num_collect",ThmC.numerals_to_Free @{thm real_num_collect}), 

	     (*"[| l is_const; m is_const |]==>l * n + m * n = (l + m) * n"*)

	 Rule.Thm ("real_num_collect_assoc_r",ThmC.numerals_to_Free @{thm real_num_collect_assoc_r}),

	     (*"[| l is_const; m is_const |] ==>  \

					\(k + m * n) + l * n = k + (l + m)*n"*)

	 Rule.Thm ("real_one_collect",ThmC.numerals_to_Free @{thm real_one_collect}),	

	     (*"m is_const ==> n + m * n = (1 + m) * n"*)

	 Rule.Thm ("real_one_collect_assoc_r",ThmC.numerals_to_Free @{thm real_one_collect_assoc_r}), 

	     (*"m is_const ==> (k + n) + m * n = k + (m + 1) * n"*)

         Rule.Eval ("Groups.plus_class.plus", eval_binop "#add_"),

	 (*MG: Reihenfolge der folgenden 2 Rule.Thm muss so bleiben, wegen

		     (a+a)+a --> a + 2*a --> 3*a and not (a+a)+a --> 2*a + a *)

         Rule.Thm ("real_mult_2_assoc_r",ThmC.numerals_to_Free @{thm real_mult_2_assoc_r}),

	     (*"(k + z1) + z1 = k + 2 * z1"*)

	 Rule.Thm ("sym_real_mult_2",ThmC.numerals_to_Free (@{thm real_mult_2} RS @{thm sym}))

	     (*"z1 + z1 = 2 * z1"*)

	], scr = Rule.Empty_Prog};

val reduce_012_ = 

  Rule_Def.Repeat{id = "reduce_012_", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty,srls = Rule_Set.Empty, calc = [], errpatts = [],

      rules = [Rule.Thm ("mult_1_left",ThmC.numerals_to_Free @{thm mult_1_left}),                 

	       (*"1 * z = z"*)

	       Rule.Thm ("mult_zero_left",ThmC.numerals_to_Free @{thm mult_zero_left}),        

	       (*"0 * z = 0"*)

	       Rule.Thm ("mult_zero_right",ThmC.numerals_to_Free @{thm mult_zero_right}),

	       (*"z * 0 = 0"*)

	       Rule.Thm ("add_0_left",ThmC.numerals_to_Free @{thm add_0_left}),

	       (*"0 + z = z"*)

	       Rule.Thm ("add_0_right",ThmC.numerals_to_Free @{thm add_0_right}),

	       (*"z + 0 = z"*) (*wegen a+b-b --> a+(1-1)*b --> a+0 --> a*)

	       (*Rule.Thm ("realpow_oneI",ThmC.numerals_to_Free @{thm realpow_oneI})*)

	       (*"?r \<up> 1 = ?r"*)

	       Rule.Thm ("division_ring_divide_zero",ThmC.numerals_to_Free @{thm division_ring_divide_zero})

	       (*"0 / ?x = 0"*)

	       ], scr = Rule.Empty_Prog};

val discard_parentheses1 = 

    Rule_Set.append_rules "discard_parentheses1" Rule_Set.empty 

	       [Rule.Thm ("sym_mult.assoc",

                      ThmC.numerals_to_Free (@{thm mult.assoc} RS @{thm sym}))

		(*"?z1.1 * (?z2.1 * ?z3.1) = ?z1.1 * ?z2.1 * ?z3.1"*)

		(*Rule.Thm ("sym_add.assoc",

                        ThmC.numerals_to_Free (@{thm add.assoc} RS @{thm sym}))*)

		(*"?z1.1 + (?z2.1 + ?z3.1) = ?z1.1 + ?z2.1 + ?z3.1"*)

		 ];

val expand_poly =

  Rule_Def.Repeat{id = "expand_poly", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty,srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      (*asm_thm = [],*)

      rules = [Rule.Thm ("distrib_right" ,ThmC.numerals_to_Free @{thm distrib_right}),

	       (*"(z1.0 + z2.0) * w = z1.0 * w + z2.0 * w"*)

	       Rule.Thm ("distrib_left",ThmC.numerals_to_Free @{thm distrib_left}),

	       (*"w * (z1.0 + z2.0) = w * z1.0 + w * z2.0"*)

	       (*Rule.Thm ("distrib_right1",ThmC.numerals_to_Free @{thm distrib_right}1),

		....... 18.3.03 undefined???*)

	       Rule.Thm ("real_plus_binom_pow2",ThmC.numerals_to_Free @{thm real_plus_binom_pow2}),

	       (*"(a + b) \<up> 2 = a \<up> 2 + 2*a*b + b \<up> 2"*)

	       Rule.Thm ("real_minus_binom_pow2_p",ThmC.numerals_to_Free @{thm real_minus_binom_pow2_p}),

	       (*"(a - b) \<up> 2 = a \<up> 2 + -2*a*b + b \<up> 2"*)

	       Rule.Thm ("real_plus_minus_binom1_p",

		    ThmC.numerals_to_Free @{thm real_plus_minus_binom1_p}),

	       (*"(a + b)*(a - b) = a \<up> 2 + -1*b \<up> 2"*)

	       Rule.Thm ("real_plus_minus_binom2_p",

		    ThmC.numerals_to_Free @{thm real_plus_minus_binom2_p}),

	       (*"(a - b)*(a + b) = a \<up> 2 + -1*b \<up> 2"*)

	       Rule.Thm ("minus_minus",ThmC.numerals_to_Free @{thm minus_minus}),

	       (*"- (- ?z) = ?z"*)

	       Rule.Thm ("real_diff_minus",ThmC.numerals_to_Free @{thm real_diff_minus}),

	       (*"a - b = a + -1 * b"*)

	       Rule.Thm ("sym_real_mult_minus1",

                     ThmC.numerals_to_Free (@{thm real_mult_minus1} RS @{thm sym}))

	       (*- ?z = "-1 * ?z"*)

	       (*Rule.Thm ("real_minus_add_distrib",

		      ThmC.numerals_to_Free @{thm real_minus_add_distrib}),*)

	       (*"- (?x + ?y) = - ?x + - ?y"*)

	       (*Rule.Thm ("real_diff_plus",ThmC.numerals_to_Free @{thm real_diff_plus})*)

	       (*"a - b = a + -b"*)

	       ], scr = Rule.Empty_Prog};

val simplify_power = 

  Rule_Def.Repeat{id = "simplify_power", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty, srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      rules = [Rule.Thm ("realpow_multI", ThmC.numerals_to_Free @{thm realpow_multI}),

	       (*"(r * s) \<up> n = r \<up> n * s \<up> n"*)

	       Rule.Thm ("sym_realpow_twoI",

                     ThmC.numerals_to_Free( @{thm realpow_twoI} RS @{thm sym})),	

	       (*"r1 * r1 = r1 \<up> 2"*)

	       Rule.Thm ("realpow_plus_1",ThmC.numerals_to_Free @{thm realpow_plus_1}),		

	       (*"r * r \<up> n = r \<up> (n + 1)"*)

	       Rule.Thm ("realpow_pow",ThmC.numerals_to_Free @{thm realpow_pow}),

	       (*"(a \<up> b) \<up> c = a \<up> (b * c)"*)

	       Rule.Thm ("sym_realpow_addI",

                     ThmC.numerals_to_Free (@{thm realpow_addI} RS @{thm sym})),

	       (*"r \<up> n * r \<up> m = r \<up> (n + m)"*)

	       Rule.Thm ("realpow_oneI",ThmC.numerals_to_Free @{thm realpow_oneI}),

	       (*"r \<up> 1 = r"*)

	       Rule.Thm ("realpow_eq_oneI",ThmC.numerals_to_Free @{thm realpow_eq_oneI})

	       (*"1 \<up> n = 1"*)

	       ], scr = Rule.Empty_Prog};

val collect_numerals = 

  Rule_Def.Repeat{id = "collect_numerals", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Atools_erls(*erls3.4.03*),srls = Rule_Set.Empty,

      calc = [("PLUS"  , ("Groups.plus_class.plus", eval_binop "#add_")), 

	      ("TIMES" , ("Groups.times_class.times", eval_binop "#mult_")),

	      ("POWER", ("Transcendental.powr", eval_binop "#power_"))

	      ], errpatts = [],

      rules = [Rule.Thm ("real_num_collect",ThmC.numerals_to_Free @{thm real_num_collect}), 

	       (*"[| l is_const; m is_const |]==>l * n + m * n = (l + m) * n"*)

	       Rule.Thm ("real_num_collect_assoc",ThmC.numerals_to_Free @{thm real_num_collect_assoc}),

	       (*"[| l is_const; m is_const |] ==>  

				l * n + (m * n + k) =  (l + m) * n + k"*)

	       Rule.Thm ("real_one_collect",ThmC.numerals_to_Free @{thm real_one_collect}),	

	       (*"m is_const ==> n + m * n = (1 + m) * n"*)

	       Rule.Thm ("real_one_collect_assoc",ThmC.numerals_to_Free @{thm real_one_collect_assoc}), 

	       (*"m is_const ==> k + (n + m * n) = k + (1 + m) * n"*)

	       Rule.Eval ("Groups.plus_class.plus", eval_binop "#add_"), 

	       Rule.Eval ("Groups.times_class.times", eval_binop "#mult_"),

	       Rule.Eval ("Transcendental.powr", eval_binop "#power_")

	       ], scr = Rule.Empty_Prog};

val reduce_012 = 

  Rule_Def.Repeat{id = "reduce_012", preconds = [], 

      rew_ord = ("dummy_ord", Rewrite_Ord.dummy_ord),

      erls = Rule_Set.empty,srls = Rule_Set.Empty,

      calc = [], errpatts = [],

      rules = [Rule.Thm ("mult_1_left",ThmC.numerals_to_Free @{thm mult_1_left}),                 

	       (*"1 * z = z"*)

	       (*Rule.Thm ("real_mult_minus1",ThmC.numerals_to_Free @{thm real_mult_minus1}),14.3.03*)

	       (*"-1 * z = - z"*)

	       Rule.Thm ("minus_mult_left", 

		    ThmC.numerals_to_Free (@{thm minus_mult_left} RS @{thm sym})),

	       (*- (?x * ?y) = "- ?x * ?y"*)

	       (*Rule.Thm ("real_minus_mult_cancel",

                       ThmC.numerals_to_Free @{thm real_minus_mult_cancel}),

	       (*"- ?x * - ?y = ?x * ?y"*)---*)

	       Rule.Thm ("mult_zero_left",ThmC.numerals_to_Free @{thm mult_zero_left}),        

	       (*"0 * z = 0"*)

	       Rule.Thm ("add_0_left",ThmC.numerals_to_Free @{thm add_0_left}),

	       (*"0 + z = z"*)

	       Rule.Thm ("right_minus",ThmC.numerals_to_Free @{thm right_minus}),

	       (*"?z + - ?z = 0"*)

	       Rule.Thm ("sym_real_mult_2",

                     ThmC.numerals_to_Free (@{thm real_mult_2} RS @{thm sym})),