src/Tools/isac/Interpret/calchead.sml
author Walther Neuper <neuper@ist.tugraz.at>
Mon, 18 Jul 2011 09:29:17 +0200
branchdecompose-isar
changeset 42092 1a6d6089e594
parent 42090 908dfde7cf75
child 42432 7dc25d1526a5
permissions -rw-r--r--
intermed: make autocalc..CompleteCalc run with x+1=2

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