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