src/HOL/Code_Numeral.thy
author haftmann
Sun Oct 08 22:28:22 2017 +0200 (20 months ago)
changeset 66815 93c6632ddf44
parent 66806 a4e82b58d833
child 66817 0b12755ccbb2
permissions -rw-r--r--
one uniform type class for parity structures
haftmann@51143
     1
(*  Title:      HOL/Code_Numeral.thy
haftmann@51143
     2
    Author:     Florian Haftmann, TU Muenchen
haftmann@51143
     3
*)
haftmann@24999
     4
wenzelm@60758
     5
section \<open>Numeric types for code generation onto target language numerals only\<close>
haftmann@24999
     6
haftmann@31205
     7
theory Code_Numeral
haftmann@51143
     8
imports Nat_Transfer Divides Lifting
haftmann@51143
     9
begin
haftmann@51143
    10
wenzelm@60758
    11
subsection \<open>Type of target language integers\<close>
haftmann@51143
    12
wenzelm@61076
    13
typedef integer = "UNIV :: int set"
haftmann@51143
    14
  morphisms int_of_integer integer_of_int ..
haftmann@51143
    15
haftmann@59487
    16
setup_lifting type_definition_integer
haftmann@51143
    17
haftmann@51143
    18
lemma integer_eq_iff:
haftmann@51143
    19
  "k = l \<longleftrightarrow> int_of_integer k = int_of_integer l"
haftmann@51143
    20
  by transfer rule
haftmann@51143
    21
haftmann@51143
    22
lemma integer_eqI:
haftmann@51143
    23
  "int_of_integer k = int_of_integer l \<Longrightarrow> k = l"
haftmann@51143
    24
  using integer_eq_iff [of k l] by simp
haftmann@51143
    25
haftmann@51143
    26
lemma int_of_integer_integer_of_int [simp]:
haftmann@51143
    27
  "int_of_integer (integer_of_int k) = k"
haftmann@51143
    28
  by transfer rule
haftmann@51143
    29
haftmann@51143
    30
lemma integer_of_int_int_of_integer [simp]:
haftmann@51143
    31
  "integer_of_int (int_of_integer k) = k"
haftmann@51143
    32
  by transfer rule
haftmann@51143
    33
haftmann@51143
    34
instantiation integer :: ring_1
haftmann@24999
    35
begin
haftmann@24999
    36
haftmann@51143
    37
lift_definition zero_integer :: integer
haftmann@51143
    38
  is "0 :: int"
haftmann@51143
    39
  .
haftmann@51143
    40
haftmann@51143
    41
declare zero_integer.rep_eq [simp]
haftmann@24999
    42
haftmann@51143
    43
lift_definition one_integer :: integer
haftmann@51143
    44
  is "1 :: int"
haftmann@51143
    45
  .
haftmann@51143
    46
haftmann@51143
    47
declare one_integer.rep_eq [simp]
haftmann@24999
    48
haftmann@51143
    49
lift_definition plus_integer :: "integer \<Rightarrow> integer \<Rightarrow> integer"
haftmann@51143
    50
  is "plus :: int \<Rightarrow> int \<Rightarrow> int"
haftmann@51143
    51
  .
haftmann@51143
    52
haftmann@51143
    53
declare plus_integer.rep_eq [simp]
haftmann@24999
    54
haftmann@51143
    55
lift_definition uminus_integer :: "integer \<Rightarrow> integer"
haftmann@51143
    56
  is "uminus :: int \<Rightarrow> int"
haftmann@51143
    57
  .
haftmann@51143
    58
haftmann@51143
    59
declare uminus_integer.rep_eq [simp]
haftmann@24999
    60
haftmann@51143
    61
lift_definition minus_integer :: "integer \<Rightarrow> integer \<Rightarrow> integer"
haftmann@51143
    62
  is "minus :: int \<Rightarrow> int \<Rightarrow> int"
haftmann@51143
    63
  .
haftmann@51143
    64
haftmann@51143
    65
declare minus_integer.rep_eq [simp]
haftmann@24999
    66
haftmann@51143
    67
lift_definition times_integer :: "integer \<Rightarrow> integer \<Rightarrow> integer"
haftmann@51143
    68
  is "times :: int \<Rightarrow> int \<Rightarrow> int"
haftmann@51143
    69
  .
haftmann@51143
    70
haftmann@51143
    71
declare times_integer.rep_eq [simp]
haftmann@28708
    72
haftmann@51143
    73
instance proof
haftmann@51143
    74
qed (transfer, simp add: algebra_simps)+
haftmann@51143
    75
haftmann@51143
    76
end
haftmann@51143
    77
haftmann@64241
    78
instance integer :: Rings.dvd ..
haftmann@64241
    79
haftmann@64241
    80
lemma [transfer_rule]:
haftmann@64241
    81
  "rel_fun pcr_integer (rel_fun pcr_integer HOL.iff) Rings.dvd Rings.dvd"
haftmann@64241
    82
  unfolding dvd_def by transfer_prover
haftmann@64241
    83
haftmann@51143
    84
lemma [transfer_rule]:
blanchet@55945
    85
  "rel_fun HOL.eq pcr_integer (of_nat :: nat \<Rightarrow> int) (of_nat :: nat \<Rightarrow> integer)"
haftmann@64178
    86
  by (rule transfer_rule_of_nat) transfer_prover+
haftmann@51143
    87
haftmann@51143
    88
lemma [transfer_rule]:
blanchet@55945
    89
  "rel_fun HOL.eq pcr_integer (\<lambda>k :: int. k :: int) (of_int :: int \<Rightarrow> integer)"
haftmann@51143
    90
proof -
blanchet@55945
    91
  have "rel_fun HOL.eq pcr_integer (of_int :: int \<Rightarrow> int) (of_int :: int \<Rightarrow> integer)"
haftmann@64178
    92
    by (rule transfer_rule_of_int) transfer_prover+
haftmann@51143
    93
  then show ?thesis by (simp add: id_def)
haftmann@24999
    94
qed
haftmann@24999
    95
haftmann@51143
    96
lemma [transfer_rule]:
blanchet@55945
    97
  "rel_fun HOL.eq pcr_integer (numeral :: num \<Rightarrow> int) (numeral :: num \<Rightarrow> integer)"
haftmann@64178
    98
  by (rule transfer_rule_numeral) transfer_prover+
haftmann@26140
    99
haftmann@51143
   100
lemma [transfer_rule]:
blanchet@55945
   101
  "rel_fun HOL.eq (rel_fun HOL.eq pcr_integer) (Num.sub :: _ \<Rightarrow> _ \<Rightarrow> int) (Num.sub :: _ \<Rightarrow> _ \<Rightarrow> integer)"
haftmann@51143
   102
  by (unfold Num.sub_def [abs_def]) transfer_prover
haftmann@51143
   103
haftmann@51143
   104
lemma int_of_integer_of_nat [simp]:
haftmann@51143
   105
  "int_of_integer (of_nat n) = of_nat n"
haftmann@51143
   106
  by transfer rule
haftmann@51143
   107
haftmann@51143
   108
lift_definition integer_of_nat :: "nat \<Rightarrow> integer"
haftmann@51143
   109
  is "of_nat :: nat \<Rightarrow> int"
haftmann@51143
   110
  .
haftmann@51143
   111
haftmann@51143
   112
lemma integer_of_nat_eq_of_nat [code]:
haftmann@51143
   113
  "integer_of_nat = of_nat"
haftmann@51143
   114
  by transfer rule
haftmann@51143
   115
haftmann@51143
   116
lemma int_of_integer_integer_of_nat [simp]:
haftmann@51143
   117
  "int_of_integer (integer_of_nat n) = of_nat n"
haftmann@51143
   118
  by transfer rule
haftmann@51143
   119
haftmann@51143
   120
lift_definition nat_of_integer :: "integer \<Rightarrow> nat"
haftmann@51143
   121
  is Int.nat
haftmann@51143
   122
  .
haftmann@26140
   123
haftmann@51143
   124
lemma nat_of_integer_of_nat [simp]:
haftmann@51143
   125
  "nat_of_integer (of_nat n) = n"
haftmann@51143
   126
  by transfer simp
haftmann@51143
   127
haftmann@51143
   128
lemma int_of_integer_of_int [simp]:
haftmann@51143
   129
  "int_of_integer (of_int k) = k"
haftmann@51143
   130
  by transfer simp
haftmann@51143
   131
haftmann@51143
   132
lemma nat_of_integer_integer_of_nat [simp]:
haftmann@51143
   133
  "nat_of_integer (integer_of_nat n) = n"
haftmann@51143
   134
  by transfer simp
haftmann@51143
   135
haftmann@51143
   136
lemma integer_of_int_eq_of_int [simp, code_abbrev]:
haftmann@51143
   137
  "integer_of_int = of_int"
haftmann@51143
   138
  by transfer (simp add: fun_eq_iff)
haftmann@26140
   139
haftmann@51143
   140
lemma of_int_integer_of [simp]:
haftmann@51143
   141
  "of_int (int_of_integer k) = (k :: integer)"
haftmann@51143
   142
  by transfer rule
haftmann@51143
   143
haftmann@51143
   144
lemma int_of_integer_numeral [simp]:
haftmann@51143
   145
  "int_of_integer (numeral k) = numeral k"
haftmann@51143
   146
  by transfer rule
haftmann@51143
   147
haftmann@51143
   148
lemma int_of_integer_sub [simp]:
haftmann@51143
   149
  "int_of_integer (Num.sub k l) = Num.sub k l"
haftmann@51143
   150
  by transfer rule
haftmann@51143
   151
haftmann@66190
   152
definition integer_of_num :: "num \<Rightarrow> integer"
haftmann@66190
   153
  where [simp]: "integer_of_num = numeral"
haftmann@61275
   154
haftmann@61275
   155
lemma integer_of_num [code]:
haftmann@66190
   156
  "integer_of_num Num.One = 1"
haftmann@66190
   157
  "integer_of_num (Num.Bit0 n) = (let k = integer_of_num n in k + k)"
haftmann@66190
   158
  "integer_of_num (Num.Bit1 n) = (let k = integer_of_num n in k + k + 1)"
haftmann@66190
   159
  by (simp_all only: integer_of_num_def numeral.simps Let_def)
haftmann@61275
   160
haftmann@61275
   161
lemma integer_of_num_triv:
haftmann@61275
   162
  "integer_of_num Num.One = 1"
haftmann@61275
   163
  "integer_of_num (Num.Bit0 Num.One) = 2"
haftmann@66190
   164
  by simp_all
haftmann@61275
   165
haftmann@64592
   166
instantiation integer :: "{linordered_idom, equal}"
haftmann@26140
   167
begin
haftmann@26140
   168
haftmann@51143
   169
lift_definition abs_integer :: "integer \<Rightarrow> integer"
haftmann@51143
   170
  is "abs :: int \<Rightarrow> int"
haftmann@51143
   171
  .
haftmann@51143
   172
haftmann@51143
   173
declare abs_integer.rep_eq [simp]
haftmann@26140
   174
haftmann@51143
   175
lift_definition sgn_integer :: "integer \<Rightarrow> integer"
haftmann@51143
   176
  is "sgn :: int \<Rightarrow> int"
haftmann@51143
   177
  .
haftmann@51143
   178
haftmann@51143
   179
declare sgn_integer.rep_eq [simp]
haftmann@51143
   180
haftmann@51143
   181
lift_definition less_eq_integer :: "integer \<Rightarrow> integer \<Rightarrow> bool"
haftmann@51143
   182
  is "less_eq :: int \<Rightarrow> int \<Rightarrow> bool"
haftmann@51143
   183
  .
haftmann@51143
   184
haftmann@64592
   185
haftmann@51143
   186
lift_definition less_integer :: "integer \<Rightarrow> integer \<Rightarrow> bool"
haftmann@51143
   187
  is "less :: int \<Rightarrow> int \<Rightarrow> bool"
haftmann@51143
   188
  .
haftmann@51143
   189
haftmann@51143
   190
lift_definition equal_integer :: "integer \<Rightarrow> integer \<Rightarrow> bool"
haftmann@51143
   191
  is "HOL.equal :: int \<Rightarrow> int \<Rightarrow> bool"
haftmann@51143
   192
  .
haftmann@51143
   193
haftmann@64592
   194
instance
haftmann@64592
   195
  by standard (transfer, simp add: algebra_simps equal less_le_not_le [symmetric] mult_strict_right_mono linear)+
haftmann@26140
   196
haftmann@26140
   197
end
haftmann@26140
   198
haftmann@51143
   199
lemma [transfer_rule]:
blanchet@55945
   200
  "rel_fun pcr_integer (rel_fun pcr_integer pcr_integer) (min :: _ \<Rightarrow> _ \<Rightarrow> int) (min :: _ \<Rightarrow> _ \<Rightarrow> integer)"
haftmann@51143
   201
  by (unfold min_def [abs_def]) transfer_prover
haftmann@51143
   202
haftmann@51143
   203
lemma [transfer_rule]:
blanchet@55945
   204
  "rel_fun pcr_integer (rel_fun pcr_integer pcr_integer) (max :: _ \<Rightarrow> _ \<Rightarrow> int) (max :: _ \<Rightarrow> _ \<Rightarrow> integer)"
haftmann@51143
   205
  by (unfold max_def [abs_def]) transfer_prover
haftmann@51143
   206
haftmann@51143
   207
lemma int_of_integer_min [simp]:
haftmann@51143
   208
  "int_of_integer (min k l) = min (int_of_integer k) (int_of_integer l)"
haftmann@51143
   209
  by transfer rule
haftmann@51143
   210
haftmann@51143
   211
lemma int_of_integer_max [simp]:
haftmann@51143
   212
  "int_of_integer (max k l) = max (int_of_integer k) (int_of_integer l)"
haftmann@51143
   213
  by transfer rule
haftmann@26140
   214
haftmann@51143
   215
lemma nat_of_integer_non_positive [simp]:
haftmann@51143
   216
  "k \<le> 0 \<Longrightarrow> nat_of_integer k = 0"
haftmann@51143
   217
  by transfer simp
haftmann@51143
   218
haftmann@51143
   219
lemma of_nat_of_integer [simp]:
haftmann@51143
   220
  "of_nat (nat_of_integer k) = max 0 k"
haftmann@51143
   221
  by transfer auto
haftmann@51143
   222
haftmann@64848
   223
instantiation integer :: normalization_semidom
haftmann@64592
   224
begin
haftmann@64592
   225
haftmann@64592
   226
lift_definition normalize_integer :: "integer \<Rightarrow> integer"
haftmann@64592
   227
  is "normalize :: int \<Rightarrow> int"
haftmann@64592
   228
  .
haftmann@64592
   229
haftmann@64592
   230
declare normalize_integer.rep_eq [simp]
haftmann@64592
   231
haftmann@64592
   232
lift_definition unit_factor_integer :: "integer \<Rightarrow> integer"
haftmann@64592
   233
  is "unit_factor :: int \<Rightarrow> int"
haftmann@64592
   234
  .
haftmann@64592
   235
haftmann@64592
   236
declare unit_factor_integer.rep_eq [simp]
haftmann@64592
   237
haftmann@64592
   238
lift_definition divide_integer :: "integer \<Rightarrow> integer \<Rightarrow> integer"
haftmann@64592
   239
  is "divide :: int \<Rightarrow> int \<Rightarrow> int"
haftmann@64592
   240
  .
haftmann@64592
   241
haftmann@64592
   242
declare divide_integer.rep_eq [simp]
haftmann@64848
   243
  
haftmann@64848
   244
instance
haftmann@64848
   245
  by (standard; transfer)
haftmann@64848
   246
    (auto simp add: mult_sgn_abs sgn_mult abs_eq_iff')
haftmann@64592
   247
haftmann@64848
   248
end
haftmann@64848
   249
haftmann@66806
   250
instantiation integer :: unique_euclidean_ring
haftmann@64848
   251
begin
haftmann@64848
   252
  
haftmann@64592
   253
lift_definition modulo_integer :: "integer \<Rightarrow> integer \<Rightarrow> integer"
haftmann@64592
   254
  is "modulo :: int \<Rightarrow> int \<Rightarrow> int"
haftmann@64592
   255
  .
haftmann@64592
   256
haftmann@64592
   257
declare modulo_integer.rep_eq [simp]
haftmann@64592
   258
haftmann@66806
   259
lift_definition euclidean_size_integer :: "integer \<Rightarrow> nat"
haftmann@66806
   260
  is "euclidean_size :: int \<Rightarrow> nat"
haftmann@66806
   261
  .
haftmann@66806
   262
haftmann@66806
   263
declare euclidean_size_integer.rep_eq [simp]
haftmann@66806
   264
haftmann@66806
   265
lift_definition uniqueness_constraint_integer :: "integer \<Rightarrow> integer \<Rightarrow> bool"
haftmann@66806
   266
  is "uniqueness_constraint :: int \<Rightarrow> int \<Rightarrow> bool"
haftmann@66806
   267
  .
haftmann@66806
   268
haftmann@66806
   269
declare uniqueness_constraint_integer.rep_eq [simp]
haftmann@66806
   270
haftmann@64592
   271
instance
haftmann@66806
   272
  by (standard; transfer)
haftmann@66806
   273
    (use mult_le_mono2 [of 1] in \<open>auto simp add: sgn_mult_abs abs_mult sgn_mult abs_mod_less sgn_mod nat_mult_distrib\<close>, rule div_eqI, simp_all)
haftmann@64592
   274
haftmann@64592
   275
end
haftmann@64592
   276
haftmann@66806
   277
lemma [code]:
haftmann@66806
   278
  "euclidean_size = nat_of_integer \<circ> abs"
haftmann@66806
   279
  by (simp add: fun_eq_iff nat_of_integer.rep_eq)
haftmann@66806
   280
haftmann@66806
   281
lemma [code]:
haftmann@66806
   282
  "uniqueness_constraint (k :: integer) l \<longleftrightarrow> unit_factor k = unit_factor l"
haftmann@66806
   283
  by (simp add: integer_eq_iff)
haftmann@66806
   284
haftmann@66815
   285
instance integer :: ring_parity
haftmann@66815
   286
  by (standard; transfer) (simp_all add: of_nat_div odd_iff_mod_2_eq_one)
haftmann@66815
   287
haftmann@66806
   288
instantiation integer :: unique_euclidean_semiring_numeral
haftmann@61275
   289
begin
haftmann@61275
   290
haftmann@61275
   291
definition divmod_integer :: "num \<Rightarrow> num \<Rightarrow> integer \<times> integer"
haftmann@61275
   292
where
haftmann@61275
   293
  divmod_integer'_def: "divmod_integer m n = (numeral m div numeral n, numeral m mod numeral n)"
haftmann@61275
   294
haftmann@61275
   295
definition divmod_step_integer :: "num \<Rightarrow> integer \<times> integer \<Rightarrow> integer \<times> integer"
haftmann@61275
   296
where
haftmann@61275
   297
  "divmod_step_integer l qr = (let (q, r) = qr
haftmann@61275
   298
    in if r \<ge> numeral l then (2 * q + 1, r - numeral l)
haftmann@61275
   299
    else (2 * q, r))"
haftmann@61275
   300
haftmann@61275
   301
instance proof
haftmann@61275
   302
  show "divmod m n = (numeral m div numeral n :: integer, numeral m mod numeral n)"
haftmann@61275
   303
    for m n by (fact divmod_integer'_def)
haftmann@61275
   304
  show "divmod_step l qr = (let (q, r) = qr
haftmann@61275
   305
    in if r \<ge> numeral l then (2 * q + 1, r - numeral l)
haftmann@61275
   306
    else (2 * q, r))" for l and qr :: "integer \<times> integer"
haftmann@61275
   307
    by (fact divmod_step_integer_def)
haftmann@61275
   308
qed (transfer,
haftmann@61275
   309
  fact le_add_diff_inverse2
haftmann@66806
   310
  unique_euclidean_semiring_numeral_class.div_less
haftmann@66806
   311
  unique_euclidean_semiring_numeral_class.mod_less
haftmann@66806
   312
  unique_euclidean_semiring_numeral_class.div_positive
haftmann@66806
   313
  unique_euclidean_semiring_numeral_class.mod_less_eq_dividend
haftmann@66806
   314
  unique_euclidean_semiring_numeral_class.pos_mod_bound
haftmann@66806
   315
  unique_euclidean_semiring_numeral_class.pos_mod_sign
haftmann@66806
   316
  unique_euclidean_semiring_numeral_class.mod_mult2_eq
haftmann@66806
   317
  unique_euclidean_semiring_numeral_class.div_mult2_eq
haftmann@66806
   318
  unique_euclidean_semiring_numeral_class.discrete)+
haftmann@61275
   319
haftmann@61275
   320
end
haftmann@61275
   321
haftmann@61275
   322
declare divmod_algorithm_code [where ?'a = integer,
haftmann@66190
   323
  folded integer_of_num_def, unfolded integer_of_num_triv, 
haftmann@61275
   324
  code]
haftmann@53069
   325
Andreas@55427
   326
lemma integer_of_nat_0: "integer_of_nat 0 = 0"
Andreas@55427
   327
by transfer simp
Andreas@55427
   328
Andreas@55427
   329
lemma integer_of_nat_1: "integer_of_nat 1 = 1"
Andreas@55427
   330
by transfer simp
Andreas@55427
   331
Andreas@55427
   332
lemma integer_of_nat_numeral:
Andreas@55427
   333
  "integer_of_nat (numeral n) = numeral n"
Andreas@55427
   334
by transfer simp
haftmann@26140
   335
wenzelm@60758
   336
subsection \<open>Code theorems for target language integers\<close>
haftmann@51143
   337
wenzelm@60758
   338
text \<open>Constructors\<close>
haftmann@26140
   339
haftmann@51143
   340
definition Pos :: "num \<Rightarrow> integer"
haftmann@51143
   341
where
haftmann@61274
   342
  [simp, code_post]: "Pos = numeral"
haftmann@51143
   343
haftmann@51143
   344
lemma [transfer_rule]:
blanchet@55945
   345
  "rel_fun HOL.eq pcr_integer numeral Pos"
haftmann@51143
   346
  by simp transfer_prover
haftmann@30245
   347
haftmann@61274
   348
lemma Pos_fold [code_unfold]:
haftmann@61274
   349
  "numeral Num.One = Pos Num.One"
haftmann@61274
   350
  "numeral (Num.Bit0 k) = Pos (Num.Bit0 k)"
haftmann@61274
   351
  "numeral (Num.Bit1 k) = Pos (Num.Bit1 k)"
haftmann@61274
   352
  by simp_all
haftmann@61274
   353
haftmann@51143
   354
definition Neg :: "num \<Rightarrow> integer"
haftmann@51143
   355
where
haftmann@54489
   356
  [simp, code_abbrev]: "Neg n = - Pos n"
haftmann@51143
   357
haftmann@51143
   358
lemma [transfer_rule]:
blanchet@55945
   359
  "rel_fun HOL.eq pcr_integer (\<lambda>n. - numeral n) Neg"
haftmann@54489
   360
  by (simp add: Neg_def [abs_def]) transfer_prover
haftmann@51143
   361
haftmann@51143
   362
code_datatype "0::integer" Pos Neg
haftmann@51143
   363
haftmann@64994
   364
  
haftmann@64994
   365
text \<open>A further pair of constructors for generated computations\<close>
haftmann@64994
   366
haftmann@64994
   367
context
haftmann@64994
   368
begin  
haftmann@64994
   369
haftmann@64994
   370
qualified definition positive :: "num \<Rightarrow> integer"
haftmann@64994
   371
  where [simp]: "positive = numeral"
haftmann@64994
   372
haftmann@64994
   373
qualified definition negative :: "num \<Rightarrow> integer"
haftmann@64994
   374
  where [simp]: "negative = uminus \<circ> numeral"
haftmann@64994
   375
haftmann@64994
   376
lemma [code_computation_unfold]:
haftmann@64994
   377
  "numeral = positive"
haftmann@64994
   378
  "Pos = positive"
haftmann@64994
   379
  "Neg = negative"
haftmann@64994
   380
  by (simp_all add: fun_eq_iff)
haftmann@64994
   381
haftmann@64994
   382
end
haftmann@64994
   383
haftmann@51143
   384
wenzelm@60758
   385
text \<open>Auxiliary operations\<close>
haftmann@51143
   386
haftmann@51143
   387
lift_definition dup :: "integer \<Rightarrow> integer"
haftmann@51143
   388
  is "\<lambda>k::int. k + k"
haftmann@51143
   389
  .
haftmann@26140
   390
haftmann@51143
   391
lemma dup_code [code]:
haftmann@51143
   392
  "dup 0 = 0"
haftmann@51143
   393
  "dup (Pos n) = Pos (Num.Bit0 n)"
haftmann@51143
   394
  "dup (Neg n) = Neg (Num.Bit0 n)"
haftmann@54489
   395
  by (transfer, simp only: numeral_Bit0 minus_add_distrib)+
haftmann@51143
   396
haftmann@51143
   397
lift_definition sub :: "num \<Rightarrow> num \<Rightarrow> integer"
haftmann@51143
   398
  is "\<lambda>m n. numeral m - numeral n :: int"
haftmann@51143
   399
  .
haftmann@26140
   400
haftmann@51143
   401
lemma sub_code [code]:
haftmann@51143
   402
  "sub Num.One Num.One = 0"
haftmann@51143
   403
  "sub (Num.Bit0 m) Num.One = Pos (Num.BitM m)"
haftmann@51143
   404
  "sub (Num.Bit1 m) Num.One = Pos (Num.Bit0 m)"
haftmann@51143
   405
  "sub Num.One (Num.Bit0 n) = Neg (Num.BitM n)"
haftmann@51143
   406
  "sub Num.One (Num.Bit1 n) = Neg (Num.Bit0 n)"
haftmann@51143
   407
  "sub (Num.Bit0 m) (Num.Bit0 n) = dup (sub m n)"
haftmann@51143
   408
  "sub (Num.Bit1 m) (Num.Bit1 n) = dup (sub m n)"
haftmann@51143
   409
  "sub (Num.Bit1 m) (Num.Bit0 n) = dup (sub m n) + 1"
haftmann@51143
   410
  "sub (Num.Bit0 m) (Num.Bit1 n) = dup (sub m n) - 1"
haftmann@51143
   411
  by (transfer, simp add: dbl_def dbl_inc_def dbl_dec_def)+
haftmann@28351
   412
haftmann@24999
   413
wenzelm@60758
   414
text \<open>Implementations\<close>
haftmann@24999
   415
haftmann@51143
   416
lemma one_integer_code [code, code_unfold]:
haftmann@51143
   417
  "1 = Pos Num.One"
haftmann@51143
   418
  by simp
haftmann@24999
   419
haftmann@51143
   420
lemma plus_integer_code [code]:
haftmann@51143
   421
  "k + 0 = (k::integer)"
haftmann@51143
   422
  "0 + l = (l::integer)"
haftmann@51143
   423
  "Pos m + Pos n = Pos (m + n)"
haftmann@51143
   424
  "Pos m + Neg n = sub m n"
haftmann@51143
   425
  "Neg m + Pos n = sub n m"
haftmann@51143
   426
  "Neg m + Neg n = Neg (m + n)"
haftmann@51143
   427
  by (transfer, simp)+
haftmann@24999
   428
haftmann@51143
   429
lemma uminus_integer_code [code]:
haftmann@51143
   430
  "uminus 0 = (0::integer)"
haftmann@51143
   431
  "uminus (Pos m) = Neg m"
haftmann@51143
   432
  "uminus (Neg m) = Pos m"
haftmann@51143
   433
  by simp_all
haftmann@28708
   434
haftmann@51143
   435
lemma minus_integer_code [code]:
haftmann@51143
   436
  "k - 0 = (k::integer)"
haftmann@51143
   437
  "0 - l = uminus (l::integer)"
haftmann@51143
   438
  "Pos m - Pos n = sub m n"
haftmann@51143
   439
  "Pos m - Neg n = Pos (m + n)"
haftmann@51143
   440
  "Neg m - Pos n = Neg (m + n)"
haftmann@51143
   441
  "Neg m - Neg n = sub n m"
haftmann@51143
   442
  by (transfer, simp)+
haftmann@46028
   443
haftmann@51143
   444
lemma abs_integer_code [code]:
haftmann@51143
   445
  "\<bar>k\<bar> = (if (k::integer) < 0 then - k else k)"
haftmann@51143
   446
  by simp
huffman@47108
   447
haftmann@51143
   448
lemma sgn_integer_code [code]:
haftmann@51143
   449
  "sgn k = (if k = 0 then 0 else if (k::integer) < 0 then - 1 else 1)"
huffman@47108
   450
  by simp
haftmann@46028
   451
haftmann@51143
   452
lemma times_integer_code [code]:
haftmann@51143
   453
  "k * 0 = (0::integer)"
haftmann@51143
   454
  "0 * l = (0::integer)"
haftmann@51143
   455
  "Pos m * Pos n = Pos (m * n)"
haftmann@51143
   456
  "Pos m * Neg n = Neg (m * n)"
haftmann@51143
   457
  "Neg m * Pos n = Neg (m * n)"
haftmann@51143
   458
  "Neg m * Neg n = Pos (m * n)"
haftmann@51143
   459
  by simp_all
haftmann@51143
   460
haftmann@64592
   461
lemma normalize_integer_code [code]:
haftmann@64592
   462
  "normalize = (abs :: integer \<Rightarrow> integer)"
haftmann@64592
   463
  by transfer simp
haftmann@64592
   464
haftmann@64592
   465
lemma unit_factor_integer_code [code]:
haftmann@64592
   466
  "unit_factor = (sgn :: integer \<Rightarrow> integer)"
haftmann@64592
   467
  by transfer simp
haftmann@64592
   468
haftmann@51143
   469
definition divmod_integer :: "integer \<Rightarrow> integer \<Rightarrow> integer \<times> integer"
haftmann@51143
   470
where
haftmann@51143
   471
  "divmod_integer k l = (k div l, k mod l)"
haftmann@51143
   472
haftmann@66801
   473
lemma fst_divmod_integer [simp]:
haftmann@51143
   474
  "fst (divmod_integer k l) = k div l"
haftmann@51143
   475
  by (simp add: divmod_integer_def)
haftmann@51143
   476
haftmann@66801
   477
lemma snd_divmod_integer [simp]:
haftmann@51143
   478
  "snd (divmod_integer k l) = k mod l"
haftmann@51143
   479
  by (simp add: divmod_integer_def)
haftmann@51143
   480
haftmann@51143
   481
definition divmod_abs :: "integer \<Rightarrow> integer \<Rightarrow> integer \<times> integer"
haftmann@51143
   482
where
haftmann@51143
   483
  "divmod_abs k l = (\<bar>k\<bar> div \<bar>l\<bar>, \<bar>k\<bar> mod \<bar>l\<bar>)"
haftmann@51143
   484
haftmann@51143
   485
lemma fst_divmod_abs [simp]:
haftmann@51143
   486
  "fst (divmod_abs k l) = \<bar>k\<bar> div \<bar>l\<bar>"
haftmann@51143
   487
  by (simp add: divmod_abs_def)
haftmann@51143
   488
haftmann@51143
   489
lemma snd_divmod_abs [simp]:
haftmann@51143
   490
  "snd (divmod_abs k l) = \<bar>k\<bar> mod \<bar>l\<bar>"
haftmann@51143
   491
  by (simp add: divmod_abs_def)
haftmann@28708
   492
haftmann@53069
   493
lemma divmod_abs_code [code]:
haftmann@53069
   494
  "divmod_abs (Pos k) (Pos l) = divmod k l"
haftmann@53069
   495
  "divmod_abs (Neg k) (Neg l) = divmod k l"
haftmann@53069
   496
  "divmod_abs (Neg k) (Pos l) = divmod k l"
haftmann@53069
   497
  "divmod_abs (Pos k) (Neg l) = divmod k l"
haftmann@51143
   498
  "divmod_abs j 0 = (0, \<bar>j\<bar>)"
haftmann@51143
   499
  "divmod_abs 0 j = (0, 0)"
haftmann@51143
   500
  by (simp_all add: prod_eq_iff)
haftmann@51143
   501
haftmann@51143
   502
lemma divmod_integer_code [code]:
haftmann@51143
   503
  "divmod_integer k l =
haftmann@51143
   504
    (if k = 0 then (0, 0) else if l = 0 then (0, k) else
haftmann@51143
   505
    (apsnd \<circ> times \<circ> sgn) l (if sgn k = sgn l
haftmann@51143
   506
      then divmod_abs k l
haftmann@51143
   507
      else (let (r, s) = divmod_abs k l in
haftmann@51143
   508
        if s = 0 then (- r, 0) else (- r - 1, \<bar>l\<bar> - s))))"
haftmann@51143
   509
proof -
haftmann@51143
   510
  have aux1: "\<And>k l::int. sgn k = sgn l \<longleftrightarrow> k = 0 \<and> l = 0 \<or> 0 < l \<and> 0 < k \<or> l < 0 \<and> k < 0"
haftmann@51143
   511
    by (auto simp add: sgn_if)
haftmann@51143
   512
  have aux2: "\<And>q::int. - int_of_integer k = int_of_integer l * q \<longleftrightarrow> int_of_integer k = int_of_integer l * - q" by auto
haftmann@51143
   513
  show ?thesis
blanchet@55414
   514
    by (simp add: prod_eq_iff integer_eq_iff case_prod_beta aux1)
haftmann@51143
   515
      (auto simp add: zdiv_zminus1_eq_if zmod_zminus1_eq_if div_minus_right mod_minus_right aux2)
haftmann@51143
   516
qed
haftmann@51143
   517
haftmann@51143
   518
lemma div_integer_code [code]:
haftmann@51143
   519
  "k div l = fst (divmod_integer k l)"
haftmann@28708
   520
  by simp
haftmann@28708
   521
haftmann@51143
   522
lemma mod_integer_code [code]:
haftmann@51143
   523
  "k mod l = snd (divmod_integer k l)"
haftmann@25767
   524
  by simp
haftmann@24999
   525
haftmann@51143
   526
lemma equal_integer_code [code]:
haftmann@51143
   527
  "HOL.equal 0 (0::integer) \<longleftrightarrow> True"
haftmann@51143
   528
  "HOL.equal 0 (Pos l) \<longleftrightarrow> False"
haftmann@51143
   529
  "HOL.equal 0 (Neg l) \<longleftrightarrow> False"
haftmann@51143
   530
  "HOL.equal (Pos k) 0 \<longleftrightarrow> False"
haftmann@51143
   531
  "HOL.equal (Pos k) (Pos l) \<longleftrightarrow> HOL.equal k l"
haftmann@51143
   532
  "HOL.equal (Pos k) (Neg l) \<longleftrightarrow> False"
haftmann@51143
   533
  "HOL.equal (Neg k) 0 \<longleftrightarrow> False"
haftmann@51143
   534
  "HOL.equal (Neg k) (Pos l) \<longleftrightarrow> False"
haftmann@51143
   535
  "HOL.equal (Neg k) (Neg l) \<longleftrightarrow> HOL.equal k l"
haftmann@51143
   536
  by (simp_all add: equal)
haftmann@51143
   537
haftmann@51143
   538
lemma equal_integer_refl [code nbe]:
haftmann@51143
   539
  "HOL.equal (k::integer) k \<longleftrightarrow> True"
haftmann@51143
   540
  by (fact equal_refl)
haftmann@31266
   541
haftmann@51143
   542
lemma less_eq_integer_code [code]:
haftmann@51143
   543
  "0 \<le> (0::integer) \<longleftrightarrow> True"
haftmann@51143
   544
  "0 \<le> Pos l \<longleftrightarrow> True"
haftmann@51143
   545
  "0 \<le> Neg l \<longleftrightarrow> False"
haftmann@51143
   546
  "Pos k \<le> 0 \<longleftrightarrow> False"
haftmann@51143
   547
  "Pos k \<le> Pos l \<longleftrightarrow> k \<le> l"
haftmann@51143
   548
  "Pos k \<le> Neg l \<longleftrightarrow> False"
haftmann@51143
   549
  "Neg k \<le> 0 \<longleftrightarrow> True"
haftmann@51143
   550
  "Neg k \<le> Pos l \<longleftrightarrow> True"
haftmann@51143
   551
  "Neg k \<le> Neg l \<longleftrightarrow> l \<le> k"
haftmann@51143
   552
  by simp_all
haftmann@51143
   553
haftmann@51143
   554
lemma less_integer_code [code]:
haftmann@51143
   555
  "0 < (0::integer) \<longleftrightarrow> False"
haftmann@51143
   556
  "0 < Pos l \<longleftrightarrow> True"
haftmann@51143
   557
  "0 < Neg l \<longleftrightarrow> False"
haftmann@51143
   558
  "Pos k < 0 \<longleftrightarrow> False"
haftmann@51143
   559
  "Pos k < Pos l \<longleftrightarrow> k < l"
haftmann@51143
   560
  "Pos k < Neg l \<longleftrightarrow> False"
haftmann@51143
   561
  "Neg k < 0 \<longleftrightarrow> True"
haftmann@51143
   562
  "Neg k < Pos l \<longleftrightarrow> True"
haftmann@51143
   563
  "Neg k < Neg l \<longleftrightarrow> l < k"
haftmann@51143
   564
  by simp_all
haftmann@26140
   565
haftmann@51143
   566
lift_definition num_of_integer :: "integer \<Rightarrow> num"
haftmann@51143
   567
  is "num_of_nat \<circ> nat"
haftmann@51143
   568
  .
haftmann@51143
   569
haftmann@51143
   570
lemma num_of_integer_code [code]:
haftmann@51143
   571
  "num_of_integer k = (if k \<le> 1 then Num.One
haftmann@51143
   572
     else let
haftmann@51143
   573
       (l, j) = divmod_integer k 2;
haftmann@51143
   574
       l' = num_of_integer l;
haftmann@51143
   575
       l'' = l' + l'
haftmann@51143
   576
     in if j = 0 then l'' else l'' + Num.One)"
haftmann@51143
   577
proof -
haftmann@51143
   578
  {
haftmann@51143
   579
    assume "int_of_integer k mod 2 = 1"
haftmann@51143
   580
    then have "nat (int_of_integer k mod 2) = nat 1" by simp
haftmann@51143
   581
    moreover assume *: "1 < int_of_integer k"
haftmann@51143
   582
    ultimately have **: "nat (int_of_integer k) mod 2 = 1" by (simp add: nat_mod_distrib)
haftmann@51143
   583
    have "num_of_nat (nat (int_of_integer k)) =
haftmann@51143
   584
      num_of_nat (2 * (nat (int_of_integer k) div 2) + nat (int_of_integer k) mod 2)"
haftmann@51143
   585
      by simp
haftmann@51143
   586
    then have "num_of_nat (nat (int_of_integer k)) =
haftmann@51143
   587
      num_of_nat (nat (int_of_integer k) div 2 + nat (int_of_integer k) div 2 + nat (int_of_integer k) mod 2)"
haftmann@51143
   588
      by (simp add: mult_2)
haftmann@51143
   589
    with ** have "num_of_nat (nat (int_of_integer k)) =
haftmann@51143
   590
      num_of_nat (nat (int_of_integer k) div 2 + nat (int_of_integer k) div 2 + 1)"
haftmann@51143
   591
      by simp
haftmann@51143
   592
  }
haftmann@51143
   593
  note aux = this
haftmann@51143
   594
  show ?thesis
blanchet@55414
   595
    by (auto simp add: num_of_integer_def nat_of_integer_def Let_def case_prod_beta
haftmann@51143
   596
      not_le integer_eq_iff less_eq_integer_def
haftmann@51143
   597
      nat_mult_distrib nat_div_distrib num_of_nat_One num_of_nat_plus_distrib
haftmann@51143
   598
       mult_2 [where 'a=nat] aux add_One)
haftmann@25918
   599
qed
haftmann@25918
   600
haftmann@51143
   601
lemma nat_of_integer_code [code]:
haftmann@51143
   602
  "nat_of_integer k = (if k \<le> 0 then 0
haftmann@51143
   603
     else let
haftmann@51143
   604
       (l, j) = divmod_integer k 2;
haftmann@51143
   605
       l' = nat_of_integer l;
haftmann@51143
   606
       l'' = l' + l'
haftmann@51143
   607
     in if j = 0 then l'' else l'' + 1)"
haftmann@33340
   608
proof -
haftmann@51143
   609
  obtain j where "k = integer_of_int j"
haftmann@51143
   610
  proof
haftmann@51143
   611
    show "k = integer_of_int (int_of_integer k)" by simp
haftmann@51143
   612
  qed
haftmann@51143
   613
  moreover have "2 * (j div 2) = j - j mod 2"
haftmann@64246
   614
    by (simp add: minus_mod_eq_mult_div [symmetric] mult.commute)
haftmann@51143
   615
  ultimately show ?thesis
haftmann@63950
   616
    by (auto simp add: split_def Let_def modulo_integer_def nat_of_integer_def not_le
haftmann@51143
   617
      nat_add_distrib [symmetric] Suc_nat_eq_nat_zadd1)
haftmann@51143
   618
      (auto simp add: mult_2 [symmetric])
haftmann@33340
   619
qed
haftmann@28708
   620
haftmann@51143
   621
lemma int_of_integer_code [code]:
haftmann@51143
   622
  "int_of_integer k = (if k < 0 then - (int_of_integer (- k))
haftmann@51143
   623
     else if k = 0 then 0
haftmann@51143
   624
     else let
haftmann@51143
   625
       (l, j) = divmod_integer k 2;
haftmann@51143
   626
       l' = 2 * int_of_integer l
haftmann@51143
   627
     in if j = 0 then l' else l' + 1)"
haftmann@64246
   628
  by (auto simp add: split_def Let_def integer_eq_iff minus_mod_eq_mult_div [symmetric])
haftmann@28708
   629
haftmann@51143
   630
lemma integer_of_int_code [code]:
haftmann@51143
   631
  "integer_of_int k = (if k < 0 then - (integer_of_int (- k))
haftmann@51143
   632
     else if k = 0 then 0
haftmann@51143
   633
     else let
haftmann@60868
   634
       l = 2 * integer_of_int (k div 2);
haftmann@60868
   635
       j = k mod 2
haftmann@60868
   636
     in if j = 0 then l else l + 1)"
haftmann@64246
   637
  by (auto simp add: split_def Let_def integer_eq_iff minus_mod_eq_mult_div [symmetric])
haftmann@51143
   638
haftmann@51143
   639
hide_const (open) Pos Neg sub dup divmod_abs
huffman@46547
   640
haftmann@28708
   641
wenzelm@60758
   642
subsection \<open>Serializer setup for target language integers\<close>
haftmann@24999
   643
haftmann@51143
   644
code_reserved Eval int Integer abs
haftmann@25767
   645
haftmann@52435
   646
code_printing
haftmann@52435
   647
  type_constructor integer \<rightharpoonup>
haftmann@52435
   648
    (SML) "IntInf.int"
haftmann@52435
   649
    and (OCaml) "Big'_int.big'_int"
haftmann@52435
   650
    and (Haskell) "Integer"
haftmann@52435
   651
    and (Scala) "BigInt"
haftmann@52435
   652
    and (Eval) "int"
haftmann@52435
   653
| class_instance integer :: equal \<rightharpoonup>
haftmann@52435
   654
    (Haskell) -
haftmann@24999
   655
haftmann@52435
   656
code_printing
haftmann@52435
   657
  constant "0::integer" \<rightharpoonup>
haftmann@58400
   658
    (SML) "!(0/ :/ IntInf.int)"
haftmann@52435
   659
    and (OCaml) "Big'_int.zero'_big'_int"
haftmann@58400
   660
    and (Haskell) "!(0/ ::/ Integer)"
haftmann@52435
   661
    and (Scala) "BigInt(0)"
huffman@47108
   662
wenzelm@60758
   663
setup \<open>
haftmann@58399
   664
  fold (fn target =>
haftmann@58399
   665
    Numeral.add_code @{const_name Code_Numeral.Pos} I Code_Printer.literal_numeral target
haftmann@58399
   666
    #> Numeral.add_code @{const_name Code_Numeral.Neg} (op ~) Code_Printer.literal_numeral target)
haftmann@58399
   667
    ["SML", "OCaml", "Haskell", "Scala"]
wenzelm@60758
   668
\<close>
haftmann@51143
   669
haftmann@52435
   670
code_printing
haftmann@52435
   671
  constant "plus :: integer \<Rightarrow> _ \<Rightarrow> _" \<rightharpoonup>
haftmann@52435
   672
    (SML) "IntInf.+ ((_), (_))"
haftmann@52435
   673
    and (OCaml) "Big'_int.add'_big'_int"
haftmann@52435
   674
    and (Haskell) infixl 6 "+"
haftmann@52435
   675
    and (Scala) infixl 7 "+"
haftmann@52435
   676
    and (Eval) infixl 8 "+"
haftmann@52435
   677
| constant "uminus :: integer \<Rightarrow> _" \<rightharpoonup>
haftmann@52435
   678
    (SML) "IntInf.~"
haftmann@52435
   679
    and (OCaml) "Big'_int.minus'_big'_int"
haftmann@52435
   680
    and (Haskell) "negate"
haftmann@52435
   681
    and (Scala) "!(- _)"
haftmann@52435
   682
    and (Eval) "~/ _"
haftmann@52435
   683
| constant "minus :: integer \<Rightarrow> _" \<rightharpoonup>
haftmann@52435
   684
    (SML) "IntInf.- ((_), (_))"
haftmann@52435
   685
    and (OCaml) "Big'_int.sub'_big'_int"
haftmann@52435
   686
    and (Haskell) infixl 6 "-"
haftmann@52435
   687
    and (Scala) infixl 7 "-"
haftmann@52435
   688
    and (Eval) infixl 8 "-"
haftmann@52435
   689
| constant Code_Numeral.dup \<rightharpoonup>
haftmann@52435
   690
    (SML) "IntInf.*/ (2,/ (_))"
haftmann@52435
   691
    and (OCaml) "Big'_int.mult'_big'_int/ (Big'_int.big'_int'_of'_int/ 2)"
haftmann@52435
   692
    and (Haskell) "!(2 * _)"
haftmann@52435
   693
    and (Scala) "!(2 * _)"
haftmann@52435
   694
    and (Eval) "!(2 * _)"
haftmann@52435
   695
| constant Code_Numeral.sub \<rightharpoonup>
haftmann@52435
   696
    (SML) "!(raise/ Fail/ \"sub\")"
haftmann@52435
   697
    and (OCaml) "failwith/ \"sub\""
haftmann@52435
   698
    and (Haskell) "error/ \"sub\""
haftmann@52435
   699
    and (Scala) "!sys.error(\"sub\")"
haftmann@52435
   700
| constant "times :: integer \<Rightarrow> _ \<Rightarrow> _" \<rightharpoonup>
haftmann@52435
   701
    (SML) "IntInf.* ((_), (_))"
haftmann@52435
   702
    and (OCaml) "Big'_int.mult'_big'_int"
haftmann@52435
   703
    and (Haskell) infixl 7 "*"
haftmann@52435
   704
    and (Scala) infixl 8 "*"
haftmann@52435
   705
    and (Eval) infixl 9 "*"
haftmann@52435
   706
| constant Code_Numeral.divmod_abs \<rightharpoonup>
haftmann@52435
   707
    (SML) "IntInf.divMod/ (IntInf.abs _,/ IntInf.abs _)"
haftmann@52435
   708
    and (OCaml) "Big'_int.quomod'_big'_int/ (Big'_int.abs'_big'_int _)/ (Big'_int.abs'_big'_int _)"
haftmann@52435
   709
    and (Haskell) "divMod/ (abs _)/ (abs _)"
haftmann@52435
   710
    and (Scala) "!((k: BigInt) => (l: BigInt) =>/ if (l == 0)/ (BigInt(0), k) else/ (k.abs '/% l.abs))"
haftmann@52435
   711
    and (Eval) "Integer.div'_mod/ (abs _)/ (abs _)"
haftmann@52435
   712
| constant "HOL.equal :: integer \<Rightarrow> _ \<Rightarrow> bool" \<rightharpoonup>
haftmann@52435
   713
    (SML) "!((_ : IntInf.int) = _)"
haftmann@52435
   714
    and (OCaml) "Big'_int.eq'_big'_int"
haftmann@52435
   715
    and (Haskell) infix 4 "=="
haftmann@52435
   716
    and (Scala) infixl 5 "=="
haftmann@52435
   717
    and (Eval) infixl 6 "="
haftmann@52435
   718
| constant "less_eq :: integer \<Rightarrow> _ \<Rightarrow> bool" \<rightharpoonup>
haftmann@52435
   719
    (SML) "IntInf.<= ((_), (_))"
haftmann@52435
   720
    and (OCaml) "Big'_int.le'_big'_int"
haftmann@52435
   721
    and (Haskell) infix 4 "<="
haftmann@52435
   722
    and (Scala) infixl 4 "<="
haftmann@52435
   723
    and (Eval) infixl 6 "<="
haftmann@52435
   724
| constant "less :: integer \<Rightarrow> _ \<Rightarrow> bool" \<rightharpoonup>
haftmann@52435
   725
    (SML) "IntInf.< ((_), (_))"
haftmann@52435
   726
    and (OCaml) "Big'_int.lt'_big'_int"
haftmann@52435
   727
    and (Haskell) infix 4 "<"
haftmann@52435
   728
    and (Scala) infixl 4 "<"
haftmann@52435
   729
    and (Eval) infixl 6 "<"
Andreas@61857
   730
| constant "abs :: integer \<Rightarrow> _" \<rightharpoonup>
Andreas@61857
   731
    (SML) "IntInf.abs"
Andreas@61857
   732
    and (OCaml) "Big'_int.abs'_big'_int"
Andreas@61857
   733
    and (Haskell) "Prelude.abs"
Andreas@61857
   734
    and (Scala) "_.abs"
Andreas@61857
   735
    and (Eval) "abs"
haftmann@51143
   736
haftmann@52435
   737
code_identifier
haftmann@52435
   738
  code_module Code_Numeral \<rightharpoonup> (SML) Arith and (OCaml) Arith and (Haskell) Arith
huffman@46547
   739
haftmann@51143
   740
wenzelm@60758
   741
subsection \<open>Type of target language naturals\<close>
haftmann@51143
   742
wenzelm@61076
   743
typedef natural = "UNIV :: nat set"
haftmann@51143
   744
  morphisms nat_of_natural natural_of_nat ..
haftmann@51143
   745
haftmann@59487
   746
setup_lifting type_definition_natural
haftmann@51143
   747
haftmann@51143
   748
lemma natural_eq_iff [termination_simp]:
haftmann@51143
   749
  "m = n \<longleftrightarrow> nat_of_natural m = nat_of_natural n"
haftmann@51143
   750
  by transfer rule
haftmann@51143
   751
haftmann@51143
   752
lemma natural_eqI:
haftmann@51143
   753
  "nat_of_natural m = nat_of_natural n \<Longrightarrow> m = n"
haftmann@51143
   754
  using natural_eq_iff [of m n] by simp
haftmann@51143
   755
haftmann@51143
   756
lemma nat_of_natural_of_nat_inverse [simp]:
haftmann@51143
   757
  "nat_of_natural (natural_of_nat n) = n"
haftmann@51143
   758
  by transfer rule
haftmann@51143
   759
haftmann@51143
   760
lemma natural_of_nat_of_natural_inverse [simp]:
haftmann@51143
   761
  "natural_of_nat (nat_of_natural n) = n"
haftmann@51143
   762
  by transfer rule
haftmann@51143
   763
haftmann@51143
   764
instantiation natural :: "{comm_monoid_diff, semiring_1}"
haftmann@51143
   765
begin
haftmann@51143
   766
haftmann@51143
   767
lift_definition zero_natural :: natural
haftmann@51143
   768
  is "0 :: nat"
haftmann@51143
   769
  .
haftmann@51143
   770
haftmann@51143
   771
declare zero_natural.rep_eq [simp]
haftmann@51143
   772
haftmann@51143
   773
lift_definition one_natural :: natural
haftmann@51143
   774
  is "1 :: nat"
haftmann@51143
   775
  .
haftmann@51143
   776
haftmann@51143
   777
declare one_natural.rep_eq [simp]
haftmann@51143
   778
haftmann@51143
   779
lift_definition plus_natural :: "natural \<Rightarrow> natural \<Rightarrow> natural"
haftmann@51143
   780
  is "plus :: nat \<Rightarrow> nat \<Rightarrow> nat"
haftmann@51143
   781
  .
haftmann@51143
   782
haftmann@51143
   783
declare plus_natural.rep_eq [simp]
haftmann@51143
   784
haftmann@51143
   785
lift_definition minus_natural :: "natural \<Rightarrow> natural \<Rightarrow> natural"
haftmann@51143
   786
  is "minus :: nat \<Rightarrow> nat \<Rightarrow> nat"
haftmann@51143
   787
  .
haftmann@51143
   788
haftmann@51143
   789
declare minus_natural.rep_eq [simp]
haftmann@51143
   790
haftmann@51143
   791
lift_definition times_natural :: "natural \<Rightarrow> natural \<Rightarrow> natural"
haftmann@51143
   792
  is "times :: nat \<Rightarrow> nat \<Rightarrow> nat"
haftmann@51143
   793
  .
haftmann@51143
   794
haftmann@51143
   795
declare times_natural.rep_eq [simp]
haftmann@51143
   796
haftmann@51143
   797
instance proof
haftmann@51143
   798
qed (transfer, simp add: algebra_simps)+
haftmann@51143
   799
haftmann@51143
   800
end
haftmann@51143
   801
haftmann@64241
   802
instance natural :: Rings.dvd ..
haftmann@64241
   803
haftmann@64241
   804
lemma [transfer_rule]:
haftmann@64241
   805
  "rel_fun pcr_natural (rel_fun pcr_natural HOL.iff) Rings.dvd Rings.dvd"
haftmann@64241
   806
  unfolding dvd_def by transfer_prover
haftmann@64241
   807
haftmann@51143
   808
lemma [transfer_rule]:
blanchet@55945
   809
  "rel_fun HOL.eq pcr_natural (\<lambda>n::nat. n) (of_nat :: nat \<Rightarrow> natural)"
haftmann@51143
   810
proof -
blanchet@55945
   811
  have "rel_fun HOL.eq pcr_natural (of_nat :: nat \<Rightarrow> nat) (of_nat :: nat \<Rightarrow> natural)"
haftmann@51143
   812
    by (unfold of_nat_def [abs_def]) transfer_prover
haftmann@51143
   813
  then show ?thesis by (simp add: id_def)
haftmann@51143
   814
qed
haftmann@51143
   815
haftmann@51143
   816
lemma [transfer_rule]:
blanchet@55945
   817
  "rel_fun HOL.eq pcr_natural (numeral :: num \<Rightarrow> nat) (numeral :: num \<Rightarrow> natural)"
haftmann@51143
   818
proof -
blanchet@55945
   819
  have "rel_fun HOL.eq pcr_natural (numeral :: num \<Rightarrow> nat) (\<lambda>n. of_nat (numeral n))"
haftmann@51143
   820
    by transfer_prover
haftmann@51143
   821
  then show ?thesis by simp
haftmann@51143
   822
qed
haftmann@51143
   823
haftmann@51143
   824
lemma nat_of_natural_of_nat [simp]:
haftmann@51143
   825
  "nat_of_natural (of_nat n) = n"
haftmann@51143
   826
  by transfer rule
haftmann@51143
   827
haftmann@51143
   828
lemma natural_of_nat_of_nat [simp, code_abbrev]:
haftmann@51143
   829
  "natural_of_nat = of_nat"
haftmann@51143
   830
  by transfer rule
haftmann@51143
   831
haftmann@51143
   832
lemma of_nat_of_natural [simp]:
haftmann@51143
   833
  "of_nat (nat_of_natural n) = n"
haftmann@51143
   834
  by transfer rule
haftmann@51143
   835
haftmann@51143
   836
lemma nat_of_natural_numeral [simp]:
haftmann@51143
   837
  "nat_of_natural (numeral k) = numeral k"
haftmann@51143
   838
  by transfer rule
haftmann@51143
   839
haftmann@64592
   840
instantiation natural :: "{linordered_semiring, equal}"
haftmann@51143
   841
begin
haftmann@51143
   842
haftmann@51143
   843
lift_definition less_eq_natural :: "natural \<Rightarrow> natural \<Rightarrow> bool"
haftmann@51143
   844
  is "less_eq :: nat \<Rightarrow> nat \<Rightarrow> bool"
haftmann@51143
   845
  .
haftmann@51143
   846
haftmann@51143
   847
declare less_eq_natural.rep_eq [termination_simp]
haftmann@51143
   848
haftmann@51143
   849
lift_definition less_natural :: "natural \<Rightarrow> natural \<Rightarrow> bool"
haftmann@51143
   850
  is "less :: nat \<Rightarrow> nat \<Rightarrow> bool"
haftmann@51143
   851
  .
haftmann@51143
   852
haftmann@51143
   853
declare less_natural.rep_eq [termination_simp]
haftmann@51143
   854
haftmann@51143
   855
lift_definition equal_natural :: "natural \<Rightarrow> natural \<Rightarrow> bool"
haftmann@51143
   856
  is "HOL.equal :: nat \<Rightarrow> nat \<Rightarrow> bool"
haftmann@51143
   857
  .
haftmann@51143
   858
haftmann@51143
   859
instance proof
haftmann@51143
   860
qed (transfer, simp add: algebra_simps equal less_le_not_le [symmetric] linear)+
haftmann@51143
   861
haftmann@24999
   862
end
haftmann@46664
   863
haftmann@51143
   864
lemma [transfer_rule]:
blanchet@55945
   865
  "rel_fun pcr_natural (rel_fun pcr_natural pcr_natural) (min :: _ \<Rightarrow> _ \<Rightarrow> nat) (min :: _ \<Rightarrow> _ \<Rightarrow> natural)"
haftmann@51143
   866
  by (unfold min_def [abs_def]) transfer_prover
haftmann@51143
   867
haftmann@51143
   868
lemma [transfer_rule]:
blanchet@55945
   869
  "rel_fun pcr_natural (rel_fun pcr_natural pcr_natural) (max :: _ \<Rightarrow> _ \<Rightarrow> nat) (max :: _ \<Rightarrow> _ \<Rightarrow> natural)"
haftmann@51143
   870
  by (unfold max_def [abs_def]) transfer_prover
haftmann@51143
   871
haftmann@51143
   872
lemma nat_of_natural_min [simp]:
haftmann@51143
   873
  "nat_of_natural (min k l) = min (nat_of_natural k) (nat_of_natural l)"
haftmann@51143
   874
  by transfer rule
haftmann@51143
   875
haftmann@51143
   876
lemma nat_of_natural_max [simp]:
haftmann@51143
   877
  "nat_of_natural (max k l) = max (nat_of_natural k) (nat_of_natural l)"
haftmann@51143
   878
  by transfer rule
haftmann@51143
   879
haftmann@66806
   880
instantiation natural :: unique_euclidean_semiring
haftmann@64592
   881
begin
haftmann@64592
   882
haftmann@64592
   883
lift_definition normalize_natural :: "natural \<Rightarrow> natural"
haftmann@64592
   884
  is "normalize :: nat \<Rightarrow> nat"
haftmann@64592
   885
  .
haftmann@64592
   886
haftmann@64592
   887
declare normalize_natural.rep_eq [simp]
haftmann@64592
   888
haftmann@64592
   889
lift_definition unit_factor_natural :: "natural \<Rightarrow> natural"
haftmann@64592
   890
  is "unit_factor :: nat \<Rightarrow> nat"
haftmann@64592
   891
  .
haftmann@64592
   892
haftmann@64592
   893
declare unit_factor_natural.rep_eq [simp]
haftmann@64592
   894
haftmann@64592
   895
lift_definition divide_natural :: "natural \<Rightarrow> natural \<Rightarrow> natural"
haftmann@64592
   896
  is "divide :: nat \<Rightarrow> nat \<Rightarrow> nat"
haftmann@64592
   897
  .
haftmann@64592
   898
haftmann@64592
   899
declare divide_natural.rep_eq [simp]
haftmann@64592
   900
haftmann@64592
   901
lift_definition modulo_natural :: "natural \<Rightarrow> natural \<Rightarrow> natural"
haftmann@64592
   902
  is "modulo :: nat \<Rightarrow> nat \<Rightarrow> nat"
haftmann@64592
   903
  .
haftmann@64592
   904
haftmann@64592
   905
declare modulo_natural.rep_eq [simp]
haftmann@64592
   906
haftmann@66806
   907
lift_definition euclidean_size_natural :: "natural \<Rightarrow> nat"
haftmann@66806
   908
  is "euclidean_size :: nat \<Rightarrow> nat"
haftmann@66806
   909
  .
haftmann@66806
   910
haftmann@66806
   911
declare euclidean_size_natural.rep_eq [simp]
haftmann@66806
   912
haftmann@66806
   913
lift_definition uniqueness_constraint_natural :: "natural \<Rightarrow> natural \<Rightarrow> bool"
haftmann@66806
   914
  is "uniqueness_constraint :: nat \<Rightarrow> nat \<Rightarrow> bool"
haftmann@66806
   915
  .
haftmann@66806
   916
haftmann@66806
   917
declare uniqueness_constraint_natural.rep_eq [simp]
haftmann@66806
   918
haftmann@64592
   919
instance
haftmann@66806
   920
  by (standard; transfer)
haftmann@66806
   921
    (auto simp add: algebra_simps unit_factor_nat_def gr0_conv_Suc)
haftmann@64592
   922
haftmann@64592
   923
end
haftmann@64592
   924
haftmann@66806
   925
lemma [code]:
haftmann@66806
   926
  "euclidean_size = nat_of_natural"
haftmann@66806
   927
  by (simp add: fun_eq_iff)
haftmann@66806
   928
haftmann@66806
   929
lemma [code]:
haftmann@66806
   930
  "uniqueness_constraint = (\<top> :: natural \<Rightarrow> natural \<Rightarrow> bool)"
haftmann@66806
   931
  by (simp add: fun_eq_iff)
haftmann@66806
   932
haftmann@66815
   933
instance natural :: semiring_parity
haftmann@66815
   934
  by (standard; transfer) (simp_all add: of_nat_div odd_iff_mod_2_eq_one)
haftmann@66815
   935
haftmann@51143
   936
lift_definition natural_of_integer :: "integer \<Rightarrow> natural"
haftmann@51143
   937
  is "nat :: int \<Rightarrow> nat"
haftmann@51143
   938
  .
haftmann@51143
   939
haftmann@51143
   940
lift_definition integer_of_natural :: "natural \<Rightarrow> integer"
haftmann@51143
   941
  is "of_nat :: nat \<Rightarrow> int"
haftmann@51143
   942
  .
haftmann@51143
   943
haftmann@51143
   944
lemma natural_of_integer_of_natural [simp]:
haftmann@51143
   945
  "natural_of_integer (integer_of_natural n) = n"
haftmann@51143
   946
  by transfer simp
haftmann@51143
   947
haftmann@51143
   948
lemma integer_of_natural_of_integer [simp]:
haftmann@51143
   949
  "integer_of_natural (natural_of_integer k) = max 0 k"
haftmann@51143
   950
  by transfer auto
haftmann@51143
   951
haftmann@51143
   952
lemma int_of_integer_of_natural [simp]:
haftmann@51143
   953
  "int_of_integer (integer_of_natural n) = of_nat (nat_of_natural n)"
haftmann@51143
   954
  by transfer rule
haftmann@51143
   955
haftmann@51143
   956
lemma integer_of_natural_of_nat [simp]:
haftmann@51143
   957
  "integer_of_natural (of_nat n) = of_nat n"
haftmann@51143
   958
  by transfer rule
haftmann@51143
   959
haftmann@51143
   960
lemma [measure_function]:
haftmann@51143
   961
  "is_measure nat_of_natural"
haftmann@51143
   962
  by (rule is_measure_trivial)
haftmann@51143
   963
haftmann@51143
   964
wenzelm@60758
   965
subsection \<open>Inductive representation of target language naturals\<close>
haftmann@51143
   966
haftmann@51143
   967
lift_definition Suc :: "natural \<Rightarrow> natural"
haftmann@51143
   968
  is Nat.Suc
haftmann@51143
   969
  .
haftmann@51143
   970
haftmann@51143
   971
declare Suc.rep_eq [simp]
haftmann@51143
   972
blanchet@58306
   973
old_rep_datatype "0::natural" Suc
haftmann@51143
   974
  by (transfer, fact nat.induct nat.inject nat.distinct)+
haftmann@51143
   975
blanchet@55416
   976
lemma natural_cases [case_names nat, cases type: natural]:
haftmann@51143
   977
  fixes m :: natural
haftmann@51143
   978
  assumes "\<And>n. m = of_nat n \<Longrightarrow> P"
haftmann@51143
   979
  shows P
haftmann@51143
   980
  using assms by transfer blast
haftmann@51143
   981
blanchet@58390
   982
lemma [simp, code]: "size_natural = nat_of_natural"
blanchet@58390
   983
proof (rule ext)
blanchet@58390
   984
  fix n
blanchet@58390
   985
  show "size_natural n = nat_of_natural n"
blanchet@58390
   986
    by (induct n) simp_all
blanchet@58390
   987
qed
blanchet@58379
   988
blanchet@58390
   989
lemma [simp, code]: "size = nat_of_natural"
blanchet@58390
   990
proof (rule ext)
blanchet@58390
   991
  fix n
blanchet@58390
   992
  show "size n = nat_of_natural n"
blanchet@58390
   993
    by (induct n) simp_all
blanchet@58390
   994
qed
blanchet@58379
   995
haftmann@51143
   996
lemma natural_decr [termination_simp]:
haftmann@51143
   997
  "n \<noteq> 0 \<Longrightarrow> nat_of_natural n - Nat.Suc 0 < nat_of_natural n"
haftmann@51143
   998
  by transfer simp
haftmann@51143
   999
blanchet@58379
  1000
lemma natural_zero_minus_one: "(0::natural) - 1 = 0"
blanchet@58379
  1001
  by (rule zero_diff)
haftmann@51143
  1002
blanchet@58379
  1003
lemma Suc_natural_minus_one: "Suc n - 1 = n"
haftmann@51143
  1004
  by transfer simp
haftmann@51143
  1005
haftmann@51143
  1006
hide_const (open) Suc
haftmann@51143
  1007
haftmann@51143
  1008
wenzelm@60758
  1009
subsection \<open>Code refinement for target language naturals\<close>
haftmann@51143
  1010
haftmann@51143
  1011
lift_definition Nat :: "integer \<Rightarrow> natural"
haftmann@51143
  1012
  is nat
haftmann@51143
  1013
  .
haftmann@51143
  1014
haftmann@51143
  1015
lemma [code_post]:
haftmann@51143
  1016
  "Nat 0 = 0"
haftmann@51143
  1017
  "Nat 1 = 1"
haftmann@51143
  1018
  "Nat (numeral k) = numeral k"
haftmann@51143
  1019
  by (transfer, simp)+
haftmann@51143
  1020
haftmann@51143
  1021
lemma [code abstype]:
haftmann@51143
  1022
  "Nat (integer_of_natural n) = n"
haftmann@51143
  1023
  by transfer simp
haftmann@51143
  1024
haftmann@63174
  1025
lemma [code]:
haftmann@63174
  1026
  "natural_of_nat n = natural_of_integer (integer_of_nat n)"
haftmann@63174
  1027
  by transfer simp
haftmann@51143
  1028
haftmann@51143
  1029
lemma [code abstract]:
haftmann@51143
  1030
  "integer_of_natural (natural_of_integer k) = max 0 k"
haftmann@51143
  1031
  by simp
haftmann@51143
  1032
haftmann@51143
  1033
lemma [code_abbrev]:
haftmann@51143
  1034
  "natural_of_integer (Code_Numeral.Pos k) = numeral k"
haftmann@51143
  1035
  by transfer simp
haftmann@51143
  1036
haftmann@51143
  1037
lemma [code abstract]:
haftmann@51143
  1038
  "integer_of_natural 0 = 0"
haftmann@51143
  1039
  by transfer simp
haftmann@51143
  1040
haftmann@51143
  1041
lemma [code abstract]:
haftmann@51143
  1042
  "integer_of_natural 1 = 1"
haftmann@51143
  1043
  by transfer simp
haftmann@51143
  1044
haftmann@51143
  1045
lemma [code abstract]:
haftmann@51143
  1046
  "integer_of_natural (Code_Numeral.Suc n) = integer_of_natural n + 1"
haftmann@51143
  1047
  by transfer simp
haftmann@51143
  1048
haftmann@51143
  1049
lemma [code]:
haftmann@51143
  1050
  "nat_of_natural = nat_of_integer \<circ> integer_of_natural"
haftmann@51143
  1051
  by transfer (simp add: fun_eq_iff)
haftmann@51143
  1052
haftmann@51143
  1053
lemma [code, code_unfold]:
blanchet@55416
  1054
  "case_natural f g n = (if n = 0 then f else g (n - 1))"
haftmann@51143
  1055
  by (cases n rule: natural.exhaust) (simp_all, simp add: Suc_def)
haftmann@51143
  1056
blanchet@55642
  1057
declare natural.rec [code del]
haftmann@51143
  1058
haftmann@51143
  1059
lemma [code abstract]:
haftmann@51143
  1060
  "integer_of_natural (m + n) = integer_of_natural m + integer_of_natural n"
haftmann@51143
  1061
  by transfer simp
haftmann@51143
  1062
haftmann@51143
  1063
lemma [code abstract]:
haftmann@51143
  1064
  "integer_of_natural (m - n) = max 0 (integer_of_natural m - integer_of_natural n)"
haftmann@51143
  1065
  by transfer simp
haftmann@51143
  1066
haftmann@51143
  1067
lemma [code abstract]:
haftmann@51143
  1068
  "integer_of_natural (m * n) = integer_of_natural m * integer_of_natural n"
haftmann@64592
  1069
  by transfer simp
haftmann@64592
  1070
haftmann@64592
  1071
lemma [code]:
haftmann@64592
  1072
  "normalize n = n" for n :: natural
haftmann@64592
  1073
  by transfer simp
haftmann@64592
  1074
haftmann@64592
  1075
lemma [code]:
haftmann@64592
  1076
  "unit_factor n = of_bool (n \<noteq> 0)" for n :: natural
haftmann@64592
  1077
proof (cases "n = 0")
haftmann@64592
  1078
  case True
haftmann@64592
  1079
  then show ?thesis
haftmann@64592
  1080
    by simp
haftmann@64592
  1081
next
haftmann@64592
  1082
  case False
haftmann@64592
  1083
  then have "unit_factor n = 1"
haftmann@64592
  1084
  proof transfer
haftmann@64592
  1085
    fix n :: nat
haftmann@64592
  1086
    assume "n \<noteq> 0"
haftmann@64592
  1087
    then obtain m where "n = Suc m"
haftmann@64592
  1088
      by (cases n) auto
haftmann@64592
  1089
    then show "unit_factor n = 1"
haftmann@64592
  1090
      by simp
haftmann@64592
  1091
  qed
haftmann@64592
  1092
  with False show ?thesis
haftmann@64592
  1093
    by simp
haftmann@64592
  1094
qed
haftmann@51143
  1095
haftmann@51143
  1096
lemma [code abstract]:
haftmann@51143
  1097
  "integer_of_natural (m div n) = integer_of_natural m div integer_of_natural n"
haftmann@51143
  1098
  by transfer (simp add: zdiv_int)
haftmann@51143
  1099
haftmann@51143
  1100
lemma [code abstract]:
haftmann@51143
  1101
  "integer_of_natural (m mod n) = integer_of_natural m mod integer_of_natural n"
haftmann@51143
  1102
  by transfer (simp add: zmod_int)
haftmann@51143
  1103
haftmann@51143
  1104
lemma [code]:
haftmann@51143
  1105
  "HOL.equal m n \<longleftrightarrow> HOL.equal (integer_of_natural m) (integer_of_natural n)"
haftmann@51143
  1106
  by transfer (simp add: equal)
haftmann@51143
  1107
blanchet@58379
  1108
lemma [code nbe]: "HOL.equal n (n::natural) \<longleftrightarrow> True"
blanchet@58379
  1109
  by (rule equal_class.equal_refl)
haftmann@51143
  1110
blanchet@58379
  1111
lemma [code]: "m \<le> n \<longleftrightarrow> integer_of_natural m \<le> integer_of_natural n"
haftmann@51143
  1112
  by transfer simp
haftmann@51143
  1113
blanchet@58379
  1114
lemma [code]: "m < n \<longleftrightarrow> integer_of_natural m < integer_of_natural n"
haftmann@51143
  1115
  by transfer simp
haftmann@51143
  1116
haftmann@51143
  1117
hide_const (open) Nat
haftmann@51143
  1118
kuncar@55736
  1119
lifting_update integer.lifting
kuncar@55736
  1120
lifting_forget integer.lifting
kuncar@55736
  1121
kuncar@55736
  1122
lifting_update natural.lifting
kuncar@55736
  1123
lifting_forget natural.lifting
haftmann@51143
  1124
haftmann@51143
  1125
code_reflect Code_Numeral
haftmann@63174
  1126
  datatypes natural
haftmann@63174
  1127
  functions "Code_Numeral.Suc" "0 :: natural" "1 :: natural"
haftmann@63174
  1128
    "plus :: natural \<Rightarrow> _" "minus :: natural \<Rightarrow> _"
haftmann@63174
  1129
    "times :: natural \<Rightarrow> _" "divide :: natural \<Rightarrow> _"
haftmann@63950
  1130
    "modulo :: natural \<Rightarrow> _"
haftmann@63174
  1131
    integer_of_natural natural_of_integer
haftmann@51143
  1132
haftmann@51143
  1133
end