test/Tools/isac/Knowledge/partial_fractions.sml
author Walther Neuper <neuper@ist.tugraz.at>
Tue, 11 Oct 2011 16:43:02 +0200
branchdecompose-isar
changeset 42310 55931ca19f4d
parent 42305 3ff2cf70f845
child 42313 f6a46e84067a
permissions -rw-r--r--
tuned
     1 (* Title:  partial fraction decomposition
     2    Author: Jan Rocnik
     3    (c) due to copyright terms
     4 *)
     5 
     6 "--------------------------------------------------------";
     7 "table of contents --------------------------------------";
     8 "--------------------------------------------------------";
     9 "----------- why helpless here ? ------------------------";
    10 "----------- why not nxt = Model_Problem here ? ---------";
    11 "--------------------------------------------------------";
    12 "--------------------------------------------------------";
    13 "--------------------------------------------------------";
    14 
    15 
    16 "----------- why helpless here ? ------------------------";
    17 "----------- why helpless here ? ------------------------";
    18 "----------- why helpless here ? ------------------------";
    19 val fmz = ["filterExpression (X z = 3 / (z - 1/4 + -1/8 * (1/(z::real))))", 
    20   "stepResponse (x[n::real]::bool)"];
    21 val (dI,pI,mI) = ("Isac", ["inverse", "Z_Transform", "SignalProcessing"], 
    22   ["SignalProcessing","Z_Transform","inverse"]);
    23 val (p,_,f,nxt,_,pt)  = CalcTreeTEST [(fmz, (dI,pI,mI))]; 
    24 val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Add_Given";
    25 val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Add_Find";
    26 val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Specify_Theory";
    27 val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Specify_Problem";
    28 val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Specify_Method";
    29 val (p,_,f,nxt,_,pt) = me nxt p [] pt; "nxt = Apply_Method";
    30 val (p,_,f,nxt,_,pt) = me nxt p [] pt; 
    31   "nxt = Rewrite (ruleZY, Inverse_Z_Transform.ruleZY) --> X z = 3 / (z - 1 / 4 + -1 / 8 * (1 / z))";
    32 val (p,_,f,nxt,_,pt) = me nxt p [] pt; 
    33   "nxt = Rewrite_Set norm_Rational --> X' z = 3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))";
    34 "~~~~~ fun me, args:"; val ((_,tac), (p:pos'), _, (pt:ptree)) = (nxt, p, [], pt);
    35 val ("ok", (_, _, ptp)) = locatetac tac (pt,p)
    36 val (pt, p) = ptp;
    37 "~~~~~ fun step, args:"; val (((ip as (_,p_)):pos'), ((ptp as (pt,p), tacis):calcstate)) = 
    38                            (p, ((pt, e_pos'),[]));
    39 val pIopt = get_pblID (pt,ip);
    40 ip = ([],Res); "false";
    41 tacis; " = []";
    42 pIopt; (* = SOME ["inverse", "Z_Transform", "SignalProcessing"]*)
    43 member op = [Pbl,Met] p_ andalso is_none (get_obj g_env pt (fst p)); "false";
    44 (*nxt_solve_ (pt,ip); "WAS isalist2list applied to NON-list 'no_meth'"
    45    THIS MEANS: replace no_meth by [no_meth] in Script.*)
    46 (*WAS val ("helpless",_) = step p ((pt, e_pos'),[]) *)
    47 (*WAS val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Empty_Tac instead SubProblem";*)
    48 
    49 "----------- why not nxt = Model_Problem here ? ---------";
    50 "----------- why not nxt = Model_Problem here ? ---------";
    51 "----------- why not nxt = Model_Problem here ? ---------";
    52