src/Tools/isac/Interpret/ctree.sml
author Walther Neuper <neuper@ist.tugraz.at>
Tue, 28 Sep 2010 09:06:56 +0200
branchisac-update-Isa09-2
changeset 38031 460c24a6a6ba
parent 38015 67ba02dffacc
child 38050 4c52ad406c20
permissions -rw-r--r--
tuned error and writeln

# raise error --> error
# writeln in atomtyp, atomty, atomt and xmlsrc
neuper@37906
     1
(* use"../ME/ctree.sml";
neuper@37906
     2
   use"ME/ctree.sml";
neuper@37906
     3
   use"ctree.sml";
neuper@37906
     4
   W.N.26.10.99
neuper@37906
     5
neuper@38015
     6
tracing (pr_ptree pr_short pt); 
neuper@37906
     7
neuper@37906
     8
val Nd ( _, ns) = pt;
neuper@37906
     9
neuper@37906
    10
*)
neuper@37906
    11
neuper@37906
    12
(*structure Ptree (**): PTREE (**) = ###### outcommented ######*)
neuper@37906
    13
signature PTREE =
neuper@37906
    14
sig
neuper@37906
    15
  type ptree
neuper@37906
    16
  type envp
neuper@37906
    17
  val e_ptree : ptree
neuper@37906
    18
  exception PTREE of string
neuper@37906
    19
  type branch
neuper@37906
    20
  type ostate
neuper@37906
    21
  type cellID
neuper@37906
    22
  type cid
neuper@37906
    23
  type posel
neuper@37906
    24
  type pos
neuper@37906
    25
  type pos'
neuper@37906
    26
  type loc
neuper@37906
    27
  type domID
neuper@37906
    28
  type pblID
neuper@37906
    29
  type metID
neuper@37906
    30
  type spec
neuper@37906
    31
  type 'a ppc
neuper@37906
    32
  type con
neuper@37906
    33
  type subs
neuper@37906
    34
  type subst
neuper@37906
    35
  type env
neuper@37906
    36
  type ets
neuper@37906
    37
  val ets2str : ets -> string
neuper@37906
    38
  type item
neuper@37906
    39
  type tac
neuper@37906
    40
  type tac_
neuper@37906
    41
  val tac_2str : tac_ -> string
neuper@37906
    42
  type safe
neuper@37906
    43
  val safe2str : safe -> string
neuper@37906
    44
neuper@37906
    45
  type meth
neuper@37906
    46
  val cappend_atomic : ptree -> pos -> loc -> cterm' -> tac
neuper@37906
    47
    -> cterm' -> ostate -> cid -> ptree * posel list * cid
neuper@37906
    48
  val cappend_form : ptree
neuper@37906
    49
    -> pos -> loc -> cterm' -> cid -> ptree * pos * cid
neuper@37906
    50
  val cappend_parent : ptree -> pos -> loc -> cterm' -> tac
neuper@37906
    51
    -> branch -> cid -> ptree * int list * cid
neuper@37906
    52
  val cappend_problem : ptree -> posel list(*FIXME*) -> loc
neuper@37906
    53
    -> cterm' list * spec -> cid -> ptree * int list * cellID list
neuper@37906
    54
  val append_result : ptree -> pos -> cterm' -> ostate -> ptree * pos
neuper@37906
    55
neuper@37906
    56
  type ppobj
neuper@37906
    57
  val g_branch : ppobj -> branch
neuper@37906
    58
  val g_cell : ppobj -> cid
neuper@37906
    59
  val g_args : ppobj -> (int * (term list)) list (*args of scr*)
neuper@37906
    60
  val g_form : ppobj -> cterm'
neuper@37906
    61
  val g_loc : ppobj -> loc
neuper@37906
    62
  val g_met : ppobj -> meth
neuper@37906
    63
  val g_domID : ppobj -> domID
neuper@37906
    64
  val g_metID : ppobj -> metID
neuper@37906
    65
  val g_model : ppobj -> cterm' ppc
neuper@37906
    66
  val g_tac : ppobj -> tac
neuper@37906
    67
  val g_origin : ppobj -> cterm' list * spec
neuper@37906
    68
  val g_ostate : ppobj -> ostate
neuper@37906
    69
  val g_pbl : ppobj -> pblID * item ppc
neuper@37906
    70
  val g_result : ppobj -> cterm'
neuper@37906
    71
  val g_spec : ppobj -> spec
neuper@37906
    72
(*  val get_all : (ppobj -> 'a) -> ptree -> 'a list
neuper@37906
    73
  val get_alls : (ppobj -> 'a) -> ptree list -> 'a list *)
neuper@37906
    74
  val get_obj : (ppobj -> 'a) -> ptree -> pos -> 'a     
neuper@37906
    75
  val gpt_cell : ptree -> cid
neuper@37906
    76
  val par_pblobj : ptree -> pos -> pos
neuper@37906
    77
  val pre_pos : pos -> pos
neuper@37906
    78
  val lev_dn : int list -> int list
neuper@37906
    79
  val lev_on : pos -> posel list
neuper@37906
    80
  val lev_pred : pos -> pos
neuper@37906
    81
  val lev_up : pos -> pos
neuper@37906
    82
(*  val pr_cell : pos -> ppobj -> string
neuper@37906
    83
  val pr_pos : int list -> string        *)
neuper@37906
    84
  val pr_ptree : (pos -> ppobj -> string) -> ptree -> string
neuper@37906
    85
  val pr_short : pos -> ppobj -> string
neuper@37906
    86
(*  val repl : 'a list -> int -> 'a -> 'a list
neuper@37906
    87
  val repl_app : 'a list -> int -> 'a -> 'a list
neuper@37906
    88
  val repl_branch : branch -> ppobj -> ppobj
neuper@37906
    89
  val repl_domID : domID -> ppobj -> ppobj
neuper@37906
    90
  val repl_form : cterm' -> ppobj -> ppobj
neuper@37906
    91
  val repl_met : item ppc -> ppobj -> ppobj
neuper@37906
    92
  val repl_metID : metID -> ppobj -> ppobj
neuper@37906
    93
  val repl_model : cterm' list -> ppobj -> ppobj
neuper@37906
    94
  val repl_tac : tac -> ppobj -> ppobj
neuper@37906
    95
  val repl_pbl : item ppc -> ppobj -> ppobj
neuper@37906
    96
  val repl_pblID : pblID -> ppobj -> ppobj
neuper@37906
    97
  val repl_result : cterm' -> ostate -> ppobj -> ppobj
neuper@37906
    98
  val repl_spec : spec -> ppobj -> ppobj
neuper@37906
    99
  val repl_subs : (string * string) list -> ppobj -> ppobj     *)
neuper@37906
   100
  val rootthy : ptree -> domID
neuper@37906
   101
(*  val test_trans : ppobj -> bool
neuper@37906
   102
  val uni__asm : (string * pos) list -> ppobj -> ppobj
neuper@37906
   103
  val uni__cid : cellID list -> ppobj -> ppobj                 *)
neuper@37906
   104
  val union_asm : ptree -> pos -> (string * pos) list -> ptree
neuper@37906
   105
  val union_cid : ptree -> pos -> cellID list -> ptree
neuper@37906
   106
  val update_branch : ptree -> pos -> branch -> ptree
neuper@37906
   107
  val update_domID : ptree -> pos -> domID -> ptree
neuper@37906
   108
  val update_met : ptree -> pos -> meth -> ptree
neuper@37906
   109
  val update_metppc : ptree -> pos -> item ppc -> ptree
neuper@37906
   110
  val update_metID : ptree -> pos -> metID -> ptree
neuper@37906
   111
  val update_tac : ptree -> pos -> tac -> ptree
neuper@37906
   112
  val update_pbl : ptree -> pos -> pblID * item ppc -> ptree
neuper@37906
   113
  val update_pblppc : ptree -> pos -> item ppc -> ptree
neuper@37906
   114
  val update_pblID : ptree -> pos -> pblID -> ptree
neuper@37906
   115
  val update_spec : ptree -> pos -> spec -> ptree
neuper@37906
   116
  val update_subs : ptree -> pos -> (string * string) list -> ptree
neuper@37906
   117
neuper@37906
   118
  val rep_pblobj : ppobj
neuper@37906
   119
    -> {branch:branch, cell:cid, env:envp, loc:loc, meth:meth, model:cterm' ppc,
neuper@37906
   120
        origin:cterm' list * spec, ostate:ostate, probl:pblID * item ppc,
neuper@37906
   121
        result:cterm', spec:spec}
neuper@37906
   122
  val rep_prfobj : ppobj
neuper@37906
   123
    -> {branch:branch, cell:cid, form:cterm', loc:loc, tac:tac,
neuper@37906
   124
        ostate:ostate, result:cterm'}
neuper@37906
   125
end 
neuper@37906
   126
neuper@37906
   127
(* -------------- 
neuper@37906
   128
structure Ptree (**): PTREE (**) =
neuper@37906
   129
struct
neuper@37906
   130
 -------------- *)
neuper@37906
   131
neuper@37906
   132
type env = (term * term) list;
neuper@37906
   133
neuper@37906
   134
   
neuper@37906
   135
datatype branch = 
neuper@37906
   136
	 NoBranch | AndB | OrB 
neuper@37906
   137
       | TransitiveB  (* FIXXXME.8.03: set branch from met in Apply_Method
neuper@37906
   138
                         FIXXXME.0402: -"- in Begin_Trans'*)
neuper@37906
   139
       | SequenceB | IntersectB | CollectB | MapB;
neuper@37906
   140
fun branch2str NoBranch = "NoBranch"
neuper@37906
   141
  | branch2str AndB = "AndB"
neuper@37906
   142
  | branch2str OrB = "OrB"
neuper@37906
   143
  | branch2str TransitiveB = "TransitiveB" 
neuper@37906
   144
  | branch2str SequenceB = "SequenceB"
neuper@37906
   145
  | branch2str IntersectB = "IntersectB"
neuper@37906
   146
  | branch2str CollectB = "CollectB"
neuper@37906
   147
  | branch2str MapB = "MapB";
neuper@37906
   148
neuper@37906
   149
datatype ostate = 
neuper@37906
   150
    Incomplete | Complete | Inconsistent(*WN041020 latter unused*);
neuper@37906
   151
fun ostate2str Incomplete = "Incomplete"
neuper@37906
   152
  | ostate2str Complete = "Complete"
neuper@37906
   153
  | ostate2str Inconsistent = "Inconsistent";
neuper@37906
   154
neuper@37906
   155
type cellID = int;     
neuper@37906
   156
type cid = cellID list;
neuper@37906
   157
neuper@37906
   158
type posel = int;     (* roundabout for (some of) nice signatures *)
neuper@37906
   159
type pos = posel list;
neuper@37906
   160
val pos2str = ints2str';
neuper@37906
   161
datatype pos_ = 
neuper@37906
   162
    Pbl    (*PblObj-position: problem-type*)
neuper@37906
   163
  | Met    (*PblObj-position: method*)
neuper@37906
   164
  | Frm    (*PblObj-position: -> Pbl in ME (not by moveDown !)
neuper@37906
   165
           | PrfObj-position: formula*)
neuper@37906
   166
  | Res    (*PblObj | PrfObj-position: result*)
neuper@37906
   167
  | Und;   (*undefined*)
neuper@37906
   168
fun pos_2str Pbl = "Pbl"
neuper@37906
   169
  | pos_2str Met = "Met"
neuper@37906
   170
  | pos_2str Frm = "Frm"
neuper@37906
   171
  | pos_2str Res = "Res"
neuper@37906
   172
  | pos_2str Und = "Und";
neuper@37906
   173
neuper@37906
   174
type pos' = pos * pos_;
neuper@37906
   175
(*WN.12.03 remembering interator (pos * pos_) for ptree 
neuper@37906
   176
	   pos : lev_on, lev_dn, lev_up, 
neuper@37906
   177
                 lev_onFrm, lev_dnRes (..see solve Apply_Method !) 
neuper@37906
   178
           pos_:
neuper@37906
   179
# generate1 sets pos_ if possible  ...?WN0502?NOT...
neuper@37906
   180
# generate1 does NOT set pos, because certain nodes can be lev_on OR lev_dn
neuper@37906
   181
                     exceptions: Begin/End_Trans
neuper@37906
   182
# thus generate(1) called in
neuper@37906
   183
.# assy, locate_gen 
neuper@37906
   184
.# nxt_solv (tac_ -cases); general case: 
neuper@37906
   185
  val pos' = case pos' of (p,Res) => (lev_on p',Res) | _ => pos'
neuper@37906
   186
# WN050220, S(604):
neuper@37906
   187
  generate1...(Rewrite(f,..,res))..(pos, pos_)
neuper@37906
   188
     cappend_atomic.................pos //////  gets f+res always!!!
neuper@37906
   189
        cut_tree....................pos, pos_ 
neuper@37906
   190
*)
neuper@37906
   191
fun pos'2str (p,p_) = pair2str (ints2str' p, pos_2str p_);
neuper@37906
   192
fun pos's2str ps = (strs2str' o (map pos'2str)) ps;
neuper@37906
   193
val e_pos' = ([],Und):pos';
neuper@37906
   194
neuper@37906
   195
fun res2str (t, ts) = pair2str (term2str t, terms2str ts);
neuper@37906
   196
fun asm2str (t, p:pos) = pair2str (term2str t, ints2str' p);
neuper@37906
   197
fun asms2str asms = (strs2str' o (map asm2str)) asms;
neuper@37906
   198
neuper@37906
   199
neuper@37906
   200
neuper@37906
   201
(*26.4.02: never used after introduction of scripts !!!
neuper@37906
   202
type loc =  loc_ *        (* + interpreter-state          *)
neuper@37906
   203
	    (loc_ * rls') (* -"- for script of the ruleset*)
neuper@37906
   204
		option;
neuper@37926
   205
val e_loc = ([],NONE):loc;
neuper@37906
   206
val ee_loc = (e_loc,e_loc);*)
neuper@37906
   207
neuper@37906
   208
neuper@37906
   209
datatype safe = Sundef | Safe | Unsafe | Helpless;
neuper@37906
   210
fun safe2str Sundef   = "Sundef"
neuper@37906
   211
  | safe2str Safe     = "Safe"
neuper@37906
   212
  | safe2str Unsafe   = "Unsafe" 
neuper@37906
   213
  | safe2str Helpless = "Helpless";
neuper@37906
   214
neuper@37906
   215
type subs = cterm' list; (*16.11.00 for FE-KE*)
neuper@37906
   216
val e_subs = ["(bdv, x)"];
neuper@37906
   217
neuper@37906
   218
(*._sub_stitution as strings of _e_qualities.*)
neuper@37906
   219
type sube = cterm' list;
neuper@37906
   220
val e_sube = []:cterm' list;
neuper@37906
   221
fun sube2str s = strs2str s;
neuper@37906
   222
neuper@37906
   223
(*._sub_stitution as _t_erms of _e_qualities.*)
neuper@37906
   224
type subte = term list;
neuper@37906
   225
val e_subte = []:term list;
neuper@37906
   226
fun subte2str ss = terms2str ss;
neuper@37906
   227
neuper@37906
   228
fun subte2sube ss = map term2str ss;
neuper@37906
   229
neuper@37906
   230
fun subst2subs s = map (pair2str o 
neuper@37929
   231
			(apfst (Syntax.string_of_term (thy2ctxt' "Isac"))) o
neuper@37929
   232
			(apsnd (Syntax.string_of_term (thy2ctxt' "Isac")))) s;
neuper@37929
   233
fun subst2subs' s = map ((apfst (Syntax.string_of_term (thy2ctxt' "Isac"))) o
neuper@37929
   234
			 (apsnd (Syntax.string_of_term (thy2ctxt' "Isac")))) s;
neuper@37906
   235
fun subs2subst thy s = map (isapair2pair o term_of o the o (parse thy)) s;
neuper@37906
   236
(*> subs2subst thy ["(bdv,x)","(err,#0)"];
neuper@37906
   237
val it =
neuper@37906
   238
  [(Free ("bdv","RealDef.real"),Free ("x","RealDef.real")),
neuper@37906
   239
   (Free ("err","RealDef.real"),Free ("#0","RealDef.real"))] 
neuper@37906
   240
   : (term * term) list*)
neuper@37906
   241
(*["bdv=x","err=0"] ---> [(bdv,x), (err,0)]*)
neuper@37906
   242
fun sube2subst thy s = map (dest_equals' o term_of o the o (parse thy)) s;
neuper@37906
   243
(* val ts = sube2subst thy ["bdv=x","err=0"];
neuper@37906
   244
   subst2str' ts;
neuper@37906
   245
   *)
neuper@37906
   246
fun sube2subte ss = map str2term ss;
neuper@37906
   247
neuper@37906
   248
neuper@37906
   249
fun isasub2subst isasub = ((map isapair2pair) o isalist2list) isasub;
neuper@37906
   250
neuper@37906
   251
neuper@37906
   252
type scrstate =       (*state for script interpreter*)
neuper@37906
   253
	 env(*stack*) (*used to instantiate tac for checking assod
neuper@37906
   254
		       12.03.noticed: e_ not updated during execution ?!?*)
neuper@37906
   255
	 * loc_       (*location of tac in script*)
neuper@37906
   256
	 * term option(*argument of curried functions*)
neuper@37906
   257
	 * term       (*value obtained by tac executed
neuper@37906
   258
		       updated also after a derivation by 'new_val'*)
neuper@37906
   259
	 * safe       (*estimation of how result will be obtained*)
neuper@37906
   260
	 * bool;      (*true = strongly .., false = weakly associated: 
neuper@37906
   261
					    only used during ass_dn/up*)
neuper@37925
   262
val e_scrstate = ([],[],NONE,e_term,Sundef,false):scrstate;
neuper@37906
   263
neuper@37906
   264
neuper@37906
   265
(*21.8.02 ---> definitions.sml for datatype scr 
neuper@37906
   266
type rrlsstate =      (*state for reverse rewriting*)
neuper@37906
   267
     (term *          (*the current formula*)
neuper@37906
   268
      rule list      (*of reverse rewrite set (#1#)*)
neuper@37906
   269
	    list *    (*may be serveral, eg. in norm_rational*)
neuper@37906
   270
      (rule *         (*Thm (+ Thm generated from Calc) resulting in ...*)
neuper@37906
   271
       (term *        (*... rewrite with ...*)
neuper@37906
   272
	term list))   (*... assumptions*)
neuper@37906
   273
	  list);      (*derivation from given term to normalform
neuper@37906
   274
		       in reverse order with sym_thm; 
neuper@37906
   275
                       (#1#) could be extracted from here #1*) --------*)
neuper@37906
   276
     
neuper@37906
   277
datatype istate =     (*interpreter state*)
neuper@37906
   278
	 Uistate                 (*undefined in modspec, in '_deriv'ation*)
neuper@37906
   279
       | ScrState of scrstate    (*for script interpreter*)
neuper@37906
   280
       | RrlsState of rrlsstate; (*for reverse rewriting*)
neuper@37925
   281
val e_istate = (ScrState ([],[],NONE,e_term,Sundef,false)):istate;
neuper@37906
   282
neuper@37906
   283
type iist = istate option * istate option;
neuper@37906
   284
(*val e_iist = (e_istate, e_istate); --- sinnlos f"ur NICHT-equality-type*) 
neuper@37906
   285
neuper@37906
   286
neuper@37906
   287
fun rta2str (r,(t,a)) = "\n("^(rule2str r)^",("^(term2str t)^", "^
neuper@37906
   288
		      (terms2str a)^"))";
neuper@37906
   289
fun istate2str Uistate = "Uistate"
neuper@37906
   290
  | istate2str (ScrState (e,l,to,t,s,b):istate) =
neuper@37906
   291
    "ScrState ("^ subst2str e ^",\n "^ 
neuper@37906
   292
    loc_2str l ^", "^ termopt2str to ^",\n "^
neuper@37906
   293
    term2str t ^", "^ safe2str s ^", "^ bool2str b ^")"
neuper@37906
   294
  | istate2str (RrlsState (t,t1,rss,rtas)) = 
neuper@37906
   295
    "RrlsState ("^(term2str t)^", "^(term2str t1)^", "^
neuper@37906
   296
    ((strs2str o (map (strs2str o (map rule2str)))) rss)^", "^
neuper@37906
   297
    ((strs2str o (map rta2str)) rtas)^")";
neuper@37925
   298
fun istates2str (NONE, NONE) = "(#NONE, #NONE)"
neuper@37925
   299
  | istates2str (NONE, SOME ist) = "(#NONE,\n#SOME "^istate2str ist^")"
neuper@37925
   300
  | istates2str (SOME ist, NONE) = "(#SOME "^istate2str ist^",\n #NONE)"
neuper@37925
   301
  | istates2str (SOME i1, SOME i2) = "(#SOME "^istate2str i1^",\n #SOME "^
neuper@37906
   302
				     istate2str i2^")";
neuper@37906
   303
neuper@37906
   304
fun new_val v (ScrState (env, loc_, topt, _, safe, bool)) =
neuper@37906
   305
    (ScrState (env, loc_, topt, v, safe, bool))
neuper@38031
   306
  | new_val _ _ = error "new_val: only for ScrState";
neuper@37906
   307
neuper@37906
   308
datatype con = land | lor;
neuper@37906
   309
neuper@37906
   310
neuper@37906
   311
type spec = 
neuper@37906
   312
     domID * (*WN.12.03: is replaced by thy from get_met ?FIXME? in:
neuper@37906
   313
	      specify (Init_Proof..), nxt_specify_init_calc,
neuper@37906
   314
	      assod (.SubProblem...), stac2tac (.SubProblem...)*)
neuper@37906
   315
     pblID * 
neuper@37906
   316
     metID;
neuper@37906
   317
fun spec2str ((dom,pbl,met)(*:spec*)) = 
neuper@37906
   318
  "(" ^ (quote dom) ^ ", " ^ (strs2str pbl) ^ 
neuper@37906
   319
  ", " ^ (strs2str met) ^ ")";
neuper@37906
   320
(*> spec2str empty_spec;
neuper@37906
   321
val it = "(\"\", [], (\"\", \"\"))" : string *)
neuper@37906
   322
val empty_spec = (e_domID,e_pblID,e_metID):spec;
neuper@37906
   323
val e_spec = empty_spec;
neuper@37906
   324
neuper@37906
   325
neuper@37906
   326
neuper@37906
   327
(*.tactics propagate the construction of the calc-tree;
neuper@37906
   328
   there are
neuper@37906
   329
   (a) 'specsteps' for the specify-phase, and others for the solve-phase
neuper@37906
   330
   (b) those of the solve-phase are 'initac's and others;
neuper@37906
   331
       initacs start with a formula different from the preceding formula.
neuper@37906
   332
   see 'type tac_' for the internal representation of tactics.*)
neuper@37906
   333
datatype tac = 
neuper@37906
   334
  Init_Proof of ((cterm' list) * spec)
neuper@37906
   335
(*'specsteps'...*)
neuper@37906
   336
| Model_Problem
neuper@37906
   337
| Refine_Problem of pblID              | Refine_Tacitly of pblID
neuper@37906
   338
neuper@37906
   339
| Add_Given of cterm'                  | Del_Given of cterm'
neuper@37906
   340
| Add_Find of cterm'                   | Del_Find of cterm'
neuper@37906
   341
| Add_Relation of cterm'               | Del_Relation of cterm'
neuper@37906
   342
neuper@37906
   343
| Specify_Theory of domID              | Specify_Problem of pblID
neuper@37906
   344
| Specify_Method of metID
neuper@37906
   345
(*...'specsteps'*)
neuper@37906
   346
| Apply_Method of metID 
neuper@37906
   347
(*.creates an 'istate' in PblObj.env; in case of 'init_form' 
neuper@37906
   348
   creates a formula at ((lev_on o lev_dn) p, Frm) and in this ppobj.'loc' 
neuper@37925
   349
   'SOME istate' (at fst of 'loc').
neuper@37906
   350
   As each step (in the solve-phase) has a resulting formula (at the front-end)
neuper@37906
   351
   Apply_Method also does the 1st step in the script (an 'initac') if there
neuper@37906
   352
   is no 'init_form' .*)
neuper@37906
   353
| Check_Postcond of pblID
neuper@37906
   354
| Free_Solve
neuper@37906
   355
neuper@37906
   356
| Rewrite_Inst of ( subs * thm')       | Rewrite of thm'
neuper@37906
   357
                                       | Rewrite_Asm of thm'
neuper@37906
   358
| Rewrite_Set_Inst of ( subs * rls')   | Rewrite_Set of rls'        
neuper@37906
   359
| Detail_Set_Inst of ( subs * rls')    | Detail_Set of rls'
neuper@37906
   360
| End_Detail  (*end of script from next_tac, 
neuper@37906
   361
                in solve: switches back to parent script WN0509 drop!*)
neuper@37906
   362
| Derive of rls' (*an input formula using rls WN0509 drop!*)
neuper@37906
   363
| Calculate of string (* plus | minus | times | cancel | pow | sqrt *)
neuper@37906
   364
| End_Ruleset
neuper@37906
   365
| Substitute of sube                   | Apply_Assumption of cterm' list
neuper@37906
   366
neuper@37906
   367
| Take of cterm'      (*an 'initac'*)
neuper@37906
   368
| Take_Inst of cterm'  
neuper@37906
   369
| Group of (con * int list ) 
neuper@37906
   370
| Subproblem of (domID * pblID) (*an 'initac'*)
neuper@37906
   371
| CAScmd of cterm'  (*6.6.02 URD: Function formula; WN0509 drop!*)
neuper@37906
   372
| End_Subproblem    (*WN0509 drop!*)
neuper@37906
   373
neuper@37906
   374
| Split_And                            | Conclude_And
neuper@37906
   375
| Split_Or                             | Conclude_Or
neuper@37906
   376
| Begin_Trans                          | End_Trans
neuper@37906
   377
| Begin_Sequ                           | End_Sequ(* substitute root.env *)
neuper@37906
   378
| Split_Intersect                      | End_Intersect
neuper@37906
   379
| Check_elementwise of cterm'          | Collect_Trues
neuper@37906
   380
| Or_to_List
neuper@37906
   381
neuper@37906
   382
| Empty_Tac (*TODO.11.6.03 ... of string: could carry msg of (Notappl msg)
neuper@37906
   383
	       in 'helpless'*)
neuper@37906
   384
| Tac of string(* eg.'repeat'*WN0509 drop!*)
neuper@37906
   385
| User                                 (*internal, for ets*WN0509 drop!*)
neuper@37906
   386
| End_Proof';(* inout*)
neuper@37906
   387
neuper@37906
   388
(* tac2str /--> library.sml: needed in dialog.sml for 'separable *)
neuper@37906
   389
fun tac2str (ma:tac) = case ma of
neuper@37906
   390
    Init_Proof (ppc, spec)  => 
neuper@37906
   391
      "Init_Proof "^(pair2str (strs2str ppc, spec2str spec))
neuper@37906
   392
  | Model_Problem           => "Model_Problem "
neuper@37906
   393
  | Refine_Tacitly pblID    => "Refine_Tacitly "^(strs2str pblID)
neuper@37906
   394
  | Refine_Problem pblID    => "Refine_Problem "^(strs2str pblID)
neuper@37906
   395
  | Add_Given cterm'        => "Add_Given "^cterm'
neuper@37906
   396
  | Del_Given cterm'        => "Del_Given "^cterm'
neuper@37906
   397
  | Add_Find cterm'         => "Add_Find "^cterm'
neuper@37906
   398
  | Del_Find cterm'         => "Del_Find "^cterm'
neuper@37906
   399
  | Add_Relation cterm'     => "Add_Relation "^cterm'
neuper@37906
   400
  | Del_Relation cterm'     => "Del_Relation "^cterm'
neuper@37906
   401
neuper@37906
   402
  | Specify_Theory domID    => "Specify_Theory "^(quote domID    )
neuper@37906
   403
  | Specify_Problem pblID   => "Specify_Problem "^(strs2str pblID )
neuper@37906
   404
  | Specify_Method metID    => "Specify_Method "^(strs2str metID)
neuper@37906
   405
  | Apply_Method metID      => "Apply_Method "^(strs2str metID)
neuper@37906
   406
  | Check_Postcond pblID    => "Check_Postcond "^(strs2str pblID)
neuper@37906
   407
  | Free_Solve              => "Free_Solve"
neuper@37906
   408
neuper@37906
   409
  | Rewrite_Inst (subs,thm')=> 
neuper@37906
   410
      "Rewrite_Inst "^(pair2str (subs2str subs, spair2str thm'))
neuper@37906
   411
  | Rewrite thm'            => "Rewrite "^(spair2str thm')
neuper@37906
   412
  | Rewrite_Asm thm'        => "Rewrite_Asm "^(spair2str thm')
neuper@37906
   413
  | Rewrite_Set_Inst (subs, rls) => 
neuper@37906
   414
      "Rewrite_Set_Inst "^(pair2str (subs2str subs, quote rls))
neuper@37906
   415
  | Rewrite_Set rls         => "Rewrite_Set "^(quote rls    )
neuper@37906
   416
  | Detail_Set rls          => "Detail_Set "^(quote rls    )
neuper@37906
   417
  | Detail_Set_Inst (subs, rls) => 
neuper@37906
   418
      "Detail_Set_Inst "^(pair2str (subs2str subs, quote rls))
neuper@37906
   419
  | End_Detail              => "End_Detail"
neuper@37906
   420
  | Derive rls'             => "Derive "^rls' 
neuper@37906
   421
  | Calculate op_           => "Calculate "^op_ 
neuper@37906
   422
  | Substitute sube         => "Substitute "^sube2str sube	     
neuper@37906
   423
  | Apply_Assumption ct's   => "Apply_Assumption "^(strs2str ct's)
neuper@37906
   424
neuper@37906
   425
  | Take cterm'             => "Take "^(quote cterm'	)
neuper@37906
   426
  | Take_Inst cterm'        => "Take_Inst "^(quote cterm' )
neuper@37906
   427
  | Group (con, ints)       => 
neuper@37906
   428
      "Group "^(pair2str (con2str con, ints2str ints))
neuper@37906
   429
  | Subproblem (domID, pblID) => 
neuper@37906
   430
      "Subproblem "^(pair2str (domID, strs2str pblID))
neuper@37906
   431
(*| Subproblem_Full (spec, cts') => 
neuper@37906
   432
      "Subproblem_Full "^(pair2str (spec2str spec, strs2str cts'))*)
neuper@37906
   433
  | End_Subproblem          => "End_Subproblem"
neuper@37906
   434
  | CAScmd cterm'           => "CAScmd "^(quote cterm')
neuper@37906
   435
neuper@37906
   436
  | Check_elementwise cterm'=> "Check_elementwise "^(quote cterm') 
neuper@37906
   437
  | Or_to_List              => "Or_to_List "
neuper@37906
   438
  | Collect_Trues           => "Collect_Trues"
neuper@37906
   439
neuper@37906
   440
  | Empty_Tac             => "Empty_Tac"
neuper@37906
   441
  | Tac string            => "Tac "^string
neuper@37906
   442
  | User                    => "User"
neuper@37906
   443
  | End_Proof'              => "tac End_Proof'"
neuper@37906
   444
  | _                       => "tac2str not impl. for ?!";
neuper@37906
   445
neuper@37906
   446
fun is_rewset (Rewrite_Set_Inst _) = true
neuper@37906
   447
  | is_rewset (Rewrite_Set _) = true 
neuper@37906
   448
  | is_rewset _ = false;
neuper@37906
   449
fun is_rewtac (Rewrite _) = true
neuper@37906
   450
  | is_rewtac (Rewrite_Inst _) = true
neuper@37906
   451
  | is_rewtac (Rewrite_Asm _) = true
neuper@37906
   452
  | is_rewtac tac = is_rewset tac;
neuper@37906
   453
neuper@37906
   454
fun tac2IDstr (ma:tac) = case ma of
neuper@37906
   455
    Model_Problem           => "Model_Problem"
neuper@37906
   456
  | Refine_Tacitly pblID    => "Refine_Tacitly"
neuper@37906
   457
  | Refine_Problem pblID    => "Refine_Problem"
neuper@37906
   458
  | Add_Given cterm'        => "Add_Given"
neuper@37906
   459
  | Del_Given cterm'        => "Del_Given"
neuper@37906
   460
  | Add_Find cterm'         => "Add_Find"
neuper@37906
   461
  | Del_Find cterm'         => "Del_Find"
neuper@37906
   462
  | Add_Relation cterm'     => "Add_Relation"
neuper@37906
   463
  | Del_Relation cterm'     => "Del_Relation"
neuper@37906
   464
neuper@37906
   465
  | Specify_Theory domID    => "Specify_Theory"
neuper@37906
   466
  | Specify_Problem pblID   => "Specify_Problem"
neuper@37906
   467
  | Specify_Method metID    => "Specify_Method"
neuper@37906
   468
  | Apply_Method metID      => "Apply_Method"
neuper@37906
   469
  | Check_Postcond pblID    => "Check_Postcond"
neuper@37906
   470
  | Free_Solve              => "Free_Solve"
neuper@37906
   471
neuper@37906
   472
  | Rewrite_Inst (subs,thm')=> "Rewrite_Inst"
neuper@37906
   473
  | Rewrite thm'            => "Rewrite"
neuper@37906
   474
  | Rewrite_Asm thm'        => "Rewrite_Asm"
neuper@37906
   475
  | Rewrite_Set_Inst (subs, rls) => "Rewrite_Set_Inst"
neuper@37906
   476
  | Rewrite_Set rls         => "Rewrite_Set"
neuper@37906
   477
  | Detail_Set rls          => "Detail_Set"
neuper@37906
   478
  | Detail_Set_Inst (subs, rls) => "Detail_Set_Inst"
neuper@37906
   479
  | Derive rls'             => "Derive "
neuper@37906
   480
  | Calculate op_           => "Calculate "
neuper@37906
   481
  | Substitute subs         => "Substitute" 
neuper@37906
   482
  | Apply_Assumption ct's   => "Apply_Assumption"
neuper@37906
   483
neuper@37906
   484
  | Take cterm'             => "Take"
neuper@37906
   485
  | Take_Inst cterm'        => "Take_Inst"
neuper@37906
   486
  | Group (con, ints)       => "Group"
neuper@37906
   487
  | Subproblem (domID, pblID) => "Subproblem"
neuper@37906
   488
  | End_Subproblem          => "End_Subproblem"
neuper@37906
   489
  | CAScmd cterm'           => "CAScmd"
neuper@37906
   490
neuper@37906
   491
  | Check_elementwise cterm'=> "Check_elementwise"
neuper@37906
   492
  | Or_to_List              => "Or_to_List "
neuper@37906
   493
  | Collect_Trues           => "Collect_Trues"
neuper@37906
   494
neuper@37906
   495
  | Empty_Tac             => "Empty_Tac"
neuper@37906
   496
  | Tac string            => "Tac "
neuper@37906
   497
  | User                    => "User"
neuper@37906
   498
  | End_Proof'              => "End_Proof'"
neuper@37906
   499
  | _                       => "tac2str not impl. for ?!";
neuper@37906
   500
neuper@37906
   501
fun rls_of (Rewrite_Set_Inst (_, rls)) = rls
neuper@37906
   502
  | rls_of (Rewrite_Set rls) = rls
neuper@38031
   503
  | rls_of tac = error ("rls_of: called with tac '"^tac2IDstr tac^"'");
neuper@37906
   504
neuper@37906
   505
fun thm_of_rew (Rewrite_Inst (subs,(thmID,_))) = 
neuper@37925
   506
    (thmID, SOME ((subs2subst (assoc_thy "Isac.thy") subs):subst))
neuper@37925
   507
  | thm_of_rew (Rewrite  (thmID,_)) = (thmID, NONE)
neuper@37925
   508
  | thm_of_rew (Rewrite_Asm (thmID,_)) = (thmID, NONE);
neuper@37906
   509
neuper@37906
   510
fun rls_of_rewset (Rewrite_Set_Inst (subs,rls)) = 
neuper@37925
   511
    (rls, SOME ((subs2subst (assoc_thy "Isac.thy") subs):subst))
neuper@37925
   512
  | rls_of_rewset (Rewrite_Set rls) = (rls, NONE)
neuper@37925
   513
  | rls_of_rewset (Detail_Set rls) = (rls, NONE)
neuper@37906
   514
  | rls_of_rewset (Detail_Set_Inst (subs, rls)) = 
neuper@37925
   515
    (rls, SOME ((subs2subst (assoc_thy "Isac.thy") subs):subst));
neuper@37906
   516
neuper@37906
   517
fun rule2tac _ (Calc (opID, thm)) = Calculate (calID2calcID opID)
neuper@37906
   518
  | rule2tac [] (Thm (thmID, thm)) = Rewrite (thmID, string_of_thmI thm)
neuper@37906
   519
  | rule2tac subst (Thm (thmID, thm)) = 
neuper@37906
   520
    Rewrite_Inst (subst2subs subst, (thmID, string_of_thmI thm))
neuper@37906
   521
  | rule2tac [] (Rls_ rls) = Rewrite_Set (id_rls rls)
neuper@37906
   522
  | rule2tac subst (Rls_ rls) = 
neuper@37906
   523
    Rewrite_Set_Inst (subst2subs subst, (id_rls rls))
neuper@37906
   524
  | rule2tac _ rule = 
neuper@38031
   525
    error ("rule2tac: called with '" ^ rule2str rule ^ "'");
neuper@37906
   526
neuper@37906
   527
type fmz_ = cterm' list;
neuper@37906
   528
neuper@37906
   529
(*.a formalization of an example containing data 
neuper@37906
   530
   sufficient for mechanically finding the solution for the example.*)
neuper@37906
   531
(*FIXME.WN051014: dont store fmz = (_,spec) in the PblObj, 
neuper@37906
   532
  this is done in origin*)
neuper@37906
   533
type fmz = fmz_ * spec;
neuper@37906
   534
val e_fmz = ([],e_spec);
neuper@37906
   535
neuper@37906
   536
(*tac_ is made from tac in applicable_in,
neuper@37906
   537
  and carries all data necessary for generate;*)
neuper@37906
   538
datatype tac_ = 
neuper@37906
   539
(* datatype tac = *)
neuper@37906
   540
  Init_Proof' of ((cterm' list) * spec)
neuper@37906
   541
                (* ori list !: code specify -> applicable*)
neuper@37906
   542
| Model_Problem' of pblID * 
neuper@37906
   543
		    itm list *  (*the 'untouched' pbl*)
neuper@37906
   544
		    itm list    (*the casually completed met*)
neuper@37906
   545
| Refine_Tacitly' of pblID *    (*input*)
neuper@37906
   546
		     pblID *    (*the refined from applicable_in*)
neuper@37906
   547
		     domID *    (*from new pbt?! filled in specify*)
neuper@37906
   548
		     metID *    (*from new pbt?! filled in specify*)
neuper@37906
   549
		     itm list   (*drop ! 9.03: remains [] for
neuper@37906
   550
                                  Model_Problem recognizing its activation*)
neuper@37906
   551
| Refine_Problem' of (pblID * (itm list * (bool * Term.term) list))
neuper@37906
   552
 (*FIXME?040215 drop: done automatically in init_proof + Subproblem'*)
neuper@37906
   553
| Add_Given'    of cterm' *
neuper@37906
   554
		   itm list (*updated with input in fun specify_additem*)
neuper@37906
   555
| Add_Find'     of cterm' *
neuper@37906
   556
		   itm list (*updated with input in fun specify_additem*)
neuper@37906
   557
| Add_Relation' of cterm' *
neuper@37906
   558
		 itm list (*updated with input in fun specify_additem*)
neuper@37906
   559
| Del_Given' of cterm'   | Del_Find' of cterm'   | Del_Relation' of cterm'
neuper@37906
   560
  (*4.00.: all..    term: in applicable_in ..? Syn ?only for FormFK?*)
neuper@37906
   561
neuper@37906
   562
| Specify_Theory' of domID              
neuper@37906
   563
| Specify_Problem' of (pblID *        (*               *)
neuper@37906
   564
		       (bool *        (* matches	     *)
neuper@37906
   565
			(itm list *   (* ppc	     *)
neuper@37906
   566
			 (bool * term) list))) (* preconditions *)
neuper@37906
   567
| Specify_Method' of metID *
neuper@37906
   568
		     ori list * (*repl. "#undef"*)
neuper@37906
   569
		     itm list   (*... updated from pbl to met*)
neuper@37906
   570
| Apply_Method' of metID * 
neuper@37906
   571
		   (term option) * (*init_form*)
neuper@37906
   572
		   istate		        
neuper@37906
   573
| Check_Postcond' of 
neuper@37906
   574
  pblID * 
neuper@37906
   575
  (term *      (*returnvalue of script in solve*)
neuper@37906
   576
   cterm' list)(*collect by get_assumptions_ in applicable_in, except if 
neuper@37906
   577
                 butlast tac is Check_elementwise: take only these asms*)
neuper@37906
   578
| Free_Solve'
neuper@37906
   579
neuper@37906
   580
| Rewrite_Inst' of theory' * rew_ord' * rls
neuper@37906
   581
		   * bool * subst * thm' * term * (term  * term list)
neuper@37906
   582
| Rewrite' of theory' * rew_ord' * rls * bool * thm' * 
neuper@37906
   583
	      term * (term * term list)
neuper@37906
   584
| Rewrite_Asm' of theory' * rew_ord' * rls * bool * thm' * 
neuper@37906
   585
  term * (term * term list)
neuper@37906
   586
| Rewrite_Set_Inst' of theory' * bool * subst * rls * 
neuper@37906
   587
		       term * (term * term list)
neuper@37906
   588
| Detail_Set_Inst' of theory' * bool * subst * rls * 
neuper@37906
   589
		      term * (term * term list)
neuper@37906
   590
| Rewrite_Set' of theory' * bool * rls * term * (term * term list)
neuper@37906
   591
| Detail_Set' of theory' * bool * rls * term * (term * term list)
neuper@37906
   592
| End_Detail' of (term * (term list)) (*see End_Trans'*)
neuper@37906
   593
| End_Ruleset' of term
neuper@37906
   594
| Derive' of rls
neuper@37906
   595
| Calculate' of theory' * string * term * (term * thm') 
neuper@37906
   596
	      (*WN.29.4.03 asm?: * term list??*)
neuper@37906
   597
| Substitute' of subte  (*the 'substitution': terms of type bool*) 
neuper@37906
   598
		 * term (*to be substituted in*)
neuper@37906
   599
		 * term (*resulting from the substitution*)
neuper@37906
   600
| Apply_Assumption' of term list * term
neuper@37906
   601
neuper@37906
   602
| Take' of term                         | Take_Inst' of term  
neuper@37906
   603
| Group' of (con * int list * term)
neuper@37906
   604
| Subproblem' of (spec * 
neuper@37906
   605
		  (ori list) * (*filled in assod Subproblem'*)
neuper@37906
   606
		  term *       (*-"-, headline of calc-head *)
neuper@37906
   607
		  fmz_ * 
neuper@37906
   608
		  term)        (*Subproblem(dom,pbl)*)  
neuper@37906
   609
| CAScmd' of term
neuper@37906
   610
| End_Subproblem' of term (*???*)
neuper@37906
   611
| Split_And' of term                    | Conclude_And' of term
neuper@37906
   612
| Split_Or' of term                     | Conclude_Or' of term
neuper@37906
   613
| Begin_Trans' of term                  | End_Trans' of (term * (term list))
neuper@37906
   614
| Begin_Sequ'                           | End_Sequ'(* substitute root.env*)
neuper@37906
   615
| Split_Intersect' of term              | End_Intersect' of term
neuper@37906
   616
| Check_elementwise' of (*special case:*)
neuper@37906
   617
  term *   (*(1)the current formula: [x=1,x=...]*)
neuper@37906
   618
  string * (*(2)the pred from Check_elementwise   *)
neuper@37906
   619
  (term *  (*(3)composed from (1) and (2): {x. pred}*)
neuper@37906
   620
   term list) (*20.5.03 assumptions*)
neuper@37906
   621
neuper@37906
   622
| Or_to_List' of term * term            (* (a | b, [a,b]) *)
neuper@37906
   623
| Collect_Trues' of term
neuper@37906
   624
neuper@37906
   625
| Empty_Tac_                          | Tac_ of  (*for dummies*)
neuper@37906
   626
                                            theory *
neuper@37906
   627
                                            string * (*form*)
neuper@37906
   628
					    string * (*in Tac*)
neuper@37906
   629
					    string   (*result of Tac".."*)
neuper@37906
   630
| User' (*internal for ets*)            | End_Proof'';(*End_Proof:inout*)
neuper@37906
   631
neuper@37906
   632
fun tac_2str ma = case ma of
neuper@37906
   633
    Init_Proof' (ppc, spec)  => 
neuper@37906
   634
      "Init_Proof' "^(pair2str (strs2str ppc, spec2str spec))
neuper@37906
   635
  | Model_Problem' (pblID,_,_)     => "Model_Problem' "^(strs2str pblID )
neuper@37906
   636
  | Refine_Tacitly'(p,prefin,domID,metID,itms)=> 
neuper@37906
   637
    "Refine_Tacitly' ("
neuper@37906
   638
    ^(strs2str p)^", "^(strs2str prefin)^", "
neuper@37906
   639
    ^domID^", "^(strs2str metID)^", pbl-itms)"
neuper@37906
   640
  | Refine_Problem' ms       => "Refine_Problem' ("^(*matchs2str ms*)"..."^")"
neuper@37906
   641
(*| Match_Problem' (pI, (ok, (itms, pre))) => 
neuper@37906
   642
    "Match_Problem' "^(spair2str (strs2str pI,
neuper@37906
   643
				  spair2str (bool2str ok,
neuper@37924
   644
					     spair2str ("itms2str_ itms", 
neuper@37906
   645
							"items2str pre"))))*)
neuper@37906
   646
  | Add_Given' cterm'        => "Add_Given' "(*^cterm'*)
neuper@37906
   647
  | Del_Given' cterm'        => "Del_Given' "(*^cterm'*)
neuper@37906
   648
  | Add_Find' cterm'         => "Add_Find' "(*^cterm'*)
neuper@37906
   649
  | Del_Find' cterm'         => "Del_Find' "(*^cterm'*)
neuper@37906
   650
  | Add_Relation' cterm'     => "Add_Relation' "(*^cterm'*)
neuper@37906
   651
  | Del_Relation' cterm'     => "Del_Relation' "(*^cterm'*)
neuper@37906
   652
neuper@37906
   653
  | Specify_Theory' domID    => "Specify_Theory' "^(quote domID    )
neuper@37906
   654
  | Specify_Problem' (pI, (ok, (itms, pre))) => 
neuper@37906
   655
    "Specify_Problem' "^(spair2str (strs2str pI,
neuper@37906
   656
				  spair2str (bool2str ok,
neuper@37924
   657
					     spair2str ("itms2str_ itms", 
neuper@37906
   658
							"items2str pre"))))
neuper@37906
   659
  | Specify_Method' (pI,oris,itms) => 
neuper@37906
   660
    "Specify_Method' ("^metID2str pI^", "^oris2str oris^", )"
neuper@37906
   661
neuper@37906
   662
  | Apply_Method' (metID,_,_)      => "Apply_Method' "^(strs2str metID)
neuper@37906
   663
  | Check_Postcond' (pblID,(scval,asm)) => 
neuper@37906
   664
      "Check_Postcond' "^(spair2str(strs2str pblID, 
neuper@37906
   665
				    spair2str (term2str scval, strs2str asm)))
neuper@37906
   666
neuper@37906
   667
  | Free_Solve'              => "Free_Solve'"
neuper@37906
   668
neuper@37906
   669
  | Rewrite_Inst' (*subs,thm'*) _ => 
neuper@37906
   670
      "Rewrite_Inst' "(*^(pair2str (subs2str subs, spair2str thm'))*)
neuper@37906
   671
  | Rewrite' thm'            => "Rewrite' "(*^(spair2str thm')*)
neuper@37906
   672
  | Rewrite_Asm' thm'        => "Rewrite_Asm' "(*^(spair2str thm')*)
neuper@37906
   673
  | Rewrite_Set_Inst' (*subs,thm'*) _ => 
neuper@37906
   674
      "Rewrite_Set_Inst' "(*^(pair2str (subs2str subs, quote rls))*)
neuper@37906
   675
  | Rewrite_Set'(thy',pasm,rls',f,(f',asm))          
neuper@37906
   676
    => "Rewrite_Set' ("^thy'^","^(bool2str pasm)^","^(id_rls rls')^","
neuper@37929
   677
    ^(Syntax.string_of_term (thy2ctxt' "Isac") f)^",("^(Syntax.string_of_term (thy2ctxt' "Isac") f')
neuper@37929
   678
    ^","^((strs2str o (map (Syntax.string_of_term (thy2ctxt' "Isac")))) asm)^"))"
neuper@37906
   679
neuper@37906
   680
  | End_Detail' _             => "End_Detail' xxx"
neuper@37906
   681
  | Detail_Set' _             => "Detail_Set' xxx"
neuper@37906
   682
  | Detail_Set_Inst' _        => "Detail_Set_Inst' xxx"
neuper@37906
   683
neuper@37906
   684
  | Derive' rls              => "Derive' "^id_rls rls
neuper@37906
   685
  | Calculate'  _            => "Calculate' "
neuper@37906
   686
  | Substitute' subs         => "Substitute' "(*^(subs2str subs)*)    
neuper@37906
   687
  | Apply_Assumption' ct's   => "Apply_Assumption' "(*^(strs2str ct's)*)
neuper@37906
   688
neuper@37906
   689
  | Take' cterm'             => "Take' "(*^(quote cterm'	)*)
neuper@37906
   690
  | Take_Inst' cterm'        => "Take_Inst' "(*^(quote cterm' )*)
neuper@37906
   691
  | Group' (con, ints, _)     => 
neuper@37906
   692
      "Group' "^(pair2str (con2str con, ints2str ints))
neuper@37906
   693
  | Subproblem' (spec, oris, _,_,pbl_form) => 
neuper@37906
   694
      "Subproblem' "(*^(pair2str (domID, strs2str ,...))*)
neuper@37906
   695
  | End_Subproblem'  _       => "End_Subproblem'"
neuper@37906
   696
  | CAScmd' cterm'           => "CAScmd' "(*^(quote cterm')*)
neuper@37906
   697
neuper@37906
   698
  | Empty_Tac_             => "Empty_Tac_"
neuper@37906
   699
  | User'                    => "User'"
neuper@37906
   700
  | Tac_ (_,form,id,result) => "Tac_ (thy,"^form^","^id^","^result^")"
neuper@37906
   701
  | _                       => "tac_2str not impl. for arg";
neuper@37906
   702
neuper@37906
   703
(*'executed tactics' (tac_s) with local environment etc.;
neuper@37906
   704
  used for continuing eval script + for generate*)
neuper@37906
   705
type ets =
neuper@37906
   706
    (loc_ *      (* of tactic in scr, tactic (weakly) associated with tac_*)
neuper@37906
   707
     (tac_ * 	 (* (for generate)  *)
neuper@37906
   708
      env *      (* with 'tactic=result' as a rule, tactic ev. _not_ ready:
neuper@37906
   709
		  for handling 'parallel let'*)
neuper@37906
   710
      env *      (* with results of (ready) tacs        *)
neuper@37906
   711
      term *     (* itr_arg of tactic, for upd. env at Repeat, Try*)
neuper@37906
   712
      term * 	 (* result value of the tac         *)
neuper@37906
   713
      safe))
neuper@37906
   714
    list;
neuper@37906
   715
val Ets = []:ets;
neuper@37906
   716
neuper@37906
   717
neuper@37906
   718
fun ets2s (l,(m,eno,env,iar,res,s)) = 
neuper@37906
   719
  "\n("^(loc_2str l)^",("^(tac_2str m)^
neuper@37906
   720
  ",\n  ens= "^(subst2str eno)^
neuper@37906
   721
  ",\n  env= "^(subst2str env)^
neuper@37929
   722
  ",\n  iar= "^(Syntax.string_of_term (thy2ctxt' "Isac") iar)^
neuper@37929
   723
  ",\n  res= "^(Syntax.string_of_term (thy2ctxt' "Isac") res)^
neuper@37906
   724
  ",\n  "^(safe2str s)^"))";
neuper@37906
   725
fun ets2str (ets:ets) = (strs2str o (map ets2s)) ets;
neuper@37906
   726
neuper@37906
   727
neuper@37906
   728
type envp =(*9.5.03: unused, delete with field in ptree.PblObj FIXXXME*)
neuper@37906
   729
   (int * term list) list * (*assoc-list: args of met*)
neuper@37906
   730
   (int * rls) list *       (*assoc-list: tacs already done ///15.9.00*)
neuper@37906
   731
   (int * ets) list *       (*assoc-list: tacs etc. already done*)
neuper@37906
   732
   (string * pos) list;     (*asms * from where*)
neuper@37906
   733
val empty_envp = ([],[],[],[]):envp; 
neuper@37906
   734
neuper@37906
   735
datatype ppobj = 
neuper@37906
   736
    PrfObj of {cell  : lrd option, (*where in form tac has been applied*)
neuper@37906
   737
	       (*^^^FIXME.WN0607 rename this field*)
neuper@37906
   738
	       form  : term,    
neuper@37906
   739
	       tac   : tac,         (* also in istate*)
neuper@37906
   740
	       loc   : istate option * istate option, (*for form, result 
neuper@37925
   741
13.8.02: (NONE,NONE) <==> e_istate ! see update_loc, get_loc*)
neuper@37906
   742
	       branch: branch,
neuper@37906
   743
	       result: term * term list,    
neuper@37906
   744
	       ostate: ostate}    (*Complete <=> result is OK*)
neuper@37906
   745
  | PblObj of {cell  : lrd option,(*unused: meaningful only for some _Prf_Obj*)
neuper@37906
   746
	       fmz   : fmz,       (*from init:FIXME never use this spec;-drop*)
neuper@37906
   747
	       origin: (ori list) * (*representation from fmz+pbt
neuper@37906
   748
                                  for efficiently adding items in probl, meth*)
neuper@37906
   749
		       spec *     (*updated by Refine_Tacitly*)
neuper@37906
   750
		       term,      (*headline of calc-head, as calculated 
neuper@37906
   751
							      initially(!)*)
neuper@37906
   752
		       (*# the origin of a root-pbl is created from fmz
neuper@37906
   753
                           (thus providing help for input to the user),
neuper@37906
   754
			 # the origin of a sub-pbl is created from the argument
neuper@37906
   755
			   -list of a script-tac 'SubProblem (spec) [arg-list]'
neuper@37906
   756
			   by 'match_ags'*)
neuper@37906
   757
	       spec  : spec,      (*explicitly input*)
neuper@37906
   758
	       probl : itm list,  (*itms explicitly input*)
neuper@37906
   759
	       meth  : itm list,  (*itms automatically added to copy of probl
neuper@37906
   760
				   TODO: input like to 'probl'*)
neuper@37906
   761
	       env   : istate option,(*for problem with initac in script*)
neuper@37906
   762
	       loc   : istate option * istate option, (*for pbl+met * result*)
neuper@37906
   763
	       branch: branch,
neuper@37906
   764
	       result: term * term list,
neuper@37906
   765
	       ostate: ostate};   (*Complete <=> result is _proven_ OK*)
neuper@37906
   766
neuper@37906
   767
(*.this tree contains isac's calculations; TODO.WN03 rename to ctree;
neuper@37906
   768
   the structure has been copied from an early version of Theorema(c);
neuper@37906
   769
   it has the disadvantage, that there is no space 
neuper@37906
   770
   for the first tactic in a script generating the first formula at (p,Frm);
neuper@37906
   771
   this trouble has been covered by 'init_form' and 'Take' so far,
neuper@37906
   772
   but it is crucial if the first tactic in a script is eg. 'Subproblem';
neuper@37906
   773
   see 'type tac ', Apply_Method.
neuper@37906
   774
.*)
neuper@37906
   775
datatype ptree = 
neuper@37906
   776
    EmptyPtree
neuper@37906
   777
  | Nd of ppobj * (ptree list);
neuper@37906
   778
val e_ptree = EmptyPtree;
neuper@37906
   779
neuper@37906
   780
fun rep_prfobj (PrfObj {cell,form,tac,loc,branch,result,ostate}) =
neuper@37906
   781
  {cell=cell,form=form,tac=tac,loc=loc,branch=branch,result=result,ostate=ostate};
neuper@37906
   782
fun rep_pblobj (PblObj {cell,origin,fmz,spec,probl,meth,env,
neuper@37906
   783
			loc,branch,result,ostate}) =
neuper@37906
   784
  {cell=cell,origin=origin,fmz=fmz,spec=spec,probl=probl,meth=meth,
neuper@37906
   785
   env=env,loc=loc,branch=branch,result=result,ostate=ostate};
neuper@37906
   786
fun is_prfobj (PrfObj _) = true
neuper@37906
   787
  | is_prfobj _ =false;
neuper@37906
   788
(*val is_prfobj' = get_obj is_prfobj; *)
neuper@37906
   789
fun is_pblobj (PblObj _) = true
neuper@37906
   790
  | is_pblobj _ = false;
neuper@37906
   791
(*val is_pblobj' = get_obj is_pblobj; 'Error: unbound constructor get_obj'*)
neuper@37906
   792
neuper@37906
   793
neuper@37906
   794
exception PTREE of string;
neuper@37906
   795
fun nth _ []      = raise PTREE "nth _ []"
neuper@37906
   796
  | nth 1 (x::xs) = x
neuper@37906
   797
  | nth n (x::xs) = nth (n-1) xs;
neuper@37906
   798
(*> nth 2 [11,22,33]; -->> val it = 22 : int*)
neuper@37906
   799
neuper@37906
   800
fun lev_up ([]:pos) = raise PTREE "lev_up []"
neuper@37906
   801
  | lev_up p = (drop_last p):pos;
neuper@37906
   802
fun lev_on ([]:pos) = raise PTREE "lev_on []"
neuper@37906
   803
  | lev_on pos = 
neuper@37906
   804
    let val len = length pos
neuper@37906
   805
    in (drop_last pos) @ [(nth len pos)+1] end;
neuper@37906
   806
fun lev_onFrm ((p,_):pos') = (lev_on p,Frm):pos'
neuper@37906
   807
  | lev_onFrm p = raise PTREE ("*** lev_onFrm: pos'="^(pos'2str p));
neuper@37906
   808
(*040216: for inform --> embed_deriv: remains on same level*)
neuper@37906
   809
fun lev_back (([],_):pos') = raise PTREE "lev_on_back: called by ([],_)"
neuper@37906
   810
  | lev_back (p,_) =
neuper@37906
   811
    if last_elem p <= 1 then (p, Frm):pos' 
neuper@37906
   812
    else ((drop_last p) @ [(nth (length p) p) - 1], Res);
neuper@37906
   813
(*.increase pos by n within a level.*)
neuper@37906
   814
fun pos_plus 0 pos = pos
neuper@37906
   815
  | pos_plus n ((p,Frm):pos') = pos_plus (n-1) (p, Res)
neuper@37906
   816
  | pos_plus n ((p,  _):pos') = pos_plus (n-1) (lev_on p, Res);
neuper@37906
   817
neuper@37906
   818
neuper@37906
   819
neuper@37906
   820
fun lev_pred ([]:pos) = raise PTREE "lev_pred []"
neuper@37906
   821
  | lev_pred (pos:pos) = 
neuper@37906
   822
    let val len = length pos
neuper@37906
   823
    in ((drop_last pos) @ [(nth len pos)-1]):pos end;
neuper@37906
   824
(*lev_pred [1,2,3];
neuper@37906
   825
val it = [1,2,2] : pos
neuper@37906
   826
> lev_pred [1];
neuper@37906
   827
val it = [0] : pos          *)
neuper@37906
   828
neuper@37906
   829
fun lev_dn p = p @ [0];
neuper@37906
   830
(*> (lev_dn o lev_on) [1,2,3];
neuper@37906
   831
val it = [1,2,4,0] : pos    *)
neuper@37906
   832
(*fun lev_dn' ((p,p_):pos') = (lev_dn p, Frm):pos'; WN.3.12.03: never used*)
neuper@37906
   833
fun lev_dnRes ((p,_):pos') = (lev_dn p, Res):pos';
neuper@37906
   834
neuper@37906
   835
(*4.4.00*)
neuper@37906
   836
fun lev_up_ ((p,Res):pos') = (lev_up p,Res):pos'
neuper@38031
   837
  | lev_up_ p' = error ("lev_up_: called for "^(pos'2str p'));
neuper@37906
   838
fun lev_dn_ ((p,_):pos') = (lev_dn p,Res):pos'
neuper@37906
   839
fun ind ((p,_):pos') = length p; (*WN050108 deprecated in favour of lev_of*)
neuper@37906
   840
fun lev_of ((p,_):pos') = length p;
neuper@37906
   841
neuper@37906
   842
neuper@37906
   843
(** convert ptree to a string **)
neuper@37906
   844
neuper@37906
   845
(* convert a pos from list to string *)
neuper@37906
   846
fun pr_pos ps = (space_implode "." (map string_of_int ps))^".   ";
neuper@37906
   847
(* show hd origin or form only *)
neuper@37906
   848
fun pr_short (p:pos) (PblObj {origin = (ori,_,_),...}) = 
neuper@37906
   849
  ((pr_pos p) ^ " ----- pblobj -----\n")
neuper@37929
   850
(*   ((((Syntax.string_of_term (thy2ctxt' "Isac")) o #4 o hd) ori)^" "^
neuper@37929
   851
    (((Syntax.string_of_term (thy2ctxt' "Isac")) o hd(*!?!*) o #5 o hd) ori))^
neuper@37906
   852
   "\n") *)
neuper@37906
   853
  | pr_short p (PrfObj {form = form,...}) =
neuper@37906
   854
  ((pr_pos p) ^ (term2str form) ^ "\n");
neuper@37906
   855
(*
neuper@37906
   856
fun pr_cell (p:pos) (PblObj {cell = c, origin = (ori,_,_),...}) = 
neuper@37906
   857
  ((ints2str c) ^"   "^ 
neuper@37929
   858
   ((((Syntax.string_of_term (thy2ctxt' "Isac")) o #4 o hd) ori)^" "^
neuper@37929
   859
    (((Syntax.string_of_term (thy2ctxt' "Isac")) o hd(*!?!*) o #5 o hd) ori))^
neuper@37906
   860
   "\n")
neuper@37906
   861
  | pr_cell p (PrfObj {cell = c, form = form,...}) =
neuper@37906
   862
  ((ints2str c) ^"   "^ (term2str form) ^ "\n");
neuper@37906
   863
*)
neuper@37906
   864
neuper@37906
   865
(* convert ptree *)
neuper@37906
   866
fun pr_ptree f pt =
neuper@37906
   867
  let
neuper@37906
   868
    fun pr_pt pfn _  EmptyPtree = ""
neuper@37906
   869
      | pr_pt pfn ps (Nd (b, [])) = pfn ps b
neuper@37906
   870
      | pr_pt pfn ps (Nd (b, ts)) = (pfn ps b)^
neuper@37906
   871
      (prts pfn (ps:pos) 1 ts)
neuper@37906
   872
    and prts pfn ps p [] = ""
neuper@37906
   873
      | prts pfn ps p (t::ts) = (pr_pt pfn (ps @ [p]) t)^
neuper@37906
   874
      (prts pfn ps (p+1) ts)
neuper@37906
   875
  in pr_pt f [] pt end;
neuper@37906
   876
(*
neuper@37906
   877
> fun prfn ps b = (pr_pos ps)^"   "^b(*TODO*)^"\n";
neuper@38006
   878
> val pt = Unsynchronized.ref EmptyPtree;
neuper@37986
   879
> pt:=Nd("root'",
neuper@37906
   880
       [Nd("xx1",[]),
neuper@37906
   881
	Nd("xx2",
neuper@37906
   882
	   [Nd("xx2.1.",[]),
neuper@37906
   883
	    Nd("xx2.2.",[])]),
neuper@37906
   884
	Nd("xx3",[])]);
neuper@38015
   885
> tracing (pr_ptree prfn (!pt));
neuper@37906
   886
*)
neuper@37906
   887
neuper@37906
   888
neuper@37906
   889
(** access the branches of ptree **)
neuper@37906
   890
neuper@37906
   891
fun ins_nth 1 e l  = e::l
neuper@37906
   892
  | ins_nth n e [] = raise PTREE "ins_nth n e []"
neuper@37906
   893
  | ins_nth n e (l::ls) = l::(ins_nth (n-1) e ls);
neuper@37906
   894
fun repl []      _ _ = raise PTREE "repl [] _ _"
neuper@37906
   895
  | repl (l::ls) 1 e = e::ls
neuper@37906
   896
  | repl (l::ls) n e = l::(repl ls (n-1) e);
neuper@37906
   897
fun repl_app ls n e = 
neuper@37906
   898
    let val lim = 1 + length ls
neuper@37906
   899
    in if n > lim then raise PTREE "repl_app: n > lim"
neuper@37906
   900
       else if n = lim then ls @ [e]
neuper@37906
   901
	    else repl ls n e end;
neuper@37906
   902
(*  
neuper@37906
   903
> repl [1,2,3] 2 22222;
neuper@37906
   904
val it = [1,22222,3] : int list
neuper@37906
   905
> repl_app [1,2,3,4] 5 5555;
neuper@37906
   906
val it = [1,2,3,4,5555] : int list
neuper@37906
   907
> repl_app [1,2,3] 2 22222;
neuper@37906
   908
val it = [1,22222,3] : int list
neuper@37906
   909
> repl_app [1] 2 22222 ;
neuper@37906
   910
val it = [1,22222] : int list
neuper@37906
   911
*)
neuper@37906
   912
neuper@37906
   913
neuper@37906
   914
(*.get from obj at pos by f : ppobj -> 'a.*)
neuper@37906
   915
fun get_obj f EmptyPtree  (_:pos)  = raise PTREE "get_obj f EmptyPtree"
neuper@37906
   916
  | get_obj f (Nd (b,  _)) []      = f b
neuper@37906
   917
  | get_obj f (Nd (b, bs)) (p::ps) =
neuper@37906
   918
(* val (f, Nd (b, bs), (p::ps)) = (I, pt, p);
neuper@37906
   919
   *)
neuper@37906
   920
  let val _ = (nth p bs) handle _ => raise PTREE ("get_obj: pos = "^
neuper@37906
   921
			   (ints2str' (p::ps))^" does not exist");
neuper@37906
   922
  in (get_obj f (nth p bs) (ps:pos)) 
neuper@37906
   923
      (*before WN050419: 'wrong type..' raised also if pos doesn't exist*)
neuper@37906
   924
    handle _ => raise PTREE (*"get_obj: at pos = "^
neuper@37906
   925
			     (ints2str' (p::ps))^" wrong type of ppobj"*)
neuper@37906
   926
			  ("get_obj: pos = "^
neuper@37906
   927
			   (ints2str' (p::ps))^" does not exist")
neuper@37906
   928
  end;
neuper@37906
   929
fun get_nd EmptyPtree _ = raise PTREE "get_nd EmptyPtree"
neuper@37906
   930
  | get_nd n [] = n
neuper@37906
   931
  | get_nd (Nd (_,nds)) (pos as p::(ps:pos)) = (get_nd (nth p nds) ps)
neuper@37906
   932
    handle _ => raise PTREE ("get_nd: not existent pos = "^(ints2str' pos));
neuper@37906
   933
neuper@37906
   934
neuper@37906
   935
(* for use by get_obj *)
neuper@37925
   936
fun g_cell   (PblObj {cell = c,...}) = NONE
neuper@37906
   937
  | g_cell   (PrfObj {cell = c,...}) = c;(*WN0607 hack for quick introduction of lrd + rewrite-at (thms, calcs)*)
neuper@37906
   938
fun g_form   (PrfObj {form = f,...}) = f
neuper@37906
   939
  | g_form   (PblObj {origin=(_,_,f),...}) = f;
neuper@37906
   940
fun g_form' (Nd (PrfObj {form = f,...}, _)) = f
neuper@37906
   941
  | g_form' (Nd (PblObj {origin=(_,_,f),...}, _)) = f;
neuper@37906
   942
(*  | g_form   _ = raise PTREE "g_form not for PblObj";*)
neuper@37906
   943
fun g_origin (PblObj {origin = ori,...}) = ori
neuper@37906
   944
  | g_origin _ = raise PTREE "g_origin not for PrfObj";
neuper@37906
   945
fun g_fmz (PblObj {fmz = f,...}) = f
neuper@37906
   946
  | g_fmz _ = raise PTREE "g_fmz not for PrfObj";
neuper@37906
   947
fun g_spec   (PblObj {spec = s,...}) = s
neuper@37906
   948
  | g_spec _   = raise PTREE "g_spec not for PrfObj";
neuper@37906
   949
fun g_pbl    (PblObj {probl = p,...}) = p
neuper@37906
   950
  | g_pbl  _   = raise PTREE "g_pbl not for PrfObj";
neuper@37906
   951
fun g_met    (PblObj {meth = p,...}) = p
neuper@37906
   952
  | g_met  _   = raise PTREE "g_met not for PrfObj";
neuper@37906
   953
fun g_domID  (PblObj {spec = (d,_,_),...}) = d
neuper@37906
   954
  | g_domID  _ = raise PTREE "g_metID not for PrfObj";
neuper@37906
   955
fun g_metID  (PblObj {spec = (_,_,m),...}) = m
neuper@37906
   956
  | g_metID  _ = raise PTREE "g_metID not for PrfObj";
neuper@37906
   957
fun g_env    (PblObj {env,...}) = env
neuper@37906
   958
  | g_env    _ = raise PTREE "g_env not for PrfObj"; 
neuper@37906
   959
fun g_loc    (PblObj {loc = l,...}) = l
neuper@37906
   960
  | g_loc    (PrfObj {loc = l,...}) = l;
neuper@37906
   961
fun g_branch (PblObj {branch = b,...}) = b
neuper@37906
   962
  | g_branch (PrfObj {branch = b,...}) = b;
neuper@37906
   963
fun g_tac  (PblObj {spec = (d,p,m),...}) = Apply_Method m
neuper@37906
   964
  | g_tac  (PrfObj {tac = m,...}) = m;
neuper@37906
   965
fun g_result (PblObj {result = r,...}) = r
neuper@37906
   966
  | g_result (PrfObj {result = r,...}) = r;
neuper@37906
   967
fun g_res (PblObj {result = (r,_),...}) = r
neuper@37906
   968
  | g_res (PrfObj {result = (r,_),...}) = r;
neuper@37906
   969
fun g_res' (Nd (PblObj {result = (r,_),...}, _)) = r
neuper@37906
   970
  | g_res' (Nd (PrfObj {result = (r,_),...}, _)) = r;
neuper@37906
   971
fun g_ostate (PblObj {ostate = r,...}) = r
neuper@37906
   972
  | g_ostate (PrfObj {ostate = r,...}) = r;
neuper@37906
   973
fun g_ostate' (Nd (PblObj {ostate = r,...}, _)) = r
neuper@37906
   974
  | g_ostate' (Nd (PrfObj {ostate = r,...}, _)) = r;
neuper@37906
   975
neuper@37925
   976
fun gpt_cell (Nd (PblObj {cell = c,...},_)) = NONE
neuper@37906
   977
  | gpt_cell (Nd (PrfObj {cell = c,...},_)) = c;
neuper@37906
   978
neuper@37906
   979
(*in CalcTree/Subproblem an 'just_created_' model is created;
neuper@37906
   980
  this is filled to 'untouched' by Model/Refine_Problem*)
neuper@37906
   981
fun just_created_ (PblObj {meth, probl, spec, ...}) = 
neuper@37906
   982
    null meth andalso null probl andalso spec = e_spec;
neuper@37906
   983
val e_origin = ([],e_spec,e_term): (ori list) * spec * term;
neuper@37906
   984
neuper@37906
   985
fun just_created (pt,(p,_):pos') =
neuper@37906
   986
    let val ppobj = get_obj I pt p
neuper@37906
   987
    in is_pblobj ppobj andalso just_created_ ppobj end;
neuper@37906
   988
neuper@37906
   989
(*.does the pos in the ctree exist ?.*)
neuper@37906
   990
fun existpt pos pt = can (get_obj I pt) pos;
neuper@37906
   991
(*.does the pos' in the ctree exist, ie. extra check for result in the node.*)
neuper@37906
   992
fun existpt' ((p,p_):pos') pt = 
neuper@37906
   993
    if can (get_obj I pt) p 
neuper@37906
   994
    then case p_ of 
neuper@37906
   995
	     Res => get_obj g_ostate pt p = Complete
neuper@37906
   996
	   | _ => true
neuper@37906
   997
    else false;
neuper@37906
   998
neuper@37906
   999
(*.is this position appropriate for calculating intermediate steps?.*)
neuper@37906
  1000
fun is_interpos ((_, Res):pos') = true
neuper@37906
  1001
  | is_interpos _ = false;
neuper@37906
  1002
neuper@37906
  1003
fun last_onlev pt pos = not (existpt (lev_on pos) pt);
neuper@37906
  1004
neuper@37906
  1005
neuper@37906
  1006
(*.find the position of the next parent which is a PblObj in ptree.*)
neuper@37906
  1007
fun par_pblobj pt ([]:pos) = ([]:pos)
neuper@37906
  1008
  | par_pblobj pt p =
neuper@37906
  1009
    let fun par pt [] = []
neuper@37906
  1010
	  | par pt p = if is_pblobj (get_obj I pt p) then p
neuper@37906
  1011
		       else par pt (lev_up p)
neuper@37906
  1012
    in par pt (lev_up p) end; 
neuper@37906
  1013
(* lev_up for hard_gen operating with pos = [...,0] *)
neuper@37906
  1014
neuper@37906
  1015
(*.find the position and the children of the next parent which is a PblObj.*)
neuper@37906
  1016
fun par_children (Nd (PblObj _, children)) ([]:pos) = (children, []:pos)
neuper@37906
  1017
  | par_children (pt as Nd (PblObj _, children)) p =
neuper@37906
  1018
    let fun par [] = (children, [])
neuper@37906
  1019
	  | par p = let val Nd (obj, children) = get_nd pt p
neuper@37906
  1020
		    in if is_pblobj obj then (children, p) else par (lev_up p)
neuper@37906
  1021
		    end;
neuper@37906
  1022
    in par (lev_up p) end; 
neuper@37906
  1023
neuper@37906
  1024
(*.get the children of a node in ptree.*)
neuper@37906
  1025
fun children (Nd (PblObj _, cn)) = cn
neuper@37906
  1026
  | children (Nd (PrfObj _, cn)) = cn;
neuper@37906
  1027
neuper@37906
  1028
neuper@37906
  1029
(*.find the next parent, which is either a PblObj (return true)
neuper@37906
  1030
  or a PrfObj with tac = Detail_Set (return false).*)
neuper@37906
  1031
(*FIXME.3.4.03:re-organize par_pbl_det after rls' --> rls*)
neuper@37906
  1032
fun par_pbl_det pt ([]:pos) = (true, []:pos, Erls)
neuper@37906
  1033
  | par_pbl_det pt p =
neuper@37906
  1034
    let fun par pt [] = (true, [], Erls)
neuper@37906
  1035
	  | par pt p = if is_pblobj (get_obj I pt p) then (true, p, Erls)
neuper@37906
  1036
		       else case get_obj g_tac pt p of
neuper@37906
  1037
				(*Detail_Set rls' => (false, p, assoc_rls rls')
neuper@37906
  1038
			      (*^^^--- before 040206 after ---vvv*)
neuper@37906
  1039
			      |*)Rewrite_Set rls' => (false, p, assoc_rls rls')
neuper@37906
  1040
			      | Rewrite_Set_Inst (_, rls') => 
neuper@37906
  1041
				(false, p, assoc_rls rls')
neuper@37906
  1042
			      | _ => par pt (lev_up p)
neuper@37906
  1043
    in par pt (lev_up p) end; 
neuper@37906
  1044
neuper@37906
  1045
neuper@37906
  1046
neuper@37906
  1047
neuper@37906
  1048
(*.get from the whole ptree by f : ppobj -> 'a.*)
neuper@37906
  1049
fun get_all f EmptyPtree   = []
neuper@37906
  1050
  | get_all f (Nd (b, [])) = [f b]
neuper@37906
  1051
  | get_all f (Nd (b, bs)) = [f b] @ (get_alls f bs)
neuper@37906
  1052
and get_alls f [] = []
neuper@37906
  1053
  | get_alls f pts = flat (map (get_all f) pts);
neuper@37906
  1054
neuper@37906
  1055
neuper@37906
  1056
(*.insert obj b into ptree at pos, ev.overwriting this pos.*)
neuper@37906
  1057
fun insert b EmptyPtree   ([]:pos)  = Nd (b, [])
neuper@37906
  1058
  | insert b EmptyPtree    _        = raise PTREE "insert b Empty _"
neuper@37906
  1059
  | insert b (Nd ( _,  _)) []       = raise PTREE "insert b _ []"
neuper@37906
  1060
  | insert b (Nd (b', bs)) (p::[])  = 
neuper@37906
  1061
     Nd (b', repl_app bs p (Nd (b,[]))) 
neuper@37906
  1062
  | insert b (Nd (b', bs)) (p::ps)  =
neuper@37906
  1063
     Nd (b', repl_app bs p (insert b (nth p bs) ps));
neuper@37906
  1064
(*
neuper@37906
  1065
> type ppobj = string;
neuper@38015
  1066
> tracing (pr_ptree prfn (!pt));
neuper@38006
  1067
  val pt = Unsynchronized.ref Empty;
neuper@37986
  1068
  pt:= insert ("root'":ppobj) EmptyPtree [];
neuper@37906
  1069
  pt:= insert ("xx1":ppobj) (!pt) [1];
neuper@37906
  1070
  pt:= insert ("xx2":ppobj) (!pt) [2];
neuper@37906
  1071
  pt:= insert ("xx3":ppobj) (!pt) [3];
neuper@37906
  1072
  pt:= insert ("xx2.1":ppobj) (!pt) [2,1];
neuper@37906
  1073
  pt:= insert ("xx2.2":ppobj) (!pt) [2,2];
neuper@37906
  1074
  pt:= insert ("xx2.1.1":ppobj) (!pt) [2,1,1];
neuper@37906
  1075
  pt:= insert ("xx2.1.2":ppobj) (!pt) [2,1,2];
neuper@37906
  1076
  pt:= insert ("xx2.1.3":ppobj) (!pt) [2,1,3];
neuper@37906
  1077
*)
neuper@37906
  1078
neuper@37906
  1079
(*.insert children to a node without children.*)
neuper@37906
  1080
(*compare: fun insert*)
neuper@37906
  1081
fun ins_chn _  EmptyPtree   (_:pos) = raise PTREE "ins_chn: EmptyPtree"
neuper@37906
  1082
  | ins_chn ns (Nd _)       []      = raise PTREE "ins_chn: pos = []"
neuper@37906
  1083
  | ins_chn ns (Nd (b, bs)) (p::[]) =
neuper@37906
  1084
    if p > length bs then raise PTREE "ins_chn: pos not existent"
neuper@37906
  1085
    else let val Nd (b', bs') = nth p bs
neuper@37906
  1086
	 in if null bs' then Nd (b, repl_app bs p (Nd (b', ns)))
neuper@37906
  1087
	    else raise PTREE "ins_chn: pos mustNOT be overwritten" end
neuper@37906
  1088
  | ins_chn ns (Nd (b, bs)) (p::ps) =
neuper@37906
  1089
     Nd (b, repl_app bs p (ins_chn ns (nth p bs) ps));
neuper@37906
  1090
neuper@37906
  1091
(* print_depth 11;ins_chn;print_depth 3; ###insert#########################*);
neuper@37906
  1092
neuper@37906
  1093
neuper@37906
  1094
(** apply f to obj at pos, f: ppobj -> ppobj **)
neuper@37906
  1095
neuper@37906
  1096
fun appl_to_node f (Nd (b,bs)) = Nd (f b, bs);
neuper@37906
  1097
fun appl_obj f EmptyPtree    []      = EmptyPtree
neuper@37906
  1098
  | appl_obj f EmptyPtree    _       = raise PTREE "appl_obj f Empty _"
neuper@37906
  1099
  | appl_obj f (Nd (b, bs)) []       = Nd (f b, bs)
neuper@37906
  1100
  | appl_obj f (Nd (b, bs)) (p::[])  = 
neuper@37906
  1101
     Nd (b, repl_app bs p (((appl_to_node f) o (nth p)) bs))
neuper@37906
  1102
  | appl_obj f (Nd (b, bs)) (p::ps)  =
neuper@37906
  1103
     Nd (b, repl_app bs p (appl_obj f (nth p bs) (ps:pos)));
neuper@37906
  1104
 
neuper@37906
  1105
(* for use by appl_obj *) 
neuper@37906
  1106
fun repl_form f (PrfObj {cell=c,form= _,tac=tac,loc=loc,
neuper@37906
  1107
			 branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1108
    PrfObj {cell=c,form= f,tac=tac,loc=loc,
neuper@37906
  1109
	    branch=branch,result=result,ostate=ostate}
neuper@37906
  1110
  | repl_form _ _ = raise PTREE "repl_form takes no PblObj";
neuper@37906
  1111
fun repl_pbl x    (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1112
			   spec=spec,probl=_,meth=meth,env=env,loc=loc,
neuper@37906
  1113
			   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1114
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=spec,probl= x,
neuper@37906
  1115
	  meth=meth,env=env,loc=loc,branch=branch,result=result,ostate=ostate}
neuper@37906
  1116
  | repl_pbl _ _ = raise PTREE "repl_pbl takes no PrfObj";
neuper@37906
  1117
fun repl_met x    (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1118
			   spec=spec,probl=probl,meth=_,env=env,loc=loc,
neuper@37906
  1119
			   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1120
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1121
	  meth= x,env=env,loc=loc,branch=branch,result=result,ostate=ostate}
neuper@37906
  1122
  | repl_met _ _ = raise PTREE "repl_pbl takes no PrfObj";
neuper@37906
  1123
neuper@37906
  1124
fun repl_spec  x    (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1125
			   spec= _,probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1126
			   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1127
  PblObj {cell=cell,origin=origin,fmz=fmz,spec= x,probl=probl,
neuper@37906
  1128
	  meth=meth,env=env,loc=loc,branch=branch,result=result,ostate=ostate}
neuper@37906
  1129
  | repl_spec  _ _ = raise PTREE "repl_domID takes no PrfObj";
neuper@37906
  1130
fun repl_domID x    (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1131
			   spec=(_,p,m),probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1132
			   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1133
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=(x,p,m),probl=probl,
neuper@37906
  1134
	  meth=meth,env=env,loc=loc,branch=branch,result=result,ostate=ostate}
neuper@37906
  1135
  | repl_domID _ _ = raise PTREE "repl_domID takes no PrfObj";
neuper@37906
  1136
fun repl_pblID x    (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1137
			   spec=(d,_,m),probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1138
			   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1139
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=(d,x,m),probl=probl,
neuper@37906
  1140
	  meth=meth,env=env,loc=loc,branch=branch,result=result,ostate=ostate}
neuper@37906
  1141
  | repl_pblID _ _ = raise PTREE "repl_pblID takes no PrfObj";
neuper@37906
  1142
fun repl_metID x (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1143
			   spec=(d,p,_),probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1144
			   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1145
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=(d,p,x),probl=probl,
neuper@37906
  1146
	  meth=meth,env=env,loc=loc,branch=branch,result=result,ostate=ostate}
neuper@37906
  1147
  | repl_metID _ _ = raise PTREE "repl_metID takes no PrfObj";
neuper@37906
  1148
neuper@37906
  1149
fun repl_result l f' s (PrfObj {cell=cell,form=form,tac=tac,loc=_,
neuper@37906
  1150
			     branch=branch,result = _ ,ostate = _}) =
neuper@37906
  1151
    PrfObj {cell=cell,form=form,tac=tac,loc= l,
neuper@37906
  1152
	    branch=branch,result = f',ostate = s}
neuper@37906
  1153
  | repl_result l f' s (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1154
			     spec=spec,probl=probl,meth=meth,env=env,loc=_,
neuper@37906
  1155
			     branch=branch,result= _ ,ostate= _}) =
neuper@37906
  1156
    PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1157
	    spec=spec,probl=probl,meth=meth,env=env,loc= l,
neuper@37906
  1158
	    branch=branch,result= f',ostate= s};
neuper@37906
  1159
neuper@37906
  1160
fun repl_tac x (PrfObj {cell=cell,form=form,tac= _,loc=loc,
neuper@37906
  1161
			  branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1162
    PrfObj {cell=cell,form=form,tac= x,loc=loc,
neuper@37906
  1163
	    branch=branch,result=result,ostate=ostate}
neuper@37906
  1164
  | repl_tac _ _ = raise PTREE "repl_tac takes no PblObj";
neuper@37906
  1165
neuper@37906
  1166
fun repl_branch b (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1167
			   spec=spec,probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1168
			   branch= _,result=result,ostate=ostate}) =
neuper@37906
  1169
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1170
	  meth=meth,env=env,loc=loc,branch= b,result=result,ostate=ostate}
neuper@37906
  1171
  | repl_branch b (PrfObj {cell=cell,form=form,tac=tac,loc=loc,
neuper@37906
  1172
			  branch= _,result=result,ostate=ostate}) =
neuper@37906
  1173
    PrfObj {cell=cell,form=form,tac=tac,loc=loc,
neuper@37906
  1174
	    branch= b,result=result,ostate=ostate};
neuper@37906
  1175
neuper@37906
  1176
fun repl_env e
neuper@37906
  1177
  (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1178
	   spec=spec,probl=probl,meth=meth,env=_,loc=loc,
neuper@37906
  1179
	   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1180
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1181
	  meth=meth,env=e,loc=loc,branch=branch,
neuper@37906
  1182
	  result=result,ostate=ostate}
neuper@37906
  1183
  | repl_env _ _ = raise PTREE "repl_ets takes no PrfObj";
neuper@37906
  1184
neuper@37906
  1185
fun repl_oris oris
neuper@37906
  1186
  (PblObj {cell=cell,origin=(_,spe,hdf),fmz=fmz,
neuper@37906
  1187
	   spec=spec,probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1188
	   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1189
  PblObj{cell=cell,origin=(oris,spe,hdf),fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1190
	  meth=meth,env=env,loc=loc,branch=branch,
neuper@37906
  1191
	  result=result,ostate=ostate}
neuper@37906
  1192
  | repl_oris _ _ = raise PTREE "repl_oris takes no PrfObj";
neuper@37906
  1193
fun repl_orispec spe
neuper@37906
  1194
  (PblObj {cell=cell,origin=(oris,_,hdf),fmz=fmz,
neuper@37906
  1195
	   spec=spec,probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1196
	   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1197
  PblObj{cell=cell,origin=(oris,spe,hdf),fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1198
	  meth=meth,env=env,loc=loc,branch=branch,
neuper@37906
  1199
	  result=result,ostate=ostate}
neuper@37906
  1200
  | repl_orispec _ _ = raise PTREE "repl_orispec takes no PrfObj";
neuper@37906
  1201
neuper@37906
  1202
fun repl_loc l (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1203
			spec=spec,probl=probl,meth=meth,env=env,loc=_,
neuper@37906
  1204
			branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1205
  PblObj {cell=cell,origin=origin,fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1206
	  meth=meth,env=env,loc=l,branch=branch,result=result,ostate=ostate}
neuper@37906
  1207
  | repl_loc l (PrfObj {cell=cell,form=form,tac=tac,loc=_,
neuper@37906
  1208
			branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1209
  PrfObj {cell=cell,form=form,tac=tac,loc= l,
neuper@37906
  1210
	  branch=branch,result=result,ostate=ostate};
neuper@37906
  1211
(*
neuper@37906
  1212
fun uni__cid cell' 
neuper@37906
  1213
  (PblObj {cell=cell,origin=origin,fmz=fmz,
neuper@37906
  1214
	   spec=spec,probl=probl,meth=meth,env=env,loc=loc,
neuper@37906
  1215
	   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1216
  PblObj {cell=cell union cell',origin=origin,fmz=fmz,spec=spec,probl=probl,
neuper@37906
  1217
	  meth=meth,env=env,loc=loc,branch=branch,
neuper@37906
  1218
	  result=result,ostate=ostate}
neuper@37906
  1219
  | uni__cid cell'
neuper@37906
  1220
  (PrfObj {cell=cell,form=form,tac=tac,loc=loc,
neuper@37906
  1221
	   branch=branch,result=result,ostate=ostate}) =
neuper@37906
  1222
  PrfObj {cell=cell union cell',form=form,tac=tac,loc=loc,
neuper@37906
  1223
	  branch=branch,result=result,ostate=ostate};
neuper@37906
  1224
*)
neuper@37906
  1225
neuper@37906
  1226
(*WN050219 put here for interpreting code for cut_tree below...*)
neuper@37906
  1227
type ocalhd =
neuper@37906
  1228
     bool *                (*ALL itms+preconds true*)
neuper@37906
  1229
     pos_ *                (*model belongs to Problem | Method*)
neuper@37906
  1230
     term *                (*header: Problem... or Cas
neuper@37906
  1231
				FIXXXME.12.03: item! for marking syntaxerrors*)
neuper@37906
  1232
     itm list *            (*model: given, find, relate*)
neuper@37906
  1233
     ((bool * term) list) *(*model: preconds*)
neuper@37906
  1234
     spec;                 (*specification*)
neuper@37906
  1235
val e_ocalhd = (false, Und, e_term, [e_itm], [(false, e_term)], e_spec);
neuper@37906
  1236
neuper@37906
  1237
datatype ptform =
neuper@37906
  1238
	 Form of term
neuper@37906
  1239
       | ModSpec of ocalhd;
neuper@37906
  1240
val e_ptform = Form e_term;
neuper@37906
  1241
val e_ptform' = ModSpec e_ocalhd;
neuper@37906
  1242
neuper@37906
  1243
neuper@37906
  1244
neuper@37906
  1245
(*.applies (snd f) to the branches at a pos if ((fst f) b),
neuper@37906
  1246
   f : (ppobj -> bool) * (int -> ptree list -> ptree list).*)
neuper@37906
  1247
neuper@37906
  1248
fun appl_branch f EmptyPtree [] = (EmptyPtree, false)
neuper@37906
  1249
  | appl_branch f EmptyPtree _  = raise PTREE "appl_branch f Empty _"
neuper@37906
  1250
  | appl_branch f (Nd ( _, _)) [] = raise PTREE "appl_branch f _ []"
neuper@37906
  1251
  | appl_branch f (Nd (b, bs)) (p::[]) = 
neuper@37906
  1252
    if (fst f) b then (Nd (b, (snd f) (p:posel) bs), true)
neuper@37906
  1253
    else (Nd (b, bs), false)
neuper@37906
  1254
  | appl_branch f (Nd (b, bs)) (p::ps) =
neuper@37906
  1255
	let val (b',bool) = appl_branch f (nth p bs) ps
neuper@37906
  1256
	in (Nd (b, repl_app bs p b'), bool) end;
neuper@37906
  1257
neuper@37906
  1258
(* for cut_level;  appl_branch(deprecated) *)
neuper@37906
  1259
fun test_trans (PrfObj{branch = Transitive,...}) = true
neuper@37906
  1260
  | test_trans (PblObj{branch = Transitive,...}) = true
neuper@37906
  1261
  | test_trans _ = false;
neuper@37906
  1262
neuper@37906
  1263
fun is_pblobj' pt (p:pos) =
neuper@37906
  1264
    let val ppobj = get_obj I pt p
neuper@37906
  1265
    in is_pblobj ppobj end;
neuper@37906
  1266
neuper@37906
  1267
neuper@37906
  1268
fun delete_result pt (p:pos) =
neuper@37925
  1269
    (appl_obj (repl_result (fst (get_obj g_loc pt p), NONE) 
neuper@37906
  1270
			   (e_term,[]) Incomplete) pt p);
neuper@37906
  1271
neuper@37906
  1272
fun del_res (PblObj {cell, fmz, origin, spec, probl, meth, 
neuper@37906
  1273
		     env, loc=(l1,_), branch, result, ostate}) =
neuper@37906
  1274
    PblObj {cell=cell,fmz=fmz,origin=origin,spec=spec,probl=probl,meth=meth,
neuper@37925
  1275
	    env=env, loc=(l1,NONE), branch=branch, result=(e_term,[]), 
neuper@37906
  1276
	    ostate=Incomplete}
neuper@37906
  1277
neuper@37906
  1278
  | del_res (PrfObj {cell, form, tac, loc=(l1,_), branch, result, ostate}) =
neuper@37925
  1279
    PrfObj {cell=cell,form=form,tac=tac, loc=(l1,NONE), branch=branch, 
neuper@37906
  1280
	    result=(e_term,[]), ostate=Incomplete};
neuper@37906
  1281
neuper@37906
  1282
neuper@37906
  1283
(*
neuper@37906
  1284
fun update_fmz  pt pos x = appl_obj (repl_fmz  x) pt pos;
neuper@37906
  1285
                                       1.00 not used anymore*)
neuper@37906
  1286
neuper@37906
  1287
(*FIXME.WN.12.03: update_X X pos pt -> pt could be chained by o (efficiency?)*)
neuper@37906
  1288
fun update_env    pt pos x = appl_obj (repl_env    x) pt pos;
neuper@37906
  1289
fun update_domID  pt pos x = appl_obj (repl_domID  x) pt pos;
neuper@37906
  1290
fun update_pblID  pt pos x = appl_obj (repl_pblID  x) pt pos;
neuper@37906
  1291
fun update_metID  pt pos x = appl_obj (repl_metID  x) pt pos;
neuper@37906
  1292
fun update_spec   pt pos x = appl_obj (repl_spec   x) pt pos;
neuper@37906
  1293
neuper@37906
  1294
fun update_pbl    pt pos x = appl_obj (repl_pbl    x) pt pos;
neuper@37906
  1295
fun update_pblppc pt pos x = appl_obj (repl_pbl    x) pt pos;
neuper@37906
  1296
neuper@37906
  1297
fun update_met    pt pos x = appl_obj (repl_met    x) pt pos;
neuper@37906
  1298
(*1.09.01 ----
neuper@37906
  1299
fun update_metppc pt pos x = 
neuper@37906
  1300
  let val {rew_ord'=od,rls'=rs,asm_thm=at,asm_rls=ar,...} =
neuper@37906
  1301
    get_obj g_met pt pos
neuper@37906
  1302
  in appl_obj (repl_met 
neuper@37906
  1303
     {rew_ord'=od,rls'=rs,asm_thm=at,asm_rls=ar,ppc=x}) 
neuper@37906
  1304
    pt pos end;*)
neuper@37906
  1305
fun update_metppc pt pos x = appl_obj (repl_met    x) pt pos;
neuper@37906
  1306
			 			   
neuper@37906
  1307
(*fun union_cid     pt pos x = appl_obj (uni__cid    x) pt pos;*)
neuper@37906
  1308
neuper@37906
  1309
fun update_branch pt pos x = appl_obj (repl_branch x) pt pos;
neuper@37906
  1310
fun update_tac  pt pos x = appl_obj (repl_tac  x) pt pos;
neuper@37906
  1311
neuper@37906
  1312
fun update_oris   pt pos x = appl_obj (repl_oris   x) pt pos;
neuper@37906
  1313
fun update_orispec   pt pos x = appl_obj (repl_orispec   x) pt pos;
neuper@37906
  1314
neuper@37906
  1315
 (*done by append_* !! 3.5.02;  ununsed WN050305 thus outcommented
neuper@37925
  1316
fun update_loc pt (p,_) (ScrState ([],[],NONE,
neuper@37906
  1317
				   Const ("empty",_),Sundef,false)) = 
neuper@37925
  1318
    appl_obj (repl_loc (NONE,NONE)) pt p
neuper@37906
  1319
  | update_loc pt (p,Res) x =  
neuper@37906
  1320
    let val (lform,_) = get_obj g_loc pt p
neuper@37925
  1321
    in appl_obj (repl_loc (lform,SOME x)) pt p end
neuper@37906
  1322
neuper@37906
  1323
  | update_loc pt (p,_) x = 
neuper@37906
  1324
    let val (_,lres) = get_obj g_loc pt p
neuper@37925
  1325
    in appl_obj (repl_loc (SOME x,lres)) pt p end;-------------*)
neuper@37906
  1326
neuper@37906
  1327
(*WN050305 for handling cut_tree in cappend_atomic -- TODO redesign !*)
neuper@37906
  1328
fun update_loc' pt p iss = appl_obj (repl_loc iss) pt p;
neuper@37906
  1329
neuper@37906
  1330
(*13.8.02---------------------------
neuper@37925
  1331
fun get_loc EmptyPtree _ = NONE
neuper@37906
  1332
  | get_loc pt (p,Res) =
neuper@37906
  1333
  let val (lfrm,lres) = get_obj g_loc pt p
neuper@37906
  1334
  in if lres = e_istate then lfrm else lres end
neuper@37906
  1335
  | get_loc pt (p,_) =
neuper@37906
  1336
  let val (lfrm,lres) = get_obj g_loc pt p
neuper@37906
  1337
  in if lfrm = e_istate then lres else lfrm end;  5.10.00: too liberal ?*)
neuper@37906
  1338
(*13.8.02: options, because istate is no equalitype any more*)
neuper@37906
  1339
fun get_loc EmptyPtree _ = e_istate
neuper@37906
  1340
  | get_loc pt (p,Res) =
neuper@37906
  1341
    (case get_obj g_loc pt p of
neuper@37925
  1342
	 (SOME i, NONE) => i
neuper@37925
  1343
       | (NONE  , NONE) => e_istate
neuper@37925
  1344
       | (_     , SOME i) => i)
neuper@37906
  1345
  | get_loc pt (p,_) =
neuper@37906
  1346
    (case get_obj g_loc pt p of
neuper@37925
  1347
	 (NONE  , SOME i) => i (*13.8.02 just copied from ^^^: too liberal ?*)
neuper@37925
  1348
       | (NONE  , NONE) => e_istate
neuper@37925
  1349
       | (SOME i, _) => i);
neuper@37906
  1350
val get_istate = get_loc; (*3.5.02*)
neuper@37906
  1351
neuper@37906
  1352
(*.collect the assumptions within a problem up to a certain position.*)
neuper@37906
  1353
type asms = (term * pos) list;(*WN0502 should be (pos' * term) list
neuper@37906
  1354
				       ...........===^===*)
neuper@37906
  1355
neuper@37906
  1356
fun get_asm (b:pos, p:pos) (Nd (PblObj {result=(_,asm),...},_)) = 
neuper@38015
  1357
    ((*tracing ("### get_asm PblObj:(b,p)= "^
neuper@37906
  1358
		(pair2str(ints2str b, ints2str p)));*)
neuper@37906
  1359
     (map (rpair b) asm):asms)
neuper@37906
  1360
  | get_asm (b, p) (Nd (PrfObj {result=(_,asm),...}, [])) = 
neuper@38015
  1361
    ((*tracing ("### get_asm PrfObj []:(b,p)= "^
neuper@37906
  1362
	      (pair2str(ints2str b, ints2str p)));*)
neuper@37906
  1363
     (map (rpair b) asm))
neuper@37906
  1364
  | get_asm (b, p:pos) (Nd (PrfObj _, nds)) = 
neuper@38015
  1365
    let (*val _= tracing ("### get_asm PrfObj nds:(b,p)= "^
neuper@37906
  1366
	      (pair2str(ints2str b, ints2str p)));*)
neuper@37906
  1367
	val levdn = 
neuper@37906
  1368
	    if p <> [] then (b @ [hd p]:pos, tl p:pos) 
neuper@37906
  1369
	    else (b @ [1], [99999]) (*_deeper_ nesting is always _before_ p*)
neuper@37906
  1370
    in gets_asm levdn 1 nds end
neuper@37906
  1371
and gets_asm _ _ [] = []
neuper@37906
  1372
  | gets_asm (b, p' as p::ps) i (nd::nds) = 
neuper@37906
  1373
    if p < i then [] 
neuper@38015
  1374
    else ((*tracing ("### gets_asm: (b,p')= "^(pair2str(ints2str b,
neuper@37906
  1375
						      ints2str p')));*)
neuper@37906
  1376
	  (get_asm (b @ [i], ps) nd) @ (gets_asm (b, p') (i + 1) nds));
neuper@37906
  1377
neuper@37906
  1378
fun get_assumptions_ (Nd (PblObj {result=(r,asm),...}, cn)) (([], _):pos') = 
neuper@37906
  1379
    if r = e_term then gets_asm ([], [99999]) 1 cn
neuper@37906
  1380
    else map (rpair []) asm
neuper@37906
  1381
  | get_assumptions_ pt (p,p_) =
neuper@37906
  1382
    let val (cn, base) = par_children pt p
neuper@37906
  1383
	val offset = drop (length base, p)
neuper@37906
  1384
	val base' = replicate (length base) 1
neuper@37906
  1385
	val offset' = case p_ of 
neuper@37906
  1386
			 Frm => let val (qs,q) = split_last offset
neuper@37906
  1387
				in qs @ [q - 1] end
neuper@37906
  1388
		       | _ => offset
neuper@38015
  1389
        (*val _= tracing ("... get_assumptions: (b,o)= "^
neuper@37906
  1390
			(pair2str(ints2str base',ints2str offset)))*)
neuper@37906
  1391
    in gets_asm (base', offset) 1 cn end;
neuper@37906
  1392
neuper@37906
  1393
neuper@37906
  1394
(*---------
neuper@37906
  1395
end
neuper@37906
  1396
neuper@37906
  1397
open Ptree;
neuper@37906
  1398
----------*)
neuper@37906
  1399
neuper@37906
  1400
(*pos of the formula on FE relative to the current pos,
neuper@37906
  1401
  which is the next writepos*)
neuper@37906
  1402
fun pre_pos ([]:pos) = []:pos
neuper@37906
  1403
  | pre_pos pp =
neuper@37906
  1404
  let val (ps,p) = split_last pp
neuper@37906
  1405
  in case p of 1 => ps | n => ps @ [n-1] end;
neuper@37906
  1406
neuper@37906
  1407
(*WN.20.5.03 ... but not used*)
neuper@37906
  1408
fun posless [] (_::_) = true
neuper@37906
  1409
  | posless (_::_) [] = false
neuper@37906
  1410
  | posless (p::ps) (q::qs) = if p = q then posless ps qs else p < q;
neuper@37906
  1411
(* posless [2,3,4] [3,4,5];
neuper@37906
  1412
true
neuper@37906
  1413
>  posless [2,3,4] [1,2,3];
neuper@37906
  1414
false
neuper@37906
  1415
>  posless [2,3] [2,3,4];
neuper@37906
  1416
true
neuper@37906
  1417
>  posless [2,3,4] [2,3];
neuper@37906
  1418
false                    
neuper@37906
  1419
>  posless [6] [6,5,2];
neuper@37906
  1420
true
neuper@37906
  1421
+++ see Isabelle/../library.ML*)
neuper@37906
  1422
neuper@37906
  1423
neuper@37906
  1424
(**.development for extracting an 'interval' from ptree.**)
neuper@37906
  1425
neuper@37906
  1426
(*version 1 stopped 8.03 in favour of get_interval with !!!move_dn
neuper@37906
  1427
  actually used (inefficient) version with move_dn: see modspec.sml*)
neuper@37906
  1428
local
neuper@37906
  1429
neuper@37906
  1430
fun hdp [] = 1     | hdp [0] = 1     | hdp x = hd x;(*start with first*)
neuper@37906
  1431
fun hdq	[] = 99999 | hdq [0] = 99999 | hdq x = hd x;(*take until last*)
neuper@37906
  1432
fun tlp [] = [0]     | tlp [_] = [0]     | tlp x = tl x;
neuper@37906
  1433
fun tlq [] = [99999] | tlq [_] = [99999] | tlq x = tl x;
neuper@37906
  1434
neuper@37906
  1435
fun getnd i (b,p) q (Nd (po, nds)) =
neuper@37906
  1436
    (if  i <= 0 then [[b]] else []) @
neuper@37906
  1437
    (getnds (i-1) true (b@[hdp p], tlp p) (tlq q)
neuper@37906
  1438
	   (take_fromto (hdp p) (hdq q) nds))
neuper@37906
  1439
neuper@37906
  1440
and getnds _ _ _ _ [] = []                         (*no children*)
neuper@37906
  1441
  | getnds i _ (b,p) q [nd] = (getnd i (b,p) q nd) (*l+r-margin*)
neuper@37906
  1442
neuper@37906
  1443
  | getnds i true (b,p) q [n1, n2] =               (*l-margin,  r-margin*)
neuper@37906
  1444
    (getnd i      (       b, p ) [99999] n1) @
neuper@37906
  1445
    (getnd ~99999 (lev_on b,[0]) q       n2)
neuper@37906
  1446
neuper@37906
  1447
  | getnds i _    (b,p) q [n1, n2] =               (*intern,  r-margin*)
neuper@37906
  1448
    (getnd i      (       b,[0]) [99999] n1) @
neuper@37906
  1449
    (getnd ~99999 (lev_on b,[0]) q       n2)
neuper@37906
  1450
neuper@37906
  1451
  | getnds i true (b,p) q (nd::(nds as _::_)) =    (*l-margin, intern*)
neuper@37906
  1452
    (getnd i             (       b, p ) [99999] nd) @
neuper@37906
  1453
    (getnds ~99999 false (lev_on b,[0]) q nds)
neuper@37906
  1454
neuper@37906
  1455
  | getnds i _ (b,p) q (nd::(nds as _::_)) =       (*intern, ...*)
neuper@37906
  1456
    (getnd i             (       b,[0]) [99999] nd) @
neuper@37906
  1457
    (getnds ~99999 false (lev_on b,[0]) q nds); 
neuper@37906
  1458
in
neuper@37906
  1459
(*get an 'interval from to' from a ptree as 'intervals f t' of respective nodes
neuper@37906
  1460
  where 'from' are pos, i.e. a key as int list, 'f' an int (to,t analoguous)
neuper@37906
  1461
(1) the 'f' are given 
neuper@37906
  1462
(1a) by 'from' if 'f' = the respective element of 'from' (left margin)
neuper@37906
  1463
(1b) -inifinity, if 'f' > the respective element of 'from' (internal node)
neuper@37906
  1464
(2) the 't' ar given
neuper@37906
  1465
(2a) by 'to' if 't' = the respective element of 'to' (right margin)
neuper@37906
  1466
(2b) inifinity, if 't' < the respective element of 'to (internal node)'
neuper@37906
  1467
the 'f' and 't' are set by hdp,... *)
neuper@37906
  1468
fun get_trace pt p q =
neuper@37906
  1469
    (flat o (getnds ((length p) -1) true ([hdp p], tlp p) (tlq q))) 
neuper@37906
  1470
	(take_fromto (hdp p) (hdq q) (children pt));
neuper@37906
  1471
end;
neuper@37906
  1472
(*WN0510 stoppde this development;
neuper@37906
  1473
 actually used (inefficient) version with move_dn: getFormulaeFromTo*)
neuper@37906
  1474
neuper@37906
  1475
neuper@37906
  1476
neuper@37906
  1477
neuper@37906
  1478
fun get_somespec ((dI,pI,mI):spec) ((dI',pI',mI'):spec) =
neuper@37906
  1479
    let val domID = if dI = e_domID
neuper@37906
  1480
		    then if dI' = e_domID 
neuper@38031
  1481
			 then error"pt_extract: no domID in probl,origin"
neuper@37906
  1482
			 else dI'
neuper@37906
  1483
		    else dI
neuper@37906
  1484
	val pblID = if pI = e_pblID
neuper@37906
  1485
		    then if pI' = e_pblID 
neuper@38031
  1486
			 then error"pt_extract: no pblID in probl,origin"
neuper@37906
  1487
			 else pI'
neuper@37906
  1488
		    else pI
neuper@37906
  1489
	val metID = if mI = e_metID
neuper@37906
  1490
		    then if pI' = e_metID 
neuper@38031
  1491
			 then error"pt_extract: no metID in probl,origin"
neuper@37906
  1492
			 else mI'
neuper@37906
  1493
		    else mI
neuper@37906
  1494
    in (domID, pblID, metID):spec end;
neuper@37906
  1495
fun get_somespec' ((dI,pI,mI):spec) ((dI',pI',mI'):spec) =
neuper@37906
  1496
    let val domID = if dI = e_domID then dI' else dI
neuper@37906
  1497
	val pblID = if pI = e_pblID then pI' else pI
neuper@37906
  1498
	val metID = if mI = e_metID then mI' else mI
neuper@37906
  1499
    in (domID, pblID, metID):spec end;
neuper@37906
  1500
neuper@37906
  1501
(*extract a formula or model from ptree for itms2itemppc or model2xml*)
neuper@37906
  1502
fun preconds2str bts = 
neuper@37906
  1503
    (strs2str o (map (linefeed o pair2str o
neuper@37906
  1504
		      (apsnd term2str) o 
neuper@37906
  1505
		      (apfst bool2str)))) bts;
neuper@37906
  1506
fun ocalhd2str ((b, p, hdf, itms, prec, spec):ocalhd) =
neuper@37906
  1507
    "("^bool2str b^", "^pos_2str p^", "^term2str hdf^
neuper@37928
  1508
    ", "^itms2str_ (thy2ctxt' "Isac") itms^
neuper@37906
  1509
    ", "^preconds2str prec^", \n"^spec2str spec^" )";
neuper@37906
  1510
neuper@37906
  1511
neuper@37906
  1512
neuper@37906
  1513
fun is_pblnd (Nd (ppobj, _)) = is_pblobj ppobj;
neuper@37906
  1514
neuper@37906
  1515
neuper@37906
  1516
(**.functions for the 'ptree iterator' as seen from the FE-Kernel interface.**)
neuper@37906
  1517
neuper@37906
  1518
(*move one step down into existing nodes of ptree; regard TransitiveB
neuper@37906
  1519
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~##################
neuper@37906
  1520
fun move_dn _ (Nd (c, ns)) ([],p_) = (*root problem*)
neuper@37906
  1521
(* val (Nd (c, ns), ([],p_)) = (pt, get_pos cI uI);
neuper@37906
  1522
   *)
neuper@37906
  1523
    if is_pblobj c 
neuper@37906
  1524
    then case p_ of (*Frm => ([], Pbl) 1.12.03
neuper@37906
  1525
		  |*) Res => raise PTREE "move_dn: end of calculation"
neuper@37906
  1526
		  | _ => if null ns (*go down from Pbl + Met*)
neuper@37906
  1527
			 then raise PTREE "move_dn: solve problem not started"
neuper@37906
  1528
			 else ([1], Frm)
neuper@37906
  1529
    else (case p_ of Res => raise PTREE "move_dn: end of (sub-)tree"
neuper@37906
  1530
		  | _ => if null ns
neuper@37906
  1531
			 then raise PTREE "move_dn: pos not existent 1"
neuper@37906
  1532
			 else ([1], Frm))
neuper@37906
  1533
neuper@37906
  1534
  (*iterate towards end of pos*)
neuper@37906
  1535
(* val (P,(Nd (_, ns)),(p::(ps as (_::_)),p_)) = ([]:pos, pt, get_pos cI uI);
neuper@37906
  1536
   val (P,(Nd (_, ns)),(p::(ps as (_::_)),p_)) = ((P@[p]),(nth p ns),(ps, p_));
neuper@37906
  1537
   *) 
neuper@37906
  1538
 | move_dn P  (Nd (_, ns)) (p::(ps as (_::_)),p_) =
neuper@37906
  1539
    if p > length ns then raise PTREE "move_dn: pos not existent 2"
neuper@37906
  1540
    else move_dn ((P@[p]): pos) (nth p ns) (ps, p_)
neuper@37906
  1541
(* val (P, (Nd (c, ns)), ([p], p_)) = ((P@[p]), (nth p ns), (ps, p_));
neuper@37906
  1542
   val (P, (Nd (c, ns)), ([p], p_)) = ([],pt,get_pos cI uI);
neuper@37906
  1543
   *)
neuper@37906
  1544
  | move_dn P (Nd (c, ns)) ([p], p_) = (*act on last element of pos*)
neuper@37906
  1545
    if p > length ns then raise PTREE "move_dn: pos not existent 3"
neuper@37906
  1546
    else if is_pblnd (nth p ns)  then
neuper@38015
  1547
	((*tracing("### move_dn: is_pblnd (nth p ns), P= "^ints2str' P^", \n"^
neuper@37906
  1548
		 "length ns= "^((string_of_int o length) ns)^
neuper@37906
  1549
		 ", p= "^string_of_int p^", p_= "^pos_2str p_);*)
neuper@37906
  1550
	 case p_ of Res => if p = length ns 
neuper@37906
  1551
			   then if g_ostate c = Complete then (P, Res)
neuper@37906
  1552
				else raise PTREE (ints2str' P^" not complete")
neuper@37906
  1553
			   (*FIXME here handle not-sequent-branches*)
neuper@37906
  1554
			   else if g_branch c = TransitiveB 
neuper@37906
  1555
				   andalso (not o is_pblnd o (nth (p+1))) ns
neuper@37906
  1556
			   then (P@[p+1], Res)
neuper@37906
  1557
			   else (P@[p+1], if is_pblnd (nth (p+1) ns) 
neuper@37906
  1558
					  then Pbl else Frm)
neuper@37906
  1559
		  | _ => if (null o children o (nth p)) ns (*go down from Pbl*)
neuper@37906
  1560
			 then raise PTREE "move_dn: solve subproblem not started"
neuper@37906
  1561
			 else (P @ [p, 1], 
neuper@37906
  1562
			       if (is_pblnd o hd o children o (nth p)) ns
neuper@37906
  1563
			       then Pbl else Frm)
neuper@37906
  1564
			      )
neuper@37906
  1565
    (* val (P, Nd (c, ns), ([p], p_)) = ([], pt, ([1], Frm));
neuper@37906
  1566
        *)
neuper@37906
  1567
    else case p_ of Frm => if (null o children o (nth p)) ns 
neuper@37906
  1568
			 (*then if g_ostate c = Complete then (P@[p],Res)*)
neuper@37906
  1569
			   then if g_ostate' (nth p ns) = Complete 
neuper@37906
  1570
				then (P@[p],Res)
neuper@37906
  1571
				else raise PTREE "move_dn: pos not existent 4"
neuper@37906
  1572
			   else (P @ [p, 1], (*go down*) 
neuper@37906
  1573
				 if (is_pblnd o hd o children o (nth p)) ns
neuper@37906
  1574
				 then Pbl else Frm)
neuper@37906
  1575
		  | Res => if p = length ns 
neuper@37906
  1576
			   then 
neuper@37906
  1577
			      if g_ostate c = Complete then (P, Res)
neuper@37906
  1578
			      else raise PTREE (ints2str' P^" not complete")
neuper@37906
  1579
			   else 
neuper@37906
  1580
			       if g_branch c = TransitiveB 
neuper@37906
  1581
				  andalso (not o is_pblnd o (nth (p+1))) ns
neuper@37906
  1582
			       then if (null o children o (nth (p+1))) ns
neuper@37906
  1583
				    then (P@[p+1], Res)
neuper@37906
  1584
				    else (P@[p+1,1], Frm)(*040221*)
neuper@37906
  1585
			       else (P@[p+1], if is_pblnd (nth (p+1) ns) 
neuper@37906
  1586
					      then Pbl else Frm); 
neuper@37906
  1587
*)
neuper@37906
  1588
(*.move one step down into existing nodes of ptree; skip Res = Frm.nxt;
neuper@37906
  1589
   move_dn at the end of the calc-tree raises PTREE.*)
neuper@37906
  1590
fun move_dn _ (Nd (c, ns)) ([],p_) = (*root problem*)
neuper@37906
  1591
    (case p_ of 
neuper@37906
  1592
	     Res => raise PTREE "move_dn: end of calculation"
neuper@37906
  1593
	   | _ => if null ns (*go down from Pbl + Met*)
neuper@37906
  1594
		  then raise PTREE "move_dn: solve problem not started"
neuper@37906
  1595
		  else ([1], Frm))
neuper@37906
  1596
  | move_dn P  (Nd (_, ns)) (p::(ps as (_::_)),p_) =(*iterate to end of pos*)
neuper@37906
  1597
    if p > length ns then raise PTREE "move_dn: pos not existent 2"
neuper@37906
  1598
    else move_dn ((P@[p]): pos) (nth p ns) (ps, p_)
neuper@37906
  1599
neuper@37906
  1600
  | move_dn P (Nd (c, ns)) ([p], p_) = (*act on last element of pos*)
neuper@37906
  1601
    if p > length ns then raise PTREE "move_dn: pos not existent 3"
neuper@37906
  1602
    else case p_ of 
neuper@37906
  1603
	     Res => 
neuper@37906
  1604
	     if p = length ns (*last Res on this level: go a level up*)
neuper@37906
  1605
	     then if g_ostate c = Complete then (P, Res)
neuper@37906
  1606
		  else raise PTREE (ints2str' P^" not complete 1")
neuper@37906
  1607
	     else (*go to the next Nd on this level, or down into the next Nd*)
neuper@37906
  1608
		 if is_pblnd (nth (p+1) ns) then (P@[p+1], Pbl)
neuper@37906
  1609
		 else 
neuper@37906
  1610
		     if g_res' (nth p ns) = g_form' (nth (p+1) ns)
neuper@37906
  1611
		     then if (null o children o (nth (p+1))) ns
neuper@37906
  1612
			  then (*take the Res if Complete*) 
neuper@37906
  1613
			      if g_ostate' (nth (p+1) ns) = Complete 
neuper@37906
  1614
			      then (P@[p+1], Res)
neuper@37906
  1615
			      else raise PTREE (ints2str' (P@[p+1])^
neuper@37906
  1616
						" not complete 2")
neuper@37906
  1617
			  else (P@[p+1,1], Frm)(*go down into the next PrfObj*)
neuper@37906
  1618
		     else (P@[p+1], Frm)(*take Frm: exists if the Nd exists*)
neuper@37906
  1619
	   | Frm => (*go down or to the Res of this Nd*)
neuper@37906
  1620
	     if (null o children o (nth p)) ns
neuper@37906
  1621
	     then if g_ostate' (nth p ns) = Complete then (P @ [p], Res)
neuper@37906
  1622
		  else raise PTREE (ints2str' (P @ [p])^" not complete 3")
neuper@37906
  1623
	     else (P @ [p, 1], Frm)
neuper@37906
  1624
	   | _ => (*is Pbl or Met*)
neuper@37906
  1625
	     if (null o children o (nth p)) ns
neuper@37906
  1626
	     then raise PTREE "move_dn:solve subproblem not startd"
neuper@37906
  1627
	     else (P @ [p, 1], 
neuper@37906
  1628
		   if (is_pblnd o hd o children o (nth p)) ns
neuper@37906
  1629
		   then Pbl else Frm);
neuper@37906
  1630
neuper@37906
  1631
neuper@37906
  1632
(*.go one level down into ptree.*)
neuper@37906
  1633
fun movelevel_dn [] (Nd (c, ns)) ([],p_) = (*root problem*)
neuper@37906
  1634
    if is_pblobj c 
neuper@37906
  1635
    then if null ns 
neuper@37906
  1636
	 then raise PTREE "solve problem not started"
neuper@37906
  1637
	 else ([1], if (is_pblnd o hd) ns then Pbl else Frm)
neuper@37906
  1638
    else raise PTREE "pos not existent 1"
neuper@37906
  1639
neuper@37906
  1640
  (*iterate towards end of pos*)
neuper@37906
  1641
  | movelevel_dn P (Nd (_, ns)) (p::(ps as (_::_)),p_) =
neuper@37906
  1642
    if p > length ns then raise PTREE "pos not existent 2"
neuper@37906
  1643
    else movelevel_dn (P@[p]) (nth p ns) (ps, p_)
neuper@37906
  1644
neuper@37906
  1645
  | movelevel_dn P (Nd (c, ns)) ([p], p_) = (*act on last element of pos*)
neuper@37906
  1646
    if p > length ns then raise PTREE "pos not existent 3" else
neuper@37906
  1647
    case p_ of Res => 
neuper@37906
  1648
	       if p = length ns 
neuper@37906
  1649
	       then raise PTREE "no children"
neuper@37906
  1650
	       else 
neuper@37906
  1651
		   if g_branch c = TransitiveB
neuper@37906
  1652
		   then if (null o children o (nth (p+1))) ns
neuper@37906
  1653
			then raise PTREE "no children"
neuper@37906
  1654
			else (P @ [p+1, 1], 
neuper@37906
  1655
			      if (is_pblnd o hd o children o (nth (p+1))) ns
neuper@37906
  1656
			      then Pbl else Frm)
neuper@37906
  1657
		   else if (null o children o (nth p)) ns
neuper@37906
  1658
		   then raise PTREE "no children"
neuper@37906
  1659
		   else (P @ [p, 1], if (is_pblnd o hd o children o (nth p)) ns
neuper@37906
  1660
				     then Pbl else Frm)
neuper@37906
  1661
	     | _ => if (null o children o (nth p)) ns 
neuper@37906
  1662
		    then raise PTREE "no children"
neuper@37906
  1663
		    else (P @ [p, 1], (*go down*)
neuper@37906
  1664
			  if (is_pblnd o hd o children o (nth p)) ns
neuper@37906
  1665
			  then Pbl else Frm);
neuper@37906
  1666
neuper@37906
  1667
neuper@37906
  1668
neuper@37906
  1669
(*.go to the previous position in ptree; regard TransitiveB.*)
neuper@37906
  1670
fun move_up _ (Nd (c, ns)) (([],p_):pos') = (*root problem*)
neuper@37906
  1671
    if is_pblobj c 
neuper@37906
  1672
    then case p_ of Res => if null ns then ([], Pbl) (*Res -> Pbl (not Met)!*)
neuper@37906
  1673
			   else ([length ns], Res)
neuper@37906
  1674
		  | _  => raise PTREE "begin of calculation"
neuper@37906
  1675
    else raise PTREE "pos not existent"
neuper@37906
  1676
neuper@37906
  1677
  | move_up P  (Nd (_, ns)) (p::(ps as (_::_)),p_) = (*iterate to end of pos*)
neuper@37906
  1678
    if p > length ns then raise PTREE "pos not existent"
neuper@37906
  1679
    else move_up (P@[p]) (nth p ns) (ps,p_)
neuper@37906
  1680
neuper@37906
  1681
  | move_up P (Nd (c, ns)) ([p], p_) = (*act on last element of pos*)
neuper@37906
  1682
    if p > length ns then raise PTREE "pos not existent"
neuper@37906
  1683
    else if is_pblnd (nth p ns)  then
neuper@37906
  1684
	case p_ of Res => 
neuper@37906
  1685
		   let val nc = (length o children o (nth p)) ns
neuper@37906
  1686
		   in if nc = 0 then (P@[p], Pbl) (*Res -> Pbl (not Met)!*)
neuper@37906
  1687
		      else (P @ [p, nc], Res) end (*go down*)
neuper@37906
  1688
		 | _ => if p = 1 then (P, Pbl) else (P@[p-1], Res) 
neuper@37906
  1689
    else case p_ of Frm => if p <> 1 then (P, Frm) 
neuper@37906
  1690
			  else if is_pblobj c then (P, Pbl) else (P, Frm)
neuper@37906
  1691
		  | Res => 
neuper@37906
  1692
		    let val nc = (length o children o (nth p)) ns
neuper@37906
  1693
		    in if nc = 0 (*cannot go down*)
neuper@37906
  1694
		       then if g_branch c = TransitiveB andalso p <> 1
neuper@37906
  1695
			    then (P@[p-1], Res) else (P@[p], Frm)
neuper@37906
  1696
		       else (P @ [p, nc], Res) end; (*go down*)
neuper@37906
  1697
neuper@37906
  1698
neuper@37906
  1699
neuper@37906
  1700
(*.go one level up in ptree; sets the position on Frm.*)
neuper@37906
  1701
fun movelevel_up _ (Nd (c, ns)) (([],p_):pos') = (*root problem*)
neuper@37906
  1702
    raise PTREE "pos not existent"
neuper@37906
  1703
neuper@37906
  1704
  (*iterate towards end of pos*)
neuper@37906
  1705
  | movelevel_up P  (Nd (_, ns)) (p::(ps as (_::_)),p_) = 
neuper@37906
  1706
    if p > length ns then raise PTREE "pos not existent"
neuper@37906
  1707
    else movelevel_up (P@[p]) (nth p ns) (ps,p_)
neuper@37906
  1708
neuper@37906
  1709
  | movelevel_up P (Nd (c, ns)) ([p], p_) = (*act on last element of pos*)
neuper@37906
  1710
    if p > length ns then raise PTREE "pos not existent"
neuper@37906
  1711
    else if is_pblobj c then (P, Pbl) else (P, Frm);
neuper@37906
  1712
neuper@37906
  1713
neuper@37906
  1714
(*.go to the next calc-head up in the calc-tree.*)
neuper@37906
  1715
fun movecalchd_up pt ((p, Res):pos') =
neuper@37906
  1716
    (par_pblobj pt p, Pbl):pos'
neuper@37906
  1717
  | movecalchd_up pt (p, _) =
neuper@37906
  1718
    if is_pblobj (get_obj I pt p) 
neuper@37906
  1719
    then (p, Pbl) else (par_pblobj pt p, Pbl);
neuper@37906
  1720
neuper@37906
  1721
(*.determine the previous pos' on the same level.*)
neuper@37906
  1722
(*WN0502 made for interSteps; _only_ works for branch TransitiveB*)
neuper@37906
  1723
fun lev_pred' pt (pos:pos' as ([],Res)) = ([],Pbl):pos'
neuper@37906
  1724
  | lev_pred' pt (pos:pos' as (p, Res)) =
neuper@37906
  1725
    let val (p', last) = split_last p
neuper@37906
  1726
    in if last = 1 
neuper@37906
  1727
       then if (is_pblobj o (get_obj I pt)) p then (p,Pbl) else (p, Frm)
neuper@37906
  1728
       else if get_obj g_res pt (p' @ [last - 1]) = get_obj g_form pt p
neuper@37906
  1729
       then (p' @ [last - 1], Res) (*TransitiveB*)
neuper@37906
  1730
       else if (is_pblobj o (get_obj I pt)) p then (p,Pbl) else (p, Frm)
neuper@37906
  1731
    end;
neuper@37906
  1732
neuper@37906
  1733
(*.determine the next pos' on the same level.*)
neuper@37906
  1734
fun lev_on' pt (([],Pbl):pos') = ([],Res):pos'
neuper@37906
  1735
  | lev_on' pt (p, Res) =
neuper@37906
  1736
    if get_obj g_res pt p = get_obj g_form pt (lev_on p)(*TransitiveB*)
neuper@37906
  1737
    then if existpt' (lev_on p, Res) pt then (lev_on p, Res)
neuper@38031
  1738
	 else error ("lev_on': (p, Res) -> (p, Res) not existent, \
neuper@37906
  1739
		      \p = "^ints2str' (lev_on p))
neuper@37906
  1740
    else (lev_on p, Frm)
neuper@37906
  1741
  | lev_on' pt (p, _) =
neuper@37906
  1742
    if existpt' (p, Res) pt then (p, Res)
neuper@38031
  1743
    else error ("lev_on': (p, Frm) -> (p, Res) not existent, \
neuper@37906
  1744
		      \p = "^ints2str' p);
neuper@37906
  1745
neuper@37906
  1746
fun exist_lev_on' pt p = (lev_on' pt p; true) handle _ => false;
neuper@37906
  1747
neuper@37906
  1748
(*.is the pos' at the last element of a calulation _AND_ can be continued.*)
neuper@37906
  1749
(* val (pt, pos as (p,p_)) = (pt, ([1],Frm));
neuper@37906
  1750
   *)
neuper@37906
  1751
fun is_curr_endof_calc pt (([],Res) : pos') = false
neuper@37906
  1752
  | is_curr_endof_calc pt (pos as (p,_)) =
neuper@37906
  1753
    not (exist_lev_on' pt pos) 
neuper@37906
  1754
    andalso get_obj g_ostate pt (lev_up p) = Incomplete;
neuper@37906
  1755
neuper@37906
  1756
neuper@37906
  1757
(**.insert into ctree and cut branches accordingly.**)
neuper@37906
  1758
  
neuper@37906
  1759
(*.get all positions of certain intervals on the ctree.*)
neuper@37906
  1760
(*OLD VERSION without move_dn; kept for occasional redesign
neuper@37906
  1761
   get all pos's to be cut in a ptree
neuper@37906
  1762
   below a pos or from a ptree list after i-th element (NO level_up).*)
neuper@37906
  1763
fun get_allpos' (_:pos, _:posel) EmptyPtree   = ([]:pos' list)
neuper@37906
  1764
  | get_allpos' (p, 1) (Nd (b, bs)) = (*p is pos of Nd*)
neuper@37906
  1765
    if g_ostate b = Incomplete 
neuper@38015
  1766
    then ((*tracing("get_allpos' (p, 1) Incomplete: p="^ints2str' p);*)
neuper@37906
  1767
	  [(p,Frm)] @ (get_allpos's (p, 1) bs)
neuper@37906
  1768
	  )
neuper@38015
  1769
    else ((*tracing("get_allpos' (p, 1) else: p="^ints2str' p);*)
neuper@37906
  1770
	  [(p,Frm)] @ (get_allpos's (p, 1) bs) @ [(p,Res)]
neuper@37906
  1771
	  )
neuper@37906
  1772
    (*WN041020 here we assume what is presented on the worksheet ?!*)
neuper@37906
  1773
  | get_allpos' (p, i) (Nd (b, bs)) = (*p is pos of Nd*)
neuper@37906
  1774
    if length bs > 0 orelse is_pblobj b
neuper@37906
  1775
    then if g_ostate b = Incomplete 
neuper@37906
  1776
	 then [(p,Frm)] @ (get_allpos's (p, 1) bs)
neuper@37906
  1777
	 else [(p,Frm)] @ (get_allpos's (p, 1) bs) @ [(p,Res)]
neuper@37906
  1778
    else 
neuper@37906
  1779
	if g_ostate b = Incomplete 
neuper@37906
  1780
	then []
neuper@37906
  1781
	else [(p,Res)]
neuper@37906
  1782
(*WN041020 here we assume what is presented on the worksheet ?!*)
neuper@37906
  1783
and get_allpos's _ [] = []
neuper@37906
  1784
  | get_allpos's (p, i) (pt::pts) = (*p is pos of parent-Nd*)
neuper@37906
  1785
    (get_allpos' (p@[i], i) pt) @ (get_allpos's (p, i+1) pts);
neuper@37906
  1786
neuper@37906
  1787
(*.get all positions of certain intervals on the ctree.*)
neuper@37906
  1788
(*NEW version WN050225*)
neuper@37906
  1789
neuper@37906
  1790
neuper@37906
  1791
(*.cut branches.*)
neuper@37906
  1792
(*before WN041019......
neuper@37906
  1793
val cut_branch = (test_trans, curry take):
neuper@37906
  1794
    (ppobj -> bool) * (int -> ptree list -> ptree list);
neuper@37906
  1795
.. formlery used for ...
neuper@37906
  1796
fun cut_tree''' _ [] = EmptyPtree
neuper@37906
  1797
  | cut_tree''' pt pos = 
neuper@37906
  1798
  let val (pt',cut) = appl_branch cut_branch pt pos
neuper@37906
  1799
  in if cut andalso length pos > 1 then cut_tree''' pt' (lev_up pos)
neuper@37906
  1800
     else pt' end;
neuper@37906
  1801
*)
neuper@37906
  1802
(*OLD version before WN050225*)
neuper@37906
  1803
(*WN050106 like cut_level, but deletes exactly 1 node --- for tests ONLY*)
neuper@37906
  1804
fun cut_level_'_ (_:pos' list) (_:pos) EmptyPtree (_:pos') =
neuper@37906
  1805
    raise PTREE "cut_level_'_ Empty _"
neuper@37906
  1806
  | cut_level_'_ _ _ (Nd ( _, _)) ([],_) = raise PTREE "cut_level_'_ _ []"
neuper@37906
  1807
  | cut_level_'_ cuts P (Nd (b, bs)) (p::[],p_) = 
neuper@37906
  1808
    if test_trans b 
neuper@37906
  1809
    then (Nd (b, drop_nth [] (p:posel, bs)),
neuper@37906
  1810
	  (*     ~~~~~~~~~~~*)
neuper@37906
  1811
	  cuts @ 
neuper@37906
  1812
	  (if p_ = Frm then [(P@[p],Res)] else ([]:pos' list)) @
neuper@37906
  1813
	  (*WN041020 here we assume what is presented on the worksheet ?!*)
neuper@37906
  1814
	  (get_allpos's (P, p+1) (drop_nth [] (p, bs))))
neuper@37906
  1815
    (*                            ~~~~~~~~~~~*)
neuper@37906
  1816
    else (Nd (b, bs), cuts)
neuper@37906
  1817
  | cut_level_'_ cuts P (Nd (b, bs)) ((p::ps),p_) =
neuper@37906
  1818
    let val (bs',cuts') = cut_level_'_ cuts P (nth p bs) (ps, p_)
neuper@37906
  1819
    in (Nd (b, repl_app bs p bs'), cuts @ cuts') end;
neuper@37906
  1820
neuper@37906
  1821
(*before WN050219*)
neuper@37906
  1822
fun cut_level (_:pos' list) (_:pos) EmptyPtree (_:pos') =
neuper@37906
  1823
    raise PTREE "cut_level EmptyPtree _"
neuper@37906
  1824
  | cut_level _ _ (Nd ( _, _)) ([],_) = raise PTREE "cut_level _ []"
neuper@37906
  1825
neuper@37906
  1826
  | cut_level cuts P (Nd (b, bs)) (p::[],p_) = 
neuper@37906
  1827
    if test_trans b 
neuper@37906
  1828
    then (Nd (b, take (p:posel, bs)),
neuper@37906
  1829
	  cuts @ 
neuper@37906
  1830
	  (if p_ = Frm andalso (*#*) g_ostate b = Complete
neuper@37906
  1831
	   then [(P@[p],Res)] else ([]:pos' list)) @
neuper@37906
  1832
	  (*WN041020 here we assume what is presented on the worksheet ?!*)
neuper@37906
  1833
	  (get_allpos's (P, p+1) (takerest (p, bs))))
neuper@37906
  1834
    else (Nd (b, bs), cuts)
neuper@37906
  1835
neuper@37906
  1836
  | cut_level cuts P (Nd (b, bs)) ((p::ps),p_) =
neuper@37906
  1837
    let val (bs',cuts') = cut_level cuts P (nth p bs) (ps, p_)
neuper@37906
  1838
    in (Nd (b, repl_app bs p bs'), cuts @ cuts') end;
neuper@37906
  1839
neuper@37906
  1840
(*OLD version before WN050219, overwritten below*)
neuper@37906
  1841
fun cut_tree _ (([],_):pos') = raise PTREE "cut_tree _ ([],_)"
neuper@37906
  1842
  | cut_tree pt (pos as ([p],_)) =
neuper@37906
  1843
    let	val (pt', cuts) = cut_level ([]:pos' list) [] pt pos
neuper@37906
  1844
    in (pt', cuts @ (if get_obj g_ostate pt [] = Incomplete 
neuper@37906
  1845
		     then [] else [([],Res)])) end
neuper@37906
  1846
  | cut_tree pt (p,p_) =
neuper@37906
  1847
    let	
neuper@37906
  1848
	fun cutfn pt cuts (p,p_) = 
neuper@37906
  1849
	    let val (pt', cuts') = cut_level [] (lev_up p) pt (p,p_)
neuper@37906
  1850
		val cuts'' = if get_obj g_ostate pt (lev_up p) = Incomplete 
neuper@37906
  1851
			     then [] else [(lev_up p, Res)]
neuper@37906
  1852
	    in if length cuts' > 0 andalso length p > 1
neuper@37906
  1853
	       then cutfn pt' (cuts @ cuts') (lev_up p, Frm(*-->(p,Res)*))
neuper@37906
  1854
	       else (pt',cuts @ cuts') end
neuper@37906
  1855
	val (pt', cuts) = cutfn pt [] (p,p_)
neuper@37906
  1856
    in (pt', cuts @ (if get_obj g_ostate pt [] = Incomplete 
neuper@37906
  1857
		     then [] else [([], Res)])) end;
neuper@37906
  1858
neuper@37906
  1859
neuper@37906
  1860
(*########/ inserted from ctreeNEW.sml \#################################**)
neuper@37906
  1861
neuper@37906
  1862
(*.get all positions in a ptree until ([],Res) or ostate=Incomplete
neuper@37906
  1863
val get_allp = fn : 
neuper@37906
  1864
  pos' list -> : accumulated, start with []
neuper@37906
  1865
  pos ->       : the offset for subtrees wrt the root
neuper@37906
  1866
  ptree ->     : (sub)tree
neuper@37906
  1867
  pos'         : initialization (the last pos' before ...)
neuper@37906
  1868
  -> pos' list : of positions in this (sub) tree (relative to the root)
neuper@37906
  1869
.*)
neuper@37906
  1870
(* val (cuts, P, pt, pos) = ([], [3], get_nd pt [3], ([], Frm):pos');
neuper@37906
  1871
   val (cuts, P, pt, pos) = ([], [2], get_nd pt [2], ([], Frm):pos');
neuper@37906
  1872
   length (children pt);
neuper@37906
  1873
   *)
neuper@37906
  1874
fun get_allp (cuts:pos' list) (P:pos, pos:pos') pt =
neuper@37906
  1875
    (let val nxt = move_dn [] pt pos (*exn if Incomplete reached*)
neuper@37906
  1876
     in if nxt <> ([],Res) 
neuper@37906
  1877
	then get_allp (cuts @ [nxt]) (P, nxt) pt
neuper@37906
  1878
	else (map (apfst (curry op@ P)) (cuts @ [nxt])): pos' list
neuper@37906
  1879
     end) handle PTREE _ => (map (apfst (curry op@ P)) cuts);
neuper@37906
  1880
neuper@37906
  1881
neuper@37906
  1882
(*the pts are assumed to be on the same level*)
neuper@37906
  1883
fun get_allps (cuts: pos' list) (P:pos) [] = cuts
neuper@37906
  1884
  | get_allps cuts P (pt::pts) =
neuper@37906
  1885
    let val below = get_allp [] (P, ([], Frm)) pt
neuper@37906
  1886
	val levfrm = 
neuper@37906
  1887
	    if is_pblnd pt 
neuper@37906
  1888
	    then (P, Pbl)::below
neuper@37906
  1889
	    else if last_elem P = 1 
neuper@37906
  1890
	    then (P, Frm)::below
neuper@37906
  1891
	    else (*Trans*) below
neuper@37906
  1892
	val levres = levfrm @ (if null below then [(P, Res)] else [])
neuper@37906
  1893
    in get_allps (cuts @ levres) (lev_on P) pts end;
neuper@37906
  1894
neuper@37906
  1895
neuper@37906
  1896
(**.these 2 funs decide on how far cut_tree goes.**)
neuper@37906
  1897
(*.shall the nodes _after_ the pos to be inserted at be deleted?.*)
neuper@37906
  1898
fun test_trans (PrfObj{branch = Transitive,...}) = true
neuper@37906
  1899
  | test_trans (PrfObj{branch = NoBranch,...}) = true
neuper@37906
  1900
  | test_trans (PblObj{branch = Transitive,...}) = true 
neuper@37906
  1901
  | test_trans (PblObj{branch = NoBranch,...}) = true 
neuper@37906
  1902
  | test_trans _ = false;
neuper@37906
  1903
(*.shall cutting be continued on the higher level(s)?
neuper@37906
  1904
   the Nd regarded will NOT be changed.*)
neuper@37906
  1905
fun cutlevup (PblObj _) = false (*for tests of LK0502*)
neuper@37906
  1906
  | cutlevup _ = true;
neuper@37906
  1907
val cutlevup = test_trans;(*WN060727 after summerterm tests.LK0502 withdrawn*)
neuper@37906
  1908
    
neuper@37906
  1909
(*cut_bottom new sml603..608
neuper@37906
  1910
cut the level at the bottom of the pos (used by cappend_...)
neuper@37906
  1911
and handle the parent in order to avoid extra case for root
neuper@37906
  1912
fn: ptree ->         : the _whole_ ptree for cut_levup
neuper@37906
  1913
    pos * posel ->   : the pos after split_last
neuper@37906
  1914
    ptree ->         : the parent of the Nd to be cut
neuper@37906
  1915
return
neuper@37906
  1916
    (ptree *         : the updated ptree
neuper@37906
  1917
     pos' list) *    : the pos's cut
neuper@37906
  1918
     bool            : cutting shall be continued on the higher level(s)
neuper@37906
  1919
*)
neuper@37906
  1920
fun cut_bottom _ (pt' as Nd (b, [])) = ((pt', []), cutlevup b)
neuper@37906
  1921
  | cut_bottom (P:pos, p:posel) (Nd (b, bs)) =
neuper@37906
  1922
    let (*divide level into 3 parts...*)
neuper@37906
  1923
	val keep = take (p - 1, bs)
neuper@37906
  1924
	val pt' as Nd (_,bs') = nth p bs
neuper@37906
  1925
	(*^^^^^_here_ will be 'insert'ed by 'append_..'*)
neuper@37906
  1926
	val (tail, tp) = (takerest (p, bs), 
neuper@37906
  1927
			  if null (takerest (p, bs)) then 0 else p + 1)
neuper@37906
  1928
	val (children, cuts) = 
neuper@37906
  1929
	    if test_trans b
neuper@37906
  1930
	    then (keep,
neuper@37906
  1931
		  (if is_pblnd pt' then [(P @ [p], Pbl)] else [])
neuper@37906
  1932
		  @ (get_allp  [] (P @ [p], (P, Frm)) pt')
neuper@37906
  1933
		  @ (get_allps [] (P @ [p+1]) tail))
neuper@37906
  1934
	    else (keep @ [(*'insert'ed by 'append_..'*)] @ tail,
neuper@37906
  1935
		  get_allp  [] (P @ [p], (P, Frm)) pt')
neuper@37906
  1936
	val (pt'', cuts) = 
neuper@37906
  1937
	    if cutlevup b
neuper@37906
  1938
	    then (Nd (del_res b, children), 
neuper@37906
  1939
		  cuts @ (if g_ostate b = Incomplete then [] else [(P,Res)]))
neuper@37906
  1940
	    else (Nd (b, children), cuts)
neuper@38015
  1941
	(*val _= tracing("####cut_bottom (P, p)="^pos2str (P @ [p])^
neuper@37906
  1942
		       ", Nd=.............................................")
neuper@37906
  1943
	val _= show_pt pt''
neuper@38015
  1944
	val _= tracing("####cut_bottom form='"^
neuper@37906
  1945
		       term2str (get_obj g_form pt'' []))
neuper@38015
  1946
	val _= tracing("####cut_bottom cuts#="^string_of_int (length cuts)^
neuper@37906
  1947
		       ", cuts="^pos's2str cuts)*)
neuper@37906
  1948
    in ((pt'', cuts:pos' list), cutlevup b) end;
neuper@37906
  1949
neuper@37906
  1950
neuper@37906
  1951
(*.go all levels from the bottom of 'pos' up to the root, 
neuper@37906
  1952
 on each level compose the children of a node and accumulate the cut Nds
neuper@37906
  1953
args
neuper@37906
  1954
   pos' list ->      : for accumulation
neuper@37906
  1955
   bool -> 	     : cutting shall be continued on the higher level(s)
neuper@37906
  1956
   ptree -> 	     : the whole ptree for 'get_nd pt P' on each level
neuper@37906
  1957
   ptree -> 	     : the Nd from the lower level for insertion at path
neuper@37906
  1958
   pos * posel ->    : pos=path split for convenience
neuper@37906
  1959
   ptree -> 	     : Nd the children of are under consideration on this call 
neuper@37906
  1960
returns		     :
neuper@37906
  1961
   ptree * pos' list : the updated parent-Nd and the pos's of the Nds cut
neuper@37906
  1962
.*)
neuper@37906
  1963
fun cut_levup (cuts:pos' list) clevup pt pt' (P:pos, p:posel) (Nd (b, bs)) =
neuper@37906
  1964
    let (*divide level into 3 parts...*)
neuper@37906
  1965
	val keep = take (p - 1, bs)
neuper@37906
  1966
	(*val pt' comes as argument from below*)
neuper@37906
  1967
	val (tail, tp) = (takerest (p, bs), 
neuper@37906
  1968
			  if null (takerest (p, bs)) then 0 else p + 1)
neuper@37906
  1969
	val (children, cuts') = 
neuper@37906
  1970
	    if clevup
neuper@37906
  1971
	    then (keep @ [pt'], get_allps [] (P @ [p+1]) tail)
neuper@37906
  1972
	    else (keep @ [pt'] @ tail, [])
neuper@37906
  1973
	val clevup' = if clevup then cutlevup b else false 
neuper@37906
  1974
	(*the first Nd with false stops cutting on all levels above*)
neuper@37906
  1975
	val (pt'', cuts') = 
neuper@37906
  1976
	    if clevup'
neuper@37906
  1977
	    then (Nd (del_res b, children), 
neuper@37906
  1978
		  cuts' @ (if g_ostate b = Incomplete then [] else [(P,Res)]))
neuper@37906
  1979
	    else (Nd (b, children), cuts')
neuper@38015
  1980
	(*val _= tracing("#####cut_levup clevup= "^bool2str clevup)
neuper@38015
  1981
	val _= tracing("#####cut_levup cutlevup b= "^bool2str (cutlevup b))
neuper@38015
  1982
	val _= tracing("#####cut_levup (P, p)="^pos2str (P @ [p])^
neuper@37906
  1983
		       ", Nd=.............................................")
neuper@37906
  1984
	val _= show_pt pt''
neuper@38015
  1985
	val _= tracing("#####cut_levup form='"^
neuper@37906
  1986
		       term2str (get_obj g_form pt'' []))
neuper@38015
  1987
	val _= tracing("#####cut_levup cuts#="^string_of_int (length cuts)^
neuper@37906
  1988
		       ", cuts="^pos's2str cuts)*)
neuper@37906
  1989
    in if null P then (pt'', (cuts @ cuts'):pos' list)
neuper@37906
  1990
       else let val (P, p) = split_last P
neuper@37906
  1991
	    in cut_levup (cuts @ cuts') clevup' pt pt'' (P, p) (get_nd pt P)
neuper@37906
  1992
	    end
neuper@37906
  1993
    end;
neuper@37906
  1994
 
neuper@37906
  1995
(*.cut nodes after and below an inserted node in the ctree;
neuper@37906
  1996
   the cuts range is limited by the predicate 'fun cutlevup'.*)
neuper@37906
  1997
fun cut_tree pt (pos,_) =
neuper@37906
  1998
    if not (existpt pos pt) 
neuper@37906
  1999
    then (pt,[]) (*appending a formula never cuts anything*)
neuper@37906
  2000
    else let val (P, p) = split_last pos
neuper@37906
  2001
	     val ((pt', cuts), clevup) = cut_bottom (P, p) (get_nd pt P)
neuper@37906
  2002
	 (*        pt' is the updated parent of the Nd to cappend_..*)
neuper@37906
  2003
	 in if null P then (pt', cuts)
neuper@37906
  2004
	    else let val (P, p) = split_last P
neuper@37906
  2005
		 in cut_levup cuts clevup pt pt' (P, p) (get_nd pt P)
neuper@37906
  2006
		 end
neuper@37906
  2007
	 end;
neuper@37906
  2008
neuper@37906
  2009
fun append_atomic p l f r f' s pt = 
neuper@38015
  2010
  let (**val _= tracing("#@append_atomic: pos ="^pos2str p)**)
neuper@37906
  2011
	val (iss, f) = if existpt p pt andalso get_obj g_tac pt p=Empty_Tac
neuper@37906
  2012
		     then (*after Take*)
neuper@37925
  2013
			 ((fst (get_obj g_loc pt p), SOME l), 
neuper@37906
  2014
			  get_obj g_form pt p) 
neuper@37925
  2015
		     else ((NONE, SOME l), f)
neuper@37925
  2016
  in insert (PrfObj {cell = NONE,
neuper@37906
  2017
		     form  = f,
neuper@37906
  2018
		     tac  = r,
neuper@37906
  2019
		     loc   = iss,
neuper@37906
  2020
		     branch= NoBranch,
neuper@37906
  2021
		     result= f',
neuper@37906
  2022
		     ostate= s}) pt p end;
neuper@37906
  2023
neuper@37906
  2024
neuper@37906
  2025
(*20.8.02: cappend_* FIXXXXME cut branches below cannot be decided here:
neuper@37906
  2026
  detail - generate - cappend: inserted, not appended !!!
neuper@37906
  2027
neuper@37906
  2028
  cut decided in applicable_in !?!
neuper@37906
  2029
*)
neuper@37906
  2030
fun cappend_atomic pt p loc f r f' s = 
neuper@37906
  2031
(* val (pt, p, loc, f, r, f', s) = 
neuper@37906
  2032
       (pt,p,l,f,Rewrite_Set_Inst (subst2subs subs',id_rls rls'),
neuper@37906
  2033
	(f',asm),Complete);
neuper@37906
  2034
   *)
neuper@38015
  2035
((*tracing("##@cappend_atomic: pos ="^pos2str p);*)
neuper@37906
  2036
  apfst (append_atomic p loc f r f' s) (cut_tree pt (p,Frm))
neuper@37906
  2037
);
neuper@37906
  2038
(*TODO.WN050305 redesign the handling of istates*)
neuper@37906
  2039
fun cappend_atomic pt p ist_res f r f' s = 
neuper@37906
  2040
    if existpt p pt andalso get_obj g_tac pt p=Empty_Tac
neuper@37906
  2041
    then (*after Take: transfer Frm and respective istate*)
neuper@37906
  2042
	let val (ist_form, f) = (get_loc pt (p,Frm), 
neuper@37906
  2043
				 get_obj g_form pt p)
neuper@37906
  2044
	    val (pt, cs) = cut_tree pt (p,Frm)
neuper@37906
  2045
	    val pt = append_atomic p e_istate f r f' s pt
neuper@37925
  2046
	    val pt = update_loc' pt p (SOME ist_form, SOME ist_res)
neuper@37906
  2047
	in (pt, cs) end
neuper@37906
  2048
    else apfst (append_atomic p ist_res f r f' s) (cut_tree pt (p,Frm));
neuper@37906
  2049
neuper@37906
  2050
neuper@37906
  2051
(* called by Take *)
neuper@37906
  2052
fun append_form p l f pt = 
neuper@38015
  2053
((*tracing("##@append_form: pos ="^pos2str p);*)
neuper@37925
  2054
  insert (PrfObj {cell = NONE,
neuper@37906
  2055
		  form  = (*if existpt p pt 
neuper@37906
  2056
		  andalso get_obj g_tac pt p = Empty_Tac 
neuper@37906
  2057
			    (*distinction from 'old' (+complete!) pobjs*)
neuper@37906
  2058
			    then get_obj g_form pt p else*) f,
neuper@37906
  2059
		  tac  = Empty_Tac,
neuper@37925
  2060
		  loc   = (SOME l, NONE),
neuper@37906
  2061
		  branch= NoBranch,
neuper@37906
  2062
		  result= (e_term,[]),
neuper@37906
  2063
		  ostate= Incomplete}) pt p
neuper@37906
  2064
);
neuper@37906
  2065
(* val (p,loc,f) = ([1], e_istate, str2term "x + 1 = 2");
neuper@37906
  2066
   val (p,loc,f) = (fst p, e_istate, str2term "-1 + x = 0");
neuper@37906
  2067
   *)
neuper@37906
  2068
fun cappend_form pt p loc f =
neuper@38015
  2069
((*tracing("##@cappend_form: pos ="^pos2str p);*)
neuper@37906
  2070
  apfst (append_form p loc f) (cut_tree pt (p,Frm))
neuper@37906
  2071
);
neuper@37906
  2072
fun cappend_form pt p loc f =
neuper@38015
  2073
let (*val _= tracing("##@cappend_form: pos ="^pos2str p)
neuper@38015
  2074
    val _= tracing("##@cappend_form before cut_tree: loc ="^istate2str loc)*)
neuper@37906
  2075
    val (pt', cs) = cut_tree pt (p,Frm)
neuper@37906
  2076
    val pt'' = append_form p loc f pt'
neuper@38015
  2077
    (*val _= tracing("##@cappend_form after append: loc ="^
neuper@37906
  2078
		   istates2str (get_obj g_loc pt'' p))*)
neuper@37906
  2079
in (pt'', cs) end;
neuper@37906
  2080
neuper@37906
  2081
neuper@37906
  2082
    
neuper@37906
  2083
fun append_result pt p l f s =
neuper@38015
  2084
((*tracing("##@append_result: pos ="^pos2str p);*)
neuper@37906
  2085
    (appl_obj (repl_result (fst (get_obj g_loc pt p),
neuper@37925
  2086
			    SOME l) f s) pt p, [])
neuper@37906
  2087
);
neuper@37906
  2088
neuper@37906
  2089
neuper@37906
  2090
(*WN041022 deprecated, still for kbtest/diffapp.sml, /systest/root-equ.sml*)
neuper@37906
  2091
fun append_parent p l f r b pt = 
neuper@38015
  2092
  let (*val _= tracing("###append_parent: pos ="^pos2str p);*)
neuper@37906
  2093
    val (ll,f) = if existpt p pt andalso get_obj g_tac pt p=Empty_Tac
neuper@37925
  2094
		  then ((fst (get_obj g_loc pt p), SOME l), 
neuper@37906
  2095
			get_obj g_form pt p) 
neuper@37925
  2096
		 else ((SOME l, NONE), f)
neuper@37906
  2097
  in insert (PrfObj 
neuper@37925
  2098
	  {cell = NONE,
neuper@37906
  2099
	   form  = f,
neuper@37906
  2100
	   tac  = r,
neuper@37906
  2101
	   loc   = ll,
neuper@37906
  2102
	   branch= b,
neuper@37906
  2103
	   result= (e_term,[]),
neuper@37906
  2104
	   ostate= Incomplete}) pt p end;
neuper@37906
  2105
fun cappend_parent pt p loc f r b =
neuper@38015
  2106
((*tracing("###cappend_parent: pos ="^pos2str p);*)
neuper@37906
  2107
  apfst (append_parent p loc f r b) (cut_tree pt (p,Und))
neuper@37906
  2108
);
neuper@37906
  2109
neuper@37906
  2110
neuper@37906
  2111
fun append_problem [] l fmz (strs,spec,hdf) _ =
neuper@38015
  2112
((*tracing("###append_problem: pos = []");*)
neuper@37906
  2113
  (Nd (PblObj 
neuper@37925
  2114
	       {cell  = NONE,
neuper@37906
  2115
		origin= (strs,spec,hdf),
neuper@37906
  2116
		fmz   = fmz,
neuper@37906
  2117
		spec  = empty_spec,
neuper@37906
  2118
		probl = []:itm list,
neuper@37906
  2119
		meth  = []:itm list,
neuper@37925
  2120
		env   = NONE,
neuper@37925
  2121
		loc   = (SOME l, NONE),
neuper@37906
  2122
		branch= TransitiveB,(*FIXXXXXME.27.8.03: for equations only*)
neuper@37906
  2123
		result= (e_term,[]),
neuper@37906
  2124
		ostate= Incomplete},[]))
neuper@37906
  2125
)
neuper@37906
  2126
  | append_problem p l fmz (strs,spec,hdf) pt =
neuper@38015
  2127
((*tracing("###append_problem: pos ="^pos2str p);*)
neuper@37906
  2128
  insert (PblObj 
neuper@37925
  2129
	  {cell  = NONE,
neuper@37906
  2130
	   origin= (strs,spec,hdf),
neuper@37906
  2131
	   fmz   = fmz,
neuper@37906
  2132
	   spec  = empty_spec,
neuper@37906
  2133
	   probl = []:itm list,
neuper@37906
  2134
	   meth  = []:itm list,
neuper@37925
  2135
	   env   = NONE,
neuper@37925
  2136
	   loc   = (SOME l, NONE),
neuper@37906
  2137
	   branch= TransitiveB,
neuper@37906
  2138
	   result= (e_term,[]),
neuper@37906
  2139
	   ostate= Incomplete}) pt p
neuper@37906
  2140
);
neuper@37906
  2141
fun cappend_problem _ [] loc fmz ori =
neuper@38015
  2142
((*tracing("###cappend_problem: pos = []");*)
neuper@37906
  2143
  (append_problem [] loc fmz ori EmptyPtree,[])
neuper@37906
  2144
)
neuper@37906
  2145
  | cappend_problem pt p loc fmz ori = 
neuper@38015
  2146
((*tracing("###cappend_problem: pos ="^pos2str p);*)
neuper@37906
  2147
  apfst (append_problem p (loc:istate) fmz ori) (cut_tree pt (p,Frm))
neuper@37906
  2148
);
neuper@37906
  2149
neuper@37906
  2150
(*.get the theory explicitly specified for the rootpbl;
neuper@37906
  2151
   thus use this function _after_ finishing specification.*)
neuper@37906
  2152
fun rootthy (Nd (PblObj {spec=(thyID, _, _),...}, _)) = assoc_thy thyID
neuper@38031
  2153
  | rootthy _ = error "rootthy";
neuper@37906
  2154