Documentation

LeanPool.NavierStokesAndEuler.Euler.TransverseEndpointEnergy

The actual endpoint energy in the activation argument. Paths are genuine initial-zero Bochner H¹ paths, and the zero-terminal correction is solved in the existing closed transverse derivative space. Symmetry, positivity and minimum energy are conclusions of the construction.

The endpoint Schur complement of a coercive quadratic form. The stationary extension is constructed by the inverse of the form on the closed zero-trace space. No stationary extension or Dirichlet-to-Neumann map is an input.

noncomputable def EulerDirichletEndpointReduction.correction {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) :
E →L[] S

The zero-trace correction obtained by a genuine coercive inverse.

Equations
Instances For
    noncomputable def EulerDirichletEndpointReduction.stationaryPart {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) :

    Subtract the solved zero-trace correction from any trial extension.

    Equations
    Instances For
      theorem EulerDirichletEndpointReduction.stationaryPart_eq {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (x : E) :
      (stationaryPart S A c hc hA) x = x - ((correction S A c hc hA) x)
      theorem EulerDirichletEndpointReduction.correction_equation {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (x : E) (v : S) :
      inner (A ((correction S A c hc hA) x)) v = inner (A x) v
      theorem EulerDirichletEndpointReduction.stationaryPart_orthogonal {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (x : E) (v : S) :
      inner (A ((stationaryPart S A c hc hA) x)) v = 0

      The constructed extension satisfies all zero-trace stationary equations.

      theorem EulerDirichletEndpointReduction.stationaryPart_sub_mem {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (x : E) :
      (stationaryPart S A c hc hA) x - x S
      theorem EulerDirichletEndpointReduction.stationaryPart_subspace {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (v : S) :
      (stationaryPart S A c hc hA) v = 0
      theorem EulerDirichletEndpointReduction.stationaryPart_eq_of_sub_mem {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (x y : E) (hxy : x - y S) :
      (stationaryPart S A c hc hA) x = (stationaryPart S A c hc hA) y

      The actual extension depends only on the terminal class of the trial.

      theorem EulerDirichletEndpointReduction.energy_split {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (hAs : (↑A).IsSymmetric) (x : E) (v : S) :
      inner (A ((stationaryPart S A c hc hA) x + v)) ((stationaryPart S A c hc hA) x + v) = inner (A ((stationaryPart S A c hc hA) x)) ((stationaryPart S A c hc hA) x) + inner (A v) v

      Energy splits orthogonally along the stationary extension and zero-trace variations.

      theorem EulerDirichletEndpointReduction.stationaryPart_minimizes {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (hAs : (↑A).IsSymmetric) (x : E) :
      inner (A ((stationaryPart S A c hc hA) x)) ((stationaryPart S A c hc hA) x) inner (A x) x

      The solved extension minimizes the actual quadratic form in its trace class.

      noncomputable def EulerDirichletEndpointReduction.endpointExtension {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) (L : U →L[] E) :

      A prescribed bounded trial lift followed by the actual stationary projection.

      Equations
      Instances For
        noncomputable def EulerDirichletEndpointReduction.endpointOperator {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) [CompleteSpace E] (L : U →L[] E) :

        The operator representing the actual stationary endpoint energy.

        Equations
        • One or more equations did not get rendered due to their size.
        Instances For
          theorem EulerDirichletEndpointReduction.endpointOperator_inner {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) [CompleteSpace E] (L : U →L[] E) (x y : U) :
          inner ((endpointOperator S A c hc hA L) x) y = inner (A ((endpointExtension S A c hc hA L) x)) ((endpointExtension S A c hc hA L) y)
          theorem EulerDirichletEndpointReduction.endpointOperator_positive {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) [CompleteSpace E] (hAs : (↑A).IsSymmetric) (L : U →L[] E) :

          Positivity is inherited from the actual displacement form.

          theorem EulerDirichletEndpointReduction.endpointOperator_trial_bound {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) [CompleteSpace E] (hAs : (↑A).IsSymmetric) (L : U →L[] E) (x : U) :
          inner ((endpointOperator S A c hc hA L) x) x inner (A (L x)) (L x)

          Any explicit admissible trial controls the endpoint quadratic form.

          theorem EulerDirichletEndpointReduction.positive_norm_le_of_quadratic {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (B : E →L[] E) (hB : B.IsPositive) (C : ) (hC : 0 C) (hb : ∀ (x : E), inner (B x) x C * x ^ 2) :

          A nonnegative quadratic-form bound gives the same operator-norm bound.

          theorem EulerDirichletEndpointReduction.endpointOperator_norm_le {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (S : Submodule E) [CompleteSpace S] (A : E →L[] E) (c : ) (hc : 0 < c) (hA : ∀ (x : E), c * x ^ 2 inner (A x) x) [CompleteSpace E] (hAs : (↑A).IsSymmetric) (L : U →L[] E) (C : ) (hC : 0 C) (hL : ∀ (x : U), inner (A (L x)) (L x) C * x ^ 2) :
          endpointOperator S A c hc hA L C

          The endpoint norm is controlled by the energy of the trial, without an inverse norm loss.

          The physical kinetic-minus-potential form on all initial-zero H¹ paths.

          Equations
          • One or more equations did not get rendered due to their size.
          Instances For
            theorem EulerTransverseEndpointEnergy.energyOperator_coercive {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (u : (EulerTimeLp.TimeLp T E)) :
            1 / 2 * u ^ 2 inner ((energyOperator T hT H) u) u

            The source upper Hessian bound and time smallness give actual coercivity, also when the terminal value is nonzero.

            theorem EulerTransverseEndpointEnergy.timeMultiplier_symmetric {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] (T : ) (hT : 0 T) (H : C((Set.Icc 0 T), E →L[] E)) (hH : ∀ (t : (Set.Icc 0 T)), (↑(H t)).IsSymmetric) :
            theorem EulerTransverseEndpointEnergy.energyOperator_symmetric {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] (T : ) (hT : 0 T) (H : C((Set.Icc 0 T), E →L[] E)) (hH : ∀ (t : (Set.Icc 0 T)), (↑(H t)).IsSymmetric) :
            noncomputable def EulerTransverseEndpointEnergy.endpointDerivative {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) :

            Solve the zero-endpoint variation problem for an explicit terminal trial lift.

            Equations
            • One or more equations did not get rendered due to their size.
            Instances For
              noncomputable def EulerTransverseEndpointEnergy.endpointDisplacement {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) :
              U →L[] C((Set.Icc 0 T), E)

              The constructed physical stationary path.

              Equations
              • One or more equations did not get rendered due to their size.
              Instances For
                noncomputable def EulerTransverseEndpointEnergy.dirichletToNeumann {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) :

                The genuine endpoint quadratic form represented by a bounded operator.

                Equations
                • One or more equations did not get rendered due to their size.
                Instances For
                  theorem EulerTransverseEndpointEnergy.endpointDisplacement_initial {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) (Y : U) :
                  ((endpointDisplacement T hT m H K hK hH hsmall L) Y) 0, = 0
                  theorem EulerTransverseEndpointEnergy.endpointDerivative_sub_mem {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) (Y : U) :
                  theorem EulerTransverseEndpointEnergy.endpointDisplacement_terminal {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) (Y : U) :
                  ((endpointDisplacement T hT m H K hK hH hsmall L) Y) T, = ((EulerInitialTimePrimitive.initialPrimitive T hT) (L Y)) T,

                  The stationary correction preserves the actual terminal trace.

                  theorem EulerTransverseEndpointEnergy.endpointDisplacement_tangent {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) (hL : ∀ (Y : U) (t : (Set.Icc 0 T)), inner (m t) (((EulerInitialTimePrimitive.initialPrimitive T hT) (L Y)) t) = 0) (Y : U) (t : (Set.Icc 0 T)) :
                  inner (m t) (((endpointDisplacement T hT m H K hK hH hsmall L) Y) t) = 0

                  A genuinely tangent trial produces a genuinely tangent stationary path.

                  theorem EulerTransverseEndpointEnergy.endpointDerivative_weak {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) (Y : U) (v : (EulerTransverseVariationalInverse.transverseDerivatives T hT m)) :

                  Every zero-endpoint transverse test satisfies the literal stationary weak equation.

                  theorem EulerTransverseEndpointEnergy.dirichletToNeumann_inner {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (L : U →L[] (EulerTimeLp.TimeLp T E)) (Y Z : U) :
                  inner ((dirichletToNeumann T hT m H K hK hH hsmall L) Y) Z = inner ((endpointDerivative T hT m H K hK hH hsmall L) Y) ((endpointDerivative T hT m H K hK hH hsmall L) Z) - inner ((EulerTimeLp.timeMultiplier T hT H) ((EulerInitialTimePrimitive.initialPrimitiveTimeLp T hT) ((endpointDerivative T hT m H K hK hH hsmall L) Y))) ((EulerInitialTimePrimitive.initialPrimitiveTimeLp T hT) ((endpointDerivative T hT m H K hK hH hsmall L) Z))

                  The endpoint pairing is the actual kinetic-minus-potential energy pairing.

                  theorem EulerTransverseEndpointEnergy.dirichletToNeumann_positive {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (hHs : ∀ (t : (Set.Icc 0 T)), (↑(H t)).IsSymmetric) (L : U →L[] (EulerTimeLp.TimeLp T E)) :
                  (dirichletToNeumann T hT m H K hK hH hsmall L).IsPositive

                  The source endpoint operator is symmetric positive semidefinite.

                  theorem EulerTransverseEndpointEnergy.dirichletToNeumann_norm_le {E : Type u_1} {U : Type u_2} [NormedAddCommGroup E] [InnerProductSpace E] [CompleteSpace E] [NormedAddCommGroup U] [InnerProductSpace U] [CompleteSpace U] (T : ) (hT : 0 T) (m : (Set.Icc 0 T)E) (H : C((Set.Icc 0 T), E →L[] E)) (K : ) (hK : 0 K) (hH : ∀ (t : (Set.Icc 0 T)) (x : E), inner ((H t) x) x K * x ^ 2) (hsmall : K * (T ^ 2 / 2) 1 / 2) (hHs : ∀ (t : (Set.Icc 0 T)), (↑(H t)).IsSymmetric) (L : U →L[] (EulerTimeLp.TimeLp T E)) (C : ) (hC : 0 C) (hL : ∀ (Y : U), L Y ^ 2 - inner ((EulerTimeLp.timeMultiplier T hT H) ((EulerInitialTimePrimitive.initialPrimitiveTimeLp T hT) (L Y))) ((EulerInitialTimePrimitive.initialPrimitiveTimeLp T hT) (L Y)) C * Y ^ 2) :
                  dirichletToNeumann T hT m H K hK hH hsmall L C

                  The endpoint norm costs only trial energy, with no inverse or frame norm factor.