test/Tools/isac/Interpret/mathengine.sml
author Walther Neuper <neuper@ist.tugraz.at>
Sun, 04 May 2014 14:02:23 +0200
changeset 55412 632d2ecab96f
parent 55279 130688f277ba
child 55415 f8d44cb86330
permissions -rw-r--r--
error (from <NEXT> on new Worksheet) related to thehier
neuper@42279
     1
(* Title:  tests for Interpret/mathengine.sml
neuper@42279
     2
   Author: Walther Neuper 2000, 2006
neuper@37906
     3
   (c) due to copyright terms
neuper@37906
     4
*)
neuper@38065
     5
"--------------------------------------------------------";
neuper@38065
     6
"table of contents --------------------------------------";
neuper@38065
     7
"--------------------------------------------------------";
neuper@41940
     8
"----------- change to parse ctxt -----------------------";
neuper@38065
     9
"----------- debugging setContext..pbl_ -----------------";
neuper@38065
    10
"----------- tryrefine ----------------------------------";
neuper@55412
    11
"---------- mini-subpbl isac-java 3992e3bd1948: initial <NEXT> ---------------";
neuper@41982
    12
"----------- fun step: Apply_Method without init_form ---";
neuper@38065
    13
"----------- fun step -----------------------------------";
neuper@38065
    14
"----------- fun autocalc -------------------------------";
neuper@38065
    15
"----------- fun autoCalculate --------------------------";
neuper@42067
    16
"----------- fun autoCalculate..CompleteCalc ------------";
t@42225
    17
"----------- search for Or_to_List ----------------------";
neuper@42279
    18
"----------- check thy in CalcTreeTEST ------------------";
neuper@38065
    19
"--------------------------------------------------------";
neuper@38065
    20
"--------------------------------------------------------";
neuper@38065
    21
"--------------------------------------------------------";
neuper@37906
    22
neuper@41940
    23
"----------- change to parse ctxt -----------------------";
neuper@41940
    24
"----------- change to parse ctxt -----------------------";
neuper@41940
    25
"----------- change to parse ctxt -----------------------";
neuper@41940
    26
"===== start calculation: from problem description (fmz) to origin";
neuper@41940
    27
val fmz = ["realTestGiven (((1+2)*4/3)^^^2)","realTestFind s"];
neuper@41940
    28
val (thyID, pblID, metID) =
neuper@41940
    29
  ("Test", ["calculate", "test"], ["Test", "test_calculate"]);
neuper@41940
    30
val (p,_,_,nxt,_,pt) = CalcTreeTEST [(fmz, (thyID, pblID, metID))];
neuper@41940
    31
"----- ";
neuper@41940
    32
(* call sequence: CalcTreeTEST 
neuper@41940
    33
                > nxt_specify_init_calc 
neuper@41940
    34
                > prep_ori
neuper@41940
    35
*)
neuper@41940
    36
val (thy, pbt) = (Thy_Info.get_theory thyID, (#ppc o get_pbt) pblID);
neuper@41940
    37
"----- in  prep_ori";
neuper@48761
    38
val ctxt = Proof_Context.init_global thy;
neuper@41940
    39
bonzai@41952
    40
val ctopts = map (parseNEW ctxt) fmz;
bonzai@41952
    41
val dts = map (split_dts o the) ctopts;
neuper@41940
    42
(*split_dts:
neuper@41940
    43
(term * term list) list
neuper@41940
    44
        formulas: e.g. ((1+2)*4/3)^^^2
neuper@41940
    45
 description: e.g. realTestGiven
neuper@41940
    46
*)
neuper@41940
    47
 prep_ori;
neuper@41940
    48
(* FROM
neuper@41940
    49
val it = fn:
neuper@41940
    50
   string list -> theory -> (string * (term * 'a)) list -> ori list
neuper@41940
    51
TO
neuper@41940
    52
val it = fn:
neuper@41940
    53
   string list -> theory -> (string * (term * 'a)) list -> (ori list * ctxt)
neuper@41940
    54
AND
neuper@41940
    55
FROM val oris = prep_ori...
neuper@41940
    56
TO   val (oris, _) = prep_ori...
neuper@41940
    57
*)
neuper@41940
    58
"----- insert ctxt in ptree";
neuper@41940
    59
(* datatype ppobj
neuper@41940
    60
FROM loc   : istate option * istate option,
neuper@41940
    61
TO   loc   : (istate * ctxt) option * (istate * ctxt) option,
neuper@41940
    62
e.g.
neuper@41940
    63
FROM val pt = Nd (PblObj
neuper@41940
    64
       {.., loc = (SOME (ScrState ([], [], NONE, Const ("empty", "'a"), Sundef, false)),
neuper@41940
    65
          NONE),
neuper@41940
    66
TO   val pt = Nd (PblObj
neuper@41940
    67
       {.., loc = 
neuper@41940
    68
        ((SOME (ScrState ([], [], NONE, Const ("empty", "'a"), Sundef, false)), ctxt),
neuper@41940
    69
          NONE),
neuper@41940
    70
*)
neuper@41940
    71
neuper@41940
    72
"===== interactive specification: from origin to specification (probl)";
neuper@41940
    73
val (p,_,_,nxt,_,pt) = me nxt p [1] pt; 
neuper@41940
    74
  (*nxt =("Add_Given", Model_Problem)*)
neuper@41940
    75
val (p,_,_,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
    76
  (*nxt = ("Add_Find",Add_Find "realTestFind s") : string * tac*)
neuper@41940
    77
"----- ";
neuper@41940
    78
(* call sequence: me Model_Problem 
neuper@41940
    79
                > me ("Add_Given", Add_Given "realTestGiven (((1 + 2) * 4 / 3) ^^^ 2)")
neuper@41940
    80
                > locatetac tac
neuper@41940
    81
                > loc_specify_
neuper@41940
    82
                > specify          GET ctxt (stored in ctree)
neuper@41940
    83
                > specify_additem
neuper@41940
    84
                > appl_add
neuper@41940
    85
neuper@41940
    86
*)
neuper@41940
    87
"----- in appl_add";
neuper@41940
    88
(* FROM appl_add thy
neuper@41940
    89
   TO   appl_add ctxt
neuper@41940
    90
   FROM parse thy str
neuper@41940
    91
   TO   parseNEW ctxt str
neuper@41940
    92
*)
neuper@41940
    93
val (p,_,_,nxt,_,pt) = me nxt p [1] pt; 
neuper@41940
    94
  (*nxt = ("Specify_Theory",Specify_Theory "Test") : string * tac*)
neuper@41940
    95
val (p,_,_,nxt,_,pt) = me nxt p [1] pt; 
neuper@41940
    96
  (*nxt = ("Specify_Problem",Specify_Problem ["calculate","test"])*)
neuper@41940
    97
neuper@41940
    98
"===== end specify: from specification (probl) to guard (method)";
neuper@41940
    99
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   100
  (*nxt = ("Specify_Method",Specify_Method ("Test","test_calculate"))*)
neuper@41940
   101
neuper@41940
   102
"===== start interpretation: from guard to environment";
neuper@41940
   103
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   104
  (*nxt = ("Apply_Method",Apply_Method ("Test","test_calculate"))*)
neuper@41940
   105
"----- ";
neuper@41940
   106
(* call sequence: me ("Apply_Method",...
neuper@41940
   107
                > locatetac
neuper@41940
   108
                > loc_solve_
neuper@41940
   109
                > solve ("Apply_Method",...
neuper@41940
   110
*)
neuper@41940
   111
val ((_,tac), ptp) = (nxt, (pt, p));
neuper@41940
   112
locatetac tac (pt,p);
neuper@41940
   113
  val (mI, m) = mk_tac'_ tac;
neuper@41940
   114
  val Appl m = applicable_in p pt m;
neuper@41940
   115
  member op = specsteps mI;
neuper@41940
   116
  loc_solve_ (mI,m) ptp;
neuper@41940
   117
    val (m, (pt, pos)) = ((mI,m), ptp);
neuper@41940
   118
    solve m (pt, pos);
bonzai@41949
   119
      val ((_, m as Apply_Method' (mI, _, _, _)), (pt, (pos as (p,_)))) = 
neuper@41940
   120
        (m, (pt, pos));
neuper@41940
   121
      val {srls,...} = get_met mI;
neuper@41940
   122
      val PblObj{meth=itms,...} = get_obj I pt p;
neuper@41940
   123
      val thy' = get_obj g_domID pt p;
neuper@41940
   124
      val thy = assoc_thy thy';
bonzai@41949
   125
      val (is as ScrState (env,_,_,_,_,_), ctxt, sc) = init_scrstate thy itms mI;
neuper@41940
   126
neuper@41940
   127
"----- go on in the calculation";
neuper@41940
   128
val (p,_,f,nxt,_,pt) = me nxt pos [1] pt;
neuper@41940
   129
  (*nxt = ("Calculate",Calculate "PLUS")*)
neuper@41940
   130
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   131
  (*nxt = ("Calculate",Calculate "TIMES")*)
neuper@41940
   132
neuper@41940
   133
"===== input a formula to be derived from previous istate";
neuper@41940
   134
"----- appendFormula TODO: first repair error";
neuper@41940
   135
  val cs = ((pt,p),[]);
neuper@41940
   136
  val ("ok", cs' as (_,_,(pt,p))) = step p cs;
neuper@41943
   137
  val ifo = (* encode "4^^^2" \<longrightarrow>*) "4^2";
neuper@41940
   138
(*
neuper@41943
   139
  val ("no derivation found", (_,_,(pt, p))) = inform cs' ((*encode*) ifo);
neuper@41940
   140
  here ^^^^^^^^^^^^^^^^^^^^^ should be "ok"
neuper@41940
   141
*)
neuper@41940
   142
neuper@41940
   143
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   144
(*nxt = ("Calculate",Calculate "DIVIDE")*)
neuper@41940
   145
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   146
(*nxt = ("Calculate",Calculate "POWER")*)
neuper@41940
   147
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   148
(*nxt = ("Check_Postcond",Check_Postcond ["calculate","test"])*)
neuper@41940
   149
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   150
(*nxt = ("End_Proof'",End_Proof')*)
neuper@41940
   151
if f = Form' (FormKF (~1,EdUndef,0,Nundef,"16")) then ()
neuper@41940
   152
else error "calculate.sml: script test_calculate changed behaviour";
neuper@41940
   153
neuper@41940
   154
"===== tactic Subproblem: from environment to origin";
neuper@41940
   155
"----- TODO";
neuper@41940
   156
neuper@37906
   157
neuper@38065
   158
"----------- debugging setContext..pbl_ -----------------";
neuper@38065
   159
"----------- debugging setContext..pbl_ -----------------";
neuper@38065
   160
"----------- debugging setContext..pbl_ -----------------";
neuper@37906
   161
states:=[];
neuper@41970
   162
CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@41970
   163
  ("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906
   164
   ["Test","squ-equ-test-subpbl1"]))];
neuper@37906
   165
Iterator 1;
neuper@37906
   166
moveActiveRoot 1; modelProblem 1;
neuper@37906
   167
neuper@37906
   168
val pos as (p,_) = ([],Pbl);
neuper@37906
   169
val guh = "pbl_equ_univ";
neuper@37906
   170
checkContext 1 pos guh;
neuper@37906
   171
val ((pt,_),_) = get_calc 1;
neuper@37906
   172
val pp = par_pblobj pt p;
neuper@37906
   173
val keID = guh2kestoreID guh;
neuper@37906
   174
case context_pbl keID pt pp of (true,["univariate", "equation"],_,_,_)=>()
neuper@38031
   175
| _ => error "mathengine.sml: context_pbl .. pbl_equ_univ checked";
neuper@37906
   176
neuper@37906
   177
case get_obj g_spec pt p of (_, ["e_pblID"], _) => ()
neuper@38031
   178
| _ => error "mathengine.sml: context_pbl .. pbl still empty";
neuper@37906
   179
setContext 1 pos guh;
neuper@37906
   180
val ((pt,_),_) = get_calc 1;
neuper@37906
   181
case get_obj g_spec pt p of (_, ["univariate", "equation"], _) => ()
neuper@38031
   182
| _ => error "mathengine.sml: context_pbl .. pbl set";
neuper@37906
   183
neuper@37906
   184
neuper@37906
   185
setContext 1 pos "met_eq_lin";
neuper@37906
   186
val ((pt,_),_) = get_calc 1;
neuper@37906
   187
case get_obj g_spec pt p of (_,  _, ["LinEq", "solve_lineq_equation"]) => ()
neuper@38031
   188
| _ => error "mathengine.sml: context_pbl .. pbl set";
neuper@37906
   189
neuper@37906
   190
neuper@38065
   191
"----------- tryrefine ----------------------------------";
neuper@38065
   192
"----------- tryrefine ----------------------------------";
neuper@38065
   193
"----------- tryrefine ----------------------------------";
neuper@37906
   194
states:=[];
neuper@37906
   195
CalcTree [(["equality (x/(x^2 - 6*x+9) - 1/(x^2 - 3*x) =1/x)", 
neuper@37906
   196
	    "solveFor x", "solutions L"],
neuper@38058
   197
	   ("RatEq",["univariate","equation"],["no_met"]))];
neuper@37906
   198
Iterator 1;
neuper@37906
   199
moveActiveRoot 1; autoCalculate 1 CompleteCalc;
neuper@37906
   200
neuper@37906
   201
refineProblem 1 ([1],Res) "pbl_equ_univ" 
neuper@37906
   202
(*gives "pbl_equ_univ_rat" correct*);
neuper@37906
   203
neuper@55279
   204
refineProblem 1 ([1],Res) (pblID2guh ["LINEAR","univariate","equation"])
neuper@38065
   205
(*gives "pbl_equ_univ_lin" incorrect*);
neuper@38065
   206
neuper@55412
   207
"---------- mini-subpbl isac-java 3992e3bd1948: initial <NEXT> ---------------";
neuper@55412
   208
"---------- mini-subpbl isac-java 3992e3bd1948: initial <NEXT> ---------------";
neuper@55412
   209
"---------- mini-subpbl isac-java 3992e3bd1948: initial <NEXT> ---------------";
neuper@55412
   210
states:=[];
neuper@55412
   211
(*this is the exact sequence of input provided by isac-java 3992e3bd1948;
neuper@55412
   212
  Surrounding ML { * ... * } are omitted.
neuper@55412
   213
  The multiple calls suggest to replicate the CalcTree in the Dialogue.
neuper@55412
   214
neuper@55412
   215
  Reproduce error: On the newly opened Worksheet click <NEXT>.
neuper@55412
   216
  The error indicates a problem with the hierarchy of theories. 
neuper@55412
   217
*)
neuper@55412
   218
CalcTree [(["equality (x+1=(2::real))", "solveFor x", "solutions L"],
neuper@55412
   219
  ("Test",["sqroot-test","univariate","equation","test"],["Test","squ-equ-test-subpbl1"]))]; 
neuper@55412
   220
Iterator 1; 
neuper@55412
   221
moveActiveRoot 1; 
neuper@55412
   222
getFormulaeFromTo 1 ([],Pbl) ([],Pbl) 0 false; 
neuper@55412
   223
getFormulaeFromTo 1 ([],Pbl) ([],Pbl) 0 false; 
neuper@55412
   224
(*completeCalcHead*)autoCalculate 1 CompleteCalcHead; 
neuper@55412
   225
(*completeCalcHead*)getActiveFormula 1 ; 
neuper@55412
   226
(*completeCalcHead*)refFormula 1 ([],Met); 
neuper@55412
   227
refFormula 1 ([],Pbl); 
neuper@55412
   228
(*<REFFORMULA>  <CALCID> 1 </CALCID>  <CALCHEAD status = "correct"> *)
neuper@55412
   229
(*fetchProposedTactic 1; 
neuper@55412
   230
<SYSERROR>  <CALCID> 1 </CALCID>  <ERROR> error in kernel </ERROR></SYSERROR>*)
neuper@55412
   231
"~~~~~ fun fetchProposedTactic, args:"; val (cI) = (1);
neuper@55412
   232
val ("ok", (tacis, _, _)) = step (get_pos cI 1) (get_calc cI)
neuper@55412
   233
val _= upd_tacis cI tacis
neuper@55412
   234
	       val (tac,_,_) = last_elem tacis;
neuper@55412
   235
(*fetchproposedtacticOK2xml cI tac (fetchErrorpatterns tac);
neuper@55412
   236
**************************************************************************
neuper@55412
   237
ERROR: app_py: not found: ["IsacKnowledge","Tools","Rulesets","list_rls"]
neuper@55412
   238
**************************************************************************)
neuper@38065
   239
neuper@41982
   240
"----------- fun step: Apply_Method without init_form ---";
neuper@41982
   241
"----------- fun step: Apply_Method without init_form ---";
neuper@41982
   242
"----------- fun step: Apply_Method without init_form ---";
neuper@41982
   243
val fmz = ["equality (x+1=(2::real))", "solveFor x","solutions L"];
neuper@41982
   244
val (dI',pI',mI') =
neuper@41982
   245
  ("Test", ["sqroot-test","univariate","equation","test"],
neuper@41982
   246
   ["Test","squ-equ-test-subpbl1"]);
neuper@41982
   247
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
neuper@41982
   248
"~~~~~ fun me, args:"; val (_,tac) = nxt;
neuper@41982
   249
val ("ok", (_, _, (pt, p))) = locatetac tac (pt,p);
neuper@41982
   250
"~~~~~ fun step, args:"; val (ip as (_,p_), ptp as (pt,p), tacis) = (p, (pt, e_pos'), []);
neuper@41982
   251
val pIopt = get_pblID (pt,ip);
neuper@41982
   252
ip = ([],Res) (*false*);
neuper@41982
   253
val SOME pI = pIopt;
neuper@55412
   254
member op = [Pbl,Met] p_ andalso is_none (get_obj g_env pt (fst p));
neuper@41982
   255
			                          (*^^^^^^^^: Apply_Method without init_form*)
neuper@38065
   256
neuper@38065
   257
"----------- fun step -----------------------------------";
neuper@38065
   258
"----------- fun step -----------------------------------";
neuper@38065
   259
"----------- fun step -----------------------------------";
neuper@38065
   260
val p = e_pos'; val c = []; 
neuper@38065
   261
val (p,_,f,nxt,_,pt) = 
neuper@38065
   262
    CalcTreeTEST 
neuper@38065
   263
        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"], 
neuper@38065
   264
          ("Integrate", ["integrate","function"], ["diff","integration"]))];
neuper@38065
   265
"----- step 1: returns tac = Model_Problem ---";
neuper@38065
   266
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   267
"----- step 2: returns tac =  ---";
neuper@38065
   268
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   269
"----- step 3: returns tac =  ---";
neuper@38065
   270
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   271
"----- step 4: returns tac =  ---";
neuper@38065
   272
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   273
"----- step 5: returns tac =  ---";
neuper@38065
   274
neuper@48895
   275
(*========== inhibit exn AK110718 ==============================================
neuper@38066
   276
2002 stops here as well: TODO review actual arguments:
neuper@38065
   277
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   278
"----- step 6: returns tac =  ---";
neuper@38065
   279
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   280
"----- step 7: returns tac =  ---";
neuper@38065
   281
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   282
"----- step 8: returns tac =  ---";
neuper@38065
   283
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@48895
   284
========== inhibit exn AK110718 ==============================================*)
neuper@38065
   285
neuper@38065
   286
neuper@38065
   287
"----------- fun autocalc -------------------------------";
neuper@38065
   288
"----------- fun autocalc -------------------------------";
neuper@38065
   289
"----------- fun autocalc -------------------------------";
neuper@38065
   290
val p = e_pos'; val c = []; 
neuper@38065
   291
val (p,_,f,nxt,_,pt) = 
neuper@38065
   292
    CalcTreeTEST 
neuper@38065
   293
        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"], 
neuper@38065
   294
          ("Integrate",["integrate","function"], ["diff","integration"]))];
neuper@38066
   295
(*-----since Model_Problem + complete_mod_ in case cas of SOME-----*
neuper@38066
   296
modeling is skipped FIXME 
neuper@38066
   297
 *-----since Model_Problem + complete_mod_ in case cas of SOME-----*)
neuper@38065
   298
tracing "----- step 1 ---";
neuper@38065
   299
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   300
tracing "----- step 2 ---";
neuper@38065
   301
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   302
tracing "----- step 3 ---";
neuper@38065
   303
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   304
tracing "----- step 4 ---";
neuper@38065
   305
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   306
tracing "----- step 5 ---";
neuper@38065
   307
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   308
tracing "----- step 6 ---";
neuper@38065
   309
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   310
tracing "----- step 7 ---";
neuper@38065
   311
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   312
tracing "----- step 8 ---";
neuper@38065
   313
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@48895
   314
(*========== inhibit exn AK110718 ==============================================
neuper@42067
   315
WN110628: Integration does not work, see Knowledge/integrate.sml
neuper@42067
   316
neuper@41929
   317
if str = "end-of-calculation" andalso 
neuper@41929
   318
   term2str (get_obj g_res pt (fst p)) = "c + x + 1 / 3 * x ^^^ 3" then ()
neuper@38065
   319
else error "mathengine.sml -- fun autocalc -- end";
neuper@48895
   320
========== inhibit exn AK110718 ==============================================*)
neuper@38065
   321
neuper@38065
   322
neuper@38065
   323
"----------- fun autoCalculate -----------------------------------";
neuper@38065
   324
"----------- fun autoCalculate -----------------------------------";
neuper@38065
   325
"----------- fun autoCalculate -----------------------------------";
neuper@38065
   326
states := [];
neuper@38065
   327
CalcTree (*ATTENTION: encode_fmz ... unlike CalcTreeTEST*)
neuper@38065
   328
    [(["functionTerm (x^2 + 1)", "integrateBy x", "antiDerivative FF"],
neuper@38065
   329
      ("Integrate", ["integrate", "function"], ["diff", "integration"]))];
neuper@38065
   330
Iterator 1;
neuper@38065
   331
moveActiveRoot 1;
neuper@38066
   332
(*-----since Model_Problem + complete_mod_ in case cas of SOME-----*
neuper@38066
   333
modeling is skipped FIXME 
neuper@38066
   334
see test/../interface -- solve_linear as rootpbl FE -- for OTHER expl:
neuper@38066
   335
 setNextTactic 1 (Add_Given "equality (1 + -1 * 2 + x = 0)");
neuper@38066
   336
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1); (*equality added*);
neuper@38066
   337
neuper@38066
   338
 fetchProposedTactic 1;
neuper@38066
   339
 setNextTactic 1 (Add_Given "solveFor x");
neuper@38066
   340
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);
neuper@38066
   341
neuper@38066
   342
 fetchProposedTactic 1;
neuper@38066
   343
 setNextTactic 1 (Add_Find "solutions L");
neuper@38066
   344
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);
neuper@38066
   345
neuper@38066
   346
 fetchProposedTactic 1;
neuper@38066
   347
 setNextTactic 1 (Specify_Theory "Test");
neuper@38066
   348
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);
neuper@38066
   349
*-----since Model_Problem + complete_mod_ in case cas of SOME-----*)
neuper@38065
   350
autoCalculate 1 (Step 1); 
neuper@38065
   351
"----- step 1 ---";
neuper@38065
   352
autoCalculate 1 (Step 1);
neuper@38065
   353
"----- step 2 ---";
neuper@38065
   354
autoCalculate 1 (Step 1);
neuper@38065
   355
"----- step 3 ---";
neuper@38065
   356
autoCalculate 1 (Step 1);
neuper@38065
   357
"----- step 4 ---";
neuper@38065
   358
autoCalculate 1 (Step 1);
neuper@38065
   359
"----- step 5 ---";
neuper@38065
   360
autoCalculate 1 (Step 1);
neuper@38065
   361
"----- step 6 ---";
neuper@38065
   362
autoCalculate 1 (Step 1);
neuper@38065
   363
"----- step 7 ---";
neuper@38065
   364
autoCalculate 1 (Step 1);
neuper@38065
   365
"----- step 8 ---";
neuper@38065
   366
autoCalculate 1 (Step 1);
neuper@38065
   367
val (ptp as (_, p), _) = get_calc 1;
neuper@38065
   368
val (Form t,_,_) = pt_extract ptp;
akargl@42108
   369
neuper@38065
   370
if term2str t = "c + x + 1 / 3 * x ^^^ 3" andalso p = ([], Res) then ()
neuper@38065
   371
else error "mathengine.sml -- fun autoCalculate -- end";
neuper@38065
   372
neuper@42067
   373
"----------- fun autoCalculate..CompleteCalc ------------";
neuper@42067
   374
"----------- fun autoCalculate..CompleteCalc ------------";
neuper@42067
   375
"----------- fun autoCalculate..CompleteCalc ------------";
neuper@42067
   376
 states:=[];
neuper@42067
   377
 CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@42067
   378
   ("Test", ["sqroot-test","univariate","equation","test"],
neuper@42067
   379
    ["Test","squ-equ-test-subpbl1"]))];
neuper@42067
   380
 Iterator 1;
neuper@42103
   381
 moveActiveRoot 1; 
neuper@42103
   382
(*autoCalculate 1 CompleteCalc; (*WAS <SYSERROR>..<ERROR> error in kernel </ERROR>*)*)
neuper@42103
   383
"~~~~~ fun autoCalculate, args:"; val (cI, auto) = (1, CompleteCalc);
neuper@42103
   384
val pold = get_pos cI 1;
neuper@42103
   385
(*autocalc [] pold (get_calc cI) auto;  (*WAS Type unification failed
neuper@42103
   386
  Type error in application: incompatible operand type
neuper@42103
   387
  Operator:  solveFor :: real \<Rightarrow> una
neuper@42103
   388
  Operand:   x :: 'a *)*)
neuper@42103
   389
"~~~~~ fun autocalc, args:"; val (c, (pos as (_,p_)), ((pt,_), _), auto) 
neuper@42103
   390
                               = ([] : pos' list, pold, (get_calc cI), auto);
neuper@42103
   391
autoord auto > 3 andalso just_created (pt, pos); (*true*)
neuper@42103
   392
val ptp = all_modspec (pt, pos);
neuper@42103
   393
"TODO all_modspec: preconds for rootpbl";
neuper@42103
   394
(*all_solve auto c ptp; (*WAS Type unification failed...*)*)
neuper@42103
   395
"~~~~~ and all_solve, args:"; val (auto, c, (ptp as (pt, pos as (p,_)))) =
neuper@42103
   396
                                                                  (auto, c, ptp);
neuper@42103
   397
    val (_,_,mI) = get_obj g_spec pt p
neuper@42103
   398
    val ctxt = get_ctxt pt pos
neuper@42103
   399
    val (ttt, c', ptp) = nxt_solv (Apply_Method' (mI, NONE, e_istate, ctxt)) (e_istate, ctxt) ptp;
neuper@42103
   400
(* complete_solve auto (c @ c') ptp; (*WAS Type unification failed...*) pos = ([], Met)*)
neuper@42103
   401
"~~~~~ fun complete_solve, args:"; val (auto, c, (ptp as (_, p as (_,p_)))) =
neuper@42103
   402
                                                           (auto, (c @ c'), ptp);
neuper@42103
   403
p = ([], Res) (*false p = ([1], Frm)*);
neuper@42103
   404
member op = [Pbl,Met] p_ (*false*);
neuper@42103
   405
val (ttt, c', ptp') = nxt_solve_ ptp; (*ttt = Rewrite_Set "norm_equation"*)
neuper@42103
   406
(* complete_solve auto (c @ c') ptp'; (*WAS Type unification failed...*) pos = ([], Met)*)
neuper@42103
   407
"~~~~~ fun complete_solve, args:"; val (auto, c, (ptp as (_, p as (_,p_)))) =
neuper@42103
   408
                                                           (auto, (c @ c'), ptp');
neuper@42103
   409
p = ([], Res) (*false p = ([1], Res)*);
neuper@42103
   410
member op = [Pbl,Met] p_ (*false*);
neuper@42103
   411
val (ttt, c', ptp') = nxt_solve_ ptp; (*ttt = Rewrite_Set "Test_simplify"*)
neuper@42103
   412
(* complete_solve auto (c @ c') ptp'; (*WAS Type unification failed...*) pos = ([], Met)*)
neuper@42103
   413
"~~~~~ fun complete_solve, args:"; val (auto, c, (ptp as (_, p as (_,p_)))) =
neuper@42103
   414
                                                           (auto, (c @ c'), ptp');
neuper@42103
   415
p = ([], Res) (*false p = ([2], Res)*);
neuper@42103
   416
member op = [Pbl,Met] p_ (*false*);
neuper@42103
   417
val ((Subproblem _, tac_, (_, is))::_, c', ptp') = nxt_solve_ ptp;
neuper@42103
   418
autoord auto < 5 (*false*);
neuper@42103
   419
(* val ptp = all_modspec ptp' (*WAS Type unification failed...*)*)
neuper@42103
   420
"~~~~~ fun all_modspec, args:"; val (pt, pos as (p,_)) = (ptp');
neuper@42103
   421
    val (PblObj {fmz = (fmz_, _), origin = (pors, spec as (dI,pI,mI), hdl), ...}) = get_obj I pt p;
neuper@42103
   422
    val thy = assoc_thy dI;
neuper@42103
   423
	    val {ppc, ...} = get_met mI;
neuper@42103
   424
(*  val (mors, ctxt) = prep_ori fmz_ thy ppc; WAS Type unification failed because
neuper@48761
   425
val ctxt' = dI |> Thy_Info.get_theory |> Proof_Context.init_global;
neuper@42103
   426
(parseNEW ctxt' "x" |> the |> type_of) = TFree ("'a",[]); 
neuper@42103
   427
                                               ^^^^^ *)
neuper@42103
   428
(*vvv from:  | assod pt _ (Subproblem'...*)
neuper@42103
   429
    val (fmz_, vals) = oris2fmz_vals pors;
neuper@42103
   430
(**)
neuper@48761
   431
    val ctxt = dI |> Thy_Info.get_theory |> Proof_Context.init_global 
neuper@42103
   432
      |> declare_constraints' vals
neuper@42103
   433
(**)
neuper@42103
   434
(*^^^ from:  | assod pt _ (Subproblem'...*)
neuper@42103
   435
val [(1, [1], "#Given", dsc_eq, [equality]),
neuper@42103
   436
     (2, [1], "#Given", dsc_solvefor, [xxx]),
neuper@42103
   437
     (3, [1], "#Find", dsc_solutions, [x_i])] = pors;
neuper@42103
   438
if term2str xxx = "x" andalso type_of xxx = HOLogic.realT then () 
neuper@42103
   439
else error "autoCalculate..CompleteCalc: SubProblem broken 1";
neuper@42103
   440
    val pt = update_pblppc pt p (map (ori2Coritm ppc) pors);
neuper@42103
   441
	    val pt = update_metppc pt p (map (ori2Coritm ppc) pors);
neuper@42103
   442
	    val pt = update_spec pt p (dI,pI,mI);
neuper@42103
   443
    val pt = update_ctxt pt p ctxt;
neuper@42103
   444
"~~~~~ return to complete_solve, args:"; val (ptp) = (pt, (p,Met));
neuper@42103
   445
val (msg, c, (pt, p)) = complete_solve auto (c @ c') ptp;
neuper@42103
   446
if msg = "end-of-calculation" andalso c = [] andalso p = ([], Res) then ()
neuper@42103
   447
else error "autoCalculate..CompleteCalc: final result";
neuper@42103
   448
if terms2strs (get_assumptions_ pt p) = 
neuper@42105
   449
  ["matches (?a = ?b) (-1 + x = 0)", (*precond of submet*) 
neuper@42104
   450
   "x = 1",                          (*result of subpbl and rootpbl*)
neuper@42104
   451
   "precond_rootmet x"]              (*precond of rootmet*)
neuper@42103
   452
then () else error "autoCalculate..CompleteCalc: ctxt at final result";
neuper@42067
   453
 
akargl@42108
   454
neuper@42394
   455
"--------- search for Or_to_List ------------------------";
neuper@42394
   456
"--------- search for Or_to_List ------------------------";
t@42225
   457
"--------- search for Or_to_List ------------------------";
t@42225
   458
val fmz = ["equality (x^^^2 + 4*x + 5 = (2::real))", "solveFor x","solutions L"];
t@42225
   459
val (dI',pI',mI') = ("Isac", ["univariate","equation"], ["no_met"])
t@42225
   460
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
t@42225
   461
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   462
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   463
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   464
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   465
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   466
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   467
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   468
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   469
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   470
val (p,_,f,nxt,_,pt) = me nxt p [] pt; (*nxt = ("Subproblem"*)
t@42225
   471
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   472
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   473
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   474
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   475
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   476
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   477
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   478
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   479
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   480
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   481
"~~~~~ fun me, args:"; val (((_,tac):tac'_), (p:pos'), (_:NEW(*remove*)), (pt:ptree)) = 
t@42225
   482
                            (nxt, p, [], pt);
t@42225
   483
val ("ok", (_, _, ptp))  = locatetac tac (pt,p)
t@42225
   484
val (pt, p) = ptp;
t@42225
   485
"~~~~~ fun step, args:"; val (((ip as (_,p_)):pos'), ((ptp as (pt,p), tacis):calcstate)) =
t@42225
   486
                          (p, ((pt, e_pos'),[]));
t@42225
   487
val pIopt = get_pblID (pt,ip);
t@42225
   488
ip = ([],Res); (*false*)
t@42225
   489
ip = p; (*false*)
t@42225
   490
member op = [Pbl,Met] p_; (*false*)
t@42225
   491
"~~~~~ fun nxt_solve_, args:"; val (ptp as (pt, pos as (p,p_))) = (pt,ip);
t@42225
   492
e_metID = get_obj g_metID pt (par_pblobj pt p); (*false*)
t@42225
   493
val thy' = get_obj g_domID pt (par_pblobj pt p);
t@42225
   494
val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
t@42225
   495
val (tac_,is,(t,_)) = next_tac (thy',srls) (pt,pos) sc is; (*WAS Empty_Tac_: tac_*)
t@42225
   496
case tac_ of Or_to_List' _ => ()
t@42225
   497
| _ => error "Or_to_List broken ?"
akargl@42108
   498
neuper@42394
   499
neuper@42279
   500
"----------- check thy in CalcTreeTEST ------------------";
neuper@42279
   501
"----------- check thy in CalcTreeTEST ------------------";
neuper@42279
   502
"----------- check thy in CalcTreeTEST ------------------";
neuper@42279
   503
"WN1109 with Inverse_Z_Transform.thy when starting tests with CalcTreeTEST \n" ^
neuper@42279
   504
"we got the message: ME_Isa: thy 'Build_Inverse_Z_Transform' not in system \n" ^
neuper@42279
   505
"Below there are the steps which found out the reason: \n" ^
neuper@42279
   506
"store_pbt mistakenly stored that theory.";
neuper@48761
   507
val ctxt = Proof_Context.init_global @{theory Isac};
neuper@42279
   508
val SOME t = parseNEW ctxt "filterExpression (X  = 3 / (z - 1/4 + -1/8 * (1/(z::real))))";
neuper@42279
   509
val SOME t = parseNEW ctxt "stepResponse (x[n::real]::bool)";
akargl@42108
   510
neuper@42279
   511
val fmz = ["filterExpression (X  = 3 / (z - 1/4 + -1/8 * (1/(z::real))))", 
neuper@42279
   512
  "stepResponse (x[n::real]::bool)"];
neuper@42405
   513
val (dI,pI,mI) = ("Isac", ["Inverse", "Z_Transform", "SignalProcessing"], 
neuper@42405
   514
  ["SignalProcessing","Z_Transform","Inverse"]);
neuper@42279
   515
neuper@42279
   516
val (p,_,fb,nxt,_,pt)  = CalcTreeTEST [(fmz, (dI,pI,mI))]; 
neuper@42279
   517
(*WAS ME_Isa: thy 'Build_Inverse_Z_Transform' not in system*)
neuper@42279
   518
neuper@42279
   519
"~~~~~ T fun CalcTreeTEST, args:"; val [(fmz, sp):fmz] = [(fmz, (dI,pI,mI))];
neuper@42279
   520
"~~~~~ fun nxt_specify_init_calc, args:"; val (fmz, (dI, pI, mI)) = (fmz, sp);
neuper@42279
   521
	      val thy = assoc_thy dI;
neuper@42279
   522
    mI = ["no_met"]; (*false*)
neuper@42279
   523
		      val (pI, (pors, pctxt), mI) = (pI, get_pbt pI |> #ppc |> prep_ori fmz thy, mI)
neuper@42279
   524
	      val {cas, ppc, thy=thy',...} = get_pbt pI; (*take dI from _refined_ pbl*)
neuper@42279
   525
	"tracing bottom up"; Context.theory_name thy' = "Build_Inverse_Z_Transform"; (*WAS true*)
neuper@42279
   526
      val dI = theory2theory' (maxthy thy thy');
neuper@42279
   527
"tracing bottom up"; dI = "Build_Inverse_Z_Transform"; (*WAS true*)
neuper@42279
   528
    cas = NONE; (*true*)
neuper@42279
   529
	      val hdl = pblterm dI pI;
neuper@42279
   530
        val (pt, _) = cappend_problem e_ptree [] (e_istate, pctxt) (fmz, (dI, pI, mI))
neuper@42279
   531
				  (pors, (dI, pI, mI), hdl)
neuper@42279
   532
        val pt = update_ctxt pt [] pctxt;
neuper@42279
   533
neuper@42279
   534
"~~~~~ fun nxt_specif, args:"; val ((tac as Model_Problem), (pt, pos as (p,p_))) = (Model_Problem, (pt, ([], Pbl)));
neuper@42279
   535
      val PblObj {origin = (oris, ospec, _), probl, spec, ctxt, ...} = get_obj I pt p;
neuper@42279
   536
"tracing bottom up"; #1(ospec) = "Build_Inverse_Z_Transform"; (*WAS true*)
neuper@42279
   537
neuper@42279
   538
"~~~~~ fun some_spec, args:"; val ((odI, opI, omI), (dI, pI, mI)) = (ospec, spec);
neuper@42279
   539
dI = e_domID; (*true*)
neuper@42279
   540
odI = "Build_Inverse_Z_Transform"; (*true*)
neuper@42279
   541
dI = "e_domID"; (*true*)
neuper@42279
   542
"~~~~~ after fun some_spec:";
neuper@42279
   543
      val (dI, pI, mI) = some_spec ospec spec;
neuper@42279
   544
"tracing bottom up"; dI = "Build_Inverse_Z_Transform"; (*WAS true*)
neuper@42279
   545
      val thy = assoc_thy dI; (*caused ME_Isa: thy 'Build_Inverse_Z_Transform' not in system*)
neuper@42279
   546