test/Tools/isac/Interpret/script.sml
author Walther Neuper <neuper@ist.tugraz.at>
Wed, 05 Dec 2012 15:56:38 +0100
changeset 48790 98df8f6dc3f9
parent 42438 31e1aa39b5cb
child 48895 35751d90365e
permissions -rw-r--r--
Test_Isac.thy works until "Interpret/mstools.sml"
neuper@41999
     1
(* Title:  test/../script.sml
neuper@41999
     2
   Author: Walther Neuper 050908
neuper@37906
     3
   (c) copyright due to lincense terms.
neuper@37906
     4
*)
neuper@37906
     5
"-----------------------------------------------------------------";
neuper@37906
     6
"table of contents -----------------------------------------------";
neuper@37906
     7
"-----------------------------------------------------------------";
neuper@37906
     8
"----------- WN0509 why does next_tac doesnt find Substitute -----";
neuper@37906
     9
"----------- WN0509 Substitute 2nd part --------------------------";
neuper@37906
    10
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@41969
    11
"----------- fun init_form, fun get_stac -------------------------";
neuper@41972
    12
"----------- x+1=2 start SubProblem 'stac2tac_ TODO: no match for ";
neuper@42283
    13
"----------- Take as 1st stac in script --------------------------";
neuper@42394
    14
"----------- 'trace_script' from outside 'fun me '----------------";
neuper@42438
    15
"----------- fun sel_rules ---------------------------------------";
neuper@42438
    16
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@37906
    17
"-----------------------------------------------------------------";
neuper@37906
    18
"-----------------------------------------------------------------";
neuper@37906
    19
"-----------------------------------------------------------------";
neuper@37906
    20
neuper@42438
    21
val thy =  @{theory Isac};
neuper@42281
    22
neuper@37906
    23
"----------- WN0509 why does next_tac doesnt find Substitute -----";
neuper@37906
    24
"----------- WN0509 why does next_tac doesnt find Substitute -----";
neuper@37906
    25
"----------- WN0509 why does next_tac doesnt find Substitute -----";
neuper@37906
    26
neuper@37906
    27
(*replace met 'IntegrierenUndKonstanteBestimmen' with this script...*)
neuper@37906
    28
val str = (*#1#*)
akargl@42169
    29
"Script BiegelinieScript                                             " ^ 
akargl@42169
    30
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^ 
akargl@42169
    31
"(rb_rb::bool list) (rm_rm::bool list) =                             " ^
akargl@42169
    32
"  (let q___q = Take (M_b v_v = q__q);                               " ^
akargl@42169
    33
"       (M1__M1::bool) = ((Substitute [v_v = 0])) q___q              " ^
akargl@42169
    34
" in True)";
akargl@42169
    35
neuper@37906
    36
val sc' = ((inst_abs thy) o term_of o the o (parse thy)) str;
neuper@37906
    37
neuper@37906
    38
val str = (*#2#*)
akargl@42169
    39
"Script BiegelinieScript                                             " ^ 
akargl@42169
    40
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^ 
akargl@42169
    41
"(rb_rb::bool list) (rm_rm::bool list) =                             " ^ 
akargl@42169
    42
"       (M1__M1::bool) = ((Substitute [v_v = 0]) @@                  " ^
akargl@42169
    43
" in True)";(*..doesnt find Substitute with ..@@ !!!!!!!!!!!!!!!!!!!!!*)
neuper@37906
    44
neuper@37906
    45
val str = (*#3#*)
akargl@42169
    46
"Script BiegelinieScript                                             " ^
akargl@42169
    47
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^
akargl@42169
    48
"(rb_rb::bool list) (rm_rm::bool list) =                             " ^
akargl@42169
    49
"  (let q___q = Take (q_q v_v = q__q);                               " ^
akargl@42169
    50
"      (M1__M1::bool) = Substitute [v_v = 0]      q___q ;            " ^
akargl@42169
    51
"       M1__m1 =        Substitute [M_b 0 = 0]  M1__M1               " ^
akargl@42169
    52
" in True)"
neuper@37906
    53
;
neuper@37906
    54
val str = (*#4#*)
akargl@42169
    55
"Script BiegelinieScript                                             " ^
akargl@42169
    56
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^
akargl@42169
    57
"(rb_rb::bool list) (rm_rm::bool list) =                             " ^
akargl@42169
    58
"  (let q___q = Take (q_q v_v = q__q);                               " ^
akargl@42169
    59
"      (M1__M1::bool) = Substitute [v_v = 0]      q___q ;            " ^
akargl@42169
    60
"       M1__M1 =        Substitute [v_v = 1]      q___q ;            " ^
akargl@42169
    61
"       M1__M1 =        Substitute [v_v = 2]      q___q ;            " ^
akargl@42169
    62
"       M1__M1 =        Substitute [v_v = 3]      q___q ;            " ^
akargl@42169
    63
"       M1__M1 =        Substitute [M_b 0 = 0]  M1__M1               " ^
akargl@42169
    64
" in True)"
neuper@37906
    65
;
neuper@37906
    66
val str = (*#5#*)
akargl@42169
    67
"Script BiegelinieScript                                             " ^
akargl@42169
    68
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^
akargl@42169
    69
"(rb_rb::bool list) (rm_rm::bool list) =                             " ^
akargl@42169
    70
"  (let q___q = Take (M_b v_v = q__q);                               " ^
akargl@42169
    71
"      (M1__M1::bool) = Substitute [v_v = 0]      q___q ;            " ^
akargl@42169
    72
"       M2__M2 = Take q___q ;                                        " ^
akargl@42169
    73
"       M2__M2 =        Substitute [v_v = 2]      q___q              " ^
akargl@42169
    74
" in True)"
neuper@37906
    75
;
akargl@42145
    76
neuper@37906
    77
val sc' = ((inst_abs thy) o term_of o the o (parse thy)) str;
neuper@37906
    78
atomty sc';
akargl@42145
    79
neuper@48790
    80
val {scr = Prog sc,...} = get_met ["IntegrierenUndKonstanteBestimmen"];
neuper@37906
    81
(*---------------------------------------------------------------------
neuper@38031
    82
if sc = sc' then () else error"script.sml, doesnt find Substitute #1";
neuper@37906
    83
---------------------------------------------------------------------*)
neuper@37906
    84
neuper@37906
    85
"----------- WN0509 Substitute 2nd part --------------------------";
neuper@37906
    86
"----------- WN0509 Substitute 2nd part --------------------------";
neuper@37906
    87
"----------- WN0509 Substitute 2nd part --------------------------";
neuper@37906
    88
(*replace met 'IntegrierenUndKonstanteBestimmen' with this script...*)
neuper@37906
    89
val str = (*Substitute ; Substitute works*)
akargl@42169
    90
"Script BiegelinieScript                                             " ^
akargl@42169
    91
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^
akargl@42169
    92
"(rb_rb::bool list) (rm_rm::bool list) =                             "^
neuper@37906
    93
(*begin after the 2nd integrate*)
akargl@42169
    94
"  (let M__M = Take (M_b v_v = q__q);                                " ^
akargl@42169
    95
"       e1__e1 = nth_nth 1 rm_rm ;                                   " ^
akargl@42169
    96
"       (x1__x1::real) = argument_in (lhs e1__e1);                   " ^
akargl@42169
    97
"       (M1__M1::bool) = Substitute [v_v = x1__x1] M__M;             " ^
akargl@42169
    98
"        M1__M1        = Substitute [e1__e1] M1__M1                  " ^
akargl@42169
    99
" in True)"
neuper@37906
   100
;
neuper@37906
   101
val sc' = ((inst_abs thy) o term_of o the o (parse thy)) str;
neuper@37906
   102
atomty sc';
akargl@42145
   103
neuper@37906
   104
val str = (*Substitute @@ Substitute does NOT work???*)
akargl@42169
   105
"Script BiegelinieScript                                             " ^
akargl@42169
   106
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real)              " ^
akargl@42169
   107
"(rb_rb::bool list) (rm_rm::bool list) =                             "^
neuper@37906
   108
(*begin after the 2nd integrate*)
akargl@42169
   109
"  (let M__M = Take (M_b v_v = q__q);                                " ^
akargl@42169
   110
"       e1__e1 = nth_nth 1 rm_rm ;                                   " ^
akargl@42169
   111
"       (x1__x1::real) = argument_in (lhs e1__e1);                   " ^
akargl@42169
   112
"       (M1__M1::bool) = ((Substitute [v_v = x1__x1]) @@             " ^
akargl@42169
   113
"                       (Substitute [e1__e1]))        M__M           " ^
akargl@42169
   114
" in True)"
neuper@37906
   115
;
neuper@37906
   116
neuper@48790
   117
val {scr = Prog sc,...} = get_met ["IntegrierenUndKonstanteBestimmen"];
neuper@37906
   118
(*---------------------------------------------------------------------
akargl@42169
   119
if sc = sc' then () else error "script.sml, doesnt find Substitute #1";
neuper@37906
   120
---------------------------------------------------------------------*)
neuper@37906
   121
val fmz = ["Traegerlaenge L",
neuper@37906
   122
	   "Streckenlast (- q_0 * x^^^2 / 2 + x * c + c_2)",
neuper@37906
   123
	   "Biegelinie y",
neuper@37906
   124
	   "RandbedingungenBiegung [y 0 = 0, y L = 0]",
neuper@37906
   125
	   "RandbedingungenMoment [M_b 0 = 0, M_b L = 0]",
neuper@37906
   126
	   "FunktionsVariable x"];
neuper@37906
   127
val (dI',pI',mI') =
neuper@38058
   128
  ("Biegelinie",["MomentBestimmte","Biegelinien"],
neuper@37906
   129
   ["IntegrierenUndKonstanteBestimmen"]);
neuper@37906
   130
val p = e_pos'; val c = [];
neuper@37906
   131
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
akargl@42169
   132
val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
akargl@42169
   133
val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
akargl@42169
   134
val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
akargl@42169
   135
val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
akargl@42169
   136
val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@37906
   137
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@37906
   138
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@37906
   139
case nxt of (_, Apply_Method ["IntegrierenUndKonstanteBestimmen"]) => ()
neuper@38031
   140
	  | _ => error "script.sml, doesnt find Substitute #2";
akargl@42169
   141
(* ERROR: caused by f2str f *)
neuper@37906
   142
trace_rewrite:=true;
neuper@42387
   143
neuper@37906
   144
val (p,_,f,nxt,_,pt) = me nxt p c pt; f2str f(*------------------------*);
neuper@37906
   145
trace_rewrite:=false;
akargl@42169
   146
  (*Exception TYPE raised:
akargl@42169
   147
    Illegal type for constant "HOL.eq" :: "[real, bool] => bool"
akargl@42169
   148
    Atools.argument'_in (Tools.lhs (ListG.nth_ (1 + -1 + -1) [])) =
akargl@42169
   149
    ListG.nth_ (1 + -1 + -1) []
akargl@42169
   150
    Exception-
akargl@42169
   151
       TYPE
akargl@42169
   152
          ("Illegal type for constant \"op =\" :: \"[real, bool] => bool\"",
akargl@42169
   153
             [],
akargl@42169
   154
             [Const ("HOL.Trueprop", "bool => prop") $
akargl@42169
   155
                   (Const ("HOL.eq", "[RealDef.real, bool] => bool") $ ... $ ...)])
akargl@42169
   156
       raised
akargl@42169
   157
    ... this was because eval_argument_in took "argument_in (lhs (M_b 0 = 0))"
akargl@42169
   158
    ie. the argument had not been simplified before          ^^^^^^^^^^^^^^^
akargl@42169
   159
    thus corrected eval_argument_in OK*)
neuper@37906
   160
val {srls,...} = get_met ["IntegrierenUndKonstanteBestimmen"];
neuper@42387
   161
val e1__e1 = str2term "NTH 1 [M_b 0 = 0, M_b L = 0]";
akargl@42169
   162
val e1__e1 = eval_listexpr_ @{theory Biegelinie} srls e1__e1; term2str e1__e1;
akargl@42169
   163
if term2str e1__e1 = "M_b 0 = 0"
akargl@42169
   164
  then ()
akargl@42169
   165
  else error"script.sml diff.beh. eval_listexpr_ nth_nth 1 [...";
neuper@37906
   166
neuper@37906
   167
(*TODO.WN050913 ??? doesnt eval_listexpr_ go into subterms ???...
akargl@42169
   168
 val x1__ = str2term"argument_in (lhs (M_b 0 = 0))";
akargl@42169
   169
 val x1__ = eval_listexpr_ Biegelinie.thy srls x1__; term2str x1__;
akargl@42169
   170
 (*no rewrite*)
akargl@42169
   171
 calculate_ Biegelinie.thy ("Tools.lhs", eval_rhs"eval_lhs_") x1__;
akargl@42169
   172
 val SOME (str, t) = eval_lhs 0 "Tools.lhs" (str2term"lhs (M_b 0 = 0)") 0;*)
neuper@37906
   173
akargl@42169
   174
val l__l = str2term "lhs (M_b 0 = 0)";
akargl@42169
   175
val l__l = eval_listexpr_ @{theory Biegelinie} srls l__l; term2str l__l;
akargl@42169
   176
val SOME (str, t) = eval_lhs 0 "Tools.lhs" (str2term "lhs (M_b 0 = 0)") 0;
neuper@42387
   177
akargl@42169
   178
val (p,_,f,nxt,_,pt) = me nxt p c pt; f2str f (*------------------------*);
neuper@42387
   179
if nxt = ("Rewrite", Rewrite ("Belastung_Querkraft", "Biegelinie.Belastung_Querkraft"))
neuper@42387
   180
then () else error "";
neuper@37906
   181
neuper@37906
   182
neuper@37906
   183
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@37906
   184
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@37906
   185
"----------- fun sel_appl_atomic_tacs ----------------------------";
akargl@42145
   186
neuper@37906
   187
states:=[];
neuper@41970
   188
CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@41970
   189
  ("Test", ["sqroot-test","univariate","equation","test"],
neuper@37906
   190
   ["Test","squ-equ-test-subpbl1"]))];
neuper@37906
   191
Iterator 1;
neuper@37906
   192
moveActiveRoot 1;
neuper@37906
   193
autoCalculate 1 CompleteCalcHead;
neuper@37906
   194
autoCalculate 1 (Step 1);
neuper@37906
   195
autoCalculate 1 (Step 1);
neuper@37906
   196
val ((pt, p), _) = get_calc 1; show_pt pt;
neuper@37906
   197
val appltacs = sel_appl_atomic_tacs pt p;
akargl@42145
   198
akargl@42169
   199
(*========== inhibit exn AK110721 ================================================
akargl@42169
   200
(*(*These SHOULD be the same*)
akargl@42145
   201
 print_depth(999);
akargl@42145
   202
 appltacs;
akargl@42145
   203
 [Rewrite ("radd_commute", "?m + ?n = ?n + ?m"),
akargl@42145
   204
    Rewrite ("radd_assoc", "?m + ?n + ?k = ?m + (?n + ?k)"),
akargl@42169
   205
    Calculate "TIMES"];*)
akargl@42145
   206
neuper@37906
   207
if appltacs = 
neuper@37906
   208
   [Rewrite ("radd_commute", "?m + ?n = ?n + ?m"),
neuper@37906
   209
    Rewrite ("radd_assoc", "?m + ?n + ?k = ?m + (?n + ?k)"),
akargl@42169
   210
    Calculate "TIMES"] then ()
neuper@38031
   211
else error "script.sml fun sel_appl_atomic_tacs diff.behav.";
akargl@42169
   212
========== inhibit exn AK110721 ================================================*)
neuper@37906
   213
neuper@37906
   214
trace_script := true;
neuper@37906
   215
trace_script := false;
neuper@37906
   216
applyTactic 1 p (hd appltacs);
neuper@37906
   217
val ((pt, p), _) = get_calc 1; show_pt pt;
neuper@37906
   218
val appltacs = sel_appl_atomic_tacs pt p;
neuper@37906
   219
akargl@42169
   220
(*(* AK110722 Debugging start*)
akargl@42169
   221
(*applyTactic 1 p (hd appltacs); WAS assy: call by ([3], Pbl)*)
akargl@42169
   222
"~~~~~ fun applyTactic, args:"; val ((cI:calcID), ip, tac) = (1, p, (hd appltacs));
akargl@42169
   223
val ((pt, _), _) = get_calc cI;
akargl@42169
   224
val p = get_pos cI 1;
akargl@42169
   225
akargl@42169
   226
locatetac;
akargl@42169
   227
locatetac tac;
akargl@42169
   228
locatetac tac;
akargl@42169
   229
akargl@42169
   230
(*locatetac tac (pt, ip); (*WAS assy: call by ([3], Pbl)*)*)
akargl@42169
   231
"~~~~~ fun locatetac, args:"; val (tac, (ptp as (pt, p))) = (tac, (pt, ip));
akargl@42169
   232
val (mI,m) = mk_tac'_ tac;
akargl@42169
   233
akargl@42169
   234
applicable_in p pt m ; (*Appl*)
akargl@42169
   235
akargl@42169
   236
member op = specsteps mI; (*false*)
akargl@42169
   237
akargl@42169
   238
(*loc_solve_ (mI,m) ptp; (*Argument: (mI, m) : string * tac Reason: Can't unify tac_ with tac*)*)
akargl@42169
   239
loc_solve_; (*without _ ??? result is -> lOc writing error ???*)
akargl@42169
   240
(*val it = fn: string * tac_ -> ptree * (pos * pos_) -> lOc_*)
akargl@42169
   241
(mI,m); (*string * tac*)
akargl@42169
   242
ptp (*ptree * pos'*);
akargl@42169
   243
akargl@42169
   244
"~~~~~ fun loc_solve_, args:"; val (m ,(pt,pos)) = ((mI,m), ptp);
akargl@42169
   245
(*val (msg, cs') = solve m (pt, pos);
akargl@42169
   246
(*Argument: m : string * tac Reason: Can't unify tac_ with tac (Different type constructors)*)*)
akargl@42169
   247
akargl@42169
   248
m;
akargl@42169
   249
(pt, pos);
akargl@42169
   250
akargl@42169
   251
"~~~~~ fun solve, args:"; val ((mI,m), (pt, po as (p,p_))) = (m, (pt, pos));
akargl@42169
   252
(*solve has still p_!!!!!!!!!!!!!!!!!!!!!!!!!*)
akargl@42169
   253
akargl@42169
   254
e_metID = get_obj g_metID pt (par_pblobj pt p); (*false*)
akargl@42169
   255
val ctxt = get_ctxt pt po;
akargl@42169
   256
akargl@42169
   257
(*generate1 (assoc_thy (get_obj g_domID pt (par_pblobj pt p))) m (e_istate, ctxt) (p,p_) pt;
akargl@42169
   258
(*Argument: m : tac Reason: Can't unify tac_ with tac (Different type constructors)*)*)
akargl@42169
   259
(assoc_thy (get_obj g_domID pt (par_pblobj pt p)));
akargl@42169
   260
assoc_thy;
akargl@42169
   261
akargl@42169
   262
(*not finished - error continues in fun generate1*)
akargl@42169
   263
(* AK110722 Debugging end*)*)
akargl@42169
   264
(*========== inhibit exn AK110721 ================================================
akargl@42169
   265
"----- WN080102 these vvv do not work, because locatetac starts the search" ^
akargl@42169
   266
 "1 stac too low";
akargl@42169
   267
(* AK110722 ERROR: assy: call by ([3], Pbl) *)
neuper@37906
   268
applyTactic 1 p (hd appltacs);
akargl@42169
   269
============ inhibit exn AK110721 ==============================================*)
akargl@42169
   270
neuper@37906
   271
autoCalculate 1 CompleteCalc;
neuper@37906
   272
neuper@41969
   273
"----------- fun init_form, fun get_stac -------------------------";
neuper@41969
   274
"----------- fun init_form, fun get_stac -------------------------";
neuper@41969
   275
"----------- fun init_form, fun get_stac -------------------------";
neuper@41969
   276
val fmz = ["equality (x+1=(2::real))", "solveFor x", "solutions L"];
neuper@41969
   277
val (dI',pI',mI') = ("Test",["sqroot-test","univariate","equation","test"],
neuper@41969
   278
  ["Test","squ-equ-test-subpbl1"]);
neuper@41969
   279
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
neuper@41969
   280
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   281
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   282
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   283
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   284
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   285
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   286
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@41969
   287
"~~~~~ fun me, args:"; val (_,tac) = nxt;
neuper@41969
   288
"~~~~~ fun locatetac, args:"; val ptp as (pt, p) = (pt, p);
neuper@41969
   289
val (mI,m) = mk_tac'_ tac;
neuper@41969
   290
val Appl m = applicable_in p pt m;
neuper@41969
   291
member op = specsteps mI; (*false*)
neuper@41969
   292
"~~~~~ fun loc_solve_, args:"; val (m, (pt, pos) ) = ((mI, m), ptp);
neuper@41969
   293
"~~~~~ fun solve, args:"; val (("Apply_Method", Apply_Method' (mI,_,_,_)), pos as (p,_)) = (m, pos);
neuper@41969
   294
val {srls, ...} = get_met mI;
neuper@41969
   295
val PblObj {meth=itms, ...} = get_obj I pt p;
neuper@41969
   296
val thy' = get_obj g_domID pt p;
neuper@41969
   297
val thy = assoc_thy thy';
neuper@41969
   298
val (is as ScrState (env,_,_,_,_,_), ctxt, sc) = init_scrstate thy itms mI;
neuper@41969
   299
val ini = init_form thy sc env;
neuper@48790
   300
"----- fun init_form, args:"; val (Prog sc) = sc;
neuper@41969
   301
"----- fun get_stac, args:"; val (y, h $ body) = (thy, sc);
neuper@41969
   302
case get_stac thy sc of SOME (Free ("e_e", _)) => ()
neuper@41969
   303
| _ => error "script: Const (?, ?) in script (as term) changed";;
neuper@41969
   304
neuper@41972
   305
neuper@41972
   306
"----------- x+1=2 start SubProblem 'stac2tac_ TODO: no match for ";
neuper@41972
   307
"----------- x+1=2 start SubProblem 'stac2tac_ TODO: no match for ";
neuper@41972
   308
"----------- x+1=2 start SubProblem 'stac2tac_ TODO: no match for ";
neuper@41972
   309
val fmz = ["equality (x+1=(2::real))", "solveFor x","solutions L"];
neuper@41972
   310
val (dI',pI',mI') =
neuper@41972
   311
  ("Test", ["sqroot-test","univariate","equation","test"],
neuper@41972
   312
   ["Test","squ-equ-test-subpbl1"]);
neuper@41972
   313
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
neuper@41972
   314
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   315
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   316
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   317
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   318
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   319
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   320
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   321
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   322
val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
neuper@41972
   323
(*val (p,_,f,nxt,_,pt) = me nxt p [1] pt; WAS: Empty_Tac, helpless*)
neuper@41972
   324
show_pt pt;
neuper@41972
   325
"~~~~~ fun me, args:"; val (_,tac) = nxt;
neuper@41973
   326
val (pt, p) = case locatetac tac (pt,p) of
neuper@41973
   327
	("ok", (_, _, ptp))  => ptp | _ => error "script.sml locatetac";
neuper@41972
   328
"~~~~~ fun step, args:"; val (ip as (_,p_), (ptp as (pt,p), tacis)) = (p, ((pt, e_pos'), []))
neuper@41999
   329
val pIopt = get_pblID (pt,ip); (*SOME ["linear", "univariate", "equation", "test"]*)
neuper@41972
   330
tacis; (*= []*)
neuper@41972
   331
member op = [Pbl,Met] p_ andalso is_none (get_obj g_env pt (fst p)); (*= false*)
neuper@41972
   332
"~~~~~ fun nxt_solve_, args:"; val (ptp as (pt, pos as (p, p_))) = (pt, ip);
neuper@41972
   333
                                 (*WAS stac2tac_ TODO: no match for SubProblem*)
neuper@41972
   334
val thy' = get_obj g_domID pt (par_pblobj pt p);
neuper@41972
   335
val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
neuper@48790
   336
"~~~~~ fun next_tac, args:"; val (thy, (ptp as (pt,(p,_)):ptree * pos'), (sc as Prog (h $ body)), 
neuper@41972
   337
	                               (ScrState (E,l,a,v,s,b), ctxt)) = ((thy',srls), (pt,pos), sc, is);
neuper@41972
   338
l; (*= [R, L, R, L, R, R]*)
neuper@41972
   339
val Appy (m', scrst as (_,_,_,v,_,_)) = nstep_up thy ptp sc E l Skip_ a v;
neuper@41972
   340
"~~~~~ dont like to go into nstep_up...";
neuper@41972
   341
val t = str2term ("SubProblem" ^ 
neuper@41972
   342
  "(Test', [linear, univariate, equation, test], [Test, solve_linear])" ^
neuper@41972
   343
  "[BOOL (-1 + x = 0), REAL x]");
neuper@41972
   344
val (tac, tac_) = stac2tac_ pt @{theory "Isac"} t; (*WAS stac2tac_ TODO: no match for SubProblem*)
neuper@41972
   345
case tac_ of 
neuper@41972
   346
  Subproblem' _ => ()
neuper@41972
   347
| _ => error "script.sml x+1=2 start SubProblem from script";
neuper@41972
   348
neuper@41972
   349
neuper@42283
   350
"----------- Take as 1st stac in script --------------------------";
neuper@42283
   351
"----------- Take as 1st stac in script --------------------------";
neuper@42283
   352
"----------- Take as 1st stac in script --------------------------";
neuper@42283
   353
val p = e_pos'; val c = []; 
neuper@42283
   354
val (p,_,f,nxt,_,pt) = 
neuper@42283
   355
    CalcTreeTEST 
neuper@42283
   356
        [(["functionTerm (x^^^2 + 1)", "integrateBy x", "antiDerivative FF"],
neuper@42283
   357
          ("Integrate", ["integrate","function"], ["diff","integration"]))];
neuper@42283
   358
val (p,_,f,nxt,_,pt) = me nxt p c pt; (*nxt = ("Tac ", ...) --> Add_Given...*)
neuper@42283
   359
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@42283
   360
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@42283
   361
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@42283
   362
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@42283
   363
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@42283
   364
val (p,_,f,nxt,_,pt) = me nxt p c pt;
neuper@42283
   365
case nxt of (_, Apply_Method ["diff", "integration"]) => ()
neuper@42283
   366
          | _ => error "integrate.sml -- me method [diff,integration] -- spec";
neuper@42283
   367
"----- step 8: returns nxt = Rewrite_Set_Inst ([\"(bdv, x)\"],\"integration\")";
neuper@42283
   368
neuper@42283
   369
"~~~~~ fun me, args:"; val (((_,tac):tac'_), (p:pos'), _, (pt:ptree)) = (nxt, p, c, pt);
neuper@42283
   370
"~~~~~ fun locatetac, args:"; val (tac, (ptp as (pt, p))) = (tac, (pt,p));
neuper@42283
   371
val (mI,m) = mk_tac'_ tac;
neuper@42283
   372
val Appl m = applicable_in p pt m;
neuper@42283
   373
member op = specsteps mI (*false*);
neuper@42283
   374
"~~~~~ fun loc_solve_, args:"; val (m, (pt,pos)) = ((mI,m), ptp);
neuper@42283
   375
"~~~~~ fun solve, args:"; val ((_, m as Apply_Method' (mI, _, _, _)), (pt, (pos as (p,_)))) = 
neuper@42283
   376
                            (m, (pt, pos));
neuper@42283
   377
val {srls, ...} = get_met mI;
neuper@42283
   378
        val PblObj {meth=itms, ...} = get_obj I pt p;
neuper@42283
   379
        val thy' = get_obj g_domID pt p;
neuper@42283
   380
        val thy = assoc_thy thy';
neuper@42283
   381
        val (is as ScrState (env,_,_,_,_,_), ctxt, sc) = init_scrstate thy itms mI;
neuper@42283
   382
        val ini = init_form thy sc env;
neuper@42283
   383
        val p = lev_dn p;
neuper@42283
   384
ini = NONE; (*true*)
neuper@42283
   385
            val (m', (is', ctxt'), _) = next_tac (thy', srls) (pt, (p, Res)) sc (is, ctxt);
neuper@42283
   386
	            val d = e_rls (*FIXME: get simplifier from domID*);
neuper@42283
   387
"~~~~~ fun locate_gen, args:"; val ((ts as (thy',srls)), (m:tac_), ((pt,p):ptree * pos'), 
neuper@48790
   388
	                     (scr as Prog (h $ body),d), (ScrState (E,l,a,v,S,b), ctxt)) = 
neuper@42283
   389
                   ((thy',srls), m', (pt,(p, Res)), (sc,d), (is', ctxt'));
neuper@42283
   390
val thy = assoc_thy thy';
neuper@42283
   391
l = [] orelse ((*init.in solve..Apply_Method...*)
neuper@42283
   392
			         (last_elem o fst) p = 0 andalso snd p = Res); (*true*)
neuper@42283
   393
"~~~~~ fun assy, args:"; val (ya, (is as (E,l,a,v,S,b),ss), 
neuper@42283
   394
                   (Const ("HOL.Let",_) $ e $ (Abs (id,T,body)))) = 
neuper@42283
   395
                 ((thy',ctxt,srls,d,Aundef), ((E,[R],a,v,S,b), [(m,EmptyMout,pt,p,[])]), body);
neuper@42283
   396
 (*check: true*) term2str e = "Take (Integral f_f D v_v)";
neuper@42283
   397
"~~~~~ fun assy, args:"; val ((thy',ctxt,sr,d,ap), (is as (E,l,a,v,S,b), (m,_,pt,(p,p_),c)::ss), t) = 
neuper@42283
   398
                           (ya, ((E , l@[L,R], a,v,S,b),ss), e);
neuper@42283
   399
val (a', STac stac) = handle_leaf "locate" thy' sr E a v t;
neuper@42283
   400
(*             val ctxt = get_ctxt pt (p,p_)
neuper@42283
   401
exception PTREE "get_obj: pos = [0] does not exist" raised 
neuper@42283
   402
(line 908 of "/usr/local/isabisac/src/Tools/isac/Interpret/ctree.sml")*)
neuper@42283
   403
neuper@42283
   404
neuper@42394
   405
"----------- 'trace_script' from outside 'fun me '----------------";
neuper@42394
   406
"----------- 'trace_script' from outside 'fun me '----------------";
neuper@42394
   407
"----------- 'trace_script' from outside 'fun me '----------------";
neuper@42394
   408
val fmz = ["equality (x+1=(2::real))", "solveFor x","solutions L"];
neuper@42394
   409
val (dI',pI',mI') =
neuper@42394
   410
   ("Test", ["sqroot-test","univariate","equation","test"],
neuper@42394
   411
    ["Test","squ-equ-test-subpbl1"]);
neuper@42394
   412
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
neuper@42394
   413
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@42394
   414
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@42394
   415
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@42394
   416
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@42394
   417
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@42394
   418
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
neuper@42394
   419
val (p,_,f,nxt,_,pt) = me nxt p [] pt; (*nxt = Apply_Method ["Test", "squ-equ-test-subpbl1"]*)
neuper@42394
   420
(*from_pblobj_or_detail': no istate = from_pblobj_or_detail' "Isac" p pt;*)
neuper@42394
   421
neuper@42394
   422
val (p,_,f,nxt,_,pt) = me nxt p [] pt; (*---------------------------------------------------*)
neuper@42394
   423
val (p'''', pt'''') = (p, pt);
neuper@42394
   424
f2str f = "x + 1 = 2";          snd nxt = Rewrite_Set "norm_equation";
neuper@42394
   425
val (p, p_) = p(* = ([1], Frm)*);
neuper@48790
   426
val (srls, (ScrState scrstate, ctxt), Prog sc) = from_pblobj_or_detail' "Isac" (p, p_) pt;
neuper@42394
   427
val (env, loc_, curry_arg, res, Safe, true) = scrstate;
neuper@42394
   428
  env2str env = "[\"\n(e_e, x + 1 = 2)\",\"\n(v_v, x)\"]";
neuper@42394
   429
  loc_ = [];
neuper@42394
   430
  curry_arg = NONE;
neuper@42394
   431
  term2str res = "??.empty";                     (* scrstate before starting interpretation *)
neuper@48790
   432
(*term2str (go loc_ sc) = "Prog Solve_root_equation e_e v_v = ...*)
neuper@42394
   433
neuper@42394
   434
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   435
val (p'''', pt'''') = (p, pt);
neuper@42394
   436
f2str f = "x + 1 + -1 * 2 = 0"; snd nxt = Rewrite_Set "Test_simplify";
neuper@42394
   437
val (p, p_) = p(* = ([1], Res)*);
neuper@48790
   438
val (srls, (ScrState scrstate, ctxt), Prog sc) = from_pblobj_or_detail' "Isac" (p, p_) pt;
neuper@42394
   439
val (env, loc_, curry_arg, res, Safe, true) = scrstate;
neuper@42394
   440
  env2str env = "[\"\n(e_e, x + 1 = 2)\",\"\n(v_v, x)\"]";
neuper@42394
   441
  loc_ = [R, L, R, L, L, R, R];
neuper@42394
   442
  curry_arg = SOME (str2term "e_e::bool");
neuper@42394
   443
  term2str res = "x + 1 + -1 * 2 = 0";(* scrstate after norm_equation, before Test_simplify *)
neuper@42394
   444
term2str (go loc_ sc) = "Rewrite_Set norm_equation False";
neuper@42394
   445
neuper@42394
   446
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   447
val (p'''', pt'''') = (p, pt);
neuper@42394
   448
f2str f = "-1 + x = 0"; snd nxt = Subproblem ("Test", ["linear", "univariate", "equation", "test"]);
neuper@42394
   449
val (p, p_) = p(* = ([2], Res)*);
neuper@48790
   450
val (srls, (ScrState scrstate, ctxt), Prog sc) = from_pblobj_or_detail' "Isac" (p, p_) pt;
neuper@42394
   451
val (env, loc_, curry_arg, res, Safe, false) = scrstate;
neuper@42394
   452
  env2str env = "[\"\n(e_e, x + 1 = 2)\",\"\n(v_v, x)\"]";
neuper@42394
   453
  loc_ = [R, L, R, L, R, R];
neuper@42394
   454
  curry_arg = SOME (str2term "e_e::bool");
neuper@42394
   455
  term2str res = "-1 + x = 0";           (* scrstate after Test_simplify, before Subproblem *)
neuper@42394
   456
term2str (go loc_ sc) = "Rewrite_Set Test_simplify False";
neuper@42394
   457
neuper@42394
   458
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   459
val (p'''', pt'''') = (p, pt);     (*f2str f = exception Match; ...is a model, not a formula*)
neuper@42394
   460
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   461
val (p'''', pt'''') = (p, pt);
neuper@42394
   462
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   463
val (p'''', pt'''') = (p, pt);
neuper@42394
   464
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   465
val (p'''', pt'''') = (p, pt);
neuper@42394
   466
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   467
val (p'''', pt'''') = (p, pt);
neuper@42394
   468
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   469
val (p'''', pt'''') = (p, pt);
neuper@42394
   470
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   471
val (p'''', pt'''') = (p, pt);
neuper@42394
   472
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*nxt = Apply_Method ["Test", "solve_linear"]*)
neuper@42394
   473
(*from_pblobj_or_detail': no istate = from_pblobj_or_detail' "Isac" p pt;*)
neuper@42394
   474
neuper@42394
   475
val (p,_,f,nxt,_,pt) = me nxt p [] pt; (*---------------------------------------------------*)
neuper@42394
   476
val (p'''', pt'''') = (p, pt);
neuper@42394
   477
f2str f = "-1 + x = 0"; snd nxt = Rewrite_Set_Inst (["(bdv, x)"], "isolate_bdv");
neuper@42394
   478
val (p, p_) = p(* = ([3, 1], Frm)*);
neuper@48790
   479
val (srls, (ScrState scrstate, ctxt), Prog sc) = from_pblobj_or_detail' "Isac" (p, p_) pt;
neuper@42394
   480
val (env, loc_, curry_arg, res, Safe, true) = scrstate;
neuper@42394
   481
  env2str env = "[\"\n(e_e, -1 + x = 0)\",\"\n(v_v, x)\"]";
neuper@42394
   482
  loc_ = [];
neuper@42394
   483
  curry_arg = NONE;
neuper@42394
   484
  term2str res = "??.empty";                     (* scrstate before starting interpretation *)
neuper@48790
   485
(*term2str (go loc_ sc) = "Prog Solve_linear e_e v_v = ...*)
neuper@42394
   486
neuper@42394
   487
val (p,_,f,nxt,_,pt) = me nxt p'''' [] pt''''; (*-------------------------------------------*)
neuper@42394
   488
val (p'''', pt'''') = (p, pt);
neuper@42394
   489
f2str f = "x = 0 + -1 * -1"; snd nxt = Rewrite_Set "Test_simplify";
neuper@42394
   490
val (p, p_) = p(* = ([3, 1], Res)*);
neuper@48790
   491
val (srls, (ScrState scrstate, ctxt), Prog sc) = from_pblobj_or_detail' "Isac" (p, p_) pt;
neuper@42394
   492
val (env, loc_, curry_arg, res, Safe, true) = scrstate;
neuper@42394
   493
  env2str env = "[\"\n(e_e, -1 + x = 0)\",\"\n(v_v, x)\"]";
neuper@42394
   494
  loc_ = [R, L, R, L, R, L, R];
neuper@42394
   495
  curry_arg = SOME (str2term "e_e::bool");
neuper@42394
   496
  term2str res = "x = 0 + -1 * -1";(* scrstate after isolate_bdv, before Test_simplify *)
neuper@42394
   497
term2str (go loc_ sc) = "Rewrite_Set_Inst [(bdv, v_v)] isolate_bdv False";
neuper@42394
   498
neuper@42438
   499
"----------- fun sel_rules ---------------------------------------";
neuper@42438
   500
"----------- fun sel_rules ---------------------------------------";
neuper@42438
   501
"----------- fun sel_rules ---------------------------------------";
neuper@42438
   502
(* compare test/../interface.sml
neuper@42438
   503
"--------- getTactic, fetchApplicableTactics ------------";
neuper@42438
   504
*)
neuper@42438
   505
 states:=[];
neuper@42438
   506
 CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@42438
   507
   ("Test", ["sqroot-test","univariate","equation","test"],
neuper@42438
   508
    ["Test","squ-equ-test-subpbl1"]))];
neuper@42438
   509
 Iterator 1;
neuper@42438
   510
 moveActiveRoot 1;
neuper@42438
   511
 autoCalculate 1 CompleteCalc;
neuper@42438
   512
 val ((pt,_),_) = get_calc 1;
neuper@42438
   513
 show_pt pt;
neuper@42394
   514
neuper@42438
   515
(*UC\label{SOLVE:MANUAL:TACTIC:listall} UC 30.3.2.2 p.175*)
neuper@42438
   516
 val tacs = sel_rules pt ([],Pbl);
neuper@42438
   517
 if tacs = [Apply_Method ["Test", "squ-equ-test-subpbl1"]] then ()
neuper@42438
   518
 else error "script.sml: diff.behav. in sel_rules ([],Pbl)";
neuper@42438
   519
neuper@42438
   520
 val tacs = sel_rules pt ([1],Res);
neuper@42438
   521
 if tacs = [Rewrite_Set "norm_equation", Rewrite_Set "Test_simplify",
neuper@42438
   522
      Subproblem ("Test", ["linear", "univariate", "equation", "test"]),
neuper@42438
   523
      Check_elementwise "Assumptions"] then ()
neuper@42438
   524
 else error "script.sml: diff.behav. in sel_rules ([1],Res)";
neuper@42438
   525
neuper@42438
   526
 val tacs = sel_rules pt ([3],Pbl);
neuper@42438
   527
 if tacs = [Apply_Method ["Test", "solve_linear"]] then ()
neuper@42438
   528
 else error "script.sml: diff.behav. in sel_rules ([3],Pbl)";
neuper@42438
   529
neuper@42438
   530
 val tacs = sel_rules pt ([3,1],Res);
neuper@42438
   531
 if tacs = [Rewrite_Set_Inst (["(bdv, x)"], "isolate_bdv"),
neuper@42438
   532
      Rewrite_Set "Test_simplify"] then ()
neuper@42438
   533
 else error "script.sml: diff.behav. in sel_rules ([3,1],Res)";
neuper@42438
   534
neuper@42438
   535
 val tacs = sel_rules pt ([3],Res);
neuper@42438
   536
 if tacs = [Rewrite_Set "norm_equation", Rewrite_Set "Test_simplify",
neuper@42438
   537
      Subproblem ("Test", ["linear", "univariate", "equation", "test"]),
neuper@42438
   538
      Check_elementwise "Assumptions"] then ()
neuper@42438
   539
 else error "script.sml: diff.behav. in sel_rules ([3],Res)";
neuper@42438
   540
neuper@42438
   541
 val tacs = (sel_rules pt ([],Res)) handle PTREE str => [Tac str];
neuper@42438
   542
 if tacs = [Tac "no tactics applicable at the end of a calculation"] then ()
neuper@42438
   543
 else error "script.sml: diff.behav. in sel_rules ([],Res)";
neuper@42438
   544
neuper@42438
   545
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@42438
   546
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@42438
   547
"----------- fun sel_appl_atomic_tacs ----------------------------";
neuper@42438
   548
 states:=[];
neuper@42438
   549
 CalcTree [(["equality (x+1=(2::real))", "solveFor x","solutions L"], 
neuper@42438
   550
   ("Test", ["sqroot-test","univariate","equation","test"],
neuper@42438
   551
    ["Test","squ-equ-test-subpbl1"]))];
neuper@42438
   552
 Iterator 1;
neuper@42438
   553
 moveActiveRoot 1;
neuper@42438
   554
 autoCalculate 1 CompleteCalc;
neuper@42438
   555
neuper@42438
   556
(*UC\label{SOLVE:MANUAL:TACTIC:listall} UC 30.3.2.2 p.175*)
neuper@42438
   557
 fetchApplicableTactics 1 99999 ([],Pbl);
neuper@42438
   558
neuper@42438
   559
 fetchApplicableTactics 1 99999 ([1],Res);
neuper@42438
   560
"~~~~~ fun fetchApplicableTactics, args:"; val (cI, (scope:int), (p:pos')) = (1, 99999, ([1],Res));
neuper@42438
   561
val ((pt, _), _) = get_calc cI;
neuper@42438
   562
(*version 1:*)
neuper@42438
   563
if sel_rules pt p = [Rewrite_Set "norm_equation", Rewrite_Set "Test_simplify",
neuper@42438
   564
  Subproblem ("Test", ["linear", "univariate", "equation", "test"]),
neuper@42438
   565
  Check_elementwise "Assumptions"] then ()
neuper@42438
   566
else error "fetchApplicableTactics ([1],Res) changed";
neuper@42438
   567
(*version 2:*)
neuper@42438
   568
(*WAS:
neuper@42438
   569
sel_appl_atomic_tacs pt p;
neuper@42438
   570
...
neuper@42438
   571
### atomic_appl_tacs: not impl. for tac='Subproblem(Test,["linear","univariate","equation","test"])' 
neuper@42438
   572
### atomic_appl_tacs: not impl. for tac = 'Check_elementwise "Assumptions"' 
neuper@42438
   573
*)
neuper@42438
   574
neuper@42438
   575
"~~~~~ fun sel_appl_atomic_tacs, args:"; val (pt, (p,p_)) = (pt, p);
neuper@42438
   576
is_spec_pos p_ = false;
neuper@42438
   577
        val pp = par_pblobj pt p
neuper@42438
   578
        val thy' = (get_obj g_domID pt pp):theory'
neuper@42438
   579
        val thy = assoc_thy thy'
neuper@42438
   580
        val metID = get_obj g_metID pt pp
neuper@42438
   581
        val metID' =
neuper@42438
   582
          if metID = e_metID 
neuper@42438
   583
          then (thd3 o snd3) (get_obj g_origin pt pp)
neuper@42438
   584
          else metID
neuper@48790
   585
        val {scr=Prog sc,srls,erls,rew_ord'=ro,...} = get_met metID'
neuper@42438
   586
        val ScrState (env,_,a,v,_,_) = get_istate pt (p,p_)
neuper@42438
   587
        val alltacs = (*we expect at least 1 stac in a script*)
neuper@42438
   588
          map ((stac2tac pt thy) o rep_stacexpr o #2 o
neuper@42438
   589
           (handle_leaf "selrul" thy' srls env a v)) (stacpbls sc)
neuper@42438
   590
        val f =
neuper@42438
   591
          case p_ of
neuper@42438
   592
              Frm => get_obj g_form pt p
neuper@42438
   593
            | Res => (fst o (get_obj g_result pt)) p
neuper@42438
   594
(*WN120611 stopped and took version 1 again for fetchApplicableTactics !
neuper@42438
   595
(distinct o flat o (map (atomic_appl_tacs thy ro erls f))) alltacs
neuper@42438
   596
...
neuper@42438
   597
### atomic_appl_tacs: not impl. for tac='Subproblem(Test,["linear","univariate","equation","test"])' 
neuper@42438
   598
### atomic_appl_tacs: not impl. for tac = 'Check_elementwise "Assumptions"' 
neuper@42438
   599
*)
neuper@42438
   600