Standard part function

Given a finite element in a non-archimedean field, the standard part function rounds it to the unique closest real number. That is, it chops off any infinitesimals.

Let K be a linearly ordered field. The subset of finite elements (i.e. those bounded by a natural number) is a ValuationSubring, which means we can construct its residue field FiniteResidueField, roughly corresponding to the finite elements quotiented by infinitesimals. This field inherits a LinearOrder instance, which makes it into an Archimedean linearly ordered field, meaning we can uniquely embed it in the reals.

Given a finite element of the field, the ArchimedeanClass.stdPart function returns the real number corresponding to this unique embedding. This function generalizes, among other things, the standard part function on Hyperreal.

TODO

Redefine Hyperreal.st in terms of ArchimedeanClass.stdPart.

References

• https://en.wikipedia.org/wiki/Standard_part_function

Finite residue field

noncomputable def ArchimedeanClass.FiniteElement (K : Type u_1) [LinearOrder K] [Field K] [IsOrderedRing K] : Type u_1

The valuation subring of elements in non-negative Archimedean classes, i.e. elements bounded by some natural number.

Equations

• ArchimedeanClass.FiniteElement K = ↥(AddValuation.toValuation (ArchimedeanClass.addValuation K)).valuationSubring

instance ArchimedeanClass.FiniteElement.instCommRing {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : CommRing (FiniteElement K)

Equations

• ArchimedeanClass.FiniteElement.instCommRing = id inferInstance

instance ArchimedeanClass.FiniteElement.instIsDomain {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : IsDomain (FiniteElement K)

instance ArchimedeanClass.FiniteElement.instValuationRing {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : ValuationRing (FiniteElement K)

instance ArchimedeanClass.FiniteElement.instLinearOrder {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : LinearOrder (FiniteElement K)

Equations

• ArchimedeanClass.FiniteElement.instLinearOrder = id inferInstance

instance ArchimedeanClass.FiniteElement.instIsStrictOrderedRing {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : IsStrictOrderedRing (FiniteElement K)

@[simp]

theorem ArchimedeanClass.FiniteElement.val_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : ↑0 = 0

@[simp]

theorem ArchimedeanClass.FiniteElement.val_one {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : ↑1 = 1

@[simp]

theorem ArchimedeanClass.FiniteElement.val_add {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x y : FiniteElement K) : ↑(x + y) = ↑x + ↑y

@[simp]

theorem ArchimedeanClass.FiniteElement.val_sub {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x y : FiniteElement K) : ↑(x - y) = ↑x - ↑y

@[simp]

theorem ArchimedeanClass.FiniteElement.val_mul {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x y : FiniteElement K) : ↑(x * y) = ↑x * ↑y

theorem ArchimedeanClass.FiniteElement.ext {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : FiniteElement K} (h : ↑x = ↑y) : x = y

theorem ArchimedeanClass.FiniteElement.ext_iff {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : FiniteElement K} : x = y ↔ ↑x = ↑y

def ArchimedeanClass.FiniteElement.mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x : K) (h : 0 ≤ mk x) : FiniteElement K

The constructor for FiniteElement.

Equations

• ArchimedeanClass.FiniteElement.mk x h = ⟨x, h⟩

@[simp]

theorem ArchimedeanClass.FiniteElement.mk_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (h : 0 ≤ mk 0) : FiniteElement.mk 0 h = 0

@[simp]

theorem ArchimedeanClass.FiniteElement.mk_one {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (h : 0 ≤ mk 1) : FiniteElement.mk 1 h = 1

@[simp]

theorem ArchimedeanClass.FiniteElement.mk_natCast {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {n : ℕ} (h : 0 ≤ mk ↑n) : FiniteElement.mk (↑n) h = ↑n

@[simp]

theorem ArchimedeanClass.FiniteElement.mk_intCast {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {n : ℤ} (h : 0 ≤ mk ↑n) : FiniteElement.mk (↑n) h = ↑n

@[simp]

theorem ArchimedeanClass.FiniteElement.mk_neg {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : K} (h : 0 ≤ mk x) : -FiniteElement.mk x h = FiniteElement.mk (-x) ⋯

theorem ArchimedeanClass.FiniteElement.not_isUnit_iff_mk_pos {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : FiniteElement K} : ¬IsUnit x ↔ 0 < mk ↑x

theorem ArchimedeanClass.FiniteElement.isUnit_iff_mk_eq_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : FiniteElement K} : IsUnit x ↔ mk ↑x = 0

def ArchimedeanClass.FiniteResidueField (K : Type u_1) [LinearOrder K] [Field K] [IsOrderedRing K] : Type u_1

The residue field of FiniteElement. This quotient inherits an order from K, which makes it into a linearly ordered Archimedean field.

Equations

• ArchimedeanClass.FiniteResidueField K = IsLocalRing.ResidueField (ArchimedeanClass.FiniteElement K)

noncomputable instance ArchimedeanClass.FiniteResidueField.instField {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : Field (FiniteResidueField K)

Equations

• ArchimedeanClass.FiniteResidueField.instField = inferInstanceAs (Field (IsLocalRing.ResidueField (ArchimedeanClass.FiniteElement K)))

theorem ArchimedeanClass.FiniteResidueField.ordConnected_preimage_mk' {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x : Quotient (Submodule.quotientRel (IsLocalRing.maximalIdeal (FiniteElement K)))) : (Quotient.mk (Submodule.quotientRel (IsLocalRing.maximalIdeal (FiniteElement K))) ⁻¹' {x}).OrdConnected

noncomputable instance ArchimedeanClass.FiniteResidueField.instLinearOrder {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : LinearOrder (FiniteResidueField K)

Equations

• ArchimedeanClass.FiniteResidueField.instLinearOrder = Quotient.instLinearOrder

def ArchimedeanClass.FiniteResidueField.mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : FiniteElement K →+*o FiniteResidueField K

The quotient map from finite elements on the field to the associated residue field.

Equations

• ArchimedeanClass.FiniteResidueField.mk = { toRingHom := IsLocalRing.residue (ArchimedeanClass.FiniteElement K), monotone' := ⋯ }

theorem ArchimedeanClass.FiniteResidueField.ind {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {motive : FiniteResidueField K → Prop} (mk : ∀ (x : FiniteElement K), motive (mk x)) (x : FiniteResidueField K) : motive x

instance ArchimedeanClass.FiniteResidueField.ordConnected_preimage_mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x : FiniteResidueField K) : (⇑mk ⁻¹' {x}).OrdConnected

theorem ArchimedeanClass.FiniteResidueField.mk_eq_mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : FiniteElement K} : mk x = mk y ↔ 0 < ArchimedeanClass.mk (↑x - ↑y)

theorem ArchimedeanClass.FiniteResidueField.mk_eq_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : FiniteElement K} : mk x = 0 ↔ 0 < ArchimedeanClass.mk ↑x

theorem ArchimedeanClass.FiniteResidueField.mk_ne_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : FiniteElement K} : mk x ≠ 0 ↔ ArchimedeanClass.mk ↑x = 0

theorem ArchimedeanClass.FiniteResidueField.mk_le_mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : FiniteElement K} : mk x ≤ mk y ↔ x ≤ y ∨ mk x = mk y

theorem ArchimedeanClass.FiniteResidueField.mk_lt_mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : FiniteElement K} : mk x < mk y ↔ x < y ∧ mk x ≠ mk y

theorem ArchimedeanClass.FiniteResidueField.lt_of_mk_lt_mk {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : FiniteElement K} (h : mk x < mk y) : x < y

instance ArchimedeanClass.FiniteResidueField.instIsOrderedRing {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : IsOrderedRing (FiniteResidueField K)

instance ArchimedeanClass.FiniteResidueField.instArchimedean {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : Archimedean (FiniteResidueField K)

def ArchimedeanClass.FiniteResidueField.ofArchimedean {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {R : Type u_2} [LinearOrder R] [CommRing R] [IsOrderedRing R] [Archimedean R] (f : R →+*o K) : R →+*o FiniteResidueField K

An embedding from an Archimedean field into K induces an embedding into FiniteResidueField K.

Equations

• One or more equations did not get rendered due to their size.

theorem ArchimedeanClass.FiniteResidueField.ofArchimedean_apply {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {R : Type u_2} [LinearOrder R] [CommRing R] [IsOrderedRing R] [Archimedean R] (f : R →+*o K) (r : R) : (ofArchimedean f) r = mk (FiniteElement.mk (f r) ⋯)

theorem ArchimedeanClass.FiniteResidueField.ofArchimedean_injective {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {R : Type u_2} [LinearOrder R] [CommRing R] [IsOrderedRing R] [Archimedean R] (f : R →+*o K) : Function.Injective ⇑(ofArchimedean f)

@[simp]

theorem ArchimedeanClass.FiniteResidueField.ofArchimedean_inj {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {R : Type u_2} [LinearOrder R] [CommRing R] [IsOrderedRing R] [Archimedean R] (f : R →+*o K) {x y : R} : (ofArchimedean f) x = (ofArchimedean f) y ↔ x = y

Standard part

noncomputable def ArchimedeanClass.stdPart {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x : K) : ℝ

The standard part of a FiniteElement is the unique real number with an infinitesimal difference.

For any infinite inputs, this function outputs a junk value of 0.

Equations

• ArchimedeanClass.stdPart x = if h : 0 ≤ ArchimedeanClass.mk x then (default.comp ArchimedeanClass.FiniteResidueField.mk) (ArchimedeanClass.FiniteElement.mk x h) else 0

theorem ArchimedeanClass.stdPart_of_mk_ne_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : K} (h : mk x ≠ 0) : stdPart x = 0

theorem ArchimedeanClass.stdPart_of_mk_nonneg {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : K} (f : FiniteResidueField K →+*o ℝ) (h : 0 ≤ mk x) : stdPart x = f (FiniteResidueField.mk (FiniteElement.mk x h))

@[simp]

theorem ArchimedeanClass.stdPart_zero {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : stdPart 0 = 0

@[simp]

theorem ArchimedeanClass.stdPart_one {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] : stdPart 1 = 1

@[simp]

theorem ArchimedeanClass.stdPart_neg {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x : K) : stdPart (-x) = -stdPart x

@[simp]

theorem ArchimedeanClass.stdPart_inv {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (x : K) : stdPart x⁻¹ = (stdPart x)⁻¹

theorem ArchimedeanClass.stdPart_add {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : K} (hx : 0 ≤ mk x) (hy : 0 ≤ mk y) : stdPart (x + y) = stdPart x + stdPart y

theorem ArchimedeanClass.stdPart_sub {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : K} (hx : 0 ≤ mk x) (hy : 0 ≤ mk y) : stdPart (x - y) = stdPart x - stdPart y

theorem ArchimedeanClass.stdPart_mul {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : K} (hx : 0 ≤ mk x) (hy : 0 ≤ mk y) : stdPart (x * y) = stdPart x * stdPart y

theorem ArchimedeanClass.stdPart_div {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x y : K} (hx : 0 ≤ mk x) (hy : 0 ≤ -mk y) : stdPart (x / y) = stdPart x / stdPart y

@[simp]

theorem ArchimedeanClass.stdPart_intCast {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (n : ℤ) : stdPart ↑n = ↑n

@[simp]

theorem ArchimedeanClass.stdPart_natCast {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (n : ℕ) : stdPart ↑n = ↑n

@[simp]

theorem ArchimedeanClass.stdPart_ofNat {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (n : ℕ) [n.AtLeastTwo] : stdPart (OfNat.ofNat n) = ↑n

@[simp]

theorem ArchimedeanClass.stdPart_ratCast {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (q : ℚ) : stdPart ↑q = ↑q

@[simp]

theorem ArchimedeanClass.stdPart_real {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] (f : ℝ →+*o K) (r : ℝ) : stdPart (f r) = r

theorem ArchimedeanClass.ofArchimedean_stdPart {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : K} (f : ℝ →+*o K) (hx : 0 ≤ mk x) : (FiniteResidueField.ofArchimedean f) (stdPart x) = FiniteResidueField.mk (FiniteElement.mk x hx)

theorem ArchimedeanClass.mk_sub_pos_iff {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : K} (f : ℝ →+*o K) {r : ℝ} (hx : 0 ≤ mk x) : 0 < mk (x - f r) ↔ stdPart x = r

The standard part of x is the unique real r such that x - r is infinitesimal.

theorem ArchimedeanClass.mk_sub_stdPart_pos {K : Type u_1} [LinearOrder K] [Field K] [IsOrderedRing K] {x : K} (f : ℝ →+*o K) (hx : 0 ≤ mk x) : 0 < mk (x - f (stdPart x))