| author | webertj | 
| Mon, 10 Jul 2006 21:02:29 +0200 | |
| changeset 20066 | b045b835cb3b | 
| parent 19336 | fb5e19d26d5e | 
| child 20853 | 3ff5a2e05810 | 
| permissions | -rw-r--r-- | 
| 12869 | 1  | 
(* Title: HOL/Tools/svc_funcs.ML  | 
2  | 
ID: $Id$  | 
|
3  | 
Author: Lawrence C Paulson  | 
|
4  | 
Copyright 1999 University of Cambridge  | 
|
5  | 
||
6  | 
Translation functions for the interface to SVC  | 
|
7  | 
||
| 
19336
 
fb5e19d26d5e
removed some illegal characters: they were crashing SML/NJ
 
paulson 
parents: 
19277 
diff
changeset
 | 
8  | 
Based upon the work of Soren T. Heilmann  | 
| 12869 | 9  | 
|
10  | 
Integers and naturals are translated as follows:  | 
|
11  | 
In a positive context, replace x<y by x+1<=y  | 
|
12  | 
In a negative context, replace x<=y by x<y+1  | 
|
13  | 
In a negative context, replace x=y by x<y+1 & y<x+1  | 
|
14  | 
Biconditionals (if-and-only-iff) are expanded if they require such translations  | 
|
15  | 
in either operand.  | 
|
16  | 
||
17  | 
For each variable of type nat, an assumption is added that it is non-negative.  | 
|
18  | 
*)  | 
|
19  | 
||
20  | 
structure Svc =  | 
|
21  | 
struct  | 
|
22  | 
val trace = ref false;  | 
|
23  | 
||
24  | 
datatype expr =  | 
|
25  | 
Buildin of string * expr list  | 
|
26  | 
| Interp of string * expr list  | 
|
27  | 
| UnInterp of string * expr list  | 
|
| 16836 | 28  | 
| FalseExpr  | 
| 12869 | 29  | 
| TrueExpr  | 
| 
15965
 
f422f8283491
Use of IntInf.int instead of int in most numeric simprocs; avoids
 
paulson 
parents: 
15574 
diff
changeset
 | 
30  | 
| Int of IntInf.int  | 
| 
 
f422f8283491
Use of IntInf.int instead of int in most numeric simprocs; avoids
 
paulson 
parents: 
15574 
diff
changeset
 | 
31  | 
| Rat of IntInf.int * IntInf.int;  | 
| 12869 | 32  | 
|
| 16836 | 33  | 
fun signedInt i =  | 
| 
15965
 
f422f8283491
Use of IntInf.int instead of int in most numeric simprocs; avoids
 
paulson 
parents: 
15574 
diff
changeset
 | 
34  | 
if i < 0 then "-" ^ IntInf.toString (~i)  | 
| 
 
f422f8283491
Use of IntInf.int instead of int in most numeric simprocs; avoids
 
paulson 
parents: 
15574 
diff
changeset
 | 
35  | 
else IntInf.toString i;  | 
| 16836 | 36  | 
|
| 12869 | 37  | 
fun is_intnat T = T = HOLogic.intT orelse T = HOLogic.natT;  | 
| 16836 | 38  | 
|
| 12869 | 39  | 
fun is_numeric T = is_intnat T orelse T = HOLogic.realT;  | 
| 16836 | 40  | 
|
| 12869 | 41  | 
fun is_numeric_op T = is_numeric (domain_type T);  | 
42  | 
||
43  | 
fun toString t =  | 
|
| 16836 | 44  | 
let fun ue (Buildin(s, l)) =  | 
45  | 
             "(" ^ s ^ (Library.foldl (fn (a, b) => a ^ " " ^ (ue b)) ("", l)) ^ ") "
 | 
|
46  | 
| ue (Interp(s, l)) =  | 
|
47  | 
             "{" ^ s ^ (Library.foldl (fn (a, b) => a ^ " " ^ (ue b)) ("", l)) ^ "} "
 | 
|
48  | 
| ue (UnInterp(s, l)) =  | 
|
49  | 
             "(" ^ s ^ (Library.foldl (fn (a, b) => a ^ " " ^ (ue b)) ("", l)) ^ ") "
 | 
|
50  | 
| ue (FalseExpr) = "FALSE "  | 
|
51  | 
| ue (TrueExpr) = "TRUE "  | 
|
52  | 
| ue (Int i) = (signedInt i) ^ " "  | 
|
53  | 
| ue (Rat(i, j)) = (signedInt i) ^ "|" ^ (signedInt j) ^ " "  | 
|
| 12869 | 54  | 
in  | 
| 16836 | 55  | 
ue t  | 
| 12869 | 56  | 
end;  | 
57  | 
||
| 16836 | 58  | 
fun valid e =  | 
59  | 
let val svc_home = getenv "SVC_HOME"  | 
|
| 12869 | 60  | 
val svc_machine = getenv "SVC_MACHINE"  | 
61  | 
val check_valid = if svc_home = ""  | 
|
| 16836 | 62  | 
then error "Environment variable SVC_HOME not set"  | 
63  | 
else if svc_machine = ""  | 
|
64  | 
then error "Environment variable SVC_MACHINE not set"  | 
|
65  | 
else svc_home ^ "/" ^ svc_machine ^ "/bin/check_valid"  | 
|
| 12869 | 66  | 
val svc_input = toString e  | 
67  | 
      val _ = if !trace then tracing ("Calling SVC:\n" ^ svc_input) else ()
 | 
|
68  | 
val svc_input_file = File.tmp_path (Path.basic "SVM_in");  | 
|
69  | 
val svc_output_file = File.tmp_path (Path.basic "SVM_out");  | 
|
70  | 
val _ = (File.write svc_input_file svc_input;  | 
|
| 16836 | 71  | 
execute (check_valid ^ " -dump-result " ^  | 
72  | 
File.shell_path svc_output_file ^  | 
|
73  | 
" " ^ File.shell_path svc_input_file ^  | 
|
74  | 
">/dev/null 2>&1"))  | 
|
| 12869 | 75  | 
val svc_output =  | 
76  | 
(case Library.try File.read svc_output_file of  | 
|
| 15531 | 77  | 
SOME out => out  | 
78  | 
| NONE => error "SVC returned no output");  | 
|
| 12869 | 79  | 
in  | 
80  | 
      if ! trace then tracing ("SVC Returns:\n" ^ svc_output)
 | 
|
81  | 
else (File.rm svc_input_file; File.rm svc_output_file);  | 
|
82  | 
String.isPrefix "VALID" svc_output  | 
|
83  | 
end  | 
|
84  | 
||
| 16836 | 85  | 
 fun fail t = raise TERM ("SVC oracle", [t]);
 | 
| 12869 | 86  | 
|
87  | 
fun apply c args =  | 
|
88  | 
let val (ts, bs) = ListPair.unzip args  | 
|
89  | 
in (list_comb(c,ts), exists I bs) end;  | 
|
90  | 
||
91  | 
(*Determining whether the biconditionals must be unfolded: if there are  | 
|
92  | 
int or nat comparisons below*)  | 
|
93  | 
val iff_tag =  | 
|
94  | 
let fun tag t =  | 
|
| 16836 | 95  | 
let val (c,ts) = strip_comb t  | 
96  | 
in case c of  | 
|
97  | 
             Const("op &", _)   => apply c (map tag ts)
 | 
|
98  | 
           | Const("op |", _)   => apply c (map tag ts)
 | 
|
99  | 
           | Const("op -->", _) => apply c (map tag ts)
 | 
|
100  | 
           | Const("Not", _)    => apply c (map tag ts)
 | 
|
101  | 
           | Const("True", _)   => (c, false)
 | 
|
102  | 
           | Const("False", _)  => (c, false)
 | 
|
103  | 
           | Const("op =", Type ("fun", [T,_])) =>
 | 
|
104  | 
if T = HOLogic.boolT then  | 
|
105  | 
(*biconditional: with int/nat comparisons below?*)  | 
|
106  | 
let val [t1,t2] = ts  | 
|
107  | 
val (u1,b1) = tag t1  | 
|
108  | 
and (u2,b2) = tag t2  | 
|
109  | 
val cname = if b1 orelse b2 then "unfold" else "keep"  | 
|
110  | 
in  | 
|
111  | 
                        (Const ("SVC_Oracle.iff_" ^ cname, dummyT) $ u1 $ u2,
 | 
|
112  | 
b1 orelse b2)  | 
|
113  | 
end  | 
|
114  | 
else (*might be numeric equality*) (t, is_intnat T)  | 
|
| 19277 | 115  | 
           | Const("Orderings.less", Type ("fun", [T,_]))  => (t, is_intnat T)
 | 
116  | 
           | Const("Orderings.less_eq", Type ("fun", [T,_])) => (t, is_intnat T)
 | 
|
| 16836 | 117  | 
| _ => (t, false)  | 
118  | 
end  | 
|
| 12869 | 119  | 
in #1 o tag end;  | 
120  | 
||
121  | 
(*Map expression e to 0<=a --> e, where "a" is the name of a nat variable*)  | 
|
| 16836 | 122  | 
fun add_nat_var (a, e) =  | 
| 12869 | 123  | 
     Buildin("=>", [Buildin("<=", [Int 0, UnInterp (a, [])]),
 | 
| 16836 | 124  | 
e]);  | 
| 12869 | 125  | 
|
126  | 
fun param_string [] = ""  | 
|
127  | 
| param_string is = "_" ^ space_implode "_" (map string_of_int is)  | 
|
128  | 
||
129  | 
(*Translate an Isabelle formula into an SVC expression  | 
|
130  | 
pos ["positive"]: true if an assumption, false if a goal*)  | 
|
131  | 
fun expr_of pos t =  | 
|
132  | 
let  | 
|
133  | 
val params = rev (rename_wrt_term t (Term.strip_all_vars t))  | 
|
134  | 
and body = Term.strip_all_body t  | 
|
135  | 
val nat_vars = ref ([] : string list)  | 
|
136  | 
(*translation of a variable: record all natural numbers*)  | 
|
137  | 
fun trans_var (a,T,is) =  | 
|
| 16836 | 138  | 
(if T = HOLogic.natT then nat_vars := (a ins_string (!nat_vars))  | 
139  | 
else ();  | 
|
| 12869 | 140  | 
UnInterp (a ^ param_string is, []))  | 
141  | 
(*A variable, perhaps applied to a series of parameters*)  | 
|
142  | 
    fun var (Free(a,T), is)      = trans_var ("F_" ^ a, T, is)
 | 
|
143  | 
| var (Var((a, 0), T), is) = trans_var (a, T, is)  | 
|
| 16836 | 144  | 
| var (Bound i, is) =  | 
| 12869 | 145  | 
let val (a,T) = List.nth (params, i)  | 
| 16836 | 146  | 
          in  trans_var ("B_" ^ a, T, is)  end
 | 
| 12869 | 147  | 
| var (t $ Bound i, is) = var(t,i::is)  | 
148  | 
(*removing a parameter from a Var: the bound var index will  | 
|
149  | 
become part of the Var's name*)  | 
|
| 16836 | 150  | 
| var (t,_) = fail t;  | 
| 12869 | 151  | 
(*translation of a literal*)  | 
152  | 
    fun lit (Const("Numeral.number_of", _) $ w) =
 | 
|
153  | 
(HOLogic.dest_binum w handle TERM _ => raise Match)  | 
|
154  | 
      | lit (Const("0", _)) = 0
 | 
|
155  | 
      | lit (Const("1", _)) = 1
 | 
|
156  | 
(*translation of a literal expression [no variables]*)  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
157  | 
    fun litExp (Const("HOL.plus", T) $ x $ y) =
 | 
| 16836 | 158  | 
if is_numeric_op T then (litExp x) + (litExp y)  | 
159  | 
else fail t  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
160  | 
      | litExp (Const("HOL.minus", T) $ x $ y) =
 | 
| 16836 | 161  | 
if is_numeric_op T then (litExp x) - (litExp y)  | 
162  | 
else fail t  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
163  | 
      | litExp (Const("HOL.times", T) $ x $ y) =
 | 
| 16836 | 164  | 
if is_numeric_op T then (litExp x) * (litExp y)  | 
165  | 
else fail t  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
166  | 
      | litExp (Const("HOL.uminus", T) $ x)   =
 | 
| 16836 | 167  | 
if is_numeric_op T then ~(litExp x)  | 
168  | 
else fail t  | 
|
169  | 
| litExp t = lit t  | 
|
170  | 
handle Match => fail t  | 
|
| 12869 | 171  | 
(*translation of a real/rational expression*)  | 
172  | 
    fun suc t = Interp("+", [Int 1, t])
 | 
|
173  | 
    fun tm (Const("Suc", T) $ x) = suc (tm x)
 | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
174  | 
      | tm (Const("HOL.plus", T) $ x $ y) =
 | 
| 16836 | 175  | 
          if is_numeric_op T then Interp("+", [tm x, tm y])
 | 
176  | 
else fail t  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
177  | 
      | tm (Const("HOL.minus", T) $ x $ y) =
 | 
| 16836 | 178  | 
if is_numeric_op T then  | 
179  | 
              Interp("+", [tm x, Interp("*", [Int ~1, tm y])])
 | 
|
180  | 
else fail t  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
181  | 
      | tm (Const("HOL.times", T) $ x $ y) =
 | 
| 16836 | 182  | 
          if is_numeric_op T then Interp("*", [tm x, tm y])
 | 
183  | 
else fail t  | 
|
184  | 
      | tm (Const("RealDef.rinv", T) $ x) =
 | 
|
185  | 
if domain_type T = HOLogic.realT then  | 
|
186  | 
Rat(1, litExp x)  | 
|
187  | 
else fail t  | 
|
| 
19233
 
77ca20b0ed77
renamed HOL + - * etc. to HOL.plus HOL.minus HOL.times etc.
 
haftmann 
parents: 
17465 
diff
changeset
 | 
188  | 
      | tm (Const("HOL.uminus", T) $ x) =
 | 
| 16836 | 189  | 
          if is_numeric_op T then Interp("*", [Int ~1, tm x])
 | 
190  | 
else fail t  | 
|
191  | 
| tm t = Int (lit t)  | 
|
192  | 
handle Match => var (t,[])  | 
|
| 12869 | 193  | 
(*translation of a formula*)  | 
| 16836 | 194  | 
    and fm pos (Const("op &", _) $ p $ q) =
 | 
195  | 
            Buildin("AND", [fm pos p, fm pos q])
 | 
|
196  | 
      | fm pos (Const("op |", _) $ p $ q) =
 | 
|
197  | 
            Buildin("OR", [fm pos p, fm pos q])
 | 
|
198  | 
      | fm pos (Const("op -->", _) $ p $ q) =
 | 
|
199  | 
            Buildin("=>", [fm (not pos) p, fm pos q])
 | 
|
200  | 
      | fm pos (Const("Not", _) $ p) =
 | 
|
201  | 
            Buildin("NOT", [fm (not pos) p])
 | 
|
| 12869 | 202  | 
      | fm pos (Const("True", _)) = TrueExpr
 | 
203  | 
      | fm pos (Const("False", _)) = FalseExpr
 | 
|
| 16836 | 204  | 
      | fm pos (Const("SVC_Oracle.iff_keep", _) $ p $ q) =
 | 
205  | 
(*polarity doesn't matter*)  | 
|
206  | 
            Buildin("=", [fm pos p, fm pos q])
 | 
|
207  | 
      | fm pos (Const("SVC_Oracle.iff_unfold", _) $ p $ q) =
 | 
|
208  | 
            Buildin("AND",   (*unfolding uses both polarities*)
 | 
|
209  | 
                         [Buildin("=>", [fm (not pos) p, fm pos q]),
 | 
|
210  | 
                          Buildin("=>", [fm (not pos) q, fm pos p])])
 | 
|
211  | 
      | fm pos (t as Const("op =", Type ("fun", [T,_])) $ x $ y) =
 | 
|
212  | 
let val tx = tm x and ty = tm y  | 
|
213  | 
in if pos orelse T = HOLogic.realT then  | 
|
214  | 
                       Buildin("=", [tx, ty])
 | 
|
215  | 
else if is_intnat T then  | 
|
216  | 
                       Buildin("AND",
 | 
|
217  | 
                                    [Buildin("<", [tx, suc ty]),
 | 
|
218  | 
                                     Buildin("<", [ty, suc tx])])
 | 
|
219  | 
else fail t  | 
|
220  | 
end  | 
|
221  | 
(*inequalities: possible types are nat, int, real*)  | 
|
| 19277 | 222  | 
      | fm pos (t as Const("Orderings.less",  Type ("fun", [T,_])) $ x $ y) =
 | 
| 16836 | 223  | 
if not pos orelse T = HOLogic.realT then  | 
224  | 
                Buildin("<", [tm x, tm y])
 | 
|
225  | 
else if is_intnat T then  | 
|
226  | 
                Buildin("<=", [suc (tm x), tm y])
 | 
|
227  | 
else fail t  | 
|
| 19277 | 228  | 
      | fm pos (t as Const("Orderings.less_eq",  Type ("fun", [T,_])) $ x $ y) =
 | 
| 16836 | 229  | 
if pos orelse T = HOLogic.realT then  | 
230  | 
                Buildin("<=", [tm x, tm y])
 | 
|
231  | 
else if is_intnat T then  | 
|
232  | 
                Buildin("<", [tm x, suc (tm y)])
 | 
|
233  | 
else fail t  | 
|
| 12869 | 234  | 
| fm pos t = var(t,[]);  | 
235  | 
(*entry point, and translation of a meta-formula*)  | 
|
236  | 
      fun mt pos ((c as Const("Trueprop", _)) $ p) = fm pos (iff_tag p)
 | 
|
| 16836 | 237  | 
        | mt pos ((c as Const("==>", _)) $ p $ q) =
 | 
238  | 
            Buildin("=>", [mt (not pos) p, mt pos q])
 | 
|
239  | 
| mt pos t = fm pos (iff_tag t) (*it might be a formula*)  | 
|
| 12869 | 240  | 
|
241  | 
val body_e = mt pos body (*evaluate now to assign into !nat_vars*)  | 
|
| 16836 | 242  | 
in  | 
243  | 
foldr add_nat_var body_e (!nat_vars)  | 
|
| 12869 | 244  | 
end;  | 
245  | 
||
246  | 
||
247  | 
(*The oracle proves the given formula t, if possible*)  | 
|
| 16836 | 248  | 
fun oracle thy t =  | 
249  | 
(conditional (! trace) (fn () =>  | 
|
250  | 
    tracing ("SVC oracle: problem is\n" ^ Sign.string_of_term thy t));
 | 
|
251  | 
if valid (expr_of false t) then t else fail t);  | 
|
| 12869 | 252  | 
|
253  | 
end;  |