1 (* Title: solve an example by interpreting a method's script
2 Author: Walther Neuper 1999
3 (c) copyright due to lincense terms.
4 1234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
5 10 20 30 40 50 60 70 80 90 100
8 fun safe (ScrState (_,_,_,_,s,_)) = s
9 | safe (RrlsState _) = Safe;
12 type tac'_ = mstID * tac; (*DG <-> ME*)
13 val e_tac'_ = ("Empty_Tac", Empty_Tac):tac'_;
15 fun mk_tac'_ m = case m of
16 Init_Proof (ppc, spec) => ("Init_Proof", Init_Proof (ppc, spec ))
17 | Model_Problem => ("Model_Problem", Model_Problem)
18 | Refine_Tacitly pblID => ("Refine_Tacitly", Refine_Tacitly pblID)
19 | Refine_Problem pblID => ("Refine_Problem", Refine_Problem pblID)
20 | Add_Given cterm' => ("Add_Given", Add_Given cterm')
21 | Del_Given cterm' => ("Del_Given", Del_Given cterm')
22 | Add_Find cterm' => ("Add_Find", Add_Find cterm')
23 | Del_Find cterm' => ("Del_Find", Del_Find cterm')
24 | Add_Relation cterm' => ("Add_Relation", Add_Relation cterm')
25 | Del_Relation cterm' => ("Del_Relation", Del_Relation cterm')
27 | Specify_Theory domID => ("Specify_Theory", Specify_Theory domID)
28 | Specify_Problem pblID => ("Specify_Problem", Specify_Problem pblID)
29 | Specify_Method metID => ("Specify_Method", Specify_Method metID)
30 | Apply_Method metID => ("Apply_Method", Apply_Method metID)
31 | Check_Postcond pblID => ("Check_Postcond", Check_Postcond pblID)
32 | Free_Solve => ("Free_Solve",Free_Solve)
34 | Rewrite_Inst (subs, thm') => ("Rewrite_Inst", Rewrite_Inst (subs, thm'))
35 | Rewrite thm' => ("Rewrite", Rewrite thm')
36 | Rewrite_Asm thm' => ("Rewrite_Asm", Rewrite_Asm thm')
37 | Rewrite_Set_Inst (subs, rls')
38 => ("Rewrite_Set_Inst", Rewrite_Set_Inst (subs, rls'))
39 | Rewrite_Set rls' => ("Rewrite_Set", Rewrite_Set rls')
40 | End_Ruleset => ("End_Ruleset", End_Ruleset)
42 | End_Detail => ("End_Detail", End_Detail)
43 | Detail_Set rls' => ("Detail_Set", Detail_Set rls')
44 | Detail_Set_Inst (s, rls') => ("Detail_Set_Inst", Detail_Set_Inst (s, rls'))
46 | Calculate op_ => ("Calculate", Calculate op_)
47 | Substitute sube => ("Substitute", Substitute sube)
48 | Apply_Assumption cts' => ("Apply_Assumption", Apply_Assumption cts')
50 | Take cterm' => ("Take", Take cterm')
51 | Take_Inst cterm' => ("Take_Inst", Take_Inst cterm')
52 | Group (con, ints) => ("Group", Group (con, ints))
53 | Subproblem (domID, pblID) => ("Subproblem", Subproblem (domID, pblID))
55 | Subproblem_Full(spec,cts')=> ("Subproblem_Full", Subproblem_Full(spec,cts'))
57 | End_Subproblem => ("End_Subproblem",End_Subproblem)
58 | CAScmd cterm' => ("CAScmd", CAScmd cterm')
60 | Split_And => ("Split_And", Split_And)
61 | Conclude_And => ("Conclude_And", Conclude_And)
62 | Split_Or => ("Split_Or", Split_Or)
63 | Conclude_Or => ("Conclude_Or", Conclude_Or)
64 | Begin_Trans => ("Begin_Trans", Begin_Trans)
65 | End_Trans => ("End_Trans", End_Trans)
66 | Begin_Sequ => ("Begin_Sequ", Begin_Sequ)
67 | End_Sequ => ("End_Sequ", Begin_Sequ)
68 | Split_Intersect => ("Split_Intersect", Split_Intersect)
69 | End_Intersect => ("End_Intersect", End_Intersect)
70 | Check_elementwise cterm' => ("Check_elementwise", Check_elementwise cterm')
71 | Or_to_List => ("Or_to_List", Or_to_List)
72 | Collect_Trues => ("Collect_Results", Collect_Trues)
74 | Empty_Tac => ("Empty_Tac",Empty_Tac)
75 | Tac string => ("Tac",Tac string)
76 | User => ("User",User)
77 | End_Proof' => ("End_Proof'",End_Proof');
80 val empty_tac'_ = (mk_tac'_ Empty_Tac):tac'_;
82 fun mk_tac ((_,m):tac'_) = m;
83 fun mk_mstID ((mI,_):tac'_) = mI;
85 fun tac'_2str ((ID,ms):tac'_) = ID ^ (tac2str ms);
86 (* TODO: tac2str, tac'_2str NOT tested *)
90 type squ = ptree; (* TODO: safe etc. *)
92 (*13.9.02--------------
93 type ctr = (loc * pos) list;
94 val ops = [("PLUS","Groups.plus_class.plus"),("MINUS","Groups.minus_class.minus"),("TIMES","Groups.times_class.times"),
95 ("cancel","cancel"),("pow","pow"),("sqrt","sqrt")];
97 @{term "PLUS"}; (*Free ("PLUS", "'a") : term*)
98 @{term "plus"}; (*Const ("Groups.plus_class.plus", "'a => 'a => 'a")*)
99 @{term "MINUS"}; (*Free ("MINUS", "'a")*)
100 @{term "minus"}; (*Const ("Groups.minus_class.minus", "'a => 'a => 'a")*)
101 @{term "TIMES"}; (*Free ("TIMES", "'a")*)
102 @{term "times"}; (*Const ("Groups.times_class.times", "'a => 'a => 'a")*)
103 @{term "CANCEL"}; (*Free ("CANCEL", "'a")*)
104 @{term "cancel"}; (*Free ("cancel", "'a")*)
105 @{term "POWER"}; (*Free ("POWER", "'a")*)
106 @{term "pow"}; (*Free ("pow", "'a")*)
107 @{term "SQRT"}; (*Free ("SQRT", "'a")*)
108 @{term "sqrt"}; (*Const ("NthRoot.sqrt", "RealDef.real => RealDef.real")*)
111 case assoc (ops, op_) of
112 SOME op' => op' | NONE => error ("op_intern: no op= "^op_);
113 -----------------------*)
117 (* use"ME/solve.sml";
120 val ttt = (term_of o the o (parse thy))"Substitute [(bdv,x)] g";
121 val ttt = (term_of o the o (parse thy))"Rewrite thmid True g";
123 Const ("Script.Rewrite'_Inst",_) $ sub $ Free (thm',_) $ Const (pa,_) $ f'
128 val specsteps = ["Init_Proof","Refine_Tacitly","Refine_Problem",
129 "Model_Problem",(*"Match_Problem",*)
130 "Add_Given","Del_Given","Add_Find","Del_Find",
131 "Add_Relation","Del_Relation",
132 "Specify_Theory","Specify_Problem","Specify_Method"];
134 "-----------------------------------------------------------------------";
137 fun step2taci ((tac_, _, pt, p, _):step) = (*FIXXME.040312: redesign step*)
138 (tac_2tac tac_, tac_, (p, get_loc pt p)):taci;
141 (*FIXME.WN050821 compare solve ... nxt_solv*)
142 (* val ("Apply_Method",Apply_Method' (mI,_))=(mI,m);
143 val (("Apply_Method",Apply_Method' (mI,_,_)),pt, pos as (p,_))=(m,pt, pos);
145 fun solve ("Apply_Method", m as Apply_Method' (mI, _, _, _)) (pt:ptree, (pos as (p,_))) =
146 let val {srls, ...} = get_met mI;
147 val PblObj {meth=itms, ...} = get_obj I pt p;
148 val thy' = get_obj g_domID pt p;
149 val thy = assoc_thy thy';
150 val (is as ScrState (env,_,_,_,_,_), ctxt, sc) = init_scrstate thy itms mI;
151 val ini = init_form thy sc env;
156 let val (pos,c,_,pt) =
157 generate1 thy (Apply_Method' (mI, SOME t, is, ctxt))
158 (is, ctxt) (lev_on p, Frm)(*implicit Take*) pt;
159 in ("ok",([(Apply_Method mI, Apply_Method' (mI, SOME t, is, ctxt),
160 ((lev_on p, Frm), (is, ctxt)))], c, (pt,pos)):calcstate')
162 | NONE => (*execute the first tac in the Script, compare solve m*)
164 val (m', (is', ctxt'), _) = next_tac (thy', srls) (pt, (p, Res)) sc (is, ctxt);
165 val d = e_rls (*FIXME: get simplifier from domID*);
167 case locate_gen (thy',srls) m' (pt,(p, Res)) (sc,d) (is', ctxt') of
168 Steps (is'', ss as (m'',f',pt',p',c')::_) =>
169 ("ok", (map step2taci ss, c', (pt',p')))
171 let val (p,ps,f,pt) =
172 generate_hard (assoc_thy "Isac") m (p,Frm) pt;
173 in ("not-found-in-script", ([(tac_2tac m, m, (pos, (is, ctxt)))], ps, (pt,p))) end
177 | solve ("Free_Solve", Free_Solve') (pt,po as (p,_)) =
178 let (*val _=tracing"###solve Free_Solve";*)
179 val p' = lev_dn_ (p,Res);
180 val pt = update_metID pt (par_pblobj pt p) e_metID;
181 in ("ok", ([(Empty_Tac, Empty_Tac_, (po, (Uistate, e_ctxt)))], [], (pt,p')))
184 | solve ("Check_Postcond", Check_Postcond' (pI,_)) (pt,(pos as (p,p_))) =
186 val pp = par_pblobj pt p
188 (case get_obj g_tac pt p of
189 Check_elementwise _ => (*collects and instantiates asms*)
190 (snd o (get_obj g_result pt)) p
191 | _ => get_assumptions_ pt (p,p_))
192 handle _ => [] (*WN.27.5.03 asms in subpbls not completely clear*)
193 val metID = get_obj g_metID pt pp;
194 val {srls=srls,scr=sc,...} = get_met metID;
195 val loc as (ScrState (E,l,a,_,_,b), ctxt) = get_loc pt (p,p_);
196 val thy' = get_obj g_domID pt pp;
197 val thy = assoc_thy thy';
198 val (_,_,(scval,scsaf)) = next_tac (thy',srls) (pt,(p,p_)) sc loc;
203 val is = ScrState (E,l,a,scval,scsaf,b)
204 val tac_ = Check_Postcond' (pI, (scval, asm))
205 val (pos,ps,f,pt) = generate1 thy tac_ (is, ctxt) (pp,Res) pt;
206 in ("ok", ([(Check_Postcond pI, tac_, ((pp,Res),(is, ctxt)))], ps,(pt,pos))) end
208 let (*resume script of parpbl, transfer value of subpbl-script*)
209 val ppp = par_pblobj pt (lev_up p);
210 val thy' = get_obj g_domID pt ppp;
211 val thy = assoc_thy thy';
212 val metID = get_obj g_metID pt ppp;
213 val sc = (#scr o get_met) metID;
214 val (ScrState (E,l,a,_,_,b), ctxt') = get_loc pt (pp(*!/p/*),Frm);
215 val ctxt'' = from_subpbl_to_caller ctxt scval ctxt'
216 val ((p,p_),ps,f,pt) =
217 generate1 thy (Check_Postcond' (pI, (scval, asm)))
218 (ScrState (E,l,a,scval,scsaf,b), ctxt'') (pp,Res) pt;
219 in ("ok", ([(Check_Postcond pI, Check_Postcond'(pI,(scval, asm)),
220 ((pp,Res), (ScrState (E,l,a,scval,scsaf,b), ctxt'')))],ps,(pt,(p,p_))))
224 | solve (_,End_Proof'') (pt, (p,p_)) =
226 ([(Empty_Tac,Empty_Tac_,(([],Res),(Uistate, e_ctxt)))],[],(pt,(p,p_))))
228 (*-----------vvvvvvvvvvv could be done by generate1 ?!?*)
229 | solve (_,End_Detail' t) (pt, (p,p_)) =
231 val pr as (p',_) = (lev_up p, Res)
232 val pp = par_pblobj pt p
233 val r = (fst o (get_obj g_result pt)) p'
234 (*Rewrite_Set* done at Detail_Set*: this result is already in ptree*)
235 val thy' = get_obj g_domID pt pp
236 val (srls, is, sc) = from_pblobj' thy' pr pt
237 val (tac_,is',_) = next_tac (thy',srls) (pt,pr) sc is
238 in ("ok", ([(End_Detail, End_Detail' t ,
239 ((p,p_), get_loc pt (p,p_)))], [], (pt,pr)))
242 | solve (mI,m) (pt, po as (p,p_)) =
243 if e_metID = get_obj g_metID pt (par_pblobj pt p)(*29.8.02: could be detail, too !!*)
246 val ctxt = get_ctxt pt po
247 val ((p,p_),ps,f,pt) =
248 generate1 (assoc_thy (get_obj g_domID pt (par_pblobj pt p)))
249 m (e_istate, ctxt) (p,p_) pt;
250 in ("no-method-specified", (*Free_Solve*)
251 ([(Empty_Tac,Empty_Tac_, ((p,p_),(Uistate, ctxt)))], ps, (pt,(p,p_))))
255 val thy' = get_obj g_domID pt (par_pblobj pt p);
256 val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
257 val d = e_rls; (*FIXME.WN0108: canon.simplifier for domain is missing: generate from domID?*)
259 case locate_gen (thy',srls) m (pt,(p,p_)) (sc,d) is of
260 Steps (is', ss as (m',f',pt',p',c')::_) =>
262 (*27.8.02:next_tac may change to other branches in pt FIXXXXME*)
263 in ("ok", (map step2taci ss, c', (pt',p'))) end
265 let val (p,ps,f,pt) =
266 generate_hard (assoc_thy "Isac") m (p,p_) pt;
267 in ("not-found-in-script",
268 ([(tac_2tac m, m, (po,is))], ps, (pt,p)))
272 (*FIXME.WN050821 compare fun solve ... fun nxt_solv*)
273 (* nxt_solv (Apply_Method' vvv FIXME: get args in applicable_in *)
274 fun nxt_solv (Apply_Method' (mI,_,_,_)) _ (pt:ptree, pos as (p,_)) =
276 val {srls,ppc,...} = get_met mI;
277 val PblObj{meth=itms,origin=(oris,_,_),probl, ...} = get_obj I pt p;
278 val itms = if itms <> [] then itms else complete_metitms oris probl [] ppc
279 val thy' = get_obj g_domID pt p;
280 val thy = assoc_thy thy';
281 val (is as ScrState (env,_,_,_,_,_), ctxt, scr) = init_scrstate thy itms mI;
282 val ini = init_form thy scr env;
287 val pos = ((lev_on o lev_dn) p, Frm)
288 val tac_ = Apply_Method' (mI, SOME t, is, ctxt);
289 val (pos,c,_,pt) = (*implicit Take*)
290 generate1 thy tac_ (is, ctxt) pos pt
291 in ([(Apply_Method mI, tac_, (pos, (is, ctxt)))], c, (pt, pos)):calcstate' end
294 val pt = update_env pt (fst pos) (SOME (is, ctxt))
295 val (tacis, c, ptp) = nxt_solve_ (pt, pos)
297 [(Apply_Method mI, Apply_Method' (mI, NONE, e_istate, ctxt), (pos, (is, ctxt)))],
302 (*TODO.WN050913 remove unnecessary code below*)
303 | nxt_solv (Check_Postcond' (pI,_)) _ (pt, pos as (p,p_)) =
305 val pp = par_pblobj pt p
307 (case get_obj g_tac pt p of
308 Check_elementwise _ => (*collects and instantiates asms*)
309 (snd o (get_obj g_result pt)) p
310 | _ => get_assumptions_ pt (p,p_))
311 handle _ => [] (*FIXME.WN030527 asms in subpbls not completely clear*)
312 val metID = get_obj g_metID pt pp;
313 val {srls=srls,scr=sc,...} = get_met metID;
314 val loc as (ScrState (E,l,a,_,_,b), ctxt) = get_loc pt (p,p_);
315 val thy' = get_obj g_domID pt pp;
316 val thy = assoc_thy thy';
317 val (_,_,(scval,scsaf)) = next_tac (thy',srls) (pt,(p,p_)) sc loc;
322 val is = ScrState (E,l,a,scval,scsaf,b)
323 val tac_ = Check_Postcond'(pI,(scval, asm))
324 val ((p,p_),ps,f,pt) =
325 generate1 thy tac_ (is, ctxt) (pp,Res) pt;
326 in ([(Check_Postcond pI, tac_, ((pp,Res), (is, ctxt)))],ps,(pt, (p,p_))) end
328 let (*resume script of parpbl, transfer value of subpbl-script*)
329 val ppp = par_pblobj pt (lev_up p);
330 val thy' = get_obj g_domID pt ppp;
331 val thy = assoc_thy thy';
332 val metID = get_obj g_metID pt ppp;
333 val {scr,...} = get_met metID;
334 val (ScrState (E,l,a,_,_,b), ctxt') = get_loc pt (pp(*!/p/*),Frm)
335 val ctxt'' = from_subpbl_to_caller ctxt scval ctxt'
336 val tac_ = Check_Postcond' (pI, (scval, asm))
337 val is = ScrState (E,l,a,scval,scsaf,b)
338 val ((p,p_),ps,f,pt) = generate1 thy tac_ (is, ctxt'') (pp, Res) pt;
339 in ([(Check_Postcond pI, tac_, ((pp, Res), (is, ctxt'')))], ps, (pt, (p,p_))) end
342 | nxt_solv (End_Proof'') _ ptp = ([], [], ptp)
344 | nxt_solv tac_ is (pt, pos as (p,p_)) =
348 (p, Met) => ((lev_on o lev_dn) p, Frm)(*begin script*)
349 | (p, Res) => (lev_on p,Res) (*somewhere in script*)
351 val (pos',c,_,pt) = generate1 (assoc_thy "Isac") tac_ is pos pt;
352 in ([(tac_2tac tac_, tac_, (pos,is))], c, (pt, pos')) end
354 (* find the next tac from the script, nxt_solv will update the ptree *)
355 and nxt_solve_ (ptp as (pt, pos as (p,p_))) =
356 if e_metID = get_obj g_metID pt (par_pblobj pt p)
357 then ([], [], (pt,(p,p_))):calcstate'
360 val thy' = get_obj g_domID pt (par_pblobj pt p);
361 val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
362 val (tac_,is,(t,_)) = next_tac (thy',srls) (pt,pos) sc is;
363 (*TODO here ^^^ return finished/helpless/ok !*)
365 End_Detail' _ => ([(End_Detail, End_Detail' (t,[(*FIXME.040215*)]),
366 (pos, is))], [], (pt, pos))
367 | _ => nxt_solv tac_ is ptp end;
369 (* says how may steps of a calculation should be done by "fun autocalc" *)
370 (*TODO.WN0512 redesign togehter with autocalc ?*)
372 Step of int (*1 do #int steps; may stop in model/specify:
373 IS VERY INEFFICIENT IN MODEL/SPECIY*)
374 | CompleteModel (*2 complete modeling
375 if model complete, finish specifying + start solving*)
376 | CompleteCalcHead (*3 complete model/specify in one go + start solving*)
377 | CompleteToSubpbl (*4 stop at the next begin of a subproblem,
378 if none, complete the actual (sub)problem*)
379 | CompleteSubpbl (*5 complete the actual (sub)problem (incl.ev.subproblems)*)
380 | CompleteCalc; (*6 complete the calculation as a whole*)
382 fun autoord (Step _ ) = 1
383 | autoord CompleteModel = 2
384 | autoord CompleteCalcHead = 3
385 | autoord CompleteToSubpbl = 4
386 | autoord CompleteSubpbl = 5
387 | autoord CompleteCalc = 6;
389 fun complete_solve auto c (ptp as (_, p as (_,p_)): ptree * pos') =
391 then ("end-of-calculation", [], ptp)
393 if member op = [Pbl,Met] p_
396 val ptp = all_modspec ptp
397 val (_, c', ptp) = all_solve auto c ptp
398 in complete_solve auto (c @ c') ptp end
400 case nxt_solve_ ptp of
401 ((Subproblem _, tac_, (_, is))::_, c', ptp') =>
403 then ("ok", c @ c', ptp)
406 val ptp = all_modspec ptp'
407 val (_, c'', ptp) = all_solve auto (c @ c') ptp
408 in complete_solve auto (c @ c'@ c'') ptp end
409 | ((Check_Postcond _, tac_, (_, is))::_, c', ptp' as (_, p')) =>
410 if autoord auto < 6 orelse p' = ([],Res)
411 then ("ok", c @ c', ptp')
412 else complete_solve auto (c @ c') ptp'
413 | ((End_Detail, _, _)::_, c', ptp') =>
415 then ("ok", c @ c', ptp')
416 else complete_solve auto (c @ c') ptp'
417 | (_, c', ptp') => complete_solve auto (c @ c') ptp'
419 and all_solve auto c (ptp as (pt, pos as (p,_)): ptree * pos') =
421 val (_,_,mI) = get_obj g_spec pt p
422 val ctxt = get_ctxt pt pos
423 val (_, c', ptp) = nxt_solv (Apply_Method' (mI, NONE, e_istate, ctxt)) (e_istate, ctxt) ptp
424 in complete_solve auto (c @ c') ptp end;
426 (* aux.fun for detailrls with Rrls, reverse rewriting *)
427 fun rul_terms_2nds _ nds t [] = nds
428 | rul_terms_2nds thy nds t ((rule, res as (t', _)) :: rts) =
429 (append_atomic [] (e_istate, e_ctxt) t (rule2tac thy [] rule) res Complete EmptyPtree) ::
430 (rul_terms_2nds thy nds t' rts);
432 (* detail steps done internally by Rewrite_Set* into ctree by use of a script *)
433 fun detailrls pt (pos as (p,p_):pos') =
435 val t = get_obj g_form pt p
436 val tac = get_obj g_tac pt p
437 val rls = (assoc_rls o rls_of) tac
438 val ctxt = get_ctxt pt pos
441 Rrls {scr = Rfuns {init_state,...},...} =>
443 val (_,_,_,rul_terms) = init_state t
444 val newnds = rul_terms_2nds (Proof_Context.theory_of ctxt) [] t rul_terms
445 val pt''' = ins_chn newnds pt p
446 in ("detailrls", pt''', (p @ [length newnds], Res):pos') end
449 val is = init_istate tac t
450 (*TODO.WN060602 ScrState (["(t_, Problem (Isac,[equation,univar]))"]
451 is wrong for simpl, but working ?!? *)
452 val tac_ = Apply_Method' (e_metID(*WN0402: see generate1 !?!*), SOME t, is, ctxt)
453 val pos' = ((lev_on o lev_dn) p, Frm)
454 val thy = assoc_thy "Isac"
455 val (_,_,_,pt') = (*implicit Take*)generate1 thy tac_ (is, ctxt) pos' pt
456 val (_,_,(pt'',_)) = complete_solve CompleteSubpbl [] (pt',pos')
457 val newnds = children (get_nd pt'' p)
458 val pt''' = ins_chn newnds pt p
459 (*complete_solve cuts branches after*)
460 in ("detailrls", pt''', (p @ [length newnds], Res):pos') end
465 (* val(mI,m)=m;val ppp=p;(*!!!*)val(p,p_)=pos;val(_,pt,_)=ppp(*!!!*);
466 get_form ((mI,m):tac'_) ((p,p_):pos') ppp;
468 fun get_form ((mI,m):tac'_) ((p,p_):pos') pt =
469 case applicable_in (p,p_) pt m of
470 Notappl e => Error' (Error_ e)
472 (* val Appl m=applicable_in (p,p_) pt m;
474 if member op = specsteps mI
475 then let val (_,_,f,_,_,_) = specify m (p,p_) [] pt
477 else let val (*_,_,f,_,_,_*)_ = solve (mI,m) (pt,(p,p_))
478 in (*f*) EmptyMout end;