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