test/Tools/isac/Interpret/mathengine.sml
author Thomas Leh <t.leh@gmx.at>
Thu, 28 Jul 2011 11:27:04 +0200
branchdecompose-isar
changeset 42225 7daeeee85596
parent 42109 cd33f1f80c8a
child 42279 e2759e250604
permissions -rw-r--r--
completed praktikum

in Test_Isac all 'use "...test ..." WITHOUT COMMENT
indicate: no outcommented errors.
other testfiles still contain "=== inhibit exn "
or contain "-.-.-isolate response.-.-."
akargl@42108
     1
akargl@42108
     2
akargl@42108
     3
neuper@37906
     4
(* test for sml/ME/mathengine.sml
neuper@37906
     5
   authors: Walther Neuper 2000, 2006
neuper@37906
     6
   (c) due to copyright terms
neuper@37906
     7
*)
neuper@38065
     8
"--------------------------------------------------------";
neuper@38065
     9
"table of contents --------------------------------------";
neuper@38065
    10
"--------------------------------------------------------";
neuper@41940
    11
"----------- change to parse ctxt -----------------------";
neuper@38065
    12
"----------- debugging setContext..pbl_ -----------------";
neuper@38065
    13
"----------- tryrefine ----------------------------------";
neuper@41982
    14
"----------- fun step: Apply_Method without init_form ---";
neuper@38065
    15
"----------- fun step -----------------------------------";
neuper@38065
    16
"----------- fun autocalc -------------------------------";
neuper@38065
    17
"----------- fun autoCalculate --------------------------";
neuper@42067
    18
"----------- fun autoCalculate..CompleteCalc ------------";
t@42225
    19
"----------- search for Or_to_List ----------------------";
neuper@38065
    20
"--------------------------------------------------------";
neuper@38065
    21
"--------------------------------------------------------";
neuper@38065
    22
"--------------------------------------------------------";
neuper@37906
    23
neuper@41940
    24
"----------- change to parse ctxt -----------------------";
neuper@41940
    25
"----------- change to parse ctxt -----------------------";
neuper@41940
    26
"----------- change to parse ctxt -----------------------";
neuper@41940
    27
"===== start calculation: from problem description (fmz) to origin";
neuper@41940
    28
val fmz = ["realTestGiven (((1+2)*4/3)^^^2)","realTestFind s"];
neuper@41940
    29
val (thyID, pblID, metID) =
neuper@41940
    30
  ("Test", ["calculate", "test"], ["Test", "test_calculate"]);
neuper@41940
    31
val (p,_,_,nxt,_,pt) = CalcTreeTEST [(fmz, (thyID, pblID, metID))];
neuper@41940
    32
"----- ";
neuper@41940
    33
(* call sequence: CalcTreeTEST 
neuper@41940
    34
                > nxt_specify_init_calc 
neuper@41940
    35
                > prep_ori
neuper@41940
    36
*)
neuper@41940
    37
val (thy, pbt) = (Thy_Info.get_theory thyID, (#ppc o get_pbt) pblID);
neuper@41940
    38
"----- in  prep_ori";
neuper@41940
    39
val ctxt = ProofContext.init_global thy;
neuper@41940
    40
bonzai@41952
    41
val ctopts = map (parseNEW ctxt) fmz;
bonzai@41952
    42
val dts = map (split_dts o the) ctopts;
neuper@41940
    43
(*split_dts:
neuper@41940
    44
(term * term list) list
neuper@41940
    45
        formulas: e.g. ((1+2)*4/3)^^^2
neuper@41940
    46
 description: e.g. realTestGiven
neuper@41940
    47
*)
neuper@41940
    48
 prep_ori;
neuper@41940
    49
(* FROM
neuper@41940
    50
val it = fn:
neuper@41940
    51
   string list -> theory -> (string * (term * 'a)) list -> ori list
neuper@41940
    52
TO
neuper@41940
    53
val it = fn:
neuper@41940
    54
   string list -> theory -> (string * (term * 'a)) list -> (ori list * ctxt)
neuper@41940
    55
AND
neuper@41940
    56
FROM val oris = prep_ori...
neuper@41940
    57
TO   val (oris, _) = prep_ori...
neuper@41940
    58
*)
neuper@41940
    59
"----- insert ctxt in ptree";
neuper@41940
    60
(* datatype ppobj
neuper@41940
    61
FROM loc   : istate option * istate option,
neuper@41940
    62
TO   loc   : (istate * ctxt) option * (istate * ctxt) option,
neuper@41940
    63
e.g.
neuper@41940
    64
FROM val pt = Nd (PblObj
neuper@41940
    65
       {.., loc = (SOME (ScrState ([], [], NONE, Const ("empty", "'a"), Sundef, false)),
neuper@41940
    66
          NONE),
neuper@41940
    67
TO   val pt = Nd (PblObj
neuper@41940
    68
       {.., loc = 
neuper@41940
    69
        ((SOME (ScrState ([], [], NONE, Const ("empty", "'a"), Sundef, false)), ctxt),
neuper@41940
    70
          NONE),
neuper@41940
    71
*)
neuper@41940
    72
neuper@41940
    73
"===== interactive specification: from origin to specification (probl)";
neuper@41940
    74
val (p,_,_,nxt,_,pt) = me nxt p [1] pt; 
neuper@41940
    75
  (*nxt =("Add_Given", Model_Problem)*)
neuper@41940
    76
val (p,_,_,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
    77
  (*nxt = ("Add_Find",Add_Find "realTestFind s") : string * tac*)
neuper@41940
    78
"----- ";
neuper@41940
    79
(* call sequence: me Model_Problem 
neuper@41940
    80
                > me ("Add_Given", Add_Given "realTestGiven (((1 + 2) * 4 / 3) ^^^ 2)")
neuper@41940
    81
                > locatetac tac
neuper@41940
    82
                > loc_specify_
neuper@41940
    83
                > specify          GET ctxt (stored in ctree)
neuper@41940
    84
                > specify_additem
neuper@41940
    85
                > appl_add
neuper@41940
    86
neuper@41940
    87
*)
neuper@41940
    88
"----- in appl_add";
neuper@41940
    89
(* FROM appl_add thy
neuper@41940
    90
   TO   appl_add ctxt
neuper@41940
    91
   FROM parse thy str
neuper@41940
    92
   TO   parseNEW ctxt str
neuper@41940
    93
*)
neuper@41940
    94
val (p,_,_,nxt,_,pt) = me nxt p [1] pt; 
neuper@41940
    95
  (*nxt = ("Specify_Theory",Specify_Theory "Test") : string * tac*)
neuper@41940
    96
val (p,_,_,nxt,_,pt) = me nxt p [1] pt; 
neuper@41940
    97
  (*nxt = ("Specify_Problem",Specify_Problem ["calculate","test"])*)
neuper@41940
    98
neuper@41940
    99
"===== end specify: from specification (probl) to guard (method)";
neuper@41940
   100
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   101
  (*nxt = ("Specify_Method",Specify_Method ("Test","test_calculate"))*)
neuper@41940
   102
neuper@41940
   103
"===== start interpretation: from guard to environment";
neuper@41940
   104
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   105
  (*nxt = ("Apply_Method",Apply_Method ("Test","test_calculate"))*)
neuper@41940
   106
"----- ";
neuper@41940
   107
(* call sequence: me ("Apply_Method",...
neuper@41940
   108
                > locatetac
neuper@41940
   109
                > loc_solve_
neuper@41940
   110
                > solve ("Apply_Method",...
neuper@41940
   111
*)
neuper@41940
   112
val ((_,tac), ptp) = (nxt, (pt, p));
neuper@41940
   113
locatetac tac (pt,p);
neuper@41940
   114
  val (mI, m) = mk_tac'_ tac;
neuper@41940
   115
  val Appl m = applicable_in p pt m;
neuper@41940
   116
  member op = specsteps mI;
neuper@41940
   117
  loc_solve_ (mI,m) ptp;
neuper@41940
   118
    val (m, (pt, pos)) = ((mI,m), ptp);
neuper@41940
   119
    solve m (pt, pos);
bonzai@41949
   120
      val ((_, m as Apply_Method' (mI, _, _, _)), (pt, (pos as (p,_)))) = 
neuper@41940
   121
        (m, (pt, pos));
neuper@41940
   122
      val {srls,...} = get_met mI;
neuper@41940
   123
      val PblObj{meth=itms,...} = get_obj I pt p;
neuper@41940
   124
      val thy' = get_obj g_domID pt p;
neuper@41940
   125
      val thy = assoc_thy thy';
bonzai@41949
   126
      val (is as ScrState (env,_,_,_,_,_), ctxt, sc) = init_scrstate thy itms mI;
neuper@41940
   127
neuper@41940
   128
"----- go on in the calculation";
neuper@41940
   129
val (p,_,f,nxt,_,pt) = me nxt pos [1] pt;
neuper@41940
   130
  (*nxt = ("Calculate",Calculate "PLUS")*)
neuper@41940
   131
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   132
  (*nxt = ("Calculate",Calculate "TIMES")*)
neuper@41940
   133
neuper@41940
   134
"===== input a formula to be derived from previous istate";
neuper@41940
   135
"----- appendFormula TODO: first repair error";
neuper@41940
   136
  val cs = ((pt,p),[]);
neuper@41940
   137
  val ("ok", cs' as (_,_,(pt,p))) = step p cs;
neuper@41943
   138
  val ifo = (* encode "4^^^2" \<longrightarrow>*) "4^2";
neuper@41940
   139
(*
neuper@41943
   140
  val ("no derivation found", (_,_,(pt, p))) = inform cs' ((*encode*) ifo);
neuper@41940
   141
  here ^^^^^^^^^^^^^^^^^^^^^ should be "ok"
neuper@41940
   142
*)
neuper@41940
   143
neuper@41940
   144
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   145
(*nxt = ("Calculate",Calculate "DIVIDE")*)
neuper@41940
   146
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   147
(*nxt = ("Calculate",Calculate "POWER")*)
neuper@41940
   148
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   149
(*nxt = ("Check_Postcond",Check_Postcond ["calculate","test"])*)
neuper@41940
   150
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41940
   151
(*nxt = ("End_Proof'",End_Proof')*)
neuper@41940
   152
if f = Form' (FormKF (~1,EdUndef,0,Nundef,"16")) then ()
neuper@41940
   153
else error "calculate.sml: script test_calculate changed behaviour";
neuper@41940
   154
neuper@41940
   155
"===== tactic Subproblem: from environment to origin";
neuper@41940
   156
"----- TODO";
neuper@41940
   157
neuper@37906
   158
neuper@38065
   159
"----------- debugging setContext..pbl_ -----------------";
neuper@38065
   160
"----------- debugging setContext..pbl_ -----------------";
neuper@38065
   161
"----------- debugging setContext..pbl_ -----------------";
neuper@37906
   162
states:=[];
neuper@41970
   163
CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@41970
   164
  ("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906
   165
   ["Test","squ-equ-test-subpbl1"]))];
neuper@37906
   166
Iterator 1;
neuper@37906
   167
moveActiveRoot 1; modelProblem 1;
neuper@37906
   168
neuper@37906
   169
val pos as (p,_) = ([],Pbl);
neuper@37906
   170
val guh = "pbl_equ_univ";
neuper@37906
   171
checkContext 1 pos guh;
neuper@37906
   172
val ((pt,_),_) = get_calc 1;
neuper@37906
   173
val pp = par_pblobj pt p;
neuper@37906
   174
val keID = guh2kestoreID guh;
neuper@37906
   175
case context_pbl keID pt pp of (true,["univariate", "equation"],_,_,_)=>()
neuper@38031
   176
| _ => error "mathengine.sml: context_pbl .. pbl_equ_univ checked";
neuper@37906
   177
neuper@37906
   178
case get_obj g_spec pt p of (_, ["e_pblID"], _) => ()
neuper@38031
   179
| _ => error "mathengine.sml: context_pbl .. pbl still empty";
neuper@37906
   180
setContext 1 pos guh;
neuper@37906
   181
val ((pt,_),_) = get_calc 1;
neuper@37906
   182
case get_obj g_spec pt p of (_, ["univariate", "equation"], _) => ()
neuper@38031
   183
| _ => error "mathengine.sml: context_pbl .. pbl set";
neuper@37906
   184
neuper@37906
   185
neuper@37906
   186
setContext 1 pos "met_eq_lin";
neuper@37906
   187
val ((pt,_),_) = get_calc 1;
neuper@37906
   188
case get_obj g_spec pt p of (_,  _, ["LinEq", "solve_lineq_equation"]) => ()
neuper@38031
   189
| _ => error "mathengine.sml: context_pbl .. pbl set";
neuper@37906
   190
neuper@37906
   191
neuper@38065
   192
"----------- tryrefine ----------------------------------";
neuper@38065
   193
"----------- tryrefine ----------------------------------";
neuper@38065
   194
"----------- tryrefine ----------------------------------";
neuper@37906
   195
states:=[];
neuper@37906
   196
CalcTree [(["equality (x/(x^2 - 6*x+9) - 1/(x^2 - 3*x) =1/x)", 
neuper@37906
   197
	    "solveFor x", "solutions L"],
neuper@38058
   198
	   ("RatEq",["univariate","equation"],["no_met"]))];
neuper@37906
   199
Iterator 1;
neuper@37906
   200
moveActiveRoot 1; autoCalculate 1 CompleteCalc;
neuper@37906
   201
neuper@37906
   202
refineProblem 1 ([1],Res) "pbl_equ_univ" 
neuper@37906
   203
(*gives "pbl_equ_univ_rat" correct*);
neuper@37906
   204
neuper@37906
   205
refineProblem 1 ([1],Res) (pblID2guh ["linear","univariate","equation"])
neuper@38065
   206
(*gives "pbl_equ_univ_lin" incorrect*);
neuper@38065
   207
neuper@38065
   208
neuper@41982
   209
"----------- fun step: Apply_Method without init_form ---";
neuper@41982
   210
"----------- fun step: Apply_Method without init_form ---";
neuper@41982
   211
"----------- fun step: Apply_Method without init_form ---";
neuper@41982
   212
val fmz = ["equality (x+1=(2::real))", "solveFor x","solutions L"];
neuper@41982
   213
val (dI',pI',mI') =
neuper@41982
   214
  ("Test", ["sqroot-test","univariate","equation","test"],
neuper@41982
   215
   ["Test","squ-equ-test-subpbl1"]);
neuper@41982
   216
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
neuper@41982
   217
"~~~~~ fun me, args:"; val (_,tac) = nxt;
neuper@41982
   218
val ("ok", (_, _, (pt, p))) = locatetac tac (pt,p);
neuper@41982
   219
"~~~~~ fun step, args:"; val (ip as (_,p_), ptp as (pt,p), tacis) = (p, (pt, e_pos'), []);
neuper@41982
   220
val pIopt = get_pblID (pt,ip);
neuper@41982
   221
ip = ([],Res) (*false*);
neuper@41982
   222
val SOME pI = pIopt;
akargl@42108
   223
(*=== inhibit exn 110718 ? =============================================================
akargl@42108
   224
neuper@41982
   225
member op = [Pbl,Met] p_ andalso is_none (get_obj g_env pt (fst p))
neuper@41982
   226
			                          (*^^^^^^^^: Apply_Method without init_form*)
akargl@42108
   227
akargl@42108
   228
===== inhibit exn 110718 ? ===========================================================*)
neuper@38065
   229
neuper@38065
   230
"----------- fun step -----------------------------------";
neuper@38065
   231
"----------- fun step -----------------------------------";
neuper@38065
   232
"----------- fun step -----------------------------------";
neuper@38065
   233
val p = e_pos'; val c = []; 
neuper@38065
   234
val (p,_,f,nxt,_,pt) = 
neuper@38065
   235
    CalcTreeTEST 
neuper@38065
   236
        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"], 
neuper@38065
   237
          ("Integrate", ["integrate","function"], ["diff","integration"]))];
neuper@38065
   238
"----- step 1: returns tac = Model_Problem ---";
neuper@38065
   239
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   240
"----- step 2: returns tac =  ---";
neuper@38065
   241
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   242
"----- step 3: returns tac =  ---";
neuper@38065
   243
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   244
"----- step 4: returns tac =  ---";
neuper@38065
   245
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   246
"----- step 5: returns tac =  ---";
neuper@38065
   247
akargl@42108
   248
(*========== inhibit exn 110718 =======================================================
neuper@38066
   249
2002 stops here as well: TODO review actual arguments:
neuper@38065
   250
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   251
"----- step 6: returns tac =  ---";
neuper@38065
   252
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   253
"----- step 7: returns tac =  ---";
neuper@38065
   254
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
neuper@38065
   255
"----- step 8: returns tac =  ---";
neuper@38065
   256
val (str, ([(tac, _, _)], _, (pt, p))) = step p ((pt, e_pos'), []);
akargl@42108
   257
============ inhibit exn 110718 =====================================================*)
neuper@38065
   258
neuper@38065
   259
neuper@38065
   260
"----------- fun autocalc -------------------------------";
neuper@38065
   261
"----------- fun autocalc -------------------------------";
neuper@38065
   262
"----------- fun autocalc -------------------------------";
neuper@38065
   263
val p = e_pos'; val c = []; 
neuper@38065
   264
val (p,_,f,nxt,_,pt) = 
neuper@38065
   265
    CalcTreeTEST 
neuper@38065
   266
        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"], 
neuper@38065
   267
          ("Integrate",["integrate","function"], ["diff","integration"]))];
neuper@38066
   268
(*-----since Model_Problem + complete_mod_ in case cas of SOME-----*
neuper@38066
   269
modeling is skipped FIXME 
neuper@38066
   270
 *-----since Model_Problem + complete_mod_ in case cas of SOME-----*)
neuper@38065
   271
tracing "----- step 1 ---";
neuper@38065
   272
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   273
tracing "----- step 2 ---";
neuper@38065
   274
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   275
tracing "----- step 3 ---";
neuper@38065
   276
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   277
tracing "----- step 4 ---";
neuper@38065
   278
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   279
tracing "----- step 5 ---";
neuper@38065
   280
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   281
tracing "----- step 6 ---";
neuper@38065
   282
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   283
tracing "----- step 7 ---";
neuper@38065
   284
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@38065
   285
tracing "----- step 8 ---";
neuper@38065
   286
val (str, cut, (pt, p)) = autocalc [] p ((pt, e_pos'),[]) (Step 1); show_pt pt;
neuper@42067
   287
(*========== inhibit exn 110628 ================================================
neuper@42067
   288
WN110628: Integration does not work, see Knowledge/integrate.sml
neuper@42067
   289
neuper@41929
   290
if str = "end-of-calculation" andalso 
neuper@41929
   291
   term2str (get_obj g_res pt (fst p)) = "c + x + 1 / 3 * x ^^^ 3" then ()
neuper@38065
   292
else error "mathengine.sml -- fun autocalc -- end";
neuper@42067
   293
============ inhibit exn 110628 ==============================================*)
neuper@38065
   294
neuper@38065
   295
neuper@38065
   296
"----------- fun autoCalculate -----------------------------------";
neuper@38065
   297
"----------- fun autoCalculate -----------------------------------";
neuper@38065
   298
"----------- fun autoCalculate -----------------------------------";
neuper@38065
   299
states := [];
neuper@38065
   300
CalcTree (*ATTENTION: encode_fmz ... unlike CalcTreeTEST*)
neuper@38065
   301
    [(["functionTerm (x^2 + 1)", "integrateBy x", "antiDerivative FF"],
neuper@38065
   302
      ("Integrate", ["integrate", "function"], ["diff", "integration"]))];
neuper@38065
   303
Iterator 1;
neuper@38065
   304
moveActiveRoot 1;
neuper@38066
   305
(*-----since Model_Problem + complete_mod_ in case cas of SOME-----*
neuper@38066
   306
modeling is skipped FIXME 
neuper@38066
   307
see test/../interface -- solve_linear as rootpbl FE -- for OTHER expl:
neuper@38066
   308
 setNextTactic 1 (Add_Given "equality (1 + -1 * 2 + x = 0)");
neuper@38066
   309
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1); (*equality added*);
neuper@38066
   310
neuper@38066
   311
 fetchProposedTactic 1;
neuper@38066
   312
 setNextTactic 1 (Add_Given "solveFor x");
neuper@38066
   313
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);
neuper@38066
   314
neuper@38066
   315
 fetchProposedTactic 1;
neuper@38066
   316
 setNextTactic 1 (Add_Find "solutions L");
neuper@38066
   317
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);
neuper@38066
   318
neuper@38066
   319
 fetchProposedTactic 1;
neuper@38066
   320
 setNextTactic 1 (Specify_Theory "Test");
neuper@38066
   321
 autoCalculate 1 (Step 1); refFormula 1 (get_pos 1 1);
neuper@38066
   322
*-----since Model_Problem + complete_mod_ in case cas of SOME-----*)
neuper@38065
   323
autoCalculate 1 (Step 1); 
neuper@38065
   324
"----- step 1 ---";
neuper@38065
   325
autoCalculate 1 (Step 1);
neuper@38065
   326
"----- step 2 ---";
neuper@38065
   327
autoCalculate 1 (Step 1);
neuper@38065
   328
"----- step 3 ---";
neuper@38065
   329
autoCalculate 1 (Step 1);
neuper@38065
   330
"----- step 4 ---";
neuper@38065
   331
autoCalculate 1 (Step 1);
neuper@38065
   332
"----- step 5 ---";
neuper@38065
   333
autoCalculate 1 (Step 1);
neuper@38065
   334
"----- step 6 ---";
neuper@38065
   335
autoCalculate 1 (Step 1);
neuper@38065
   336
"----- step 7 ---";
neuper@38065
   337
autoCalculate 1 (Step 1);
neuper@38065
   338
"----- step 8 ---";
neuper@38065
   339
autoCalculate 1 (Step 1);
neuper@38065
   340
val (ptp as (_, p), _) = get_calc 1;
neuper@38065
   341
val (Form t,_,_) = pt_extract ptp;
akargl@42108
   342
neuper@38065
   343
if term2str t = "c + x + 1 / 3 * x ^^^ 3" andalso p = ([], Res) then ()
neuper@38065
   344
else error "mathengine.sml -- fun autoCalculate -- end";
neuper@38065
   345
neuper@42067
   346
"----------- fun autoCalculate..CompleteCalc ------------";
neuper@42067
   347
"----------- fun autoCalculate..CompleteCalc ------------";
neuper@42067
   348
"----------- fun autoCalculate..CompleteCalc ------------";
neuper@42067
   349
 states:=[];
neuper@42067
   350
 CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@42067
   351
   ("Test", ["sqroot-test","univariate","equation","test"],
neuper@42067
   352
    ["Test","squ-equ-test-subpbl1"]))];
neuper@42067
   353
 Iterator 1;
neuper@42103
   354
 moveActiveRoot 1; 
neuper@42103
   355
(*autoCalculate 1 CompleteCalc; (*WAS <SYSERROR>..<ERROR> error in kernel </ERROR>*)*)
neuper@42103
   356
"~~~~~ fun autoCalculate, args:"; val (cI, auto) = (1, CompleteCalc);
neuper@42103
   357
val pold = get_pos cI 1;
neuper@42103
   358
(*autocalc [] pold (get_calc cI) auto;  (*WAS Type unification failed
neuper@42103
   359
  Type error in application: incompatible operand type
neuper@42103
   360
  Operator:  solveFor :: real \<Rightarrow> una
neuper@42103
   361
  Operand:   x :: 'a *)*)
neuper@42103
   362
"~~~~~ fun autocalc, args:"; val (c, (pos as (_,p_)), ((pt,_), _), auto) 
neuper@42103
   363
                               = ([] : pos' list, pold, (get_calc cI), auto);
neuper@42103
   364
autoord auto > 3 andalso just_created (pt, pos); (*true*)
neuper@42103
   365
val ptp = all_modspec (pt, pos);
neuper@42103
   366
"TODO all_modspec: preconds for rootpbl";
neuper@42103
   367
(*all_solve auto c ptp; (*WAS Type unification failed...*)*)
neuper@42103
   368
"~~~~~ and all_solve, args:"; val (auto, c, (ptp as (pt, pos as (p,_)))) =
neuper@42103
   369
                                                                  (auto, c, ptp);
neuper@42103
   370
    val (_,_,mI) = get_obj g_spec pt p
neuper@42103
   371
    val ctxt = get_ctxt pt pos
neuper@42103
   372
    val (ttt, c', ptp) = nxt_solv (Apply_Method' (mI, NONE, e_istate, ctxt)) (e_istate, ctxt) ptp;
neuper@42103
   373
(* complete_solve auto (c @ c') ptp; (*WAS Type unification failed...*) pos = ([], Met)*)
neuper@42103
   374
"~~~~~ fun complete_solve, args:"; val (auto, c, (ptp as (_, p as (_,p_)))) =
neuper@42103
   375
                                                           (auto, (c @ c'), ptp);
neuper@42103
   376
p = ([], Res) (*false p = ([1], Frm)*);
neuper@42103
   377
member op = [Pbl,Met] p_ (*false*);
neuper@42103
   378
val (ttt, c', ptp') = nxt_solve_ ptp; (*ttt = Rewrite_Set "norm_equation"*)
neuper@42103
   379
(* complete_solve auto (c @ c') ptp'; (*WAS Type unification failed...*) pos = ([], Met)*)
neuper@42103
   380
"~~~~~ fun complete_solve, args:"; val (auto, c, (ptp as (_, p as (_,p_)))) =
neuper@42103
   381
                                                           (auto, (c @ c'), ptp');
neuper@42103
   382
p = ([], Res) (*false p = ([1], Res)*);
neuper@42103
   383
member op = [Pbl,Met] p_ (*false*);
neuper@42103
   384
val (ttt, c', ptp') = nxt_solve_ ptp; (*ttt = Rewrite_Set "Test_simplify"*)
neuper@42103
   385
(* complete_solve auto (c @ c') ptp'; (*WAS Type unification failed...*) pos = ([], Met)*)
neuper@42103
   386
"~~~~~ fun complete_solve, args:"; val (auto, c, (ptp as (_, p as (_,p_)))) =
neuper@42103
   387
                                                           (auto, (c @ c'), ptp');
neuper@42103
   388
p = ([], Res) (*false p = ([2], Res)*);
neuper@42103
   389
member op = [Pbl,Met] p_ (*false*);
neuper@42103
   390
val ((Subproblem _, tac_, (_, is))::_, c', ptp') = nxt_solve_ ptp;
neuper@42103
   391
autoord auto < 5 (*false*);
neuper@42103
   392
(* val ptp = all_modspec ptp' (*WAS Type unification failed...*)*)
neuper@42103
   393
"~~~~~ fun all_modspec, args:"; val (pt, pos as (p,_)) = (ptp');
neuper@42103
   394
    val (PblObj {fmz = (fmz_, _), origin = (pors, spec as (dI,pI,mI), hdl), ...}) = get_obj I pt p;
neuper@42103
   395
    val thy = assoc_thy dI;
neuper@42103
   396
	    val {ppc, ...} = get_met mI;
neuper@42103
   397
(*  val (mors, ctxt) = prep_ori fmz_ thy ppc; WAS Type unification failed because
neuper@42103
   398
val ctxt' = dI |> Thy_Info.get_theory |> ProofContext.init_global;
neuper@42103
   399
(parseNEW ctxt' "x" |> the |> type_of) = TFree ("'a",[]); 
neuper@42103
   400
                                               ^^^^^ *)
neuper@42103
   401
(*vvv from:  | assod pt _ (Subproblem'...*)
neuper@42103
   402
    val (fmz_, vals) = oris2fmz_vals pors;
neuper@42103
   403
(**)
neuper@42103
   404
    val ctxt = dI |> Thy_Info.get_theory |> ProofContext.init_global 
neuper@42103
   405
      |> declare_constraints' vals
neuper@42103
   406
(**)
neuper@42103
   407
(*^^^ from:  | assod pt _ (Subproblem'...*)
neuper@42103
   408
val [(1, [1], "#Given", dsc_eq, [equality]),
neuper@42103
   409
     (2, [1], "#Given", dsc_solvefor, [xxx]),
neuper@42103
   410
     (3, [1], "#Find", dsc_solutions, [x_i])] = pors;
neuper@42103
   411
if term2str xxx = "x" andalso type_of xxx = HOLogic.realT then () 
neuper@42103
   412
else error "autoCalculate..CompleteCalc: SubProblem broken 1";
neuper@42103
   413
    val pt = update_pblppc pt p (map (ori2Coritm ppc) pors);
neuper@42103
   414
	    val pt = update_metppc pt p (map (ori2Coritm ppc) pors);
neuper@42103
   415
	    val pt = update_spec pt p (dI,pI,mI);
neuper@42103
   416
    val pt = update_ctxt pt p ctxt;
neuper@42103
   417
"~~~~~ return to complete_solve, args:"; val (ptp) = (pt, (p,Met));
neuper@42103
   418
val (msg, c, (pt, p)) = complete_solve auto (c @ c') ptp;
neuper@42103
   419
if msg = "end-of-calculation" andalso c = [] andalso p = ([], Res) then ()
neuper@42103
   420
else error "autoCalculate..CompleteCalc: final result";
neuper@42103
   421
if terms2strs (get_assumptions_ pt p) = 
neuper@42105
   422
  ["matches (?a = ?b) (-1 + x = 0)", (*precond of submet*) 
neuper@42104
   423
   "x = 1",                          (*result of subpbl and rootpbl*)
neuper@42104
   424
   "precond_rootmet x"]              (*precond of rootmet*)
neuper@42103
   425
then () else error "autoCalculate..CompleteCalc: ctxt at final result";
neuper@42067
   426
 
akargl@42108
   427
t@42225
   428
"----------- search for Or_to_List ----------------------";
t@42225
   429
"----------- search for Or_to_List ----------------------";
t@42225
   430
"--------- search for Or_to_List ------------------------";
t@42225
   431
val fmz = ["equality (x^^^2 + 4*x + 5 = (2::real))", "solveFor x","solutions L"];
t@42225
   432
val (dI',pI',mI') = ("Isac", ["univariate","equation"], ["no_met"])
t@42225
   433
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
t@42225
   434
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   435
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   436
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   437
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   438
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   439
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   440
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   441
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   442
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   443
val (p,_,f,nxt,_,pt) = me nxt p [] pt; (*nxt = ("Subproblem"*)
t@42225
   444
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   445
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   446
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   447
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   448
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   449
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   450
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   451
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   452
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   453
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
t@42225
   454
"~~~~~ fun me, args:"; val (((_,tac):tac'_), (p:pos'), (_:NEW(*remove*)), (pt:ptree)) = 
t@42225
   455
                            (nxt, p, [], pt);
t@42225
   456
val ("ok", (_, _, ptp))  = locatetac tac (pt,p)
t@42225
   457
val (pt, p) = ptp;
t@42225
   458
"~~~~~ fun step, args:"; val (((ip as (_,p_)):pos'), ((ptp as (pt,p), tacis):calcstate)) =
t@42225
   459
                          (p, ((pt, e_pos'),[]));
t@42225
   460
val pIopt = get_pblID (pt,ip);
t@42225
   461
ip = ([],Res); (*false*)
t@42225
   462
ip = p; (*false*)
t@42225
   463
member op = [Pbl,Met] p_; (*false*)
t@42225
   464
"~~~~~ fun nxt_solve_, args:"; val (ptp as (pt, pos as (p,p_))) = (pt,ip);
t@42225
   465
e_metID = get_obj g_metID pt (par_pblobj pt p); (*false*)
t@42225
   466
val thy' = get_obj g_domID pt (par_pblobj pt p);
t@42225
   467
val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
t@42225
   468
val (tac_,is,(t,_)) = next_tac (thy',srls) (pt,pos) sc is; (*WAS Empty_Tac_: tac_*)
t@42225
   469
case tac_ of Or_to_List' _ => ()
t@42225
   470
| _ => error "Or_to_List broken ?"
akargl@42108
   471
akargl@42108
   472