src/Pure/sorts.ML
author wenzelm
Fri Oct 07 21:16:48 2016 +0200 (2016-10-07)
changeset 64092 95469c544b82
parent 62538 85ebb645b1a3
child 68295 781a98696638
permissions -rw-r--r--
accept obscure timezone used in 2011;
wenzelm@2956
     1
(*  Title:      Pure/sorts.ML
wenzelm@2956
     2
    Author:     Markus Wenzel and Stefan Berghofer, TU Muenchen
wenzelm@2956
     3
wenzelm@19514
     4
The order-sorted algebra of type classes.
wenzelm@19529
     5
wenzelm@19529
     6
Classes denote (possibly empty) collections of types that are
wenzelm@19529
     7
partially ordered by class inclusion. They are represented
wenzelm@19529
     8
symbolically by strings.
wenzelm@19529
     9
wenzelm@19529
    10
Sorts are intersections of finitely many classes. They are represented
wenzelm@19529
    11
by lists of classes.  Normal forms of sorts are sorted lists of
wenzelm@19529
    12
minimal classes (wrt. current class inclusion).
wenzelm@2956
    13
*)
wenzelm@2956
    14
wenzelm@2956
    15
signature SORTS =
wenzelm@2956
    16
sig
wenzelm@39687
    17
  val make: sort list -> sort Ord_List.T
wenzelm@39687
    18
  val subset: sort Ord_List.T * sort Ord_List.T -> bool
wenzelm@39687
    19
  val union: sort Ord_List.T -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    20
  val subtract: sort Ord_List.T -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    21
  val remove_sort: sort -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    22
  val insert_sort: sort -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    23
  val insert_typ: typ -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    24
  val insert_typs: typ list -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    25
  val insert_term: term -> sort Ord_List.T -> sort Ord_List.T
wenzelm@39687
    26
  val insert_terms: term list -> sort Ord_List.T -> sort Ord_List.T
wenzelm@19645
    27
  type algebra
wenzelm@36328
    28
  val classes_of: algebra -> serial Graph.T
wenzelm@37248
    29
  val arities_of: algebra -> (class * sort list) list Symtab.table
wenzelm@21933
    30
  val all_classes: algebra -> class list
wenzelm@19645
    31
  val super_classes: algebra -> class -> class list
wenzelm@19645
    32
  val class_less: algebra -> class * class -> bool
wenzelm@19645
    33
  val class_le: algebra -> class * class -> bool
wenzelm@19645
    34
  val sort_eq: algebra -> sort * sort -> bool
wenzelm@19645
    35
  val sort_le: algebra -> sort * sort -> bool
wenzelm@19645
    36
  val sorts_le: algebra -> sort list * sort list -> bool
wenzelm@19645
    37
  val inter_sort: algebra -> sort * sort -> sort
wenzelm@24732
    38
  val minimize_sort: algebra -> sort -> sort
wenzelm@24732
    39
  val complete_sort: algebra -> sort -> sort
wenzelm@39687
    40
  val minimal_sorts: algebra -> sort list -> sort Ord_List.T
wenzelm@61262
    41
  val add_class: Context.generic -> class * class list -> algebra -> algebra
wenzelm@61262
    42
  val add_classrel: Context.generic -> class * class -> algebra -> algebra
wenzelm@61262
    43
  val add_arities: Context.generic -> string * (class * sort list) list -> algebra -> algebra
wenzelm@19645
    44
  val empty_algebra: algebra
wenzelm@61262
    45
  val merge_algebra: Context.generic -> algebra * algebra -> algebra
wenzelm@61262
    46
  val subalgebra: Context.generic -> (class -> bool) -> (class * string -> sort list option)
haftmann@22181
    47
    -> algebra -> (sort -> sort) * algebra
wenzelm@19578
    48
  type class_error
wenzelm@61262
    49
  val class_error: Context.generic -> class_error -> string
wenzelm@19578
    50
  exception CLASS_ERROR of class_error
wenzelm@48272
    51
  val has_instance: algebra -> string -> sort -> bool
wenzelm@19645
    52
  val mg_domain: algebra -> string -> sort -> sort list   (*exception CLASS_ERROR*)
haftmann@28665
    53
  val meet_sort: algebra -> typ * sort
haftmann@28665
    54
    -> sort Vartab.table -> sort Vartab.table   (*exception CLASS_ERROR*)
haftmann@28665
    55
  val meet_sort_typ: algebra -> typ * sort -> typ -> typ   (*exception CLASS_ERROR*)
wenzelm@19645
    56
  val of_sort: algebra -> typ * sort -> bool
wenzelm@32791
    57
  val of_sort_derivation: algebra ->
haftmann@62538
    58
    {class_relation: typ -> bool -> 'a * class -> class -> 'a,
wenzelm@36102
    59
     type_constructor: string * typ list -> ('a * class) list list -> class -> 'a,
wenzelm@22570
    60
     type_variable: typ -> ('a * class) list} ->
wenzelm@19584
    61
    typ * sort -> 'a list   (*exception CLASS_ERROR*)
wenzelm@35961
    62
  val classrel_derivation: algebra ->
wenzelm@35961
    63
    ('a * class -> class -> 'a) -> 'a * class -> class -> 'a  (*exception CLASS_ERROR*)
wenzelm@31946
    64
  val witness_sorts: algebra -> string list -> (typ * sort) list -> sort list -> (typ * sort) list
wenzelm@2956
    65
end;
wenzelm@2956
    66
wenzelm@20573
    67
structure Sorts: SORTS =
wenzelm@2956
    68
struct
wenzelm@2956
    69
wenzelm@19514
    70
wenzelm@19529
    71
(** ordered lists of sorts **)
wenzelm@14782
    72
wenzelm@39687
    73
val make = Ord_List.make Term_Ord.sort_ord;
wenzelm@39687
    74
val subset = Ord_List.subset Term_Ord.sort_ord;
wenzelm@39687
    75
val union = Ord_List.union Term_Ord.sort_ord;
wenzelm@39687
    76
val subtract = Ord_List.subtract Term_Ord.sort_ord;
wenzelm@14782
    77
wenzelm@39687
    78
val remove_sort = Ord_List.remove Term_Ord.sort_ord;
wenzelm@39687
    79
val insert_sort = Ord_List.insert Term_Ord.sort_ord;
wenzelm@14782
    80
wenzelm@16598
    81
fun insert_typ (TFree (_, S)) Ss = insert_sort S Ss
wenzelm@16598
    82
  | insert_typ (TVar (_, S)) Ss = insert_sort S Ss
wenzelm@16598
    83
  | insert_typ (Type (_, Ts)) Ss = insert_typs Ts Ss
wenzelm@16598
    84
and insert_typs [] Ss = Ss
wenzelm@16598
    85
  | insert_typs (T :: Ts) Ss = insert_typs Ts (insert_typ T Ss);
wenzelm@14782
    86
wenzelm@16598
    87
fun insert_term (Const (_, T)) Ss = insert_typ T Ss
wenzelm@16598
    88
  | insert_term (Free (_, T)) Ss = insert_typ T Ss
wenzelm@16598
    89
  | insert_term (Var (_, T)) Ss = insert_typ T Ss
wenzelm@16598
    90
  | insert_term (Bound _) Ss = Ss
wenzelm@16598
    91
  | insert_term (Abs (_, T, t)) Ss = insert_term t (insert_typ T Ss)
wenzelm@16598
    92
  | insert_term (t $ u) Ss = insert_term t (insert_term u Ss);
wenzelm@14782
    93
wenzelm@16598
    94
fun insert_terms [] Ss = Ss
wenzelm@16598
    95
  | insert_terms (t :: ts) Ss = insert_terms ts (insert_term t Ss);
wenzelm@14782
    96
wenzelm@14782
    97
wenzelm@19529
    98
wenzelm@19529
    99
(** order-sorted algebra **)
wenzelm@2956
   100
wenzelm@2956
   101
(*
wenzelm@14782
   102
  classes: graph representing class declarations together with proper
wenzelm@14782
   103
    subclass relation, which needs to be transitive and acyclic.
wenzelm@2956
   104
wenzelm@14782
   105
  arities: table of association lists of all type arities; (t, ars)
wenzelm@19531
   106
    means that type constructor t has the arities ars; an element
wenzelm@37248
   107
    (c, Ss) of ars represents the arity t::(Ss)c.  "Coregularity" of
wenzelm@37248
   108
    the arities structure requires that for any two declarations
wenzelm@37248
   109
    t::(Ss1)c1 and t::(Ss2)c2 such that c1 <= c2 holds Ss1 <= Ss2.
wenzelm@2956
   110
*)
wenzelm@2956
   111
wenzelm@19645
   112
datatype algebra = Algebra of
wenzelm@20573
   113
 {classes: serial Graph.T,
wenzelm@37248
   114
  arities: (class * sort list) list Symtab.table};
wenzelm@19645
   115
wenzelm@36328
   116
fun classes_of (Algebra {classes, ...}) = classes;
wenzelm@36328
   117
fun arities_of (Algebra {arities, ...}) = arities;
wenzelm@19645
   118
wenzelm@19645
   119
fun make_algebra (classes, arities) =
wenzelm@19645
   120
  Algebra {classes = classes, arities = arities};
wenzelm@19645
   121
wenzelm@19645
   122
fun map_classes f (Algebra {classes, arities}) = make_algebra (f classes, arities);
wenzelm@19645
   123
fun map_arities f (Algebra {classes, arities}) = make_algebra (classes, f arities);
wenzelm@19645
   124
wenzelm@19645
   125
wenzelm@19645
   126
(* classes *)
wenzelm@19645
   127
wenzelm@21933
   128
fun all_classes (Algebra {classes, ...}) = Graph.all_preds classes (Graph.maximals classes);
wenzelm@21933
   129
wenzelm@44338
   130
val super_classes = Graph.immediate_succs o classes_of;
wenzelm@2956
   131
wenzelm@2956
   132
wenzelm@19529
   133
(* class relations *)
wenzelm@2956
   134
wenzelm@19645
   135
val class_less = Graph.is_edge o classes_of;
wenzelm@19645
   136
fun class_le algebra (c1, c2) = c1 = c2 orelse class_less algebra (c1, c2);
wenzelm@2956
   137
wenzelm@2956
   138
wenzelm@19529
   139
(* sort relations *)
wenzelm@2956
   140
wenzelm@19645
   141
fun sort_le algebra (S1, S2) =
wenzelm@23585
   142
  S1 = S2 orelse forall (fn c2 => exists (fn c1 => class_le algebra (c1, c2)) S1) S2;
wenzelm@2956
   143
wenzelm@19645
   144
fun sorts_le algebra (Ss1, Ss2) =
wenzelm@19645
   145
  ListPair.all (sort_le algebra) (Ss1, Ss2);
wenzelm@2956
   146
wenzelm@19645
   147
fun sort_eq algebra (S1, S2) =
wenzelm@19645
   148
  sort_le algebra (S1, S2) andalso sort_le algebra (S2, S1);
wenzelm@2956
   149
wenzelm@2956
   150
wenzelm@19529
   151
(* intersection *)
wenzelm@2956
   152
wenzelm@19645
   153
fun inter_class algebra c S =
wenzelm@2956
   154
  let
wenzelm@2956
   155
    fun intr [] = [c]
wenzelm@2956
   156
      | intr (S' as c' :: c's) =
wenzelm@19645
   157
          if class_le algebra (c', c) then S'
wenzelm@19645
   158
          else if class_le algebra (c, c') then intr c's
wenzelm@2956
   159
          else c' :: intr c's
wenzelm@2956
   160
  in intr S end;
wenzelm@2956
   161
wenzelm@19645
   162
fun inter_sort algebra (S1, S2) =
wenzelm@19645
   163
  sort_strings (fold (inter_class algebra) S1 S2);
wenzelm@2956
   164
wenzelm@2956
   165
wenzelm@24732
   166
(* normal forms *)
wenzelm@2956
   167
wenzelm@24732
   168
fun minimize_sort _ [] = []
wenzelm@24732
   169
  | minimize_sort _ (S as [_]) = S
wenzelm@24732
   170
  | minimize_sort algebra S =
wenzelm@19645
   171
      filter (fn c => not (exists (fn c' => class_less algebra (c', c)) S)) S
wenzelm@19529
   172
      |> sort_distinct string_ord;
wenzelm@2990
   173
wenzelm@24732
   174
fun complete_sort algebra =
wenzelm@24732
   175
  Graph.all_succs (classes_of algebra) o minimize_sort algebra;
wenzelm@24732
   176
wenzelm@28623
   177
fun minimal_sorts algebra raw_sorts =
wenzelm@28623
   178
  let
wenzelm@28623
   179
    fun le S1 S2 = sort_le algebra (S1, S2);
wenzelm@28623
   180
    val sorts = make (map (minimize_sort algebra) raw_sorts);
wenzelm@28623
   181
  in sorts |> filter_out (fn S => exists (fn S' => le S' S andalso not (le S S')) sorts) end;
wenzelm@28623
   182
wenzelm@2990
   183
wenzelm@2956
   184
wenzelm@19529
   185
(** build algebras **)
wenzelm@19514
   186
wenzelm@19514
   187
(* classes *)
wenzelm@19514
   188
wenzelm@23655
   189
fun err_dup_class c = error ("Duplicate declaration of class: " ^ quote c);
wenzelm@19514
   190
wenzelm@61264
   191
fun err_cyclic_classes context css =
wenzelm@19514
   192
  error (cat_lines (map (fn cs =>
wenzelm@61264
   193
    "Cycle in class relation: " ^ Syntax.string_of_classrel (Syntax.init_pretty context) cs) css));
wenzelm@19514
   194
wenzelm@61264
   195
fun add_class context (c, cs) = map_classes (fn classes =>
wenzelm@19514
   196
  let
wenzelm@20573
   197
    val classes' = classes |> Graph.new_node (c, serial ())
wenzelm@23655
   198
      handle Graph.DUP dup => err_dup_class dup;
wenzelm@19514
   199
    val classes'' = classes' |> fold Graph.add_edge_trans_acyclic (map (pair c) cs)
wenzelm@61264
   200
      handle Graph.CYCLES css => err_cyclic_classes context css;
wenzelm@19645
   201
  in classes'' end);
wenzelm@19514
   202
wenzelm@19514
   203
wenzelm@19514
   204
(* arities *)
wenzelm@19514
   205
wenzelm@19514
   206
local
wenzelm@19514
   207
wenzelm@19514
   208
fun for_classes _ NONE = ""
wenzelm@42387
   209
  | for_classes ctxt (SOME (c1, c2)) = " for classes " ^ Syntax.string_of_classrel ctxt [c1, c2];
wenzelm@19514
   210
wenzelm@61264
   211
fun err_conflict context t cc (c, Ss) (c', Ss') =
wenzelm@61264
   212
  let val ctxt = Syntax.init_pretty context in
wenzelm@47005
   213
    error ("Conflict of type arities" ^ for_classes ctxt cc ^ ":\n  " ^
wenzelm@47005
   214
      Syntax.string_of_arity ctxt (t, Ss, [c]) ^ " and\n  " ^
wenzelm@47005
   215
      Syntax.string_of_arity ctxt (t, Ss', [c']))
wenzelm@47005
   216
  end;
wenzelm@19514
   217
wenzelm@61264
   218
fun coregular context algebra t (c, Ss) ars =
wenzelm@19514
   219
  let
wenzelm@37248
   220
    fun conflict (c', Ss') =
wenzelm@19645
   221
      if class_le algebra (c, c') andalso not (sorts_le algebra (Ss, Ss')) then
wenzelm@19514
   222
        SOME ((c, c'), (c', Ss'))
wenzelm@19645
   223
      else if class_le algebra (c', c) andalso not (sorts_le algebra (Ss', Ss)) then
wenzelm@19514
   224
        SOME ((c', c), (c', Ss'))
wenzelm@19514
   225
      else NONE;
wenzelm@19514
   226
  in
wenzelm@19514
   227
    (case get_first conflict ars of
wenzelm@61264
   228
      SOME ((c1, c2), (c', Ss')) => err_conflict context t (SOME (c1, c2)) (c, Ss) (c', Ss')
wenzelm@37248
   229
    | NONE => (c, Ss) :: ars)
wenzelm@19514
   230
  end;
wenzelm@19514
   231
wenzelm@37248
   232
fun complete algebra (c, Ss) = map (rpair Ss) (c :: super_classes algebra c);
wenzelm@19645
   233
wenzelm@61264
   234
fun insert context algebra t (c, Ss) ars =
wenzelm@19514
   235
  (case AList.lookup (op =) ars c of
wenzelm@61264
   236
    NONE => coregular context algebra t (c, Ss) ars
wenzelm@37248
   237
  | SOME Ss' =>
wenzelm@19645
   238
      if sorts_le algebra (Ss, Ss') then ars
wenzelm@37248
   239
      else if sorts_le algebra (Ss', Ss)
wenzelm@61264
   240
      then coregular context algebra t (c, Ss) (remove (op =) (c, Ss') ars)
wenzelm@61264
   241
      else err_conflict context t NONE (c, Ss) (c, Ss'));
wenzelm@19514
   242
wenzelm@35975
   243
in
wenzelm@35975
   244
wenzelm@61264
   245
fun insert_ars context algebra t = fold_rev (insert context algebra t);
wenzelm@35975
   246
wenzelm@61264
   247
fun insert_complete_ars context algebra (t, ars) arities =
wenzelm@19645
   248
  let val ars' =
wenzelm@19645
   249
    Symtab.lookup_list arities t
wenzelm@61264
   250
    |> fold_rev (insert_ars context algebra t) (map (complete algebra) ars);
wenzelm@19645
   251
  in Symtab.update (t, ars') arities end;
wenzelm@19514
   252
wenzelm@61264
   253
fun add_arities context arg algebra =
wenzelm@61264
   254
  algebra |> map_arities (insert_complete_ars context algebra arg);
wenzelm@19514
   255
wenzelm@61264
   256
fun add_arities_table context algebra =
wenzelm@61264
   257
  Symtab.fold (fn (t, ars) => insert_complete_ars context algebra (t, ars));
wenzelm@19514
   258
wenzelm@19514
   259
end;
wenzelm@19514
   260
wenzelm@19529
   261
wenzelm@19645
   262
(* classrel *)
wenzelm@19645
   263
wenzelm@61264
   264
fun rebuild_arities context algebra = algebra |> map_arities (fn arities =>
wenzelm@19645
   265
  Symtab.empty
wenzelm@61264
   266
  |> add_arities_table context algebra arities);
wenzelm@19645
   267
wenzelm@61264
   268
fun add_classrel context rel = rebuild_arities context o map_classes (fn classes =>
wenzelm@19645
   269
  classes |> Graph.add_edge_trans_acyclic rel
wenzelm@61264
   270
    handle Graph.CYCLES css => err_cyclic_classes context css);
wenzelm@19645
   271
wenzelm@19645
   272
wenzelm@19645
   273
(* empty and merge *)
wenzelm@19645
   274
wenzelm@19645
   275
val empty_algebra = make_algebra (Graph.empty, Symtab.empty);
wenzelm@19645
   276
wenzelm@61264
   277
fun merge_algebra context
wenzelm@19645
   278
   (Algebra {classes = classes1, arities = arities1},
wenzelm@19645
   279
    Algebra {classes = classes2, arities = arities2}) =
wenzelm@19645
   280
  let
wenzelm@19645
   281
    val classes' = Graph.merge_trans_acyclic (op =) (classes1, classes2)
wenzelm@23655
   282
      handle Graph.DUP c => err_dup_class c
wenzelm@61264
   283
        | Graph.CYCLES css => err_cyclic_classes context css;
wenzelm@19645
   284
    val algebra0 = make_algebra (classes', Symtab.empty);
wenzelm@35975
   285
    val arities' =
wenzelm@35975
   286
      (case (pointer_eq (classes1, classes2), pointer_eq (arities1, arities2)) of
wenzelm@35975
   287
        (true, true) => arities1
wenzelm@35975
   288
      | (true, false) =>  (*no completion*)
wenzelm@35975
   289
          (arities1, arities2) |> Symtab.join (fn t => fn (ars1, ars2) =>
wenzelm@35975
   290
            if pointer_eq (ars1, ars2) then raise Symtab.SAME
wenzelm@61264
   291
            else insert_ars context algebra0 t ars2 ars1)
wenzelm@35975
   292
      | (false, true) =>  (*unary completion*)
wenzelm@35975
   293
          Symtab.empty
wenzelm@61264
   294
          |> add_arities_table context algebra0 arities1
wenzelm@35975
   295
      | (false, false) => (*binary completion*)
wenzelm@35975
   296
          Symtab.empty
wenzelm@61264
   297
          |> add_arities_table context algebra0 arities1
wenzelm@61264
   298
          |> add_arities_table context algebra0 arities2);
wenzelm@19645
   299
  in make_algebra (classes', arities') end;
wenzelm@19645
   300
wenzelm@21933
   301
wenzelm@37248
   302
(* algebra projections *)  (* FIXME potentially violates abstract type integrity *)
haftmann@28922
   303
wenzelm@61264
   304
fun subalgebra context P sargs (algebra as Algebra {classes, arities}) =
haftmann@19952
   305
  let
wenzelm@24732
   306
    val restrict_sort = minimize_sort algebra o filter P o Graph.all_succs classes;
wenzelm@37248
   307
    fun restrict_arity t (c, Ss) =
wenzelm@37248
   308
      if P c then
wenzelm@37248
   309
        (case sargs (c, t) of
wenzelm@37248
   310
          SOME sorts =>
wenzelm@37248
   311
            SOME (c, Ss |> map2 (curry (inter_sort algebra)) sorts |> map restrict_sort)
wenzelm@37248
   312
        | NONE => NONE)
haftmann@22181
   313
      else NONE;
wenzelm@46614
   314
    val classes' = classes |> Graph.restrict P;
haftmann@39020
   315
    val arities' = arities |> Symtab.map (map_filter o restrict_arity);
wenzelm@61264
   316
  in (restrict_sort, rebuild_arities context (make_algebra (classes', arities'))) end;
haftmann@20465
   317
wenzelm@19645
   318
wenzelm@19529
   319
wenzelm@19529
   320
(** sorts of types **)
wenzelm@19529
   321
wenzelm@35961
   322
(* errors -- performance tuning via delayed message composition *)
wenzelm@19578
   323
wenzelm@26639
   324
datatype class_error =
wenzelm@36105
   325
  No_Classrel of class * class |
wenzelm@36105
   326
  No_Arity of string * class |
wenzelm@36105
   327
  No_Subsort of sort * sort;
wenzelm@19529
   328
wenzelm@61264
   329
fun class_error context =
wenzelm@61264
   330
  let val ctxt = Syntax.init_pretty context in
wenzelm@47005
   331
    fn No_Classrel (c1, c2) => "No class relation " ^ Syntax.string_of_classrel ctxt [c1, c2]
wenzelm@47005
   332
     | No_Arity (a, c) => "No type arity " ^ Syntax.string_of_arity ctxt (a, [], [c])
wenzelm@47005
   333
     | No_Subsort (S1, S2) =>
wenzelm@47005
   334
        "Cannot derive subsort relation " ^
wenzelm@47005
   335
          Syntax.string_of_sort ctxt S1 ^ " < " ^ Syntax.string_of_sort ctxt S2
wenzelm@47005
   336
  end;
wenzelm@19529
   337
wenzelm@19578
   338
exception CLASS_ERROR of class_error;
wenzelm@19578
   339
wenzelm@19578
   340
wenzelm@48272
   341
(* instances *)
wenzelm@48272
   342
wenzelm@48272
   343
fun has_instance algebra a =
wenzelm@48272
   344
  forall (AList.defined (op =) (Symtab.lookup_list (arities_of algebra) a));
wenzelm@19529
   345
wenzelm@19645
   346
fun mg_domain algebra a S =
wenzelm@19529
   347
  let
wenzelm@48272
   348
    val ars = Symtab.lookup_list (arities_of algebra) a;
wenzelm@19529
   349
    fun dom c =
wenzelm@48272
   350
      (case AList.lookup (op =) ars c of
wenzelm@36105
   351
        NONE => raise CLASS_ERROR (No_Arity (a, c))
wenzelm@37248
   352
      | SOME Ss => Ss);
wenzelm@19645
   353
    fun dom_inter c Ss = ListPair.map (inter_sort algebra) (dom c, Ss);
wenzelm@19529
   354
  in
wenzelm@19529
   355
    (case S of
wenzelm@19529
   356
      [] => raise Fail "Unknown domain of empty intersection"
wenzelm@19529
   357
    | c :: cs => fold dom_inter cs (dom c))
wenzelm@19529
   358
  end;
wenzelm@19529
   359
wenzelm@19529
   360
wenzelm@26639
   361
(* meet_sort *)
wenzelm@26639
   362
wenzelm@26639
   363
fun meet_sort algebra =
wenzelm@26639
   364
  let
wenzelm@26639
   365
    fun inters S S' = inter_sort algebra (S, S');
wenzelm@26639
   366
    fun meet _ [] = I
wenzelm@26639
   367
      | meet (TFree (_, S)) S' =
wenzelm@26639
   368
          if sort_le algebra (S, S') then I
wenzelm@36105
   369
          else raise CLASS_ERROR (No_Subsort (S, S'))
wenzelm@26639
   370
      | meet (TVar (v, S)) S' =
wenzelm@26639
   371
          if sort_le algebra (S, S') then I
wenzelm@26639
   372
          else Vartab.map_default (v, S) (inters S')
wenzelm@26639
   373
      | meet (Type (a, Ts)) S = fold2 meet Ts (mg_domain algebra a S);
wenzelm@26639
   374
  in uncurry meet end;
wenzelm@26639
   375
haftmann@28665
   376
fun meet_sort_typ algebra (T, S) =
wenzelm@37248
   377
  let val tab = meet_sort algebra (T, S) Vartab.empty;
wenzelm@37248
   378
  in Term.map_type_tvar (fn (v, _) => TVar (v, (the o Vartab.lookup tab) v)) end;
haftmann@28665
   379
wenzelm@26639
   380
wenzelm@19529
   381
(* of_sort *)
wenzelm@19529
   382
wenzelm@19645
   383
fun of_sort algebra =
wenzelm@19529
   384
  let
wenzelm@19529
   385
    fun ofS (_, []) = true
wenzelm@19645
   386
      | ofS (TFree (_, S), S') = sort_le algebra (S, S')
wenzelm@19645
   387
      | ofS (TVar (_, S), S') = sort_le algebra (S, S')
wenzelm@19529
   388
      | ofS (Type (a, Ts), S) =
wenzelm@19645
   389
          let val Ss = mg_domain algebra a S in
wenzelm@19529
   390
            ListPair.all ofS (Ts, Ss)
wenzelm@19578
   391
          end handle CLASS_ERROR _ => false;
wenzelm@19529
   392
  in ofS end;
wenzelm@19529
   393
wenzelm@19529
   394
haftmann@27498
   395
(* animating derivations *)
haftmann@27498
   396
wenzelm@32791
   397
fun of_sort_derivation algebra {class_relation, type_constructor, type_variable} =
wenzelm@19529
   398
  let
wenzelm@27555
   399
    val arities = arities_of algebra;
wenzelm@19578
   400
wenzelm@36104
   401
    fun weaken T D1 S2 =
wenzelm@36104
   402
      let val S1 = map snd D1 in
wenzelm@36104
   403
        if S1 = S2 then map fst D1
wenzelm@36104
   404
        else
wenzelm@36104
   405
          S2 |> map (fn c2 =>
haftmann@62538
   406
            (case D1 |> filter (fn (_, c1) => class_le algebra (c1, c2)) of
haftmann@62538
   407
              [d1] => class_relation T true d1 c2
haftmann@62538
   408
            | (d1 :: _ :: _) => class_relation T false d1 c2
haftmann@62538
   409
            | [] => raise CLASS_ERROR (No_Subsort (S1, S2))))
wenzelm@36104
   410
      end;
wenzelm@19529
   411
wenzelm@36103
   412
    fun derive (_, []) = []
wenzelm@44338
   413
      | derive (Type (a, Us), S) =
wenzelm@19529
   414
          let
wenzelm@19645
   415
            val Ss = mg_domain algebra a S;
wenzelm@36103
   416
            val dom = map2 (fn U => fn S => derive (U, S) ~~ S) Us Ss;
wenzelm@19529
   417
          in
wenzelm@19529
   418
            S |> map (fn c =>
wenzelm@19529
   419
              let
wenzelm@37248
   420
                val Ss' = the (AList.lookup (op =) (Symtab.lookup_list arities a) c);
wenzelm@36102
   421
                val dom' = map (fn ((U, d), S') => weaken U d S' ~~ S') ((Us ~~ dom) ~~ Ss');
wenzelm@37248
   422
              in type_constructor (a, Us) dom' c end)
wenzelm@19529
   423
          end
wenzelm@36103
   424
      | derive (T, S) = weaken T (type_variable T) S;
wenzelm@36103
   425
  in derive end;
wenzelm@19529
   426
wenzelm@35961
   427
fun classrel_derivation algebra class_relation =
wenzelm@35961
   428
  let
wenzelm@35961
   429
    fun path (x, c1 :: c2 :: cs) = path (class_relation (x, c1) c2, c2 :: cs)
wenzelm@35961
   430
      | path (x, _) = x;
wenzelm@35961
   431
  in
wenzelm@35961
   432
    fn (x, c1) => fn c2 =>
wenzelm@35961
   433
      (case Graph.irreducible_paths (classes_of algebra) (c1, c2) of
wenzelm@36105
   434
        [] => raise CLASS_ERROR (No_Classrel (c1, c2))
wenzelm@35961
   435
      | cs :: _ => path (x, cs))
wenzelm@35961
   436
  end;
wenzelm@35961
   437
wenzelm@19529
   438
wenzelm@19529
   439
(* witness_sorts *)
wenzelm@19529
   440
wenzelm@19645
   441
fun witness_sorts algebra types hyps sorts =
wenzelm@19529
   442
  let
wenzelm@19645
   443
    fun le S1 S2 = sort_le algebra (S1, S2);
wenzelm@31946
   444
    fun get S2 (T, S1) = if le S1 S2 then SOME (T, S2) else NONE;
wenzelm@19645
   445
    fun mg_dom t S = SOME (mg_domain algebra t S) handle CLASS_ERROR _ => NONE;
wenzelm@19529
   446
wenzelm@19578
   447
    fun witn_sort _ [] solved_failed = (SOME (propT, []), solved_failed)
wenzelm@19578
   448
      | witn_sort path S (solved, failed) =
wenzelm@19578
   449
          if exists (le S) failed then (NONE, (solved, failed))
wenzelm@19529
   450
          else
wenzelm@31946
   451
            (case get_first (get S) solved of
wenzelm@19578
   452
              SOME w => (SOME w, (solved, failed))
wenzelm@19529
   453
            | NONE =>
wenzelm@31946
   454
                (case get_first (get S) hyps of
wenzelm@19578
   455
                  SOME w => (SOME w, (w :: solved, failed))
wenzelm@19584
   456
                | NONE => witn_types path types S (solved, failed)))
wenzelm@19529
   457
wenzelm@19578
   458
    and witn_sorts path x = fold_map (witn_sort path) x
wenzelm@19529
   459
wenzelm@19578
   460
    and witn_types _ [] S (solved, failed) = (NONE, (solved, S :: failed))
wenzelm@19578
   461
      | witn_types path (t :: ts) S solved_failed =
wenzelm@19529
   462
          (case mg_dom t S of
wenzelm@19529
   463
            SOME SS =>
wenzelm@19529
   464
              (*do not descend into stronger args (achieving termination)*)
wenzelm@19529
   465
              if exists (fn D => le D S orelse exists (le D) path) SS then
wenzelm@19578
   466
                witn_types path ts S solved_failed
wenzelm@19529
   467
              else
wenzelm@19578
   468
                let val (ws, (solved', failed')) = witn_sorts (S :: path) SS solved_failed in
wenzelm@19529
   469
                  if forall is_some ws then
wenzelm@19529
   470
                    let val w = (Type (t, map (#1 o the) ws), S)
wenzelm@19578
   471
                    in (SOME w, (w :: solved', failed')) end
wenzelm@19578
   472
                  else witn_types path ts S (solved', failed')
wenzelm@19529
   473
                end
wenzelm@19578
   474
          | NONE => witn_types path ts S solved_failed);
wenzelm@19529
   475
wenzelm@19584
   476
  in map_filter I (#1 (witn_sorts [] sorts ([], []))) end;
wenzelm@19529
   477
wenzelm@19514
   478
end;