The upper half plane and its automorphisms

This file defines UpperHalfPlane to be the upper half plane in ℂ.

We furthermore equip it with the structure of a GLPos 2 ℝ action by fractional linear transformations.

We define the notation ℍ for the upper half plane available in the locale UpperHalfPlane so as not to conflict with the quaternions.

def UpperHalfPlane : Type

The open upper half plane

Equations

• UpperHalfPlane = { point : ℂ // 0 < point.im }

def UpperHalfPlane.termℍ : Lean.ParserDescr

The open upper half plane

Equations

• UpperHalfPlane.termℍ = Lean.ParserDescr.node `UpperHalfPlane.termℍ 1024 (Lean.ParserDescr.symbol "ℍ")

def UpperHalfPlane.«term↑ₘ_» : Lean.ParserDescr

The coercion first into an element of GL(2, ℝ)⁺, then GL(2, ℝ) and finally a 2 × 2 matrix.

This notation is scoped in namespace UpperHalfPlane.

Equations

• UpperHalfPlane.«term↑ₘ_» = Lean.ParserDescr.node `UpperHalfPlane.«term↑ₘ_» 1024 (Lean.ParserDescr.binary `andthen (Lean.ParserDescr.symbol "↑ₘ") (Lean.ParserDescr.cat `term 1024))

instance UpperHalfPlane.instCoeFun : CoeFun ↥(Matrix.GLPos (Fin 2) ℝ) fun (x : ↥(Matrix.GLPos (Fin 2) ℝ)) => Fin 2 → Fin 2 → ℝ

Equations

• UpperHalfPlane.instCoeFun = { coe := fun (A : ↥(Matrix.GLPos (Fin 2) ℝ)) => ↑↑A }

def UpperHalfPlane.«term↑ₘ[_]_» : Lean.ParserDescr

The coercion into an element of GL(2, R) and finally a 2 × 2 matrix over R. This is similar to ↑ₘ, but without positivity requirements, and allows the user to specify the ring R, which can be useful to help Lean elaborate correctly.

This notation is scoped in namespace UpperHalfPlane.

Equations

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

def UpperHalfPlane.coe (z : UpperHalfPlane) : ℂ

Canonical embedding of the upper half-plane into ℂ.

Equations

• ↑z = ↑z

instance UpperHalfPlane.instCoeOutComplex : CoeOut UpperHalfPlane ℂ

Equations

• UpperHalfPlane.instCoeOutComplex = { coe := UpperHalfPlane.coe }

instance UpperHalfPlane.instInhabited : Inhabited UpperHalfPlane

Equations

• UpperHalfPlane.instInhabited = { default := ⟨Complex.I, UpperHalfPlane.instInhabited.proof_1⟩ }

theorem UpperHalfPlane.ext {a : UpperHalfPlane} {b : UpperHalfPlane} (h : ↑a = ↑b) : a = b

@[simp]

theorem UpperHalfPlane.ext_iff' {a : UpperHalfPlane} {b : UpperHalfPlane} : ↑a = ↑b ↔ a = b

instance UpperHalfPlane.canLift : CanLift ℂ UpperHalfPlane UpperHalfPlane.coe fun (z : ℂ) => 0 < z.im

Equations

• UpperHalfPlane.canLift = ⋯

def UpperHalfPlane.im (z : UpperHalfPlane) : ℝ

Imaginary part

Equations

• z.im = (↑z).im

def UpperHalfPlane.re (z : UpperHalfPlane) : ℝ

Real part

Equations

• z.re = (↑z).re

theorem UpperHalfPlane.ext' {a : UpperHalfPlane} {b : UpperHalfPlane} (hre : a.re = b.re) (him : a.im = b.im) : a = b

Extensionality lemma in terms of UpperHalfPlane.re and UpperHalfPlane.im.

def UpperHalfPlane.mk (z : ℂ) (h : 0 < z.im) : UpperHalfPlane

Constructor for UpperHalfPlane. It is useful if ⟨z, h⟩ makes Lean use a wrong typeclass instance.

Equations

• UpperHalfPlane.mk z h = ⟨z, h⟩

@[simp]

theorem UpperHalfPlane.coe_im (z : UpperHalfPlane) : (↑z).im = z.im

@[simp]

theorem UpperHalfPlane.coe_re (z : UpperHalfPlane) : (↑z).re = z.re

@[simp]

theorem UpperHalfPlane.mk_re (z : ℂ) (h : 0 < z.im) : (UpperHalfPlane.mk z h).re = z.re

@[simp]

theorem UpperHalfPlane.mk_im (z : ℂ) (h : 0 < z.im) : (UpperHalfPlane.mk z h).im = z.im

@[simp]

theorem UpperHalfPlane.coe_mk (z : ℂ) (h : 0 < z.im) : ↑(UpperHalfPlane.mk z h) = z

@[simp]

theorem UpperHalfPlane.mk_coe (z : UpperHalfPlane) (h : optParam (0 < (↑z).im) ⋯) : UpperHalfPlane.mk (↑z) h = z

theorem UpperHalfPlane.re_add_im (z : UpperHalfPlane) : ↑z.re + ↑z.im * Complex.I = ↑z

theorem UpperHalfPlane.im_pos (z : UpperHalfPlane) : 0 < z.im

theorem UpperHalfPlane.im_ne_zero (z : UpperHalfPlane) : z.im ≠ 0

theorem UpperHalfPlane.ne_zero (z : UpperHalfPlane) : ↑z ≠ 0

def UpperHalfPlane.I : UpperHalfPlane

Define I := √-1 as an element on the upper half plane.

Equations

• UpperHalfPlane.I = ⟨Complex.I, UpperHalfPlane.I.proof_1⟩

@[simp]

theorem UpperHalfPlane.I_im : UpperHalfPlane.I.im = 1

@[simp]

theorem UpperHalfPlane.I_re : UpperHalfPlane.I.re = 0

@[simp]

theorem UpperHalfPlane.coe_I : ↑UpperHalfPlane.I = Complex.I

def Mathlib.Meta.Positivity.evalUpperHalfPlaneIm : Mathlib.Meta.Positivity.PositivityExt

Extension for the positivity tactic: UpperHalfPlane.im.

def Mathlib.Meta.Positivity.evalUpperHalfPlaneCoe : Mathlib.Meta.Positivity.PositivityExt

Extension for the positivity tactic: UpperHalfPlane.coe.

theorem UpperHalfPlane.normSq_pos (z : UpperHalfPlane) : 0 < Complex.normSq ↑z

theorem UpperHalfPlane.normSq_ne_zero (z : UpperHalfPlane) : Complex.normSq ↑z ≠ 0

theorem UpperHalfPlane.im_inv_neg_coe_pos (z : UpperHalfPlane) : 0 < (-↑z)⁻¹.im

def UpperHalfPlane.num (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : ℂ

Numerator of the formula for a fractional linear transformation

Equations

• UpperHalfPlane.num g z = ↑(↑↑g 0 0) * ↑z + ↑(↑↑g 0 1)

def UpperHalfPlane.denom (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : ℂ

Denominator of the formula for a fractional linear transformation

Equations

• UpperHalfPlane.denom g z = ↑(↑↑g 1 0) * ↑z + ↑(↑↑g 1 1)

theorem UpperHalfPlane.linear_ne_zero (cd : Fin 2 → ℝ) (z : UpperHalfPlane) (h : cd ≠ 0) : ↑(cd 0) * ↑z + ↑(cd 1) ≠ 0

theorem UpperHalfPlane.denom_ne_zero (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : UpperHalfPlane.denom g z ≠ 0

theorem UpperHalfPlane.normSq_denom_pos (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : 0 < Complex.normSq (UpperHalfPlane.denom g z)

theorem UpperHalfPlane.normSq_denom_ne_zero (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : Complex.normSq (UpperHalfPlane.denom g z) ≠ 0

def UpperHalfPlane.smulAux' (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : ℂ

Fractional linear transformation, also known as the Moebius transformation

Equations

• UpperHalfPlane.smulAux' g z = UpperHalfPlane.num g z / UpperHalfPlane.denom g z

theorem UpperHalfPlane.smulAux'_im (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : (UpperHalfPlane.smulAux' g z).im = (↑↑g).det * z.im / Complex.normSq (UpperHalfPlane.denom g z)

def UpperHalfPlane.smulAux (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : UpperHalfPlane

Fractional linear transformation, also known as the Moebius transformation

Equations

• UpperHalfPlane.smulAux g z = UpperHalfPlane.mk (UpperHalfPlane.smulAux' g z) ⋯

theorem UpperHalfPlane.denom_cocycle (x : ↥(Matrix.GLPos (Fin 2) ℝ)) (y : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : UpperHalfPlane.denom (x * y) z = UpperHalfPlane.denom x (UpperHalfPlane.smulAux y z) * UpperHalfPlane.denom y z

theorem UpperHalfPlane.mul_smul' (x : ↥(Matrix.GLPos (Fin 2) ℝ)) (y : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : UpperHalfPlane.smulAux (x * y) z = UpperHalfPlane.smulAux x (UpperHalfPlane.smulAux y z)

instance UpperHalfPlane.instMulActionSubtypeGeneralLinearGroupFinOfNatNatRealMemSubgroupGLPos : MulAction (↥(Matrix.GLPos (Fin 2) ℝ)) UpperHalfPlane

The action of GLPos 2 ℝ on the upper half-plane by fractional linear transformations.

Equations

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

instance UpperHalfPlane.SLAction {R : Type u_1} [CommRing R] [Algebra R ℝ] : MulAction (Matrix.SpecialLinearGroup (Fin 2) R) UpperHalfPlane

Equations

• UpperHalfPlane.SLAction = MulAction.compHom UpperHalfPlane (Matrix.SpecialLinearGroup.toGLPos.comp (Matrix.SpecialLinearGroup.map (algebraMap R ℝ)))

def UpperHalfPlane.ModularGroup.coe' : Matrix.SpecialLinearGroup (Fin 2) ℤ → ↥(Matrix.GLPos (Fin 2) ℝ)

Canonical embedding of SL(2, ℤ) into GL(2, ℝ)⁺.

Equations

• ↑g = Matrix.SpecialLinearGroup.toGLPos ((Matrix.SpecialLinearGroup.map (Int.castRingHom ℝ)) g)

instance UpperHalfPlane.ModularGroup.instCoeSpecialLinearGroupFinOfNatNatIntSubtypeGeneralLinearGroupRealMemSubgroupGLPos : Coe (Matrix.SpecialLinearGroup (Fin 2) ℤ) ↥(Matrix.GLPos (Fin 2) ℝ)

Equations

• UpperHalfPlane.ModularGroup.instCoeSpecialLinearGroupFinOfNatNatIntSubtypeGeneralLinearGroupRealMemSubgroupGLPos = { coe := UpperHalfPlane.ModularGroup.coe' }

@[simp]

theorem UpperHalfPlane.ModularGroup.coe'_apply_complex {g : Matrix.SpecialLinearGroup (Fin 2) ℤ} {i : Fin 2} {j : Fin 2} : ↑(↑↑↑g i j) = ↑(↑g i j)

@[simp]

theorem UpperHalfPlane.ModularGroup.det_coe' {g : Matrix.SpecialLinearGroup (Fin 2) ℤ} : (↑↑↑g).det = 1

theorem UpperHalfPlane.ModularGroup.coe_one : ↑1 = 1

instance UpperHalfPlane.ModularGroup.SLOnGLPos : SMul (Matrix.SpecialLinearGroup (Fin 2) ℤ) ↥(Matrix.GLPos (Fin 2) ℝ)

Equations

• UpperHalfPlane.ModularGroup.SLOnGLPos = { smul := fun (s : Matrix.SpecialLinearGroup (Fin 2) ℤ) (g : ↥(Matrix.GLPos (Fin 2) ℝ)) => ↑s * g }

theorem UpperHalfPlane.ModularGroup.SLOnGLPos_smul_apply (s : Matrix.SpecialLinearGroup (Fin 2) ℤ) (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : (s • g) • z = (↑s * g) • z

instance UpperHalfPlane.ModularGroup.SL_to_GL_tower : IsScalarTower (Matrix.SpecialLinearGroup (Fin 2) ℤ) (↥(Matrix.GLPos (Fin 2) ℝ)) UpperHalfPlane

Equations

• UpperHalfPlane.ModularGroup.SL_to_GL_tower = ⋯

theorem UpperHalfPlane.specialLinearGroup_apply {R : Type u_1} [CommRing R] [Algebra R ℝ] (g : Matrix.SpecialLinearGroup (Fin 2) R) (z : UpperHalfPlane) : g • z = UpperHalfPlane.mk ((↑((algebraMap R ℝ) (↑g 0 0)) * ↑z + ↑((algebraMap R ℝ) (↑g 0 1))) / (↑((algebraMap R ℝ) (↑g 1 0)) * ↑z + ↑((algebraMap R ℝ) (↑g 1 1)))) ⋯

@[simp]

theorem UpperHalfPlane.coe_smul (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : ↑(g • z) = UpperHalfPlane.num g z / UpperHalfPlane.denom g z

@[simp]

theorem UpperHalfPlane.re_smul (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : (g • z).re = (UpperHalfPlane.num g z / UpperHalfPlane.denom g z).re

theorem UpperHalfPlane.im_smul (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : (g • z).im = (UpperHalfPlane.num g z / UpperHalfPlane.denom g z).im

theorem UpperHalfPlane.im_smul_eq_div_normSq (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : (g • z).im = (↑↑g).det * z.im / Complex.normSq (UpperHalfPlane.denom g z)

theorem UpperHalfPlane.c_mul_im_sq_le_normSq_denom (z : UpperHalfPlane) (g : Matrix.SpecialLinearGroup (Fin 2) ℝ) : (↑g 1 0 * z.im) ^ 2 ≤ Complex.normSq (UpperHalfPlane.denom (Matrix.SpecialLinearGroup.toGLPos g) z)

@[simp]

theorem UpperHalfPlane.neg_smul (g : ↥(Matrix.GLPos (Fin 2) ℝ)) (z : UpperHalfPlane) : -g • z = g • z

@[simp]

theorem UpperHalfPlane.ModularGroup.sl_moeb (A : Matrix.SpecialLinearGroup (Fin 2) ℤ) (z : UpperHalfPlane) : A • z = ↑A • z

@[simp]

theorem UpperHalfPlane.ModularGroup.SL_neg_smul (g : Matrix.SpecialLinearGroup (Fin 2) ℤ) (z : UpperHalfPlane) : -g • z = g • z

theorem UpperHalfPlane.ModularGroup.im_smul_eq_div_normSq (g : Matrix.SpecialLinearGroup (Fin 2) ℤ) (z : UpperHalfPlane) : (g • z).im = z.im / Complex.normSq (UpperHalfPlane.denom (↑g) z)

theorem UpperHalfPlane.ModularGroup.denom_apply (g : Matrix.SpecialLinearGroup (Fin 2) ℤ) (z : UpperHalfPlane) : UpperHalfPlane.denom (↑g) z = ↑(↑g 1 0) * ↑z + ↑(↑g 1 1)

instance UpperHalfPlane.posRealAction : MulAction { x : ℝ // 0 < x } UpperHalfPlane

Equations

• UpperHalfPlane.posRealAction = MulAction.mk UpperHalfPlane.posRealAction.proof_2 UpperHalfPlane.posRealAction.proof_3

@[simp]

theorem UpperHalfPlane.coe_pos_real_smul (x : { x : ℝ // 0 < x }) (z : UpperHalfPlane) : ↑(x • z) = ↑x • ↑z

@[simp]

theorem UpperHalfPlane.pos_real_im (x : { x : ℝ // 0 < x }) (z : UpperHalfPlane) : (x • z).im = ↑x * z.im

@[simp]

theorem UpperHalfPlane.pos_real_re (x : { x : ℝ // 0 < x }) (z : UpperHalfPlane) : (x • z).re = ↑x * z.re

instance UpperHalfPlane.instAddActionReal : AddAction ℝ UpperHalfPlane

Equations

• UpperHalfPlane.instAddActionReal = AddAction.mk UpperHalfPlane.instAddActionReal.proof_2 UpperHalfPlane.instAddActionReal.proof_3

@[simp]

theorem UpperHalfPlane.coe_vadd (x : ℝ) (z : UpperHalfPlane) : ↑(x +ᵥ z) = ↑x + ↑z

@[simp]

theorem UpperHalfPlane.vadd_re (x : ℝ) (z : UpperHalfPlane) : (x +ᵥ z).re = x + z.re

@[simp]

theorem UpperHalfPlane.vadd_im (x : ℝ) (z : UpperHalfPlane) : (x +ᵥ z).im = z.im

theorem UpperHalfPlane.modular_S_smul (z : UpperHalfPlane) : ModularGroup.S • z = UpperHalfPlane.mk (-↑z)⁻¹ ⋯

theorem UpperHalfPlane.modular_T_zpow_smul (z : UpperHalfPlane) (n : ℤ) : ModularGroup.T ^ n • z = ↑n +ᵥ z

theorem UpperHalfPlane.modular_T_smul (z : UpperHalfPlane) : ModularGroup.T • z = 1 +ᵥ z

theorem UpperHalfPlane.exists_SL2_smul_eq_of_apply_zero_one_eq_zero (g : Matrix.SpecialLinearGroup (Fin 2) ℝ) (hc : ↑g 1 0 = 0) : ∃ (u : { x : ℝ // 0 < x }) (v : ℝ), (fun (x : UpperHalfPlane) => g • x) = (fun (x : UpperHalfPlane) => v +ᵥ x) ∘ fun (x : UpperHalfPlane) => u • x

theorem UpperHalfPlane.exists_SL2_smul_eq_of_apply_zero_one_ne_zero (g : Matrix.SpecialLinearGroup (Fin 2) ℝ) (hc : ↑g 1 0 ≠ 0) : ∃ (u : { x : ℝ // 0 < x }) (v : ℝ) (w : ℝ), (fun (x : UpperHalfPlane) => g • x) = (fun (x : UpperHalfPlane) => w +ᵥ x) ∘ (fun (x : UpperHalfPlane) => ModularGroup.S • x) ∘ (fun (x : UpperHalfPlane) => v +ᵥ x) ∘ fun (x : UpperHalfPlane) => u • x