wneuper/isa: test/Tools/isac/Interpret/error-fill-pattern.sml@64eecda7a8fb (annotated)

walther@59860	1	(* Title: "Interpret/error-fill-pattern.sml"
akargl@42176	2	Author: Walther Neuper 060225,
neuper@37906	3	(c) due to copyright terms
walther@59860	4
walther@59860	5	Strange ERROR "Undefined fact: "all_left""
neuper@37906	6	*)
neuper@37906	7
neuper@37906	8	"-----------------------------------------------------------------";
neuper@37906	9	"table of contents -----------------------------------------------";
neuper@37906	10	"-----------------------------------------------------------------";
neuper@37906	11	"appendForm with miniscript with mini-subpbl:";
neuper@37906	12	"--------- appendFormula: on Res + equ_nrls ----------------------";
neuper@37906	13	"--------- appendFormula: on Frm + equ_nrls ----------------------";
neuper@37906	14	"--------- appendFormula: on Res + NO deriv ----------------------";
neuper@37906	15	"--------- appendFormula: on Res + late deriv [x = 3 + -2]---///--";
neuper@37906	16	"replaceForm with miniscript with mini-subpbl:";
neuper@37906	17	"--------- replaceFormula: on Res + = ----------------------------";
neuper@37906	18	"--------- replaceFormula: on Res + = 1st Nd ---------------------";
neuper@37906	19	"--------- replaceFormula: on Frm + = 1st Nd ---------------------";
neuper@37906	20	"--------- replaceFormula: cut calculation -----------------------";
neuper@37906	21	"--------- replaceFormula: cut calculation -----------------------";
neuper@37906	22	(* 040307 copied from informtest.sml ... old versions
neuper@37906	23	"--------- maximum-example, UC: Modeling / modifyCalcHead --------";*)
neuper@37906	24	"--------- syntax error ------------------------------------------";
neuper@37906	25	"CAS-command:";
neuper@37906	26	"--------- CAS-command on ([],Pbl) -------------------------------";
neuper@37906	27	"--------- CAS-command on ([],Pbl) FE-interface ------------------";
walther@59772	28	"--------- locate_input_term [rational,simplification] ----------------------";
neuper@37906	29	"--------- Take as 1st tac, start with <NEW> (CAS input) ---------";
neuper@37906	30	"--------- Take as 1st tac, start from exp -----------------------";
walther@59845	31	"--------- implicit_take, start with <NEW> (CAS input) ---------------";
neuper@42423	32	"--------- build fun check_err_patt ------------------------------";
neuper@42426	33	"--------- build fun check_err_patt ?bdv -------------------------";
walther@59844	34	"--------- build fun Error_Pattern.check_for ------------------------";
walther@59844	35	"--------- embed fun Error_Pattern.check_for ------------------------";
neuper@48895	36	"--------- build fun get_fillpats --------------------------------";
neuper@42430	37	"--------- embed fun find_fillpatterns ---------------------------";
neuper@48895	38	"--------- build fun is_exactly_equal, inputFillFormula ----------";
wneuper@59262	39	"--------- fun appl_adds -----------------------------------------";
wneuper@59264	40	"--------- fun concat_deriv --------------------------------------";
wneuper@59264	41	"--------- handle an input formula -------------------------------";
wneuper@59264	42	"--------- fun dropwhile' ----------------------------------------";
neuper@37906	43	"-----------------------------------------------------------------";
neuper@37906	44	"-----------------------------------------------------------------";
neuper@37906	45	"-----------------------------------------------------------------";
neuper@37906	46
neuper@37906	47
neuper@37906	48	"--------- appendFormula: on Res + equ_nrls ----------------------";
neuper@37906	49	"--------- appendFormula: on Res + equ_nrls ----------------------";
neuper@37906	50	"--------- appendFormula: on Res + equ_nrls ----------------------";
neuper@48790	51	val Prog sc = (#scr o get_met) ["Test","squ-equ-test-subpbl1"];
neuper@37906	52	(writeln o term2str) sc;
neuper@48790	53	val Prog sc = (#scr o get_met) ["Test","solve_linear"];
neuper@37906	54	(writeln o term2str) sc;
neuper@37906	55
s1210629013@55445	56	reset_states ();
neuper@41970	57	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	58	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	59	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	60	Iterator 1; moveActiveRoot 1;
wneuper@59248	61	autoCalculate 1 CompleteCalcHead;
walther@59747	62	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	63	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
neuper@37906	64
wneuper@59123	65	appendFormula 1 "-2 * 1 + (1 + x) = 0" (\|> Future.join); refFormula 1 (get_pos 1 1);
neuper@37906	66	val ((pt,_),_) = get_calc 1;
wneuper@59279	67	val str = pr_ctree pr_short pt;
neuper@48895	68	if str =
neuper@48895	69	(". ----- pblobj -----\n" ^
neuper@48895	70	"1. x + 1 = 2\n" ^
neuper@48895	71	"2. x + 1 + -1 * 2 = 0\n" ^
neuper@48895	72	"2.1. x + 1 + -1 * 2 = 0\n" ^
neuper@48895	73	"2.2. 1 + x + -1 * 2 = 0\n" ^
neuper@48895	74	"2.3. 1 + (x + -1 * 2) = 0\n" ^
neuper@48895	75	"2.4. 1 + (x + -2) = 0\n" ^
neuper@48895	76	"2.5. 1 + (x + -2 * 1) = 0\n" ^
neuper@48895	77	"2.6. 1 + x + -2 * 1 = 0\n" ) then ()
neuper@48895	78	else error "inform.sml: diff.behav.appendFormula: on Res + equ 1";
neuper@37906	79
neuper@37906	80	moveDown 1 ([ ],Pbl); refFormula 1 ([1],Frm); (x + 1 = 2)
neuper@37906	81	moveDown 1 ([1],Frm); refFormula 1 ([1],Res); (x + 1 + -1 2 = 0*)
neuper@37906	82
neuper@37906	83	(the seven steps of detailed derivation)
neuper@37906	84	moveDown 1 ([1 ],Res); refFormula 1 ([2,1],Frm);
neuper@37906	85	moveDown 1 ([2,1],Frm); refFormula 1 ([2,1],Res);
neuper@37906	86	moveDown 1 ([2,1],Res); refFormula 1 ([2,2],Res);
neuper@37906	87	moveDown 1 ([2,2],Res); refFormula 1 ([2,3],Res);
neuper@37906	88	moveDown 1 ([2,3],Res); refFormula 1 ([2,4],Res);
neuper@37906	89	moveDown 1 ([2,4],Res); refFormula 1 ([2,5],Res);
neuper@37906	90	moveDown 1 ([2,5],Res); refFormula 1 ([2,6],Res);
neuper@37906	91	val ((pt,_),_) = get_calc 1;
neuper@37906	92	if "-2 * 1 + (1 + x) = 0" = term2str (fst (get_obj g_result pt [2,6])) then()
neuper@38031	93	else error "inform.sml: diff.behav.appendFormula: on Res + equ 2";
neuper@37906	94
neuper@37906	95	fetchProposedTactic 1; (takes Iterator 1 _1_)
neuper@48895	96	(* <ERROR> error in kernel </ERROR> ALREADY IN 2009-2*)
neuper@48895	97	(*========== inhibit exn WN1130701 broken at Isabelle2002 --> 2009-2 =============
neuper@37906	98	val (_,(tac,_,_)::_) = get_calc 1;
neuper@37906	99	if tac = Rewrite_Set "Test_simplify" then ()
neuper@38031	100	else error "inform.sml: diff.behav.appendFormula: on Res + equ 3";
neuper@48895	101	============ inhibit exn WN1130701 broken at Isabelle2002 --> 2009-2 ============*)
akargl@42209	102
wneuper@59248	103	autoCalculate 1 CompleteCalc;
neuper@37906	104	val ((pt,_),_) = get_calc 1;
neuper@37906	105	if "[x = 1]" = term2str (fst (get_obj g_result pt [])) then ()
neuper@38031	106	else error "inform.sml: diff.behav.appendFormula: on Res + equ 4";
wneuper@59248	107	(* autoCalculate 1 CompleteCalc;
neuper@37906	108	val ((pt,p),_) = get_calc 1;
neuper@37906	109	(writeln o istates2str) (get_obj g_loc pt [ ]);
neuper@37906	110	(writeln o istates2str) (get_obj g_loc pt [1]);
neuper@37906	111	(writeln o istates2str) (get_obj g_loc pt [2]);
neuper@37906	112	(writeln o istates2str) (get_obj g_loc pt [3]);
neuper@37906	113	(writeln o istates2str) (get_obj g_loc pt [3,1]);
neuper@37906	114	(writeln o istates2str) (get_obj g_loc pt [3,2]);
neuper@37906	115	(writeln o istates2str) (get_obj g_loc pt [4]);
neuper@37906	116
neuper@37906	117	*)
neuper@37906	118	"----------------------------------------------------------";
akargl@42176	119
wneuper@59592	120	val fod = make_deriv (@{theory "Isac_Knowledge"}) Atools_erls
walther@59852	121	((#rules o Rule_Set.rep) Test_simplify)
akargl@42176	122	(sqrt_right false (@{theory "Pure"})) NONE
neuper@37906	123	(str2term "x + 1 + -1 * 2 = 0");
neuper@37906	124	(writeln o trtas2str) fod;
neuper@37906	125
wneuper@59592	126	val ifod = make_deriv (@{theory "Isac_Knowledge"}) Atools_erls
walther@59852	127	((#rules o Rule_Set.rep) Test_simplify)
akargl@42176	128	(sqrt_right false (@{theory "Pure"})) NONE
neuper@37906	129	(str2term "-2 * 1 + (1 + x) = 0");
neuper@37906	130	(writeln o trtas2str) ifod;
akargl@42176	131	fun equal (_,_,(t1, _)) (_,_,(t2, _)) = t1 = t2;
neuper@37906	132	val (fod', rifod') = dropwhile' equal (rev fod) (rev ifod);
neuper@37906	133	val der = fod' @ (map rev_deriv' rifod');
neuper@37906	134	(writeln o trtas2str) der;
neuper@37906	135	"----------------------------------------------------------";
wneuper@59252	136	DEconstrCalcTree 1;
neuper@37906	137
neuper@37906	138	"--------- appendFormula: on Frm + equ_nrls ----------------------";
neuper@37906	139	"--------- appendFormula: on Frm + equ_nrls ----------------------";
neuper@37906	140	"--------- appendFormula: on Frm + equ_nrls ----------------------";
s1210629013@55445	141	reset_states ();
neuper@41970	142	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	143	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	144	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	145	Iterator 1; moveActiveRoot 1;
wneuper@59248	146	autoCalculate 1 CompleteCalcHead;
walther@59747	147	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1) (x + 1 = 2);
wneuper@59123	148	appendFormula 1 "2+ -1 + x = 2" (\|> Future.join); refFormula 1 (get_pos 1 1);
neuper@37906	149
neuper@37906	150	moveDown 1 ([],Pbl); refFormula 1 ([1],Frm) (x + 1 = 2);
neuper@37906	151
neuper@37906	152	moveDown 1 ([1 ],Frm); refFormula 1 ([1,1],Frm);
neuper@37906	153	moveDown 1 ([1,1],Frm); refFormula 1 ([1,1],Res);
neuper@37906	154	moveDown 1 ([1,1],Res); refFormula 1 ([1,2],Res);
neuper@37906	155	moveDown 1 ([1,2],Res); refFormula 1 ([1,3],Res);
neuper@37906	156	moveDown 1 ([1,3],Res); refFormula 1 ([1,4],Res);
neuper@37906	157	moveDown 1 ([1,4],Res); refFormula 1 ([1,5],Res);
neuper@37906	158	moveDown 1 ([1,5],Res); refFormula 1 ([1,6],Res);
neuper@37906	159	val ((pt,_),_) = get_calc 1;
neuper@37906	160	if "2 + -1 + x = 2" = term2str (fst (get_obj g_result pt [1,6])) then()
neuper@38031	161	else error "inform.sml: diff.behav.appendFormula: on Frm + equ 1";
neuper@37906	162
neuper@37906	163	fetchProposedTactic 1; (takes Iterator 1 _1_)
neuper@37906	164	val (_,(tac,_,_)::_) = get_calc 1;
wneuper@59253	165	case tac of Rewrite_Set "norm_equation" => ()
wneuper@59253	166	\| _ => error "inform.sml: diff.behav.appendFormula: on Frm + equ 2";
wneuper@59248	167	autoCalculate 1 CompleteCalc;
neuper@37906	168	val ((pt,_),_) = get_calc 1;
neuper@37906	169	if "[x = 1]" = term2str (fst (get_obj g_result pt [])) then ()
neuper@38031	170	else error "inform.sml: diff.behav.appendFormula: on Frm + equ 3";
wneuper@59252	171	DEconstrCalcTree 1;
neuper@37906	172
neuper@37906	173	"--------- appendFormula: on Res + NO deriv ----------------------";
neuper@37906	174	"--------- appendFormula: on Res + NO deriv ----------------------";
neuper@37906	175	"--------- appendFormula: on Res + NO deriv ----------------------";
s1210629013@55445	176	reset_states ();
neuper@41970	177	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	178	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	179	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	180	Iterator 1; moveActiveRoot 1;
wneuper@59248	181	autoCalculate 1 CompleteCalcHead;
walther@59747	182	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	183	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
neuper@37906	184
wneuper@59123	185	appendFormula 1 "x = 2" (\|> Future.join);
neuper@37906	186	val ((pt,p),_) = get_calc 1;
wneuper@59279	187	val str = pr_ctree pr_short pt;
neuper@37906	188	writeln str;
neuper@37906	189	if str = ". ----- pblobj -----\n1. x + 1 = 2\n" andalso p = ([1], Res)
neuper@37906	190	then ()
neuper@38031	191	else error "inform.sml: diff.behav.appendFormula: Res + NOder 1";
neuper@37906	192
neuper@37906	193	fetchProposedTactic 1;
neuper@37906	194	val (_,(tac,_,_)::_) = get_calc 1;
wneuper@59253	195	case tac of Rewrite_Set "Test_simplify" => ()
wneuper@59253	196	\| _ => error "inform.sml: diff.behav.appendFormula: Res + NOder 2";
wneuper@59248	197	autoCalculate 1 CompleteCalc;
neuper@37906	198	val ((pt,_),_) = get_calc 1;
neuper@37906	199	if "[x = 1]" = term2str (fst (get_obj g_result pt [])) then ()
neuper@38031	200	else error "inform.sml: diff.behav.appendFormula: on Frm + equ 3";
wneuper@59252	201	DEconstrCalcTree 1;
neuper@37906	202
neuper@37906	203	"--------- appendFormula: on Res + late deriv --------------------";
neuper@37906	204	"--------- appendFormula: on Res + late deriv --------------------";
neuper@37906	205	"--------- appendFormula: on Res + late deriv --------------------";
walther@59794	206	(*cp with "fun me" to test/../lucas-interpreter.sml:
walther@59794	207	re-build: fun locate_input_term ---------------------------------------------------";
walther@59794	208	*)
s1210629013@55445	209	reset_states ();
neuper@41970	210	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	211	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	212	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	213	Iterator 1; moveActiveRoot 1;
wneuper@59248	214	autoCalculate 1 CompleteCalcHead;
walther@59747	215	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	216	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
neuper@37906	217
wneuper@59123	218	appendFormula 1 "x = 1" (\|> Future.join);
neuper@37906	219	val ((pt,p),_) = get_calc 1;
wneuper@59279	220	val str = pr_ctree pr_short pt;
neuper@37906	221	writeln str;
neuper@37906	222	if str = ". ----- pblobj -----\n1. x + 1 = 2\n2. x + 1 + -1 * 2 = 0\n3. ----- pblobj -----\n3.1. -1 + x = 0\n3.2. x = 0 + -1 * -1\n3.2.1. x = 0 + -1 * -1\n3.2.2. x = 0 + 1\n" andalso p = ([3,2], Res)
neuper@37906	223	then () (finds 1 step too early: ([3,2], Res) "x = 1" also by script !!!)
neuper@38031	224	else error "inform.sml: diff.behav.appendFormula: Res + late d 1";
akargl@42176	225
neuper@37906	226	fetchProposedTactic 1;
neuper@37906	227	val (_,(tac,_,_)::_) = get_calc 1;
wneuper@59253	228	case tac of Check_Postcond ["LINEAR", "univariate", "equation", "test"] => ()
wneuper@59253	229	\| _ => error "inform.sml: diff.behav.appendFormula: Res + late d 2";
wneuper@59248	230	autoCalculate 1 CompleteCalc;
neuper@37906	231	val ((pt,_),_) = get_calc 1;
neuper@37906	232	if "[x = 1]" = term2str (fst (get_obj g_result pt [])) then ()
neuper@38031	233	else error "inform.sml: diff.behav.appendFormula: Res + late d 3";
wneuper@59252	234	DEconstrCalcTree 1;
neuper@37906	235
neuper@37906	236	"--------- appendFormula: on Res + late deriv [x = 3 + -2]---///--";
neuper@37906	237	"--------- appendFormula: on Res + late deriv [x = 3 + -2]---///--";
neuper@37906	238	"--------- appendFormula: on Res + late deriv [x = 3 + -2]---///--";
s1210629013@55445	239	reset_states ();
neuper@41970	240	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	241	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	242	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	243	Iterator 1; moveActiveRoot 1;
wneuper@59248	244	autoCalculate 1 CompleteCalcHead;
walther@59747	245	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	246	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
wneuper@59123	247	appendFormula 1 "[x = 3 + -21]" (\|> Future.join*);
neuper@37906	248	val ((pt,p),_) = get_calc 1;
wneuper@59279	249	val str = pr_ctree pr_short pt;
neuper@37906	250	writeln str;
neuper@37906	251	if str=". ----- pblobj -----\n1. x + 1 = 2\n2. x + 1 + -1 * 2 = 0\n3. ----- pblobj -----\n3.1. -1 + x = 0\n3.2. x = 0 + -1 * -1\n4. [x = 1]\n4.1. [x = 1]\n4.2. [x = -2 + 3]\n4.3. [x = 3 + -2]\n" then ()
neuper@38031	252	else error "inform.sml: diff.behav.appendFormula: Res + latEE 1";
wneuper@59248	253	autoCalculate 1 CompleteCalc;
neuper@37906	254	val ((pt,p),_) = get_calc 1;
neuper@37906	255	if "[x = 3 + -2 * 1]" = term2str (fst (get_obj g_result pt [])) then ()
neuper@37906	256	(* ~~~~~~~~~~ simplify as last step in any script ?!*)
neuper@38031	257	else error "inform.sml: diff.behav.appendFormula: Res + latEE 2";
wneuper@59252	258	DEconstrCalcTree 1;
neuper@37906	259
neuper@37906	260	"--------- replaceFormula: on Res + = ----------------------------";
neuper@37906	261	"--------- replaceFormula: on Res + = ----------------------------";
neuper@37906	262	"--------- replaceFormula: on Res + = ----------------------------";
s1210629013@55445	263	reset_states ();
neuper@41970	264	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	265	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	266	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	267	Iterator 1; moveActiveRoot 1;
wneuper@59248	268	autoCalculate 1 CompleteCalcHead;
walther@59747	269	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	270	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
walther@59747	271	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(-1 + x);
neuper@37906	272
neuper@37906	273	replaceFormula 1 "-2 * 1 + (1 + x) = 0"; refFormula 1 (get_pos 1 1);
neuper@37906	274	val ((pt,_),_) = get_calc 1;
wneuper@59279	275	val str = pr_ctree pr_short pt;
akargl@42176	276
neuper@48895	277	(* before AK110725 this was
neuper@48895	278	". ----- pblobj -----\n
neuper@48895	279	1. x + 1 = 2\n
neuper@48895	280	2. x + 1 + -1 * 2 = 0\n
neuper@48895	281	2.1. x + 1 + -1 * 2 = 0\n
neuper@48895	282	2.2. 1 + x + -1 * 2 = 0\n
neuper@48895	283	2.3. 1 + (x + -1 * 2) = 0\n
neuper@48895	284	2.4. 1 + (x + -2) = 0\n
neuper@48895	285	2.5. 1 + (x + -2 * 1) = 0\n
neuper@48895	286	2.6. 1 + x + -2 * 1 = 0\n";
neuper@48895	287	*)
neuper@48895	288	if str =
neuper@48895	289	". ----- pblobj -----\n"^
neuper@48895	290	"1. x + 1 = 2\n"^
neuper@48895	291	"2. x + 1 + -1 * 2 = 0\n"^
neuper@48895	292	"2.1. x + 1 + -1 * 2 = 0\n"^
neuper@48895	293	"2.2. 1 + x + -1 * 2 = 0\n"^
neuper@48895	294	"2.3. 1 + (x + -1 * 2) = 0\n"^
neuper@48895	295	"2.4. 1 + (x + -2) = 0\n"^
neuper@48895	296	"2.5. 1 + (x + -2 * 1) = 0\n"^
neuper@48895	297	"2.6. 1 + x + -2 * 1 = 0\n" then()
neuper@48895	298	else error "inform.sml: diff.behav.replaceFormula: on Res += 1";
akargl@42176	299
wneuper@59248	300	autoCalculate 1 CompleteCalc;
akargl@42176	301	val ((pt,pos as (p,_)),_) = get_calc 1;
akargl@42176	302	if pos = ([],Res) andalso "[x = 1]" = (term2str o fst) (get_obj g_result pt p) then()
neuper@38031	303	else error "inform.sml: diff.behav.replaceFormula: on Res + = 2";
wneuper@59252	304	DEconstrCalcTree 1;
neuper@37906	305
neuper@37906	306	"--------- replaceFormula: on Res + = 1st Nd ---------------------";
neuper@37906	307	"--------- replaceFormula: on Res + = 1st Nd ---------------------";
neuper@37906	308	"--------- replaceFormula: on Res + = 1st Nd ---------------------";
s1210629013@55445	309	reset_states ();
neuper@41970	310	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	311	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	312	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	313	Iterator 1; moveActiveRoot 1;
wneuper@59248	314	autoCalculate 1 CompleteCalcHead;
walther@59747	315	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	316	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
neuper@37906	317
neuper@37906	318	replaceFormula 1 "x + 1 = 4 + -2"; refFormula 1 (get_pos 1 1);
neuper@37906	319	val ((pt,_),_) = get_calc 1;
wneuper@59279	320	val str = pr_ctree pr_short pt;
neuper@37906	321	writeln str;
akargl@42176	322	if str= ". ----- pblobj -----\n1. x + 1 = 2\n1.1. x + 1 = 2\n1.2. 1 + x = 2\n1.3. 1 + x = -2 + 4\n1.4. x + 1 = -2 + 4\n" then ()
neuper@38031	323	else error "inform.sml: diff.behav.replaceFormula: on Res 1 + = 1";
wneuper@59248	324	autoCalculate 1 CompleteCalc;
akargl@42176	325	val ((pt,pos as (p,_)),_) = get_calc 1;
akargl@42176	326	if pos = ([],Res) andalso "[x = 1]" = (term2str o fst)(get_obj g_result pt p) then()
neuper@38031	327	else error "inform.sml: diff.behav.replaceFormula: on Res + = 2";
wneuper@59252	328	DEconstrCalcTree 1;
neuper@37906	329
neuper@37906	330	"--------- replaceFormula: on Frm + = 1st Nd ---------------------";
neuper@37906	331	"--------- replaceFormula: on Frm + = 1st Nd ---------------------";
neuper@37906	332	"--------- replaceFormula: on Frm + = 1st Nd ---------------------";
s1210629013@55445	333	reset_states ();
neuper@41970	334	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	335	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	336	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	337	Iterator 1; moveActiveRoot 1;
wneuper@59248	338	autoCalculate 1 CompleteCalcHead;
walther@59747	339	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
neuper@37906	340
neuper@37906	341	replaceFormula 1 "x + 1 = 4 + -2"; refFormula 1 (get_pos 1 1);
neuper@37906	342	val ((pt,_),_) = get_calc 1;
wneuper@59279	343	val str = pr_ctree pr_short pt;
neuper@37906	344	writeln str;
akargl@42176	345	if str= ". ----- pblobj -----\n1. x + 1 = 2\n1.1. x + 1 = 2\n1.2. 1 + x = 2\n1.3. 1 + x = -2 + 4\n1.4. x + 1 = -2 + 4\n" then ()
neuper@38031	346	else error "inform.sml: diff.behav.replaceFormula: on Frm 1 + = 1";
wneuper@59248	347	autoCalculate 1 CompleteCalc;
akargl@42176	348	val ((pt,pos as (p,_)),_) = get_calc 1;
akargl@42176	349	if pos = ([],Res) andalso "[x = 1]" = (term2str o fst)(get_obj g_result pt p) then()
neuper@38031	350	else error "inform.sml: diff.behav.replaceFormula: on Frm 1 + = 2";
wneuper@59252	351	DEconstrCalcTree 1;
neuper@37906	352
neuper@37906	353	"--------- replaceFormula: cut calculation -----------------------";
neuper@37906	354	"--------- replaceFormula: cut calculation -----------------------";
neuper@37906	355	"--------- replaceFormula: cut calculation -----------------------";
s1210629013@55445	356	reset_states ();
neuper@41970	357	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	358	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	359	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	360	Iterator 1; moveActiveRoot 1;
wneuper@59248	361	autoCalculate 1 CompleteCalc;
neuper@37906	362	moveActiveRoot 1; moveActiveDown 1;
neuper@37906	363	if get_pos 1 1 = ([1], Frm) then ()
neuper@38031	364	else error "inform.sml: diff.behav. cut calculation 1";
neuper@37906	365
neuper@37906	366	replaceFormula 1 "x + 1 = 4 + -2"; refFormula 1 (get_pos 1 1);
neuper@37906	367	val ((pt,p),_) = get_calc 1;
wneuper@59279	368	val str = pr_ctree pr_short pt;
neuper@37906	369	writeln str;
neuper@37906	370	if p = ([1], Res) then ()
neuper@38031	371	else error "inform.sml: diff.behav. cut calculation 2";
neuper@37906	372
neuper@37906	373
neuper@37906	374	(* 040307 copied from informtest.sml; ... old version
neuper@37906	375	"---------------- maximum-example, UC: Modeling / modifyCalcHead -";
neuper@37906	376	"---------------- maximum-example, UC: Modeling / modifyCalcHead -";
neuper@37906	377	"---------------- maximum-example, UC: Modeling / modifyCalcHead -";
neuper@37906	378
neuper@37906	379	val p = ([],Pbl);
neuper@37906	380	val elems = ["fixedValues [r=Arbfix]","maximum A","valuesFor [a,b]",
neuper@37906	381	"relations [A=a*b, (a/2)^^^2 + (b/2)^^^2 = r^^^2]",
neuper@37906	382	"relations [A=a*b, (a/2)^^^2 + (b/2)^^^2 = r^^^2]",
wneuper@59582	383	"relations [A=ab, a/2=rsin alpha, b/2=r*cos (alpha::real)]",
neuper@37906	384	(^^^ these are the elements for the root-problem (in variants))
neuper@37906	385	(vvv these are elements required for subproblems)
neuper@37906	386	"boundVariable a","boundVariable b","boundVariable alpha",
neuper@37906	387	"interval {x::real. 0 <= x & x <= 2*r}",
neuper@37906	388	"interval {x::real. 0 <= x & x <= 2*r}",
neuper@37906	389	"interval {x::real. 0 <= x & x <= pi}",
neuper@37906	390	"errorBound (eps=(0::real))"]
neuper@37906	391	(specifying is not interesting for this example)
neuper@38058	392	val spec = ("DiffApp", ["maximum_of","function"],
neuper@37906	393	["DiffApp","max_by_calculus"]);
neuper@37906	394	(the empty model with descriptions for user-guidance by Model_Problem)
neuper@37906	395	val empty_model = [Given ["fixedValues []"],
neuper@37906	396	Find ["maximum", "valuesFor"],
neuper@37906	397	Relate ["relations []"]];
neuper@37906	398
neuper@37906	399
neuper@37906	400	(!!!!!!!!!!!!!!!!! DON'T USE me FOR FINDING nxt !!!!!!!!!!!!!!!!!!)
neuper@37906	401	val (p,_,f,nxt,_,pt) = CalcTreeTEST [(elems, spec)];
neuper@37906	402	(val nxt = ("Model_Problem", ...)
neuper@37924	403	val pbl = get_obj g_pbl pt (fst p); (writeln o (itms2str_ ctxt)) pbl;
neuper@37906	404
neuper@37906	405	val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@37906	406	(nxt = Add_Given "fixedValues [r = Arbfix]")
neuper@37924	407	val pbl = get_obj g_pbl pt (fst p); (writeln o (itms2str_ ctxt)) pbl;
neuper@37906	408	(*[
neuper@37906	409	(0 ,[] ,false ,#Given ,Inc fixedValues [] ,(??.empty, [])),
neuper@37906	410	(0 ,[] ,false ,#Find ,Inc maximum ,(??.empty, [])),
neuper@37906	411	(0 ,[] ,false ,#Find ,Inc valuesFor ,(??.empty, [])),
neuper@37906	412	(0 ,[] ,false ,#Relate ,Inc relations [] ,(??.empty, []))]*)
neuper@37906	413
neuper@37906	414	(the empty CalcHead is checked w.r.t the model and re-established as such)
neuper@37906	415	val (b,pt,ocalhd) = input_icalhd pt (p,"", empty_model, Pbl, e_spec);
neuper@37924	416	val pbl = get_obj g_pbl pt (fst p); (writeln o (itms2str_ ctxt)) pbl;
neuper@38031	417	if ocalhd2str ocalhd = "(Pbl, ??.empty, [\n(0 ,[] ,false ,#Given ,Inc fixedValues [] ,(??.empty, [])),\n(0 ,[] ,false ,#Find ,Inc maximum ,(??.empty, [])),\n(0 ,[] ,false ,#Find ,Inc valuesFor ,(??.empty, [])),\n(0 ,[] ,false ,#Relate ,Inc relations [] ,(??.empty, []))], [], \n(\"e_domID\", [\"e_pblID\"], [\"e_metID\"]) )" then () else error "informtest.sml: diff.behav. max 1";
neuper@37906	418
neuper@37906	419	(there is one input to the model (could be more))
neuper@37906	420	val (b,pt,ocalhd) =
neuper@37906	421	input_icalhd pt (p,"", [Given ["fixedValues [r=Arbfix]"],
neuper@37906	422	Find ["maximum", "valuesFor"],
neuper@37906	423	Relate ["relations"]], Pbl, e_spec);
neuper@37924	424	val pbl = get_obj g_pbl pt (fst p); (writeln o (itms2str_ ctxt)) pbl;
neuper@37906	425	if ocalhd2str ocalhd = "(Pbl, ??.empty, [\n(1 ,[1,2,3] ,true ,#Given ,Cor fixedValues [r = Arbfix] ,(fix_, [[r = Arbfix]])),\n(0 ,[] ,false ,#Find ,Inc maximum ,(??.empty, [])),\n(0 ,[] ,false ,#Find ,Inc valuesFor ,(??.empty, [])),\n(0 ,[] ,false ,#Relate ,Inc relations [] ,(??.empty, []))], [], \n(\"e_domID\", [\"e_pblID\"], [\"e_metID\"]) )" then ()
neuper@38031	426	else error "informtest.sml: diff.behav. max 2";
neuper@37906	427
neuper@37906	428	(*this input is complete in variant 3, but the ME doesn't recognize FIXXXXME
neuper@37906	429	val (b,pt''''',ocalhd) =
neuper@37906	430	input_icalhd pt (p,"", [Given ["fixedValues [r=Arbfix]"],
neuper@37906	431	Find ["maximum A", "valuesFor [a,b]"],
neuper@37906	432	Relate ["relations [A=ab, a/2=rsin alpha, \
neuper@37906	433	\b/2=r*cos alpha]"]], Pbl, e_spec);
neuper@37924	434	val pbl = get_obj g_pbl pt''''' (fst p); (writeln o (itms2str_ ctxt)) pbl;
neuper@37906	435	if ocalhd2str ocalhd = ------------^^^^^^^^^^ missing !!!*)
neuper@37906	436
neuper@37906	437	(this input is complete in variant 1 (variant 3 does not work yet))
neuper@37906	438	val (b,pt''''',ocalhd) =
neuper@37906	439	input_icalhd pt (p,"", [Given ["fixedValues [r=Arbfix]"],
neuper@37906	440	Find ["maximum A", "valuesFor [a,b]"],
neuper@37906	441	Relate ["relations [A=a*b, \
neuper@37906	442	\(a/2)^^^2 + (b/2)^^^2 = r^^^2]"]],
neuper@37906	443	Pbl, e_spec);
neuper@37924	444	val pbl = get_obj g_pbl pt''''' (fst p); (writeln o (itms2str_ ctxt)) pbl;
neuper@37906	445
neuper@37906	446	modifycalcheadOK2xml 111 (bool2str b) ocalhd;
neuper@37906	447	*)
wneuper@59252	448	DEconstrCalcTree 1;
neuper@37906	449
neuper@37906	450	"--------- syntax error ------------------------------------------";
neuper@37906	451	"--------- syntax error ------------------------------------------";
neuper@37906	452	"--------- syntax error ------------------------------------------";
s1210629013@55445	453	reset_states ();
neuper@41970	454	CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"],
neuper@41970	455	("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906	456	["Test","squ-equ-test-subpbl1"]))];
neuper@37906	457	Iterator 1; moveActiveRoot 1;
wneuper@59248	458	autoCalculate 1 CompleteCalcHead;
walther@59747	459	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 = 2)
walther@59747	460	autoCalculate 1 (Steps 1); refFormula 1 (get_pos 1 1);(x + 1 + -1 2 = 0*);
neuper@37906	461
neuper@37906	462	appendFormula 1 " x - "; (<ERROR> syntax error in ' x - ' </ERROR>)
neuper@37906	463	val ((pt,_),_) = get_calc 1;
wneuper@59279	464	val str = pr_ctree pr_short pt;
neuper@37906	465	writeln str;
neuper@37906	466	if str = ". ----- pblobj -----\n1. x + 1 = 2\n" then ()
neuper@38031	467	else error "inform.sml: diff.behav.appendFormula: syntax error";
wneuper@59252	468	DEconstrCalcTree 1;
neuper@37906	469
neuper@37906	470	"--------- CAS-command on ([],Pbl) -------------------------------";
neuper@37906	471	"--------- CAS-command on ([],Pbl) -------------------------------";
neuper@37906	472	"--------- CAS-command on ([],Pbl) -------------------------------";
neuper@37906	473	val (p,_,f,nxt,_,pt) =
neuper@37906	474	CalcTreeTEST [([], ("e_domID", ["e_pblID"], ["e_metID"]))];
neuper@37906	475	val ifo = "solve(x+1=2,x)";
walther@59798	476	val (_,(_,c,(pt,p))) = Step_Solve.by_term (pt,p) "solve(x+1=2,x)";
neuper@37906	477	show_pt pt;
walther@59749	478	val nxt = (Apply_Method ["Test","squ-equ-test-subpbl1"]);
neuper@37906	479	val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@37906	480	if p = ([1], Frm) andalso f2str f = "x + 1 = 2" then ()
neuper@38031	481	else error "inform.sml: diff.behav. CAScmd ([],Pbl)";
wneuper@59252	482	DEconstrCalcTree 1;
neuper@37906	483
neuper@37906	484	"--------- CAS-command on ([],Pbl) FE-interface ------------------";
neuper@37906	485	"--------- CAS-command on ([],Pbl) FE-interface ------------------";
neuper@37906	486	"--------- CAS-command on ([],Pbl) FE-interface ------------------";
s1210629013@55445	487	reset_states ();
neuper@37906	488	CalcTree [([], ("e_domID", ["e_pblID"], ["e_metID"]))];
neuper@37906	489	Iterator 1;
neuper@37906	490	moveActiveRoot 1;
neuper@37906	491	replaceFormula 1 "solve(x+1=2,x)";
wneuper@59248	492	autoCalculate 1 CompleteCalc;
neuper@37906	493	val ((pt,p),_) = get_calc 1;
neuper@37906	494	show_pt pt;
neuper@37906	495	if p = ([], Res) then ()
neuper@38031	496	else error "inform.sml: diff.behav. CAScmd ([],Pbl) FE-interface";
wneuper@59252	497	DEconstrCalcTree 1;
neuper@37906	498
wneuper@59562	499	"--------- inform [rational,simplification] ----------------------";
wneuper@59562	500	"--------- inform [rational,simplification] ----------------------";
wneuper@59562	501	"--------- inform [rational,simplification] ----------------------";
s1210629013@55445	502	reset_states ();
neuper@52105	503	CalcTree [(["Term (a * x / (b * x) + c * x / (d * x) + e / f)", "normalform N"],
neuper@52105	504	("Rational", ["rational", "simplification"], ["simplification", "of_rationals"]))];
neuper@37906	505	Iterator 1; moveActiveRoot 1;
wneuper@59248	506	autoCalculate 1 CompleteCalcHead;
neuper@52105	507
neuper@52105	508	"--- (-1) give a preview on the calculation without any input";
neuper@52105	509	(*
wneuper@59248	510	autoCalculate 1 CompleteCalc;
neuper@52105	511	val ((pt, p), _) = get_calc 1;
neuper@52105	512	show_pt pt;
neuper@52105	513	[
neuper@52105	514	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	515	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	516	(([1], Res), a / b + c / d + e / f), <--- (1) input arbitrary
neuper@52105	517	(([2], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	518	(([3], Res), (a * (d * f) + b * (c * f) + b * (d * e)) / (b * (d * f))),
neuper@52105	519	(([4], Res), (a * d * f + b * c * f + b * d * e) / (b * d * f)), <--- (2) input next
neuper@52105	520	(([], Res), (a * d * f + b * c * f + b * d * e) / (b * d * f))] <--- (3) is also final result
neuper@52105	521	EXAMPLE NOT OPTIMAL
neuper@52105	522	*)
neuper@52105	523	"--- (0) user input as the first step does not work, thus impdo at least 1 step";
walther@59747	524	autoCalculate 1 (Steps 1);
walther@59747	525	autoCalculate 1 (Steps 1);
neuper@52105	526	val ((pt, p), _) = get_calc 1;
neuper@52105	527	(*show_pt pt;
neuper@52105	528	[
neuper@52105	529	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	530	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	531	(([1], Res), a / b + c / d + e / f)]
neuper@52105	532	*)
neuper@52105	533	"--- (1) input an arbitrary next formula";
wneuper@59123	534	appendFormula 1 "((a * d) + (c * b)) / (b * d) + e / f" (\|> Future.join);
neuper@52105	535	val ((pt, p), _) = get_calc 1;
neuper@52105	536	(*show_pt pt;
neuper@52105	537	[
neuper@52105	538	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	539	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	540	(([1], Res), a / b + c / d + e / f),
neuper@52105	541	(([2,1], Frm), a / b + c / d + e / f),
neuper@52105	542	(([2,1], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	543	(([2,2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	544	(([2], Res), (a * d + c * b) / (b * d) + e / f)]
neuper@52105	545	*)
neuper@37906	546	val ((pt,p),_) = get_calc 1;
neuper@52105	547	if p = ([2], Res) andalso (length o children o (get_nd pt)) (fst p) = 2 then ()
neuper@38031	548	else error ("inform.sml: [rational,simplification] 1");
neuper@37906	549
neuper@52105	550	"--- (2) input the next formula that would be presented by mat-engine";
neuper@52105	551	(* generate a preview:
walther@59747	552	autoCalculate 1 (Steps 1);
neuper@52105	553	val ((pt, p), _) = get_calc 1;
neuper@52105	554	show_pt pt;
neuper@52105	555	[
neuper@52105	556	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	557	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	558	(([1], Res), a / b + c / d + e / f),
neuper@52105	559	(([2,1], Frm), a / b + c / d + e / f),
neuper@52105	560	(([2,1], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	561	(([2,2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	562	(([2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	563	(([3], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f))] <--- input this
neuper@52105	564	*)
wneuper@59123	565	appendFormula 1 "(b * d * e + b * c * f + a * d * f) / (b * d * f)" (\|> Future.join);
neuper@52105	566	val ((pt, p), _) = get_calc 1;
neuper@52105	567	(*show_pt pt;
neuper@52105	568	[
neuper@52105	569	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	570	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	571	(([1], Res), a / b + c / d + e / f),
neuper@52105	572	(([2,1], Frm), a / b + c / d + e / f),
neuper@52105	573	(([2,1], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	574	(([2,2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	575	(([2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	576	(([3], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f))]
neuper@52105	577	*)
neuper@52105	578	if p = ([3], Res) andalso (length o children o (get_nd pt)) (fst p) = 0 then ()
neuper@38031	579	else error ("inform.sml: [rational,simplification] 2");
neuper@37906	580
neuper@52105	581	"--- (3) input the exact final result";
wneuper@59123	582	appendFormula 1 "(b * d * e + b * c * f + a * d * f) / (b * d * f)" (\|> Future.join);
neuper@52105	583	val ((pt, p), _) = get_calc 1;
neuper@52105	584	(*show_pt pt;
neuper@52105	585	[
neuper@52105	586	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	587	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	588	(([1], Res), a / b + c / d + e / f),
neuper@52105	589	(([2,1], Frm), a / b + c / d + e / f),
neuper@52105	590	(([2,1], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	591	(([2,2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	592	(([2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	593	(([3], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	594	(([4,1], Frm), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	595	(([4,1], Res), (a * (d * f) + b * (c * f) + b * (d * e)) / (b * (d * f))),
neuper@52105	596	(([4,2], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	597	(([4], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f))]
neuper@52105	598	*)
neuper@52105	599	if p = ([4], Res) andalso (length o children o (get_nd pt)) (fst p) = 2 then ()
neuper@38031	600	else error ("inform.sml: [rational,simplification] 3");
neuper@52105	601
neuper@52105	602	"--- (4) finish the calculation + check the postcondition (in the future)";
wneuper@59248	603	autoCalculate 1 CompleteCalc;
neuper@52105	604	val ((pt, p), _) = get_calc 1;
neuper@52105	605	val (t, asm) = get_obj g_result pt [];
neuper@52105	606	if term2str t = "(a * d * f + b * c * f + b * d * e) / (b * d * f)" andalso
walther@59719	607	terms2str asm = "[\"b * d * f \<noteq> 0\",\"d \<noteq> 0\",\"b \<noteq> 0\",\"a * x / (b * x) + c * x / (d * x) + e / f is_ratpolyexp\"]"
neuper@52105	608	then () else error "inform [rational,simplification] changed at end";
neuper@52105	609	(*show_pt pt;
neuper@52105	610	[
neuper@52105	611	(([], Frm), Simplify (a * x / (b * x) + c * x / (d * x) + e / f)),
neuper@52105	612	(([1], Frm), a * x / (b * x) + c * x / (d * x) + e / f),
neuper@52105	613	(([1], Res), a / b + c / d + e / f),
neuper@52105	614	(([2,1], Frm), a / b + c / d + e / f),
neuper@52105	615	(([2,1], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	616	(([2,2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	617	(([2], Res), (a * d + c * b) / (b * d) + e / f),
neuper@52105	618	(([3], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	619	(([4,1], Frm), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	620	(([4,1], Res), (a * (d * f) + b * (c * f) + b * (d * e)) / (b * (d * f))),
neuper@52105	621	(([4,2], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	622	(([4], Res), (b * d * e + b * c * f + a * d * f) / (b * d * f)),
neuper@52105	623	(([5], Res), (a * (d * f) + b * (c * f) + b * (d * e)) / (b * (d * f))),
neuper@52105	624	(([6], Res), (a * d * f + b * c * f + b * d * e) / (b * d * f)),
neuper@52105	625	(([], Res), (a * d * f + b * c * f + b * d * e) / (b * d * f))]
neuper@52105	626	*)
wneuper@59252	627	DEconstrCalcTree 1;
neuper@37906	628
neuper@37906	629	"--------- Take as 1st tac, start with <NEW> (CAS input) ---------";
neuper@37906	630	"--------- Take as 1st tac, start with <NEW> (CAS input) ---------";
neuper@37906	631	"--------- Take as 1st tac, start with <NEW> (CAS input) ---------";
neuper@37906	632	val t = str2term "Diff (x^^^2 + x + 1, x)";
neuper@37906	633	case t of Const ("Diff.Diff", _) $ _ => ()
neuper@37906	634	\| _ => raise
neuper@37906	635	error "diff.sml behav.changed for CAS Diff (..., x)";
neuper@37906	636	atomty t;
neuper@37906	637	"-----------------------------------------------------------------";
s1210629013@55445	638	(1>)reset_states ();
neuper@37906	639	(2>)CalcTree [([], ("e_domID", ["e_pblID"], ["e_metID"]))];
neuper@37906	640	(3>)Iterator 1;moveActiveRoot 1;
neuper@37906	641	"----- here the Headline has been finished";
neuper@37906	642	(4>)moveActiveFormula 1 ([],Pbl);
neuper@37906	643	(5>)replaceFormula 1 "Diff (x^2 + x + 1, x)";
neuper@37906	644	val ((pt,_),_) = get_calc 1;
neuper@37906	645	val PblObj {probl, meth, spec, fmz, env, loc, ...} = get_obj I pt [];
neuper@37926	646	val NONE = env;
neuper@37926	647	val (SOME istate, NONE) = loc;
wneuper@59348	648	(default_print_depth 5;)
neuper@37906	649	writeln"-----------------------------------------------------------";
neuper@37906	650	spec;
neuper@37924	651	writeln (itms2str_ ctxt probl);
neuper@37924	652	writeln (itms2str_ ctxt meth);
walther@59844	653	writeln (Istate.string_of (fst istate));
neuper@37906	654
neuper@37906	655	refFormula 1 ([],Pbl) (--> correct CalcHead);
neuper@37906	656	(081016 NOT necessary (but leave it in Java):)
wneuper@59248	657	(6>)(completeCalcHead)autoCalculate 1 CompleteCalcHead;
neuper@37906	658	"----- here the CalcHead has been completed --- ONCE MORE ?????";
neuper@37906	659
neuper@37906	660	(*difference II*)
neuper@37906	661	val ((pt,p),_) = get_calc 1;
neuper@37906	662	(val p = ([], Pbl))
neuper@37906	663	val PblObj {probl, meth, spec, fmz, env, loc, ...} = get_obj I pt [];
neuper@37926	664	val NONE = env;
neuper@37926	665	val (SOME istate, NONE) = loc;
walther@59844	666	(default_print_depth 5;) writeln (Istate.string_of (fst istate)); (default_print_depth 3;)
wneuper@59583	667	(*Pstate ([],
neuper@37926	668	[], NONE,
neuper@37906	669	??.empty, Sundef, false)*)
wneuper@59348	670	(default_print_depth 5;) spec; (default_print_depth 3;)
wneuper@59592	671	(*("Isac_Knowledge",
neuper@37906	672	["derivative_of", "function"],
neuper@37906	673	["diff", "differentiate_on_R"]) : spec*)
neuper@37924	674	writeln (itms2str_ ctxt probl);
neuper@37906	675	(*[
neuper@37906	676	(1 ,[1] ,true ,#Given ,Cor functionTerm (x ^^^ 2 + x + 1) ,(f_, [x ^^^ 2 + x + 1])),
neuper@37906	677	(2 ,[1] ,true ,#Given ,Cor differentiateFor x ,(v_, [x])),
akargl@42176	678	(3 ,[1] ,true ,#Find ,Cor derivative f_'_f ,(f_'_f, [f_'_f]))]*)
neuper@37924	679	writeln (itms2str_ ctxt meth);
neuper@37906	680	(*[
neuper@37906	681	(1 ,[1] ,true ,#Given ,Cor functionTerm (x ^^^ 2 + x + 1) ,(f_, [x ^^^ 2 + x + 1])),
neuper@37906	682	(2 ,[1] ,true ,#Given ,Cor differentiateFor x ,(v_, [x])),
akargl@42176	683	(3 ,[1] ,true ,#Find ,Cor derivative f_'_f ,(f_'_f, [f_'_f]))]*)
neuper@37906	684	writeln"-----------------------------------------------------------";
neuper@37906	685	(7>)fetchProposedTactic 1 (--> Apply_Method);
neuper@37906	686	(WN081028 fixed <ERROR> helpless </ERROR> by inform returning ...(.,Met))
wneuper@59248	687	autoCalculate 1 CompleteCalc;
neuper@37906	688	val ((pt,p),_) = get_calc 1;
neuper@37906	689	val Form res = (#1 o pt_extract) (pt, ([],Res));
neuper@37906	690	show_pt pt;
neuper@37906	691	if p = ([], Res) andalso term2str res = "1 + 2 * x" then ()
neuper@38031	692	else error "diff.sml behav.changed for Diff (x^2 + x + 1, x)";
wneuper@59252	693	DEconstrCalcTree 1;
neuper@37906	694
neuper@37906	695	"--------- Take as 1st tac, start from exp -----------------------";
neuper@37906	696	"--------- Take as 1st tac, start from exp -----------------------";
neuper@37906	697	"--------- Take as 1st tac, start from exp -----------------------";
neuper@37906	698	(*the following input is copied from BridgeLog Java <==> SML,
neuper@37906	699	omitting unnecessary inputs*)
s1210629013@55445	700	(1>)reset_states ();
wneuper@59592	701	(2>)CalcTree [(["functionTerm (x^2 + x + 1)", "differentiateFor x", "derivative f_'_f"],("Isac_Knowledge",["derivative_of","function"],["diff","differentiate_on_R"]))];
neuper@37906	702	(3>)Iterator 1; moveActiveRoot 1;
neuper@37906	703
wneuper@59248	704	(6>)(completeCalcHead)autoCalculate 1 CompleteCalcHead;
neuper@37906	705	(*difference II*)
neuper@37906	706	val ((pt,_),_) = get_calc 1;
neuper@37906	707	val PblObj {probl, meth, spec, fmz, env, loc, ...} = get_obj I pt [];
neuper@37926	708	val NONE = env;
neuper@37926	709	val (SOME istate, NONE) = loc;
walther@59844	710	(default_print_depth 5;) writeln (Istate.string_of (fst istate)); (default_print_depth 3;)
wneuper@59583	711	(*Pstate ([],
neuper@37926	712	[], NONE,
neuper@37906	713	??.empty, Sundef, false)*)
wneuper@59348	714	(default_print_depth 5;) spec; (default_print_depth 3;)
wneuper@59592	715	(*("Isac_Knowledge",
neuper@37906	716	["derivative_of", "function"],
neuper@37906	717	["diff", "differentiate_on_R"]) : spec*)
neuper@37924	718	writeln (itms2str_ ctxt probl);
neuper@37906	719	(*[
neuper@37906	720	(1 ,[1] ,true ,#Given ,Cor functionTerm (x ^^^ 2 + x + 1) ,(f_, [x ^^^ 2 + x + 1])),
neuper@37906	721	(2 ,[1] ,true ,#Given ,Cor differentiateFor x ,(v_, [x])),
akargl@42176	722	(3 ,[1] ,true ,#Find ,Cor derivative f_'_f ,(f_'_f, [f_'_f]))]*)
neuper@37924	723	writeln (itms2str_ ctxt meth);
neuper@37906	724	(*[
neuper@37906	725	(1 ,[1] ,true ,#Given ,Cor functionTerm (x ^^^ 2 + x + 1) ,(f_, [x ^^^ 2 + x + 1])),
neuper@37906	726	(2 ,[1] ,true ,#Given ,Cor differentiateFor x ,(v_, [x])),
akargl@42176	727	(3 ,[1] ,true ,#Find ,Cor derivative f_'_f ,(f_'_f, [f_'_f]))]*)
neuper@37906	728	writeln"-----------------------------------------------------------";
neuper@37906	729	(7>)fetchProposedTactic 1 (--> Apply_Method);
walther@59747	730	autoCalculate 1 (Steps 1);
neuper@37906	731	val ((pt,p),_) = get_calc 1;
neuper@37906	732	val Form res = (#1 o pt_extract) (pt, p);
neuper@37906	733	if term2str res = "d_d x (x ^^^ 2 + x + 1)" then ()
neuper@38031	734	else error "diff.sml Diff (x^2 + x + 1, x) from exp";
wneuper@59252	735	DEconstrCalcTree 1;
neuper@37906	736
walther@59845	737	"--------- implicit_take, start with <NEW> (CAS input) ---------------";
walther@59845	738	"--------- implicit_take, start with <NEW> (CAS input) ---------------";
walther@59845	739	"--------- implicit_take, start with <NEW> (CAS input) ---------------";
s1210629013@55445	740	reset_states ();
neuper@37906	741	CalcTree [([], ("e_domID", ["e_pblID"], ["e_metID"]))];
neuper@37906	742	(*[[from sml: > @@@@@begin@@@@@
neuper@37906	743	[[from sml: 1
neuper@37906	744	[[from sml: <CALCTREE>
neuper@37906	745	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	746	[[from sml: </CALCTREE>
neuper@37906	747	[[from sml: @@@@@end@@@@@*)
neuper@37906	748	Iterator 1;
neuper@37906	749	(*[[from sml: > @@@@@begin@@@@@
neuper@37906	750	[[from sml: 1
neuper@37906	751	[[from sml: <ADDUSER>
neuper@37906	752	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	753	[[from sml: <USERID> 1 </USERID>
neuper@37906	754	[[from sml: </ADDUSER>
neuper@37906	755	[[from sml: @@@@@end@@@@@*)
neuper@37906	756	moveActiveRoot 1;
neuper@37906	757	(*[[from sml: > @@@@@begin@@@@@
neuper@37906	758	[[from sml: 1
neuper@37906	759	[[from sml: <CALCITERATOR>
neuper@37906	760	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	761	[[from sml: <POSITION>
neuper@37906	762	[[from sml: <INTLIST>
neuper@37906	763	[[from sml: </INTLIST>
neuper@37906	764	[[from sml: <POS> Pbl </POS>
neuper@37906	765	[[from sml: </POSITION>
neuper@37906	766	[[from sml: </CALCITERATOR>
neuper@37906	767	[[from sml: @@@@@end@@@@@*)
neuper@37906	768	getFormulaeFromTo 1 ([],Pbl) ([],Pbl) 0 false;
neuper@37906	769	(*[[from sml: > @@@@@begin@@@@@ STILL CORRECT
neuper@37906	770	[[from sml: 1
neuper@37906	771	[[from sml: <GETELEMENTSFROMTO>
neuper@37906	772	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	773	[[from sml: <FORMHEADS>
neuper@37906	774	[[from sml: <CALCFORMULA>
neuper@37906	775	[[from sml: <POSITION>
neuper@37906	776	[[from sml: <INTLIST>
neuper@37906	777	[[from sml: </INTLIST>
neuper@37906	778	[[from sml: <POS> Pbl </POS>
neuper@37906	779	[[from sml: </POSITION>
neuper@37906	780	[[from sml: <FORMULA>
neuper@37906	781	[[from sml: <MATHML>
neuper@37906	782	[[from sml: <ISA> ________________________________________________ </ISA>
neuper@37906	783	[[from sml: </MATHML>
neuper@37906	784	[[from sml:
neuper@37906	785	[[from sml: </FORMULA>
neuper@37906	786	[[from sml: </CALCFORMULA>
neuper@37906	787	[[from sml: </FORMHEADS>
neuper@37906	788	[[from sml: </GETELEMENTSFROMTO>
neuper@37906	789	[[from sml: @@@@@end@@@@@*)
neuper@37906	790	refFormula 1 ([],Pbl);
neuper@37906	791	(*[[from sml: > @@@@@begin@@@@@ STILL CORRECT
neuper@37906	792	[[from sml: 1
neuper@37906	793	[[from sml: <REFFORMULA>
neuper@37906	794	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	795	[[from sml: <CALCHEAD status = "incorrect">
neuper@37906	796	[[from sml: <POSITION>
neuper@37906	797	[[from sml: <INTLIST>
neuper@37906	798	[[from sml: </INTLIST>
neuper@37906	799	[[from sml: <POS> Pbl </POS>
neuper@37906	800	[[from sml: </POSITION>
neuper@37906	801	[[from sml: <HEAD>
neuper@37906	802	[[from sml: <MATHML>
neuper@37906	803	[[from sml: <ISA> Problem (e_domID, [e_pblID]) </ISA>
neuper@37906	804	[[from sml: </MATHML>
neuper@37906	805	[[from sml: </HEAD>
neuper@37906	806	[[from sml: <MODEL>
neuper@37906	807	[[from sml: <GIVEN> </GIVEN>
neuper@37906	808	[[from sml: <WHERE> </WHERE>
neuper@37906	809	[[from sml: <FIND> </FIND>
neuper@37906	810	[[from sml: <RELATE> </RELATE>
neuper@37906	811	[[from sml: </MODEL>
neuper@37906	812	[[from sml: <BELONGSTO> Pbl </BELONGSTO>
neuper@37906	813	[[from sml: <SPECIFICATION>
neuper@37906	814	[[from sml: <THEORYID> e_domID </THEORYID>
neuper@37906	815	[[from sml: <PROBLEMID>
neuper@37906	816	[[from sml: <STRINGLIST>
neuper@37906	817	[[from sml: <STRING> e_pblID </STRING>
neuper@37906	818	[[from sml: </STRINGLIST>
neuper@37906	819	[[from sml: </PROBLEMID>
neuper@37906	820	[[from sml: <METHODID>
neuper@37906	821	[[from sml: <STRINGLIST>
neuper@37906	822	[[from sml: <STRING> e_metID </STRING>
neuper@37906	823	[[from sml: </STRINGLIST>
neuper@37906	824	[[from sml: </METHODID>
neuper@37906	825	[[from sml: </SPECIFICATION>
neuper@37906	826	[[from sml: </CALCHEAD>
neuper@37906	827	[[from sml: </REFFORMULA>
neuper@37906	828	[[from sml: @@@@@end@@@@@*)
neuper@37906	829	moveActiveFormula 1 ([],Pbl);
neuper@37906	830	(*[[from sml: > @@@@@begin@@@@@
neuper@37906	831	[[from sml: 1
neuper@37906	832	[[from sml: <CALCITERATOR>
neuper@37906	833	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	834	[[from sml: <POSITION>
neuper@37906	835	[[from sml: <INTLIST>
neuper@37906	836	[[from sml: </INTLIST>
neuper@37906	837	[[from sml: <POS> Pbl </POS>
neuper@37906	838	[[from sml: </POSITION>
neuper@37906	839	[[from sml: </CALCITERATOR>
neuper@37906	840	[[from sml: @@@@@end@@@@@*)
neuper@37906	841	replaceFormula 1 "Simplify (1+2)";
neuper@37906	842	(*[[from sml: > @@@@@begin@@@@@
neuper@37906	843	[[from sml: 1
neuper@37906	844	[[from sml: <REPLACEFORMULA>
neuper@37906	845	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	846	[[from sml: <CALCCHANGED>
neuper@37906	847	[[from sml: <UNCHANGED>
neuper@37906	848	[[from sml: <INTLIST>
neuper@37906	849	[[from sml: </INTLIST>
neuper@37906	850	[[from sml: <POS> Pbl </POS>
neuper@37906	851	[[from sml: </UNCHANGED>
neuper@37906	852	[[from sml: <DELETED>
neuper@37906	853	[[from sml: <INTLIST>
neuper@37906	854	[[from sml: </INTLIST>
neuper@37906	855	[[from sml: <POS> Pbl </POS>
neuper@37906	856	[[from sml: </DELETED>
neuper@37906	857	[[from sml: <GENERATED>
neuper@37906	858	[[from sml: <INTLIST>
neuper@37906	859	[[from sml: </INTLIST>
neuper@37906	860	[[from sml: <POS> Met </POS> DIFFERENCE: Pbl
neuper@37906	861	[[from sml: </GENERATED>
neuper@37906	862	[[from sml: </CALCCHANGED>
neuper@37906	863	[[from sml: </REPLACEFORMULA>
neuper@37906	864	[[from sml: @@@@@end@@@@@*)
neuper@37906	865	getFormulaeFromTo 1 ([],Pbl) ([],Pbl) 0 false(* DIFFERENCE: Pbl*);
neuper@37906	866	(*@@@@@begin@@@@@
neuper@37906	867	1
neuper@37906	868	<GETELEMENTSFROMTO>
neuper@37906	869	<CALCID> 1 </CALCID>
neuper@37906	870	<FORMHEADS>
neuper@37906	871	<CALCFORMULA>
neuper@37906	872	<POSITION>
neuper@37906	873	<INTLIST>
neuper@37906	874	</INTLIST>
neuper@37906	875	<POS> Pbl </POS>
neuper@37906	876	</POSITION>
neuper@37906	877	<FORMULA>
neuper@37906	878	<MATHML>
neuper@37906	879	<ISA> Simplify (1 + 2) </ISA> WORKS !!!!!
neuper@37906	880	</MATHML>
neuper@37906	881	</FORMULA>
neuper@37906	882	</CALCFORMULA>
neuper@37906	883	</FORMHEADS>
neuper@37906	884	</GETELEMENTSFROMTO>
neuper@37906	885	@@@@@end@@@@@*)
neuper@37906	886	getFormulaeFromTo 1 ([],Pbl) ([],Met) 0 false;
neuper@37906	887	(*[[from sml: > @@@@@begin@@@@@
neuper@37906	888	[[from sml: 1
neuper@37906	889	[[from sml: <SYSERROR>
neuper@37906	890	[[from sml: <CALCID> 1 </CALCID>
neuper@37906	891	[[from sml: <ERROR> error in getFormulaeFromTo </ERROR>
neuper@37906	892	[[from sml: </SYSERROR>
neuper@37906	893	[[from sml: @@@@@end@@@@@*)
neuper@37906	894	(step into getFormulaeFromTo --- bug corrected...)
neuper@42423	895
neuper@42423	896	"--------- build fun check_err_patt ------------------------------";
neuper@42423	897	"--------- build fun check_err_patt ------------------------------";
neuper@42423	898	"--------- build fun check_err_patt ------------------------------";
neuper@42426	899	val subst = [(str2term "bdv", str2term "x")]: subst;
neuper@42423	900	val rls = norm_Rational
neuper@42423	901	val pat = parse_patt @{theory} "(?a + ?b)/(?b + ?c) = ?a / ?c";
neuper@42423	902	val (res, inf) = (str2term "(2 + 3)/(3 + 4)", str2term "2 / 4");
neuper@42423	903	val (res, inf) = (str2term "(2 + 3)/(3 + 4)", str2term "1 / 2");
neuper@42423	904
neuper@42426	905	val (res', _, _, rewritten) = (rewritten: the lhs of the pattern matches in res)
walther@59852	906	rew_sub thy 1 [] e_rew_ord Rule_Set.empty false [] (HOLogic.Trueprop $ pat) res;
neuper@42423	907	if rewritten then NONE else SOME "e_errpatID";
neuper@42423	908
neuper@42423	909	val norm_res = case rewrite_set_ (Isac()) false rls res' of
neuper@42423	910	NONE => res'
neuper@42423	911	\| SOME (norm_res, _) => norm_res
neuper@42423	912
neuper@42423	913	val norm_inf = case rewrite_set_ (Isac()) false rls inf of
neuper@42423	914	NONE => inf
neuper@42423	915	\| SOME (norm_inf, _) => norm_inf;
neuper@42423	916
neuper@42423	917	res' = inf;
neuper@42423	918	norm_res = norm_inf;
neuper@42423	919
neuper@42423	920	val pat = parse_patt @{theory} "(?a + ?b)/?a = ?b";
neuper@42423	921	val (res, inf) = (str2term "(2 + 3)/2", str2term "3");
neuper@42426	922	if check_err_patt (res, inf) subst ("errpatID", pat) rls = SOME "errpatID"
neuper@42423	923	then () else error "error patt example1 changed";
neuper@42423	924
neuper@42423	925	val pat = parse_patt @{theory} "(?a + ?b)/(?a + ?c) = ?b / ?c";
neuper@42423	926	val (res, inf) = (str2term "(2 + 3)/(2 + 4)", str2term "3 / 4");
neuper@42426	927	if check_err_patt (res, inf) subst ("errpatID", pat) rls = SOME "errpatID"
neuper@42423	928	then () else error "error patt example2 changed";
neuper@42423	929
neuper@42423	930	val pat = parse_patt @{theory} "(?a + ?b)/(?b + ?c) = ?a / ?c";
neuper@42423	931	val (res, inf) = (str2term "(2 + 3)/(3 + 4)", str2term "2 / 4");
neuper@42426	932	if check_err_patt (res, inf) subst ("errpatID", pat) rls = SOME "errpatID"
neuper@42423	933	then () else error "error patt example3 changed";
neuper@42423	934
neuper@42423	935	val inf = str2term "1 / 2";
neuper@42426	936	if check_err_patt (res, inf) subst ("errpatID", pat) rls = SOME "errpatID"
neuper@42423	937	then () else error "error patt example3 changed";
neuper@42423	938
neuper@42426	939	"--------- build fun check_err_patt ?bdv -------------------------";
neuper@42426	940	"--------- build fun check_err_patt ?bdv -------------------------";
neuper@42426	941	"--------- build fun check_err_patt ?bdv -------------------------";
neuper@42426	942	val subst = [(str2term "bdv", str2term "x")]: subst;
neuper@42426	943	val t = str2term "d_d x (x ^^^ 2 + sin (x ^^^ 4))";
neuper@42426	944	val SOME (t, _) = rewrite_set_inst_ thy false subst norm_diff t;
neuper@42426	945	if term2str t = "2 * x + cos (x ^^^ 4) * 4 * x ^^^ 3" then ()
neuper@42426	946	else error "build fun check_err_patt ?bdv changed 1";
neuper@42423	947
neuper@42426	948	val rls = norm_diff
neuper@42426	949	val pat = parse_patt @{theory} "d_d ?bdv (sin ?u) = cos (d_d ?bdv ?u)";
neuper@42426	950	val (res, inf) = (str2term "2 * x + d_d x (sin (x ^^^ 4))", str2term "2 * x + cos (4 * x ^^^ 3)");
neuper@42426	951
neuper@42426	952	val (res', _, _, rewritten) = (rewritten: the lhs of the pattern matches in res)
walther@59852	953	rew_sub thy 1 subst e_rew_ord Rule_Set.empty false [] (HOLogic.Trueprop $ pat) res;
neuper@42426	954	if term2str res' = "2 * x + cos (d_d x (x ^^^ 4))" andalso rewritten then ()
neuper@42426	955	else error "build fun check_err_patt ?bdv changed 2";
neuper@42426	956
neuper@42426	957	val norm_res = case rewrite_set_inst_ (Isac()) false subst rls res' of
neuper@42426	958	NONE => res'
neuper@42426	959	\| SOME (norm_res, _) => norm_res;
neuper@42426	960	if term2str norm_res = "2 * x + cos (4 * x ^^^ 3)" then ()
neuper@42426	961	else error "build fun check_err_patt ?bdv changed 3";
neuper@42426	962
neuper@42426	963	val norm_inf = case rewrite_set_inst_ (Isac()) false subst rls inf of
neuper@42426	964	NONE => inf
neuper@42426	965	\| SOME (norm_inf, _) => norm_inf;
neuper@42426	966	if term2str norm_inf = "2 * x + cos (4 * x ^^^ 3)" then ()
neuper@42426	967	else error "build fun check_err_patt ?bdv changed 4";
neuper@42426	968
neuper@42426	969	res' = inf;
neuper@42426	970	if norm_res = norm_inf then ()
neuper@42426	971	else error "build fun check_err_patt ?bdv changed 5";
neuper@42426	972
neuper@42426	973	if check_err_patt (res, inf) (subst: subst) ("errpatID": errpatID, pat) rls = SOME "errpatID"
neuper@42426	974	then () else error "error patt example1 changed";
neuper@42426	975
walther@59844	976	"--------- build fun Error_Pattern.check_for ------------------------";
walther@59844	977	"--------- build fun Error_Pattern.check_for ------------------------";
walther@59844	978	"--------- build fun Error_Pattern.check_for ------------------------";
neuper@42428	979	val (res, inf) =
neuper@42428	980	(str2term "d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))",
neuper@42428	981	str2term "d_d x (x ^^^ 2) + cos (4 * x ^^^ 3)");
neuper@48790	982	val {errpats, nrls = rls, scr = Prog prog, ...} = get_met ["diff", "differentiate_on_R"]
neuper@42426	983
neuper@42428	984	val env = [(str2term "v_v", str2term "x")];
neuper@42428	985	val errpats =
neuper@42428	986	[e_errpat, (generalised for testing)
neuper@42428	987	("chain-rule-diff-both",
neuper@42428	988	[parse_patt @{theory} "d_d ?bdv (sin ?u) = cos (d_d ?bdv ?u)",
neuper@42428	989	parse_patt @{theory} "d_d ?bdv (cos ?u) = - sin (d_d ?bdv ?u)",
neuper@42428	990	parse_patt @{theory} "d_d ?bdv (?u ^^^ ?n) = ?n * ?u ^^^ (?n - 1)",
neuper@42428	991	parse_patt @{theory} "d_d ?bdv (LogExp.ln ?u) = 1 / ?u",
neuper@42428	992	parse_patt @{theory} "d_d ?bdv (LogExp.ln ?u) = 1 / d_d ?bdv ?u"],
neuper@42428	993	[@{thm diff_sin_chain}, @{thm diff_cos_chain}, @{thm diff_pow_chain},
neuper@42428	994	@{thm diff_ln_chain}, @{thm diff_exp_chain}])]: errpat list;
walther@59844	995	case Error_Pattern.check_for (res, inf) (prog, env) (errpats, rls) of SOME _ => ()
walther@59844	996	\| NONE => error "Error_Pattern.check_for broken";
wneuper@59252	997	DEconstrCalcTree 1;
neuper@42428	998
walther@59844	999	"--------- embed fun Error_Pattern.check_for ------------------------";
walther@59844	1000	"--------- embed fun Error_Pattern.check_for ------------------------";
walther@59844	1001	"--------- embed fun Error_Pattern.check_for ------------------------";
wneuper@59562	1002	reset_states ();
neuper@42428	1003	CalcTree
neuper@42428	1004	[(["functionTerm (x ^ 2 + sin (x ^ 4))", "differentiateFor x", "derivative f_f'"],
wneuper@59592	1005	("Isac_Knowledge", ["derivative_of","function"], ["diff","differentiate_on_R"]))];
neuper@42428	1006	Iterator 1;
neuper@42428	1007	moveActiveRoot 1;
wneuper@59248	1008	autoCalculate 1 CompleteCalcHead;
walther@59747	1009	autoCalculate 1 (Steps 1);
walther@59747	1010	autoCalculate 1 (Steps 1);(([1], Res), d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)))
walther@59747	1011	(autoCalculate 1 (Steps 1);([2], Res), d_d x (x ^^^ 2) + cos (x ^^^ 4) d_d x (x ^^^ 4)*)
neuper@42428	1012
wneuper@59562	1013	"~~~~~ fun appendFormula , args:"; val (cI, (ifo:cterm')) = (1, "d_d x (x ^ 2) + cos (4 * x ^ 3)");
wneuper@59562	1014	"~~~~~ fun appendFormula' , args:"; val (cI, (ifo: Rule.cterm')) = (cI, ifo);
wneuper@59562	1015	val cs = get_calc cI
wneuper@59562	1016	val pos = get_pos cI 1;
walther@59844	1017	(+)if pos = ([1], Res) then () else error "inform with (positive) Error_Pattern.check_for broken 1";
walther@59798	1018	val ("ok", cs' as (_, _, ptp)) = (case) Step.do_next pos cs (of);
walther@59798	1019	(case) Step_Solve.by_term ptp (encode ifo) (of); (ERROR WAS: "no derivation found")
walther@59797	1020	"~~~~~ fun Step_Solve.by_term , args:"; val (((next_)cs as (_, _, (pt, pos as (p, _))): Chead.calcstate'), istr)
wneuper@59562	1021	= (cs', (encode ifo));
wneuper@59562	1022	val ctxt = get_ctxt pt pos (see TODO.thy)
wneuper@59592	1023	val SOME f_in = (case) TermC.parse (Celem.assoc_thy "Isac_Knowledge") istr (of);
wneuper@59562	1024	val f_in = Thm.term_of f_in
wneuper@59562	1025	val pos_pred = lev_back' pos
wneuper@59562	1026	val f_pred = Ctree.get_curr_formula (pt, pos_pred);
wneuper@59562	1027	(if) f_pred = f_in; (else)
walther@59822	1028	val NONE = (case) In_Chead.cas_input f_in (of);
walther@59795	1029	(old )val (_, _, metID) = get_obj g_spec pt (par_pblobj pt p)
walther@59795	1030	(old)val {scr = prog, ...} = Specify.get_met metID
walther@59807	1031	(old)val istate = get_istate_LI pt pos
walther@59795	1032	(old)val ctxt = get_ctxt pt pos
walther@59795	1033	( old)
walther@59795	1034	val LI.Not_Derivable =
walther@59795	1035	(case) LI.locate_input_term (pt, pos) f_in (of);
wneuper@59562	1036	val pp = Ctree.par_pblobj pt p
wneuper@59562	1037	val (errpats, nrls, prog) = case Specify.get_met (Ctree.get_obj Ctree.g_metID pt pp) of
wneuper@59562	1038	{errpats, nrls, scr = Rule.Prog prog, ...} => (errpats, nrls, prog)
wneuper@59562	1039	\| _ => error "inform: uncovered case of get_met"
wneuper@59562	1040	;
wneuper@59562	1041	(+)if errpats2str errpats = "[(\"chain-rule-diff-both\",\n[\"d_d ?bdv (sin ?u) = cos (d_d ?bdv ?u)\",\"d_d ?bdv (cos ?u) = - sin (d_d ?bdv ?u)\",\"d_d ?bdv (?u ^^^ ?n) = ?n * ?u ^^^ (?n - 1)\",\"d_d ?bdv (LogExp.ln ?u) = 1 / ?u\",\"d_d ?bdv (LogExp.ln ?u) = 1 / d_d ?bdv ?u\"],\n[\"d_d ?bdv (sin ?u) = cos ?u * d_d ?bdv ?u\",\"d_d ?bdv (cos ?u) = - sin ?u * d_d ?bdv ?u\",\"d_d ?bdv (?u ^^^ ?n) = ?n * ?u ^^^ (?n - 1) * d_d ?bdv ?u\",\"d_d ?bdv (LogExp.ln ?u) = d_d ?bdv ?u / ?u\",\"d_d ?bdv (E_ ^^^ ?u) = E_ ^^^ ?u * d_d ?x ?u\"]]"
walther@59844	1042	(+)then () else error "inform with (positive) Error_Pattern.check_for broken 3";
neuper@42428	1043
walther@59807	1044	val {env, ...} = Ctree.get_istate_LI pt pos \|> Istate.the_pstate
wneuper@59562	1045	;
wneuper@59562	1046	(+)if term2str f_pred = "d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" andalso
wneuper@59562	1047	(+) term2str f_in = "d_d x (x ^^^ 2) + cos (4 * x ^^^ 3)"
walther@59844	1048	(+)then () else error "inform with (positive) Error_Pattern.check_for broken 2";
neuper@42428	1049
walther@59844	1050	val SOME "chain-rule-diff-both" = (case) Error_Pattern.Error_Pattern.check_for (f_pred, f_in) (prog, env) (errpats, nrls) (of);
neuper@42428	1051
neuper@42428	1052	"--- final check:";
walther@59798	1053	(+)val (_, _, ptp') = cs';
walther@59798	1054	case Step_Solve.by_term ptp' (encode ifo) of
neuper@42428	1055	("error pattern #chain-rule-diff-both#", calcstate') => ()
walther@59844	1056	\| _ => error "inform with (positive) Error_Pattern.check_for broken"
neuper@42428	1057
neuper@48895	1058
neuper@42430	1059	"--------- embed fun find_fillpatterns ---------------------------";
neuper@42430	1060	"--------- embed fun find_fillpatterns ---------------------------";
neuper@42430	1061	"--------- embed fun find_fillpatterns ---------------------------";
s1210629013@55445	1062	reset_states ();
neuper@42430	1063	CalcTree
neuper@42430	1064	[(["functionTerm (x ^ 2 + sin (x ^ 4))", "differentiateFor x", "derivative f_f'"],
wneuper@59592	1065	("Isac_Knowledge", ["derivative_of","function"], ["diff","differentiate_on_R"]))];
neuper@42430	1066	Iterator 1;
neuper@42430	1067	moveActiveRoot 1;
wneuper@59248	1068	autoCalculate 1 CompleteCalcHead;
walther@59747	1069	autoCalculate 1 (Steps 1);
walther@59747	1070	autoCalculate 1 (Steps 1);(([1], Res), d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)))
wneuper@59123	1071	appendFormula 1 "d_d x (x ^ 2) + cos (4 * x ^ 3)" (\|> Future.join);
neuper@42430	1072	(<CALCMESSAGE> error pattern #chain-rule-diff-both# </CALCMESSAGE>)
neuper@42432	1073	(*or
neuper@42432	1074	<CALCMESSAGE> no derivation found </CALCMESSAGE>*)
neuper@42430	1075
neuper@42433	1076	"~~~~~ fun findFillpatterns, args:"; val (cI, errpatID) = (1, "chain-rule-diff-both");
neuper@42430	1077	val ((pt, _), _) = get_calc cI
neuper@42430	1078	val pos = get_pos cI 1;
walther@59685	1079	"~~~~~ fun find_fillpatterns , args:"; val ((pt, pos as (p, _)), errpatID) = ((pt, pos), errpatID);
walther@59685	1080	val f_curr = Ctree.get_curr_formula (pt, pos);
walther@59685	1081	val pp = Ctree.par_pblobj pt p
walther@59685	1082	val (errpats, prog) = case Specify.get_met (Ctree.get_obj Ctree.g_metID pt pp) of
walther@59685	1083	{errpats, scr = Rule.Prog prog, ...} => (errpats, prog)
walther@59685	1084	\| _ => error "find_fillpatterns: uncovered case of get_met"
walther@59807	1085	val {env, ...} = Ctree.get_istate_LI pt pos \|> Istate.the_pstate
walther@59685	1086	val subst = Rtools.get_bdv_subst prog env
neuper@42430	1087	val errpatthms = errpats
neuper@42430	1088	\|> filter ((curry op = errpatID) o (#1: errpat -> errpatID))
neuper@42430	1089	\|> map (#3: errpat -> thm list)
neuper@42430	1090	\|> flat;
neuper@42430	1091
neuper@42430	1092	case map (get_fillpats subst f_curr errpatID) errpatthms \|> flat of
neuper@42433	1093	("fill-d_d-arg", tm, thm, subs_opt) :: _ => if term2str tm =
neuper@42432	1094	"d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)) =\nd_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x ?_dummy_1"
neuper@42432	1095	then () else error "find_fillpatterns changed 1a"
neuper@42432	1096	\| _ => error "find_fillpatterns changed 1b"
neuper@42430	1097
neuper@42430	1098	"~~~~~ fun get_fillpats, args:"; val (subst, form, errpatID, thm) =
neuper@42432	1099	(subst, f_curr, errpatID, hd (simulate beginning of "map") errpatthms);
neuper@42430	1100	val thmDeriv = Thm.get_name_hint thm
neuper@42430	1101	val (part, thyID) = thy_containing_thm thmDeriv
neuper@42430	1102	val theID = [part, thyID, "Theorems", thmID_of_derivation_name thmDeriv]
neuper@42430	1103	val Hthm {fillpats, ...} = get_the theID
neuper@42433	1104	val some = map (get_fillform subst (thm, form) errpatID) fillpats;
neuper@42430	1105
neuper@42430	1106	case some \|> filter is_some \|> map the of
neuper@42433	1107	("fill-d_d-arg", tm, thm, subsopt) :: _ => if term2str tm =
neuper@42432	1108	"d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)) =\nd_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x ?_dummy_1"
neuper@42432	1109	then () else error "find_fillpatterns changed 2a"
neuper@42432	1110	\| _ => error "find_fillpatterns changed 2b"
neuper@42430	1111
neuper@42433	1112	"~~~~~ fun get_fillform, args:";
neuper@42433	1113	val ((subs_opt, subst), (thm, form), errpatID, (fillpatID, pat, erpaID)) =
neuper@42433	1114	(subst, (thm, form), errpatID, hd (simulate beginning of "map") fillpats);
neuper@42433	1115	val (form', _, _, rewritten) =
walther@59852	1116	rew_sub (Isac()) 1 subst e_rew_ord Rule_Set.empty false [] (HOLogic.Trueprop $ pat) form;
neuper@42432	1117
neuper@42432	1118	if term2str form' = "d_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x ?_dummy_1" then ()
neuper@42432	1119	else error "find_fillpatterns changed 3";
neuper@42432	1120
neuper@42433	1121	"~~~~~ to findFillpatterns return val:"; val (fillpats) =
neuper@42433	1122	(map (get_fillpats (subs_opt, subst) f_curr errpatID) errpatthms \|> flat) (only from "hd errpatthms");
neuper@42432	1123
neuper@42450	1124	"vvv--- dropped this code WN120730";
neuper@42433	1125	val msg = "fill patterns " ^
neuper@42433	1126	((map ((apsnd term2str) o quad2pair) fillpats) \|> map pair2str_ \|> strs2str_);
neuper@42450	1127	msg =
neuper@42432	1128	"fill patterns #fill-d_d-arg#d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" ^
neuper@42432	1129	" =\nd_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x ?_dummy_1" ^
neuper@42432	1130	"#fill-both-args#d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" ^
neuper@42432	1131	" =\nd_d x (x ^^^ 2) + cos ?_dummy_2 * d_d x ?_dummy_3" ^
neuper@42432	1132	"#fill-d_d#d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" ^
neuper@42432	1133	" =\nd_d x (x ^^^ 2) + cos (x ^^^ 4) * ?_dummy_1 x (x ^^^ 4)" ^
neuper@42432	1134	"#fill-inner-deriv#d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" ^
neuper@42432	1135	" =\nd_d x (x ^^^ 2) + cos (x ^^^ 4) * ?_dummy_1" ^
neuper@42450	1136	"#fill-all#d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)) = d_d x (x ^^^ 2) + ?_dummy_1#";
neuper@42450	1137	"^^^--- dropped this code WN120730";
neuper@42450	1138
neuper@42450	1139	if (map #1 fillpats) =
neuper@42450	1140	["fill-d_d-arg", "fill-both-args", "fill-d_d", "fill-inner-deriv", "fill-all"]
neuper@42450	1141	then () else error "find_fillpatterns changed 4b";
wneuper@59252	1142	DEconstrCalcTree 1;
neuper@42450	1143
neuper@42437	1144	"--------- build fun is_exactly_equal, inputFillFormula ----------";
neuper@42437	1145	"--------- build fun is_exactly_equal, inputFillFormula ----------";
neuper@42437	1146	"--------- build fun is_exactly_equal, inputFillFormula ----------";
s1210629013@55445	1147	reset_states ();
neuper@42437	1148	CalcTree
neuper@42437	1149	[(["functionTerm (x ^ 2 + sin (x ^ 4))", "differentiateFor x", "derivative f_f'"],
wneuper@59592	1150	("Isac_Knowledge", ["derivative_of","function"], ["diff","differentiate_on_R"]))];
neuper@42437	1151	Iterator 1;
neuper@42437	1152	moveActiveRoot 1;
wneuper@59248	1153	autoCalculate 1 CompleteCalcHead;
walther@59747	1154	autoCalculate 1 (Steps 1);
walther@59747	1155	autoCalculate 1 (Steps 1);(([1], Res), d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)))
wneuper@59123	1156	appendFormula 1 "d_d x (x ^ 2) + cos (4 * x ^ 3)" (\|> Future.join); (<<<<<<<=========================)
neuper@42437	1157	(* the check for errpat is maximally liberal (whole term modulo "nrls" from "type met"),
neuper@42437	1158	would recognize "cos (4 * x ^ (4 - 1)) + 2 * x" as well.
neuper@42437	1159	results in <CALCMESSAGE> error pattern #chain-rule-diff-both# </CALCMESSAGE>
neuper@42437	1160	instead of <CALCMESSAGE> no derivation found </CALCMESSAGE> *)
neuper@42437	1161	val ((pt,pos), _) = get_calc 1;
neuper@42437	1162	val p = get_pos 1 1;
neuper@42437	1163	val (Form f, _, asms) = pt_extract (pt, p);
neuper@42437	1164
wneuper@59253	1165	if p = ([1], Res) andalso term2str f = "d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" then
wneuper@59497	1166	case get_obj g_tac pt (fst p) of Rewrite_Inst (["(''bdv'', x)"],
wneuper@59253	1167	("diff_sum", thm)) =>
wneuper@59253	1168	if (term2str o Thm.prop_of) thm = "d_d ?bdv (?u + ?v) = d_d ?bdv ?u + d_d ?bdv ?v" then ()
wneuper@59253	1169	else error "embed fun get_fillform changed 11"
wneuper@59253	1170	\| _ => error "embed fun get_fillform changed 12"
wneuper@59253	1171	else error "embed fun get_fillform changed 13";
neuper@42437	1172
neuper@42437	1173	findFillpatterns 1 "chain-rule-diff-both"; (<<<<<<<=================================)
neuper@42437	1174	(*<CALCMESSAGE> fill patterns #fill-d_d-arg#d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4)) =
neuper@42437	1175	d_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x ?_dummy_1#fill-both-args#...#...#... *)
neuper@42437	1176	val ((pt,pos),_) = get_calc 1;
neuper@42437	1177	val p = get_pos 1 1;
neuper@42437	1178
neuper@42437	1179	val (Form f, _, asms) = pt_extract (pt, p);
wneuper@59253	1180	if p = ([1], Res) andalso term2str f = "d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))" then
wneuper@59497	1181	case get_obj g_tac pt (fst p) of Rewrite_Inst (["(''bdv'', x)"],
wneuper@59253	1182	("diff_sum", thm)) =>
wneuper@59253	1183	if (term2str o Thm.prop_of) thm = "d_d ?bdv (?u + ?v) = d_d ?bdv ?u + d_d ?bdv ?v" then ()
wneuper@59253	1184	else error "embed fun get_fillform changed 21"
wneuper@59253	1185	\| _ => error "embed fun get_fillform changed 22"
wneuper@59253	1186	else error "embed fun get_fillform changed 23";
neuper@42437	1187
neuper@42437	1188	requestFillformula 1 ("chain-rule-diff-both", "fill-both-args");(<<<<<<<============)
neuper@42437	1189	(<AUTOCALC> ([1], Res) ([2], Res) ([2], Res) </AUTOCALC>)
neuper@42437	1190	val ((pt,pos),_) = get_calc 1;
neuper@42437	1191	val p = get_pos 1 1;
neuper@42437	1192	val (Form f, _, asms) = pt_extract (pt, p);
wneuper@59253	1193	if p = ([1], Res) andalso existpt [2] pt
wneuper@59253	1194	andalso term2str f = "d_d x (x ^^^ 2) + d_d x (sin (x ^^^ 4))"
wneuper@59497	1195	then case get_obj g_tac pt (fst p) of Rewrite_Inst (["(''bdv'', x)"],
wneuper@59253	1196	("diff_sum", thm)) =>
wneuper@59253	1197	if (term2str o Thm.prop_of) thm = "d_d ?bdv (?u + ?v) = d_d ?bdv ?u + d_d ?bdv ?v" then ()
wneuper@59253	1198	else error "embed fun get_fillform changed 31"
wneuper@59253	1199	\| _ => error "embed fun get_fillform changed 32"
wneuper@59253	1200	else error "embed fun get_fillform changed 33";
neuper@42437	1201
neuper@42437	1202	(* input a formula which exactly fills the gaps in a "fillformula"
neuper@42437	1203	presented to the learner immediately before by "requestFillformula (errpatID, fillpatID)":
neuper@42437	1204	errpatID: lhs of the respective thm = lhs of fillformula with fillpatID.
neuper@42437	1205	the respective thm is in the ctree ................
neuper@42437	1206	*)
neuper@42437	1207	"~~~~~ fun inputFillFormula, args:"; val (cI, ifo) =
neuper@42437	1208	(1, "d_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x (x ^^^ 4)");
neuper@42437	1209	val ((pt, _), _) = get_calc cI
neuper@42437	1210	val pos = get_pos cI 1;
neuper@42437	1211
neuper@42437	1212	"~~~~~ fun is_exactly_equal, args:"; val ((pt, pos as (p, p_)), istr) = ((pt, pos), ifo);
wneuper@59592	1213	val SOME ifo = parseNEW (assoc_thy "Isac_Knowledge" \|> thy2ctxt) istr
neuper@42437	1214	val p' = lev_on p;
neuper@42437	1215	val tac = get_obj g_tac pt p';
wneuper@59497	1216	val Rewrite_Inst ([bbb as "(''bdv'', x)"], ("diff_sin_chain", ttt)) = tac;
wneuper@59253	1217	if (term2str o Thm.prop_of) ttt = "d_d ?bdv (sin ?u) = cos ?u * d_d ?bdv ?u" then ()
wneuper@59252	1218	else error "inputFillFormula changed 10";
neuper@42437	1219	val Appl rew = applicable_in pos pt tac;
neuper@42437	1220	val Rewrite_Inst' (_, _, _, _, _, _, _, (res, _)) = rew;
neuper@42437	1221
neuper@42437	1222	"~~~~~ to inputFillFormula return val:"; val ("ok", tac) = ("ok", tac);
walther@59804	1223	val ("ok", (_, c, ptp as (_,p'))) = Step.by_tactic tac (pt, pos);
neuper@42437	1224	upd_calc cI (ptp, []);
neuper@42437	1225	upd_ipos cI 1 p';
neuper@42437	1226	autocalculateOK2xml cI pos (if null c then p' else last_elem c) p';
neuper@42437	1227
neuper@42437	1228	"~~~~~ final check:";
neuper@42437	1229	val ((pt, _),_) = get_calc 1;
neuper@42437	1230	val p = get_pos 1 1;
neuper@42437	1231	val (Form f, _, asms) = pt_extract (pt, p);
wneuper@59253	1232	if p = ([2], Res) andalso term2str f = "d_d x (x ^^^ 2) + cos (x ^^^ 4) * d_d x (x ^^^ 4)"
wneuper@59497	1233	then case get_obj g_tac pt (fst p) of Rewrite_Inst (["(''bdv'', x)"],
wneuper@59253	1234	("diff_sin_chain", thm)) =>
wneuper@59253	1235	if (term2str o Thm.prop_of) thm = "d_d ?bdv (sin ?u) = cos ?u * d_d ?bdv ?u" then ()
wneuper@59253	1236	else error "inputFillFormula changed 111"
wneuper@59253	1237	\| _ => error "inputFillFormula changed 112"
wneuper@59253	1238	else error "inputFillFormula changed 113";
neuper@42437	1239
wneuper@59262	1240	"--------- fun appl_adds -----------------------------------------";
wneuper@59262	1241	"--------- fun appl_adds -----------------------------------------";
wneuper@59262	1242	"--------- fun appl_adds -----------------------------------------";
wneuper@59262	1243	(* val (dI, oris, ppc, pbt, selct::ss) =
wneuper@59262	1244	(dI, pors, probl, ppc, map itms2fstr probl);
wneuper@59262	1245	...vvv
wneuper@59262	1246	*)
wneuper@59262	1247	(* val (dI, oris, ppc, pbt, (selct::ss))=
wneuper@59262	1248	(#1 (some_spec ospec spec), oris, []:itm list,
wneuper@59262	1249	((#ppc o get_pbt) (#2 (some_spec ospec spec))),(imodel2fstr imodel));
wneuper@59262	1250	val iii = appl_adds dI oris ppc pbt (selct::ss);
wneuper@59262	1251	tracing(itms2str_ thy iii);
wneuper@59262	1252
wneuper@59262	1253	val itm = appl_add' dI oris ppc pbt selct;
wneuper@59262	1254	val ppc = insert_ppc' itm ppc;
wneuper@59262	1255
wneuper@59262	1256	val _::selct::ss = (selct::ss);
wneuper@59262	1257	val itm = appl_add' dI oris ppc pbt selct;
wneuper@59262	1258	val ppc = insert_ppc' itm ppc;
wneuper@59262	1259
wneuper@59262	1260	val _::selct::ss = (selct::ss);
wneuper@59262	1261	val itm = appl_add' dI oris ppc pbt selct;
wneuper@59262	1262	val ppc = insert_ppc' itm ppc;
wneuper@59262	1263	tracing(itms2str_ thy ppc);
wneuper@59262	1264
wneuper@59262	1265	val _::selct::ss = (selct::ss);
wneuper@59262	1266	val itm = appl_add' dI oris ppc pbt selct;
wneuper@59262	1267	val ppc = insert_ppc' itm ppc;
wneuper@59262	1268	*)
wneuper@59264	1269	"--------- fun concat_deriv --------------------------------------";
wneuper@59264	1270	"--------- fun concat_deriv --------------------------------------";
wneuper@59264	1271	"--------- fun concat_deriv --------------------------------------";
wneuper@59264	1272	(*
wneuper@59264	1273	val ({rew_ord, erls, rules,...}, fo, ifo) =
walther@59852	1274	(Rule_Set.rep Test_simplify, str2term "x+1+ -12=0", str2term "-21+(x+1)=0");
wneuper@59264	1275	(tracing o trtas2str) fod';
wneuper@59264	1276	> ["
wneuper@59264	1277	(x + 1 + -1 * 2 = 0, Thm ("radd_commute","?m + ?n = ?n + ?m"), (-1 * 2 + (x + 1) = 0, []))","
wneuper@59264	1278	(-1 * 2 + (x + 1) = 0, Thm ("radd_commute","?m + ?n = ?n + ?m"), (-1 * 2 + (1 + x) = 0, []))","
wneuper@59264	1279	(-1 * 2 + (1 + x) = 0, Thm ("radd_left_commute","?x + (?y + ?z) = ?y + (?x + ?z)"), (1 + (-1 * 2 + x) = 0, []))","
wneuper@59264	1280	(1 + (-1 * 2 + x) = 0, Thm ("#mult_Float ((~1,0), (0,0)) __ ((2,0), (0,0))","-1 * 2 = -2"), (1 + (-2 + x) = 0, []))"]
wneuper@59264	1281	val it = () : unit
wneuper@59264	1282	(tracing o trtas2str) (map rev_deriv' rifod');
wneuper@59264	1283	> ["
wneuper@59264	1284	(1 + (-2 + x) = 0, Thm ("sym_#mult_Float ((~2,0), (0,0)) __ ((1,0), (0,0))","-2 = -2 * 1"), (1 + (-2 * 1 + x) = 0, []))","
wneuper@59264	1285	(1 + (-2 * 1 + x) = 0, Thm ("sym_radd_left_commute","?y + (?x + ?z) = ?x + (?y + ?z)"), (-2 * 1 + (1 + x) = 0, []))","
wneuper@59264	1286	(-2 * 1 + (1 + x) = 0, Thm ("sym_radd_commute","?n + ?m = ?m + ?n"), (-2 * 1 + (x + 1) = 0, []))"]
wneuper@59264	1287	val it = () : unit
wneuper@59264	1288	*)
wneuper@59264	1289	"--------- handle an input formula -------------------------------";
wneuper@59264	1290	"--------- handle an input formula -------------------------------";
wneuper@59264	1291	"--------- handle an input formula -------------------------------";
wneuper@59264	1292	(*
wneuper@59264	1293	Untersuchung zur Formeleingabe (appendFormula, replaceFormla) zu einer Anregung von Alan Krempler:
wneuper@59264	1294	Welche RICHTIGEN Formeln koennen NICHT abgeleitet werden,
wneuper@59264	1295	wenn Abteilungen nur auf gleichem Level gesucht werden ?
wneuper@59264	1296	WN.040216
wneuper@59264	1297
wneuper@59264	1298	Beispiele zum Equationsolver von Richard Lang aus /src/sml/kbtest/rlang.sml
wneuper@59264	1299
wneuper@59264	1300	------------------------------------------------------------------------------
wneuper@59264	1301	"Schalk I s.87 Bsp 52a ((5*x)/(x - 2) - x/(x+2)=4)";
wneuper@59264	1302	------------------------------------------------------------------------------
wneuper@59264	1303	1. "5 * x / (x - 2) - x / (x + 2) = 4"
wneuper@59264	1304	...
wneuper@59367	1305	4. "12 * x + 4 * x ^^^ 2 = 4 * (-4 + x ^^^ 2)",Subproblem["normalise", "poly"..
wneuper@59264	1306	...
wneuper@59264	1307	4.3. "16 + 12 * x = 0", Subproblem["degree_1", "polynomial", "univariate"..
wneuper@59264	1308	...
wneuper@59264	1309	4.3.3. "[x = -4 / 3]")), Check_elementwise "Assumptions"
wneuper@59264	1310	...
wneuper@59264	1311	"[x = -4 / 3]"
wneuper@59264	1312	------------------------------------------------------------------------------
wneuper@59264	1313	(1)..(6): keine 'richtige' Eingabe kann abgeleitet werden, die einen Summanden auf die andere Seite verschiebt [Ableitung ware in 4.3.n]
wneuper@59264	1314
wneuper@59264	1315	(4.1)..(4.3): keine 'richtige' Eingabe kann abgeleitet werden, die einen Summanden auf die andere Seite verschiebt [Ableitung ware in 4.3.n]
wneuper@59264	1316	------------------------------------------------------------------------------
wneuper@59264	1317
wneuper@59264	1318
wneuper@59264	1319	------------------------------------------------------------------------------
wneuper@59264	1320	"Schalk I s.87 Bsp 55b (x/(x^^^2 - 6x+9) - 1/(x^^^2 - 3x) =1/x)";
wneuper@59264	1321	------------------------------------------------------------------------------
wneuper@59264	1322	1. "x / (x ^^^ 2 - 6 * x + 9) - 1 / (x ^^^ 2 - 3 * x) = 1 / x"
wneuper@59264	1323	...
wneuper@59264	1324	4. "(3 + (-1 * x + x ^^^ 2)) * x = 1 * (9 * x + (x ^^^ 3 + -6 * x ^^^ 2))"
wneuper@59367	1325	Subproblem["normalise", "polynomial", "univariate"..
wneuper@59264	1326	...
wneuper@59264	1327	4.4. "-6 * x + 5 * x ^^^ 2 = 0", Subproblem["bdv_only", "degree_2", "poly"..
wneuper@59264	1328	...
wneuper@59264	1329	4.4.4. "[x = 0, x = 6 / 5]", Check_elementwise "Assumptions"
wneuper@59264	1330	4.4.5. "[x = 0, x = 6 / 5]"
wneuper@59264	1331	...
wneuper@59264	1332	5. "[x = 0, x = 6 / 5]", Check_elementwise "Assumptions"
wneuper@59264	1333	"[x = 6 / 5]"
wneuper@59264	1334	------------------------------------------------------------------------------
wneuper@59264	1335	(1)..(4): keine 'richtige' Eingabe kann abgeleitet werden, die einen Summanden auf die andere Seite schiebt [Ableitung waere in 4.4.x]
wneuper@59264	1336
wneuper@59264	1337	(4.1)..(4.4.5): keine 'richtige' Eingabe kann abgeleitet werden, die dem Ergebnis "[x = 6 / 5]" aequivalent ist [Ableitung waere in 5.]
wneuper@59264	1338	------------------------------------------------------------------------------
wneuper@59264	1339
wneuper@59264	1340
wneuper@59264	1341	------------------------------------------------------------------------------
wneuper@59264	1342	"Schalk II s.56 Bsp 73b (sqrt(x+1)+sqrt(4x+4)=sqrt(9x+9))";
wneuper@59264	1343	------------------------------------------------------------------------------
wneuper@59264	1344	1. "sqrt (x + 1) + sqrt (4 * x + 4) = sqrt (9 * x + 9)"
wneuper@59264	1345	...
wneuper@59264	1346	6. "13 + 13 * x + -2 * sqrt ((4 + 4 * x) * (9 + 9 * x)) = 1 + x"
wneuper@59489	1347	Subproblem["sq", "rootX", "univariate", "equation"]
wneuper@59264	1348	...
wneuper@59264	1349	6.6. "144 + 288 * x + 144 * x ^^^ 2 = 144 + x ^^^ 2 + 288 * x + 143 * x ^^^ 2"
wneuper@59367	1350	Subproblem["normalise", "polynomial", "univariate", "equation"]
wneuper@59264	1351	...
wneuper@59264	1352	6.6.3 "0 = 0" Subproblem["degree_0", "polynomial", "univariate", "equation"]
wneuper@59264	1353	... Or_to_List
wneuper@59264	1354	6.6.3.2 "UniversalList"
wneuper@59264	1355	------------------------------------------------------------------------------
wneuper@59264	1356	(1)..(6): keine 'richtige' Eingabe kann abgeleitet werden, die eine der Wurzeln auf die andere Seite verschieb [Ableitung ware in 6.6.n]
wneuper@59264	1357
wneuper@59264	1358	(6.1)..(6.3): keine 'richtige' Eingabe kann abgeleitet werden, die einen Summanden auf die andere Seite verschiebt [Ableitung ware in 6.6.n]
wneuper@59264	1359	------------------------------------------------------------------------------
wneuper@59264	1360	*)
wneuper@59264	1361	(sh. comments auf 498)
wneuper@59264	1362	"--------- fun dropwhile' ----------------------------------------";
wneuper@59264	1363	"--------- fun dropwhile' ----------------------------------------";
wneuper@59264	1364	"--------- fun dropwhile' ----------------------------------------";
wneuper@59264	1365	(*
wneuper@59264	1366	fun equal a b = a=b;
wneuper@59264	1367	val foder = [0,1,2,3,4,5]; val ifoder = [11,12,3,4,5];
wneuper@59264	1368	val r_foder = rev foder; val r_ifoder = rev ifoder;
wneuper@59264	1369	dropwhile' equal r_foder r_ifoder;
wneuper@59264	1370	> vval it = ([0, 1, 2, 3], [3, 12, 11]) : int list * int list
wneuper@59264	1371
wneuper@59264	1372	val foder = [3,4,5]; val ifoder = [11,12,3,4,5];
wneuper@59264	1373	val r_foder = rev foder; val r_ifoder = rev ifoder;
wneuper@59264	1374	dropwhile' equal r_foder r_ifoder;
wneuper@59264	1375	> val it = ([3], [3, 12, 11]) : int list * int list
wneuper@59264	1376
wneuper@59264	1377	val foder = [5]; val ifoder = [11,12,3,4,5];
wneuper@59264	1378	val r_foder = rev foder; val r_ifoder = rev ifoder;
wneuper@59264	1379	dropwhile' equal r_foder r_ifoder;
wneuper@59264	1380	> val it = ([5], [5, 4, 3, 12, 11]) : int list * int list
wneuper@59264	1381
wneuper@59264	1382	val foder = [10,11,12,13,14,15]; val ifoder = [11,12,3,4,5];
wneuper@59264	1383	val r_foder = rev foder; val r_ifoder = rev ifoder;
wneuper@59264	1384	dropwhile' equal r_foder r_ifoder;
wneuper@59264	1385	> *** dropwhile': did not start with equal elements*)

author	Walther Neuper <walther.neuper@jku.at>
	Thu, 09 Apr 2020 12:03:14 +0200
changeset 59860	64eecda7a8fb
parent 59852	test/Tools/isac/Interpret/error-pattern.sml@ea7e6679080e
child 59861	65ec9f679c3f
permissions	-rw-r--r--