(* Title: test for rewtools.sml

   authors: Walther Neuper 2000, 2006

*)

"--------------------------------------------------------";

"--------------------------------------------------------";

"table of contents --------------------------------------";

"--------------------------------------------------------";

"----------- fun make_isac ------------------------------";

"----------- fun collect_isab_thys ----------------------";

"----------- fun thy_containing_thm ---------------------";

"----------- fun thy_containing_rls ---------------------";

"----------- fun thy_containing_cal ---------------------";

"----------- initContext Thy_ Integration-demo ----------";

"----------- initContext..Thy_, fun context_thm ---------";

"----------- initContext..Thy_, fun context_rls ---------";

"----------- checkContext..Thy_, fun context_thy --------";

"----------- checkContext..Thy_, fun context_rls --------";

"----------- checkContext..Thy_ on last formula ---------"; 

"----------- fun guh2theID ------------------------------";

"----------- debugging on Tests/solve_linear_as_rootpbl -";

"--------------------------------------------------------";

"----------- fun string_of_thmI for_[.]_) ---------------";

"----------- (sym_real_minus_eq_cancel, (?b1 = ?a1) ..._[";

"--------------------------------------------------------";

"----------- fun filter_appl_rews -----------------------";

"----------- fun is_contained_in ------------------------";

"--------------------------------------------------------";

"--------------------------------------------------------";

"----------- fun make_isac ------------------------------";

"----------- fun make_isac ------------------------------";

"----------- fun make_isac ------------------------------";

"----- rlsthmsNOTisac: in rls requested from Isabelle ---";

map #1 (rlsthmsNOTisac : (string * thm) list) = (*WN101011, Isabelle2009-2*)

["refl", "o_apply", "del_base", "del_rec", "LENGTH_CONS", "LENGTH_NIL",

 "list_diff_def", "take_Nil", "take_Cons", "if_True", "if_False",

 "sym_real_mult_minus1", "left_distrib", "right_distrib",

 "sym_realpow_twoI", "sym_realpow_addI", "mult_1_right",

 "sym_real_mult_2", "mult_1_left", "mult_zero_left", "mult_zero_right",

 "add_0_left", "add_0_right", "divide_zero_left", "sym_real_mult_assoc",

 "real_le_refl", "minus_minus", "real_mult_commute", "real_mult_assoc",

 "add_commute", "add_left_commute", "add_assoc", "minus_mult_left",

 "right_minus", "sym_add_assoc", "left_diff_distrib",

 "right_diff_distrib", "minus_divide_left", "times_divide_eq_right",

 "times_divide_eq_left", "divide_divide_eq_left",

 "divide_divide_eq_right", "divide_1", "add_divide_distrib",

 "sym_rmult_assoc"]; 

if length rlsthmsNOTisac > 40 then () 

else error "rewtools.sml: some rlsthmsNOTisac dropped";

"----- FIXME: rlsthmsNOTisac DOES contain thms from isac ---";

map (Thm.derivation_name o #2) rlsthmsNOTisac;

"----- FIXME: sym_* in rlsthmsNOTisac DOES contain thms from isac ---";

val symthms = [("real_mult_minus1", @{thm real_mult_minus1}),

               ("realpow_twoI", @{thm realpow_twoI}),

               ("realpow_addI", @{thm realpow_addI}),

               ("real_mult_2", @{thm real_mult_2}),

               ("real_mult_assoc", @{thm real_mult_assoc}),

               ("add_assoc", @{thm add_assoc}),

               ("rmult_assoc", @{thm rmult_assoc})];

map (Thm.derivation_name o #2) symthms;

"----------- fun collect_isab_thys ----------------------";

"----------- fun collect_isab_thys ----------------------";

"----------- fun collect_isab_thys ----------------------";

val ancestors = Theory.ancestors_of @{theory "Complex_Main"};

val isabthms = (flat o (map Theory.axioms_of)) ancestors;

val isacrules = (flat o (map (thms_of_rls o #2 o #2))) (! ruleset');

(*thms from rulesets*)

(*=== inhibit exn AK110725 =============================================================

val isacrlsthms = ((map rep_thm_G') o flat o 

		 (map (PureThy.all_thms_of_rls o #2 o #2))) (!ruleset');

length isacrlsthms;

(*takes a few seconds...

val isacrlsthms = gen_distinct eq_thmI isacrlsthms;

length isacrlsthms;

"----- theorems used in isac's rulesets -----vvvvvvvvvvvvvvvvvvvvv";

print_depth 99; map #1 isacrlsthms; print_depth 3;

"----- theorems used in isac's rulesets -----^^^^^^^^^^^^^^^^^^^^^";

...*)

(!theory');

map #2 (!theory');

map (PureThy.all_thms_of o #2) (!theory');

val isacthms = (flat o (map (PureThy.all_thms_of o #2))) (!theory');

(*takes a few seconds...

val rlsthmsNOTisac = gen_diff eq_thmI (isacrlsthms, isacthms);

length rlsthmsNOTisac;

"----- theorems from rulesets NOT def. in isac -----vvvvvvvvvvvvvv";

print_depth 99; map #1 rlsthmsNOTisac; print_depth 3;

"----- theorems from rulesets NOT def. in isac -----^^^^^^^^^^^^^^";

...*)

"----- for 'fun make_isab_thm_thy'";

inter eq_thmI isacrlsthms (PureThy.all_thms_of (nth 1 ancestors));

inter eq_thmI;

(inter eq_thmI);

(inter eq_thmI) isacrlsthms;

(*takes a few seconds...

curry (inter eq_thmI) isacrlsthms (PureThy.all_thms_of (nth 9 ancestors));

val thy = (nth 52 ancestors);

val sec = ((inter eq_thmI) isacrlsthms o PureThy.all_thms_of) (nth 52 ancestors);

length (PureThy.all_thms_of (nth 9 ancestors));

length sec;

...*)

(*takes 1 minute...

print_depth 99; 

map ((inter eq_thmI) rlsthmsNOTisac o PureThy.all_thms_of) ancestors;

print_depth 3;

*)

(*takes some seconds...

val isab_thm_thy = (flat o (map (make_isab_thm_thy rlsthmsNOTisac)))

		       ((#ancestors o rep_theory) first_isac_thy);

print_depth 99; isab_thm_thy; print_depth 3;

*)

=== inhibit exn AK110725 =============================================================*)

"----------- fun thy_containing_thm ---------------------";

"----------- fun thy_containing_thm ---------------------";

"----------- fun thy_containing_thm ---------------------";

val (str, (thy', thy)) = ("real_diff_minus",("Root", @{theory Root}));

if thy_contains_thm str ("XXX",thy) then ()

  else error "rewtools.sml: NOT thy_contains_thm (real_diff_minus,(Root.thy,.";

(rev (!theory'));

dropuntil (curry op = thy');

dropuntil ((curry op = thy') o (#1:theory' * theory -> theory'));

val startsearch = dropuntil ((curry op = thy') o 

			     (#1:theory' * theory -> theory')) 

			    (rev (!theory'));

if thy_containing_thm thy' str = ("IsacKnowledge", "Root") then ()

  else error "rewtools.sml: NOT thy_containin_thm (real_diff_minus,(Root.thy,.";

"----- search the same theorem somewhere further below in the list";

if thy_contains_thm str ("XXX",@{theory Poly}) then ()

  else error "rewtools.sml: NOT thy_contains_thm (real_diff_minus,(Poly.thy,.";

(*=== inhibit exn AK110725 =============================================================

(* AK110725 "LineEq" can't be "Poly" - change "LineEq" to "Poly"??? *)

if thy_containing_thm "LinEq" str = ("IsacKnowledge", "Poly") then ()

  else error "rewtools.sml: NOT thy_containing_thm (real_diff_minus,(Poly.thy,.";

=== inhibit exn AK110725 =============================================================*)

"----- second test -------------------------------";

show_thes();

(*args vor thy_containing_thm...*)

val (thy',str) = ("Test", "radd_commute");

val startsearch = dropuntil ((curry op= thy') o 

				     (#1:theory' * theory -> theory')) 

				    (rev (!theory'));

length (!theory');

length startsearch;

if thy_containing_thm thy' str = ("IsacKnowledge", "Test") then ()

else error "rewtools.sml: diff.behav. in \

		 \thy_containing_thm Test radd_commute";

"----------- fun thy_containing_rls ---------------------";

"----------- fun thy_containing_rls ---------------------";

"----------- fun thy_containing_rls ---------------------";

val thy' = "Biegelinie";

val dropthys = takewhile [] (not o (curry op= thy') o 

			     (#1:theory' * theory -> theory')) 

			 (rev (!theory'));

if length (!theory') <> length dropthys then ()

else error "rewtools.sml: diff.behav. in thy_containing_rls 1";

val dropthy's = map (get_thy o (#1 : (theory' * theory) -> theory'))

		    dropthys;

print_depth 99; dropthy's; print_depth 3;

(*WN100819=======================================================

WN100819 THIS DATE MUST BE WRONG: at this date changeset 37934:56f10b13005e

finished update ME/calchead.sml + pushed updates over all sml+test

"Isac" mem dropthy's;

op mem ("Isac", dropthy's);

(op mem) o swap;

((op mem) o swap) (dropthy's, "Isac");

curry ((op mem) o swap);

curry ((op mem) o swap) dropthy's "Isac";

val startsearch = filter_out ((curry ((op mem) o swap) dropthy's) o	 

			      ((#1 o #2) : rls' * (theory' * rls) -> theory'))

			     (rev (!ruleset'));

print_depth 99; map (#1 o #2) startsearch; print_depth 3;

if length (!ruleset') <> length startsearch then ()

else error "rewtools.sml: diff.behav. in thy_containing_rls 2";

val rls' = "norm_Poly";

case assoc (startsearch, rls') of

    SOME (thy', _) => thyID2theory' thy'

  | _ => error ("thy_containing_rls : rls '"^str^

			  "' not in !rulset' und thy '"^thy'^"'");

if thy_containing_rls thy' rls' = ("IsacKnowledge", "Poly") then ()

else error "rewtools.sml: diff.behav. in thy_containing_rls 3";

"----------- fun thy_containing_cal ---------------------";

"----------- fun thy_containing_cal ---------------------";

"----------- fun thy_containing_cal ---------------------";

val thy' = "Atools";

val dropthys = takewhile [] (not o (curry op= thy') o 

			     (#1:theory' * theory -> theory')) 

			 (rev (!theory'));

length dropthys <> length (!theory');

val dropthy's = map (get_thy o (#1 : (theory' * theory) -> theory'))

		    dropthys;

(rev (!calclist'));

map #1 (rev (!calclist'));

map (#1 : calc -> string) (rev (!calclist'));

val startsearch = filter_out ((curry ((op mem) o swap) dropthy's) o

			      (#1 : calc -> string)) (rev (!calclist'));

WN100819 =====================================================*)

"----------- initContext Thy_ Integration-demo ----------";

"----------- initContext Thy_ Integration-demo ----------";

"----------- initContext Thy_ Integration-demo ----------";

states:=[];

CalcTree

[(["functionTerm (2 * x)","integrateBy x","antiDerivative FF"], 

  ("Integrate",["integrate","function"],

  ["diff","integration"]))];

Iterator 1;

moveActiveRoot 1;

autoCalculate 1 CompleteCalc;

interSteps 1 ([1],Res);

interSteps 1 ([1,1],Res);

val ((pt,p),_) = get_calc 1; show_pt pt;

if existpt' ([1,1,1], Frm) pt then ()

else error "integrate.sml: interSteps on Rewrite_Set_Inst 1";

initContext  1 Thy_ ([1,1,1], Frm);

"----------- initContext..Thy_, fun context_thm ---------";

"----------- initContext..Thy_, fun context_thm ---------";

"----------- initContext..Thy_, fun context_thm ---------";

states:=[]; (*start of calculation, return No.1*)

CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 

  ("Test", ["sqroot-test","univariate","equation","test"],

   ["Test","squ-equ-test-subpbl1"]))];

Iterator 1; moveActiveRoot 1;

autoCalculate 1 CompleteCalc;

"----- no thy-context at result -----";

val p = ([], Res);

initContext 1 Thy_ p;

interSteps 1 ([2], Res);

interSteps 1 ([3,1], Res);

val ((pt,_),_) = get_calc 1; show_pt pt;

val p = ([2,4], Res);

val tac = Rewrite ("radd_left_commute","");

initContext 1 Thy_ p;

(*Res->Res, Rewrite "radd_left_commute 1 + (-2 + x) = 0 -> -2 + (1 + x) = 0

  --- in initContext..Thy_ ---*)

val ContThm {thm,result,...} = context_thy (pt,p) tac;

if thm = "thy_isac_Test-thm-radd_left_commute" 

   andalso term2str result = "-2 + (1 + x) = 0" then ()

else error"rewtools.sml initContext..Th_ thy_Test-thm-radd_left_commute";

val p = ([3,1,1], Frm);

val tac = Rewrite_Inst (["(bdv, x)"],("risolate_bdv_add",""));

initContext 1 Thy_ p;

(*Frm->Res, Rewrite_Inst "risolate_bdv_add"  -1 + x = 0 -> x = 0 + -1 * -1

  --- in initContext..Thy_ ---*)

val ContThmInst {thm,result,...} = context_thy (pt,p) tac;

if thm =  "thy_isac_Test-thm-risolate_bdv_add"

   andalso term2str result = "x = 0 + -1 * -1" then ()

else error "rewtools.sml initContext..Th_ thy_Test-thm-risolate_bdv_add";

initContext 1 Thy_ ([2,5], Res);

(*Res->Res, Calculate "plus"  -2 + (1 + x) = 0 -> -1 + x = 0

  --- in initContext..Thy_ ---*)

"----------- initContext..Thy_, fun context_rls ---------";

"----------- initContext..Thy_, fun context_rls ---------";

"----------- initContext..Thy_, fun context_rls ---------";

(*using pt from above*)

val p = ([1], Res);

val tac = Rewrite_Set "Test_simplify";

initContext 1 Thy_ p;

(*Res->Res, Rewrite_Set "Test_simplify" x + 1 + -1 * 2 = 0 -> -1 + x = 0

  --- in initContext..Thy_ ---*)

val ContRls {rls,result,...} = context_thy (pt,p) tac;

if rls = "thy_isac_Test-rls-Test_simplify" 

   andalso term2str result = "-1 + x = 0" then ()

else error "rewtools.sml initContext..Th_ thy_Test-thm-risolate_bdv_add";

val p = ([3,1], Frm);

val tac = Rewrite_Set_Inst (["(bdv, x)"],"isolate_bdv");

initContext 1 Thy_ p;

(*Frm->Res, Rewrite_Set_Inst "isolate_bdv" -1 + x = 0 ->  x = 0 + -1 * -1

  --- in initContext..Thy_ ---*)

val ContRlsInst {rls,result,...} = context_thy (pt,p) tac;

if rls =  "thy_isac_Test-rls-isolate_bdv"

   andalso term2str result = "x = 0 + -1 * -1" then ()

else error "rewtools.sml initContext..Th_ thy_Test-thm-risolate_bdv_add";

"----------- checkContext..Thy_, fun context_thy --------";

"----------- checkContext..Thy_, fun context_thy --------";

"----------- checkContext..Thy_, fun context_thy --------";

(*using pt from above*)

val p = ([2,4], Res);

val tac = Rewrite ("radd_left_commute","");

checkContext 1 p "thy_Test-thm-radd_left_commute";

(* radd_left_commute: 1 + (-2 + x) = 0 -> -2 + (1 + x) = 0

  --- in checkContext..Thy_ ---*)

val ContThm {thm,result,...} = context_thy (pt,p) tac;

if thm =  "thy_isac_Test-thm-radd_left_commute"

   andalso term2str result = "-2 + (1 + x) = 0" then ()

else error "rewtools.sml checkContext.._ thy_Test-thm-radd_left_commute";

val p = ([3,1,1], Frm);

val tac = Rewrite_Inst (["(bdv,x)"],("risolate_bdv_add",""));

checkContext 1 p "thy_Test-thm-risolate_bdv_add";

(* risolate_bdv_add:  -1 + x = 0 -> x = 0 + -1 * -1

  --- in checkContext..Thy_ ---*)

val ContThmInst {thm,result,...} = context_thy (pt,p) tac;

if thm =  "thy_isac_Test-thm-risolate_bdv_add"

   andalso term2str result = "x = 0 + -1 * -1" then ()

else error "rewtools.sml checkContext..T_ thy_Test-thm-risolate_bdv_add";

val p = ([2,5], Res);

val tac = Calculate "plus";

(*checkContext..Thy_ 1 ([2,5], Res);*)

(*FIXXXME #######################vvvvv kestoreID !!!!!!!!!!!!!!!!*)

checkContext 1 p ;

(* Calculate "plus"  -2 + (1 + x) = 0 -> -1 + x = 0

  --- in checkContext..Thy_ ---*)

"----------- checkContext..Thy_, fun context_rls --------";

"----------- checkContext..Thy_, fun context_rls --------";

"----------- checkContext..Thy_, fun context_rls --------";

(*using pt from above*)

show_pt pt;

val p = ([1], Res);

val tac = Rewrite_Set "Test_simplify";

checkContext 1 p "thy_isac_Test-rls-Test_simplify";

(*Res->Res, Rewrite_Set "Test_simplify" x + 1 + -1 * 2 = 0 -> -1 + x = 0

  --- in checkContext..Thy_ ---*)

val ContRls {rls,result,...} = context_thy (pt,p) tac;

if rls = "thy_isac_Test-rls-Test_simplify" 

   andalso term2str result = "-1 + x = 0" then ()

else error "rewtools.sml checkContext..Thy_ thy_Test-rls-Test_simplify";

val p = ([3,1], Frm);

val tac = Rewrite_Set_Inst (["(bdv, x)"],"isolate_bdv");

checkContext 1 p "thy_Test-rls-isolate_bdv";

val ContRlsInst {rls,result,...} = context_thy (pt,p) tac;

if rls = "thy_isac_Test-rls-isolate_bdv" 

   andalso term2str result = "x = 0 + -1 * -1" then ()

else error "rewtools.sml checkContext..Thy_ thy_Test-thm-isolate_bdv";

"----------- checkContext..Thy_ on last formula ---------"; 

"----------- checkContext..Thy_ on last formula ---------"; 

"----------- checkContext..Thy_ on last formula ---------"; 

states:=[]; (*start of calculation, return No.1*)

CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 

  ("Test", ["sqroot-test","univariate","equation","test"],

   ["Test","squ-equ-test-subpbl1"]))];

Iterator 1; moveActiveRoot 1;

autoCalculate 1 CompleteCalcHead;

autoCalculate 1 (Step 1);

val (ptp as (pt,p), tacis) = get_calc 1; show_pt pt;

initContext 1 Thy_ ([1], Frm);

checkContext 1 ([1], Frm) "thy_isac_Test-thm-radd_left_commute";

autoCalculate 1 (Step 1);

val (ptp as (pt,p), tacis) = get_calc 1; show_pt pt;

initContext 1 Thy_ ([1], Res);

checkContext 1 ([1], Res) "thy_isac_Test-rls-Test_simplify";

"----------- fun guh2theID ------------------------------";

"----------- fun guh2theID ------------------------------";

"----------- fun guh2theID ------------------------------";

val guh = "thy_scri_ListG-thm-zip_Nil";

print_depth 999;

take_fromto 1 4 (Symbol.explode guh);

take_fromto 5 9 (Symbol.explode guh);

val rest = takerest (9,(Symbol.explode guh)); 

separate "-" rest;

space_implode "-" rest;

commas rest;

val delim = "-";

val thyID = takewhile [] (not o (curry op= delim)) rest;

val rest' = dropuntil (curry op= delim) rest;

val setc = take_fromto 1 5 rest';

val xstr = takerest (5, rest');

if guh2theID guh = ["IsacScripts", "ListG", "Theorems", "zip_Nil"] then ()

else error "rewtools.sml: guh2theID thy_scri_ListG-thm-zip_Nil changed";

"----------- debugging on Tests/solve_linear_as_rootpbl -";

"----------- debugging on Tests/solve_linear_as_rootpbl -";

"----------- debugging on Tests/solve_linear_as_rootpbl -";

"----- initContext -----";

states:=[];

CalcTree 

    [(["equality (1+-1*2+x=(0::real))", "solveFor x", "solutions L"],

      ("Test",

       ["linear","univariate","equation","test"],

       ["Test","solve_linear"]))];

Iterator 1; moveActiveRoot 1;

autoCalculate 1 CompleteCalcHead;

autoCalculate 1 (Step 1); val (ptp as (pt,p), tacis) = get_calc 1; show_pt pt;

if is_curr_endof_calc pt ([1],Frm) then ()

else error "rewtools.sml is_curr_endof_calc ([1],Frm)";

autoCalculate 1 (Step 1); val (ptp as (pt,p), tacis) = get_calc 1; show_pt pt;

if not (is_curr_endof_calc pt ([1],Frm)) then ()

else error "rewtools.sml is_curr_endof_calc ([1],Frm) not";

if is_curr_endof_calc pt ([1],Res) then ()

else error "rewtools.sml is_curr_endof_calc ([1],Res)";

initContext 1 Thy_ ([1],Res);

"----- checkContext -----";

states:=[];

CalcTree 

    [(["equality (1+-1*2+x=(0::real))", "solveFor x", "solutions L"],

      ("Test",

       ["linear","univariate","equation","test"],

       ["Test","solve_linear"]))];

Iterator 1; moveActiveRoot 1;

autoCalculate 1 CompleteCalc;

interSteps 1 ([1],Res);

val (ptp as (pt,p), tacis) = get_calc 1; show_pt pt;

checkContext 1 ([1],Res) "thy_isac_Test-rls-Test_simplify";

interSteps 1 ([2],Res);

val (ptp as (pt,p), tacis) = get_calc 1; show_pt pt;

checkContext 1 ([2,1],Res) "thy_isac_Test-rls-Test_simplify";

checkContext 1 ([2,2],Res) "thy_isac_Test-rls-Test_simplify";

"----------- fun string_of_thmI for_[.]_) ---------------";

"----------- fun string_of_thmI for_[.]_) ---------------";

"----------- fun string_of_thmI for_[.]_) ---------------";

"----- these work ?!?";

(*========== inhibit exn 110718 ================================================

(* ERROR: constructor real_minus_eq_cancel has not been declared *)

val th = sym_thm real_minus_eq_cancel;

val Th = sym_Thm (Thm ("real_minus_eq_cancel", real_minus_eq_cancel));

val th'= mk_thm ( @{theory "Isac"} ) ((de_quote o string_of_thm) real_minus_eq_cancel);

val th'= mk_thm Biegelinie.thy((de_quote o string_of_thm) real_minus_eq_cancel);

"----- DIFFERENCE TO ABOVE ?!?: this is still ok, when called in next_tac...";

val ((pt,(p,p_)), _) = get_calc 1; show_pt pt;

val Appl (Rewrite' (_,_,_,_,thm',_,_)) =

    applicable_in (p,p_) pt (Rewrite ("sym_real_minus_eq_cancel",""));

"- compose stac as done in | appy (*leaf*) by handle_leaf";

val (th, sr, E, a, v, t) = 

    ("Biegelinie", 

     (#srls o get_met) ["IntegrierenUndKonstanteBestimmen"],

     [(str2term "q__::bool", str2term "q x = q_0")], 

     SOME (str2term "q x = q_0"),

     str2term "q__::bool", 

     str2term "(Rewrite sym_real_minus_eq_cancel False) (q__::bool)");

val (a', STac stac) = handle_leaf "next  " th sr E a v t;

term2str stac;

"--- but this \"m\" is already corrupted";

val (m,_) = stac2tac_ EmptyPtree (assoc_thy th) stac;

"- because in assoc_thm'...";

val (thy, (thmid, ct')) = (Biegelinie.thy, ("sym_real_minus_eq_cancel",""));

val "s"::"y"::"m"::"_"::id = explode thmid;

((num_str o (get_thm thy)) (implode id)) RS sym;

((get_thm thy) (implode id)) RS sym;

"... this creates [.]";

((get_thm thy) (implode id));

"... while this has _NO_ [.]";

"----- thus we repair the [.] in string_of_thmI...";

val thm = ((num_str o (get_thm thy)) (implode id)) RS sym;

if string_of_thmI thm = "(?b1 = ?a1) = (- ?b1 = - ?a1)" then ()

else error ("rewtools.sml: string_of_thmI " ^ string_of_thm thm ^

		  " = " ^ string_of_thmI thm);

========== inhibit exn 110718 ================================================*)

"----------- (sym_real_minus_eq_cancel, (?b1 = ?a1) ..._[";

"----------- (sym_real_minus_eq_cancel, (?b1 = ?a1) ..._[";

"----------- (sym_real_minus_eq_cancel, (?b1 = ?a1) ..._[";

states:=[];

CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",

	     "RandbedingungenBiegung [y 0 = 0, y L = 0]",

	     "RandbedingungenMoment [M_b 0 = 0, M_b L = 0]",

	     "FunktionsVariable x"],

	    ("Biegelinie",

	     ["MomentBestimmte","Biegelinien"],

	     ["IntegrierenUndKonstanteBestimmen"]))];

moveActiveRoot 1;

autoCalculate 1 CompleteCalcHead;

autoCalculate 1 (Step 1) (*Apply_Method*);

autoCalculate 1 (Step 1) (*->GENERATED ([1], Frm)*);

(*========== inhibit exn 110718 ================================================

"--- this was corrupted before 'fun string_of_thmI'";

val ((pt,(p,p_)), _) = get_calc 1; show_pt pt;

if get_obj g_tac pt p = Rewrite ("sym_real_minus_eq_cancel", 

				 "(?b1 = ?a1) = (- ?b1 = - ?a1)") then ()

else error "rewtools.sml: string_of_thmI ?!?";

getTactic 1 ([1],Frm);

"----------- fun filter_appl_rews -----------------------";

"----------- fun filter_appl_rews -----------------------";

"----------- fun filter_appl_rews -----------------------";

val f = str2term "a + z + 2*a + 3*z + 5 + 6";

val thy = assoc_thy "Isac";

val subst = [(*TODO.WN071231 test Rewrite_Inst*)];

val rls = Test_simplify;

(* val rls = rls_p_33;      filter_appl_rews  ---> length 2

   val rls = norm_Poly;     filter_appl_rews  ---> length 1

   *)

if filter_appl_rews thy subst f rls =

   [Rewrite ("radd_commute", "?m + ?n = ?n + ?m"),

    Rewrite ("radd_assoc", "?m + ?n + ?k = ?m + (?n + ?k)"),

    Calculate "plus"] then () 

else error "rewtools.sml filter_appl_rews a + z + 2*a + 3*z + 5 + 6";

============ inhibit exn 110718 ==============================================*)

"----------- fun is_contained_in ------------------------";

"----------- fun is_contained_in ------------------------";

"----------- fun is_contained_in ------------------------";

val r1 = Thm ("real_diff_minus",num_str @{thm real_diff_minus});

if contains_rule r1 Test_simplify then ()

else error "rewtools.sml contains_rule Thm";

val r1 = Calc ("Groups.plus_class.plus", eval_binop "#add_");

if contains_rule r1 Test_simplify then ()

else error "rewtools.sml contains_rule Calc";

val r1 = Calc ("Groups.minus_class.minus", eval_binop "#add_");

if not (contains_rule r1 Test_simplify) then ()

else error "rewtools.sml contains_rule Calc";