src/Pure/sorts.ML
author haftmann
Fri Dec 29 12:11:03 2006 +0100 (2006-12-29)
changeset 21926 1091904ddb19
parent 20573 c945a208e7f8
child 21933 819ef284720b
permissions -rw-r--r--
``classes`` now returns classes in topological order
wenzelm@2956
     1
(*  Title:      Pure/sorts.ML
wenzelm@2956
     2
    ID:         $Id$
wenzelm@2956
     3
    Author:     Markus Wenzel and Stefan Berghofer, TU Muenchen
wenzelm@2956
     4
wenzelm@19514
     5
The order-sorted algebra of type classes.
wenzelm@19529
     6
wenzelm@19529
     7
Classes denote (possibly empty) collections of types that are
wenzelm@19529
     8
partially ordered by class inclusion. They are represented
wenzelm@19529
     9
symbolically by strings.
wenzelm@19529
    10
wenzelm@19529
    11
Sorts are intersections of finitely many classes. They are represented
wenzelm@19529
    12
by lists of classes.  Normal forms of sorts are sorted lists of
wenzelm@19529
    13
minimal classes (wrt. current class inclusion).
wenzelm@2956
    14
*)
wenzelm@2956
    15
wenzelm@2956
    16
signature SORTS =
wenzelm@2956
    17
sig
wenzelm@16598
    18
  val eq_set: sort list * sort list -> bool
wenzelm@16598
    19
  val union: sort list -> sort list -> sort list
wenzelm@16598
    20
  val subtract: sort list -> sort list -> sort list
wenzelm@19463
    21
  val remove_sort: sort -> sort list -> sort list
wenzelm@16598
    22
  val insert_sort: sort -> sort list -> sort list
wenzelm@16598
    23
  val insert_typ: typ -> sort list -> sort list
wenzelm@16598
    24
  val insert_typs: typ list -> sort list -> sort list
wenzelm@16598
    25
  val insert_term: term -> sort list -> sort list
wenzelm@16598
    26
  val insert_terms: term list -> sort list -> sort list
wenzelm@19645
    27
  type algebra
wenzelm@19645
    28
  val rep_algebra: algebra ->
wenzelm@20573
    29
   {classes: serial Graph.T,
wenzelm@19645
    30
    arities: (class * (class * sort list)) list Symtab.table}
wenzelm@19645
    31
  val classes: algebra -> class list
wenzelm@19645
    32
  val super_classes: algebra -> class -> class list
wenzelm@19645
    33
  val class_less: algebra -> class * class -> bool
wenzelm@19645
    34
  val class_le: algebra -> class * class -> bool
wenzelm@19645
    35
  val sort_eq: algebra -> sort * sort -> bool
wenzelm@19645
    36
  val sort_le: algebra -> sort * sort -> bool
wenzelm@19645
    37
  val sorts_le: algebra -> sort list * sort list -> bool
wenzelm@19645
    38
  val inter_sort: algebra -> sort * sort -> sort
wenzelm@19645
    39
  val certify_class: algebra -> class -> class    (*exception TYPE*)
wenzelm@19645
    40
  val certify_sort: algebra -> sort -> sort       (*exception TYPE*)
wenzelm@19645
    41
  val add_class: Pretty.pp -> class * class list -> algebra -> algebra
wenzelm@19645
    42
  val add_classrel: Pretty.pp -> class * class -> algebra -> algebra
wenzelm@19645
    43
  val add_arities: Pretty.pp -> string * (class * sort list) list -> algebra -> algebra
wenzelm@19645
    44
  val empty_algebra: algebra
wenzelm@19645
    45
  val merge_algebra: Pretty.pp -> algebra * algebra -> algebra
haftmann@20465
    46
  val project_algebra: Pretty.pp -> (class -> bool) -> algebra -> (sort -> sort) * algebra
wenzelm@19578
    47
  type class_error
wenzelm@19578
    48
  val class_error: Pretty.pp -> class_error -> 'a
wenzelm@19578
    49
  exception CLASS_ERROR of class_error
wenzelm@19645
    50
  val mg_domain: algebra -> string -> sort -> sort list   (*exception CLASS_ERROR*)
wenzelm@19645
    51
  val of_sort: algebra -> typ * sort -> bool
wenzelm@19645
    52
  val of_sort_derivation: Pretty.pp -> algebra ->
wenzelm@19529
    53
    {classrel: 'a * class -> class -> 'a,
wenzelm@19529
    54
     constructor: string -> ('a * class) list list -> class -> 'a,
wenzelm@19584
    55
     variable: typ -> ('a * class) list} ->
wenzelm@19584
    56
    typ * sort -> 'a list   (*exception CLASS_ERROR*)
wenzelm@19645
    57
  val witness_sorts: algebra -> string list -> sort list -> sort list -> (typ * sort) list
wenzelm@2956
    58
end;
wenzelm@2956
    59
wenzelm@20573
    60
structure Sorts: SORTS =
wenzelm@2956
    61
struct
wenzelm@2956
    62
wenzelm@19514
    63
wenzelm@19529
    64
(** ordered lists of sorts **)
wenzelm@14782
    65
wenzelm@16598
    66
val eq_set = OrdList.eq_set Term.sort_ord;
wenzelm@16598
    67
val op union = OrdList.union Term.sort_ord;
wenzelm@16598
    68
val subtract = OrdList.subtract Term.sort_ord;
wenzelm@14782
    69
wenzelm@19463
    70
val remove_sort = OrdList.remove Term.sort_ord;
wenzelm@16598
    71
val insert_sort = OrdList.insert Term.sort_ord;
wenzelm@14782
    72
wenzelm@16598
    73
fun insert_typ (TFree (_, S)) Ss = insert_sort S Ss
wenzelm@16598
    74
  | insert_typ (TVar (_, S)) Ss = insert_sort S Ss
wenzelm@16598
    75
  | insert_typ (Type (_, Ts)) Ss = insert_typs Ts Ss
wenzelm@16598
    76
and insert_typs [] Ss = Ss
wenzelm@16598
    77
  | insert_typs (T :: Ts) Ss = insert_typs Ts (insert_typ T Ss);
wenzelm@14782
    78
wenzelm@16598
    79
fun insert_term (Const (_, T)) Ss = insert_typ T Ss
wenzelm@16598
    80
  | insert_term (Free (_, T)) Ss = insert_typ T Ss
wenzelm@16598
    81
  | insert_term (Var (_, T)) Ss = insert_typ T Ss
wenzelm@16598
    82
  | insert_term (Bound _) Ss = Ss
wenzelm@16598
    83
  | insert_term (Abs (_, T, t)) Ss = insert_term t (insert_typ T Ss)
wenzelm@16598
    84
  | insert_term (t $ u) Ss = insert_term t (insert_term u Ss);
wenzelm@14782
    85
wenzelm@16598
    86
fun insert_terms [] Ss = Ss
wenzelm@16598
    87
  | insert_terms (t :: ts) Ss = insert_terms ts (insert_term t Ss);
wenzelm@14782
    88
wenzelm@14782
    89
wenzelm@19529
    90
wenzelm@19529
    91
(** order-sorted algebra **)
wenzelm@2956
    92
wenzelm@2956
    93
(*
wenzelm@14782
    94
  classes: graph representing class declarations together with proper
wenzelm@14782
    95
    subclass relation, which needs to be transitive and acyclic.
wenzelm@2956
    96
wenzelm@14782
    97
  arities: table of association lists of all type arities; (t, ars)
wenzelm@19531
    98
    means that type constructor t has the arities ars; an element
wenzelm@19531
    99
    (c, (c0, Ss)) of ars represents the arity t::(Ss)c being derived
wenzelm@19531
   100
    via c0 <= c.  "Coregularity" of the arities structure requires
wenzelm@19531
   101
    that for any two declarations t::(Ss1)c1 and t::(Ss2)c2 such that
wenzelm@19531
   102
    c1 <= c2 holds Ss1 <= Ss2.
wenzelm@2956
   103
*)
wenzelm@2956
   104
wenzelm@19645
   105
datatype algebra = Algebra of
wenzelm@20573
   106
 {classes: serial Graph.T,
wenzelm@19645
   107
  arities: (class * (class * sort list)) list Symtab.table};
wenzelm@19645
   108
wenzelm@19645
   109
fun rep_algebra (Algebra args) = args;
wenzelm@19645
   110
wenzelm@19645
   111
val classes_of = #classes o rep_algebra;
wenzelm@19645
   112
val arities_of = #arities o rep_algebra;
wenzelm@19645
   113
wenzelm@19645
   114
fun make_algebra (classes, arities) =
wenzelm@19645
   115
  Algebra {classes = classes, arities = arities};
wenzelm@19645
   116
wenzelm@19645
   117
fun map_classes f (Algebra {classes, arities}) = make_algebra (f classes, arities);
wenzelm@19645
   118
fun map_arities f (Algebra {classes, arities}) = make_algebra (classes, f arities);
wenzelm@19645
   119
wenzelm@19645
   120
wenzelm@19645
   121
(* classes *)
wenzelm@19645
   122
haftmann@21926
   123
val classes = flat o rev o Graph.strong_conn o classes_of;
haftmann@21926
   124
  (*order allows for left-to-right traversal*)
wenzelm@19645
   125
val super_classes = Graph.imm_succs o classes_of;
wenzelm@2956
   126
wenzelm@2956
   127
wenzelm@19529
   128
(* class relations *)
wenzelm@2956
   129
wenzelm@19645
   130
val class_less = Graph.is_edge o classes_of;
wenzelm@19645
   131
fun class_le algebra (c1, c2) = c1 = c2 orelse class_less algebra (c1, c2);
wenzelm@2956
   132
wenzelm@2956
   133
wenzelm@19529
   134
(* sort relations *)
wenzelm@2956
   135
wenzelm@19645
   136
fun sort_le algebra (S1, S2) =
wenzelm@19645
   137
  forall (fn c2 => exists (fn c1 => class_le algebra (c1, c2)) S1) S2;
wenzelm@2956
   138
wenzelm@19645
   139
fun sorts_le algebra (Ss1, Ss2) =
wenzelm@19645
   140
  ListPair.all (sort_le algebra) (Ss1, Ss2);
wenzelm@2956
   141
wenzelm@19645
   142
fun sort_eq algebra (S1, S2) =
wenzelm@19645
   143
  sort_le algebra (S1, S2) andalso sort_le algebra (S2, S1);
wenzelm@2956
   144
wenzelm@2956
   145
wenzelm@19529
   146
(* intersection *)
wenzelm@2956
   147
wenzelm@19645
   148
fun inter_class algebra c S =
wenzelm@2956
   149
  let
wenzelm@2956
   150
    fun intr [] = [c]
wenzelm@2956
   151
      | intr (S' as c' :: c's) =
wenzelm@19645
   152
          if class_le algebra (c', c) then S'
wenzelm@19645
   153
          else if class_le algebra (c, c') then intr c's
wenzelm@2956
   154
          else c' :: intr c's
wenzelm@2956
   155
  in intr S end;
wenzelm@2956
   156
wenzelm@19645
   157
fun inter_sort algebra (S1, S2) =
wenzelm@19645
   158
  sort_strings (fold (inter_class algebra) S1 S2);
wenzelm@2956
   159
wenzelm@2956
   160
wenzelm@19645
   161
(* normal form *)
wenzelm@2956
   162
wenzelm@19529
   163
fun norm_sort _ [] = []
wenzelm@19529
   164
  | norm_sort _ (S as [_]) = S
wenzelm@19645
   165
  | norm_sort algebra S =
wenzelm@19645
   166
      filter (fn c => not (exists (fn c' => class_less algebra (c', c)) S)) S
wenzelm@19529
   167
      |> sort_distinct string_ord;
wenzelm@2990
   168
wenzelm@2990
   169
wenzelm@19645
   170
(* certify *)
wenzelm@19645
   171
wenzelm@19645
   172
fun certify_class algebra c =
wenzelm@19645
   173
  if can (Graph.get_node (classes_of algebra)) c then c
wenzelm@19645
   174
  else raise TYPE ("Undeclared class: " ^ quote c, [], []);
wenzelm@19645
   175
wenzelm@19645
   176
fun certify_sort classes = norm_sort classes o map (certify_class classes);
wenzelm@19645
   177
wenzelm@19645
   178
wenzelm@2956
   179
wenzelm@19529
   180
(** build algebras **)
wenzelm@19514
   181
wenzelm@19514
   182
(* classes *)
wenzelm@19514
   183
wenzelm@19514
   184
fun err_dup_classes cs =
wenzelm@19514
   185
  error ("Duplicate declaration of class(es): " ^ commas_quote cs);
wenzelm@19514
   186
wenzelm@19514
   187
fun err_cyclic_classes pp css =
wenzelm@19514
   188
  error (cat_lines (map (fn cs =>
wenzelm@19514
   189
    "Cycle in class relation: " ^ Pretty.string_of_classrel pp cs) css));
wenzelm@19514
   190
wenzelm@19645
   191
fun add_class pp (c, cs) = map_classes (fn classes =>
wenzelm@19514
   192
  let
wenzelm@20573
   193
    val classes' = classes |> Graph.new_node (c, serial ())
wenzelm@19514
   194
      handle Graph.DUP dup => err_dup_classes [dup];
wenzelm@19514
   195
    val classes'' = classes' |> fold Graph.add_edge_trans_acyclic (map (pair c) cs)
wenzelm@19514
   196
      handle Graph.CYCLES css => err_cyclic_classes pp css;
wenzelm@19645
   197
  in classes'' end);
wenzelm@19514
   198
wenzelm@19514
   199
wenzelm@19514
   200
(* arities *)
wenzelm@19514
   201
wenzelm@19514
   202
local
wenzelm@19514
   203
wenzelm@19514
   204
fun for_classes _ NONE = ""
wenzelm@19514
   205
  | for_classes pp (SOME (c1, c2)) =
wenzelm@19514
   206
      " for classes " ^ Pretty.string_of_classrel pp [c1, c2];
wenzelm@19514
   207
wenzelm@19514
   208
fun err_conflict pp t cc (c, Ss) (c', Ss') =
wenzelm@19514
   209
  error ("Conflict of type arities" ^ for_classes pp cc ^ ":\n  " ^
wenzelm@19514
   210
    Pretty.string_of_arity pp (t, Ss, [c]) ^ " and\n  " ^
wenzelm@19514
   211
    Pretty.string_of_arity pp (t, Ss', [c']));
wenzelm@19514
   212
wenzelm@19645
   213
fun coregular pp algebra t (c, (c0, Ss)) ars =
wenzelm@19514
   214
  let
wenzelm@19524
   215
    fun conflict (c', (_, Ss')) =
wenzelm@19645
   216
      if class_le algebra (c, c') andalso not (sorts_le algebra (Ss, Ss')) then
wenzelm@19514
   217
        SOME ((c, c'), (c', Ss'))
wenzelm@19645
   218
      else if class_le algebra (c', c) andalso not (sorts_le algebra (Ss', Ss)) then
wenzelm@19514
   219
        SOME ((c', c), (c', Ss'))
wenzelm@19514
   220
      else NONE;
wenzelm@19514
   221
  in
wenzelm@19514
   222
    (case get_first conflict ars of
wenzelm@19514
   223
      SOME ((c1, c2), (c', Ss')) => err_conflict pp t (SOME (c1, c2)) (c, Ss) (c', Ss')
wenzelm@19524
   224
    | NONE => (c, (c0, Ss)) :: ars)
wenzelm@19514
   225
  end;
wenzelm@19514
   226
wenzelm@19645
   227
fun complete algebra (c0, Ss) = map (rpair (c0, Ss)) (c0 :: super_classes algebra c0);
wenzelm@19645
   228
wenzelm@19645
   229
fun insert pp algebra t (c, (c0, Ss)) ars =
wenzelm@19514
   230
  (case AList.lookup (op =) ars c of
wenzelm@19645
   231
    NONE => coregular pp algebra t (c, (c0, Ss)) ars
wenzelm@19524
   232
  | SOME (_, Ss') =>
wenzelm@19645
   233
      if sorts_le algebra (Ss, Ss') then ars
wenzelm@19645
   234
      else if sorts_le algebra (Ss', Ss) then
wenzelm@19645
   235
        coregular pp algebra t (c, (c0, Ss))
wenzelm@19524
   236
          (filter_out (fn (c'', (_, Ss'')) => c = c'' andalso Ss'' = Ss') ars)
wenzelm@19514
   237
      else err_conflict pp t NONE (c, Ss) (c, Ss'));
wenzelm@19514
   238
wenzelm@19645
   239
fun insert_ars pp algebra (t, ars) arities =
wenzelm@19645
   240
  let val ars' =
wenzelm@19645
   241
    Symtab.lookup_list arities t
wenzelm@19645
   242
    |> fold_rev (fold_rev (insert pp algebra t)) (map (complete algebra) ars)
wenzelm@19645
   243
  in Symtab.update (t, ars') arities end;
wenzelm@19514
   244
wenzelm@19514
   245
in
wenzelm@19514
   246
wenzelm@19645
   247
fun add_arities pp arg algebra = algebra |> map_arities (insert_ars pp algebra arg);
wenzelm@19514
   248
wenzelm@19645
   249
fun add_arities_table pp algebra =
wenzelm@19645
   250
  Symtab.fold (fn (t, ars) => insert_ars pp algebra (t, map snd ars));
wenzelm@19514
   251
wenzelm@19514
   252
end;
wenzelm@19514
   253
wenzelm@19529
   254
wenzelm@19645
   255
(* classrel *)
wenzelm@19645
   256
wenzelm@19645
   257
fun rebuild_arities pp algebra = algebra |> map_arities (fn arities =>
wenzelm@19645
   258
  Symtab.empty
wenzelm@19645
   259
  |> add_arities_table pp algebra arities);
wenzelm@19645
   260
wenzelm@19645
   261
fun add_classrel pp rel = rebuild_arities pp o map_classes (fn classes =>
wenzelm@19645
   262
  classes |> Graph.add_edge_trans_acyclic rel
wenzelm@19645
   263
    handle Graph.CYCLES css => err_cyclic_classes pp css);
wenzelm@19645
   264
wenzelm@19645
   265
wenzelm@19645
   266
(* empty and merge *)
wenzelm@19645
   267
wenzelm@19645
   268
val empty_algebra = make_algebra (Graph.empty, Symtab.empty);
wenzelm@19645
   269
wenzelm@19645
   270
fun merge_algebra pp
wenzelm@19645
   271
   (Algebra {classes = classes1, arities = arities1},
wenzelm@19645
   272
    Algebra {classes = classes2, arities = arities2}) =
wenzelm@19645
   273
  let
wenzelm@19645
   274
    val classes' = Graph.merge_trans_acyclic (op =) (classes1, classes2)
wenzelm@19645
   275
      handle Graph.DUPS cs => err_dup_classes cs
wenzelm@19645
   276
          | Graph.CYCLES css => err_cyclic_classes pp css;
wenzelm@19645
   277
    val algebra0 = make_algebra (classes', Symtab.empty);
wenzelm@19645
   278
    val arities' = Symtab.empty
wenzelm@19645
   279
      |> add_arities_table pp algebra0 arities1
wenzelm@19645
   280
      |> add_arities_table pp algebra0 arities2;
wenzelm@19645
   281
  in make_algebra (classes', arities') end;
wenzelm@19645
   282
wenzelm@19977
   283
fun project_algebra pp proj (algebra as Algebra {classes, arities}) =
haftmann@19952
   284
  let
haftmann@20465
   285
    val proj_sort = norm_sort algebra o filter proj o Graph.all_succs classes;
haftmann@20465
   286
    fun proj_arity (c, (_, Ss)) =
haftmann@20465
   287
      if proj c then SOME (c, (c, map proj_sort Ss)) else NONE;
haftmann@20454
   288
    val classes' = classes |> Graph.project proj;
wenzelm@19977
   289
    val arities' = arities |> (Symtab.map o map_filter) proj_arity;
haftmann@20465
   290
  in (proj_sort, rebuild_arities pp (make_algebra (classes', arities'))) end;
haftmann@20465
   291
wenzelm@19645
   292
wenzelm@19529
   293
wenzelm@19529
   294
(** sorts of types **)
wenzelm@19529
   295
wenzelm@19578
   296
(* errors *)
wenzelm@19578
   297
wenzelm@19578
   298
datatype class_error = NoClassrel of class * class | NoArity of string * class;
wenzelm@19529
   299
wenzelm@19578
   300
fun class_error pp (NoClassrel (c1, c2)) =
wenzelm@19578
   301
      error ("No class relation " ^ Pretty.string_of_classrel pp [c1, c2])
wenzelm@19578
   302
  | class_error pp (NoArity (a, c)) =
wenzelm@19578
   303
      error ("No type arity " ^ Pretty.string_of_arity pp (a, [], [c]));
wenzelm@19529
   304
wenzelm@19578
   305
exception CLASS_ERROR of class_error;
wenzelm@19578
   306
wenzelm@19578
   307
wenzelm@19578
   308
(* mg_domain *)
wenzelm@19529
   309
wenzelm@19645
   310
fun mg_domain algebra a S =
wenzelm@19529
   311
  let
wenzelm@19645
   312
    val arities = arities_of algebra;
wenzelm@19529
   313
    fun dom c =
wenzelm@19529
   314
      (case AList.lookup (op =) (Symtab.lookup_list arities a) c of
wenzelm@19578
   315
        NONE => raise CLASS_ERROR (NoArity (a, c))
wenzelm@19529
   316
      | SOME (_, Ss) => Ss);
wenzelm@19645
   317
    fun dom_inter c Ss = ListPair.map (inter_sort algebra) (dom c, Ss);
wenzelm@19529
   318
  in
wenzelm@19529
   319
    (case S of
wenzelm@19529
   320
      [] => raise Fail "Unknown domain of empty intersection"
wenzelm@19529
   321
    | c :: cs => fold dom_inter cs (dom c))
wenzelm@19529
   322
  end;
wenzelm@19529
   323
wenzelm@19529
   324
wenzelm@19529
   325
(* of_sort *)
wenzelm@19529
   326
wenzelm@19645
   327
fun of_sort algebra =
wenzelm@19529
   328
  let
wenzelm@19529
   329
    fun ofS (_, []) = true
wenzelm@19645
   330
      | ofS (TFree (_, S), S') = sort_le algebra (S, S')
wenzelm@19645
   331
      | ofS (TVar (_, S), S') = sort_le algebra (S, S')
wenzelm@19529
   332
      | ofS (Type (a, Ts), S) =
wenzelm@19645
   333
          let val Ss = mg_domain algebra a S in
wenzelm@19529
   334
            ListPair.all ofS (Ts, Ss)
wenzelm@19578
   335
          end handle CLASS_ERROR _ => false;
wenzelm@19529
   336
  in ofS end;
wenzelm@19529
   337
wenzelm@19529
   338
wenzelm@19529
   339
(* of_sort_derivation *)
wenzelm@19529
   340
wenzelm@19645
   341
fun of_sort_derivation pp algebra {classrel, constructor, variable} =
wenzelm@19529
   342
  let
wenzelm@19645
   343
    val {classes, arities} = rep_algebra algebra;
haftmann@19952
   344
    fun weaken_path (x, c1 :: c2 :: cs) =
haftmann@19952
   345
          weaken_path (classrel (x, c1) c2, c2 :: cs)
wenzelm@19578
   346
      | weaken_path (x, _) = x;
wenzelm@19578
   347
    fun weaken (x, c1) c2 =
wenzelm@19578
   348
      (case Graph.irreducible_paths classes (c1, c2) of
wenzelm@19578
   349
        [] => raise CLASS_ERROR (NoClassrel (c1, c2))
wenzelm@19578
   350
      | cs :: _ => weaken_path (x, cs));
wenzelm@19578
   351
wenzelm@19529
   352
    fun weakens S1 S2 = S2 |> map (fn c2 =>
wenzelm@19645
   353
      (case S1 |> find_first (fn (_, c1) => class_le algebra (c1, c2)) of
wenzelm@19529
   354
        SOME d1 => weaken d1 c2
wenzelm@19529
   355
      | NONE => error ("Cannot derive subsort relation " ^
wenzelm@19529
   356
          Pretty.string_of_sort pp (map #2 S1) ^ " < " ^ Pretty.string_of_sort pp S2)));
wenzelm@19529
   357
wenzelm@19529
   358
    fun derive _ [] = []
wenzelm@19529
   359
      | derive (Type (a, Ts)) S =
wenzelm@19529
   360
          let
wenzelm@19645
   361
            val Ss = mg_domain algebra a S;
wenzelm@19529
   362
            val dom = map2 (fn T => fn S => derive T S ~~ S) Ts Ss;
wenzelm@19529
   363
          in
wenzelm@19529
   364
            S |> map (fn c =>
wenzelm@19529
   365
              let
wenzelm@19529
   366
                val (c0, Ss') = the (AList.lookup (op =) (Symtab.lookup_list arities a) c);
wenzelm@19529
   367
                val dom' = map2 (fn d => fn S' => weakens d S' ~~ S') dom Ss';
wenzelm@19529
   368
              in weaken (constructor a dom' c0, c0) c end)
wenzelm@19529
   369
          end
wenzelm@19529
   370
      | derive T S = weakens (variable T) S;
wenzelm@19529
   371
  in uncurry derive end;
wenzelm@19529
   372
wenzelm@19529
   373
wenzelm@19529
   374
(* witness_sorts *)
wenzelm@19529
   375
wenzelm@19645
   376
fun witness_sorts algebra types hyps sorts =
wenzelm@19529
   377
  let
wenzelm@19645
   378
    fun le S1 S2 = sort_le algebra (S1, S2);
wenzelm@19529
   379
    fun get_solved S2 (T, S1) = if le S1 S2 then SOME (T, S2) else NONE;
wenzelm@19529
   380
    fun get_hyp S2 S1 = if le S1 S2 then SOME (TFree ("'hyp", S1), S2) else NONE;
wenzelm@19645
   381
    fun mg_dom t S = SOME (mg_domain algebra t S) handle CLASS_ERROR _ => NONE;
wenzelm@19529
   382
wenzelm@19578
   383
    fun witn_sort _ [] solved_failed = (SOME (propT, []), solved_failed)
wenzelm@19578
   384
      | witn_sort path S (solved, failed) =
wenzelm@19578
   385
          if exists (le S) failed then (NONE, (solved, failed))
wenzelm@19529
   386
          else
wenzelm@19529
   387
            (case get_first (get_solved S) solved of
wenzelm@19578
   388
              SOME w => (SOME w, (solved, failed))
wenzelm@19529
   389
            | NONE =>
wenzelm@19529
   390
                (case get_first (get_hyp S) hyps of
wenzelm@19578
   391
                  SOME w => (SOME w, (w :: solved, failed))
wenzelm@19584
   392
                | NONE => witn_types path types S (solved, failed)))
wenzelm@19529
   393
wenzelm@19578
   394
    and witn_sorts path x = fold_map (witn_sort path) x
wenzelm@19529
   395
wenzelm@19578
   396
    and witn_types _ [] S (solved, failed) = (NONE, (solved, S :: failed))
wenzelm@19578
   397
      | witn_types path (t :: ts) S solved_failed =
wenzelm@19529
   398
          (case mg_dom t S of
wenzelm@19529
   399
            SOME SS =>
wenzelm@19529
   400
              (*do not descend into stronger args (achieving termination)*)
wenzelm@19529
   401
              if exists (fn D => le D S orelse exists (le D) path) SS then
wenzelm@19578
   402
                witn_types path ts S solved_failed
wenzelm@19529
   403
              else
wenzelm@19578
   404
                let val (ws, (solved', failed')) = witn_sorts (S :: path) SS solved_failed in
wenzelm@19529
   405
                  if forall is_some ws then
wenzelm@19529
   406
                    let val w = (Type (t, map (#1 o the) ws), S)
wenzelm@19578
   407
                    in (SOME w, (w :: solved', failed')) end
wenzelm@19578
   408
                  else witn_types path ts S (solved', failed')
wenzelm@19529
   409
                end
wenzelm@19578
   410
          | NONE => witn_types path ts S solved_failed);
wenzelm@19529
   411
wenzelm@19584
   412
  in map_filter I (#1 (witn_sorts [] sorts ([], []))) end;
wenzelm@19529
   413
wenzelm@19514
   414
end;