(* test for sml/ME/mathengine.sml

   authors: Walther Neuper 2000, 2006

   (c) due to copyright terms

*)

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

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

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

"----------- debugging setContext..pbl_ -----------------";

"----------- tryrefine ----------------------------------";

"----------- fun step -----------------------------------";

"----------- fun autocalc -------------------------------";

"----------- fun autoCalculate --------------------------";

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

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

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

(*=== inhibit exn ?=============================================================

"----------- debugging setContext..pbl_ -----------------";

"----------- debugging setContext..pbl_ -----------------";

"----------- debugging setContext..pbl_ -----------------";

states:=[];

CalcTree

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

  ("Test", 

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

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

Iterator 1;

moveActiveRoot 1; modelProblem 1;

val pos as (p,_) = ([],Pbl);

val guh = "pbl_equ_univ";

checkContext 1 pos guh;

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

val pp = par_pblobj pt p;

val keID = guh2kestoreID guh;

case context_pbl keID pt pp of (true,["univariate", "equation"],_,_,_)=>()

| _ => error "mathengine.sml: context_pbl .. pbl_equ_univ checked";

case get_obj g_spec pt p of (_, ["e_pblID"], _) => ()

| _ => error "mathengine.sml: context_pbl .. pbl still empty";

setContext 1 pos guh;

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

case get_obj g_spec pt p of (_, ["univariate", "equation"], _) => ()

| _ => error "mathengine.sml: context_pbl .. pbl set";

setContext 1 pos "met_eq_lin";

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

case get_obj g_spec pt p of (_,  _, ["LinEq", "solve_lineq_equation"]) => ()

| _ => error "mathengine.sml: context_pbl .. pbl set";

"----------- tryrefine ----------------------------------";

"----------- tryrefine ----------------------------------";

"----------- tryrefine ----------------------------------";

states:=[];

CalcTree [(["equality (x/(x^2 - 6*x+9) - 1/(x^2 - 3*x) =1/x)", 

	    "solveFor x", "solutions L"],

	   ("RatEq",["univariate","equation"],["no_met"]))];

Iterator 1;

moveActiveRoot 1; autoCalculate 1 CompleteCalc;

refineProblem 1 ([1],Res) "pbl_equ_univ" 

(*gives "pbl_equ_univ_rat" correct*);

refineProblem 1 ([1],Res) (pblID2guh ["linear","univariate","equation"])

(*gives "pbl_equ_univ_lin" incorrect*);

===== inhibit exn ?===========================================================*)

"----------- fun step -----------------------------------";

"----------- fun step -----------------------------------";

"----------- fun step -----------------------------------";

val p = e_pos'; val c = []; 

val (p,_,f,nxt,_,pt) = 

    CalcTreeTEST 

        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"], 

          ("Integrate", ["integrate","function"], ["diff","integration"]))];

"----- step 1: returns tac = Model_Problem ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 2: returns tac =  ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 3: returns tac =  ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 4: returns tac =  ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 5: returns tac =  ---";

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

2002 stops here as well: TODO review actual arguments:

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 6: returns tac =  ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 7: returns tac =  ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

"----- step 8: returns tac =  ---";

val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);

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

"----------- fun autocalc -------------------------------";

"----------- fun autocalc -------------------------------";

"----------- fun autocalc -------------------------------";

val p = e_pos'; val c = []; 

val (p,_,f,nxt,_,pt) = 

    CalcTreeTEST 

        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"], 

          ("Integrate",["integrate","function"], ["diff","integration"]))];

(*-----since Model_Problem + complete_mod_ in case cas of SOME-----*

modeling is skipped FIXME 

 *-----since Model_Problem + complete_mod_ in case cas of SOME-----*)

tracing "----- step 1 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 2 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 3 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 4 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 5 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 6 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 7 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

tracing "----- step 8 ---";

val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;

if str = "end-of-calculation" then ()

else error "mathengine.sml -- fun autocalc -- end";

"----------- fun autoCalculate -----------------------------------";

"----------- fun autoCalculate -----------------------------------";

"----------- fun autoCalculate -----------------------------------";

states := [];

CalcTree (*ATTENTION: encode_fmz ... unlike CalcTreeTEST*)

    [(["functionTerm (x^2 + 1)", "integrateBy x", "antiDerivative FF"],

      ("Integrate", ["integrate", "function"], ["diff", "integration"]))];

Iterator 1;

moveActiveRoot 1;

(*-----since Model_Problem + complete_mod_ in case cas of SOME-----*

modeling is skipped FIXME 

see test/../interface -- solve_linear as rootpbl FE -- for OTHER expl:

 setNextTactic 1 (Add_Given "equality (1 + -1 * 2 + x = 0)");

 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1); (*equality added*);

 fetchProposedTactic 1;

 setNextTactic 1 (Add_Given "solveFor x");

 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);

 fetchProposedTactic 1;

 setNextTactic 1 (Add_Find "solutions L");

 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);

 fetchProposedTactic 1;

 setNextTactic 1 (Specify_Theory "Test");

 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);

*-----since Model_Problem + complete_mod_ in case cas of SOME-----*)

autoCalculate 1 (Step 1); 

"----- step 1 ---";

autoCalculate 1 (Step 1);

"----- step 2 ---";

autoCalculate 1 (Step 1);

"----- step 3 ---";

autoCalculate 1 (Step 1);

"----- step 4 ---";

autoCalculate 1 (Step 1);

"----- step 5 ---";

autoCalculate 1 (Step 1);

"----- step 6 ---";

autoCalculate 1 (Step 1);

"----- step 7 ---";

autoCalculate 1 (Step 1);

"----- step 8 ---";

autoCalculate 1 (Step 1);

val (ptp as (_, p), _) = get_calc 1;

val (Form t,_,_) = pt_extract ptp;

if term2str t = "c + x + 1 / 3 * x ^^^ 3" andalso p = ([], Res) then ()

else error "mathengine.sml -- fun autoCalculate -- end";

(*-.-.-.-.-.-isolate response.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.

-.-.-.-.-.-.-isolate response.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.*)