src/Tools/isac/Interpret/calchead.sml
author Walther Neuper <neuper@ist.tugraz.at>
Thu, 05 May 2011 14:21:54 +0200
branchdecompose-isar
changeset 41976 792c59bf54d4
parent 41975 61f358925792
child 41977 a3ce4017f41d
permissions -rw-r--r--
intermed. ctxt over all minisubpbl x+1=2

overview in test:
--- all handling ctxt in minimsubpbl x+1=2 ---
ctxt initialised in prep_ori; the latter is
* not handled correctly in input_icalhd, match_pbl, refin'
* still unclear for subproblems in script
     1 (* Title:  Specify-phase: specifying and modeling a problem or a subproblem. The
     2            most important types are declared in mstools.sml.
     3    Author: Walther Neuper 991122, Mathias Lehnfeld
     4    (c) due to copyright terms
     5 
     6 12345678901234567890123456789012345678901234567890123456789012345678901234567890
     7         10        20        30        40        50        60        70        80
     8 *)
     9 
    10 (* TODO interne Funktionen aus sig entfernen *)
    11 signature CALC_HEAD =
    12   sig
    13     datatype additm = Add of SpecifyTools.itm | Err of string
    14     val all_dsc_in : SpecifyTools.itm_ list -> Term.term list
    15     val all_modspec : ptree * pos' -> ptree * pos'
    16     datatype appl = Appl of tac_ | Notappl of string
    17     val appl_add :
    18        Proof.context ->
    19        string ->
    20        SpecifyTools.ori list ->
    21        SpecifyTools.itm list ->
    22        (string * (Term.term * Term.term)) list -> cterm' -> additm
    23     type calcstate
    24     type calcstate'
    25     val chk_vars : term ppc -> string * Term.term list
    26     val chktyp :
    27        theory -> int * term list * term list -> term
    28     val chktyps :
    29        theory -> term list * term list -> term list
    30     val complete_metitms :
    31    SpecifyTools.ori list ->
    32    SpecifyTools.itm list ->
    33    SpecifyTools.itm list -> pat list -> SpecifyTools.itm list
    34     val complete_mod_ : ori list * pat list * pat list * itm list ->
    35 			itm list * itm list
    36     val complete_mod : ptree * pos' -> ptree * (pos * pos_)
    37     val complete_spec : ptree * pos' -> ptree * pos'
    38     val cpy_nam :
    39        pat list -> preori list -> pat -> preori
    40     val e_calcstate : calcstate
    41     val e_calcstate' : calcstate'
    42     val eq1 : ''a -> 'b * (''a * 'c) -> bool
    43     val eq3 :
    44        ''a -> Term.term -> 'b * 'c * 'd * ''a * SpecifyTools.itm_ -> bool
    45     val eq4 : ''a -> 'b * ''a list * 'c * 'd * 'e -> bool
    46     val eq5 :
    47        'a * 'b * 'c * 'd * SpecifyTools.itm_ ->
    48        'e * 'f * 'g * Term.term * 'h -> bool
    49     val eq_dsc : SpecifyTools.itm * SpecifyTools.itm -> bool
    50     val eq_pos' : ''a * pos_ -> ''a * pos_ -> bool
    51     val f_mout : theory -> mout -> Term.term
    52     val filter_outs :
    53        SpecifyTools.ori list ->
    54        SpecifyTools.itm list -> SpecifyTools.ori list
    55     val filter_pbt :
    56        SpecifyTools.ori list ->
    57        ('a * (Term.term * 'b)) list -> SpecifyTools.ori list
    58     val foldl1 : ('a * 'a -> 'a) -> 'a list -> 'a
    59     val foldr1 : ('a * 'a -> 'a) -> 'a list -> 'a
    60     val form : 'a -> ptree -> (string * ('a * pos_) * Term.term) list
    61     val formres : 'a -> ptree -> (string * ('a * pos_) * Term.term) list
    62     val gen_ins' : ('a * 'a -> bool) -> 'a * 'a list -> 'a list
    63     val get_formress :
    64        (string * (pos * pos_) * Term.term) list list ->
    65        pos -> ptree list -> (string * (pos * pos_) * Term.term) list
    66     val get_forms :
    67        (string * (pos * pos_) * Term.term) list list ->
    68        posel list -> ptree list -> (string * (pos * pos_) * Term.term) list
    69     val get_interval : pos' -> pos' -> int -> ptree -> (pos' * term) list
    70     val get_ocalhd : ptree * pos' -> ocalhd
    71     val get_spec_form : tac_ -> pos' -> ptree -> mout
    72     val geti_ct :
    73        theory ->
    74        SpecifyTools.ori -> SpecifyTools.itm -> string * cterm'
    75     val getr_ct : theory -> SpecifyTools.ori -> string * cterm'
    76     val has_list_type : Term.term -> bool
    77     val header : pos_ -> pblID -> metID -> pblmet
    78     val insert_ppc :
    79        theory ->
    80        int * SpecifyTools.vats * bool * string * SpecifyTools.itm_ ->
    81        SpecifyTools.itm list -> SpecifyTools.itm list
    82     val insert_ppc' :
    83        SpecifyTools.itm -> SpecifyTools.itm list -> SpecifyTools.itm list
    84     val is_complete_mod : ptree * pos' -> bool
    85     val is_complete_mod_ : SpecifyTools.itm list -> bool
    86     val is_complete_modspec : ptree * pos' -> bool
    87     val is_complete_spec : ptree * pos' -> bool
    88     val is_copy_named : 'a * ('b * Term.term) -> bool
    89     val is_copy_named_idstr : string -> bool
    90     val is_error : SpecifyTools.itm_ -> bool
    91     val is_field_correct : ''a -> ''b -> (''a * ''b list) list -> bool
    92     val is_known :
    93        Proof.context ->
    94        string ->
    95        SpecifyTools.ori list ->
    96        Term.term -> string * SpecifyTools.ori * Term.term list
    97     val is_list_type : Term.typ -> bool
    98     val is_notyet_input :
    99        Proof.context ->
   100        SpecifyTools.itm list ->
   101        Term.term list ->
   102        SpecifyTools.ori ->
   103        ('a * (Term.term * Term.term)) list -> string * SpecifyTools.itm
   104     val is_parsed : SpecifyTools.itm_ -> bool
   105     val is_untouched : SpecifyTools.itm -> bool
   106     val matc :
   107        theory ->
   108        pat list ->
   109        Term.term list ->
   110        (int list * string * Term.term * Term.term list) list ->
   111        (int list * string * Term.term * Term.term list) list
   112     val match_ags :
   113        theory -> pat list -> Term.term list -> SpecifyTools.ori list
   114     val maxl : int list -> int
   115     val match_ags_msg : string list -> Term.term -> Term.term list -> unit
   116     val memI : ''a list -> ''a -> bool
   117     val mk_additem : string -> cterm' -> tac
   118     val mk_delete : theory -> string -> SpecifyTools.itm_ -> tac
   119     val mtc :
   120        theory -> pat -> Term.term -> SpecifyTools.preori option
   121     val nxt_add :
   122        theory ->
   123        SpecifyTools.ori list ->
   124        (string * (Term.term * 'a)) list ->
   125        SpecifyTools.itm list -> (string * cterm') option
   126     val nxt_model_pbl : tac_ -> ptree * (int list * pos_) -> tac_
   127     val nxt_spec :
   128        pos_ ->
   129        bool ->
   130        SpecifyTools.ori list ->
   131        spec ->
   132        SpecifyTools.itm list * SpecifyTools.itm list ->
   133        (string * (Term.term * 'a)) list * (string * (Term.term * 'b)) list ->
   134        spec -> pos_ * tac
   135     val nxt_specif : tac -> ptree * (int list * pos_) -> calcstate'
   136     val nxt_specif_additem :
   137        string -> cterm' -> ptree * (int list * pos_) -> calcstate'
   138     val nxt_specify_init_calc : fmz -> calcstate
   139     val ocalhd_complete :
   140        SpecifyTools.itm list ->
   141        (bool * Term.term) list -> domID * pblID * metID -> bool
   142     val ori2Coritm :
   143 	pat list -> ori -> itm
   144     val ori_2itm :
   145        SpecifyTools.itm_ ->
   146        Term.term -> Term.term list -> SpecifyTools.ori -> SpecifyTools.itm
   147     val overwrite_ppc :
   148        theory ->
   149        int * SpecifyTools.vats * bool * string * SpecifyTools.itm_ ->
   150        SpecifyTools.itm list ->
   151        (int * SpecifyTools.vats * bool * string * SpecifyTools.itm_) list
   152     val parse_ok : SpecifyTools.itm_ list -> bool
   153     val posform2str : pos' * ptform -> string
   154     val posforms2str : (pos' * ptform) list -> string
   155     val posterms2str : (pos' * term) list -> string (*tests only*)
   156     val ppc135list : 'a SpecifyTools.ppc -> 'a list
   157     val ppc2list : 'a SpecifyTools.ppc -> 'a list
   158     val pt_extract :
   159        ptree * (int list * pos_) ->
   160        ptform * tac option * Term.term list
   161     val pt_form : ppobj -> ptform
   162     val pt_model : ppobj -> pos_ -> ptform
   163     val reset_calchead : ptree * pos' -> ptree * pos'
   164     val seek_oridts :
   165        Proof.context ->
   166        string ->
   167        Term.term * Term.term list ->
   168        (int * SpecifyTools.vats * string * Term.term * Term.term list) list
   169        -> string * SpecifyTools.ori * Term.term list
   170     val seek_orits :
   171        Proof.context ->
   172        string ->
   173        Term.term list ->
   174        (int * SpecifyTools.vats * string * Term.term * Term.term list) list
   175        -> string * SpecifyTools.ori * Term.term list
   176     val seek_ppc :
   177        int -> SpecifyTools.itm list -> SpecifyTools.itm option
   178     val show_pt : ptree -> unit
   179     val some_spec : spec -> spec -> spec
   180     val specify :
   181        tac_ ->
   182        pos' ->
   183        cid ->
   184        ptree ->
   185        (posel list * pos_) * ((posel list * pos_) * istate) * mout * tac *
   186        safe * ptree
   187     val specify_additem :
   188        string ->
   189        cterm' * 'a ->
   190        int list * pos_ ->
   191        'b ->
   192        ptree ->
   193        (pos * pos_) * ((pos * pos_) * istate) * mout * tac * safe * ptree
   194     val tag_form : theory -> term * term -> term
   195     val test_types : Proof.context -> Term.term * Term.term list -> string
   196     val typeless : Term.term -> Term.term
   197     val unbound_ppc : term SpecifyTools.ppc -> Term.term list
   198     val vals_of_oris : SpecifyTools.ori list -> Term.term list
   199     val variants_in : Term.term list -> int
   200     val vars_of_pbl_ : ('a * ('b * Term.term)) list -> Term.term list
   201     val vars_of_pbl_' : ('a * ('b * Term.term)) list -> Term.term list
   202   end
   203  
   204 
   205 
   206 
   207 
   208 structure CalcHead (**): CALC_HEAD(**) =
   209 
   210 struct
   211 
   212 (* datatypes *)
   213 
   214 (*.the state wich is stored after each step of calculation; it contains
   215    the calc-state and a list of [tac,istate](="tacis") to be applied next.
   216    the last_elem tacis is the first to apply to the calc-state and
   217    the (only) one shown to the front-end as the 'proposed tac'.
   218    the calc-state resulting from the application of tacis is not stored,
   219    because the tacis hold enough information for efficiently rebuilding
   220    this state just by "fun generate ".*)
   221 type calcstate = 
   222      (ptree * pos') *    (*the calc-state to which the tacis could be applied*)
   223      (taci list);        (*ev. several (hidden) steps; 
   224                            in REVERSE order: first tac_ to apply is last_elem*)
   225 val e_calcstate = ((EmptyPtree, e_pos'), [e_taci]):calcstate;
   226 
   227 (*the state used during one calculation within the mathengine; it contains
   228   a list of [tac,istate](="tacis") which generated the the calc-state;
   229   while this state's tacis are extended by each (internal) step,
   230   the calc-state is used for creating new nodes in the calc-tree
   231   (eg. applicable_in requires several particular nodes of the calc-tree)
   232   and then replaced by the the newly created;
   233   on leave of the mathengine the resuing calc-state is dropped anyway,
   234   because the tacis hold enought information for efficiently rebuilding
   235   this state just by "fun generate ".*)
   236 type calcstate' = 
   237      taci list *        (*cas. several (hidden) steps; 
   238                           in REVERSE order: first tac_ to apply is last_elem*)
   239      pos' list *        (*a "continuous" sequence of pos',
   240 			 deleted by application of taci list*)     
   241      (ptree * pos');    (*the calc-state resulting from the application of tacis*)
   242 val e_calcstate' = ([e_taci], [e_pos'], (EmptyPtree, e_pos')):calcstate';
   243 
   244 (*FIXXXME.WN020430 intermediate hack for fun ass_up*)
   245 fun f_mout thy (Form' (FormKF (_,_,_,_,f))) = (term_of o the o (parse thy)) f
   246   | f_mout thy _ = error "f_mout: not called with formula";
   247 
   248 
   249 (*.is the calchead complete ?.*)
   250 fun ocalhd_complete (its: itm list) (pre: (bool * term) list) (dI,pI,mI) = 
   251     foldl and_ (true, map #3 its) andalso 
   252     foldl and_ (true, map #1 pre) andalso 
   253     dI<>e_domID andalso pI<>e_pblID andalso mI<>e_metID;
   254 
   255 
   256 (* make a term 'typeless' for comparing with another 'typeless' term;
   257    'type-less' usually is illtyped                                  *)
   258 fun typeless (Const(s,_)) = (Const(s,e_type)) 
   259   | typeless (Free(s,_)) = (Free(s,e_type))
   260   | typeless (Var(n,_)) = (Var(n,e_type))
   261   | typeless (Bound i) = (Bound i)
   262   | typeless (Abs(s,_,t)) = Abs(s,e_type, typeless t)
   263   | typeless (t1 $ t2) = (typeless t1) $ (typeless t2);
   264 (*
   265 > val (SOME ct) = parse thy "max_relation (A=#2*a*b - a^^^#2)";
   266 > val (_,t1) = split_dsc_t hs (term_of ct);
   267 > val (SOME ct) = parse thy "A=#2*a*b - a^^^#2";
   268 > val (_,t2) = split_dsc_t hs (term_of ct);
   269 > typeless t1 = typeless t2;
   270 val it = true : bool
   271 *)
   272 
   273 
   274 
   275 (*.to an input (d,ts) find the according ori and insert the ts.*)
   276 (*WN.11.03: + dont take first inter<>[]*)
   277 fun seek_oridts ctxt sel (d,ts) [] =
   278     ("seek_oridts: input ('" ^
   279         (Print_Mode.setmp [] (Syntax.string_of_term ctxt) (comp_dts (d,ts))) ^
   280         "') not found in oris (typed)",
   281       (0, [], sel, d, ts),
   282       [])
   283   | seek_oridts ctxt sel (d,ts) ((id, vat, sel', d', ts')::oris) =
   284     if sel = sel' andalso d = d' andalso (inter op = ts ts') <> [] then
   285       ("", (id, vat, sel, d, inter op = ts ts'), ts') else
   286       seek_oridts ctxt sel (d, ts) oris;
   287 
   288 (*.to an input (_,ts) find the according ori and insert the ts.*)
   289 fun seek_orits ctxt sel ts [] = 
   290   ("seek_orits: input (_, '" ^ strs2str (map (Print_Mode.setmp [](Syntax.string_of_term 
   291                                            ctxt)) ts) ^
   292    "') not found in oris (typed)", e_ori_, [])
   293   | seek_orits ctxt sel ts ((id,vat,sel',d,ts')::(oris:ori list)) =
   294     if sel = sel' andalso (inter op = ts ts') <> [] 
   295     then if sel = sel' 
   296 	 then ("",
   297                (id,vat,sel,d, inter op = ts ts'):ori, 
   298                ts')
   299 	 else (((strs2str' o map (Print_Mode.setmp [] (Syntax.string_of_term
   300                                                            ctxt))) ts)
   301                ^ " not for " ^ sel, 
   302                e_ori_, 
   303                [])
   304     else seek_orits ctxt sel ts oris;
   305 (* false
   306 > val ((id,vat,sel',d,ts')::(ori':ori)) = ori;
   307 > seek_orits thy sel ts [(id,vat,sel',d,ts')];
   308 uncaught exception TYPE
   309 > seek_orits thy sel ts [];
   310 uncaught exception TYPE
   311 *)
   312 
   313 (*find_first item with #1 equal to id*)
   314 fun seek_ppc id [] = NONE
   315   | seek_ppc id (p::(ppc:itm list)) =
   316     if id = #1 p then SOME p else seek_ppc id ppc;
   317 
   318 
   319 
   320 (*---------------------------------------------(3) nach ptyps.sml 23.3.02*)
   321 
   322 
   323 datatype appl = Appl of tac_ | Notappl of string;
   324 
   325 fun ppc2list ({Given=gis,Where=whs,Find=fis,
   326 	       With=wis,Relate=res}: 'a ppc) =
   327   gis @ whs @ fis @ wis @ res;
   328 fun ppc135list ({Given=gis,Find=fis,Relate=res,...}: 'a ppc) =
   329   gis @ fis @ res;
   330 
   331 
   332 
   333 
   334 (* get the number of variants in a problem in 'original',
   335    assumes equal descriptions in immediate sequence    *)
   336 fun variants_in ts =
   337   let fun eq(x,y) = head_of x = head_of y;
   338     fun cnt eq [] y n = ([n],[])
   339       | cnt eq (x::xs) y n = if eq(x,y) then cnt eq xs y (n+1)
   340 			     else ([n], x::xs);
   341     fun coll eq  xs [] = xs
   342       | coll eq  xs (y::ys) = 
   343       let val (n,ys') = cnt eq (y::ys) y 0;
   344       in if ys' = [] then xs @ n else coll eq  (xs @ n) ys' end;
   345     val vts = subtract op = [1] (distinct (coll eq [] ts));
   346   in case vts of [] => 1 | [n] => n
   347       | _ => error "different variants in formalization" end;
   348 (*
   349 > cnt (op=) [2,2,2,4,5,5,5,5,5] 2 0;
   350 val it = ([3],[4,5,5,5,5,5]) : int list * int list
   351 > coll (op=) [] [1,2,2,2,4,5,5,5,5,5];
   352 val it = [1,3,1,5] : int list
   353 *)
   354 
   355 fun is_list_type (Type("List.list",_)) = true
   356   | is_list_type _ = false;
   357 (* fun destr (Type(str,sort)) = (str,sort);
   358 > val (SOME ct) = parse thy "lll::real list";
   359 > val ty = (#T o rep_cterm) ct;
   360 > is_list_type ty;
   361 val it = true : bool 
   362 > destr ty;
   363 val it = ("List.list",["RealDef.real"]) : string * typ list
   364 > atomty ((#t o rep_cterm) ct);
   365 *** -------------
   366 *** Free ( lll, real list)
   367 val it = () : unit
   368  
   369 > val (SOME ct) = parse thy "[lll::real]";
   370 > val ty = (#T o rep_cterm) ct;
   371 > is_list_type ty;
   372 val it = true : bool 
   373 > destr ty;
   374 val it = ("List.list",["'a"]) : string * typ list
   375 > atomty ((#t o rep_cterm) ct);
   376 *** -------------
   377 *** Const ( List.list.Cons, [real, real list] => real list)
   378 ***   Free ( lll, real)
   379 ***   Const ( List.list.Nil, real list) 
   380 
   381 > val (SOME ct) = parse thy "lll";
   382 > val ty = (#T o rep_cterm) ct;
   383 > is_list_type ty;
   384 val it = false : bool  *)
   385 
   386 
   387 fun has_list_type (Free(_,T)) = is_list_type T
   388   | has_list_type _ = false;
   389 (*
   390 > val (SOME ct) = parse thy "lll::real list";
   391 > has_list_type (term_of ct);
   392 val it = true : bool
   393 > val (SOME ct) = parse thy "[lll::real]";
   394 > has_list_type (term_of ct);
   395 val it = false : bool *)
   396 
   397 fun is_parsed (Syn _) = false
   398   | is_parsed _ = true;
   399 fun parse_ok its = foldl and_ (true, map is_parsed its);
   400 
   401 fun all_dsc_in itm_s =
   402   let    
   403     fun d_in (Cor ((d,_),_)) = [d]
   404       | d_in (Syn c) = []
   405       | d_in (Typ c) = []
   406       | d_in (Inc ((d,_),_)) = [d]
   407       | d_in (Sup (d,_)) = [d]
   408       | d_in (Mis (d,_)) = [d];
   409   in (flat o (map d_in)) itm_s end;  
   410 
   411 (* 30.1.00 ---
   412 fun is_Syn (Syn _) = true
   413   | is_Syn (Typ _) = true
   414   | is_Syn _ = false;
   415  --- *)
   416 fun is_error (Cor (_,ts)) = false
   417   | is_error (Sup (_,ts)) = false
   418   | is_error (Inc (_,ts)) = false
   419   | is_error (Mis (_,ts)) = false
   420   | is_error _ = true;
   421 
   422 (* 30.1.00 ---
   423 fun ct_in (Syn (c)) = c
   424   | ct_in (Typ (c)) = c
   425   | ct_in _ = error "ct_in called for Cor .. Sup";
   426  --- *)
   427 
   428 (*#############################################################*)
   429 (*#############################################################*)
   430 (* vvv--- aus nnewcode.sml am 30.1.00 ---vvv *)
   431 
   432 
   433 (* testdaten besorgen:
   434    use"test-coil-kernel.sml";
   435    val (PblObj{origin=(oris,_,_),meth={ppc=itms,...},...}) = 
   436         get_obj I pt p;
   437   *)
   438 
   439 (* given oris, ppc, 
   440    variant V: oris union ppc => int, id ID: oris union ppc => int
   441 
   442    ppc is_complete == 
   443      EX vt:V. ALL r:oris --> EX i:ppc. ID r = ID i  &  complete i
   444 
   445    and
   446      @vt = max sum(i : ppc) V i
   447 *)
   448 
   449 
   450 
   451 (*
   452 > ((vts_cnt (vts_in itms))) itms;
   453 
   454 
   455 
   456 ---^^--test 10.3.
   457 > val vts = vts_in itms;
   458 val vts = [1,2,3] : int list
   459 > val nvts = vts_cnt vts itms;
   460 val nvts = [(1,6),(2,5),(3,7)] : (int * int) list
   461 > val mx = max2 nvts;
   462 val mx = (3,7) : int * int
   463 > val v = max_vt itms;
   464 val v = 3 : int
   465 --------------------------
   466 > 
   467 *)
   468 
   469 (*.get the first term in ts from ori.*)
   470 (* val (_,_,fd,d,ts) = hd miss;
   471    *)
   472 fun getr_ct thy ((_, _, fd, d, ts) : ori) =
   473     (fd, (Print_Mode.setmp [] (Syntax.string_of_term (thy2ctxt thy)) o 
   474           comp_dts) (d,[hd ts]) : cterm');
   475 (* val t = comp_dts (d,[hd ts]);
   476    *)
   477 
   478 (* get a term from ori, notyet input in itm.
   479    the term from ori is thrown back to a string in order to reuse
   480    machinery for immediate input by the user. *)
   481 fun geti_ct thy ((_, _, _, d, ts) : ori) ((_, _, _, fd, itm_) : itm) =  
   482     (fd, Print_Mode.setmp [] ((Syntax.string_of_term (thy2ctxt thy)) o comp_dts) 
   483                           (d, subtract op = (ts_in itm_) ts) : cterm');
   484 fun geti_ct thy ((_, _, _, d, ts) : ori) ((_, _, _, fd, itm_) : itm) =  
   485     (fd, Print_Mode.setmp [] ((Syntax.string_of_term (thy2ctxt thy)) o comp_dts) 
   486                           (d, subtract op = (ts_in itm_) ts) : cterm');
   487 
   488 (* in FE dsc, not dat: this is in itms ...*)
   489 fun is_untouched ((_,_,false,_,Inc((_,[]),_)):itm) = true
   490   | is_untouched _ = false;
   491 
   492 
   493 (* select an item in oris, notyet input in itms 
   494    (precondition: in itms are only Cor, Sup, Inc) *)
   495 local infix mem;
   496 fun x mem [] = false
   497   | x mem (y :: ys) = x = y orelse x mem ys;
   498 in 
   499 fun nxt_add thy ([]:ori list) pbt itms = (*root (only) ori...fmz=[]*)
   500   let
   501     fun test_d d ((i,_,_,_,itm_):itm) = (d = (d_in itm_)) andalso i<>0; 
   502     fun is_elem itms (f,(d,t)) = 
   503       case find_first (test_d d) itms of 
   504 	SOME _ => true | NONE => false;
   505   in case filter_out (is_elem itms) pbt of
   506 (* val ((f,(d,_))::itms) = filter_out (is_elem itms) pbt;
   507    *)
   508     (f,(d,_))::itms => 
   509     SOME (f : string, 
   510           ((Print_Mode.setmp [] (Syntax.string_of_term (thy2ctxt thy))) o
   511            comp_dts) (d, []) : cterm')
   512   | _ => NONE end
   513 
   514 (* val (thy,itms) = (assoc_thy (if dI=e_domID then dI' else dI),pbl);
   515    *)
   516   | nxt_add thy oris pbt itms =
   517   let
   518     fun testr_vt v ori = (curry (op mem) v) (#2 (ori:ori))
   519       andalso (#3 ori) <>"#undef";
   520     fun testi_vt v itm = (curry (op mem) v) (#2 (itm:itm));
   521     fun test_id ids r = curry (op mem) (#1 (r:ori)) ids;
   522     fun test_subset (itm:itm) ((_,_,_,d,ts):ori) = 
   523 	(d_in (#5 itm)) = d andalso subset op = (ts_in (#5 itm), ts);
   524     fun false_and_not_Sup((i,v,false,f,Sup _):itm) = false
   525       | false_and_not_Sup (i,v,false,f, _) = true
   526       | false_and_not_Sup  _ = false;
   527 
   528     val v = if itms = [] then 1 else max_vt itms;
   529     val vors = if v = 0 then oris else filter (testr_vt v) oris;(*oris..vat*)
   530     val vits = if v = 0 then itms (*because of dsc without dat*)
   531 	       else filter (testi_vt v) itms;                   (*itms..vat*)
   532     val icl = filter false_and_not_Sup vits; (* incomplete *)
   533   in if icl = [] 
   534      then case filter_out (test_id (map #1 vits)) vors of
   535 	      [] => NONE
   536 	    (* val miss = filter_out (test_id (map #1 vits)) vors;
   537 	       *)
   538 	    | miss => SOME (getr_ct thy (hd miss))
   539      else
   540 	 case find_first (test_subset (hd icl)) vors of
   541 	     (* val SOME ori = find_first (test_subset (hd icl)) vors;
   542 	      *)
   543 	     NONE => error "nxt_add: EX itm. not(dat(itm)<=dat(ori))"
   544 	   | SOME ori => SOME (geti_ct thy ori (hd icl))
   545   end
   546 end;
   547 
   548 
   549 
   550 fun mk_delete thy "#Given"  itm_ = Del_Given   (itm_out thy itm_)
   551   | mk_delete thy "#Find"   itm_ = Del_Find    (itm_out thy itm_)
   552   | mk_delete thy "#Relate" itm_ = Del_Relation(itm_out thy itm_)
   553   | mk_delete thy str _ = 
   554   error ("mk_delete: called with field '"^str^"'");
   555 fun mk_additem "#Given" ct = Add_Given ct
   556   | mk_additem "#Find"  ct = Add_Find ct    
   557   | mk_additem "#Relate"ct = Add_Relation ct
   558   | mk_additem str _ = 
   559   error ("mk_additem: called with field '"^str^"'");
   560 
   561 
   562 (* determine the next step of specification;
   563    not done here: Refine_Tacitly (otherwise *** unknown method: (..., no_met))
   564    eg. in rootpbl 'no_met': 
   565 args:
   566   preok          predicates are _all_ ok (and problem matches completely)
   567   oris           immediately from formalization 
   568   (dI',pI',mI')  specification coming from author/parent-problem
   569   (pbl,          item lists specified by user
   570    met)          -"-, tacitly completed by copy_probl
   571   (dI,pI,mI)     specification explicitly done by the user
   572   (pbt, mpc)     problem type, guard of method
   573 *)
   574 fun nxt_spec Pbl preok (oris : ori list) ((dI', pI', mI') : spec)
   575   ((pbl : itm list), (met : itm list)) (pbt, mpc) ((dI, pI, mI) : spec) = 
   576     ((*tracing"### nxt_spec Pbl";*)
   577      if dI' = e_domID andalso dI = e_domID then (Pbl, Specify_Theory dI')
   578      else if pI' = e_pblID andalso pI = e_pblID then (Pbl, Specify_Problem pI')
   579      else case find_first (is_error o #5) (pbl:itm list) of
   580 	      SOME (_, _, _, fd, itm_) => 
   581 	      (Pbl, mk_delete 
   582 	                (assoc_thy (if dI = e_domID then dI' else dI)) fd itm_)
   583 	    | NONE => 
   584 	      ((*tracing"### nxt_spec is_error NONE";*)
   585 	       case nxt_add (assoc_thy (if dI = e_domID then dI' else dI)) 
   586 		            oris pbt pbl of
   587 	           SOME (fd,ct') => ((*tracing"### nxt_spec nxt_add SOME";*)
   588 				     (Pbl, mk_additem fd ct'))
   589 	         | NONE => (*pbl-items complete*)
   590 	           if not preok then (Pbl, Refine_Problem pI')
   591 	           else
   592 		       if dI = e_domID then (Pbl, Specify_Theory dI')
   593 		       else if pI = e_pblID then (Pbl, Specify_Problem pI')
   594 		       else if mI = e_metID then (Pbl, Specify_Method mI')
   595 		       else
   596 			   case find_first (is_error o #5) met of
   597 			       SOME (_,_,_,fd,itm_) => 
   598 			       (Met, mk_delete (assoc_thy dI) fd itm_)
   599 			     | NONE => 
   600 			       (case nxt_add (assoc_thy dI) oris mpc met of
   601 				    SOME (fd,ct') => (*30.8.01: pre?!?*)
   602 				    (Met, mk_additem fd ct')
   603 				  | NONE => 
   604 				    ((*Solv 3.4.00*)Met, Apply_Method mI))))
   605 
   606   | nxt_spec Met preok oris (dI',pI',mI') (pbl, met) (pbt,mpc) (dI,pI,mI) = 
   607   ((*tracing"### nxt_spec Met"; *)
   608    case find_first (is_error o #5) met of
   609      SOME (_,_,_,fd,itm_) => 
   610 	 (Met, mk_delete (assoc_thy (if dI=e_domID then dI' else dI)) fd itm_)
   611    | NONE => 
   612        case nxt_add (assoc_thy (if dI=e_domID then dI' else dI))oris mpc met of
   613 	 SOME (fd,ct') => (Met, mk_additem fd ct')
   614        | NONE => 
   615 	   ((*tracing"### nxt_spec Met: nxt_add NONE";*)
   616 	    if dI = e_domID then (Met, Specify_Theory dI')
   617 	    else if pI = e_pblID then (Met, Specify_Problem pI')
   618 		 else if not preok then (Met, Specify_Method mI)
   619 		      else (Met, Apply_Method mI)));
   620 	  
   621 (* di_ pI_ mI_ pos_
   622 val itms = [(1,[1],true,"#Find",Cor(e_term,[e_term])):itm,
   623 	    (2,[2],true,"#Find",Syn("empty"))];
   624 *)
   625 
   626 fun is_field_correct sel d dscpbt =
   627   case assoc (dscpbt, sel) of
   628     NONE => false
   629   | SOME ds => member op = ds d;
   630 
   631 (*. update the itm_ already input, all..from ori .*)
   632 (* val (id,vt,fd,d,ts) = (i,v,f,d,ts\\ts');
   633    *)
   634 fun ori_2itm itm_ pid all ((id,vt,fd,d,ts):ori) = 
   635   let 
   636     val ts' = union op = (ts_in itm_) ts;
   637     val pval = pbl_ids' d ts'
   638 	(*WN.9.5.03: FIXXXME [#0, epsilon]
   639 	  here would upd_penv be called for [#0, epsilon] etc. *)
   640     val complete = if eq_set op = (ts', all) then true else false;
   641   in case itm_ of
   642     (Cor _) => 
   643 	(if fd = "#undef" then (id,vt,complete,fd,Sup(d,ts')) 
   644 	 else (id,vt,complete,fd,Cor((d,ts'),(pid, pval)))):itm
   645   | (Syn c)     => error ("ori_2itm wants to overwrite "^c)
   646   | (Typ c)     => error ("ori_2itm wants to overwrite "^c)
   647   | (Inc _) => if complete
   648 	       then (id,vt,true ,fd, Cor ((d,ts'),(pid, pval)))
   649 	       else (id,vt,false,fd, Inc ((d,ts'),(pid, pval)))
   650   | (Sup ((*_,_*)d,ts')) => (*4.9.01 lost env*)
   651 	 (*if fd = "#undef" then*) (id,vt,complete,fd,Sup(d,ts'))
   652 	 (*else (id,vt,complete,fd,Cor((d,ts'),e))*)
   653 (* 28.1.00: not completely clear ---^^^ etc.*)
   654 (* 4.9.01: Mis just copied---vvv *)
   655   | (Mis _) => if complete
   656 		     then (id,vt,true ,fd, Cor ((d,ts'),(pid, pval)))
   657 		     else (id,vt,false,fd, Inc ((d,ts'),(pid, pval)))
   658   end;
   659 
   660 
   661 fun eq1 d (_,(d',_)) = (d = d');
   662 fun eq3 f d (_,_,_,f',itm_) = f = f' andalso d = (d_in itm_); 
   663 
   664 
   665 (* 'all' ts from ori; ts is the input; (ori carries rest of info)
   666    9.01: this + ori_2itm is _VERY UNCLEAR_ ? overhead ?
   667    pval: value for problem-environment _NOT_ checked for 'inter' --
   668    -- FIXXME.WN.11.03 the generation of penv has to go to insert_ppc
   669   (as it has been done for input_icalhd+insert_ppc' in 11.03)*)
   670 (*. is_input ori itms <=> 
   671     EX itm. (1) ori(field,dsc) = itm(field,dsc) & (2..4)
   672             (2) ori(ts) subset itm(ts)        --- Err "already input"       
   673 	    (3) ori(ts) inter itm(ts) = empty --- new: ori(ts)
   674 	    (4) -"- <> empty                  --- new: ori(ts) \\ inter .*)
   675 (* val(itms,(i,v,f,d,ts)) = (ppc,ori');
   676    *)
   677 fun is_notyet_input ctxt (itms:itm list) all ((i,v,f,d,ts):ori) pbt =
   678   case find_first (eq1 d) pbt of
   679       SOME (_,(_,pid)) =>(* val SOME (_,(_,pid)) = find_first (eq1 d) pbt;
   680                             val SOME (_,_,_,_,itm_)=find_first (eq3 f d) itms;
   681 			   *)
   682       (case find_first (eq3 f d) itms of
   683 	   SOME (_,_,_,_,itm_) =>
   684 	   let 
   685 	       val ts' = inter op = (ts_in itm_) ts;
   686 	   in if subset op = (ts, ts') 
   687 	      then (((strs2str' o 
   688 		      map (Print_Mode.setmp [] (Syntax.string_of_term ctxt))) ts') ^
   689 		    " already input", e_itm)                            (*2*)
   690 	      else ("", 
   691                     ori_2itm itm_ pid all (i,v,f,d,
   692                                                subtract op = ts' ts))   (*3,4*)
   693 	   end
   694 	 | NONE => ("", ori_2itm (Inc ((e_term,[]),(pid,[]))) 
   695 				 pid all (i,v,f,d,ts))                  (*1*)
   696 	)
   697     | NONE => ("", ori_2itm (Sup (d,ts)) 
   698 			      e_term all (i,v,f,d,ts));
   699 
   700 fun test_types ctxt (d,ts) =
   701   let 
   702     (*val s = !show_types; val _ = show_types:= true;*)
   703     val opt = (try comp_dts) (d,ts);
   704     val msg = case opt of 
   705       SOME _ => "" 
   706     | NONE => ((Print_Mode.setmp [] (Syntax.string_of_term ctxt) d) ^
   707                " " ^
   708 	       ((strs2str' o map (Print_Mode.setmp [] (Syntax.string_of_term
   709                                                            ctxt))) ts)
   710 	     ^ " is illtyped");
   711     (*val _ = show_types:= s*)
   712   in msg end;
   713 
   714 
   715 
   716 fun maxl [] = error "maxl of []"
   717   | maxl (y::ys) =
   718   let fun mx x [] = x
   719 	| mx x (y::ys) = if x < (y:int) then mx y ys else mx x ys
   720   in mx y ys end;
   721 
   722 
   723 (*. is the input term t known in oris ? 
   724     give feedback on all(?) strange input;
   725     return _all_ terms already input to this item (e.g. valuesFor a,b) .*)
   726 (*WN.11.03: from lists*)
   727 fun is_known ctxt sel ori t =
   728 (* val (ori,t)=(oris,term_of ct);
   729    *)
   730   let
   731     val _ = tracing ("RM is_known: t=" ^ term2str t);
   732     val ots = (distinct o flat o (map #5)) (ori:ori list);
   733     val oids = ((map (fst o dest_Free)) o distinct o 
   734 		flat o (map vars)) ots;
   735     val (d,ts(*,pval*)) = split_dts t;
   736     val ids = map (fst o dest_Free) 
   737       ((distinct o (flat o (map vars))) ts);
   738   in if (subtract op = oids ids) <> []
   739      then (("identifiers "^(strs2str' (subtract op = oids ids))^
   740 	    " not in example"), e_ori_, [])
   741      else 
   742 	 if d = e_term 
   743 	 then 
   744 	     if not (subset op = (map typeless ts, map typeless ots))
   745 	     then (("terms '"^
   746 		    ((strs2str' o 
   747                       (map (Print_Mode.setmp [] (Syntax.string_of_term 
   748 					             ctxt)))) ts)^
   749 		    "' not in example (typeless)"), e_ori_, [])
   750 	     else (case seek_orits ctxt sel ts ori of
   751 		       ("", ori_ as (_,_,_,d,ts), all) =>
   752 		       (case test_types ctxt (d,ts) of
   753 			    "" => ("", ori_, all)
   754 			  | msg => (msg, e_ori_, []))
   755 		     | (msg,_,_) => (msg, e_ori_, []))
   756 	 else 
   757 	     if member op = (map #4 ori) d
   758 	     then seek_oridts ctxt sel (d,ts) ori
   759 	     else ((Print_Mode.setmp [] (Syntax.string_of_term ctxt) d) ^
   760 		   " not in example", (0, [], sel, d, ts), [])
   761   end;
   762 
   763 
   764 (*. for return-value of appl_add .*)
   765 datatype additm =
   766 	 Add of itm
   767        | Err of string;    (*error-message*)
   768 
   769 
   770 (** add an item to the model; check wrt. oris and pbt **)
   771 (* in contrary to oris<>[] below, this part handles user-input
   772    extremely acceptive, i.e. accept input instead error-msg *)
   773 fun appl_add ctxt sel [] ppc pbt ct =
   774       let
   775         val i = 1 + (if ppc = [] then 0 else map #1 ppc |> maxl);
   776         in 
   777           case parseNEW ctxt ct of
   778             NONE => Add (i, [], false, sel, Syn ct)
   779           | SOME t =>
   780               let val (d, ts) = split_dts t;
   781               in 
   782                 if d = e_term 
   783                 then Add (i, [], false, sel, Mis (dsc_unknown, hd ts)) 
   784                 else
   785                   (case find_first (eq1 d) pbt of
   786                      NONE => Add (i, [], true, sel, Sup (d,ts))
   787                    | SOME (f, (_, id)) =>
   788                        let fun eq2 d (i, _, _, _, itm_) =
   789                              d = (d_in itm_) andalso i <> 0;
   790                        in case find_first (eq2 d) ppc of
   791                             NONE => Add (i, [], true, f,Cor ((d, ts), (id, pbl_ids' d ts)))
   792                           | SOME (i', _, _, _, itm_) => 
   793                               if is_list_dsc d 
   794                               then 
   795                                 let
   796                                   val in_itm = ts_in itm_;
   797                                   val ts' = union op = ts in_itm;
   798                                   val i'' = if in_itm = [] then i else i';
   799                                 in Add (i'', [], true, f, Cor ((d, ts'), (id, pbl_ids' d ts'))) end
   800                               else Add (i', [], true, f, Cor ((d, ts),(id, pbl_ids' d ts)))
   801                        end)
   802               end
   803       end
   804   | appl_add ctxt sel oris ppc pbt str =
   805       case parseNEW ctxt str of
   806         NONE => Err ("appl_add: syntax error in '" ^ str ^ "'")
   807       | SOME t => 
   808           case is_known ctxt sel oris t of
   809             ("", ori', all) => 
   810               (case is_notyet_input ctxt ppc all ori' pbt of
   811                  ("", itm) => Add itm
   812                | (msg, _) => Err ("[error] appl_add: is_notyet_input: " ^ msg))
   813           | (msg, _, _) => Err ("[error] appl_add: is_known: " ^ msg);
   814 
   815 
   816 (** make oris from args of the stac SubProblem and from pbt **)
   817 (* can this formal argument (of a model-pattern) be omitted in the arg-list
   818    of a SubProblem ? see ME/ptyps.sml 'type met ' *)
   819 fun is_copy_named_idstr str =
   820     case (rev o Symbol.explode) str of
   821 	"'":: _ ::"'"::_ => (tracing ((strs2str o (rev o Symbol.explode)) 
   822                                   "is_copy_named_idstr: " ^ str ^ " T"); true)
   823       | _ => (tracing ((strs2str o (rev o Symbol.explode)) 
   824                                   "is_copy_named_idstr: " ^ str ^ " F"); false);
   825 fun is_copy_named (_, (_, t)) = (is_copy_named_idstr o free2str) t;
   826 
   827 (*.should this formal argument (of a model-pattern) create a new identifier?.*)
   828 fun is_copy_named_generating_idstr str =
   829     if is_copy_named_idstr str
   830     then case (rev o Symbol.explode) str of
   831 	"'"::"'"::"'"::_ => false
   832       | _ => true
   833     else false;
   834 fun is_copy_named_generating (_, (_, t)) = 
   835     (is_copy_named_generating_idstr o free2str) t;
   836 
   837 (*.split type-wrapper from scr-arg and build part of an ori;
   838    an type-error is reported immediately, raises an exn, 
   839    subsequent handling of exn provides 2nd part of error message.*)
   840 fun mtc thy ((str, (dsc, _)):pat) (ty $ var) =
   841     (* val (thy, (str, (dsc, _)), (ty $ var)) =
   842 	   (thy,  p,               a);
   843        *)
   844     (cterm_of thy (dsc $ var);(*type check*)
   845      SOME ((([1], str, dsc, (*[var]*)
   846 	    split_dts' (dsc, var))): preori)(*:ori without leading #*))
   847     handle e as TYPE _ => 
   848 	   (tracing (dashs 70 ^ "\n"
   849 	^"*** ERROR while creating the items for the model of the ->problem\n"
   850 	^"*** from the ->stac with ->typeconstructor in arglist:\n"
   851 	^"*** item (->description ->value): "^term2str dsc^" "^term2str var^"\n"
   852 	^"*** description: "^(term_detail2str dsc)
   853 	^"*** value: "^(term_detail2str var)
   854 	^"*** typeconstructor in script: "^(term_detail2str ty)
   855 	^"*** checked by theory: "^(theory2str thy)^"\n"
   856 	^"*** " ^ dots 66);
   857 	    (*OldGoals.print_exn e; raises exn again*)
   858             writeln (PolyML.makestring e);
   859             reraise e;
   860 	(*raise ERROR "actual args do not match formal args";FIXXXME.WN100916*)
   861 	NONE);
   862 
   863 (*.match each pat of the model-pattern with an actual argument;
   864    precondition: copy-named vars are filtered out.*)
   865 fun matc thy ([]:pat list)  _  (oris:preori list) = oris
   866   | matc thy pbt [] _ =
   867     (tracing (dashs 70);
   868      error ("actual arg(s) missing for '"^pats2str pbt
   869 		 ^"' i.e. should be 'copy-named' by '*_._'"))
   870   | matc thy ((p as (s,(d,t)))::pbt) (a::ags) oris =
   871     (* val (thy, ((p as (s,(d,t)))::pbt), (a::ags), oris) =
   872 	   (thy,  pbt',                    ags,     []);
   873        (*recursion..*)
   874        val (thy, ((p as (s,(d,t)))::pbt), (a::ags), oris) =
   875 	   (thy,  pbt,                     ags,    (oris @ [ori]));
   876        *)
   877     (*del?..*)if (is_copy_named_idstr o free2str) t then oris
   878     else(*..del?*) let val opt = mtc thy p a;  
   879 	 in case opt of
   880 		(* val SOME ori = mtc thy p a;
   881 		   *)
   882 		SOME ori => matc thy pbt ags (oris @ [ori])
   883 	      | NONE => [](*WN050903 skipped by exn handled in match_ags*)
   884 	 end; 
   885 (* run subp-rooteq.sml until Init_Proof before ...
   886 > val Nd (PblObj {origin=(oris,_,_),...},_) = pt;(*from test/subp-rooteq.sml*)
   887 > fun xxxfortest (_,a,b,c,d) = (a,b,c,d);val oris = map xxxfortest oris;
   888 
   889  other vars as in mtc ..
   890 > matc thy (drop_last pbt) ags [];
   891 val it = ([[1],"#Given",Const #,[#]),(0,[#],"#Given",Const #,[#])],2)*)
   892 
   893 
   894 (* generate a new variable "x_i" name from a related given one "x"
   895    by use of oris relating "v_i_" (is_copy_named!) to "v_"
   896    e.g. (v_, x) & (v_i_, ?) --> (v_i_, x_i) *)
   897 
   898 (* generate a new variable "x_i" name from a related given one "x"
   899    by use of oris relating "v_v'i'" (is_copy_named!) to "v_v"
   900    e.g. (v_v, x) & (v_v'i', ?) --> (v_v'i', x_i),
   901    but leave is_copy_named_generating as is, e.t. ss''' *)
   902 fun cpy_nam (pbt:pat list) (oris:preori list) (p as (field,(dsc,t)):pat) =
   903   (if is_copy_named_generating p
   904    then (*WN051014 kept strange old code ...*)
   905        let fun sel (_,_,d,ts) = comp_ts (d, ts) 
   906 	   val cy' = (implode o (drop_last_n 3) o Symbol.explode o free2str) t
   907 	   val ext = (last_elem o drop_last o Symbol.explode o free2str) t
   908 	   val vars' = map (free2str o snd o snd) pbt (*cpy-nam filtered_out*)
   909 	   val vals = map sel oris
   910 	   val cy_ext = (free2str o the) (assoc (vars'~~vals, cy'))^"_"^ext
   911        in ([1], field, dsc, [mk_free (type_of t) cy_ext]):preori end
   912    else ([1], field, dsc, [t])
   913 	)
   914   handle _ => error ("cpy_nam: for "^(term2str t));
   915 
   916 (*.match the actual arguments of a SubProblem with a model-pattern
   917    and create an ori list (in root-pbl created from formalization).
   918    expects ags:pats = 1:1, while copy-named are filtered out of pats;
   919    if no 1:1 the exn raised by matc/mtc and handled at call.
   920    copy-named pats are appended in order to get them into the model-items.*)
   921 fun match_ags thy (pbt:pat list) ags =
   922     let fun flattup (i,(var,bool,str,itm_)) = (i,var,bool,str,itm_);
   923 	val pbt' = filter_out is_copy_named pbt;
   924 	val cy = filter is_copy_named pbt;
   925 	val oris' = matc thy pbt' ags [];
   926 	val cy' = map (cpy_nam pbt' oris') cy;
   927 	val ors = add_id (oris' @ cy'); 
   928     (*appended in order to get ^^^^^ into the model-items*)
   929     in (map flattup ors):ori list end;
   930 
   931 (*.report part of the error-msg which is not available in match_args.*)
   932 fun match_ags_msg pI stac ags =
   933     let (*val s = !show_types
   934 	val _ = show_types:= true*)
   935 	val pats = (#ppc o get_pbt) pI
   936 	val msg = (dots 70^"\n"
   937 		 ^"*** problem "^strs2str pI^" has the ...\n"
   938 		 ^"*** model-pattern "^pats2str pats^"\n"
   939 		 ^"*** stac   '"^term2str stac^"' has the ...\n"
   940 		 ^"*** arg-list "^terms2str ags^"\n"
   941 		 ^dashs 70)
   942 	(*WN100921 ^^^^ expect TYPE errormsg below; lost with Isa09*)
   943 	(*val _ = show_types:= s*)
   944     in tracing msg end;
   945 
   946 
   947 (*get the variables out of a pbl_; FIXME.WN.0311: is_copy_named ...obscure!!!*)
   948 fun vars_of_pbl_ pbl_ = 
   949     let fun var_of_pbl_ (gfr,(dsc,t)) = t
   950     in ((map var_of_pbl_) o (filter_out is_copy_named)) pbl_ end;
   951 fun vars_of_pbl_' pbl_ = 
   952     let fun var_of_pbl_ (gfr,(dsc,t)) = t:term
   953     in ((map var_of_pbl_)(* o (filter_out is_copy_named)*)) pbl_ end;
   954 
   955 fun overwrite_ppc thy itm ppc =
   956   let 
   957     fun repl ppc' (_,_,_,_,itm_) [] =
   958       error ("overwrite_ppc: " ^ (itm_2str_ (thy2ctxt thy) itm_) ^ 
   959                    " not found")
   960       | repl ppc' itm (p::ppc) =
   961 	if (#1 itm) = (#1 (p:itm)) then ppc' @ [itm] @ ppc
   962 	else repl (ppc' @ [p]) itm ppc
   963   in repl [] itm ppc end;
   964 
   965 (*10.3.00: insert the already compiled itm into model;
   966    ev. filter_out  untouched (in FE: (0,...)) item related to insert-item *)
   967 (* val ppc=pbl;
   968    *)
   969 fun insert_ppc thy itm ppc =
   970     let 
   971 	fun eq_untouched d ((0,_,_,_,itm_):itm) = (d = d_in itm_)
   972 	  | eq_untouched _ _ = false;
   973 	    val ppc' = 
   974 		(
   975 		 (*tracing("### insert_ppc: itm= "^(itm2str_ itm));*)       
   976 		 case seek_ppc (#1 itm) ppc of
   977 		     (* val SOME xxx = seek_ppc (#1 itm) ppc;
   978 		        *)
   979 		     SOME _ => (*itm updated in is_notyet_input WN.11.03*)
   980 		     overwrite_ppc thy itm ppc
   981 		   | NONE => (ppc @ [itm]));
   982     in filter_out (eq_untouched ((d_in o #5) itm)) ppc' end;
   983 
   984 (*from Isabelle/src/Pure/library.ML, _appends_ a new element*)
   985 fun gen_ins' eq (x, xs) = if gen_mem eq (x, xs) then xs else xs @ [x];
   986 
   987 fun eq_dsc ((_,_,_,_,itm_):itm, (_,_,_,_,iitm_):itm) = 
   988     (d_in itm_) = (d_in iitm_);
   989 (*insert_ppc = insert_ppc' for appl_add', input_icalhd 11.03,
   990     handles superfluous items carelessly*)
   991 fun insert_ppc' itm itms = gen_ins' eq_dsc (itm, itms);
   992 (* val eee = op=;
   993  > gen_ins' eee (4,[1,3,5,7]);
   994 val it = [1, 3, 5, 7, 4] : int list*)
   995 
   996 
   997 (*. output the headline to a ppc .*)
   998 fun header p_ pI mI =
   999     case p_ of Pbl => Problem (if pI = e_pblID then [] else pI) 
  1000 	     | Met => Method mI
  1001 	     | pos => error ("header called with "^ pos_2str pos);
  1002 
  1003 
  1004 fun specify_additem sel (ct,_) (p,Met) c pt = 
  1005     let
  1006       val (PblObj{meth=met,origin=(oris,(dI',pI',mI'),_),
  1007 		  probl=pbl,spec=(dI,pI,mI),...}) = get_obj I pt p;
  1008       val thy = if dI = e_domID then assoc_thy dI' else assoc_thy dI;
  1009     (*val ppt = if pI = e_pblID then get_pbt pI' else get_pbt pI;*)
  1010       val cpI = if pI = e_pblID then pI' else pI;
  1011       val cmI = if mI = e_metID then mI' else mI;
  1012       val {ppc,pre,prls,...} = get_met cmI;
  1013       val ctxt = get_ctxt pt (p,Met);
  1014     in case appl_add ctxt sel oris met ppc ct of
  1015       Add itm (*..union old input *) =>
  1016 	let (* val Add itm = appl_add thy sel oris met (#ppc (get_met cmI)) ct;
  1017                *)
  1018 	  val met' = insert_ppc thy itm met;
  1019 	  (*val pt' = update_met pt p met';*)
  1020 	  val ((p,Met),_,_,pt') = 
  1021 	      generate1 thy (case sel of
  1022 				 "#Given" => Add_Given' (ct, met')
  1023 			       | "#Find"  => Add_Find'  (ct, met')
  1024 			       | "#Relate"=> Add_Relation'(ct, met')) 
  1025 			(Uistate, e_ctxt) (p,Met) pt
  1026 	  val pre' = check_preconds thy prls pre met'
  1027 	  val pb = foldl and_ (true, map fst pre')
  1028 	  (*val _=tracing("@@@ specify_additem: Met Add before nxt_spec")*)
  1029 	  val (p_,nxt) =
  1030 	    nxt_spec Met pb oris (dI',pI',mI') (pbl,met') 
  1031 	    ((#ppc o get_pbt) cpI,ppc) (dI,pI,mI);
  1032 	in ((p,p_), ((p,p_),Uistate),
  1033 	    Form' (PpcKF (0,EdUndef,(length p),Nundef,
  1034 			  (Method cmI, itms2itemppc thy met' pre'))),
  1035 	    nxt,Safe,pt') end
  1036     | Err msg =>
  1037 	  let val pre' = check_preconds thy prls pre met
  1038 	      val pb = foldl and_ (true, map fst pre')
  1039 	    (*val _=tracing("@@@ specify_additem: Met Err before nxt_spec")*)
  1040 	      val (p_,nxt) =
  1041 	    nxt_spec Met pb oris (dI',pI',mI') (pbl,met) 
  1042 	    ((#ppc o get_pbt) cpI,(#ppc o get_met) cmI) (dI,pI,mI);
  1043 	  in ((p,p_), ((p,p_),Uistate), Error' (Error_ msg), nxt, Safe,pt) end
  1044     end
  1045 (* val (p,_) = p;
  1046    *)
  1047 | specify_additem sel (ct,_) (pos as (p,_(*Frm, Pbl*))) c pt = 
  1048     let
  1049       val (PblObj{meth=met,origin=(oris,(dI',pI',mI'),_),
  1050 		  probl=pbl,spec=(dI,pI,mI),...}) = get_obj I pt p;
  1051       val thy = if dI = e_domID then assoc_thy dI' else assoc_thy dI;
  1052       val cpI = if pI = e_pblID then pI' else pI;
  1053       val cmI = if mI = e_metID then mI' else mI;
  1054       val {ppc,where_,prls,...} = get_pbt cpI;
  1055       val ctxt = get_ctxt pt pos;
  1056     in case appl_add ctxt sel oris pbl ppc ct of
  1057       Add itm (*..union old input *) =>
  1058       (* val Add itm = appl_add thy sel oris pbl ppc ct;
  1059          *)
  1060 	let
  1061 	    (*val _= tracing("###specify_additem: itm= "^(itm2str_ itm));*)
  1062 	  val pbl' = insert_ppc thy itm pbl
  1063 	  val ((p,Pbl),_,_,pt') = 
  1064 	      generate1 thy (case sel of
  1065 				 "#Given" => Add_Given' (ct, pbl')
  1066 			       | "#Find"  => Add_Find'  (ct, pbl')
  1067 			       | "#Relate"=> Add_Relation'(ct, pbl')) 
  1068 			(Uistate, e_ctxt) (p,Pbl) pt
  1069 	  val pre = check_preconds thy prls where_ pbl'
  1070 	  val pb = foldl and_ (true, map fst pre)
  1071 	(*val _=tracing("@@@ specify_additem: Pbl Add before nxt_spec")*)
  1072 	  val (p_,nxt) =
  1073 	    nxt_spec Pbl pb oris (dI',pI',mI') (pbl',met) 
  1074 		     (ppc,(#ppc o get_met) cmI) (dI,pI,mI);
  1075 	  val ppc = if p_= Pbl then pbl' else met;
  1076 	in ((p,p_), ((p,p_),Uistate),
  1077 	    Form' (PpcKF (0,EdUndef,(length p),Nundef,
  1078 			  (header p_ pI cmI,
  1079 			   itms2itemppc thy ppc pre))), nxt,Safe,pt') end
  1080 
  1081     | Err msg =>
  1082 	  let val pre = check_preconds thy prls where_ pbl
  1083 	      val pb = foldl and_ (true, map fst pre)
  1084 	    (*val _=tracing("@@@ specify_additem: Pbl Err before nxt_spec")*)
  1085 	      val (p_,nxt) =
  1086 	    nxt_spec Pbl pb oris (dI',pI',mI') (pbl,met) 
  1087 	    (ppc,(#ppc o get_met) cmI) (dI,pI,mI);
  1088 	  in ((p,p_), ((p,p_),Uistate), Error' (Error_ msg), nxt, Safe,pt) end
  1089     end;
  1090 
  1091 fun specify (Init_Proof' (fmz,(dI',pI',mI'))) (_:pos') (_:cid) (_:ptree)= 
  1092   let          (* either """"""""""""""" all empty or complete *)
  1093     val thy = assoc_thy dI';
  1094     val (oris, ctxt) = 
  1095       if dI' = e_domID orelse pI' = e_pblID 
  1096       then ([], e_ctxt)
  1097 	    else pI' |> #ppc o get_pbt |> prep_ori fmz thy;
  1098     val (pt,c) = cappend_problem e_ptree [] (e_istate, ctxt) (fmz,(dI',pI',mI'))
  1099 			(oris, (dI',pI',mI'),e_term);
  1100     val {ppc, prls, where_, ...} = get_pbt pI'
  1101     val (pbl, pre, pb) = ([], [], false)
  1102   in 
  1103     case mI' of
  1104 	    ["no_met"] => 
  1105 	      (([], Pbl), (([], Pbl), Uistate),
  1106 	        Form' (PpcKF (0, EdUndef, (length []), Nundef,
  1107 			      (Problem [], itms2itemppc (assoc_thy dI') pbl pre))),
  1108 	        Refine_Tacitly pI', Safe, pt)
  1109      | _ => 
  1110 	      (([], Pbl), (([], Pbl), Uistate),
  1111 	        Form' (PpcKF (0, EdUndef, (length []), Nundef,
  1112 			      (Problem [], itms2itemppc (assoc_thy dI') pbl pre))),
  1113 	        Model_Problem, Safe, pt)
  1114   end
  1115 
  1116     (*ONLY for STARTING modeling phase*)
  1117   | specify (Model_Problem' (_,pbl,met)) (pos as (p,p_)) c pt =
  1118       let 
  1119         val (PblObj{origin=(oris,(dI',pI',mI'),_), spec=(dI,_,_),...}) = get_obj I pt p
  1120         val thy' = if dI = e_domID then dI' else dI
  1121         val thy = assoc_thy thy'
  1122         val {ppc,prls,where_,...} = get_pbt pI'
  1123         val pre = check_preconds thy prls where_ pbl
  1124         val pb = foldl and_ (true, map fst pre)
  1125         val ((p,_),_,_,pt) = generate1 thy (Model_Problem'([],pbl,met)) (Uistate, e_ctxt) pos pt
  1126         val (_,nxt) = nxt_spec Pbl pb oris (dI',pI',mI') (pbl,met) 
  1127 		      (ppc,(#ppc o get_met) mI') (dI',pI',mI');
  1128       in ((p, Pbl), ((p, p_), Uistate), 
  1129           Form' (PpcKF (0, EdUndef, length p, Nundef,
  1130 		       (Problem pI', itms2itemppc (assoc_thy dI') pbl pre))),
  1131           nxt, Safe, pt)
  1132       end
  1133 
  1134     (*. called only if no_met is specified .*)     
  1135   | specify (Refine_Tacitly' (pI,pIre,_,_,_)) (pos as (p,_)) c pt =
  1136   let (* val Refine_Tacitly' (pI,pIre,_,_,_) = m;
  1137          *)
  1138     val (PblObj{origin=(oris,(dI',pI',mI'),_), meth=met, ...}) = 
  1139 	get_obj I pt p;
  1140     val {prls,met,ppc,thy,where_,...} = get_pbt pIre
  1141     (*val pbl = init_pbl ppc --- Model_Problem recognizes probl=[]*)
  1142     (*val pt = update_pbl pt p pbl;
  1143     val pt = update_orispec pt p 
  1144 		(string_of_thy thy, pIre, 
  1145 		 if length met = 0 then e_metID else hd met);*)
  1146     val (domID, metID) = (string_of_thy thy, 
  1147 		      if length met = 0 then e_metID else hd met)
  1148     val ((p,_),_,_,pt) = 
  1149 	generate1 thy (Refine_Tacitly'(pI,pIre,domID,metID,(*pbl*)[])) 
  1150 		  (Uistate, e_ctxt) pos pt
  1151     (*val pre = check_preconds thy prls where_ pbl
  1152     val pb = foldl and_ (true, map fst pre)*)
  1153     val (pbl, pre, pb) = ([], [], false)
  1154   in ((p,Pbl), (pos,Uistate),
  1155       Form' (PpcKF (0,EdUndef,(length p),Nundef,
  1156 		    (Problem pIre, itms2itemppc (assoc_thy dI') pbl pre))),
  1157       Model_Problem, Safe, pt) end
  1158 
  1159   | specify (Refine_Problem' (rfd as (pI,_))) pos c pt =
  1160     let val (pos,_,_,pt) = generate1 (assoc_thy "Isac") 
  1161 				     (Refine_Problem' rfd) (Uistate, e_ctxt) pos pt
  1162     in (pos(*p,Pbl*), (pos(*p,Pbl*),Uistate), Problems (RefinedKF rfd), 
  1163 	Model_Problem, Safe, pt) end
  1164 
  1165 (* val (Specify_Problem' (pI, (ok, (itms, pre)))) = nxt; val (p,_) = p;
  1166    val (Specify_Problem' (pI, (ok, (itms, pre)))) = m; val (p,_) = p;
  1167    *)
  1168   | specify (Specify_Problem' (pI, (ok, (itms, pre)))) (pos as (p,_)) c pt =
  1169   let val (PblObj {origin=(oris,(dI',pI',mI'),_), spec=(dI,_,mI), 
  1170 		   meth=met, ...}) = get_obj I pt p;
  1171     (*val pt = update_pbl pt p itms;
  1172     val pt = update_pblID pt p pI;*)
  1173     val thy = assoc_thy dI
  1174     val ((p,Pbl),_,_,pt)= 
  1175 	generate1 thy (Specify_Problem' (pI, (ok, (itms, pre)))) (Uistate, e_ctxt) pos pt
  1176     val dI'' = assoc_thy (if dI=e_domID then dI' else dI);
  1177     val mI'' = if mI=e_metID then mI' else mI;
  1178   (*val _=tracing("@@@ specify (Specify_Problem) before nxt_spec")*)
  1179     val (_,nxt) = nxt_spec Pbl ok oris (dI',pI',mI') (itms, met) 
  1180 		((#ppc o get_pbt) pI,(#ppc o get_met) mI'') (dI,pI,mI);
  1181   in ((p,Pbl), (pos,Uistate),
  1182       Form' (PpcKF (0,EdUndef,(length p),Nundef,
  1183 		    (Problem pI, itms2itemppc dI'' itms pre))),
  1184       nxt, Safe, pt) end    
  1185 (* val Specify_Method' mID = nxt; val (p,_) = p;
  1186    val Specify_Method' mID = m;
  1187    specify (Specify_Method' mID) (p,p_) c pt;
  1188    *)
  1189   | specify (Specify_Method' (mID,_,_)) (pos as (p,_)) c pt =
  1190   let val (PblObj {origin=(oris,(dI',pI',mI'),_), probl=pbl, spec=(dI,pI,mI), 
  1191 		   meth=met, ...}) = get_obj I pt p;
  1192     val {ppc,pre,prls,...} = get_met mID
  1193     val thy = assoc_thy dI
  1194     val oris = add_field' thy ppc oris;
  1195     (*val pt = update_oris pt p oris; 20.3.02: repl. "#undef"*)
  1196     val dI'' = if dI=e_domID then dI' else dI;
  1197     val pI'' = if pI = e_pblID then pI' else pI;
  1198     val met = if met=[] then pbl else met;
  1199     val (ok, (itms, pre')) = match_itms_oris thy met (ppc,pre,prls ) oris;
  1200     (*val pt = update_met pt p itms;
  1201     val pt = update_metID pt p mID*)
  1202     val (pos,_,_,pt)= 
  1203 	generate1 thy (Specify_Method' (mID, oris, itms)) (Uistate, e_ctxt) pos pt
  1204     (*val _=tracing("@@@ specify (Specify_Method) before nxt_spec")*)
  1205     val (_,nxt) = nxt_spec Met (*ok*)true oris (dI',pI',mI') (pbl, itms) 
  1206 		((#ppc o get_pbt) pI'',ppc) (dI'',pI'',mID);
  1207   in (pos, (pos,Uistate),
  1208       Form' (PpcKF (0,EdUndef,(length p),Nundef,
  1209 		    (Method mID, itms2itemppc (assoc_thy dI'') itms pre'))),
  1210       nxt, Safe, pt) end    
  1211 (* val Add_Find' ct = nxt; val sel = "#Find"; 
  1212    *)
  1213   | specify (Add_Given' ct) p c pt = specify_additem "#Given" ct p c pt
  1214   | specify (Add_Find'  ct) p c pt = specify_additem "#Find"  ct p c pt
  1215   | specify (Add_Relation' ct) p c pt=specify_additem"#Relate"ct p c pt
  1216 (* val Specify_Theory' domID = m;
  1217    val (Specify_Theory' domID, (p,p_)) = (m, pos);
  1218    *)
  1219   | specify (Specify_Theory' domID) (pos as (p,p_)) c pt =
  1220     let val p_ = case p_ of Met => Met | _ => Pbl
  1221       val thy = assoc_thy domID;
  1222       val (PblObj{origin=(oris,(dI',pI',mI'),_), meth=met,
  1223 		  probl=pbl, spec=(dI,pI,mI),...}) = get_obj I pt p;
  1224       val mppc = case p_ of Met => met | _ => pbl;
  1225       val cpI = if pI = e_pblID then pI' else pI;
  1226       val {prls=per,ppc,where_=pwh,...} = get_pbt cpI
  1227       val cmI = if mI = e_metID then mI' else mI;
  1228       val {prls=mer,ppc=mpc,pre=mwh,...} = get_met cmI
  1229       val pre = 
  1230 	  case p_ of
  1231 	      Met => (check_preconds thy mer mwh met)
  1232 	    | _ => (check_preconds thy per pwh pbl)
  1233       val pb = foldl and_ (true, map fst pre)
  1234     in if domID = dI
  1235        then let 
  1236 	 (*val _=tracing("@@@ specify (Specify_Theory) THEN before nxt_spec")*)
  1237            val (p_,nxt) = nxt_spec p_ pb oris (dI',pI',mI') 
  1238 				   (pbl,met) (ppc,mpc) (dI,pI,mI);
  1239 	      in ((p,p_), (pos,Uistate), 
  1240 		  Form'(PpcKF (0,EdUndef,(length p), Nundef,
  1241 			       (header p_ pI cmI, itms2itemppc thy mppc pre))),
  1242 		  nxt,Safe,pt) end
  1243        else (*FIXME: check ppc wrt. (new!) domID ..? still parsable?*)
  1244 	 let 
  1245 	   (*val pt = update_domID pt p domID;11.8.03*)
  1246 	   val ((p,p_),_,_,pt) = generate1 thy (Specify_Theory' domID) 
  1247 					   (Uistate, e_ctxt) (p,p_) pt
  1248 	 (*val _=tracing("@@@ specify (Specify_Theory) ELSE before nxt_spec")*)
  1249 	   val (p_,nxt) = nxt_spec p_ pb oris (dI',pI',mI') (pbl,met) 
  1250 				   (ppc,mpc) (domID,pI,mI);
  1251 	 in ((p,p_), (pos,Uistate), 
  1252 	     Form' (PpcKF (0, EdUndef, (length p),Nundef,
  1253 			   (header p_ pI cmI, itms2itemppc thy mppc pre))),
  1254 	     nxt, Safe,pt) end
  1255     end
  1256 (* itms2itemppc thy [](*mpc*) pre
  1257    *)
  1258   | specify m' _ _ _ = 
  1259     error ("specify: not impl. for "^tac_2str m');
  1260 
  1261 (* val (sel, Add_Given ct, ptp as (pt,(p,Pbl))) = ("#Given", tac, ptp);
  1262    val (sel, Add_Find  ct, ptp as (pt,(p,Pbl))) = ("#Find", tac, ptp);
  1263    *)
  1264 fun nxt_specif_additem sel ct (ptp as (pt, (p, Pbl))) = 
  1265     let
  1266       val (PblObj{meth=met,origin=(oris,(dI',pI',_),_),
  1267 		  probl=pbl,spec=(dI,pI,_),...}) = get_obj I pt p;
  1268       val thy = if dI = e_domID then assoc_thy dI' else assoc_thy dI;
  1269       val cpI = if pI = e_pblID then pI' else pI;
  1270       val ctxt = get_ctxt pt (p, Pbl);
  1271     in case appl_add ctxt sel oris pbl ((#ppc o get_pbt) cpI) ct of
  1272 	   Add itm (*..union old input *) =>
  1273 (* val Add itm = appl_add thy sel oris pbl ppc ct;
  1274    *)
  1275 	   let
  1276 	       (*val _=tracing("###nxt_specif_additem: itm= "^(itm2str_ itm));*)
  1277 	       val pbl' = insert_ppc thy itm pbl
  1278 	       val (tac,tac_) = 
  1279 		   case sel of
  1280 		       "#Given" => (Add_Given    ct, Add_Given'   (ct, pbl'))
  1281 		     | "#Find"  => (Add_Find     ct, Add_Find'    (ct, pbl'))
  1282 		     | "#Relate"=> (Add_Relation ct, Add_Relation'(ct, pbl'))
  1283 	       val ((p,Pbl),c,_,pt') = 
  1284 		   generate1 thy tac_ (Uistate, e_ctxt) (p,Pbl) pt
  1285 	   in ([(tac,tac_,((p,Pbl),(Uistate, e_ctxt)))], c, (pt',(p,Pbl))):calcstate' end
  1286 	       
  1287 	 | Err msg => 
  1288 	   (*TODO.WN03 pass error-msgs to the frontend..
  1289              FIXME ..and dont abuse a tactic for that purpose*)
  1290 	   ([(Tac msg,
  1291 	      Tac_ (Thy_Info.get_theory "Pure", msg,msg,msg),
  1292 	      (e_pos', (e_istate, e_ctxt)))], [], ptp) 
  1293     end
  1294 
  1295 (* val sel = "#Find"; val (p,_) = p; val Add_Find' ct = nxt;
  1296    val (_,_,f,nxt',_,pt')= nxt_specif_additem sel ct (p,Met) c pt;
  1297   *)
  1298   | nxt_specif_additem sel ct (ptp as (pt,(p,Met))) = 
  1299     let
  1300       val (PblObj{meth=met,origin=(oris,(dI',pI',mI'),_),
  1301 		  probl=pbl,spec=(dI,pI,mI),...}) = get_obj I pt p;
  1302       val thy = if dI = e_domID then assoc_thy dI' else assoc_thy dI;
  1303       val cmI = if mI = e_metID then mI' else mI;
  1304       val ctxt = get_ctxt pt (p,Met);
  1305     in case appl_add ctxt sel oris met ((#ppc o get_met) cmI) ct of
  1306       Add itm (*..union old input *) =>
  1307 	let (* val Add itm = appl_add thy sel oris met (#ppc (get_met cmI)) ct;
  1308                *)
  1309 	  val met' = insert_ppc thy itm met;
  1310 	  val (tac,tac_) = 
  1311 	      case sel of
  1312 		  "#Given" => (Add_Given    ct, Add_Given'   (ct, met'))
  1313 		| "#Find"  => (Add_Find     ct, Add_Find'    (ct, met'))
  1314 		| "#Relate"=> (Add_Relation ct, Add_Relation'(ct, met'))
  1315 	  val ((p,Met),c,_,pt') = 
  1316 	      generate1 thy tac_ (Uistate, e_ctxt) (p,Met) pt
  1317 	in ([(tac,tac_,((p,Met), (Uistate, e_ctxt)))], c, (pt',(p,Met))) end
  1318 
  1319     | Err msg => ([(*tacis*)], [], ptp) 
  1320     (*nxt_me collects tacis until not hide; here just no progress*)
  1321     end;
  1322 
  1323 (* ori
  1324 val (msg,itm) = appl_add thy sel oris ppc ct;
  1325 val (Cor(d,ts)) = #5 itm;
  1326 map (atomty) ts;
  1327 
  1328 pre
  1329 *)
  1330 fun ori2Coritm pbt ((i,v,f,d,ts):ori) =
  1331     (i,v,true,f, Cor ((d,ts),(((snd o snd o the o (find_first (eq1 d))) pbt) 
  1332 			      handle _ => error ("ori2Coritm: dsc "^
  1333 						term2str d^
  1334 						"in ori, but not in pbt")
  1335 			      ,ts))):itm;
  1336 fun ori2Coritm (pbt:pat list) ((i,v,f,d,ts):ori) =
  1337     ((i,v,true,f, Cor ((d,ts),((snd o snd o the o 
  1338 			       (find_first (eq1 d))) pbt,ts))):itm)
  1339     handle _ => (*dsc in oris, but not in pbl pat list: keep this dsc*)
  1340     ((i,v,true,f, Cor ((d,ts),(d,ts))):itm);
  1341 
  1342 
  1343 (*filter out oris which have same description in itms*)
  1344 fun filter_outs oris [] = oris
  1345   | filter_outs oris (i::itms) = 
  1346     let val ors = filter_out ((curry op= ((d_in o #5) (i:itm))) o 
  1347 			      (#4:ori -> term)) oris;
  1348     in filter_outs ors itms end;
  1349 
  1350 fun memI a b = member op = a b;
  1351 (*filter oris which are in pbt, too*)
  1352 fun filter_pbt oris pbt =
  1353     let val dscs = map (fst o snd) pbt
  1354     in filter ((memI dscs) o (#4: ori -> term)) oris end;
  1355 
  1356 (*.combine itms from pbl + met and complete them wrt. pbt.*)
  1357 (*FIXXXME.WN031205 complete_metitms doesnt handle incorrect itms !*)
  1358 local infix mem;
  1359 fun x mem [] = false
  1360   | x mem (y :: ys) = x = y orelse x mem ys;
  1361 in 
  1362 fun complete_metitms (oris:ori list) (pits:itm list) (mits:itm list) met = 
  1363 (* val met = (#ppc o get_met) ["DiffApp","max_by_calculus"];
  1364    *)
  1365     let val vat = max_vt pits;
  1366         val itms = pits @ 
  1367 		   (filter ((curry (op mem) vat) o (#2:itm -> int list)) mits);
  1368 	val ors = filter ((curry (op mem) vat) o (#2:ori -> int list)) oris;
  1369         val os = filter_outs ors itms;
  1370     (*WN.12.03?: does _NOT_ add itms from met ?!*)
  1371     in itms @ (map (ori2Coritm met) os) end
  1372 end;
  1373 
  1374 
  1375 
  1376 (*.complete model and guard of a calc-head .*)
  1377 local infix mem;
  1378 fun x mem [] = false
  1379   | x mem (y :: ys) = x = y orelse x mem ys;
  1380 in 
  1381 fun complete_mod_ (oris, mpc, ppc, probl) =
  1382     let	val pits = filter_out ((curry op= false) o (#3: itm -> bool)) probl
  1383 	val vat = if probl = [] then 1 else max_vt probl
  1384 	val pors = filter ((curry (op mem) vat) o (#2:ori -> int list)) oris
  1385 	val pors = filter_outs pors pits (*which are in pbl already*)
  1386         val pors = (filter_pbt pors ppc) (*which are in pbt, too*)
  1387 
  1388 	val pits = pits @ (map (ori2Coritm ppc) pors)
  1389 	val mits = complete_metitms oris pits [] mpc
  1390     in (pits, mits) end
  1391 end;
  1392 
  1393 fun some_spec ((odI, opI, omI):spec) ((dI, pI, mI):spec) =
  1394     (if dI = e_domID then odI else dI,
  1395      if pI = e_pblID then opI else pI,
  1396      if mI = e_metID then omI else mI):spec;
  1397 
  1398 
  1399 (*.find a next applicable tac (for calcstate) and update ptree
  1400  (for ev. finding several more tacs due to hide).*)
  1401 (*FIXXXME: unify ... fun nxt_specif = nxt_spec + applicable_in + specify !!*)
  1402 (*WN.24.10.03        ~~~~~~~~~~~~~~   -> tac     -> tac_      -> -"- as arg*)
  1403 (*WN.24.10.03        fun nxt_solv   = ...................................??*)
  1404 fun nxt_specif (tac as Model_Problem) (pt, pos as (p,p_)) =
  1405   let
  1406     val PblObj {origin = (oris, ospec, _), probl, spec, ...} = get_obj I pt p
  1407     val (dI, pI, mI) = some_spec ospec spec
  1408     val thy = assoc_thy dI
  1409     val mpc = (#ppc o get_met) mI (* just for reuse complete_mod_ *)
  1410     val {cas, ppc, ...} = get_pbt pI
  1411     val pbl = init_pbl ppc (* fill in descriptions *)
  1412     (*----------------if you think, this should be done by the Dialog 
  1413      in the java front-end, search there for WN060225-modelProblem----*)
  1414     val (pbl, met) = 
  1415       case cas of NONE => (pbl, [])
  1416 			| _ => complete_mod_ (oris, mpc, ppc, probl)
  1417     (*----------------------------------------------------------------*)
  1418     val tac_ = Model_Problem' (pI, pbl, met)
  1419     val (pos,c,_,pt) = generate1 thy tac_ (Uistate, e_ctxt) pos pt
  1420   in ([(tac,tac_, (pos, (Uistate, e_ctxt)))], c, (pt,pos)):calcstate' end
  1421 
  1422   | nxt_specif (Add_Given ct) ptp = nxt_specif_additem "#Given" ct ptp
  1423   | nxt_specif (Add_Find  ct) ptp = nxt_specif_additem "#Find"  ct ptp
  1424   | nxt_specif (Add_Relation ct) ptp = nxt_specif_additem"#Relate" ct ptp
  1425 
  1426     (*. called only if no_met is specified .*)     
  1427   | nxt_specif (Refine_Tacitly pI) (ptp as (pt, pos as (p,_))) =
  1428     let val (PblObj {origin = (oris, (dI,_,_),_), ...}) = get_obj I pt p
  1429 	val opt = refine_ori oris pI
  1430     in case opt of
  1431 	   SOME pI' => 
  1432 	   let val {met,ppc,...} = get_pbt pI'
  1433 	       val pbl = init_pbl ppc
  1434 	       (*val pt = update_pbl pt p pbl ..done by Model_Problem*)
  1435 	       val mI = if length met = 0 then e_metID else hd met
  1436                val thy = assoc_thy dI
  1437 	       val (pos,c,_,pt) = 
  1438 		   generate1 thy (Refine_Tacitly' (pI,pI',dI,mI,(*pbl*)[])) 
  1439 			     (Uistate, e_ctxt) pos pt
  1440 	   in ([(Refine_Tacitly pI, Refine_Tacitly' (pI,pI',dI,mI,(*pbl*)[]),
  1441 		 (pos, (Uistate, e_ctxt)))], c, (pt,pos)) end
  1442 	 | NONE => ([], [], ptp)
  1443     end
  1444 
  1445   | nxt_specif (Refine_Problem pI) (ptp as (pt, pos as (p,_))) =
  1446     let val (PblObj {origin=(_,(dI,_,_),_),spec=(dI',_,_),
  1447 		     probl, ...}) = get_obj I pt p
  1448 	val thy = if dI' = e_domID then dI else dI'
  1449     in case refine_pbl (assoc_thy thy) pI probl of
  1450 	   NONE => ([], [], ptp)
  1451 	 | SOME (rfd as (pI',_)) => 
  1452 	   let val (pos,c,_,pt) = 
  1453 		   generate1 (assoc_thy thy) 
  1454 			     (Refine_Problem' rfd) (Uistate, e_ctxt) pos pt
  1455 	    in ([(Refine_Problem pI, Refine_Problem' rfd,
  1456 			    (pos, (Uistate, e_ctxt)))], c, (pt,pos)) end
  1457     end
  1458 
  1459   | nxt_specif (Specify_Problem pI) (pt, pos as (p,_)) =
  1460     let val (PblObj {origin=(oris,(dI,_,_),_),spec=(dI',pI',_),
  1461 		     probl, ...}) = get_obj I pt p;
  1462 	val thy = assoc_thy (if dI' = e_domID then dI else dI');
  1463         val {ppc,where_,prls,...} = get_pbt pI
  1464 	val pbl as (_,(itms,_)) = 
  1465 	    if pI'=e_pblID andalso pI=e_pblID
  1466 	    then (false, (init_pbl ppc, []))
  1467 	    else match_itms_oris thy probl (ppc,where_,prls) oris(*FIXXXXXME?*)
  1468 	(*FIXXXME~~~~~~~~~~~~~~~: take pbl and compare with new pI WN.8.03*)
  1469 	val ((p,Pbl),c,_,pt)= 
  1470 	    generate1 thy (Specify_Problem' (pI, pbl)) (Uistate, e_ctxt) pos pt
  1471     in ([(Specify_Problem pI, Specify_Problem' (pI, pbl),
  1472 		    (pos,(Uistate, e_ctxt)))], c, (pt,pos)) end
  1473 
  1474   (*transfers oris (not required in pbl) to met-model for script-env
  1475     FIXME.WN.8.03: application of several mIDs to SAME model?*)
  1476   | nxt_specif (Specify_Method mID) (ptp as (pt, pos as (p,_))) = 
  1477   let val (PblObj {origin=(oris,(dI',pI',mI'),_), probl=pbl, spec=(dI,pI,mI), 
  1478 		   meth=met, ...}) = get_obj I pt p;
  1479     val {ppc,pre,prls,...} = get_met mID
  1480     val thy = assoc_thy dI
  1481     val oris = add_field' thy ppc oris;
  1482     val dI'' = if dI=e_domID then dI' else dI;
  1483     val pI'' = if pI = e_pblID then pI' else pI;
  1484     val met = if met=[] then pbl else met;(*WN0602 what if more itms in met?*)
  1485     val (ok, (itms, pre')) = match_itms_oris thy met (ppc,pre,prls ) oris;
  1486     val (pos,c,_,pt)= 
  1487 	generate1 thy (Specify_Method' (mID, oris, itms)) (Uistate, e_ctxt) pos pt
  1488   in ([(Specify_Method mID, Specify_Method' (mID, oris, itms),
  1489 		  (pos,(Uistate, e_ctxt)))], c, (pt,pos)) end    
  1490 
  1491   | nxt_specif (Specify_Theory dI) (pt, pos as (p,Pbl)) =
  1492     let val (dI',_,_) = get_obj g_spec pt p
  1493 	val (pos,c,_,pt) = 
  1494 	    generate1 (assoc_thy "Isac") (Specify_Theory' dI) 
  1495 		      (Uistate, e_ctxt) pos pt
  1496     in  (*FIXXXME: check if pbl can still be parsed*)
  1497 	([(Specify_Theory dI, Specify_Theory' dI, (pos,(Uistate, e_ctxt)))], c,
  1498 	 (pt, pos)) end
  1499 
  1500   | nxt_specif (Specify_Theory dI) (pt, pos as (p,Met)) =
  1501     let val (dI',_,_) = get_obj g_spec pt p
  1502 	val (pos,c,_,pt) = 
  1503 	    generate1 (assoc_thy "Isac") (Specify_Theory' dI) 
  1504 		      (Uistate, e_ctxt) pos pt
  1505     in  (*FIXXXME: check if met can still be parsed*)
  1506 	([(Specify_Theory dI, Specify_Theory' dI, (pos,(Uistate, e_ctxt)))], c,
  1507 	 (pt, pos)) end
  1508 
  1509   | nxt_specif m' _ = 
  1510     error ("nxt_specif: not impl. for "^tac2str m');
  1511 
  1512 (*.get the values from oris; handle the term list w.r.t. penv.*)
  1513 
  1514 local infix mem;
  1515 fun x mem [] = false
  1516   | x mem (y :: ys) = x = y orelse x mem ys;
  1517 in 
  1518 fun vals_of_oris oris =
  1519     ((map (mkval' o (#5:ori -> term list))) o 
  1520      (filter ((curry (op mem) 1) o (#2:ori -> int list)))) oris
  1521 end;
  1522 
  1523 
  1524 
  1525 (* create a calc-tree with oris via an cas.refined pbl *)
  1526 fun nxt_specify_init_calc (([], (dI, pI, mI)) : fmz) =
  1527       if pI <> [] 
  1528       then (* from pbl-browser or from CAS cmd with pI=[e_pblID] *)
  1529 	      let 
  1530           val {cas, met, ppc, thy, ...} = get_pbt pI
  1531 	        val dI = if dI = "" then theory2theory' thy else dI
  1532 	        val thy = assoc_thy dI
  1533 	        val mI = if mI = [] then hd met else mI
  1534 	        val hdl = case cas of NONE => pblterm dI pI | SOME t => t
  1535 	        val (pt,_) = cappend_problem e_ptree [] (Uistate, e_ctxt) ([], (dI,pI,mI))
  1536 					  ([], (dI,pI,mI), hdl)
  1537 	        val pt = update_spec pt [] (dI,pI,mI)
  1538 	        val pits = init_pbl' ppc
  1539 	        val pt = update_pbl pt [] pits
  1540 	      in ((pt, ([] ,Pbl)), []) : calcstate end
  1541       else 
  1542         if mI <> [] 
  1543         then (* from met-browser *)
  1544 	        let 
  1545             val {ppc,...} = get_met mI
  1546 	          val dI = if dI = "" then "Isac" else dI
  1547 	          val thy = assoc_thy dI;
  1548 	          val (pt,_) = cappend_problem e_ptree [] (e_istate, e_ctxt) ([], (dI,pI,mI))
  1549 					    ([], (dI,pI,mI), e_term(*FIXME met*));
  1550 	          val pt = update_spec pt [] (dI,pI,mI);
  1551 	          val mits = init_pbl' ppc;
  1552 	          val pt = update_met pt [] mits;
  1553 	        in ((pt, ([], Met)), []) : calcstate end
  1554         else (* new example, pepare for interactive modeling *)
  1555 	        let val (pt,_) = cappend_problem e_ptree [] (e_istate, e_ctxt) ([], e_spec)
  1556 					  ([], e_spec, e_term)
  1557 	        in ((pt, ([], Pbl)), []) : calcstate end
  1558 
  1559   | nxt_specify_init_calc (fmz : fmz_, (dI, pI, mI) : spec) = 
  1560       let           (* both """"""""""""""""""""""""" either empty or complete *)
  1561 	      val thy = assoc_thy dI
  1562 	      val (pI, (pors, pctxt), mI) = 
  1563 	        if mI = ["no_met"] 
  1564 	        then 
  1565             let 
  1566               val (pors, pctxt) = get_pbt pI |> #ppc |> prep_ori fmz thy;
  1567 		          val pI' = refine_ori' pors pI;
  1568 		        in (pI', (pors(*refinement over models with diff.precond only*), pctxt),
  1569 		           (hd o #met o get_pbt) pI') end
  1570 	        else (pI, get_pbt pI |> #ppc |> prep_ori fmz thy, mI)
  1571 	      val {cas, ppc, thy=thy',...} = get_pbt pI (*take dI from _refined_ pbl*)
  1572 	      val dI = theory2theory' (maxthy thy thy');
  1573 	      val hdl = 
  1574           case cas of
  1575 		        NONE => pblterm dI pI
  1576 		      | SOME t => subst_atomic ((vars_of_pbl_' ppc) ~~~ vals_of_oris pors) t;
  1577         val (pt, _) = cappend_problem e_ptree [] (e_istate, pctxt) (fmz, (dI, pI, mI))
  1578 				  (pors, (dI, pI, mI), hdl)
  1579       in ((pt, ([], Pbl)), 
  1580         fst3 (nxt_specif Model_Problem (pt, ([], Pbl)))) : calcstate
  1581       end;
  1582 
  1583 
  1584 
  1585 (*18.12.99*)
  1586 fun get_spec_form (m:tac_) ((p,p_):pos') (pt:ptree) = 
  1587 (*  case appl_spec p pt m of           /// 19.1.00
  1588     Notappl e => Error' (Error_ e)
  1589   | Appl => 
  1590 *)    let val (_,_,f,_,_,_) = specify m (p,p_) [] pt
  1591       in f end;
  1592 
  1593 
  1594 (*fun tag_form thy (formal, given) = cterm_of thy
  1595 	      (((head_of o term_of) given) $ (term_of formal)); WN100819*)
  1596 fun tag_form thy (formal, given) =
  1597     (let val gf = (head_of given) $ formal;
  1598          val _ = cterm_of thy gf
  1599      in gf end)
  1600     handle _ =>
  1601            error ("calchead.tag_form: " ^ 
  1602                   Print_Mode.setmp [] (Syntax.string_of_term
  1603                                            (thy2ctxt thy)) given ^ " .. " ^
  1604                   Print_Mode.setmp [] (Syntax.string_of_term
  1605                                            (thy2ctxt thy)) formal ^
  1606                   " ..types do not match");
  1607 (* val formal = (the o (parse thy)) "[R::real]";
  1608 > val given = (the o (parse thy)) "fixed_values (cs::real list)";
  1609 > tag_form thy (formal, given);
  1610 val it = "fixed_values [R]" : cterm
  1611 *)
  1612 
  1613 fun chktyp thy (n, fs, gs) = 
  1614   ((writeln o (Print_Mode.setmp [] (Syntax.string_of_term (thy2ctxt thy))) o
  1615     (nth n)) fs;
  1616    (writeln o (Print_Mode.setmp [] (Syntax.string_of_term (thy2ctxt thy))) o
  1617     (nth n)) gs;
  1618    tag_form thy (nth n fs, nth n gs));
  1619 fun chktyps thy (fs, gs) = map (tag_form thy) (fs ~~ gs);
  1620 
  1621 (* #####################################################
  1622    find the failing item:
  1623 > val n = 2;
  1624 > val tag__form = chktyp (n,formals,givens);
  1625 > (type_of o term_of o (nth n)) formals; 
  1626 > (type_of o term_of o (nth n)) givens;
  1627 > atomty ((term_of o (nth n)) formals);
  1628 > atomty ((term_of o (nth n)) givens);
  1629 > atomty (term_of tag__form);
  1630 > use_thy"isa-98-1-HOL-plus/knowl-base/DiffAppl";
  1631  ##################################################### *)
  1632 
  1633 (* #####################################################
  1634    testdata setup
  1635 val origin = ["sqrt(9+4*x)=sqrt x + sqrt(5+x)","x::rat","(+0)"];
  1636 val formals = map (the o (parse thy)) origin;
  1637 
  1638 val given  = ["equation (lhs=rhs)",
  1639 	     "bound_variable bdv",   (* TODO type *) 
  1640 	     "error_bound apx"];
  1641 val where_ = ["e is_root_equation_in bdv",
  1642 	      "bdv is_var",
  1643 	      "apx is_const_expr"];
  1644 val find   = ["L::rat set"];
  1645 val with_  = ["L = {bdv. || ((%x. lhs) bdv) - ((%x. rhs) bdv) || < apx}"];
  1646 val chkpbl = map (the o (parse thy)) (given @ where_ @ find @ with_);
  1647 val givens = map (the o (parse thy)) given;
  1648 
  1649 val tag__forms = chktyps (formals, givens);
  1650 map ((atomty) o term_of) tag__forms;
  1651  ##################################################### *)
  1652 
  1653 
  1654 (* check pbltypes, announces one failure a time *)
  1655 (*fun chk_vars ctppc = 
  1656   let val {Given=gi,Where=wh,Find=fi,With=wi,Relate=re} = 
  1657     appc flat (mappc (vars o term_of) ctppc)
  1658   in if (wh\\gi) <> [] then ("wh\\gi",wh\\gi)
  1659      else if (re\\(gi union fi)) <> [] 
  1660 	    then ("re\\(gi union fi)",re\\(gi union fi))
  1661 	  else ("ok",[]) end;*)
  1662 fun chk_vars ctppc = 
  1663   let val {Given=gi,Where=wh,Find=fi,With=wi,Relate=re} = 
  1664           appc flat (mappc vars ctppc)
  1665       val chked = subtract op = gi wh
  1666   in if chked <> [] then ("wh\\gi", chked)
  1667      else let val chked = subtract op = (union op = gi fi) re
  1668           in if chked  <> []
  1669 	     then ("re\\(gi union fi)", chked)
  1670 	     else ("ok", []) 
  1671           end
  1672   end;
  1673 
  1674 (* check a new pbltype: variables (Free) unbound by given, find*) 
  1675 fun unbound_ppc ctppc =
  1676   let val {Given=gi,Find=fi,Relate=re,...} = 
  1677     appc flat (mappc vars ctppc)
  1678   in distinct (*re\\(gi union fi)*) 
  1679               (subtract op = (union op = gi fi) re) end;
  1680 (*
  1681 > val org = {Given=["[R=(R::real)]"],Where=[],
  1682 	   Find=["[A::real]"],With=[],
  1683 	   Relate=["[A=#2*a*b - a^^^#2, (a//#2)^^^#2 + (b//#2)^^^#2 = R^^^#2]"]
  1684 	   }:string ppc;
  1685 > val ctppc = mappc (the o (parse thy)) org;
  1686 > unbound_ppc ctppc;
  1687 val it = [("a","RealDef.real"),("b","RealDef.real")] : (string * typ) list
  1688 *)
  1689 
  1690 
  1691 (* f, a binary operator, is nested rightassociative *)
  1692 fun foldr1 f xs =
  1693   let
  1694     fun fld f (x::[]) = x
  1695       | fld f (x::x'::[]) = f (x',x)
  1696       | fld f (x::x'::xs) = f (fld f (x'::xs),x);
  1697   in ((fld f) o rev) xs end;
  1698 (*
  1699 > val (SOME ct) = parse thy "[a=b,c=d,e=f]";
  1700 > val ces = map (cterm_of thy) (isalist2list (term_of ct));
  1701 > val conj = foldr1 HOLogic.mk_conj (isalist2list (term_of ct));
  1702 > cterm_of thy conj;
  1703 val it = "(a = b & c = d) & e = f" : cterm
  1704 *)
  1705 
  1706 (* f, a binary operator, is nested leftassociative *)
  1707 fun foldl1 f (x::[]) = x
  1708   | foldl1 f (x::x'::[]) = f (x,x')
  1709   | foldl1 f (x::x'::xs) = f (x,foldl1 f (x'::xs));
  1710 (*
  1711 > val (SOME ct) = parse thy "[a=b,c=d,e=f,g=h]";
  1712 > val ces = map (cterm_of thy) (isalist2list (term_of ct));
  1713 > val conj = foldl1 HOLogic.mk_conj (isalist2list (term_of ct));
  1714 > cterm_of thy conj;
  1715 val it = "a = b & c = d & e = f & g = h" : cterm
  1716 *)
  1717 
  1718 
  1719 (* called only once, if a Subproblem has been located in the script*)
  1720 fun nxt_model_pbl (Subproblem'((_,pblID,metID),_,_,_,_)) ptp =
  1721 (* val (Subproblem'((_,pblID,metID),_,_,_,_),ptp) = (m', (pt,(p,p_)));
  1722    *)
  1723     (case metID of
  1724 	 ["no_met"] => 
  1725 	 (snd3 o hd o fst3) (nxt_specif (Refine_Tacitly pblID) ptp)
  1726        | _ => (snd3 o hd o fst3) (nxt_specif Model_Problem ptp))
  1727   (*all stored in tac_ itms     ^^^^^^^^^^*)
  1728   | nxt_model_pbl tac_ _ = 
  1729     error ("nxt_model_pbl: called by tac= "^tac_2str tac_);
  1730 (* run subp_rooteq.sml ''
  1731    until nxt=("Subproblem",Subproblem ("SqRoot",["univariate","equation"]))
  1732 > val (_, (Subproblem'((_,pblID,metID),_,_,_,_),_,_,_,_,_)) =
  1733       (last_elem o drop_last) ets'';
  1734 > val mst = (last_elem o drop_last) ets'';
  1735 > nxt_model_pbl mst;
  1736 val it = Refine_Tacitly ["univariate","equation"] : tac
  1737 *)
  1738 
  1739 (*fun eq1 d (_,(d',_)) = (d = d'); ---modspec.sml*)
  1740 fun eq4 v (_,vts,_,_,_) = member op = vts v;
  1741 fun eq5 (_,_,_,_,itm_) (_,_,_,d,_) = d_in itm_ = d;
  1742 
  1743  
  1744 (*.complete _NON_empty calc-head for autocalc (sub-)pbl from oris
  1745   + met from fmz; assumes pos on PblObj, meth = [].*)
  1746 fun complete_mod (pt, pos as (p, p_):pos') =
  1747 (* val (pt, (p, _)) = (pt, p);
  1748    val (pt, (p, _)) = (pt, pos);
  1749    *)
  1750     let val _= if p_ <> Pbl 
  1751 	       then tracing("###complete_mod: only impl.for Pbl, called with "^
  1752 			    pos'2str pos) else ()
  1753 	val (PblObj{origin=(oris, ospec, hdl), probl, spec,...}) =
  1754 	    get_obj I pt p
  1755 	val (dI,pI,mI) = some_spec ospec spec
  1756 	val mpc = (#ppc o get_met) mI
  1757 	val ppc = (#ppc o get_pbt) pI
  1758 	val (pits, mits) = complete_mod_ (oris, mpc, ppc, probl)
  1759         val pt = update_pblppc pt p pits
  1760 	val pt = update_metppc pt p mits
  1761     in (pt, (p,Met):pos') end
  1762 ;
  1763 (*| complete_mod (pt, pos as (p, Met):pos') =
  1764     error ("###complete_mod: only impl.for Pbl, called with "^
  1765 		 pos'2str pos);*)
  1766 
  1767 (*.complete _EMPTY_ calc-head for autocalc (sub-)pbl from oris(+met from fmz);
  1768    oris and spec (incl. pbl-refinement) given from init_calc or SubProblem .*)
  1769 fun all_modspec (pt, (p,_):pos') =
  1770   let 
  1771     val (PblObj{fmz=(fmz_,_), origin=(pors, spec as (dI,pI,mI), hdl),
  1772 		  ...}) = get_obj I pt p;
  1773     val thy = assoc_thy dI;
  1774 	  val {ppc,...} = get_met mI;
  1775 	  val mors = prep_ori fmz_ thy ppc |> #1;
  1776     val pt = update_pblppc pt p (map (ori2Coritm ppc) pors);
  1777 	  val pt = update_metppc pt p (map (ori2Coritm ppc) mors);
  1778 	  val pt = update_spec pt p (dI,pI,mI);
  1779   in (pt, (p,Met): pos') end;
  1780 
  1781 (*WN.12.03: use in nxt_spec, too ? what about variants ???*)
  1782 fun is_complete_mod_ ([]: itm list) = false
  1783   | is_complete_mod_ itms = 
  1784     foldl and_ (true, (map #3 itms));
  1785 
  1786 fun is_complete_mod (pt, pos as (p, Pbl): pos') =
  1787       if (is_pblobj o (get_obj I pt)) p 
  1788       then (is_complete_mod_ o (get_obj g_pbl pt)) p
  1789       else error ("is_complete_mod: called by PrfObj at "^pos'2str pos)
  1790   | is_complete_mod (pt, pos as (p, Met)) = 
  1791       if (is_pblobj o (get_obj I pt)) p 
  1792       then (is_complete_mod_ o (get_obj g_met pt)) p
  1793       else error ("is_complete_mod: called by PrfObj at "^pos'2str pos)
  1794   | is_complete_mod (_, pos) =
  1795       error ("is_complete_mod called by " ^ pos'2str pos ^
  1796 		    " (should be Pbl or Met)");
  1797 
  1798 (*.have (thy, pbl, met) _all_ been specified explicitly ?.*)
  1799 fun is_complete_spec (pt, pos as (p,_): pos') = 
  1800     if (not o is_pblobj o (get_obj I pt)) p 
  1801     then error ("is_complete_spec: called by PrfObj at "^pos'2str pos)
  1802     else let val (dI,pI,mI) = get_obj g_spec pt p
  1803 	 in dI<>e_domID andalso pI<>e_pblID andalso mI<>e_metID end;
  1804 (*.complete empty items in specification from origin (pbl, met ev.refined);
  1805   assumes 'is_complete_mod'.*)
  1806 fun complete_spec (pt, pos as (p,_): pos') = 
  1807     let val PblObj {origin = (_,ospec,_), spec,...} = get_obj I pt p
  1808 	val pt = update_spec pt p (some_spec ospec spec)
  1809     in (pt, pos) end;
  1810 
  1811 fun is_complete_modspec ptp = 
  1812     is_complete_mod ptp andalso is_complete_spec ptp;
  1813 
  1814 
  1815 
  1816 
  1817 fun pt_model (PblObj {meth,spec,origin=(_,spec',hdl),...}) Met =
  1818 (* val ((PblObj {meth,spec,origin=(_,spec',hdl),...}), Met) = (ppobj, p_);
  1819    *)
  1820     let val (_,_,metID) = get_somespec' spec spec'
  1821 	val pre = 
  1822 	    if metID = e_metID then []
  1823 	    else let val {prls,pre=where_,...} = get_met metID
  1824 		     val pre = check_preconds' prls where_ meth 0
  1825 		 in pre end
  1826 	val allcorrect = is_complete_mod_ meth
  1827 			 andalso foldl and_ (true, (map #1 pre))
  1828     in ModSpec (allcorrect, Met, hdl, meth, pre, spec) end
  1829   | pt_model (PblObj {probl,spec,origin=(_,spec',hdl),...}) _(*Frm,Pbl*) =
  1830 (* val ((PblObj {probl,spec,origin=(_,spec',hdl),...}),_) = (ppobj, p_);
  1831    *)
  1832     let val (_,pI,_) = get_somespec' spec spec'
  1833 	val pre =
  1834 	    if pI = e_pblID then []
  1835 	    else let val {prls,where_,cas,...} = get_pbt pI
  1836 		     val pre = check_preconds' prls where_ probl 0
  1837 		 in pre end
  1838 	val allcorrect = is_complete_mod_ probl
  1839 			 andalso foldl and_ (true, (map #1 pre))
  1840     in ModSpec (allcorrect, Pbl, hdl, probl, pre, spec) end;
  1841 
  1842 
  1843 fun pt_form (PrfObj {form,...}) = Form form
  1844   | pt_form (PblObj {probl,spec,origin=(_,spec',_),...}) =
  1845     let val (dI, pI, _) = get_somespec' spec spec'
  1846 	val {cas,...} = get_pbt pI
  1847     in case cas of
  1848 	   NONE => Form (pblterm dI pI)
  1849 	 | SOME t => Form (subst_atomic (mk_env probl) t)
  1850     end;
  1851 (*vvv takes the tac _generating_ the formula=result, asm ok....
  1852 fun pt_result (PrfObj {result=(t,asm), tac,...}) = 
  1853     (Form t, 
  1854      if null asm then NONE else SOME asm, 
  1855      SOME tac)
  1856   | pt_result (PblObj {result=(t,asm), origin = (_,ospec,_), spec,...}) =
  1857     let val (_,_,metID) = some_spec ospec spec
  1858     in (Form t, 
  1859 	if null asm then NONE else SOME asm, 
  1860 	if metID = e_metID then NONE else SOME (Apply_Method metID)) end;
  1861 -------------------------------------------------------------------------*)
  1862 
  1863 
  1864 (*.pt_extract returns
  1865       # the formula at pos
  1866       # the tactic applied to this formula
  1867       # the list of assumptions generated at this formula
  1868 	(by application of another tac to the preceding formula !)
  1869    pos is assumed to come from the frontend, ie. generated by moveDown.*)
  1870 (*cannot be in ctree.sml, because ModSpec has to be calculated*)
  1871 fun pt_extract (pt,([],Res)) =
  1872 (* val (pt,([],Res)) = ptp;
  1873    *)
  1874     (* ML_TODO 110417 get assumptions from ctxt !? *)
  1875     let val (f, asm) = get_obj g_result pt []
  1876     in (Form f, NONE, asm) end
  1877 (* val p = [3,2];
  1878    *)
  1879   | pt_extract (pt,(p,Res)) =
  1880 (* val (pt,(p,Res)) = ptp;
  1881    *)
  1882     let val (f, asm) = get_obj g_result pt p
  1883 	val tac = if last_onlev pt p
  1884 		  then if is_pblobj' pt (lev_up p)
  1885 		       then let val (PblObj{spec=(_,pI,_),...}) = 
  1886 				    get_obj I pt (lev_up p)
  1887 			    in if pI = e_pblID then NONE 
  1888 			       else SOME (Check_Postcond pI) end
  1889 		       else SOME End_Trans (*WN0502 TODO for other branches*)
  1890 		  else let val p' = lev_on p
  1891 		       in if is_pblobj' pt p'
  1892 			  then let val (PblObj{origin = (_,(dI,pI,_),_),...}) =
  1893 				       get_obj I pt p'
  1894 			       in SOME (Subproblem (dI, pI)) end
  1895 			  else if f = get_obj g_form pt p'
  1896 			  then SOME (get_obj g_tac pt p')
  1897 			  (*because this Frm          ~~~is not on worksheet*)
  1898 			  else SOME (Take (term2str (get_obj g_form pt p')))
  1899 		       end
  1900     in (Form f, tac, asm) end
  1901 	
  1902   | pt_extract (pt, pos as (p,p_(*Frm,Pbl*))) =
  1903 (* val (pt, pos as (p,p_(*Frm,Pbl*))) = ptp;
  1904    val (pt, pos as (p,p_(*Frm,Pbl*))) = (pt, p);
  1905    *)
  1906     let val ppobj = get_obj I pt p
  1907 	val f = if is_pblobj ppobj then pt_model ppobj p_
  1908 		else get_obj pt_form pt p
  1909 	val tac = g_tac ppobj
  1910     in (f, SOME tac, []) end;
  1911 
  1912 
  1913 (**. get the formula from a ctree-node:
  1914  take form+res from PblObj and 1.PrfObj and (PrfObj after PblObj)
  1915  take res from all other PrfObj's .**)
  1916 (*designed for interSteps, outcommented 04 in favour of calcChangedEvent*)
  1917 fun formres p (Nd (PblObj {origin = (_,_, h), result = (r, _),...}, _)) =
  1918     [("headline", (p, Frm), h), 
  1919      ("stepform", (p, Res), r)]
  1920   | formres p (Nd (PrfObj {form, result = (r, _),...}, _)) = 
  1921     [("stepform", (p, Frm), form), 
  1922      ("stepform", (p, Res), r)];
  1923 
  1924 fun form p (Nd (PrfObj {result = (r, _),...}, _)) = 
  1925     [("stepform", (p, Res), r)]
  1926 
  1927 (*assumes to take whole level, in particular hd -- for use in interSteps*)
  1928 fun get_formress fs p [] = flat fs
  1929   | get_formress fs p (nd::nds) =
  1930     (* start with   'form+res'       and continue with trying 'res' only*)
  1931     get_forms (fs @ [formres p nd]) (lev_on p) nds
  1932 and get_forms fs p [] = flat fs
  1933   | get_forms fs p (nd::nds) =
  1934     if is_pblnd nd
  1935     (* start again with      'form+res' ///ugly repeat with Check_elementwise
  1936     then get_formress (fs @ [formres p nd]) (lev_on p) nds                   *)
  1937     then get_forms    (fs @ [formres p nd]) (lev_on p) nds
  1938     (* continue with trying 'res' only*)
  1939     else get_forms    (fs @ [form    p nd]) (lev_on p) nds;
  1940 
  1941 (**.get an 'interval' 'from' 'to' of formulae from a ptree.**)
  1942 (*WN050219 made robust against _'to' below or after Complete nodes
  1943 	   by handling exn caused by move_dn*)
  1944 (*WN0401 this functionality belongs to ctree.sml, 
  1945 but fetching a calc_head requires calculations defined in modspec.sml
  1946 transfer to ME/me.sml !!!
  1947 WN051224 ^^^ doesnt hold any longer, since only the headline of a calc_head
  1948 is returned !!!!!!!!!!!!!
  1949 *)
  1950 fun eq_pos' (p1,Frm) (p2,Frm) = p1 = p2
  1951   | eq_pos' (p1,Res) (p2,Res) = p1 = p2
  1952   | eq_pos' (p1,Pbl) (p2,p2_) = p1 = p2 andalso (case p2_ of
  1953 						     Pbl => true
  1954 						   | Met => true
  1955 						   | _ => false)
  1956   | eq_pos' (p1,Met) (p2,p2_) = p1 = p2 andalso (case p2_ of
  1957 						     Pbl => true
  1958 						   | Met => true
  1959 						   | _ => false)
  1960   | eq_pos' _ _ = false;
  1961 
  1962 (*.get an 'interval' from the ctree; 'interval' is w.r.t. the 
  1963    total ordering Position#compareTo(Position p) in the java-code
  1964 val get_interval = fn
  1965     : pos' ->     : from is "move_up 1st-element" to return
  1966       pos' -> 	  : to the last element to be returned; from < to
  1967       int -> 	  : level: 0 gets the flattest sub-tree possible
  1968 			   >999 gets the deepest sub-tree possible
  1969       ptree -> 	  : 
  1970       (pos' * 	  : of the formula
  1971        Term.term) : the formula
  1972 	  list
  1973 .*)
  1974 fun get_interval from to level pt =
  1975 (* val (from,level) = (f,lev);
  1976    val (from, to, level) = (([3, 2, 1], Res), ([],Res), 9999);
  1977    *)
  1978     let fun get_inter c (from:pos') (to:pos') lev pt =
  1979 (* val (c, from, to, lev) = ([], from, to, level);
  1980    ------for recursion.......
  1981    val (c, from:pos', to:pos') = (c @ [(from, f)], move_dn [] pt from, to);
  1982    *)
  1983 	    if eq_pos' from to orelse from = ([],Res)
  1984 	    (*orelse ... avoids Exception- PTREE "end of calculation" raised,
  1985 	     if 'to' has values NOT generated by move_dn, see systest/me.sml
  1986              TODO.WN0501: introduce an order on pos' and check "from > to"..
  1987              ...there is an order in Java! 
  1988              WN051224 the hack got worse with returning term instead ptform*)
  1989 	    then let val (f,_,_) = pt_extract (pt, from)
  1990 		 in case f of
  1991 			ModSpec (_,_,headline,_,_,_) => c @ [(from, headline)] 
  1992 		      | Form t => c @ [(from, t)]
  1993 		 end
  1994 	    else 
  1995 		if lev < lev_of from
  1996 		then (get_inter c (move_dn [] pt from) to lev pt)
  1997 		     handle (PTREE _(*from move_dn too far*)) => c
  1998 		else let val (f,_,_) = pt_extract (pt, from)
  1999 			 val term = case f of
  2000 					ModSpec (_,_,headline,_,_,_)=> headline
  2001 				      | Form t => t
  2002 		     in (get_inter (c @ [(from, term)]) 
  2003 				   (move_dn [] pt from) to lev pt)
  2004 			handle (PTREE _(*from move_dn too far*)) 
  2005 			       => c @ [(from, term)] end
  2006     in get_inter [] from to level pt end;
  2007 
  2008 (*for tests*)
  2009 fun posform2str (pos:pos', form) =
  2010     "("^ pos'2str pos ^", "^
  2011     (case form of 
  2012 	 Form f => term2str f
  2013        | ModSpec c => term2str (#3 c(*the headline*)))
  2014     ^")";
  2015 fun posforms2str pfs = (strs2str' o (map (curry op ^ "\n")) o 
  2016 			(map posform2str)) pfs;
  2017 fun posterm2str (pos:pos', t) =
  2018     "("^ pos'2str pos ^", "^term2str t^")";
  2019 fun posterms2str pfs = (strs2str' o (map (curry op ^ "\n")) o 
  2020 			(map posterm2str)) pfs;
  2021 
  2022 
  2023 (*WN050225 omits the last step, if pt is incomplete*)
  2024 fun show_pt pt = 
  2025     tracing (posterms2str (get_interval ([],Frm) ([],Res) 99999 pt));
  2026 
  2027 (*.get a calchead from a PblObj-node in the ctree; 
  2028    preconditions must be calculated.*)
  2029 fun get_ocalhd (pt, pos' as (p,Pbl):pos') = 
  2030     let val PblObj {origin = (oris, ospec, hdf'), spec, probl,...} = 
  2031 	    get_obj I pt p
  2032 	val {prls,where_,...} = get_pbt (#2 (some_spec ospec spec))
  2033 	val pre = check_preconds (assoc_thy"Isac") prls where_ probl
  2034     in (ocalhd_complete probl pre spec, Pbl, hdf', probl, pre, spec):ocalhd end
  2035 | get_ocalhd (pt, pos' as (p,Met):pos') = 
  2036     let val PblObj {fmz = fmz as (fmz_,_), origin = (oris, ospec, hdf'), 
  2037 		    spec, meth,...} = 
  2038 	    get_obj I pt p
  2039 	val {prls,pre,...} = get_met (#3 (some_spec ospec spec))
  2040 	val pre = check_preconds (assoc_thy"Isac") prls pre meth
  2041     in (ocalhd_complete meth pre spec, Met, hdf', meth, pre, spec):ocalhd end;
  2042 
  2043 (*.at the activeFormula set the Model, the Guard and the Specification 
  2044    to empty and return a CalcHead;
  2045    the 'origin' remains (for reconstructing all that).*)
  2046 fun reset_calchead (pt, pos' as (p,_):pos') = 
  2047     let val PblObj {origin = (_, _, hdf'),...} = get_obj I pt p
  2048 	val pt = update_pbl pt p []
  2049 	val pt = update_met pt p []
  2050 	val pt = update_spec pt p e_spec
  2051     in (pt, (p,Pbl):pos') end;
  2052 
  2053 end
  2054 
  2055 open CalcHead;