src/Tools/isac/Interpret/solve.sml
author Walther Neuper <neuper@ist.tugraz.at>
Thu, 23 Sep 2010 14:49:23 +0200
branchisac-update-Isa09-2
changeset 38014 3e11e3c2dc42
parent 37980 c0a9d6bdc1d6
child 38015 67ba02dffacc
permissions -rw-r--r--
updated "op +", "op -", "op *". "HOL.divide" in src & test

find . -type f -exec sed -i s/"\"op +\""/"\"Groups.plus_class.plus\""/g {} \;
find . -type f -exec sed -i s/"\"op -\""/"\"Groups.minus_class.minus\""/g {} \;
find . -type f -exec sed -i s/"\"op *\""/"\"Groups.times_class.times\""/g {} \;
find . -type f -exec sed -i s/"\"HOL.divide\""/"\"Rings.inverse_class.divide\""/g {} \;
     1 (* solve an example by interpreting a method's script
     2    (c) Walther Neuper 1999
     3 
     4 use"ME/solve.sml";
     5 use"solve.sml";
     6 *)
     7 
     8 fun safe (ScrState (_,_,_,_,s,_)) = s
     9   | safe (RrlsState _) = Safe;
    10 
    11 type mstID = string;
    12 type tac'_ = mstID * tac; (*DG <-> ME*)
    13 val e_tac'_ = ("Empty_Tac", Empty_Tac):tac'_;
    14 
    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') 
    26 
    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)
    33 		    
    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)
    41 
    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'))
    45 
    46 | Calculate op_             => ("Calculate", Calculate op_)
    47 | Substitute sube           => ("Substitute", Substitute sube) 
    48 | Apply_Assumption cts'	    => ("Apply_Assumption", Apply_Assumption cts')
    49 
    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)) 
    54 (*
    55 | Subproblem_Full(spec,cts')=> ("Subproblem_Full", Subproblem_Full(spec,cts')) 
    56 *)
    57 | End_Subproblem            => ("End_Subproblem",End_Subproblem)
    58 | CAScmd cterm'		    => ("CAScmd", CAScmd cterm')
    59 			    
    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) 
    73 			    
    74 | Empty_Tac               => ("Empty_Tac",Empty_Tac)
    75 | Tac string              => ("Tac",Tac string)
    76 | User                      => ("User",User)
    77 | End_Proof'                => ("End_Proof'",End_Proof'); 
    78 
    79 (*Detail*)
    80 val empty_tac'_ = (mk_tac'_ Empty_Tac):tac'_;
    81 
    82 fun mk_tac ((_,m):tac'_) = m; 
    83 fun mk_mstID ((mI,_):tac'_) = mI;
    84 
    85 fun tac'_2str ((ID,ms):tac'_) = ID ^ (tac2str ms);
    86 (* TODO: tac2str, tac'_2str NOT tested *)
    87 
    88 
    89 
    90 type squ = ptree; (* TODO: safe etc. *)
    91 
    92 (*13.9.02--------------
    93 type ctr = (loc * pos) list;
    94 val ops = [("PLUS","Groups.plus_class.plus"),("MINUS","Groups.minus_class.minus"),("TIMES","op *"),
    95 	   ("cancel","cancel"),("pow","pow"),("sqrt","sqrt")];
    96 ML {* 
    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")*)
   109 *}
   110 fun op_intern op_ =
   111   case assoc (ops, op_) of
   112     SOME op' => op' | NONE => raise error ("op_intern: no op= "^op_);
   113 -----------------------*)
   114 
   115 
   116 
   117 (* use"ME/solve.sml";
   118    use"solve.sml";
   119 
   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";
   122 
   123   Const ("Script.Rewrite'_Inst",_) $ sub $ Free (thm',_) $ Const (pa,_) $ f'
   124    *)
   125 
   126 
   127 
   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"];
   133 
   134 "-----------------------------------------------------------------------";
   135 
   136 
   137 fun step2taci ((tac_, _, pt, p, _):step) = (*FIXXME.040312: redesign step*)
   138     (tac_2tac tac_, tac_, (p, get_istate pt p)):taci;
   139 
   140 
   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);
   144    *)
   145 fun solve ("Apply_Method", m as Apply_Method' (mI, _, _)) 
   146 	  (pt:ptree, (pos as (p,_))) =
   147   let val {srls,...} = get_met mI;
   148     val PblObj{meth=itms,...} = get_obj I pt p;
   149     val thy' = get_obj g_domID pt p;
   150     val thy = assoc_thy thy';
   151     val (is as ScrState (env,_,_,_,_,_), sc) = init_scrstate thy itms mI;
   152     val ini = init_form thy sc env;
   153     val p = lev_dn p;
   154   in 
   155       case ini of
   156 	  SOME t => (* val SOME t = ini; 
   157 	             *)
   158 	  let val (pos,c,_,pt) = 
   159 		  generate1 thy (Apply_Method' (mI, SOME t, is))
   160 			    is (lev_on p, Frm)(*implicit Take*) pt;
   161 	  in ("ok",([(Apply_Method mI, Apply_Method' (mI, SOME t, is), 
   162 		      ((lev_on p, Frm), is))], c, (pt,pos)):calcstate') 
   163 	  end	      
   164 	| NONE => (*execute the first tac in the Script, compare solve m*)
   165 	  let val (m', is', _) = next_tac (thy', srls) (pt, (p, Res)) sc is;
   166 	      val d = e_rls (*FIXME: get simplifier from domID*);
   167 	  in 
   168 	      case locate_gen (thy',srls) m' (pt,(p, Res))(sc,d) is' of 
   169 		  Steps (is'', ss as (m'',f',pt',p',c')::_) =>
   170 (* val Steps (is'', ss as (m'',f',pt',p',c')::_) =
   171        locate_gen (thy',srls) m'  (pt,(p,Res)) (sc,d) is';
   172  *)
   173 		  ("ok", (map step2taci ss, c', (pt',p')))
   174 		| NotLocatable =>  
   175 		  let val (p,ps,f,pt) = 
   176 			  generate_hard (assoc_thy "Isac.thy") m (p,Frm) pt;
   177 		  in ("not-found-in-script",
   178 		      ([(tac_2tac m, m, (pos, is))], ps, (pt,p))) end
   179     (*just-before------------------------------------------------------
   180 	      ("ok",([(Apply_Method mI,Apply_Method'(mI,NONE,e_istate),
   181 		       (pos, is))],
   182 		     [], (update_env pt (fst pos) (SOME is),pos)))
   183      -----------------------------------------------------------------*)
   184 	  end
   185   end
   186 
   187   | solve ("Free_Solve", Free_Solve')  (pt,po as (p,_)) =
   188   let (*val _=writeln"###solve Free_Solve";*)
   189     val p' = lev_dn_ (p,Res);
   190     val pt = update_metID pt (par_pblobj pt p) e_metID;
   191   in ("ok", ((*(p',Uistate,Empty_Tac_), EmptyMout, Empty_Tac, Unsafe,*)
   192       [(Empty_Tac, Empty_Tac_, (po, Uistate))], [], (pt,p'))) end
   193 
   194 (* val (("Check_Postcond",Check_Postcond' (pI,_)), (pt,(pos as (p,p_)))) =
   195        (  m,                                       (pt, pos));
   196    *)
   197   | solve ("Check_Postcond",Check_Postcond' (pI,_)) (pt,(pos as (p,p_))) =
   198     let (*val _=writeln"###solve Check_Postcond";*)
   199       val pp = par_pblobj pt p
   200       val asm = (case get_obj g_tac pt p of
   201 		    Check_elementwise _ => (*collects and instantiates asms*)
   202 		    (snd o (get_obj g_result pt)) p
   203 		  | _ => ((map fst) o (get_assumptions_ pt)) (p,p_))
   204 	  handle _ => [] (*WN.27.5.03 asms in subpbls not completely clear*)
   205       val metID = get_obj g_metID pt pp;
   206       val {srls=srls,scr=sc,...} = get_met metID;
   207       val is as ScrState (E,l,a,_,_,b) = get_istate pt (p,p_); 
   208      (*val _= writeln("### solve Check_postc, subpbl pos= "^(pos'2str (p,p_)));
   209       val _= writeln("### solve Check_postc, is= "^(istate2str is));*)
   210       val thy' = get_obj g_domID pt pp;
   211       val thy = assoc_thy thy';
   212       val (_,_,(scval,scsaf)) = next_tac (thy',srls) (pt,(p,p_)) sc is;
   213       (*val _= writeln("### solve Check_postc, scval= "^(term2str scval));*)
   214 
   215     in if pp = [] then
   216 	   let val is = ScrState (E,l,a,scval,scsaf,b)
   217 	       val tac_ = Check_Postcond'(pI,(scval, map term2str asm))
   218 	       val (pos,ps,f,pt) = generate1 thy tac_ is (pp,Res) pt;
   219 	   in ("ok", ((*(([],Res),is,End_Proof''), f, End_Proof', scsaf,*)
   220 	       [(Check_Postcond pI, tac_, ((pp,Res),is))], ps,(pt,pos))) end
   221        else
   222         let
   223 	  (*resume script of parpbl, transfer value of subpbl-script*)
   224         val ppp = par_pblobj pt (lev_up p);
   225 	val thy' = get_obj g_domID pt ppp;
   226         val thy = assoc_thy thy';
   227 	val metID = get_obj g_metID pt ppp;
   228         val sc = (#scr o get_met) metID;
   229         val is as ScrState (E,l,a,_,_,b) = get_istate pt (pp(*!/p/*),Frm); 
   230      (*val _=writeln("### solve Check_postc, parpbl pos= "^(pos'2str(pp,Frm)));
   231   	val _=writeln("### solve Check_postc, is(pt)= "^(istate2str is));
   232   	val _=writeln("### solve Check_postc, is'= "^
   233 		      (istate2str (E,l,a,scval,scsaf,b)));*)
   234         val ((p,p_),ps,f,pt) = 
   235 	    generate1 thy (Check_Postcond' (pI, (scval, map term2str asm)))
   236 		(ScrState (E,l,a,scval,scsaf,b)) (pp,Res) pt;
   237 	(*val _=writeln("### solve Check_postc, is(pt')= "^
   238 		      (istate2str (get_istate pt ([3],Res))));
   239 	val (nx,is',_) = next_tac (thy',srls) (pt,(p,p_)) sc 
   240 				(ScrState (E,l,a,scval,scsaf,b));*)
   241        in ("ok",(*((pp,Res),is',nx), f, tac_2tac nx, scsaf,*)
   242 	   ([(Check_Postcond pI, Check_Postcond'(pI,(scval, map term2str asm)),
   243 	      ((pp,Res), ScrState (E,l,a,scval,scsaf,b)))],ps,(pt,(p,p_))))
   244 	end
   245     end
   246 (* val (msg, cs') = 
   247     ("ok",([(Check_Postcond pI,Check_Postcond'(pI, (scval, map term2str asm))),
   248 	    ((pp,Res),(ScrState (E,l,a,scval,scsaf,b)))], (pt,(p,p_))));
   249     val (_,(pt',p')) = cs';
   250    (writeln o istate2str) (get_istate pt' p');
   251    (term2str o fst) (get_obj g_result pt' (fst p'));
   252    *)
   253 
   254 (* writeln(istate2str(get_istate pt (p,p_)));
   255    *)
   256   | solve (_,End_Proof'') (pt, (p,p_)) =
   257       ("end-proof",
   258        ((*(([],Res),Uistate,Empty_Tac_), EmptyMout, Empty_Tac, Safe,*)
   259        [(Empty_Tac,Empty_Tac_,(([],Res),Uistate))],[],(pt,(p,p_))))
   260 
   261 (*-----------vvvvvvvvvvv could be done by generate1 ?!?*)
   262   | solve (_,End_Detail' t) (pt,(p,p_)) =
   263     let val pr as (p',_) = (lev_up p, Res)
   264 	val pp = par_pblobj pt p
   265 	val r = (fst o (get_obj g_result pt)) p' 
   266 	(*Rewrite_Set* done at Detail_Set*: this result is already in ptree*)
   267 	val thy' = get_obj g_domID pt pp
   268 	val (srls, is, sc) = from_pblobj' thy' pr pt
   269 	val (tac_,is',_) = next_tac (thy',srls)  (pt,pr) sc is
   270     in ("ok", ((*((pp,Frm(*???*)),is,tac_), 
   271 	Form' (FormKF (~1, EdUndef, length p', Nundef, term2str r)),
   272 	tac_2tac tac_, Sundef,*)
   273 	[(End_Detail, End_Detail' t , 
   274 	  ((p,p_), get_istate pt (p,p_)))], [], (pt,pr))) end
   275 
   276   | solve (mI,m) (pt, po as (p,p_)) =
   277 (* val ((mI,m), (pt, po as (p,p_))) = (m, (pt, pos));
   278    *)
   279     if e_metID = get_obj g_metID pt (par_pblobj pt p)(*29.8.02:
   280 						      could be detail, too !!*)
   281     then let val ((p,p_),ps,f,pt) = 
   282 		 generate1 (assoc_thy (get_obj g_domID pt (par_pblobj pt p))) 
   283 			   m e_istate (p,p_) pt;
   284 	 in ("no-method-specified", (*Free_Solve*)
   285 	     ((*((p,p_),Uistate,Empty_Tac_),f, Empty_Tac, Unsafe,*)
   286 	     [(Empty_Tac,Empty_Tac_, ((p,p_),Uistate))], ps, (pt,(p,p_)))) end
   287     else
   288 	let 
   289 	    val thy' = get_obj g_domID pt (par_pblobj pt p);
   290 	    val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
   291 (*val _= writeln("### solve, before locate_gen p="^(pos'2str(p,p_)));*)
   292 		val d = e_rls; (*FIXME: canon.simplifier for domain is missing
   293 				8.01: generate from domID?*)
   294 	in case locate_gen (thy',srls) m  (pt,(p,p_)) (sc,d) is of 
   295 	       Steps (is', ss as (m',f',pt',p',c')::_) =>
   296 (* val Steps (is', ss as (m',f',pt',p',c')::_) =
   297        locate_gen (thy',srls) m  (pt,(p,p_)) (sc,d) is;
   298  *)
   299 	       let (*val _= writeln("### solve, after locate_gen: is= ")
   300 		       val _= writeln(istate2str is')*)
   301 		   (*val nxt_ = 
   302 		       case p' of (*change from solve to model subpbl*)
   303 			   (_,Pbl) => nxt_model_pbl m' (pt',p')
   304 			 | _ => fst3 (next_tac (thy',srls) (pt',p') sc is');*) 
   305 	       (*27.8.02:next_tac may change to other branches in pt FIXXXXME*)
   306 	       in ("ok", ((*(p',is',nxt_), f', tac_2tac nxt_, safe is',*)
   307 		   map step2taci ss, c', (pt',p'))) end
   308 	     | NotLocatable =>  
   309 	       let val (p,ps,f,pt) = 
   310 		       generate_hard (assoc_thy "Isac.thy") m (p,p_) pt;
   311 	       in ("not-found-in-script",
   312 		   ((*(p,Uistate,Empty_Tac_),f, Empty_Tac, Unsafe,*) 
   313 		   [(tac_2tac m, m, (po,is))], ps, (pt,p))) end
   314 	end;
   315 
   316 
   317 (*FIXME.WN050821 compare solve ... nxt_solv*)
   318 (* nxt_solv (Apply_Method'     vvv FIXME: get args in applicable_in *)
   319 fun nxt_solv (Apply_Method' (mI,_,_)) _ (pt:ptree, pos as (p,_)) =
   320 (* val ((Apply_Method' (mI,_,_)),             _,    (pt:ptree, pos as (p,_))) =
   321        ((Apply_Method' (mI, NONE, e_istate)), e_istate, ptp);
   322    *)
   323   let val {srls,ppc,...} = get_met mI;
   324     val PblObj{meth=itms,origin=(oris,_,_),probl,...} = get_obj I pt p;
   325     val itms = if itms <> [] then itms
   326 	       else complete_metitms oris probl [] ppc
   327     val thy' = get_obj g_domID pt p;
   328     val thy = assoc_thy thy';
   329     val (is as ScrState (env,_,_,_,_,_), scr) = init_scrstate thy itms mI;
   330     val ini = init_form thy scr env;
   331   in 
   332     case ini of
   333     SOME t => (* val SOME t = ini; 
   334 	         *)
   335     let val pos = ((lev_on o lev_dn) p, Frm)
   336 	val tac_ = Apply_Method' (mI, SOME t, is);
   337 	val (pos,c,_,pt) = (*implicit Take*)
   338 	    generate1 thy tac_ is pos pt
   339       (*val _= ("### nxt_solv Apply_Method, pos= "^pos'2str (lev_on p,Frm));*)
   340     in ([(Apply_Method mI, tac_, (pos, is))], c, (pt, pos)):calcstate' end
   341   | NONE =>
   342     let val pt = update_env pt (fst pos) (SOME is)
   343 	val (tacis, c, ptp) = nxt_solve_ (pt, pos)
   344     in (tacis @ 
   345 	[(Apply_Method mI, Apply_Method' (mI, NONE, e_istate), (pos, is))],
   346 	c, ptp) end
   347   end
   348 (* val ("Check_Postcond",Check_Postcond' (pI,_)) = (mI,m);
   349    val (Check_Postcond' (pI,_), _, (pt, pos as (p,p_))) = 
   350        (tac_,                  is,  ptp);
   351    *)
   352   (*TODO.WN050913 remove unnecessary code below*)
   353   | nxt_solv (Check_Postcond' (pI,_)) _ (pt, pos as (p,p_))  =
   354     let (*val _=writeln"###solve Check_Postcond";*)
   355       val pp = par_pblobj pt p
   356       val asm = (case get_obj g_tac pt p of
   357 		    Check_elementwise _ => (*collects and instantiates asms*)
   358 		    (snd o (get_obj g_result pt)) p
   359 		  | _ => ((map fst) o (get_assumptions_ pt)) (p,p_))
   360 	  handle _ => [] (*WN.27.5.03 asms in subpbls not completely clear*)
   361       val metID = get_obj g_metID pt pp;
   362       val {srls=srls,scr=sc,...} = get_met metID;
   363       val is as ScrState (E,l,a,_,_,b) = get_istate pt (p,p_); 
   364      (*val _= writeln("### solve Check_postc, subpbl pos= "^(pos'2str (p,p_)));
   365       val _= writeln("### solve Check_postc, is= "^(istate2str is));*)
   366       val thy' = get_obj g_domID pt pp;
   367       val thy = assoc_thy thy';
   368       val (_,_,(scval,scsaf)) = next_tac (thy',srls) (pt,(p,p_)) sc is;
   369       (*val _= writeln("### solve Check_postc, scval= "^(term2str scval));*)
   370     in if pp = [] then 
   371 	   let val is = ScrState (E,l,a,scval,scsaf,b)
   372 	       val tac_ = Check_Postcond'(pI,(scval, map term2str asm))
   373            (*val _= writeln"### nxt_solv2 Apply_Method: stored is =";
   374                val _= writeln(istate2str is);*)
   375 	       val ((p,p_),ps,f,pt) = 
   376 		   generate1 thy tac_ is (pp,Res) pt;
   377 	   in ([(Check_Postcond pI, tac_, ((pp,Res), is))],ps,(pt, (p,p_))) end
   378        else
   379         let
   380 	  (*resume script of parpbl, transfer value of subpbl-script*)
   381         val ppp = par_pblobj pt (lev_up p);
   382 	val thy' = get_obj g_domID pt ppp;
   383         val thy = assoc_thy thy';
   384 	val metID = get_obj g_metID pt ppp;
   385 	val {scr,...} = get_met metID;
   386         val is as ScrState (E,l,a,_,_,b) = get_istate pt (pp(*!/p/*),Frm)
   387         val tac_ = Check_Postcond' (pI, (scval, map term2str asm))
   388 	val is = ScrState (E,l,a,scval,scsaf,b)
   389     (*val _= writeln"### nxt_solv3 Apply_Method: stored is =";
   390         val _= writeln(istate2str is);*)
   391         val ((p,p_),ps,f,pt) = generate1 thy tac_ is (pp, Res) pt;
   392 	(*val (nx,is',_) = next_tac (thy',srls) (pt,(p,p_)) scr is;WN050913*)
   393        in ([(Check_Postcond pI, tac_, ((pp, Res), is))], ps, (pt, (p,p_))) end
   394     end
   395 (* writeln(istate2str(get_istate pt (p,p_)));
   396    *)
   397 
   398 (*.start interpreter and do one rewrite.*)
   399 (* val (_,Detail_Set'(thy',rls,t)) = (mI,m); val p = (p,p_);
   400    solve ("",Detail_Set'(thy', rls, t)) p pt;
   401   | nxt_solv (Detail_Set'(thy', rls, t)) _ (pt, p) = **********
   402 ---> Frontend/sml.sml
   403 
   404   | nxt_solv (End_Detail' t) _ (pt, (p,p_)) = **********
   405     let val pr as (p',_) = (lev_up p, Res)
   406 	val pp = par_pblobj pt p
   407 	val r = (fst o (get_obj g_result pt)) p' 
   408 	(*Rewrite_Set* done at Detail_Set*: this result is already in ptree*)
   409 	val thy' = get_obj g_domID pt pp
   410 	val (srls, is, sc) = from_pblobj' thy' pr pt
   411 	val (tac_,is',_) = next_tac (thy',srls)  (pt,pr) sc is
   412     in (pr, ((pp,Frm(*???*)),is,tac_), 
   413 	Form' (FormKF (~1, EdUndef, length p', Nundef, term2str r)),
   414 	tac_2tac tac_, Sundef, pt) end
   415 *)
   416   | nxt_solv (End_Proof'') _ ptp = ([], [], ptp)
   417 
   418   | nxt_solv tac_ is (pt, pos as (p,p_)) =
   419 (* val (pt, pos as (p,p_)) = ptp;
   420    *)
   421     let val pos = case pos of
   422 		      (p, Met) => ((lev_on o lev_dn) p, Frm)(*begin script*)
   423 		    | (p, Res) => (lev_on p,Res) (*somewhere in script*)
   424 		    | _ => pos  (*somewhere in script*)
   425     (*val _= writeln"### nxt_solv4 Apply_Method: stored is =";
   426         val _= writeln(istate2str is);*)
   427 	val (pos',c,_,pt) = generate1 (assoc_thy "Isac.thy") tac_ is pos pt;
   428     in ([(tac_2tac tac_, tac_, (pos,is))], c, (pt, pos')) end
   429 
   430 
   431   (*(p,p_), (([],Res),Uistate,Empty_Tac_), EmptyMout, Empty_Tac, Safe, pt*)
   432 
   433 
   434 (*.find the next tac from the script, nxt_solv will update the ptree.*)
   435 (* val (ptp as (pt,pos as (p,p_))) = ptp';
   436    val (ptp as (pt, pos as (p,p_))) = ptp'';
   437    val (ptp as (pt, pos as (p,p_))) = ptp;
   438    val (ptp as (pt, pos as (p,p_))) = (pt,ip);
   439    val (ptp as (pt, pos as (p,p_))) = (pt, pos);
   440    *)
   441 and nxt_solve_ (ptp as (pt, pos as (p,p_))) =
   442     if e_metID = get_obj g_metID pt (par_pblobj pt p)
   443     then ([], [], (pt,(p,p_))):calcstate'
   444     else let val thy' = get_obj g_domID pt (par_pblobj pt p);
   445 	     val (srls, is, sc) = from_pblobj_or_detail' thy' (p,p_) pt;
   446 	     val (tac_,is,(t,_)) = next_tac (thy',srls) (pt,pos) sc is;
   447 	 (*TODO here ^^^  return finished/helpless/ok !*)
   448 	 (* val (tac_',is',(t',_)) = next_tac (thy',srls) (pt,pos) sc is;
   449 	    *)
   450 	 in case tac_ of
   451 		End_Detail' _ => ([(End_Detail, 
   452 				    End_Detail' (t,[(*FIXME.040215*)]), 
   453 				    (pos, is))], [], (pt, pos))
   454 	      | _ => nxt_solv tac_ is ptp end;
   455 
   456 (*.says how may steps of a calculation should be done by "fun autocalc".*)
   457 (*TODO.WN0512 redesign togehter with autocalc ?*)
   458 datatype auto = 
   459   Step of int      (*1 do #int steps; may stop in model/specify:
   460 		     IS VERY INEFFICIENT IN MODEL/SPECIY*)
   461 | CompleteModel    (*2 complete modeling
   462                      if model complete, finish specifying + start solving*)
   463 | CompleteCalcHead (*3 complete model/specify in one go + start solving*)
   464 | CompleteToSubpbl (*4 stop at the next begin of a subproblem,
   465                      if none, complete the actual (sub)problem*)
   466 | CompleteSubpbl   (*5 complete the actual (sub)problem (incl.ev.subproblems)*)
   467 | CompleteCalc;    (*6 complete the calculation as a whole*)	
   468 fun autoord (Step _ ) = 1
   469   | autoord CompleteModel = 2
   470   | autoord CompleteCalcHead = 3
   471   | autoord CompleteToSubpbl = 4
   472   | autoord CompleteSubpbl = 5
   473   | autoord CompleteCalc = 6;
   474 
   475 (* val (auto, c, (ptp as (_, p))) = (auto, (c@c'), ptp);
   476    *)
   477 fun complete_solve auto c (ptp as (_, p): ptree * pos') =
   478     if p = ([], Res) then ("end-of-calculation", [], ptp) else
   479     case nxt_solve_ ptp of
   480 	((Subproblem _, tac_, (_, is))::_, c', ptp') =>
   481 (* val ptp' = ptp''';
   482    *)
   483 	if autoord auto < 5 then ("ok", c@c', ptp)
   484 	else let val ptp = all_modspec ptp';
   485 	         val (_, c'', ptp) = all_solve auto (c@c') ptp;
   486 	     in complete_solve auto (c@c'@c'') ptp end
   487       | ((Check_Postcond _, tac_, (_, is))::_, c', ptp' as (_, p')) =>
   488 	if autoord auto < 6 orelse p' = ([],Res) then ("ok", c@c', ptp')
   489 	else complete_solve auto (c@c') ptp'
   490       | ((End_Detail, _, _)::_, c', ptp') => 
   491 	if autoord auto < 6 then ("ok", c@c', ptp')
   492 	else complete_solve auto (c@c') ptp'
   493       | (_, c', ptp') => complete_solve auto (c@c') ptp'
   494 (* val (tacis, c', ptp') = nxt_solve_ ptp;
   495    val (tacis, c', ptp'') = nxt_solve_ ptp';
   496    val (tacis, c', ptp''') = nxt_solve_ ptp'';
   497    val (tacis, c', ptp'''') = nxt_solve_ ptp''';
   498    val (tacis, c', ptp''''') = nxt_solve_ ptp'''';
   499    *)
   500 and all_solve auto c (ptp as (pt, (p,_)): ptree * pos') = 
   501 (* val (ptp as (pt, (p,_))) = ptp;
   502    val (ptp as (pt, (p,_))) = ptp';
   503    val (ptp as (pt, (p,_))) = (pt, pos);
   504    *)
   505     let val (_,_,mI) = get_obj g_spec pt p;
   506         val (_, c', ptp) = nxt_solv (Apply_Method' (mI, NONE, e_istate))
   507 				e_istate ptp;
   508     in complete_solve auto (c@c') ptp end;
   509 (*@@@ vvv @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@*)
   510 fun complete_solve auto c (ptp as (_, p as (_,p_)): ptree * pos') =
   511     if p = ([], Res) then ("end-of-calculation", [], ptp) else
   512     if member op = [Pbl,Met] p_
   513     then let val ptp = all_modspec ptp
   514 	     val (_, c', ptp) = all_solve auto c ptp
   515 	 in complete_solve auto (c@c') ptp end
   516     else case nxt_solve_ ptp of
   517 	     ((Subproblem _, tac_, (_, is))::_, c', ptp') =>
   518 	     if autoord auto < 5 then ("ok", c@c', ptp)
   519 	     else let val ptp = all_modspec ptp'
   520 		      val (_, c'', ptp) = all_solve auto (c@c') ptp
   521 		  in complete_solve auto (c@c'@c'') ptp end
   522 	   | ((Check_Postcond _, tac_, (_, is))::_, c', ptp' as (_, p')) =>
   523 	     if autoord auto < 6 orelse p' = ([],Res) then ("ok", c@c', ptp')
   524 	     else complete_solve auto (c@c') ptp'
   525 	   | ((End_Detail, _, _)::_, c', ptp') => 
   526 	     if autoord auto < 6 then ("ok", c@c', ptp')
   527 	     else complete_solve auto (c@c') ptp'
   528 	   | (_, c', ptp') => complete_solve auto (c@c') ptp'
   529 and all_solve auto c (ptp as (pt, (p,_)): ptree * pos') = 
   530     let val (_,_,mI) = get_obj g_spec pt p
   531         val (_, c', ptp) = nxt_solv (Apply_Method' (mI, NONE, e_istate))
   532 				    e_istate ptp
   533     in complete_solve auto (c@c') ptp end;
   534 
   535 (*.aux.fun for detailrls with Rrls, reverse rewriting.*)
   536 (* val (nds, t, ((rule, (t', asm)) :: rts)) = ([], t, rul_terms);
   537    *)
   538 fun rul_terms_2nds nds t [] = nds
   539   | rul_terms_2nds nds t ((rule, res as (t', _)) :: rts) =
   540     (append_atomic [] e_istate t (rule2tac [] rule) res Complete EmptyPtree) ::
   541     (rul_terms_2nds nds t' rts);
   542 
   543 
   544 (*. detail steps done internally by Rewrite_Set* 
   545     into ctree by use of a script .*)
   546 (* val (pt, (p,p_)) = (pt, pos);
   547    *)
   548 fun detailrls pt ((p,p_):pos') = 
   549     let val t = get_obj g_form pt p
   550 	val tac = get_obj g_tac pt p
   551 	val rls = (assoc_rls o rls_of) tac
   552     in case rls of
   553 (* val Rrls {scr = Rfuns {init_state,...},...} = rls;
   554    *)
   555 	   Rrls {scr = Rfuns {init_state,...},...} => 
   556 	   let val (_,_,_,rul_terms) = init_state t
   557 	       val newnds = rul_terms_2nds [] t rul_terms
   558 	       val pt''' = ins_chn newnds pt p 
   559 	   in ("detailrls", pt''', (p @ [length newnds], Res):pos') end
   560 	 | _ =>
   561 	   let val is = init_istate tac t
   562 	(*TODO.WN060602 ScrState (["(t_, Problem (Isac,[equation,univar]))"]
   563 				      is wrong for simpl, but working ?!? *)
   564 	       val tac_ = Apply_Method' (e_metID(*WN0402: see generate1 !?!*), 
   565 					 SOME t, is)
   566 	       val pos' = ((lev_on o lev_dn) p, Frm)
   567 	       val thy = assoc_thy "Isac.thy"
   568 	       val (_,_,_,pt') = (*implicit Take*)generate1 thy tac_ is pos' pt
   569 	       val (_,_,(pt'',_)) = complete_solve CompleteSubpbl [] (pt',pos')
   570 	       val newnds = children (get_nd pt'' p)
   571 	       val pt''' = ins_chn newnds pt p 
   572 	   (*complete_solve cuts branches after*)
   573 	   in ("detailrls", pt'''(*, get_formress [] ((lev_on o lev_dn) p)cn*),
   574 	       (p @ [length newnds], Res):pos') end
   575     end;
   576 
   577 
   578 
   579 (* val(mI,m)=m;val ppp=p;(*!!!*)val(p,p_)=pos;val(_,pt,_)=ppp(*!!!*);
   580    get_form ((mI,m):tac'_) ((p,p_):pos') ppp;
   581    *)
   582 fun get_form ((mI,m):tac'_) ((p,p_):pos') pt = 
   583   case applicable_in (p,p_) pt m of
   584     Notappl e => Error' (Error_ e)
   585   | Appl m => 
   586       (* val Appl m=applicable_in (p,p_) pt m;
   587          *)
   588       if member op = specsteps mI
   589 	then let val (_,_,f,_,_,_) = specify m (p,p_) [] pt
   590 	     in f end
   591       else let val (*_,_,f,_,_,_*)_ = solve (mI,m) (pt,(p,p_))
   592 	   in (*f*) EmptyMout end;
   593