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