src/HOL/Tools/Predicate_Compile/code_prolog.ML
author bulwahn
Mon Sep 27 12:23:00 2010 +0200 (2010-09-27)
changeset 39724 ada0cd4900c1
parent 39546 bfe10da7d764
child 39761 c2a76ec6e2d9
permissions -rw-r--r--
adding further tracing messages; tuned
bulwahn@38073
     1
(*  Title:      HOL/Tools/Predicate_Compile/code_prolog.ML
bulwahn@38073
     2
    Author:     Lukas Bulwahn, TU Muenchen
bulwahn@38073
     3
bulwahn@38073
     4
Prototype of an code generator for logic programming languages (a.k.a. Prolog)
bulwahn@38073
     5
*)
bulwahn@38073
     6
bulwahn@38073
     7
signature CODE_PROLOG =
bulwahn@38073
     8
sig
bulwahn@38792
     9
  datatype prolog_system = SWI_PROLOG | YAP
bulwahn@38792
    10
  type code_options =
bulwahn@38947
    11
    {ensure_groundness : bool,
bulwahn@38947
    12
     limited_types : (typ * int) list,
bulwahn@38959
    13
     limited_predicates : (string list * int) list,
bulwahn@38947
    14
     replacing : ((string * string) * string) list,
bulwahn@39461
    15
     manual_reorder : ((string * int) * int list) list}
bulwahn@39461
    16
  val set_ensure_groundness : code_options -> code_options
bulwahn@39461
    17
  val map_limit_predicates : ((string list * int) list -> (string list * int) list)
bulwahn@39461
    18
    -> code_options -> code_options
bulwahn@38950
    19
  val code_options_of : theory -> code_options 
bulwahn@38950
    20
  val map_code_options : (code_options -> code_options) -> theory -> theory
bulwahn@39461
    21
  
bulwahn@38113
    22
  datatype arith_op = Plus | Minus
bulwahn@38112
    23
  datatype prol_term = Var of string | Cons of string | AppF of string * prol_term list
bulwahn@38113
    24
    | Number of int | ArithOp of arith_op * prol_term list;
bulwahn@38113
    25
  datatype prem = Conj of prem list
bulwahn@38113
    26
    | Rel of string * prol_term list | NotRel of string * prol_term list
bulwahn@38113
    27
    | Eq of prol_term * prol_term | NotEq of prol_term * prol_term
bulwahn@38727
    28
    | ArithEq of prol_term * prol_term | NotArithEq of prol_term * prol_term
bulwahn@38727
    29
    | Ground of string * typ;
bulwahn@38727
    30
bulwahn@38079
    31
  type clause = ((string * prol_term list) * prem);
bulwahn@38073
    32
  type logic_program = clause list;
bulwahn@38079
    33
  type constant_table = (string * string) list
bulwahn@39461
    34
  
bulwahn@39461
    35
  val generate : Predicate_Compile_Aux.mode option * bool ->
bulwahn@39461
    36
    Proof.context -> string -> (logic_program * constant_table)
bulwahn@38079
    37
  val write_program : logic_program -> string
bulwahn@39464
    38
  val run : (Time.time * prolog_system) -> logic_program -> (string * prol_term list) -> string list -> int option -> prol_term list list
bulwahn@39461
    39
  
bulwahn@39461
    40
  val quickcheck : Proof.context -> term -> int -> term list option * (bool list * bool)
bulwahn@38733
    41
bulwahn@38079
    42
  val trace : bool Unsynchronized.ref
bulwahn@38947
    43
  
bulwahn@38947
    44
  val replace : ((string * string) * string) -> logic_program -> logic_program
bulwahn@38073
    45
end;
bulwahn@38073
    46
bulwahn@38073
    47
structure Code_Prolog : CODE_PROLOG =
bulwahn@38073
    48
struct
bulwahn@38073
    49
bulwahn@38079
    50
(* diagnostic tracing *)
bulwahn@38079
    51
bulwahn@38079
    52
val trace = Unsynchronized.ref false
bulwahn@38079
    53
bulwahn@38079
    54
fun tracing s = if !trace then Output.tracing s else () 
bulwahn@38727
    55
bulwahn@38727
    56
(* code generation options *)
bulwahn@38727
    57
bulwahn@38727
    58
bulwahn@38792
    59
type code_options =
bulwahn@38947
    60
  {ensure_groundness : bool,
bulwahn@38947
    61
   limited_types : (typ * int) list,
bulwahn@38959
    62
   limited_predicates : (string list * int) list,
bulwahn@38947
    63
   replacing : ((string * string) * string) list,
bulwahn@39461
    64
   manual_reorder : ((string * int) * int list) list}
bulwahn@39461
    65
bulwahn@38792
    66
bulwahn@39461
    67
fun set_ensure_groundness {ensure_groundness, limited_types, limited_predicates,
bulwahn@39461
    68
  replacing, manual_reorder} =
bulwahn@39461
    69
  {ensure_groundness = true, limited_types = limited_types,
bulwahn@39461
    70
   limited_predicates = limited_predicates, replacing = replacing,
bulwahn@39461
    71
   manual_reorder = manual_reorder}
bulwahn@39461
    72
bulwahn@39461
    73
fun map_limit_predicates f {ensure_groundness, limited_types, limited_predicates,
bulwahn@39461
    74
  replacing, manual_reorder} =
bulwahn@39461
    75
  {ensure_groundness = ensure_groundness, limited_types = limited_types,
bulwahn@39461
    76
   limited_predicates = f limited_predicates, replacing = replacing,
bulwahn@39461
    77
   manual_reorder = manual_reorder}
bulwahn@39461
    78
  
bulwahn@38950
    79
structure Options = Theory_Data
bulwahn@38950
    80
(
bulwahn@38950
    81
  type T = code_options
bulwahn@38950
    82
  val empty = {ensure_groundness = false,
bulwahn@39461
    83
    limited_types = [], limited_predicates = [], replacing = [], manual_reorder = []}
bulwahn@38950
    84
  val extend = I;
bulwahn@38950
    85
  fun merge
bulwahn@38950
    86
    ({ensure_groundness = ensure_groundness1, limited_types = limited_types1,
bulwahn@38960
    87
      limited_predicates = limited_predicates1, replacing = replacing1,
bulwahn@39461
    88
      manual_reorder = manual_reorder1},
bulwahn@38950
    89
     {ensure_groundness = ensure_groundness2, limited_types = limited_types2,
bulwahn@38960
    90
      limited_predicates = limited_predicates2, replacing = replacing2,
bulwahn@39461
    91
      manual_reorder = manual_reorder2}) =
bulwahn@38950
    92
    {ensure_groundness = ensure_groundness1 orelse ensure_groundness2,
bulwahn@38950
    93
     limited_types = AList.merge (op =) (K true) (limited_types1, limited_types2),
bulwahn@38950
    94
     limited_predicates = AList.merge (op =) (K true) (limited_predicates1, limited_predicates2),
bulwahn@38960
    95
     manual_reorder = AList.merge (op =) (K true) (manual_reorder1, manual_reorder2),
bulwahn@39461
    96
     replacing = Library.merge (op =) (replacing1, replacing2)};
bulwahn@38950
    97
);
bulwahn@38950
    98
bulwahn@38950
    99
val code_options_of = Options.get
bulwahn@38950
   100
bulwahn@38950
   101
val map_code_options = Options.map
bulwahn@38727
   102
bulwahn@39461
   103
(* system configuration *)
bulwahn@39461
   104
bulwahn@39461
   105
datatype prolog_system = SWI_PROLOG | YAP
bulwahn@39461
   106
bulwahn@39462
   107
fun string_of_system SWI_PROLOG = "swiprolog"
bulwahn@39462
   108
  | string_of_system YAP = "yap"
bulwahn@39462
   109
bulwahn@39461
   110
type system_configuration = {timeout : Time.time, prolog_system : prolog_system}
bulwahn@39462
   111
                                                
bulwahn@39461
   112
structure System_Config = Generic_Data
bulwahn@39461
   113
(
bulwahn@39461
   114
  type T = system_configuration
bulwahn@39461
   115
  val empty = {timeout = Time.fromSeconds 10, prolog_system = SWI_PROLOG}
bulwahn@39461
   116
  val extend = I;
bulwahn@39461
   117
  fun merge ({timeout = timeout1, prolog_system = prolog_system1},
bulwahn@39461
   118
        {timeout = timeout2, prolog_system = prolog_system2}) =
bulwahn@39461
   119
    {timeout = timeout1, prolog_system = prolog_system1}
bulwahn@39461
   120
)
bulwahn@39461
   121
bulwahn@38073
   122
(* general string functions *)
bulwahn@38073
   123
bulwahn@38073
   124
val first_upper = implode o nth_map 0 Symbol.to_ascii_upper o explode;
bulwahn@38073
   125
val first_lower = implode o nth_map 0 Symbol.to_ascii_lower o explode;
bulwahn@38073
   126
bulwahn@38073
   127
(* internal program representation *)
bulwahn@38073
   128
bulwahn@38113
   129
datatype arith_op = Plus | Minus
bulwahn@38113
   130
bulwahn@38112
   131
datatype prol_term = Var of string | Cons of string | AppF of string * prol_term list
bulwahn@38113
   132
  | Number of int | ArithOp of arith_op * prol_term list;
bulwahn@38113
   133
bulwahn@38735
   134
fun dest_Var (Var v) = v
bulwahn@38735
   135
bulwahn@38735
   136
fun add_vars (Var v) = insert (op =) v
bulwahn@38735
   137
  | add_vars (ArithOp (_, ts)) = fold add_vars ts
bulwahn@38735
   138
  | add_vars (AppF (_, ts)) = fold add_vars ts
bulwahn@38735
   139
  | add_vars _ = I
bulwahn@38735
   140
bulwahn@38735
   141
fun map_vars f (Var v) = Var (f v)
bulwahn@38735
   142
  | map_vars f (ArithOp (opr, ts)) = ArithOp (opr, map (map_vars f) ts)
bulwahn@38735
   143
  | map_vars f (AppF (fs, ts)) = AppF (fs, map (map_vars f) ts)
bulwahn@38735
   144
  | map_vars f t = t
bulwahn@38735
   145
  
bulwahn@38728
   146
fun maybe_AppF (c, []) = Cons c
bulwahn@38728
   147
  | maybe_AppF (c, xs) = AppF (c, xs)
bulwahn@38728
   148
bulwahn@38113
   149
fun is_Var (Var _) = true
bulwahn@38113
   150
  | is_Var _ = false
bulwahn@38113
   151
bulwahn@38113
   152
fun is_arith_term (Var _) = true
bulwahn@38113
   153
  | is_arith_term (Number _) = true
bulwahn@38113
   154
  | is_arith_term (ArithOp (_, operands)) = forall is_arith_term operands
bulwahn@38113
   155
  | is_arith_term _ = false
bulwahn@38073
   156
bulwahn@38081
   157
fun string_of_prol_term (Var s) = "Var " ^ s
bulwahn@38075
   158
  | string_of_prol_term (Cons s) = "Cons " ^ s
bulwahn@38075
   159
  | string_of_prol_term (AppF (f, args)) = f ^ "(" ^ commas (map string_of_prol_term args) ^ ")" 
bulwahn@38112
   160
  | string_of_prol_term (Number n) = "Number " ^ string_of_int n
bulwahn@38075
   161
bulwahn@38113
   162
datatype prem = Conj of prem list
bulwahn@38113
   163
  | Rel of string * prol_term list | NotRel of string * prol_term list
bulwahn@38113
   164
  | Eq of prol_term * prol_term | NotEq of prol_term * prol_term
bulwahn@38727
   165
  | ArithEq of prol_term * prol_term | NotArithEq of prol_term * prol_term
bulwahn@38727
   166
  | Ground of string * typ;
bulwahn@38735
   167
bulwahn@38073
   168
fun dest_Rel (Rel (c, ts)) = (c, ts)
bulwahn@38735
   169
bulwahn@38735
   170
fun map_term_prem f (Conj prems) = Conj (map (map_term_prem f) prems)
bulwahn@38735
   171
  | map_term_prem f (Rel (r, ts)) = Rel (r, map f ts)
bulwahn@38735
   172
  | map_term_prem f (NotRel (r, ts)) = NotRel (r, map f ts)
bulwahn@38735
   173
  | map_term_prem f (Eq (l, r)) = Eq (f l, f r)
bulwahn@38735
   174
  | map_term_prem f (NotEq (l, r)) = NotEq (f l, f r)
bulwahn@38735
   175
  | map_term_prem f (ArithEq (l, r)) = ArithEq (f l, f r)
bulwahn@38735
   176
  | map_term_prem f (NotArithEq (l, r)) = NotArithEq (f l, f r)
bulwahn@38735
   177
  | map_term_prem f (Ground (v, T)) = Ground (dest_Var (f (Var v)), T)
bulwahn@38735
   178
bulwahn@38735
   179
fun fold_prem_terms f (Conj prems) = fold (fold_prem_terms f) prems
bulwahn@38735
   180
  | fold_prem_terms f (Rel (_, ts)) = fold f ts
bulwahn@38735
   181
  | fold_prem_terms f (NotRel (_, ts)) = fold f ts
bulwahn@38735
   182
  | fold_prem_terms f (Eq (l, r)) = f l #> f r
bulwahn@38735
   183
  | fold_prem_terms f (NotEq (l, r)) = f l #> f r
bulwahn@38735
   184
  | fold_prem_terms f (ArithEq (l, r)) = f l #> f r
bulwahn@38735
   185
  | fold_prem_terms f (NotArithEq (l, r)) = f l #> f r
bulwahn@38735
   186
  | fold_prem_terms f (Ground (v, T)) = f (Var v)
bulwahn@38735
   187
  
bulwahn@38079
   188
type clause = ((string * prol_term list) * prem);
bulwahn@38073
   189
bulwahn@38073
   190
type logic_program = clause list;
bulwahn@38073
   191
bulwahn@38073
   192
(* translation from introduction rules to internal representation *)
bulwahn@38073
   193
bulwahn@38958
   194
fun mk_conform f empty avoid name =
bulwahn@38956
   195
  let
bulwahn@38956
   196
    fun dest_Char (Symbol.Char c) = c
bulwahn@38956
   197
    val name' = space_implode "" (map (dest_Char o Symbol.decode)
bulwahn@38956
   198
      (filter (fn s => Symbol.is_ascii_letter s orelse Symbol.is_ascii_digit s)
bulwahn@38956
   199
        (Symbol.explode name)))
bulwahn@38958
   200
    val name'' = f (if name' = "" then empty else name')
bulwahn@38958
   201
  in (if member (op =) avoid name'' then Name.variant avoid name'' else name'') end
bulwahn@38956
   202
bulwahn@38079
   203
(** constant table **)
bulwahn@38079
   204
bulwahn@38079
   205
type constant_table = (string * string) list
bulwahn@38079
   206
bulwahn@38079
   207
fun declare_consts consts constant_table =
bulwahn@38956
   208
  let
bulwahn@38956
   209
    fun update' c table =
bulwahn@38956
   210
      if AList.defined (op =) table c then table else
bulwahn@38956
   211
        let
bulwahn@38958
   212
          val c' = mk_conform first_lower "pred" (map snd table) (Long_Name.base_name c)
bulwahn@38956
   213
        in
bulwahn@38956
   214
          AList.update (op =) (c, c') table
bulwahn@38956
   215
        end
bulwahn@38956
   216
  in
bulwahn@38956
   217
    fold update' consts constant_table
bulwahn@38956
   218
  end
bulwahn@38079
   219
  
bulwahn@38079
   220
fun translate_const constant_table c =
bulwahn@38079
   221
  case AList.lookup (op =) constant_table c of
bulwahn@38079
   222
    SOME c' => c'
bulwahn@38079
   223
  | NONE => error ("No such constant: " ^ c)
bulwahn@38073
   224
bulwahn@38079
   225
fun inv_lookup _ [] _ = NONE
bulwahn@38079
   226
  | inv_lookup eq ((key, value)::xs) value' =
bulwahn@38079
   227
      if eq (value', value) then SOME key
bulwahn@38079
   228
      else inv_lookup eq xs value';
bulwahn@38079
   229
bulwahn@38079
   230
fun restore_const constant_table c =
bulwahn@38079
   231
  case inv_lookup (op =) constant_table c of
bulwahn@38079
   232
    SOME c' => c'
bulwahn@38079
   233
  | NONE => error ("No constant corresponding to "  ^ c)
bulwahn@38727
   234
bulwahn@38079
   235
(** translation of terms, literals, premises, and clauses **)
bulwahn@38079
   236
bulwahn@38113
   237
fun translate_arith_const @{const_name "Groups.plus_class.plus"} = SOME Plus
bulwahn@38113
   238
  | translate_arith_const @{const_name "Groups.minus_class.minus"} = SOME Minus
bulwahn@38113
   239
  | translate_arith_const _ = NONE
bulwahn@38113
   240
bulwahn@38734
   241
fun mk_nat_term constant_table n =
bulwahn@38734
   242
  let
bulwahn@38734
   243
    val zero = translate_const constant_table @{const_name "Groups.zero_class.zero"}
bulwahn@38734
   244
    val Suc = translate_const constant_table @{const_name "Suc"}
bulwahn@38734
   245
  in funpow n (fn t => AppF (Suc, [t])) (Cons zero) end
bulwahn@38734
   246
bulwahn@38079
   247
fun translate_term ctxt constant_table t =
bulwahn@38112
   248
  case try HOLogic.dest_number t of
bulwahn@38112
   249
    SOME (@{typ "int"}, n) => Number n
bulwahn@38734
   250
  | SOME (@{typ "nat"}, n) => mk_nat_term constant_table n
bulwahn@38112
   251
  | NONE =>
bulwahn@38112
   252
      (case strip_comb t of
bulwahn@38112
   253
        (Free (v, T), []) => Var v 
bulwahn@38112
   254
      | (Const (c, _), []) => Cons (translate_const constant_table c)
bulwahn@38112
   255
      | (Const (c, _), args) =>
bulwahn@38113
   256
        (case translate_arith_const c of
bulwahn@38113
   257
          SOME aop => ArithOp (aop, map (translate_term ctxt constant_table) args)
bulwahn@38113
   258
        | NONE =>                                                             
bulwahn@38113
   259
            AppF (translate_const constant_table c, map (translate_term ctxt constant_table) args))
bulwahn@38112
   260
      | _ => error ("illegal term for translation: " ^ Syntax.string_of_term ctxt t))
bulwahn@38073
   261
bulwahn@38079
   262
fun translate_literal ctxt constant_table t =
bulwahn@38073
   263
  case strip_comb t of
haftmann@38864
   264
    (Const (@{const_name HOL.eq}, _), [l, r]) =>
bulwahn@38113
   265
      let
bulwahn@38113
   266
        val l' = translate_term ctxt constant_table l
bulwahn@38113
   267
        val r' = translate_term ctxt constant_table r
bulwahn@38113
   268
      in
bulwahn@38728
   269
        (if is_Var l' andalso is_arith_term r' andalso not (is_Var r') then ArithEq else Eq) (l', r')
bulwahn@38113
   270
      end
bulwahn@38079
   271
  | (Const (c, _), args) =>
bulwahn@38079
   272
      Rel (translate_const constant_table c, map (translate_term ctxt constant_table) args)
bulwahn@38073
   273
  | _ => error ("illegal literal for translation: " ^ Syntax.string_of_term ctxt t)
bulwahn@38073
   274
bulwahn@38073
   275
fun NegRel_of (Rel lit) = NotRel lit
bulwahn@38073
   276
  | NegRel_of (Eq eq) = NotEq eq
bulwahn@38113
   277
  | NegRel_of (ArithEq eq) = NotArithEq eq
bulwahn@38113
   278
bulwahn@38727
   279
fun mk_groundness_prems t = map Ground (Term.add_frees t [])
bulwahn@38727
   280
  
bulwahn@38792
   281
fun translate_prem ensure_groundness ctxt constant_table t =  
bulwahn@38073
   282
    case try HOLogic.dest_not t of
bulwahn@38727
   283
      SOME t =>
bulwahn@38792
   284
        if ensure_groundness then
bulwahn@38727
   285
          Conj (mk_groundness_prems t @ [NegRel_of (translate_literal ctxt constant_table t)])
bulwahn@38727
   286
        else
bulwahn@38727
   287
          NegRel_of (translate_literal ctxt constant_table t)
bulwahn@38079
   288
    | NONE => translate_literal ctxt constant_table t
bulwahn@38114
   289
    
bulwahn@38114
   290
fun imp_prems_conv cv ct =
bulwahn@38114
   291
  case Thm.term_of ct of
bulwahn@38114
   292
    Const ("==>", _) $ _ $ _ => Conv.combination_conv (Conv.arg_conv cv) (imp_prems_conv cv) ct
bulwahn@38114
   293
  | _ => Conv.all_conv ct
bulwahn@38114
   294
bulwahn@38114
   295
fun Trueprop_conv cv ct =
bulwahn@38114
   296
  case Thm.term_of ct of
haftmann@38558
   297
    Const (@{const_name Trueprop}, _) $ _ => Conv.arg_conv cv ct  
bulwahn@38114
   298
  | _ => raise Fail "Trueprop_conv"
bulwahn@38114
   299
bulwahn@38114
   300
fun preprocess_intro thy rule =
bulwahn@38114
   301
  Conv.fconv_rule
bulwahn@38114
   302
    (imp_prems_conv
bulwahn@38114
   303
      (Trueprop_conv (Conv.try_conv (Conv.rewr_conv @{thm Predicate.eq_is_eq}))))
bulwahn@38114
   304
    (Thm.transfer thy rule)
bulwahn@38114
   305
bulwahn@38792
   306
fun translate_intros ensure_groundness ctxt gr const constant_table =
bulwahn@38073
   307
  let
bulwahn@38114
   308
    val intros = map (preprocess_intro (ProofContext.theory_of ctxt)) (Graph.get_node gr const)
bulwahn@38073
   309
    val (intros', ctxt') = Variable.import_terms true (map prop_of intros) ctxt
bulwahn@38079
   310
    val constant_table' = declare_consts (fold Term.add_const_names intros' []) constant_table
bulwahn@38073
   311
    fun translate_intro intro =
bulwahn@38073
   312
      let
bulwahn@38073
   313
        val head = HOLogic.dest_Trueprop (Logic.strip_imp_concl intro)
bulwahn@38727
   314
        val prems = map HOLogic.dest_Trueprop (Logic.strip_imp_prems intro)
bulwahn@38792
   315
        val prems' = Conj (map (translate_prem ensure_groundness ctxt' constant_table') prems)
bulwahn@38079
   316
        val clause = (dest_Rel (translate_literal ctxt' constant_table' head), prems')
bulwahn@38073
   317
      in clause end
bulwahn@39724
   318
  in
bulwahn@39724
   319
    (map translate_intro intros', constant_table')
bulwahn@39724
   320
  end
bulwahn@38073
   321
bulwahn@38731
   322
fun depending_preds_of (key, intros) =
bulwahn@38731
   323
  fold Term.add_const_names (map Thm.prop_of intros) []
bulwahn@38731
   324
bulwahn@38731
   325
fun add_edges edges_of key G =
bulwahn@38731
   326
  let
bulwahn@38731
   327
    fun extend' key (G, visited) = 
bulwahn@38731
   328
      case try (Graph.get_node G) key of
bulwahn@38731
   329
          SOME v =>
bulwahn@38731
   330
            let
bulwahn@38731
   331
              val new_edges = filter (fn k => is_some (try (Graph.get_node G) k)) (edges_of (key, v))
bulwahn@38731
   332
              val (G', visited') = fold extend'
bulwahn@38731
   333
                (subtract (op =) (key :: visited) new_edges) (G, key :: visited)
bulwahn@38731
   334
            in
bulwahn@38731
   335
              (fold (Graph.add_edge o (pair key)) new_edges G', visited')
bulwahn@38731
   336
            end
bulwahn@38731
   337
        | NONE => (G, visited)
bulwahn@38731
   338
  in
bulwahn@38731
   339
    fst (extend' key (G, []))
bulwahn@38731
   340
  end
bulwahn@38731
   341
bulwahn@39183
   342
fun print_intros ctxt gr consts =
bulwahn@39183
   343
  tracing (cat_lines (map (fn const =>
bulwahn@39183
   344
    "Constant " ^ const ^ "has intros:\n" ^
bulwahn@39183
   345
    cat_lines (map (Display.string_of_thm ctxt) (Graph.get_node gr const))) consts))
bulwahn@39461
   346
bulwahn@39461
   347
(* translation of moded predicates *)
bulwahn@39461
   348
bulwahn@39461
   349
(** generating graph of moded predicates **)
bulwahn@39461
   350
bulwahn@39461
   351
(* could be moved to Predicate_Compile_Core *)
bulwahn@39461
   352
fun requires_modes polarity cls =
bulwahn@39461
   353
  let
bulwahn@39461
   354
    fun req_mode_of pol (t, derivation) =
bulwahn@39461
   355
      (case fst (strip_comb t) of
bulwahn@39461
   356
        Const (c, _) => SOME (c, (pol, Predicate_Compile_Core.head_mode_of derivation))
bulwahn@39461
   357
      | _ => NONE)
bulwahn@39461
   358
    fun req (Predicate_Compile_Aux.Prem t, derivation) = req_mode_of polarity (t, derivation)
bulwahn@39461
   359
      | req (Predicate_Compile_Aux.Negprem t, derivation) = req_mode_of (not polarity) (t, derivation)
bulwahn@39461
   360
      | req _ = NONE
bulwahn@39461
   361
  in      
bulwahn@39461
   362
    maps (fn (_, prems) => map_filter req prems) cls
bulwahn@39461
   363
  end
bulwahn@39461
   364
 
bulwahn@39461
   365
structure Mode_Graph = Graph(type key = string * (bool * Predicate_Compile_Aux.mode)
bulwahn@39461
   366
  val ord = prod_ord fast_string_ord (prod_ord bool_ord Predicate_Compile_Aux.mode_ord));
bulwahn@39461
   367
bulwahn@39461
   368
fun mk_moded_clauses_graph ctxt scc gr =
bulwahn@39461
   369
  let
bulwahn@39461
   370
    val options = Predicate_Compile_Aux.default_options
bulwahn@39461
   371
    val mode_analysis_options =
bulwahn@39461
   372
      {use_random = true, reorder_premises = true, infer_pos_and_neg_modes = true}
bulwahn@39461
   373
    fun infer prednames (gr, (pos_modes, neg_modes, random)) =
bulwahn@39461
   374
      let
bulwahn@39461
   375
        val (lookup_modes, lookup_neg_modes, needs_random) =
bulwahn@39461
   376
          ((fn s => the (AList.lookup (op =) pos_modes s)),
bulwahn@39461
   377
           (fn s => the (AList.lookup (op =) neg_modes s)),
bulwahn@39461
   378
           (fn s => member (op =) (the (AList.lookup (op =) random s))))
bulwahn@39461
   379
        val (preds, all_vs, param_vs, all_modes, clauses) =
bulwahn@39461
   380
          Predicate_Compile_Core.prepare_intrs options ctxt prednames
bulwahn@39461
   381
            (maps (Predicate_Compile_Core.intros_of ctxt) prednames)
bulwahn@39461
   382
        val ((moded_clauses, random'), _) =
bulwahn@39461
   383
          Predicate_Compile_Core.infer_modes mode_analysis_options options 
bulwahn@39461
   384
            (lookup_modes, lookup_neg_modes, needs_random) ctxt preds all_modes param_vs clauses
bulwahn@39461
   385
        val modes = map (fn (p, mps) => (p, map fst mps)) moded_clauses
bulwahn@39461
   386
        val pos_modes' = map (apsnd (map_filter (fn (true, m) => SOME m | _ => NONE))) modes
bulwahn@39461
   387
        val neg_modes' = map (apsnd (map_filter (fn (false, m) => SOME m | _ => NONE))) modes
bulwahn@39461
   388
        val _ = tracing ("Inferred modes:\n" ^
bulwahn@39461
   389
          cat_lines (map (fn (s, ms) => s ^ ": " ^ commas (map
bulwahn@39461
   390
            (fn (p, m) => Predicate_Compile_Aux.string_of_mode m ^ (if p then "pos" else "neg")) ms)) modes))
bulwahn@39461
   391
        val gr' = gr
bulwahn@39461
   392
          |> fold (fn (p, mps) => fold (fn (mode, cls) =>
bulwahn@39461
   393
                Mode_Graph.new_node ((p, mode), cls)) mps)
bulwahn@39461
   394
            moded_clauses
bulwahn@39461
   395
          |> fold (fn (p, mps) => fold (fn (mode, cls) => fold (fn req =>
bulwahn@39461
   396
              Mode_Graph.add_edge ((p, mode), req)) (requires_modes (fst mode) cls)) mps)
bulwahn@39461
   397
            moded_clauses
bulwahn@39461
   398
      in
bulwahn@39461
   399
        (gr', (AList.merge (op =) (op =) (pos_modes, pos_modes'),
bulwahn@39461
   400
          AList.merge (op =) (op =) (neg_modes, neg_modes'),
bulwahn@39461
   401
          AList.merge (op =) (op =) (random, random')))
bulwahn@39461
   402
      end
bulwahn@39461
   403
  in  
bulwahn@39461
   404
    fst (fold infer (rev scc) (Mode_Graph.empty, ([], [], []))) 
bulwahn@39461
   405
  end
bulwahn@39461
   406
bulwahn@39461
   407
fun declare_moded_predicate moded_preds table =
bulwahn@39461
   408
  let
bulwahn@39461
   409
    fun update' (p as (pred, (pol, mode))) table =
bulwahn@39461
   410
      if AList.defined (op =) table p then table else
bulwahn@39461
   411
        let
bulwahn@39461
   412
          val name = Long_Name.base_name pred ^ (if pol then "p" else "n")
bulwahn@39461
   413
            ^ Predicate_Compile_Aux.ascii_string_of_mode mode
bulwahn@39461
   414
          val p' = mk_conform first_lower "pred" (map snd table) name
bulwahn@39461
   415
        in
bulwahn@39461
   416
          AList.update (op =) (p, p') table
bulwahn@39461
   417
        end
bulwahn@39461
   418
  in
bulwahn@39461
   419
    fold update' moded_preds table
bulwahn@39461
   420
  end
bulwahn@39461
   421
bulwahn@39461
   422
fun mk_program ctxt moded_gr moded_preds (prog, (moded_pred_table, constant_table)) =
bulwahn@39461
   423
  let
bulwahn@39461
   424
    val moded_pred_table' = declare_moded_predicate moded_preds moded_pred_table
bulwahn@39461
   425
    fun mk_literal pol derivation constant_table' t =
bulwahn@39461
   426
      let
bulwahn@39461
   427
        val (p, args) = strip_comb t
bulwahn@39461
   428
        val mode = Predicate_Compile_Core.head_mode_of derivation 
bulwahn@39461
   429
        val name = fst (dest_Const p)
bulwahn@39461
   430
        
bulwahn@39461
   431
        val p' = the (AList.lookup (op =) moded_pred_table' (name, (pol, mode)))
bulwahn@39461
   432
        val args' = map (translate_term ctxt constant_table') args
bulwahn@39461
   433
      in
bulwahn@39461
   434
        Rel (p', args')
bulwahn@39461
   435
      end
bulwahn@39461
   436
    fun mk_prem pol (indprem, derivation) constant_table =
bulwahn@39461
   437
      case indprem of
bulwahn@39461
   438
        Predicate_Compile_Aux.Generator (s, T) => (Ground (s, T), constant_table)
bulwahn@39461
   439
      | _ =>
bulwahn@39461
   440
        declare_consts (Term.add_const_names (Predicate_Compile_Aux.dest_indprem indprem) []) constant_table
bulwahn@39461
   441
        |> (fn constant_table' =>
bulwahn@39461
   442
          (case indprem of Predicate_Compile_Aux.Negprem t =>
bulwahn@39461
   443
            NegRel_of (mk_literal (not pol) derivation constant_table' t)
bulwahn@39461
   444
          | _ =>
bulwahn@39461
   445
            mk_literal pol derivation constant_table' (Predicate_Compile_Aux.dest_indprem indprem), constant_table'))
bulwahn@39461
   446
    fun mk_clause pred_name pol (ts, prems) (prog, constant_table) =
bulwahn@39461
   447
    let
bulwahn@39461
   448
      val constant_table' = declare_consts (fold Term.add_const_names ts []) constant_table
bulwahn@39461
   449
      val args = map (translate_term ctxt constant_table') ts
bulwahn@39461
   450
      val (prems', constant_table'') = fold_map (mk_prem pol) prems constant_table'
bulwahn@39461
   451
    in
bulwahn@39461
   452
      (((pred_name, args), Conj prems') :: prog, constant_table'')
bulwahn@39461
   453
    end
bulwahn@39461
   454
    fun mk_clauses (pred, mode as (pol, _)) =
bulwahn@39461
   455
      let
bulwahn@39461
   456
        val clauses = Mode_Graph.get_node moded_gr (pred, mode)
bulwahn@39461
   457
        val pred_name = the (AList.lookup (op =) moded_pred_table' (pred, mode))
bulwahn@39461
   458
      in
bulwahn@39461
   459
        fold (mk_clause pred_name pol) clauses
bulwahn@39461
   460
      end
bulwahn@39461
   461
  in
bulwahn@39461
   462
    apsnd (pair moded_pred_table') (fold mk_clauses moded_preds (prog, constant_table))
bulwahn@39461
   463
  end
bulwahn@39461
   464
bulwahn@39461
   465
fun generate (use_modes, ensure_groundness) ctxt const =
bulwahn@38073
   466
  let 
bulwahn@38731
   467
    fun strong_conn_of gr keys =
bulwahn@38073
   468
      Graph.strong_conn (Graph.subgraph (member (op =) (Graph.all_succs gr keys)) gr)
bulwahn@38732
   469
    val gr = Predicate_Compile_Core.intros_graph_of ctxt
bulwahn@38731
   470
    val gr' = add_edges depending_preds_of const gr
bulwahn@38731
   471
    val scc = strong_conn_of gr' [const]
bulwahn@39461
   472
    val initial_constant_table = 
bulwahn@39461
   473
      declare_consts [@{const_name "Groups.zero_class.zero"}, @{const_name "Suc"}] []
bulwahn@38073
   474
  in
bulwahn@39461
   475
    case use_modes of
bulwahn@39461
   476
      SOME mode =>
bulwahn@39461
   477
        let
bulwahn@39461
   478
          val moded_gr = mk_moded_clauses_graph ctxt scc gr
bulwahn@39461
   479
          val moded_gr' = Mode_Graph.subgraph
bulwahn@39461
   480
            (member (op =) (Mode_Graph.all_succs moded_gr [(const, (true, mode))])) moded_gr
bulwahn@39461
   481
          val scc = Mode_Graph.strong_conn moded_gr' 
bulwahn@39461
   482
        in
bulwahn@39461
   483
          apfst rev (apsnd snd
bulwahn@39461
   484
            (fold (mk_program ctxt moded_gr') (rev scc) ([], ([], initial_constant_table))))
bulwahn@39461
   485
        end
bulwahn@39461
   486
      | NONE =>
bulwahn@39461
   487
        let 
bulwahn@39461
   488
          val _ = print_intros ctxt gr (flat scc)
bulwahn@39461
   489
          val constant_table = declare_consts (flat scc) initial_constant_table
bulwahn@39461
   490
        in
bulwahn@39461
   491
          apfst flat (fold_map (translate_intros ensure_groundness ctxt gr) (flat scc) constant_table)
bulwahn@39461
   492
        end
bulwahn@38073
   493
  end
bulwahn@38727
   494
  
bulwahn@38789
   495
(* implementation for fully enumerating predicates and
bulwahn@38789
   496
  for size-limited predicates for enumerating the values of a datatype upto a specific size *)
bulwahn@38073
   497
bulwahn@38727
   498
fun add_ground_typ (Conj prems) = fold add_ground_typ prems
bulwahn@38727
   499
  | add_ground_typ (Ground (_, T)) = insert (op =) T
bulwahn@38727
   500
  | add_ground_typ _ = I
bulwahn@38073
   501
bulwahn@38728
   502
fun mk_relname (Type (Tcon, Targs)) =
bulwahn@38728
   503
  first_lower (Long_Name.base_name Tcon) ^ space_implode "_" (map mk_relname Targs)
bulwahn@38728
   504
  | mk_relname _ = raise Fail "unexpected type"
bulwahn@38728
   505
bulwahn@38789
   506
fun mk_lim_relname T = "lim_" ^  mk_relname T
bulwahn@38789
   507
bulwahn@38728
   508
(* This is copied from "pat_completeness.ML" *)
bulwahn@38728
   509
fun inst_constrs_of thy (T as Type (name, _)) =
bulwahn@38728
   510
  map (fn (Cn,CT) =>
bulwahn@38728
   511
    Envir.subst_term_types (Sign.typ_match thy (body_type CT, T) Vartab.empty) (Const (Cn, CT)))
bulwahn@38728
   512
    (the (Datatype.get_constrs thy name))
bulwahn@38728
   513
  | inst_constrs_of thy T = raise TYPE ("inst_constrs_of", [T], [])
bulwahn@38789
   514
bulwahn@38789
   515
fun is_recursive_constr T (Const (constr_name, T')) = member (op =) (binder_types T') T
bulwahn@38728
   516
  
bulwahn@38789
   517
fun mk_ground_impl ctxt limited_types (T as Type (Tcon, Targs)) (seen, constant_table) =
bulwahn@38728
   518
  if member (op =) seen T then ([], (seen, constant_table))
bulwahn@38728
   519
  else
bulwahn@38728
   520
    let
bulwahn@38789
   521
      val (limited, size) = case AList.lookup (op =) limited_types T of
bulwahn@38789
   522
        SOME s => (true, s)
bulwahn@38789
   523
      | NONE => (false, 0)      
bulwahn@38789
   524
      val rel_name = (if limited then mk_lim_relname else mk_relname) T
bulwahn@38789
   525
      fun mk_impl (Const (constr_name, cT), recursive) (seen, constant_table) =
bulwahn@38727
   526
        let
bulwahn@38727
   527
          val constant_table' = declare_consts [constr_name] constant_table
bulwahn@38789
   528
          val Ts = binder_types cT
bulwahn@38728
   529
          val (rec_clauses, (seen', constant_table'')) =
bulwahn@38789
   530
            fold_map (mk_ground_impl ctxt limited_types) Ts (seen, constant_table')
bulwahn@38789
   531
          val vars = map (fn i => Var ("x" ^ string_of_int i)) (1 upto (length Ts))
bulwahn@38789
   532
          val lim_var =
bulwahn@38789
   533
            if limited then
bulwahn@38789
   534
              if recursive then [AppF ("suc", [Var "Lim"])]              
bulwahn@38789
   535
              else [Var "Lim"]
bulwahn@38789
   536
            else [] 
bulwahn@38789
   537
          fun mk_prem v T' =
bulwahn@38789
   538
            if limited andalso T' = T then Rel (mk_lim_relname T', [Var "Lim", v])
bulwahn@38789
   539
            else Rel (mk_relname T', [v])
bulwahn@38728
   540
          val clause =
bulwahn@38789
   541
            ((rel_name, lim_var @ [maybe_AppF (translate_const constant_table'' constr_name, vars)]),
bulwahn@38789
   542
             Conj (map2 mk_prem vars Ts))
bulwahn@38727
   543
        in
bulwahn@38728
   544
          (clause :: flat rec_clauses, (seen', constant_table''))
bulwahn@38727
   545
        end
bulwahn@38728
   546
      val constrs = inst_constrs_of (ProofContext.theory_of ctxt) T
bulwahn@38789
   547
      val constrs' = (constrs ~~ map (is_recursive_constr T) constrs)
bulwahn@38789
   548
        |> (fn cs => filter_out snd cs @ filter snd cs)
bulwahn@38789
   549
      val (clauses, constant_table') =
bulwahn@38789
   550
        apfst flat (fold_map mk_impl constrs' (T :: seen, constant_table))
bulwahn@38789
   551
      val size_term = funpow size (fn t => AppF ("suc", [t])) (Cons "zero")
bulwahn@38789
   552
    in
bulwahn@38789
   553
      ((if limited then
bulwahn@38789
   554
        cons ((mk_relname T, [Var "x"]), Rel (mk_lim_relname T, [size_term, Var "x"]))
bulwahn@38789
   555
      else I) clauses, constant_table')
bulwahn@38789
   556
    end
bulwahn@38789
   557
 | mk_ground_impl ctxt _ T (seen, constant_table) =
bulwahn@38728
   558
   raise Fail ("unexpected type :" ^ Syntax.string_of_typ ctxt T)
bulwahn@38728
   559
bulwahn@38727
   560
fun replace_ground (Conj prems) = Conj (map replace_ground prems)
bulwahn@38728
   561
  | replace_ground (Ground (x, T)) =
bulwahn@38728
   562
    Rel (mk_relname T, [Var x])  
bulwahn@38727
   563
  | replace_ground p = p
bulwahn@38727
   564
  
bulwahn@38789
   565
fun add_ground_predicates ctxt limited_types (p, constant_table) =
bulwahn@38727
   566
  let
bulwahn@38727
   567
    val ground_typs = fold (add_ground_typ o snd) p []
bulwahn@38789
   568
    val (grs, (_, constant_table')) = fold_map (mk_ground_impl ctxt limited_types) ground_typs ([], constant_table)
bulwahn@38727
   569
    val p' = map (apsnd replace_ground) p
bulwahn@38073
   570
  in
bulwahn@38727
   571
    ((flat grs) @ p', constant_table')
bulwahn@38073
   572
  end
bulwahn@38789
   573
bulwahn@38947
   574
(* make depth-limited version of predicate *)
bulwahn@38947
   575
bulwahn@38947
   576
fun mk_lim_rel_name rel_name = "lim_" ^ rel_name
bulwahn@38947
   577
bulwahn@38959
   578
fun mk_depth_limited rel_names ((rel_name, ts), prem) =
bulwahn@38947
   579
  let
bulwahn@38947
   580
    fun has_positive_recursive_prems (Conj prems) = exists has_positive_recursive_prems prems
bulwahn@38959
   581
      | has_positive_recursive_prems (Rel (rel, ts)) = member (op =) rel_names rel
bulwahn@38947
   582
      | has_positive_recursive_prems _ = false
bulwahn@38947
   583
    fun mk_lim_prem (Conj prems) = Conj (map mk_lim_prem prems)
bulwahn@38947
   584
      | mk_lim_prem (p as Rel (rel, ts)) =
bulwahn@38959
   585
        if member (op =) rel_names rel then Rel (mk_lim_rel_name rel, Var "Lim" :: ts) else p
bulwahn@38947
   586
      | mk_lim_prem p = p
bulwahn@38947
   587
  in
bulwahn@38947
   588
    if has_positive_recursive_prems prem then
bulwahn@38947
   589
      ((mk_lim_rel_name rel_name, (AppF ("suc", [Var "Lim"]))  :: ts), mk_lim_prem prem)
bulwahn@38947
   590
    else
bulwahn@38947
   591
      ((mk_lim_rel_name rel_name, (Var "Lim") :: ts), prem)
bulwahn@38947
   592
  end
bulwahn@38947
   593
bulwahn@38947
   594
fun add_limited_predicates limited_predicates =
bulwahn@38956
   595
  let                                     
bulwahn@38959
   596
    fun add (rel_names, limit) (p, constant_table) = 
bulwahn@38947
   597
      let
bulwahn@38959
   598
        val clauses = filter (fn ((rel, _), _) => member (op =) rel_names rel) p
bulwahn@38959
   599
        val clauses' = map (mk_depth_limited rel_names) clauses
bulwahn@38947
   600
        fun nat_term_of n = funpow n (fn t => AppF ("suc", [t])) (Cons "zero")
bulwahn@38959
   601
        fun mk_entry_clause rel_name =
bulwahn@38959
   602
          let
bulwahn@38959
   603
            val nargs = length (snd (fst
bulwahn@38959
   604
              (the (find_first (fn ((rel, _), _) => rel = rel_name) clauses))))
bulwahn@38959
   605
            val vars = map (fn i => Var ("x" ^ string_of_int i)) (1 upto nargs)        
bulwahn@38959
   606
          in
bulwahn@38959
   607
            (("limited_" ^ rel_name, vars), Rel ("lim_" ^ rel_name, nat_term_of limit :: vars))
bulwahn@38959
   608
          end
bulwahn@38959
   609
      in (p @ (map mk_entry_clause rel_names) @ clauses', constant_table) end
bulwahn@38947
   610
  in
bulwahn@38947
   611
    fold add limited_predicates
bulwahn@38947
   612
  end
bulwahn@38947
   613
bulwahn@38947
   614
bulwahn@38947
   615
(* replace predicates in clauses *)
bulwahn@38947
   616
bulwahn@38947
   617
(* replace (A, B, C) p = replace A by B in clauses of C *)
bulwahn@38947
   618
fun replace ((from, to), location) p =
bulwahn@38947
   619
  let
bulwahn@38947
   620
    fun replace_prem (Conj prems) = Conj (map replace_prem prems)
bulwahn@38947
   621
      | replace_prem (r as Rel (rel, ts)) =
bulwahn@38947
   622
          if rel = from then Rel (to, ts) else r
bulwahn@38947
   623
      | replace_prem r = r
bulwahn@38947
   624
  in
bulwahn@38947
   625
    map (fn ((rel, args), prem) => ((rel, args), (if rel = location then replace_prem else I) prem)) p
bulwahn@38947
   626
  end
bulwahn@38947
   627
bulwahn@38960
   628
  
bulwahn@38960
   629
(* reorder manually : reorder premises of ith clause of predicate p by a permutation perm *)
bulwahn@38947
   630
bulwahn@38960
   631
fun reorder_manually reorder p =
bulwahn@38960
   632
  let
bulwahn@38960
   633
    fun reorder' (clause as ((rel, args), prem)) seen =
bulwahn@38960
   634
      let
bulwahn@38960
   635
        val seen' = AList.map_default (op =) (rel, 0) (fn x => x + 1) seen
bulwahn@38960
   636
        val i = the (AList.lookup (op =) seen' rel)
bulwahn@38960
   637
        val perm = AList.lookup (op =) reorder (rel, i)
bulwahn@38960
   638
        val prem' = (case perm of 
bulwahn@38960
   639
          SOME p => (case prem of Conj prems => Conj (map (nth prems) p) | _ => prem)
bulwahn@38960
   640
        | NONE => prem)
bulwahn@38960
   641
      in (((rel, args), prem'), seen') end
bulwahn@38960
   642
  in
bulwahn@38960
   643
    fst (fold_map reorder' p [])
bulwahn@38960
   644
  end
bulwahn@39462
   645
bulwahn@38735
   646
(* rename variables to prolog-friendly names *)
bulwahn@38735
   647
bulwahn@38735
   648
fun rename_vars_term renaming = map_vars (fn v => the (AList.lookup (op =) renaming v))
bulwahn@38735
   649
bulwahn@38735
   650
fun rename_vars_prem renaming = map_term_prem (rename_vars_term renaming)
bulwahn@38735
   651
bulwahn@38735
   652
fun is_prolog_conform v =
bulwahn@38735
   653
  forall (fn s => Symbol.is_ascii_letter s orelse Symbol.is_ascii_digit s) (Symbol.explode v)
bulwahn@38735
   654
  
bulwahn@38735
   655
fun mk_renaming v renaming =
bulwahn@38958
   656
  (v, mk_conform first_upper "Var" (map snd renaming) v) :: renaming
bulwahn@38735
   657
bulwahn@38735
   658
fun rename_vars_clause ((rel, args), prem) =
bulwahn@38735
   659
  let
bulwahn@38735
   660
    val vars = fold_prem_terms add_vars prem (fold add_vars args [])
bulwahn@38735
   661
    val renaming = fold mk_renaming vars []
bulwahn@38735
   662
  in ((rel, map (rename_vars_term renaming) args), rename_vars_prem renaming prem) end
bulwahn@38735
   663
  
bulwahn@38735
   664
val rename_vars_program = map rename_vars_clause
bulwahn@38956
   665
bulwahn@39542
   666
(* post processing of generated prolog program *)
bulwahn@39542
   667
bulwahn@39542
   668
fun post_process ctxt options (p, constant_table) =
bulwahn@39542
   669
  (p, constant_table)
bulwahn@39542
   670
  |> (if #ensure_groundness options then
bulwahn@39542
   671
        add_ground_predicates ctxt (#limited_types options)
bulwahn@39542
   672
      else I)
bulwahn@39724
   673
  |> tap (fn _ => tracing "Adding limited predicates...")
bulwahn@39542
   674
  |> add_limited_predicates (#limited_predicates options)
bulwahn@39724
   675
  |> tap (fn _ => tracing "Replacing predicates...")
bulwahn@39542
   676
  |> apfst (fold replace (#replacing options))
bulwahn@39542
   677
  |> apfst (reorder_manually (#manual_reorder options))
bulwahn@39542
   678
  |> apfst rename_vars_program
bulwahn@39542
   679
bulwahn@38073
   680
(* code printer *)
bulwahn@38073
   681
bulwahn@38113
   682
fun write_arith_op Plus = "+"
bulwahn@38113
   683
  | write_arith_op Minus = "-"
bulwahn@38113
   684
bulwahn@38735
   685
fun write_term (Var v) = v
bulwahn@38079
   686
  | write_term (Cons c) = c
bulwahn@38113
   687
  | write_term (AppF (f, args)) = f ^ "(" ^ space_implode ", " (map write_term args) ^ ")"
bulwahn@38113
   688
  | write_term (ArithOp (oper, [a1, a2])) = write_term a1 ^ " " ^ write_arith_op oper ^ " " ^ write_term a2
bulwahn@38112
   689
  | write_term (Number n) = string_of_int n
bulwahn@38073
   690
bulwahn@38073
   691
fun write_rel (pred, args) =
bulwahn@38073
   692
  pred ^ "(" ^ space_implode ", " (map write_term args) ^ ")" 
bulwahn@38073
   693
bulwahn@38073
   694
fun write_prem (Conj prems) = space_implode ", " (map write_prem prems)
bulwahn@38073
   695
  | write_prem (Rel p) = write_rel p  
bulwahn@38073
   696
  | write_prem (NotRel p) = "not(" ^ write_rel p ^ ")"
bulwahn@38073
   697
  | write_prem (Eq (l, r)) = write_term l ^ " = " ^ write_term r
bulwahn@38073
   698
  | write_prem (NotEq (l, r)) = write_term l ^ " \\= " ^ write_term r
bulwahn@38113
   699
  | write_prem (ArithEq (l, r)) = write_term l ^ " is " ^ write_term r
bulwahn@38113
   700
  | write_prem (NotArithEq (l, r)) = write_term l ^ " =\\= " ^ write_term r
bulwahn@39461
   701
  | write_prem _ = raise Fail "Not a valid prolog premise"
bulwahn@38073
   702
bulwahn@38073
   703
fun write_clause (head, prem) =
bulwahn@38073
   704
  write_rel head ^ (if prem = Conj [] then "." else " :- " ^ write_prem prem ^ ".")
bulwahn@38073
   705
bulwahn@38073
   706
fun write_program p =
bulwahn@38073
   707
  cat_lines (map write_clause p) 
bulwahn@38073
   708
bulwahn@38790
   709
(* query templates *)
bulwahn@38078
   710
bulwahn@38792
   711
(** query and prelude for swi-prolog **)
bulwahn@38792
   712
bulwahn@39464
   713
fun swi_prolog_query_first (rel, args) vnames =
bulwahn@39464
   714
  "eval :- once("  ^ rel ^ "(" ^ space_implode ", " (map write_term args) ^ ")),\n" ^
bulwahn@38082
   715
  "writef('" ^ space_implode ";" (map (fn v => v ^ " = %w") vnames) ^
bulwahn@38082
   716
  "\\n', [" ^ space_implode ", " vnames ^ "]).\n"
bulwahn@38077
   717
  
bulwahn@39464
   718
fun swi_prolog_query_firstn n (rel, args) vnames =
bulwahn@38077
   719
  "eval :- findnsols(" ^ string_of_int n ^ ", (" ^ space_implode ", " vnames ^ "), " ^
bulwahn@39464
   720
    rel ^ "(" ^ space_implode ", " (map write_term args) ^ "), Sols), writelist(Sols).\n" ^
bulwahn@38077
   721
    "writelist([]).\n" ^
bulwahn@39546
   722
    "writelist([(" ^ space_implode ", " vnames ^ ")|SolutionTail]) :- " ^
bulwahn@38079
   723
    "writef('" ^ space_implode ";" (map (fn v => v ^ " = %w") vnames) ^
bulwahn@39546
   724
    "\\n', [" ^ space_implode ", " vnames ^ "]), writelist(SolutionTail).\n"
bulwahn@38077
   725
  
bulwahn@38792
   726
val swi_prolog_prelude =
bulwahn@38073
   727
  "#!/usr/bin/swipl -q -t main -f\n\n" ^
bulwahn@38077
   728
  ":- use_module(library('dialect/ciao/aggregates')).\n" ^
bulwahn@38729
   729
  ":- style_check(-singleton).\n" ^
bulwahn@38729
   730
  ":- style_check(-discontiguous).\n" ^ 	
bulwahn@38729
   731
  ":- style_check(-atom).\n\n" ^
bulwahn@38073
   732
  "main :- catch(eval, E, (print_message(error, E), fail)), halt.\n" ^
bulwahn@38073
   733
  "main :- halt(1).\n"
bulwahn@38075
   734
bulwahn@38792
   735
(** query and prelude for yap **)
bulwahn@38792
   736
bulwahn@39464
   737
fun yap_query_first (rel, args) vnames =
bulwahn@39464
   738
  "eval :- once(" ^ rel ^ "(" ^ space_implode ", " (map write_term args) ^ ")),\n" ^
bulwahn@38792
   739
  "format('" ^ space_implode ";" (map (fn v => v ^ " = ~w") vnames) ^
bulwahn@38792
   740
  "\\n', [" ^ space_implode ", " vnames ^ "]).\n"
bulwahn@38792
   741
bulwahn@38792
   742
val yap_prelude =
bulwahn@38792
   743
  "#!/usr/bin/yap -L\n\n" ^
bulwahn@38792
   744
  ":- initialization(eval).\n"
bulwahn@38792
   745
bulwahn@38792
   746
(* system-dependent query, prelude and invocation *)
bulwahn@38792
   747
bulwahn@38792
   748
fun query system nsols = 
bulwahn@38792
   749
  case system of
bulwahn@38792
   750
    SWI_PROLOG =>
bulwahn@38792
   751
      (case nsols of NONE => swi_prolog_query_first | SOME n => swi_prolog_query_firstn n)
bulwahn@38792
   752
  | YAP =>
bulwahn@38792
   753
      case nsols of NONE => yap_query_first | SOME n =>
bulwahn@38792
   754
        error "No support for querying multiple solutions in the prolog system yap"
bulwahn@38792
   755
bulwahn@38792
   756
fun prelude system =
bulwahn@38792
   757
  case system of SWI_PROLOG => swi_prolog_prelude | YAP => yap_prelude
bulwahn@38792
   758
bulwahn@38792
   759
fun invoke system file_name =
bulwahn@38792
   760
  let
bulwahn@39462
   761
    val env_var =
bulwahn@39462
   762
      (case system of SWI_PROLOG => "EXEC_SWIPL"| YAP => "EXEC_YAP")
bulwahn@39462
   763
    val prog = getenv env_var
bulwahn@38792
   764
    val cmd =
bulwahn@39462
   765
      case system of SWI_PROLOG => prog ^ " -f " | YAP => prog ^ " -L "
bulwahn@39462
   766
  in
bulwahn@39462
   767
    if prog = "" then
bulwahn@39462
   768
      (warning (env_var ^ " not set; could not execute code for " ^ string_of_system system); "")
bulwahn@39462
   769
    else fst (bash_output (cmd ^ file_name))
bulwahn@39462
   770
  end
bulwahn@38792
   771
bulwahn@38075
   772
(* parsing prolog solution *)
bulwahn@38790
   773
bulwahn@38115
   774
val scan_number =
bulwahn@38115
   775
  Scan.many1 Symbol.is_ascii_digit
bulwahn@38075
   776
bulwahn@38075
   777
val scan_atom =
bulwahn@38728
   778
  Scan.many1 (fn s => Symbol.is_ascii_lower s orelse Symbol.is_ascii_digit s orelse Symbol.is_ascii_quasi s)
bulwahn@38075
   779
bulwahn@38075
   780
val scan_var =
bulwahn@38078
   781
  Scan.many1
bulwahn@38078
   782
    (fn s => Symbol.is_ascii_upper s orelse Symbol.is_ascii_digit s orelse Symbol.is_ascii_quasi s)
bulwahn@38075
   783
bulwahn@38076
   784
val scan_ident =
bulwahn@38076
   785
  Scan.repeat (Scan.one
bulwahn@38076
   786
    (fn s => Symbol.is_ascii_letter s orelse Symbol.is_ascii_digit s orelse Symbol.is_ascii_quasi s))
bulwahn@38076
   787
bulwahn@38075
   788
fun dest_Char (Symbol.Char s) = s
bulwahn@38075
   789
bulwahn@38075
   790
val string_of = concat o map (dest_Char o Symbol.decode)
bulwahn@38075
   791
bulwahn@38076
   792
val is_atom_ident = forall Symbol.is_ascii_lower
bulwahn@38076
   793
bulwahn@38076
   794
val is_var_ident =
bulwahn@38076
   795
  forall (fn s => Symbol.is_ascii_upper s orelse Symbol.is_ascii_digit s orelse Symbol.is_ascii_quasi s)
bulwahn@38078
   796
bulwahn@38115
   797
fun int_of_symbol_list xs = fold (fn x => fn s => s * 10 + (ord x - ord "0")) xs 0
bulwahn@38115
   798
bulwahn@38078
   799
fun scan_terms xs = (((scan_term --| $$ ",") ::: scan_terms)
bulwahn@38078
   800
  || (scan_term >> single)) xs
bulwahn@38078
   801
and scan_term xs =
bulwahn@38115
   802
  ((scan_number >> (Number o int_of_symbol_list))
bulwahn@38115
   803
  || (scan_var >> (Var o string_of))
bulwahn@38078
   804
  || ((scan_atom -- ($$ "(" |-- scan_terms --| $$ ")"))
bulwahn@38079
   805
    >> (fn (f, ts) => AppF (string_of f, ts)))
bulwahn@38078
   806
  || (scan_atom >> (Cons o string_of))) xs
bulwahn@38079
   807
bulwahn@38075
   808
val parse_term = fst o Scan.finite Symbol.stopper
bulwahn@38077
   809
    (Scan.error (!! (fn _ => raise Fail "parsing prolog output failed")) scan_term)
bulwahn@38075
   810
  o explode
bulwahn@38075
   811
  
bulwahn@38079
   812
fun parse_solutions sol =
bulwahn@38075
   813
  let
bulwahn@38077
   814
    fun dest_eq s = case space_explode "=" s of
bulwahn@38075
   815
        (l :: r :: []) => parse_term (unprefix " " r)
bulwahn@38078
   816
      | _ => raise Fail "unexpected equation in prolog output"
bulwahn@38079
   817
    fun parse_solution s = map dest_eq (space_explode ";" s)
bulwahn@38961
   818
    val sols = case space_explode "\n" sol of [] => [] | s => fst (split_last s)  
bulwahn@38075
   819
  in
bulwahn@38961
   820
    map parse_solution sols
bulwahn@38075
   821
  end 
bulwahn@38073
   822
  
bulwahn@38073
   823
(* calling external interpreter and getting results *)
bulwahn@38073
   824
bulwahn@39464
   825
fun run (timeout, system) p (query_rel, args) vnames nsols =
bulwahn@38073
   826
  let
bulwahn@39464
   827
    val renaming = fold mk_renaming (fold add_vars args vnames) [] 
bulwahn@38735
   828
    val vnames' = map (fn v => the (AList.lookup (op =) renaming v)) vnames
bulwahn@39464
   829
    val args' = map (rename_vars_term renaming) args
bulwahn@39542
   830
    val prog = prelude system ^ query system nsols (query_rel, args') vnames' ^ write_program p
bulwahn@38079
   831
    val _ = tracing ("Generated prolog program:\n" ^ prog)
bulwahn@39187
   832
    val solution = TimeLimit.timeLimit timeout (fn prog => Cache_IO.with_tmp_file "prolog_file" (fn prolog_file =>
bulwahn@39187
   833
      (File.write prolog_file prog; invoke system (Path.implode prolog_file)))) prog
bulwahn@38079
   834
    val _ = tracing ("Prolog returned solution(s):\n" ^ solution)
bulwahn@38079
   835
    val tss = parse_solutions solution
bulwahn@38073
   836
  in
bulwahn@38079
   837
    tss
bulwahn@38073
   838
  end
bulwahn@38073
   839
bulwahn@38790
   840
(* restoring types in terms *)
bulwahn@38075
   841
bulwahn@38081
   842
fun restore_term ctxt constant_table (Var s, T) = Free (s, T)
bulwahn@38115
   843
  | restore_term ctxt constant_table (Number n, @{typ "int"}) = HOLogic.mk_number @{typ "int"} n
bulwahn@38115
   844
  | restore_term ctxt constant_table (Number n, _) = raise (Fail "unexpected type for number") 
bulwahn@38079
   845
  | restore_term ctxt constant_table (Cons s, T) = Const (restore_const constant_table s, T)
bulwahn@38079
   846
  | restore_term ctxt constant_table (AppF (f, args), T) =
bulwahn@38079
   847
    let
bulwahn@38079
   848
      val thy = ProofContext.theory_of ctxt
bulwahn@38079
   849
      val c = restore_const constant_table f
bulwahn@38079
   850
      val cT = Sign.the_const_type thy c
bulwahn@38079
   851
      val (argsT, resT) = strip_type cT
bulwahn@38079
   852
      val subst = Sign.typ_match thy (resT, T) Vartab.empty
bulwahn@38079
   853
      val argsT' = map (Envir.subst_type subst) argsT
bulwahn@38079
   854
    in
bulwahn@38079
   855
      list_comb (Const (c, Envir.subst_type subst cT),
bulwahn@38079
   856
        map (restore_term ctxt constant_table) (args ~~ argsT'))
bulwahn@38079
   857
    end
bulwahn@38079
   858
bulwahn@39465
   859
    
bulwahn@39465
   860
(* restore numerals in natural numbers *)
bulwahn@39465
   861
bulwahn@39465
   862
fun restore_nat_numerals t =
bulwahn@39465
   863
  if fastype_of t = @{typ nat} andalso is_some (try HOLogic.dest_nat t) then
bulwahn@39465
   864
    HOLogic.mk_number @{typ nat} (HOLogic.dest_nat t)
bulwahn@39465
   865
  else
bulwahn@39465
   866
    (case t of
bulwahn@39465
   867
        t1 $ t2 => restore_nat_numerals t1 $ restore_nat_numerals t2
bulwahn@39465
   868
      | t => t)
bulwahn@39465
   869
  
bulwahn@38790
   870
(* values command *)
bulwahn@38790
   871
bulwahn@38790
   872
val preprocess_options = Predicate_Compile_Aux.Options {
bulwahn@38790
   873
  expected_modes = NONE,
bulwahn@39383
   874
  proposed_modes = [],
bulwahn@38790
   875
  proposed_names = [],
bulwahn@38790
   876
  show_steps = false,
bulwahn@38790
   877
  show_intermediate_results = false,
bulwahn@38790
   878
  show_proof_trace = false,
bulwahn@38790
   879
  show_modes = false,
bulwahn@38790
   880
  show_mode_inference = false,
bulwahn@38790
   881
  show_compilation = false,
bulwahn@38790
   882
  show_caught_failures = false,
bulwahn@39383
   883
  show_invalid_clauses = false,
bulwahn@38790
   884
  skip_proof = true,
bulwahn@38790
   885
  no_topmost_reordering = false,
bulwahn@38790
   886
  function_flattening = true,
bulwahn@38790
   887
  specialise = false,
bulwahn@38790
   888
  fail_safe_function_flattening = false,
bulwahn@38790
   889
  no_higher_order_predicate = [],
bulwahn@38790
   890
  inductify = false,
bulwahn@38790
   891
  detect_switches = true,
bulwahn@38790
   892
  compilation = Predicate_Compile_Aux.Pred
bulwahn@38790
   893
}
bulwahn@38790
   894
bulwahn@38075
   895
fun values ctxt soln t_compr =
bulwahn@38075
   896
  let
bulwahn@38950
   897
    val options = code_options_of (ProofContext.theory_of ctxt)
bulwahn@38075
   898
    val split = case t_compr of (Const (@{const_name Collect}, _) $ t) => t
bulwahn@38075
   899
      | _ => error ("Not a set comprehension: " ^ Syntax.string_of_term ctxt t_compr);
bulwahn@38075
   900
    val (body, Ts, fp) = HOLogic.strip_psplits split;
bulwahn@38075
   901
    val output_names = Name.variant_list (Term.add_free_names body [])
bulwahn@38075
   902
      (map (fn i => "x" ^ string_of_int i) (1 upto length Ts))
bulwahn@38080
   903
    val output_frees = rev (map2 (curry Free) output_names Ts)
bulwahn@38075
   904
    val body = subst_bounds (output_frees, body)
bulwahn@38075
   905
    val (pred as Const (name, T), all_args) =
bulwahn@38075
   906
      case strip_comb body of
bulwahn@38075
   907
        (Const (name, T), all_args) => (Const (name, T), all_args)
bulwahn@38075
   908
      | (head, _) => error ("Not a constant: " ^ Syntax.string_of_term ctxt head)
bulwahn@38732
   909
    val _ = tracing "Preprocessing specification..."
bulwahn@38732
   910
    val T = Sign.the_const_type (ProofContext.theory_of ctxt) name
bulwahn@38732
   911
    val t = Const (name, T)
wenzelm@38755
   912
    val thy' =
wenzelm@38755
   913
      Theory.copy (ProofContext.theory_of ctxt)
wenzelm@38755
   914
      |> Predicate_Compile.preprocess preprocess_options t
wenzelm@38755
   915
    val ctxt' = ProofContext.init_global thy'
bulwahn@38079
   916
    val _ = tracing "Generating prolog program..."
bulwahn@39461
   917
    val (p, constant_table) = generate (NONE, #ensure_groundness options) ctxt' name (* FIXME *)
bulwahn@39542
   918
      |> post_process ctxt' options
bulwahn@39464
   919
    val constant_table' = declare_consts (fold Term.add_const_names all_args []) constant_table
bulwahn@39464
   920
    val args' = map (translate_term ctxt constant_table') all_args
bulwahn@38079
   921
    val _ = tracing "Running prolog program..."
bulwahn@39461
   922
    val system_config = System_Config.get (Context.Proof ctxt)
bulwahn@39461
   923
    val tss = run (#timeout system_config, #prolog_system system_config)
bulwahn@39464
   924
      p (translate_const constant_table' name, args') output_names soln
bulwahn@38079
   925
    val _ = tracing "Restoring terms..."
bulwahn@38115
   926
    val empty = Const("Orderings.bot_class.bot", fastype_of t_compr)
bulwahn@38115
   927
    fun mk_insert x S =
bulwahn@38115
   928
      Const (@{const_name "Set.insert"}, fastype_of x --> fastype_of S --> fastype_of S) $ x $ S 
bulwahn@38115
   929
    fun mk_set_compr in_insert [] xs =
bulwahn@38115
   930
       rev ((Free ("...", fastype_of t_compr)) ::
bulwahn@38115
   931
        (if null in_insert then xs else (fold mk_insert in_insert empty) :: xs))
bulwahn@38115
   932
      | mk_set_compr in_insert (t :: ts) xs =
bulwahn@38115
   933
        let
bulwahn@38115
   934
          val frees = Term.add_frees t []
bulwahn@38115
   935
        in
bulwahn@38115
   936
          if null frees then
bulwahn@38115
   937
            mk_set_compr (t :: in_insert) ts xs
bulwahn@38115
   938
          else
bulwahn@38115
   939
            let
wenzelm@38755
   940
              val uu as (uuN, uuT) = singleton (Variable.variant_frees ctxt' [t]) ("uu", fastype_of t)
bulwahn@38115
   941
              val set_compr =
bulwahn@38115
   942
                HOLogic.mk_Collect (uuN, uuT, fold (fn (s, T) => fn t => HOLogic.mk_exists (s, T, t))
bulwahn@38115
   943
                  frees (HOLogic.mk_conj (HOLogic.mk_eq (Free uu, t), @{term "True"})))
bulwahn@38115
   944
            in
bulwahn@38729
   945
              mk_set_compr [] ts
bulwahn@38729
   946
                (set_compr :: (if null in_insert then xs else (fold mk_insert in_insert empty) :: xs))  
bulwahn@38115
   947
            end
bulwahn@38115
   948
        end
bulwahn@38075
   949
  in
bulwahn@38115
   950
      foldl1 (HOLogic.mk_binop @{const_name sup}) (mk_set_compr []
bulwahn@39465
   951
        (map (fn ts => HOLogic.mk_tuple 
bulwahn@39465
   952
          (map (restore_nat_numerals o restore_term ctxt' constant_table) (ts ~~ Ts))) tss) [])
bulwahn@38075
   953
  end
bulwahn@38075
   954
bulwahn@38075
   955
fun values_cmd print_modes soln raw_t state =
bulwahn@38075
   956
  let
bulwahn@38075
   957
    val ctxt = Toplevel.context_of state
bulwahn@38075
   958
    val t = Syntax.read_term ctxt raw_t
bulwahn@38075
   959
    val t' = values ctxt soln t
bulwahn@38075
   960
    val ty' = Term.type_of t'
bulwahn@38075
   961
    val ctxt' = Variable.auto_fixes t' ctxt
bulwahn@38115
   962
    val _ = tracing "Printing terms..."
bulwahn@38075
   963
    val p = Print_Mode.with_modes print_modes (fn () =>
bulwahn@38075
   964
      Pretty.block [Pretty.quote (Syntax.pretty_term ctxt' t'), Pretty.fbrk,
bulwahn@38075
   965
        Pretty.str "::", Pretty.brk 1, Pretty.quote (Syntax.pretty_typ ctxt' ty')]) ();
bulwahn@38075
   966
  in Pretty.writeln p end;
bulwahn@38075
   967
bulwahn@38075
   968
bulwahn@38075
   969
(* renewing the values command for Prolog queries *)
bulwahn@38075
   970
bulwahn@38075
   971
val opt_print_modes =
bulwahn@38075
   972
  Scan.optional (Parse.$$$ "(" |-- Parse.!!! (Scan.repeat1 Parse.xname --| Parse.$$$ ")")) [];
bulwahn@38075
   973
bulwahn@38075
   974
val _ = Outer_Syntax.improper_command "values" "enumerate and print comprehensions" Keyword.diag
bulwahn@38077
   975
  (opt_print_modes -- Scan.optional (Parse.nat >> SOME) NONE -- Parse.term
bulwahn@38075
   976
   >> (fn ((print_modes, soln), t) => Toplevel.keep
haftmann@38504
   977
        (values_cmd print_modes soln t))); (*FIXME does not preserve the previous functionality*)
bulwahn@38075
   978
bulwahn@38733
   979
(* quickcheck generator *)
bulwahn@38733
   980
bulwahn@39541
   981
(* FIXME: a small clone of Predicate_Compile_Quickcheck - maybe refactor out commons *)
bulwahn@38733
   982
bulwahn@39461
   983
fun quickcheck ctxt t size =
bulwahn@38733
   984
  let
bulwahn@38950
   985
    val options = code_options_of (ProofContext.theory_of ctxt)
wenzelm@38755
   986
    val thy = Theory.copy (ProofContext.theory_of ctxt)
bulwahn@39541
   987
    val ((((full_constname, constT), vs'), intro), thy1) =
bulwahn@39541
   988
      Predicate_Compile_Aux.define_quickcheck_predicate t thy
bulwahn@38733
   989
    val thy2 = Context.theory_map (Predicate_Compile_Alternative_Defs.add_thm intro) thy1
bulwahn@39541
   990
    val thy3 = Predicate_Compile.preprocess preprocess_options (Const (full_constname, constT)) thy2
wenzelm@38755
   991
    val ctxt' = ProofContext.init_global thy3
bulwahn@38733
   992
    val _ = tracing "Generating prolog program..."
bulwahn@39461
   993
    val (p, constant_table) = generate (NONE, true) ctxt' full_constname
bulwahn@39542
   994
      |> post_process ctxt' (set_ensure_groundness options)
bulwahn@38733
   995
    val _ = tracing "Running prolog program..."
bulwahn@39461
   996
    val system_config = System_Config.get (Context.Proof ctxt)
bulwahn@39461
   997
    val tss = run (#timeout system_config, #prolog_system system_config)
bulwahn@39464
   998
      p (translate_const constant_table full_constname, map (Var o fst) vs') (map fst vs') (SOME 1)
bulwahn@38733
   999
    val _ = tracing "Restoring terms..."
bulwahn@38961
  1000
    val res =
bulwahn@38961
  1001
      case tss of
bulwahn@39541
  1002
        [ts] => SOME (map (restore_term ctxt' constant_table) (ts ~~ map snd vs'))
bulwahn@38961
  1003
      | _ => NONE
bulwahn@38733
  1004
    val empty_report = ([], false)
bulwahn@38733
  1005
  in
bulwahn@38733
  1006
    (res, empty_report)
bulwahn@38733
  1007
  end; 
bulwahn@38732
  1008
bulwahn@38073
  1009
end;