Documentation

TauCeti.Analysis.ODE.TrajectoryOperator

The bounded first-order trajectory operator #

The operator D_A = d/ds + A(s) acts from the first-order weak-derivative graph to Lᵖ for any essentially bounded, strongly measurable family of bounded operators on a real Banach space. We use the existing WeakDerivStep over the identity map on Lᵖ, so its norm controls both the value and its weak derivative. No new Sobolev carrier is introduced.

WeakDerivStep.trajectoryOperator is a continuous linear map. Its operator norm is at most 1 + ‖A‖, and changing the coefficient changes the operator by at most the L∞ distance between coefficients. Solving D_A u = f is equivalent to the weak ordinary differential equation u' = f - A u. For differentiable representatives with locally integrable derivative, the operator agrees almost everywhere with the classical expression u' + A u. In particular, for a path in W^{1,2}(ℝ, F) with a C¹ representative φ, the squared norms of the path, of its weak derivative and of its image under D_A are ∫ ‖φ‖², ∫ ‖φ'‖² and ∫ ‖φ' + A φ‖² (TauCeti.WeakDerivStep.norm_prev_sq_eq_integral, TauCeti.WeakDerivStep.norm_weakFDeriv_sq_eq_integral and TauCeti.WeakDerivStep.norm_trajectoryOperator_sq_eq_integral).

On the line, for a path with a C¹ representative φ, the squared norm of D_A u is the integral of ‖φ' + A φ‖² (TauCeti.WeakDerivStep.norm_trajectoryOperator_sq_eq_integral).

These statements supply the bounded realization to which coercive estimates and compact-error criteria can be applied. They do not assert Fredholmness, asymptotic invertibility, or density of smooth paths. Neither symmetry nor finite dimensionality of the target is required.

References #

The coefficient action is Mathlib's ContinuousLinearMap.holderL, applied to evaluation; the domain is the closed weak-derivative graph of TauCeti.Analysis.Sobolev.GraphStep.Basic.

The first-order operator u ↦ u' + A u on the weak-derivative graph over Lᵖ. The coefficient is an L∞ operator field, so the definition depends only on its almost everywhere equivalence class.

Instances For
    @[simp]

    The operator is the sum of the weak derivative in direction 1 and the coefficient action.

    theorem TauCeti.WeakDerivStep.trajectoryOperator_apply_ae {F : Type u_1} [NormedAddCommGroup F] [NormedSpace ℝ F] [CompleteSpace F] {mu : MeasureTheory.Measure ℝ} [MeasureTheory.IsFiniteMeasureOnCompacts mu] {Omega : TopologicalSpace.Opens ℝ} {p : ENNReal} [Fact (1 ≤ p)] (A : ↥(MeasureTheory.Lp (F →L[ℝ] F) ⊤ (mu.restrict ↑Omega))) (u : ↥(WeakDerivStep mu Omega p (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))))) :
    ↑↑((trajectoryOperator A) u) =ᵐ[mu.restrict ↑Omega] fun (s : ℝ) => (↑↑(weakFDeriv (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) u) s) 1 + (↑↑A s) (↑↑(prev (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) u) s)

    On representatives, the bounded trajectory operator has the expected pointwise formula almost everywhere.

    The operator norm of the first-order trajectory operator is at most 1 + ‖A‖.

    The weak derivative of a path is controlled by its image under the trajectory operator and its Lᵖ norm: ‖u'‖ ≤ ‖u' + A u‖ + ‖A‖ ‖u‖.

    Subtracting two trajectory operators cancels the derivative part and leaves multiplication by the difference of their coefficients.

    An L∞ perturbation of the coefficient changes the trajectory operator by at most the same amount in operator norm.

    The trajectory operator depends 1-Lipschitz continuously on its L∞ coefficient, for the operator norm.

    theorem TauCeti.WeakDerivStep.trajectoryOperator_eq_iff_hasWeakLineDerivOn {F : Type u_1} [NormedAddCommGroup F] [NormedSpace ℝ F] [CompleteSpace F] {mu : MeasureTheory.Measure ℝ} [MeasureTheory.IsFiniteMeasureOnCompacts mu] {Omega : TopologicalSpace.Opens ℝ} {p : ENNReal} [Fact (1 ≤ p)] (A : ↥(MeasureTheory.Lp (F →L[ℝ] F) ⊤ (mu.restrict ↑Omega))) (u : ↥(WeakDerivStep mu Omega p (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))))) (f : ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) :
    (trajectoryOperator A) u = f ↔ HasWeakLineDerivOn mu Omega (↑↑(prev (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) u)) (fun (s : ℝ) => ↑↑f s - (↑↑A s) (↑↑(prev (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) u) s)) 1

    Solving the bounded trajectory equation is equivalent to solving the weak ordinary differential equation u' = f - A u.

    theorem TauCeti.WeakDerivStep.exists_prev_eq_and_trajectoryOperator_eq_of_hasWeakLineDerivOn {F : Type u_1} [NormedAddCommGroup F] [NormedSpace ℝ F] [CompleteSpace F] {mu : MeasureTheory.Measure ℝ} [MeasureTheory.IsFiniteMeasureOnCompacts mu] {Omega : TopologicalSpace.Opens ℝ} {p : ENNReal} [Fact (1 ≤ p)] (A : ↥(MeasureTheory.Lp (F →L[ℝ] F) ⊤ (mu.restrict ↑Omega))) {x f : ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))} (h : HasWeakLineDerivOn mu Omega (↑↑x) (fun (s : ℝ) => ↑↑f s - (↑↑A s) (↑↑x s)) 1) :
    ∃ (u : ↥(WeakDerivStep mu Omega p (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))))), prev (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) u = x ∧ (trajectoryOperator A) u = f

    Weak solutions are values of graph elements. If an Lᵖ path x has weak derivative f - A x, then x is the value of an element u of the weak-derivative graph with D_A u = f. This is the converse of trajectoryOperator_eq_iff_hasWeakLineDerivOn: there the element u is given, here it is constructed from x.

    theorem TauCeti.WeakDerivStep.trajectoryOperator_apply_ae_eq_deriv_add {F : Type u_1} [NormedAddCommGroup F] [NormedSpace ℝ F] [CompleteSpace F] {mu : MeasureTheory.Measure ℝ} [mu.IsAddHaarMeasure] {Omega : TopologicalSpace.Opens ℝ} {p : ENNReal} [Fact (1 ≤ p)] (A : ↥(MeasureTheory.Lp (F →L[ℝ] F) ⊤ (mu.restrict ↑Omega))) (u : ↥(WeakDerivStep mu Omega p (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))))) {v : ℝ → F} (hv : ↑↑(prev (ContinuousLinearMap.id ℝ ↥(MeasureTheory.Lp F p (mu.restrict ↑Omega))) u) =ᵐ[mu.restrict ↑Omega] v) (hdv : MeasureTheory.LocallyIntegrableOn (fderiv ℝ v) (↑Omega) mu) (hvd : ∀ s ∈ ↑Omega, DifferentiableAt ℝ v s) :
    ↑↑((trajectoryOperator A) u) =ᵐ[mu.restrict ↑Omega] fun (s : ℝ) => deriv v s + (↑↑A s) (v s)

    For a differentiable representative whose derivative is locally integrable, the Sobolev trajectory operator agrees almost everywhere with the classical expression v' + A v.

    Norms of a path with a C¹ representative #

    On a path with a C¹ representative φ, the squared norm of D_A is ∫ ‖φ' + A φ‖².