Compare commits

..

No commits in common. "6a12258093f067e486abcb53e378629d5a2ae946" and "83f9aa5d725edbd84fb73222f021bb131b55a014" have entirely different histories.

3 changed files with 157 additions and 283 deletions

View File

@ -1,6 +1,6 @@
import Nevanlinna.complexHarmonic import Nevanlinna.complexHarmonic
import Nevanlinna.holomorphicAt import Nevanlinna.holomorphicAt
import Nevanlinna.holomorphic_primitive2 import Nevanlinna.holomorphic_primitive
import Nevanlinna.mathlibAddOn import Nevanlinna.mathlibAddOn
@ -108,175 +108,131 @@ function.
theorem harmonic_is_realOfHolomorphic theorem harmonic_is_realOfHolomorphic
{f : } {f : }
{z : } (hf : ∀ z, HarmonicAt f z) :
{R : } ∃ F : , (∀ z, HolomorphicAt F z) ∧ (Complex.reCLM ∘ F = f) := by
(hR : 0 < R)
(hf : ∀ x ∈ Metric.ball z R, HarmonicAt f x) :
∃ F : , (∀ x ∈ Metric.ball z R, HolomorphicAt F x) ∧ (Set.EqOn (Complex.reCLM ∘ F) f (Metric.ball z R)) := by
let f_1 : := Complex.ofRealCLM ∘ (partialDeriv 1 f) let f_1 : := Complex.ofRealCLM ∘ (partialDeriv 1 f)
let f_I : := Complex.ofRealCLM ∘ (partialDeriv Complex.I f) let f_I : := Complex.ofRealCLM ∘ (partialDeriv Complex.I f)
let g : := f_1 - Complex.I • f_I let g : := f_1 - Complex.I • f_I
have contDiffOn_if_contDiffAt have reg₂f : ContDiff 2 f := by
{f' : } apply contDiff_iff_contDiffAt.mpr
{z' : } intro z
{R' : } exact (hf z).1
{n' : }
(hf' : ∀ x ∈ Metric.ball z' R', ContDiffAt n' f' x) :
ContDiffOn n' f' (Metric.ball z' R') := by
intro z hz
apply ContDiffAt.contDiffWithinAt
exact hf' z hz
have reg₂f : ContDiffOn 2 f (Metric.ball z R) := by have reg₁f_1 : ContDiff 1 f_1 := by
apply contDiffOn_if_contDiffAt apply contDiff_iff_contDiffAt.mpr
intro x hx intro z
exact (hf x hx).1
have contDiffOn_if_contDiffAt'
{f' : }
{z' : }
{R' : }
{n' : }
(hf' : ∀ x ∈ Metric.ball z' R', ContDiffAt n' f' x) :
ContDiffOn n' f' (Metric.ball z' R') := by
intro z hz
apply ContDiffAt.contDiffWithinAt
exact hf' z hz
have reg₁f_1 : ContDiffOn 1 f_1 (Metric.ball z R) := by
apply contDiffOn_if_contDiffAt'
intro z hz
dsimp [f_1] dsimp [f_1]
apply ContDiffAt.continuousLinearMap_comp apply ContDiffAt.continuousLinearMap_comp
exact partialDeriv_contDiffAt (hf z hz).1 1 exact partialDeriv_contDiffAt (hf z).1 1
have reg₁f_I : ContDiffOn 1 f_I (Metric.ball z R) := by have reg₁f_I : ContDiff 1 f_I := by
apply contDiffOn_if_contDiffAt' apply contDiff_iff_contDiffAt.mpr
intro z hz intro z
dsimp [f_I] dsimp [f_I]
apply ContDiffAt.continuousLinearMap_comp apply ContDiffAt.continuousLinearMap_comp
exact partialDeriv_contDiffAt (hf z hz).1 Complex.I exact partialDeriv_contDiffAt (hf z).1 Complex.I
have reg₁g : ContDiffOn 1 g (Metric.ball z R) := by have reg₁g : ContDiff 1 g := by
dsimp [g] dsimp [g]
apply ContDiffOn.sub apply ContDiff.sub
exact reg₁f_1 exact reg₁f_1
have : Complex.I • f_I = fun x ↦ Complex.I • f_I x := by rfl apply ContDiff.const_smul'
rw [this]
apply ContDiffOn.const_smul
exact reg₁f_I exact reg₁f_I
have reg₁ : DifferentiableOn g (Metric.ball z R) := by have reg₁ : Differentiable g := by
intro x hx intro z
apply DifferentiableAt.differentiableWithinAt
apply CauchyRiemann₇.2 apply CauchyRiemann₇.2
constructor constructor
· apply DifferentiableWithinAt.differentiableAt (reg₁g.differentiableOn le_rfl x hx) · apply Differentiable.differentiableAt
apply IsOpen.mem_nhds Metric.isOpen_ball hx apply ContDiff.differentiable
exact reg₁g
rfl
· dsimp [g] · dsimp [g]
rw [partialDeriv_sub₂_differentiableAt, partialDeriv_sub₂_differentiableAt] rw [partialDeriv_sub₂, partialDeriv_sub₂]
simp
dsimp [f_1, f_I] dsimp [f_1, f_I]
rw [partialDeriv_smul'₂, partialDeriv_smul'₂] rw [partialDeriv_smul'₂, partialDeriv_smul'₂]
rw [partialDeriv_compContLinAt, partialDeriv_compContLinAt] rw [partialDeriv_compContLin, partialDeriv_compContLin, partialDeriv_compContLin, partialDeriv_compContLin]
simp
rw [partialDeriv_compContLinAt, partialDeriv_compContLinAt]
rw [mul_sub] rw [mul_sub]
simp simp
rw [← mul_assoc] rw [← mul_assoc]
simp simp
rw [add_comm, sub_eq_add_neg] rw [add_comm, sub_eq_add_neg]
congr 1 congr 1
· rw [partialDeriv_commOn _ reg₂f Complex.I 1] · rw [partialDeriv_comm reg₂f Complex.I 1]
exact hx · let A := Filter.EventuallyEq.eq_of_nhds (hf z).2
exact Metric.isOpen_ball
· let A := Filter.EventuallyEq.eq_of_nhds (hf x hx).2
simp at A simp at A
unfold Complex.laplace at A unfold Complex.laplace at A
conv => conv =>
right right
right right
rw [← sub_zero (partialDeriv 1 (partialDeriv 1 f) x)] rw [← sub_zero (partialDeriv 1 (partialDeriv 1 f) z)]
rw [← A] rw [← A]
simp simp
--DifferentiableAt (partialDeriv _ f) --Differentiable (partialDeriv _ f)
repeat repeat
apply ContDiffAt.differentiableAt apply ContDiff.differentiable
apply partialDeriv_contDiffAt (hf x hx).1 apply contDiff_iff_contDiffAt.mpr
exact fun w ↦ partialDeriv_contDiffAt (hf w).1 _
apply le_rfl apply le_rfl
-- DifferentiableAt f_1 x
apply (reg₁f_1.differentiableOn le_rfl).differentiableAt
apply IsOpen.mem_nhds Metric.isOpen_ball hx
-- DifferentiableAt (Complex.I • f_I)
have : Complex.I • f_I = fun x ↦ Complex.I • f_I x := by rfl
rw [this]
apply DifferentiableAt.const_smul
apply (reg₁f_I.differentiableOn le_rfl).differentiableAt
apply IsOpen.mem_nhds Metric.isOpen_ball hx
-- Differentiable f_1 -- Differentiable f_1
apply (reg₁f_1.differentiableOn le_rfl).differentiableAt exact reg₁f_1.differentiable le_rfl
apply IsOpen.mem_nhds Metric.isOpen_ball hx
-- Differentiable (Complex.I • f_I) -- Differentiable (Complex.I • f_I)
have : Complex.I • f_I = fun x ↦ Complex.I • f_I x := by rfl apply Differentiable.const_smul'
rw [this] exact reg₁f_I.differentiable le_rfl
apply DifferentiableAt.const_smul -- Differentiable f_1
apply (reg₁f_I.differentiableOn le_rfl).differentiableAt exact reg₁f_1.differentiable le_rfl
apply IsOpen.mem_nhds Metric.isOpen_ball hx -- Differentiable (Complex.I • f_I)
apply Differentiable.const_smul'
exact reg₁f_I.differentiable le_rfl
let F := fun z' ↦ (primitive z g) z' + f z let F := fun z ↦ (primitive 0 g) z + f 0
have regF : DifferentiableOn F (Metric.ball z R) := by have regF : Differentiable F := by
apply DifferentiableOn.add apply Differentiable.add
apply primitive_differentiableOn reg₁ apply primitive_differentiable reg₁
simp simp
have pF'' : ∀ x ∈ Metric.ball z R, (fderiv F x) = ContinuousLinearMap.lsmul (g x) := by have pF'' : ∀ x, (fderiv F x) = ContinuousLinearMap.lsmul (g x) := by
intro x hx intro x
have : DifferentiableAt F x := by rw [DifferentiableAt.fderiv_restrictScalars (regF x)]
apply (regF x hx).differentiableAt
apply IsOpen.mem_nhds Metric.isOpen_ball hx
rw [DifferentiableAt.fderiv_restrictScalars this]
dsimp [F] dsimp [F]
rw [fderiv_add_const] rw [fderiv_add_const]
rw [primitive_fderiv' reg₁ x hx] rw [primitive_fderiv']
exact rfl exact rfl
exact reg₁
use F use F
constructor constructor
· -- ∀ x ∈ Metric.ball z R, HolomorphicAt F x · -- ∀ (z : ), HolomorphicAt F z
intro x hx intro z
apply HolomorphicAt_iff.2 apply HolomorphicAt_iff.2
use Metric.ball z R use Set.univ
constructor constructor
· exact Metric.isOpen_ball · exact isOpen_const
· constructor · constructor
· assumption · simp
· intro w hw · intro w _
apply (regF w hw).differentiableAt exact regF w
apply IsOpen.mem_nhds Metric.isOpen_ball hw · -- (F z).re = f z
· -- Set.EqOn (⇑Complex.reCLM ∘ F) f (Metric.ball z R) have A := reg₂f.differentiable one_le_two
have : DifferentiableOn (Complex.reCLM ∘ F) (Metric.ball z R) := by have B : Differentiable (Complex.reCLM ∘ F) := by
apply DifferentiableOn.comp apply Differentiable.comp
apply Differentiable.differentiableOn exact ContinuousLinearMap.differentiable Complex.reCLM
apply ContinuousLinearMap.differentiable Complex.reCLM exact Differentiable.restrictScalars regF
apply regF.restrictScalars have C : (F 0).re = f 0 := by
exact Set.mapsTo'.mpr fun ⦃a⦄ _ => hR dsimp [F]
have hz : z ∈ Metric.ball z R := by exact Metric.mem_ball_self hR rw [primitive_zeroAtBasepoint]
apply Convex.eqOn_of_fderivWithin_eq _ this _ _ _ hz _ simp
exact convex_ball z R
apply reg₂f.differentiableOn one_le_two
apply IsOpen.uniqueDiffOn Metric.isOpen_ball
intro x hx apply eq_of_fderiv_eq B A _ 0 C
rw [fderivWithin_eq_fderiv, fderivWithin_eq_fderiv]
intro x
rw [fderiv.comp] rw [fderiv.comp]
simp simp
apply ContinuousLinearMap.ext apply ContinuousLinearMap.ext
@ -291,24 +247,7 @@ theorem harmonic_is_realOfHolomorphic
rw [(fderiv f x).map_smul, (fderiv f x).map_smul] rw [(fderiv f x).map_smul, (fderiv f x).map_smul]
rw [smul_eq_mul, smul_eq_mul] rw [smul_eq_mul, smul_eq_mul]
ring ring
-- DifferentiableAt (⇑Complex.reCLM) (F x)
assumption fun_prop
exact ContinuousLinearMap.differentiableAt Complex.reCLM -- DifferentiableAt F x
apply (regF.restrictScalars x hx).differentiableAt exact regF.restrictScalars x
apply IsOpen.mem_nhds Metric.isOpen_ball hx
apply IsOpen.uniqueDiffOn Metric.isOpen_ball
assumption
apply (reg₂f.differentiableOn one_le_two).differentiableAt
apply IsOpen.mem_nhds Metric.isOpen_ball hx
apply IsOpen.uniqueDiffOn Metric.isOpen_ball
assumption
-- DifferentiableAt (⇑Complex.reCLM ∘ F) x
apply DifferentiableAt.comp
apply Differentiable.differentiableAt
exact ContinuousLinearMap.differentiable Complex.reCLM
apply (regF.restrictScalars x hx).differentiableAt
apply IsOpen.mem_nhds Metric.isOpen_ball hx
--
dsimp [F]
rw [primitive_zeroAtBasepoint]
simp

View File

@ -49,6 +49,9 @@ theorem primitive_fderivAtBasepointZero
rw [this] rw [this]
obtain ⟨s, h₁s, h₂s⟩ : ∃ s ⊆ f⁻¹' Metric.ball (f 0) (c / (4 : )), IsOpen s ∧ 0 ∈ s := by obtain ⟨s, h₁s, h₂s⟩ : ∃ s ⊆ f⁻¹' Metric.ball (f 0) (c / (4 : )), IsOpen s ∧ 0 ∈ s := by
have B : Metric.ball (f 0) (c / 4) ∈ nhds (f 0) := by
apply Metric.ball_mem_nhds (f 0)
linarith
apply eventually_nhds_iff.mp apply eventually_nhds_iff.mp
apply continuousAt_def.1 apply continuousAt_def.1
apply Continuous.continuousAt apply Continuous.continuousAt
@ -121,106 +124,74 @@ theorem primitive_fderivAtBasepointZero
apply h₁s apply h₁s
exact h₂ε.1 hy exact h₂ε.1 hy
have t₀ {r : } (hr : r ∈ Metric.ball 0 ε) : IntervalIntegrable (fun x => f { re := x, im := 0 }) MeasureTheory.volume 0 r := by
apply ContinuousOn.intervalIntegrable
apply ContinuousOn.comp
exact hf
have : (fun x => ({ re := x, im := 0 } : )) = Complex.ofRealLI := by rfl
rw [this]
apply Continuous.continuousOn
continuity
intro x hx
apply h₂ε.2
simp
constructor
· simp
calc |x|
_ < ε := by
sorry
· simpa
have intervalComputation_uIcc {x' y' : } (h : x' ∈ Set.uIcc 0 y') : |x'| ≤ |y'| := by have t₁ {r : } (hr : r ∈ Metric.ball 0 ε) : IntervalIntegrable (fun _ => f 0) MeasureTheory.volume 0 r := by
let A := h.1 apply ContinuousOn.intervalIntegrable
let B := h.2 apply ContinuousOn.comp
rcases le_total 0 y' with hy | hy apply hf
· simp [hy] at A fun_prop
simp [hy] at B intro x hx
rwa [abs_of_nonneg A, abs_of_nonneg hy] simpa
· simp [hy] at A
simp [hy] at B
rw [abs_of_nonpos hy]
rw [abs_of_nonpos]
linarith [h.1]
exact B
have t₂ {a b : } : IntervalIntegrable (fun x_1 => f { re := a, im := x_1 }) MeasureTheory.volume 0 b := by
apply Continuous.intervalIntegrable
apply Continuous.comp hf
have : (Complex.mk a) = (fun x => Complex.I • Complex.ofRealCLM x + { re := a, im := 0 }) := by
funext x
apply Complex.ext
rw [Complex.add_re]
simp
simp
rw [this]
apply Continuous.add
continuity
fun_prop
have t₃ {a : } : IntervalIntegrable (fun _ => f 0) MeasureTheory.volume 0 a := by
apply Continuous.intervalIntegrable
apply Continuous.comp
exact hf
fun_prop
have {A B C D :E} : (A + B) - (C + D) = (A - C) + (B - D) := by
abel
conv =>
left
intro x
left
arg 1
rw [this]
rw [← smul_sub]
rw [← intervalIntegral.integral_sub t₀ t₁]
rw [← intervalIntegral.integral_sub t₂ t₃]
rw [Filter.eventually_iff_exists_mem] rw [Filter.eventually_iff_exists_mem]
use Metric.ball 0 (ε / (4 : ))
use Metric.ball 0 (ε / (4 : ))
constructor constructor
· apply Metric.ball_mem_nhds 0 · apply Metric.ball_mem_nhds 0
linarith linarith
· intro y hy · intro y hy
have {A B C D :E} : (A + B) - (C + D) = (A - C) + (B - D) := by
abel
rw [this]
rw [← smul_sub]
have t₀ : IntervalIntegrable (fun x => f { re := x, im := 0 }) MeasureTheory.volume 0 y.re := by
apply ContinuousOn.intervalIntegrable
apply ContinuousOn.comp
exact hf
have : (fun x => ({ re := x, im := 0 } : )) = Complex.ofRealLI := by rfl
rw [this]
apply Continuous.continuousOn
continuity
intro x hx
apply h₂ε.2
simp
constructor
· simp
calc |x|
_ ≤ |y.re| := by apply intervalComputation_uIcc hx
_ ≤ Complex.abs y := by exact Complex.abs_re_le_abs y
_ < ε / 4 := by simp at hy; assumption
_ < ε := by linarith
· simpa
have t₁ : IntervalIntegrable (fun _ => f 0) MeasureTheory.volume 0 y.re := by
apply ContinuousOn.intervalIntegrable
apply ContinuousOn.comp
apply hf
fun_prop
intro x _
simpa
rw [← intervalIntegral.integral_sub t₀ t₁]
have t₂ : IntervalIntegrable (fun x_1 => f { re := y.re, im := x_1 }) MeasureTheory.volume 0 y.im := by
apply ContinuousOn.intervalIntegrable
apply ContinuousOn.comp
exact hf
have : (Complex.mk y.re) = (fun x => Complex.I • Complex.ofRealCLM x + { re := y.re, im := 0 }) := by
funext x
apply Complex.ext
rw [Complex.add_re]
simp
simp
rw [this]
apply ContinuousOn.add
apply Continuous.continuousOn
continuity
fun_prop
intro x hx
apply h₂ε.2
constructor
· simp
calc |y.re|
_ ≤ Complex.abs y := by exact Complex.abs_re_le_abs y
_ < ε / 4 := by simp at hy; assumption
_ < ε := by linarith
· simp
calc |x|
_ ≤ |y.im| := by apply intervalComputation_uIcc hx
_ ≤ Complex.abs y := by exact Complex.abs_im_le_abs y
_ < ε / 4 := by simp at hy; assumption
_ < ε := by linarith
have t₃ : IntervalIntegrable (fun _ => f 0) MeasureTheory.volume 0 y.im := by
apply ContinuousOn.intervalIntegrable
apply ContinuousOn.comp
exact hf
fun_prop
intro x _
apply h₂ε.2
simp
constructor
· simpa
· simpa
rw [← intervalIntegral.integral_sub t₂ t₃]
have h₁y : |y.re| < ε / 4 := by have h₁y : |y.re| < ε / 4 := by
calc |y.re| calc |y.re|
_ ≤ Complex.abs y := by apply Complex.abs_re_le_abs _ ≤ Complex.abs y := by apply Complex.abs_re_le_abs
@ -781,11 +752,12 @@ theorem primitive_hasDerivAt
apply hasDerivAt_const apply hasDerivAt_const
theorem primitive_differentiableOn
theorem primitive_differentiable
{E : Type u} [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E] {E : Type u} [NormedAddCommGroup E] [NormedSpace E] [CompleteSpace E]
{f : → E} {f : → E}
{z₀ : } (z₀ : )
{R : } (R : )
(hf : DifferentiableOn f (Metric.ball z₀ R)) (hf : DifferentiableOn f (Metric.ball z₀ R))
: :
DifferentiableOn (primitive z₀ f) (Metric.ball z₀ R) := by DifferentiableOn (primitive z₀ f) (Metric.ball z₀ R) := by

View File

@ -38,10 +38,10 @@ theorem partialDeriv_eventuallyEq
exact fun v => rfl exact fun v => rfl
theorem partialDeriv_smul₁ theorem partialDeriv_smul₁
{f : E → F} {f : E → F}
{a : 𝕜} {a : 𝕜}
{v : E} : {v : E} :
partialDeriv 𝕜 (a • v) f = a • partialDeriv 𝕜 v f := by partialDeriv 𝕜 (a • v) f = a • partialDeriv 𝕜 v f := by
unfold partialDeriv unfold partialDeriv
conv => conv =>
@ -102,6 +102,20 @@ theorem partialDeriv_add₂ {f₁ f₂ : E → F} (h₁ : Differentiable 𝕜 f
simp simp
theorem partialDeriv_sub₂ {f₁ f₂ : E → F} (h₁ : Differentiable 𝕜 f₁) (h₂ : Differentiable 𝕜 f₂) : ∀ v : E, partialDeriv 𝕜 v (f₁ - f₂) = (partialDeriv 𝕜 v f₁) - (partialDeriv 𝕜 v f₂) := by
unfold partialDeriv
intro v
have : f₁ - f₂ = fun y ↦ f₁ y - f₂ y := by rfl
rw [this]
conv =>
left
intro w
left
rw [fderiv_sub (h₁ w) (h₂ w)]
funext w
simp
theorem partialDeriv_add₂_differentiableAt theorem partialDeriv_add₂_differentiableAt
{f₁ f₂ : E → F} {f₁ f₂ : E → F}
{v : E} {v : E}
@ -139,57 +153,6 @@ theorem partialDeriv_add₂_contDiffAt
exact (hf₂' x (Set.mem_of_mem_inter_right hx)).differentiableAt exact (hf₂' x (Set.mem_of_mem_inter_right hx)).differentiableAt
theorem partialDeriv_sub₂ {f₁ f₂ : E → F} (h₁ : Differentiable 𝕜 f₁) (h₂ : Differentiable 𝕜 f₂) : ∀ v : E, partialDeriv 𝕜 v (f₁ - f₂) = (partialDeriv 𝕜 v f₁) - (partialDeriv 𝕜 v f₂) := by
unfold partialDeriv
intro v
have : f₁ - f₂ = fun y ↦ f₁ y - f₂ y := by rfl
rw [this]
conv =>
left
intro w
left
rw [fderiv_sub (h₁ w) (h₂ w)]
funext w
simp
theorem partialDeriv_sub₂_differentiableAt
{f₁ f₂ : E → F}
{v : E}
{x : E}
(h₁ : DifferentiableAt 𝕜 f₁ x)
(h₂ : DifferentiableAt 𝕜 f₂ x) :
partialDeriv 𝕜 v (f₁ - f₂) x = (partialDeriv 𝕜 v f₁) x - (partialDeriv 𝕜 v f₂) x := by
unfold partialDeriv
have : f₁ - f₂ = fun y ↦ f₁ y - f₂ y := by rfl
rw [this]
rw [fderiv_sub h₁ h₂]
rfl
theorem partialDeriv_sub₂_contDiffAt
{f₁ f₂ : E → F}
{v : E}
{x : E}
(h₁ : ContDiffAt 𝕜 1 f₁ x)
(h₂ : ContDiffAt 𝕜 1 f₂ x) :
partialDeriv 𝕜 v (f₁ - f₂) =ᶠ[nhds x] (partialDeriv 𝕜 v f₁) - (partialDeriv 𝕜 v f₂) := by
obtain ⟨f₁', u₁, hu₁, _, hf₁'⟩ := contDiffAt_one_iff.1 h₁
obtain ⟨f₂', u₂, hu₂, _, hf₂'⟩ := contDiffAt_one_iff.1 h₂
apply Filter.eventuallyEq_iff_exists_mem.2
use u₁ ∩ u₂
constructor
· exact Filter.inter_mem hu₁ hu₂
· intro x hx
simp
apply partialDeriv_sub₂_differentiableAt 𝕜
exact (hf₁' x (Set.mem_of_mem_inter_left hx)).differentiableAt
exact (hf₂' x (Set.mem_of_mem_inter_right hx)).differentiableAt
theorem partialDeriv_compContLin theorem partialDeriv_compContLin
{f : E → F} {f : E → F}
{l : F →L[𝕜] G} {l : F →L[𝕜] G}