src/Tools/isac/CalcElements/calcelems.sml
author Walther Neuper <walther.neuper@jku.at>
Wed, 15 Jan 2020 11:47:38 +0100
changeset 59767 c4acd312bd53
parent 59728 f47a69ee4504
child 59771 cdf6417d27b6
permissions -rw-r--r--
preps for IJCAR paper
     1 (* ~~/src/Tools/isac/calcelems.sml
     2    elements of calculations.
     3    they are partially held in association lists as ref's for
     4    switching language levels (meta-string, object-values).                                  
     5    in order to keep these ref's during re-evaluation of code,                       
     6    they are defined here at the beginning of the code.
     7    Author: Walther Neuper 2003
     8    (c) copyright due to lincense terms
     9 *)
    10 
    11 signature CALC_ELEMENT =
    12   sig
    13     type cas_elem
    14     type pbt
    15     type ptyps
    16     type metID
    17     type pblID
    18     type mets
    19     type met
    20     datatype 'a ptyp = Ptyp of string * 'a list * 'a ptyp list
    21 
    22     type authors
    23     type guh
    24 
    25     type fillpat
    26     datatype thydata
    27       = Hcal of {calc: Rule.calc, coursedesign: authors, guh: guh, mathauthors: authors}
    28       | Hord of {coursedesign: authors, guh: guh, mathauthors: authors, ord: Rule.subst -> term * term -> bool}
    29       | Hrls of {coursedesign: authors, guh: guh, mathauthors: authors, thy_rls: Rule.thyID * Rule.rls}
    30       | Hthm of {coursedesign: authors, fillpats: fillpat list, guh: guh, mathauthors: authors, thm: thm}
    31       | Html of {coursedesign: authors, guh: guh, html: string, mathauthors: authors}
    32     type theID
    33     type rlss_elem
    34     val merge_rlss: rlss_elem list * rlss_elem list -> rlss_elem list
    35     val rls_eq: (''a * ('b * 'c)) * (''a * ('d * 'e)) -> bool
    36     type spec
    37     val cas_eq: cas_elem * cas_elem -> bool
    38     val e_Ptyp: pbt ptyp
    39     val merge_ptyps: 'a ptyp list * 'a ptyp list -> 'a ptyp list
    40     val check_guhs_unique: bool Unsynchronized.ref
    41     val check_pblguh_unique: guh -> pbt ptyp list -> unit
    42     val insrt: pblID -> 'a -> string list -> 'a ptyp list -> 'a ptyp list
    43     val e_Mets: met ptyp
    44     val check_metguh_unique: guh -> met ptyp list -> unit
    45     val add_thydata: string list * string list -> thydata -> thydata ptyp list -> thydata ptyp list
    46     val get_py: 'a ptyp list -> pblID -> string list -> 'a
    47     val update_hthm: thydata -> fillpat list -> thydata
    48     val update_ptyps: string list -> string list -> 'a -> 'a ptyp list -> 'a ptyp list
    49     val part2guh: theID -> guh
    50     val spec2str: string * string list * string list -> string
    51     val linefeed: string -> string
    52     val pbts2str: pbt list -> string
    53     val thes2str: thydata list -> string
    54     val theID2str: string list -> string
    55     val the2str: thydata -> string
    56     val trace_calc: bool Unsynchronized.ref
    57     eqtype thmID
    58     type thm'
    59     val trace_rewrite: bool Unsynchronized.ref
    60     val depth: int Unsynchronized.ref
    61     val assoc_thy: Rule.theory' -> theory
    62     type thm''
    63     val metID2str: string list -> string
    64     val e_pblID: pblID
    65     val e_metID: metID
    66     val empty_spec: spec
    67     val e_spec: spec
    68     datatype ketype = Exp_ | Met_ | Pbl_ | Thy_
    69     type kestoreID
    70     val app_py: 'a ptyp list -> ('a ptyp -> 'b) -> pblID -> string list -> 'b
    71     val ketype2str: ketype -> string
    72     val coll_pblguhs: pbt ptyp list -> guh list
    73     val coll_metguhs: met ptyp list -> guh list
    74     type pat
    75     val pats2str: pat list -> string
    76     val maxthy: theory -> theory -> theory
    77     eqtype filename
    78     val lim_deriv: int Unsynchronized.ref
    79     val id_of_thm: Rule.rule -> string
    80     val isabthys: unit -> theory list
    81     val thyID_of_derivation_name: string -> string
    82     val partID': Rule.theory' -> string
    83     val thm2guh: string * Rule.thyID -> thmID -> guh
    84     val thmID_of_derivation_name: string -> string
    85     val rls2guh: string * Rule.thyID -> Rule.rls' -> guh
    86     val theID2guh: theID -> guh
    87     eqtype fillpatID
    88     type pbt_ = string * (term * term)
    89     eqtype xml
    90     val cal2guh: string * Rule.thyID -> string -> guh
    91     val ketype2str': ketype -> string
    92     val str2ketype': string -> ketype
    93     val thmID_of_derivation_name': thm -> string
    94     eqtype filepath
    95     val theID2thyID: theID -> Rule.thyID
    96     val thy2guh: theID -> guh
    97     val thypart2guh: theID -> guh
    98     val ord2guh: string * Rule.theory' -> string -> string
    99     val update_hrls: thydata -> Rule.errpatID list -> thydata
   100     eqtype iterID
   101     eqtype calcID
   102     val thm''_of_thm: thm -> thm''
   103     val thm_of_thm: Rule.rule -> thm
   104 (* ---- for tests only: shifted from below to remove the Warning "unused" at fun.def. --------- *)
   105     val pats2str' : pat list -> string
   106     val insert_fillpats: thydata ptyp list -> (pblID * fillpat list) list -> thydata ptyp list ->
   107       thydata ptyp list
   108 (*/-------------------------------------------------------- ! aktivate for Test_Isac BEGIN ---\* )
   109     val knowthys: unit -> theory list
   110     val e_pbt: pbt
   111     val e_met: met
   112 ( *\--- ! aktivate for Test_Isac END ----------------------------------------------------------/*)
   113 
   114 (*----- unused code, kept as hints to design ideas ---------------------------------------------*)
   115 val overwritelthy: theory -> (Rule.rls' * (string * Rule.rls)) list * (Rule.rls' * Rule.rls) list ->
   116   (Rule.rls' * (string * Rule.rls)) list  end
   117 
   118 (**)
   119 structure Celem(**): CALC_ELEMENT(**) =
   120 struct
   121 (**)
   122 
   123 val linefeed = (curry op^) "\n"; (* ?\<longrightarrow> libraryC ?*)
   124 type authors = string list;
   125 
   126 type iterID = int;
   127 type calcID = int;
   128 
   129 (* TODO CLEANUP Thm:
   130 Thm (string, thm): (a) needs string to identify sym_thmID for handling in front-end;
   131                    (b) investigate if ""RS sym" attaches a [.]" still occurs: string_of_thmI
   132 thmID            : type for data from user input + program
   133 thmDeriv         : type for thy_hierarchy ONLY
   134 obsolete types   : thm' (SEE "ad thm'"), thm''. 
   135 revise funs      : id_of_thm, thm_of_thm, rep_thm_G', eq_thmI, eq_thmI', thm''_of_thm thm.
   136 activate         : thmID_of_derivation_name'
   137 *)
   138 type thmID = string;    (* identifier for a thm (the shortest possible identifier)       *)
   139 type thmDeriv = string; (* WN120524 deprecated
   140   thyID ^"."^ xxx ^"."^ thmID, see fun thmID_of_derivation_name 
   141   WN120524: dont use Thm.derivation_name, this is destroyed by num_str;
   142   Thm.get_name_hint survives num_str and seems perfectly reliable *)
   143 
   144 type thm' = thmID * Rule.cterm';(*WN060610 deprecated in favour of thm''; WN180324: used in TODO review:
   145 val thm'2xml : int -> Celem.thm' -> Celem.xml
   146 val assoc_thm': theory -> Celem.thm' -> thm
   147 | Calculate' of Rule.theory' * string * term * (term * Celem.thm')
   148 *)
   149 (* tricky combination of (string, term) for theorems in Isac:
   150   * case 1 general: frontend + lucin, e.g. applicable_in..Rewrite: (thmID, _) --> (thmID, thm)
   151     by Global_Theory.get_thm, special cases ("add_commute",..) see convert_metaview_to_thmid.
   152   * case 2 "sym_..": Global_Theory.get_thm..RS sym
   153   * case 3 ad-hoc thm "#..." mk_thm from ad-hoc term (numerals only) in calculate_:
   154     from applicable_in..Calculate: opstr --calculate_/adhoc_thm--> (thmID, thm)
   155 *)
   156 type thm'' = thmID * thm; (* only for transport via libisabelle isac-java <--- ME *)
   157 
   158 (*.a 'guh'='globally unique handle' is a string unique for each element
   159    of isac's KEStore and persistent over time
   160    (in particular under shifts within the respective hierarchy);
   161    specialty for thys:
   162    # guh NOT resistant agains shifts from one thy to another
   163    (which is the price for Isabelle's design: thy's overwrite ids of subthy's)
   164    # requirement for matchTheory: induce guh from tac + current thy
   165    (see 'fun thy_containing_thm', 'fun thy_containing_rls' etc.)
   166    TODO: introduce to pbl, met.*)
   167 type guh = string;
   168 
   169 type xml = string; (* rm together with old code replaced by XML.tree *)
   170 
   171 
   172 (* for (at least) 2 kinds of access:
   173   (1) given an errpatID, find the respective fillpats (e.g. in fun find_fill_pats)
   174   (2) given a thm, find respective fillpats *)
   175 type fillpatID = string
   176 type fillpat =
   177   fillpatID   (* DESIGN ?TODO: give an order w.r.t difficulty ? *)
   178   * term      (* the pattern with fill-in gaps                  *)
   179   * Rule.errpatID; (* which the fillpat would be a help for          
   180                  DESIGN ?TODO: list for several patterns ?      *)
   181 
   182 
   183 (*  WN0509 discussion:
   184 #############################################################################
   185 #   How to manage theorys in subproblems wrt. the requirement,              #
   186 #   that scripts should be re-usable ?                                      #
   187 #############################################################################
   188 
   189     eg. 'Program Solve_rat_equation' calls 'SubProblem (RatEq',..'
   190     which would not allow to 'solve (y'' = -M_b / EI, M_b)' by this script
   191     because Biegelinie.thy is subthy of RatEq.thy and thus Biegelinie.M_b
   192     is unknown in RatEq.thy and M_b cannot be parsed into the scripts guard
   193     (see match_ags).
   194 
   195     Preliminary solution:
   196     # the thy in 'SubProblem (thy', pbl, arglist)' is not taken automatically,
   197     # instead the 'maxthy (rootthy pt) thy' is taken for each subpbl
   198     # however, a thy specified by the user in the rootpbl may lead to
   199       errors in far-off subpbls (which are not yet reported properly !!!)
   200       and interactively specifiying thys in subpbl is not very relevant.
   201 
   202     Other solutions possible:
   203     # always parse and type-check with Thy_Info_get_theory "Isac_Knowledge"
   204       (rejected due to the vague idea eg. to re-use equations for R in C etc.)
   205     # regard the subthy-relation in specifying thys of subpbls
   206     # specifically handle 'SubProblem (undefined, pbl, arglist)'
   207     # ???
   208 *)
   209 
   210 fun id_of_thm (Rule.Thm (id, _)) = id  (* TODO re-arrange code for rule2str *)
   211   | id_of_thm _ = raise ERROR ("id_of_thm: uncovered case " (* ^ rule2str r *))
   212 fun thm_of_thm (Rule.Thm (_, thm)) = thm  (* TODO re-arrange code for rule2str *)
   213   | thm_of_thm _ = raise ERROR ("thm_of_thm: uncovered case " (* ^ rule2str r *))
   214 
   215 fun thmID_of_derivation_name dn = last_elem (space_explode "." dn);
   216 fun thmID_of_derivation_name' thm = (thmID_of_derivation_name o Thm.get_name_hint) thm
   217 fun thyID_of_derivation_name dn = hd (space_explode "." dn);
   218 fun thm''_of_thm thm = (thmID_of_derivation_name' thm, thm) : thm''
   219 
   220 
   221 
   222 (*the key into the hierarchy ob theory elements*)
   223 type theID = string list;
   224 val theID2str = strs2str; (*theID eg. is ["IsacKnowledge", "Test", "Rulesets", "ac_plus_times"]*)
   225 fun theID2thyID theID =
   226   if length theID >= 3 then (last_elem o (drop_last_n 2)) theID
   227   else error ("theID2thyID called with " ^ theID2str theID);
   228 
   229 (*the key into the hierarchy ob problems*)
   230 type pblID = string list; (* domID :: ...*)
   231 val e_pblID = ["e_pblID"];
   232 
   233 (*the key into the hierarchy ob methods*)
   234 type metID = string list;
   235 type spec = Rule.domID * pblID * metID;
   236 fun spec2str (dom, pbl, met) = 
   237   "(" ^ quote dom  ^ ", " ^ strs2str pbl ^ ", " ^ strs2str met ^ ")";
   238 val e_metID = ["e_metID"];
   239 val metID2str = strs2str;
   240 val empty_spec = (Rule.e_domID, e_pblID, e_metID);
   241 val e_spec = empty_spec;
   242 
   243 (* association list with cas-commands, for generating a complete calc-head *)
   244 type generate_fn = 
   245   (term list ->    (* the arguments of the cas-command, eg. (x+1=2, x) *)
   246     (term *        (* description of an element                        *)
   247       term list)   (* value of the element (always put into a list)    *)
   248   list)            (* of elements in the formalization                 *)
   249 type cas_elem = 
   250   (term *          (* cas-command, eg. 'solve'                         *)
   251     (spec *        (* theory, problem, method                          *)
   252     generate_fn))
   253 fun cas_eq ((t1, (_, _)) : cas_elem, (t2, (_, _)) : cas_elem) = t1 = t2
   254 
   255 (*either theID or pblID or metID*)
   256 type kestoreID = string list;
   257 
   258 (* for distinction of contexts WN130621: disambiguate with Isabelle's Context *)
   259 datatype ketype = Exp_ | Thy_ | Pbl_ | Met_;
   260 fun ketype2str Exp_ = "Exp_"
   261   | ketype2str Thy_ = "Thy_"
   262   | ketype2str Pbl_ = "Pbl_"
   263   | ketype2str Met_ = "Met_";
   264 fun ketype2str' Exp_ = "Example"
   265   | ketype2str' Thy_ = "Theory"
   266   | ketype2str' Pbl_ = "Problem"
   267   | ketype2str' Met_ = "Method";
   268 (* for conversion from XML *)
   269 fun str2ketype' "exp" = Exp_
   270   | str2ketype' "thy" = Thy_
   271   | str2ketype' "pbl" = Pbl_
   272   | str2ketype' "met" = Met_
   273   | str2ketype' str = raise ERROR ("str2ketype': WRONG arg = " ^ str)
   274 
   275 (* A tree for storing data defined in different theories 
   276   for access from the Interpreter and from dialogue authoring 
   277   using a string list as key.
   278   'a is for pbt | met | thydata; after WN030424 naming "pbt" became inappropriate *)
   279 datatype 'a ptyp =
   280 	Ptyp of string *  (* element of the key                                                        *)
   281 		'a list *       (* several pbts with different domIDs/thy TODO: select by subthy (isaref.p.69)
   282 			                 presently only _ONE_ elem FOR ALL KINDS OF CONTENT pbt | met | thydata    *)
   283 		('a ptyp) list; (* the children nodes                                                        *)
   284 
   285 (* datatype for collecting thydata for hierarchy *)
   286 (*WN060720 more consistent naming would be 'type thyelem' or 'thelem'*)
   287 datatype thydata =
   288   Html of {guh: guh, coursedesign: authors, mathauthors: authors, html: string}
   289 | Hthm of {guh: guh, coursedesign: authors, mathauthors: authors, fillpats: fillpat list,
   290    thm: thm} (* here no sym_thm, thus no thmID required *)
   291 | Hrls of {guh: guh, coursedesign: authors, mathauthors: authors, thy_rls: (Rule.thyID * Rule.rls)}
   292 | Hcal of {guh: guh, coursedesign: authors, mathauthors: authors, calc: Rule.calc}
   293 | Hord of {guh: guh, coursedesign: authors, mathauthors: authors,
   294     ord: (Rule.subst -> (term * term) -> bool)};
   295 fun the2str (Html {guh, ...}) = guh
   296   | the2str (Hthm {guh, ...}) = guh
   297   | the2str (Hrls {guh, ...}) = guh
   298   | the2str (Hcal {guh, ...}) = guh
   299   | the2str (Hord {guh, ...}) = guh
   300 fun thes2str thes = map the2str thes |> list2str;
   301 
   302 (* notes on thehier concerning sym_thmID theorems (created in derivations, reverse rewriting)
   303      (a): thehier does not contain sym_thmID theorems
   304      (b): lookup for sym_thmID directly from Isabelle using sym_thm
   305           (within math-engine NO lookup in thehier -- within java in *.xml only!)
   306 TODO (c): export from thehier to xml
   307 TODO (c1)   creates one entry for "thmID" (and NONE for "sym_thmID") in the hierarchy
   308 TODO (c2)   creates 2 files "thy_*-thm-thmID.xml" and "thy_*-thm-sym_thmID.xml"
   309 TODO (d): 1 entry in the MiniBrowser's hierarchy (generated from xml)
   310           stands for both, "thmID" and "sym_thmID" 
   311 TODO (d1)   lookup from calctxt          
   312 TODO (d1)   lookup from from rule set in MiniBrowser *)
   313 type thehier = (thydata ptyp) list;
   314 (* required to determine sequence of main nodes of thehier in KEStore.thy *)
   315 fun part2guh [str] = (case str of
   316 	  "Isabelle" => "thy_isab_" ^ str ^ "-part" : guh
   317       | "IsacScripts" => "thy_scri_" ^ str ^ "-part"
   318       | "IsacKnowledge" => "thy_isac_" ^ str ^ "-part"
   319       | str => raise ERROR ("thy2guh: called with \""^ str ^"\""))
   320   | part2guh theID = raise ERROR ("part2guh called with theID = \"" ^ theID2str theID ^ "'");
   321 
   322 fun thy2guh [part, thyID] = (case part of
   323       "Isabelle" => "thy_isab_" ^ thyID
   324     | "IsacScripts" => "thy_scri_" ^ thyID
   325     | "IsacKnowledge" => "thy_isac_" ^ thyID
   326     | str => raise ERROR ("thy2guh: called with \"" ^ str ^ "\""))
   327   | thy2guh theID = raise ERROR ("thy2guh called with \"" ^ strs2str' theID ^ "\"");
   328 			
   329 fun thypart2guh ([part, thyID, thypart] : theID) = (case part of
   330       "Isabelle" => "thy_isab_" ^ thyID ^ "-" ^ thypart : guh
   331     | "IsacScripts" => "thy_scri_" ^ thyID ^ "-" ^ thypart
   332     | "IsacKnowledge" => "thy_isac_" ^ thyID ^ "-" ^ thypart
   333     | str => raise ERROR ("thypart2guh: called with '" ^ str ^ "'"))
   334   | thypart2guh strs = raise ERROR ("thypart2guh called with \"" ^ strs2str' strs ^ "\"");
   335 
   336   
   337 (* convert the data got via contextToThy to a globally unique handle.
   338    there is another way to get the guh: get out of the 'theID' in the hierarchy *)
   339 fun thm2guh (isa, thyID) thmID = case isa of
   340     "Isabelle" => "thy_isab_" ^ Rule.theory'2thyID thyID ^ "-thm-" ^ strip_thy thmID : guh
   341   | "IsacKnowledge" => "thy_isac_" ^ Rule.theory'2thyID thyID ^ "-thm-" ^ strip_thy thmID
   342   | "IsacScripts" => "thy_scri_" ^ Rule.theory'2thyID thyID ^ "-thm-" ^ strip_thy thmID
   343   | _ => raise ERROR
   344     ("thm2guh called with (isa, thyID) = (" ^ isa ^ ", " ^ thyID ^ ") for thm = \"" ^ thmID ^ "\"");
   345 
   346 fun rls2guh (isa, thyID) rls' = case isa of
   347     "Isabelle" => "thy_isab_" ^ Rule.theory'2thyID thyID ^ "-rls-" ^ rls' : guh
   348   | "IsacKnowledge" => "thy_isac_" ^ Rule.theory'2thyID thyID ^ "-rls-" ^ rls'
   349   | "IsacScripts" => "thy_scri_" ^ Rule.theory'2thyID thyID ^ "-rls-" ^ rls'
   350   | _ => raise ERROR
   351     ("rls2guh called with (isa, thyID) = (" ^ isa ^ ", " ^ thyID ^ ") for rls = \"" ^ rls' ^ "\"");
   352 			  
   353 fun cal2guh (isa, thyID) calID = case isa of
   354     "Isabelle" => "thy_isab_" ^ Rule.theory'2thyID thyID ^ "-cal-" ^ calID : guh
   355   | "IsacKnowledge" => "thy_isac_" ^ Rule.theory'2thyID thyID ^ "-cal-" ^ calID
   356   | "IsacScripts" => "thy_scri_" ^ Rule.theory'2thyID thyID ^ "-cal-" ^ calID
   357   | _ => raise ERROR
   358     ("cal2guh called with (isa, thyID) = (" ^ isa ^ ", " ^ thyID ^ ") for cal = \"" ^ calID ^ "\"");
   359 			  
   360 fun ord2guh (isa, thyID) rew_ord' = case isa of
   361     "Isabelle" => "thy_isab_" ^ Rule.theory'2thyID thyID ^ "-ord-" ^ rew_ord' : guh
   362   | "IsacKnowledge" => "thy_isac_" ^ Rule.theory'2thyID thyID ^ "-ord-" ^ rew_ord'
   363   | "IsacScripts" => "thy_scri_" ^ Rule.theory'2thyID thyID ^ "-ord-" ^ rew_ord'
   364   | _ => raise ERROR
   365     ("ord2guh called with (isa, thyID) = (" ^ isa ^ ", " ^ thyID ^ ") for ord = \"" ^ rew_ord' ^ "\"");
   366 
   367 (* not only for thydata, but also for thy's etc *)
   368 (* TODO
   369 fun theID2guh theID = case length theID of
   370     0 => error ("theID2guh: called with theID = " ^ strs2str' theID)
   371   | 1 => part2guh theID
   372   | 2 => thy2guh theID
   373   | 3 => thypart2guh theID
   374   | 4 => 
   375     let val [isa, thyID, typ, elemID] = theID
   376     in case typ of
   377         "Theorems" => thm2guh (isa, thyID) elemID
   378       | "Rulesets" => rls2guh (isa, thyID) elemID
   379       | "Calculations" => cal2guh (isa, thyID) elemID
   380       | "Orders" => ord2guh (isa, thyID) elemID
   381       | "Theorems" => thy2guh [isa, thyID]
   382       | str => raise ERROR ("theID2guh: called with theID = " ^ strs2str' theID)
   383     end
   384   | n => raise ERROR ("theID2guh called with theID = " ^ strs2str' theID);
   385 *)
   386 (* not only for thydata, but also for thy's etc *)
   387 fun theID2guh [] = raise ERROR ("theID2guh: called with []")
   388   | theID2guh [str] = part2guh [str]
   389   | theID2guh [s1, s2] = thy2guh [s1, s2]
   390   | theID2guh [s1, s2, s3] = thypart2guh [s1, s2, s3]
   391   | theID2guh (strs as [isa, thyID, typ, elemID]) = (case typ of
   392       "Theorems" => thm2guh (isa, thyID) elemID
   393     | "Rulesets" => rls2guh (isa, thyID) elemID
   394     | "Calculations" => cal2guh (isa, thyID) elemID
   395     | "Orders" => ord2guh (isa, thyID) elemID
   396     | _ => raise ERROR ("theID2guh: called with theID = " ^ strs2str' strs))
   397   | theID2guh strs = raise ERROR ("theID2guh called with theID = " ^ strs2str' strs);
   398 
   399 type filepath = string;
   400 type filename = string;
   401 
   402 
   403 
   404 (* datastructure for KEStore_Elems, intermediate for thehier *)
   405 type rlss_elem = 
   406   (Rule.rls' *     (* identifier unique within Isac *)
   407   (Rule.theory' *  (* just for assignment in thehier, not appropriate for parsing etc *)
   408     Rule.rls))     (* ((#id o rep_rls) rls) = rls'   by coding discipline *)
   409 fun rls_eq ((id1, (_, _)), (id2, (_, _))) = id1 = id2
   410 
   411 fun insert_merge_rls (re as (id, (thyID, r1)) : rlss_elem) ys = 
   412     case get_index (fn y => if curry rls_eq re y then SOME y else NONE) ys of
   413       NONE => re :: ys
   414     | SOME (i, (_, (_, r2))) => 
   415       let
   416         val r12 = Rule.merge_rls id r1 r2
   417       in list_update ys i (id, (thyID, r12)) end
   418 fun merge_rlss (s1, s2) = fold insert_merge_rls s1 s2;
   419 
   420 
   421 fun assoc_thy thy =
   422     if thy = "e_domID"
   423     then (Rule.Thy_Info_get_theory "Base_Tools") (*lower bound of Knowledge*)
   424     else (Rule.Thy_Info_get_theory thy) handle _ => error ("ME_Isa: thy \"" ^ thy ^ "\" not in system");
   425 
   426 (* overwrite an element in an association list and pair it with a thyID
   427    in order to create the thy_hierarchy;
   428    overwrites existing rls' even if they are defined in a different thy;
   429    this is related to assoc_rls, TODO.WN060120: assoc_rew_ord, assoc_calc *)
   430 (* WN060120 ...these are NOT compatible to "fun assoc_thm'" in that
   431    they do NOT handle overlays by re-using an identifier in different thys;
   432    "thyID.rlsID" would be a good solution, if the "." would be possible
   433    in scripts...
   434    actually a hack to get alltogether run again with minimal effort *)
   435 fun insthy thy' (rls', rls) = (rls', (thy', rls));
   436 fun overwritelthy thy (al, bl: (Rule.rls' * Rule.rls) list) =
   437     let val bl' = map (insthy ((get_thy o Rule.theory2theory') thy)) bl
   438     in overwritel (al, bl') end;
   439 
   440 fun maxthy thy1 thy2 = if Context.subthy (thy1, thy2) then thy2 else thy1;
   441 
   442 
   443 (* trace internal steps of isac's numeral calculations *)
   444 val trace_calc = Unsynchronized.ref false;
   445 (* trace internal steps of isac's rewriter *)
   446 val trace_rewrite = Unsynchronized.ref false;
   447 (* depth of recursion in traces of the rewriter, if trace_rewrite:=true *)
   448 val depth = Unsynchronized.ref 99999;
   449 (* no of rewrites exceeding this int -> NO rewrite *)
   450 val lim_deriv = Unsynchronized.ref 100;
   451 (* switch for checking guhs unique before storing a pbl or met;
   452    set true at startup (done at begin of ROOT.ML)
   453    set false for editing IsacKnowledge (done at end of ROOT.ML) *)
   454 val check_guhs_unique = Unsynchronized.ref true;
   455 
   456 (* the pattern for an item of a problems model or a methods guard *)
   457 type pat =
   458   (string *     (* field               *)
   459 	  (term *     (* description         *)
   460 	     term))   (* id | arbitrary term *);
   461 fun pat2str ((field, (dsc, id)) : pat) = 
   462   pair2str (field, pair2str (Rule.term2str dsc, Rule.term2str id))
   463 fun pats2str pats = (strs2str o (map pat2str)) pats
   464 fun pat2str' ((field, (dsc, id)) : pat) = 
   465   pair2str (field, pair2str (Rule.term2str dsc, Rule.term2str id)) ^ "\n"
   466 fun pats2str' pats = (strs2str o (map pat2str')) pats
   467 
   468 (* types for problems models (TODO rename to specification models) *)
   469 type pbt_ =
   470   (string *   (* field "#Given",..*)(*deprecated due to 'type pat'*)
   471     (term *   (* description      *)
   472       term)); (* id | struct-var  *)
   473 type pbt = 
   474   {guh : guh,         (* unique within this isac-knowledge                               *)
   475   mathauthors : string list, (* copyright                                                *)
   476   init : pblID,       (* to start refinement with                                        *)
   477   thy  : theory,      (* which allows to compile that pbt
   478                         TODO: search generalized for subthy (ref.p.69*)
   479                         (*^^^ WN050912 NOT used during application of the problem,
   480                         because applied terms may be from 'subthy' as well as from super;
   481                         thus we take 'maxthy'; see match_ags !                           *)
   482   cas : term option,  (* 'CAS-command'                                                   *)
   483   prls : Rule.rls,    (* for preds in where_                                             *)
   484   where_ : term list, (* where - predicates                                              *)
   485   ppc : pat list,     (* this is the model-pattern; 
   486                          it contains "#Given","#Where","#Find","#Relate"-patterns
   487                          for constraints on identifiers see "fun cpy_nam"                *)
   488   met : metID list}   (* methods solving the pbt                                         *)
   489 
   490 val e_pbt = {guh = "pbl_empty", mathauthors = [], init = e_pblID, thy = Thy_Info.get_theory "Pure",
   491   cas = NONE, prls = Rule.Erls, where_ = [], ppc = [], met = []} : pbt
   492 fun pbt2str ({cas = cas', guh = guh', init = init', mathauthors = ma', met = met', ppc = ppc',
   493       prls = prls', thy = thy', where_ = w'} : pbt)
   494     = "{cas = " ^ (Rule.termopt2str cas') ^  ", guh = \"" ^ guh'  ^ "\", init = "
   495       ^ (strs2str init') ^ ", mathauthors = " ^ (strs2str ma' |> quote) ^ ", met = "
   496       ^ (strslist2strs met') ^ ", ppc = " ^ pats2str ppc' ^ ", prls = "
   497       ^ (Rule.rls2str prls' |> quote) ^ ", thy = {" ^ (Rule.theory2str thy') ^ "}, where_ = "
   498       ^ (Rule.terms2str w') ^ "}" |> linefeed;
   499 fun pbts2str pbts = map pbt2str pbts |> list2str;
   500 
   501 val e_Ptyp = Ptyp ("e_pblID", [e_pbt], [])
   502 type ptyps = (pbt ptyp) list
   503 
   504 fun coll_pblguhs pbls =
   505   let
   506     fun node coll (Ptyp (_, [n], ns)) = [(#guh : pbt -> guh) n] @ (nodes coll ns)
   507       | node _ _ = raise ERROR "coll_pblguhs - node"
   508 	  and nodes coll [] = coll
   509       | nodes coll (n :: ns) = (node coll n) @ (nodes coll ns);
   510   in nodes [] pbls end;
   511 fun check_pblguh_unique guh pbls =
   512   if member op = (coll_pblguhs pbls) guh
   513   then error ("check_guh_unique failed with \""^ guh ^"\";\n"^
   514 	      "use \"sort_pblguhs()\" for a list of guhs;\n"^
   515 	      "consider setting \"check_guhs_unique := false\"")
   516   else ();
   517 
   518 fun insrt _ pbt [k] [] = [Ptyp (k, [pbt], [])]
   519   | insrt d pbt [k] ((Ptyp (k', [p], ps)) :: pys) =
   520     ((*writeln ("### insert 1: ks = " ^ strs2str [k] ^ "    k'= " ^ k');*)
   521     if k = k'
   522     then ((Ptyp (k', [pbt], ps)) :: pys)
   523     else  ((Ptyp (k', [p], ps)) :: (insrt d pbt [k] pys))
   524     )			 
   525   | insrt d pbt (k::ks) ((Ptyp (k', [p], ps)) :: pys) =
   526     ((*writeln ("### insert 2: ks = "^(strs2str (k::ks))^"    k'= "^k');*)
   527     if k = k'
   528     then ((Ptyp (k', [p], insrt d pbt ks ps)) :: pys)
   529     else 
   530       if length pys = 0
   531       then error ("insert: not found " ^ (strs2str (d : pblID)))
   532       else ((Ptyp (k', [p], ps)) :: (insrt d pbt (k :: ks) pys))
   533     )
   534   | insrt _ _ _ _ = raise ERROR "";
   535 
   536 fun update_ptyps ID _ _ [] =
   537     error ("update_ptyps: " ^ strs2str' ID ^ " does not exist")
   538   | update_ptyps ID [i] data ((py as Ptyp (key, _, pys)) :: pyss) =
   539     if i = key
   540     then 
   541       if length pys = 0
   542       then ((Ptyp (key, [data], [])) :: pyss)
   543       else error ("update_ptyps: " ^ strs2str' ID ^ " has descendants")
   544     else py :: update_ptyps ID [i] data pyss
   545   | update_ptyps ID (i :: is) data ((py as Ptyp (key, d, pys)) :: pyss) =
   546     if i = key
   547     then ((Ptyp (key,  d, update_ptyps ID is data pys)) :: pyss)
   548     else (py :: (update_ptyps ID (i :: is) data pyss))
   549   | update_ptyps _ _ _ _ = raise ERROR "update_ptyps called with undef pattern.";
   550 
   551 (* this function only works wrt. the way Isabelle evaluates Theories and is not a general merge
   552   function for trees / ptyps *)
   553 fun merge_ptyps ([], pt) = pt
   554   | merge_ptyps (pt, []) = pt
   555   | merge_ptyps ((x' as Ptyp (k, _, ps)) :: xs, (xs' as Ptyp (k', y, ps') :: ys)) =
   556       if k = k'
   557       then Ptyp (k, y, merge_ptyps (ps, ps')) :: merge_ptyps (xs, ys)
   558       else x' :: merge_ptyps (xs, xs');
   559 
   560 (* data for methods stored in 'methods'-database*)
   561 type met = 
   562      {guh        : guh,             (* unique within this isac-knowledge             *)
   563       mathauthors: string list,     (* copyright                                     *)
   564       init       : pblID,           (* WN060721 introduced mistakenly--TODO.REMOVE!  *)
   565       rew_ord'   : Rule.rew_ord',   (* for rules in Detail                           
   566 			                                 TODO.WN0509 store fun itself, see 'type pbt'  *)
   567       erls       : Rule.rls,        (* the eval_rls for cond. in rules FIXME "rls'   
   568 				                               instead erls in "fun prep_met"                *)
   569       srls       : Rule.rls,        (* for evaluating list expressions in scr        *)
   570       prls       : Rule.rls,        (* for evaluating predicates in modelpattern     *)
   571       crls       : Rule.rls,        (* for check_elementwise, ie. formulae in calc.  *)
   572       nrls       : Rule.rls,        (* canonical simplifier specific for this met    *)
   573       errpats    : Rule.errpat list,(* error patterns expected in this method        *)
   574       calc       : Rule.calc list,  (* Theory_Data in fun prep_met                   *)
   575       (*branch   : TransitiveB set in append_problem at generation ob pblobj         *)
   576       ppc        : pat list,        (* items in given, find, relate;
   577 	      items (in "#Find") which need not occur in the arg-list of a SubProblem
   578         are 'copy-named' with an identifier "*'.'".
   579         copy-named items are 'generating' if they are NOT "*'''" ?WN120516??
   580         see ME/calchead.sml 'fun is_copy_named'.                                     *)
   581       pre        : term list,       (* preconditions in where                        *)
   582       scr        : Rule.program     (* progam, empty as @{thm refl} or Rfuns         *)
   583 	   };
   584 val e_met = {guh = "met_empty", mathauthors = [], init = e_metID, rew_ord' = "e_rew_ord'",
   585 	erls = Rule.e_rls, srls = Rule.e_rls, prls = Rule.e_rls, calc = [], crls = Rule.e_rls,
   586 	errpats = [], nrls = Rule.e_rls, ppc = [], pre = [], scr = Rule.EmptyScr};
   587 val e_Mets = Ptyp ("e_metID", [e_met],[]);
   588 
   589 type mets = (met ptyp) list;
   590 fun coll_metguhs mets =
   591   let
   592     fun node coll (Ptyp (_, [n], ns)) = [(#guh : met -> guh) n] @ (nodes coll ns)
   593       | node _ _ = raise ERROR "coll_pblguhs - node"
   594   	and nodes coll [] = coll
   595   	  | nodes coll (n :: ns) = (node coll n) @ (nodes coll ns);
   596     in nodes [] mets end;
   597 fun check_metguh_unique (guh:guh) (mets: (met ptyp) list) =
   598     if member op = (coll_metguhs mets) guh
   599     then raise ERROR ("check_guh_unique failed with \"" ^ guh ^"\";\n"^
   600 		(*"use \"sort_metguhs()\" for a list of guhs;\n" ^        ...evaluates to [] ?!?*)
   601 		  "consider setting \"check_guhs_unique := false\"")
   602     else ();
   603 
   604 fun Html_default exist = (Html {guh = theID2guh exist, 
   605   coursedesign = ["isac team 2006"], mathauthors = [], html = ""})
   606 
   607 fun fill_parents (_, [i]) thydata = Ptyp (i, [thydata], [])
   608   | fill_parents (exist, i :: is) thydata =
   609     Ptyp (i, [Html_default (exist @ [i])], [fill_parents (exist @ [i], is) thydata])
   610   | fill_parents _ _ = raise ERROR "Html_default: avoid ML warning: Matches are not exhaustive"
   611 
   612 fun add_thydata (exist, is) thydata [] = [fill_parents (exist, is) thydata]
   613   | add_thydata (exist, [i]) data (pys as (py as Ptyp (key, _, _)) :: pyss) = 
   614     if i = key
   615     then pys (* preserve existing thydata *) 
   616     else py :: add_thydata (exist, [i]) data pyss
   617   | add_thydata (exist, iss as (i :: is)) data ((py as Ptyp (key, d, pys)) :: pyss) = 
   618     if i = key
   619     then       
   620       if length pys = 0
   621       then Ptyp (key, d, [fill_parents (exist @ [i], is) data]) :: pyss
   622       else Ptyp (key, d, add_thydata (exist @ [i], is) data pys) :: pyss
   623     else py :: add_thydata (exist, iss) data pyss
   624   | add_thydata _ _ _ = raise ERROR "add_thydata: avoid ML warning: Matches are not exhaustive"
   625 
   626 fun update_hthm (Hthm {guh, coursedesign, mathauthors, thm, ...}) fillpats' =
   627   Hthm {guh = guh, coursedesign = coursedesign, mathauthors = mathauthors,
   628     fillpats = fillpats', thm = thm}
   629   | update_hthm _ _ = raise ERROR "update_hthm: wrong arguments";
   630 
   631 (* for dialog-authoring *)
   632 fun update_hrls (Hrls {guh, coursedesign, mathauthors, thy_rls = (thyID, rls)}) errpatIDs =
   633     let
   634       val rls' = 
   635         case rls of
   636           Rule.Rls {id, preconds, rew_ord, erls, srls, calc, rules, scr, ...}
   637           => Rule.Rls {id = id, preconds = preconds, rew_ord = rew_ord, erls = erls, srls = srls,
   638                calc = calc, rules = rules, scr = scr, errpatts = errpatIDs}
   639         | Rule.Seq {id, preconds, rew_ord, erls, srls, calc, rules, scr, ...}
   640           => Rule.Seq {id = id, preconds = preconds, rew_ord = rew_ord, erls = erls, srls = srls,
   641                calc = calc, rules = rules, scr = scr, errpatts = errpatIDs}
   642         | Rule.Rrls {id, prepat, rew_ord, erls, calc, scr, ...}
   643           => Rule.Rrls {id = id, prepat = prepat, rew_ord = rew_ord, erls = erls, calc = calc,
   644                scr = scr, errpatts = errpatIDs}
   645         | Erls => Erls
   646     in
   647       Hrls {guh = guh, coursedesign = coursedesign, mathauthors = mathauthors,
   648         thy_rls = (thyID, rls')}
   649     end
   650   | update_hrls _ _ = raise ERROR "update_hrls: wrong arguments";
   651 
   652 fun app_py p f (d:pblID) (k(*:pblRD*)) =
   653   let
   654     fun py_err _ = raise ERROR ("app_py: not found: " ^ strs2str d);
   655     fun app_py' _ [] = py_err ()
   656       | app_py' [] _ = py_err ()
   657       | app_py' [k0] ((p' as Ptyp (k', _, _  )) :: ps) =
   658         if k0 = k' then f p' else app_py' [k0] ps
   659       | app_py' (k' as (k0 :: ks)) (Ptyp (k'', _, ps) :: ps') =
   660         if k0 = k'' then app_py' ks ps else app_py' k' ps';
   661   in app_py' k p end;
   662 fun get_py p =
   663   let
   664     fun extract_py (Ptyp (_, [py], _)) = py
   665       | extract_py _ = raise ERROR ("extract_py: Ptyp has wrong format.");
   666   in app_py p extract_py end;
   667 
   668 fun (*KEStore_Elems.*)insert_fillpats th fis = (* for tests bypassing setup KEStore_Elems *)
   669   let
   670     fun update_elem th (theID, fillpats) =
   671       let
   672         val hthm = get_py th theID theID
   673         val hthm' = update_hthm hthm fillpats
   674           handle ERROR _ => error ("insert_fillpats: " ^ strs2str theID ^ "must address a theorem")
   675       in update_ptyps theID theID hthm' end
   676   in fold (update_elem th) fis end
   677 
   678 (* group the theories defined in Isac, compare Build_Thydata:
   679   section "Get and group the theories defined in Isac" *) 
   680 fun isabthys () = (*["Complex_Main", "Taylor", .., "Pure"]*)
   681   let
   682     val allthys = Theory.ancestors_of (Rule.Thy_Info_get_theory "Build_Thydata")
   683   in
   684     drop ((find_index (curry Context.eq_thy (Thy_Info.get_theory "Complex_Main")) allthys), allthys)
   685   end
   686 fun knowthys () = (*["Isac_Knowledge", .., "Descript"]*)
   687   let
   688     fun isacthys () = (* ["Isac_Knowledge", .., "KEStore"] without Build_Isac thys: "Interpret" etc *)
   689       let
   690         val allthys = filter_out (member Context.eq_thy
   691           [(*Thy_Info_get_theory "ProgLang",*) Rule.Thy_Info_get_theory "Interpret", 
   692             Rule.Thy_Info_get_theory "MathEngine", Rule.Thy_Info_get_theory "BridgeLibisabelle"]) 
   693           (Theory.ancestors_of (Rule.Thy_Info_get_theory "Build_Thydata"))
   694       in
   695         take ((find_index (curry Context.eq_thy (Thy_Info.get_theory "Complex_Main")) allthys), 
   696         allthys)
   697       end
   698     val isacthys' = isacthys ()
   699     val proglang_parent = Rule.Thy_Info_get_theory "ProgLang"
   700   in
   701     take ((find_index (curry Context.eq_thy proglang_parent) isacthys'), isacthys')
   702   end
   703 
   704 fun progthys () = (*["Isac_Knowledge", .., "Descript"]*)
   705   let
   706     fun isacthys () = (* ["Isac_Knowledge", .., "KEStore"] without Build_Isac thys: "Interpret" etc *)
   707       let                                                        
   708         val allthys = filter_out (member Context.eq_thy
   709           [(*Thy_Info_get_theory "ProgLang",*) Rule.Thy_Info_get_theory "Interpret", 
   710             Rule.Thy_Info_get_theory "MathEngine", Rule.Thy_Info_get_theory "BridgeLibisabelle"]) 
   711           (Theory.ancestors_of (Rule.Thy_Info_get_theory "Build_Thydata"))
   712       in
   713         take ((find_index (curry Context.eq_thy (Thy_Info.get_theory "Complex_Main")) allthys), 
   714         allthys)
   715       end
   716     val isacthys' = isacthys ()
   717     val proglang_parent = Rule.Thy_Info_get_theory "ProgLang"
   718   in
   719     drop ((find_index (curry Context.eq_thy proglang_parent) isacthys') + 1(*ProgLang*), isacthys')
   720   end
   721 
   722 fun partID thy = 
   723   if member Context.eq_thy (knowthys ()) thy then "IsacKnowledge"
   724   else if member Context.eq_thy (progthys ()) thy then "IsacScripts"
   725   else if member Context.eq_thy (isabthys ()) thy then "Isabelle"
   726   else error ("closure of thys in Isac is broken by " ^ Rule.string_of_thy thy)
   727 fun partID' thy' = partID (Rule.Thy_Info_get_theory thy')
   728 
   729 end (*struct*)
   730