TheoremBase

prNpr_N is Lipschitz from H to RmR^m, so w inherits (semi)continuity from Z (clause 1). For clauses 2-3, with a sign σ: w^∓_δ = w ∓ δh on V; moving x̂ along high modes leaves w fixed and h splits V-orthogonally, so ψ = −φ∘T ± δh∘Π_N has a penalised extremum, and the penalised-maximum lemma gives x̂ ∈ D(A) and Dφ(x̂) ± δAx̂ vanishing off Γ_N. Restricting φ ± δh to the slice ι_N z + (x̂ − Π_N x̂) gives a Euclidean test function at prNpr_N x̂; the Galerkin viscosity inequality translates term by term into the δ-shifted inequality for F^♯_N at y = x̂.

Proof

Each result cited below is universally quantified over the data in its own statement.

The notation of The Wick-Square Problem on the Sobolev Triple of Order Two: Standing Notation, Hilbert Triples: Standing Notation and Background and The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential is in force, and elementary facts about real numbers, finite sums, reindexing of finite sums along bijections, absolute values and squares are those carried by The Real Numbers: Standing Notation and Background. The items Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple, The δ\delta-Envelopes uδ−u^-_\delta and uδ+u^+_\delta of a Function on an Open Subset of a Hilbert Triple, Second-Order Equation Operator on an Open Subset of a Hilbert Triple and Its δ\delta-Shifts, Basic Properties of the δ\delta-Envelopes: Semicontinuity, Duality, Closed Superlevel Sets, Bounds and Monotonicity, First- and Second-Order Conditions at a Local Extremum of a C2C^2 Function Penalised by hh on the Small Space and The Penalty Function h=12∣⋅∣V2h=\tfrac12|\cdot|_V^2 of a Hilbert Triple: Expansion Identities, Closed Sublevel Sets, Density and Local Bounds are stated in the setting of Hilbert Triples: Standing Notation and Background; their notation (HH, VV, AA, D(A)D(A), hh, C2(H)C^{2}(H), Sym(H)\mathrm{Sym}(H), Sym(V)\mathrm{Sym}(V), Y∣VY|_{V}, IVI_{V}) is the notation fixed for the present triple by Hilbert Triples: Standing Notation and Background, and we apply them with U=HU=H, so that V∩U=VV\cap U=V and W=D(A)W=D(A), as recorded in Hilbert Triples: Standing Notation and Background §open-sets.

Since ΓN\Gamma_{N} is nonempty, m≥1m\ge1. For i∈[m]i\in[m] we write ki=κN(i)k_{i}=\kappa_{N}(i) and Ei=ekiE_{i}=e_{k_{i}}. Since κN\kappa_{N} is a bijection from [m][m] onto ΓN\Gamma_{N} (The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential §cutoff), ∑k∈ΓNf(k)=∑i=1mf(ki)\sum_{k\in\Gamma_{N}}f(k)=\sum_{i=1}^{m}f(k_{i}) for every real-valued ff on ΓN\Gamma_{N}, and ki≠kjk_{i}\neq k_{j} for i≠ji\neq j. Recall 1≤μk1\le\mu_{k} for every mode kk (The Wick-Square Problem on the Torus: Standing Notation §modes).

Step 1 (Coordinate identities). Every unit family eke_{k} lies in H−1=D(A)⊆V⊆HH^{-1}=D(A)\subseteq V\subseteq H (The Wick-Square Problem on the Torus: Standing Notation §units, The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §domain, Hilbert Triples: Standing Notation and Background §triple). By The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §pairings, for every z∈Hz\in H and every mode kk,

⟨z,ek⟩H=z(k)μk3,and, if z∈V,⟨z,ek⟩V=z(k)μk2.(1)\langle z,e_{k}\rangle_{H}=\frac{z(k)}{\mu_{k}^{3}},\qquad\text{and, if }z\in V,\qquad \langle z,e_{k}\rangle_{V}=\frac{z(k)}{\mu_{k}^{2}}.\qquad(1)

Taking z=ekz=e_{k} gives ∣ek∣H2=μk−3≤1|e_{k}|_{H}^{2}=\mu_{k}^{-3}\le1, hence ∣ek∣H≤1|e_{k}|_{H}\le1; and ⟨el,ek⟩V=el(k)/μk2\langle e_{l},e_{k}\rangle_{V}=e_{l}(k)/\mu_{k}^{2}, which is μk−2\mu_{k}^{-2} for l=kl=k and 00 for l≠kl\neq k.

(a) By The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential §coordinates, ιNζ=∑i=1mζiEi\iota_{N}\zeta=\sum_{i=1}^{m}\zeta_{i}E_{i} has value ζi\zeta_{i} at kik_{i} and 00 off ΓN\Gamma_{N}, so prNιNζ=ζ\mathrm{pr}_{N}\iota_{N}\zeta=\zeta for ζ∈Rm\zeta\in\mathbb{R}^{m}; and for every family xx, ΠNx=ιN(prNx)=∑i=1mx(ki)Ei\Pi_{N}x=\iota_{N}(\mathrm{pr}_{N}x)=\sum_{i=1}^{m}x(k_{i})E_{i} agrees with xx at each kik_{i} and vanishes off ΓN\Gamma_{N}. Hence prNΠN=prN\mathrm{pr}_{N}\Pi_{N}=\mathrm{pr}_{N} and ΠNΠN=ΠN\Pi_{N}\Pi_{N}=\Pi_{N}, and prN\mathrm{pr}_{N}, ιN\iota_{N}, ΠN\Pi_{N} are linear, the operations on families being pointwise. For x∈Hx\in H, ΠNx∈H−1⊆V\Pi_{N}x\in H^{-1}\subseteq V by The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §data. For a mode k∉ΓNk\notin\Gamma_{N}, eke_{k} vanishes at every kik_{i}, so ΠNek=0H\Pi_{N}e_{k}=0_{H}.

(b) Let ζ,ζ′∈Rm\zeta,\zeta'\in\mathbb{R}^{m} and let c∈Hc\in H vanish at every kik_{i}. By Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations and (1),

⟨ιNζ,c⟩H=∑i=1mζic(ki)μki3=0,⟨ιNζ,ιNζ′⟩V=∑i=1mζiζi′μki2,\langle\iota_{N}\zeta,c\rangle_{H}=\sum_{i=1}^{m}\zeta_{i}\frac{c(k_{i})}{\mu_{k_{i}}^{3}}=0,\qquad \langle\iota_{N}\zeta,\iota_{N}\zeta'\rangle_{V}=\sum_{i=1}^{m}\frac{\zeta_{i}\zeta'_{i}}{\mu_{k_{i}}^{2}},

and, if moreover c∈Vc\in V, ⟨ιNζ,c⟩V=∑iζic(ki)/μki2=0\langle\iota_{N}\zeta,c\rangle_{V}=\sum_{i}\zeta_{i}c(k_{i})/\mu_{k_{i}}^{2}=0. With The Penalty Function h=12∣⋅∣V2h=\tfrac12|\cdot|_V^2 of a Hilbert Triple: Expansion Identities, Closed Sublevel Sets, Density and Local Bounds §expansion-v, applied to the points ιNζ\iota_{N}\zeta and ιNζ+c\iota_{N}\zeta+c of VV, this gives, for c∈Vc\in V vanishing at every kik_{i},

h(ιNζ+c)=h(ιNζ)+h(c),h(ιNζ)=12∑i=1mζi2μki2.(2)h(\iota_{N}\zeta+c)=h(\iota_{N}\zeta)+h(c),\qquad h(\iota_{N}\zeta)=\tfrac12\sum_{i=1}^{m}\frac{\zeta_{i}^{2}}{\mu_{k_{i}}^{2}}.\qquad(2)

(c) Let a,a′∈Ha,a'\in H. By Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations and (1), ⟨ΠNa,a′⟩H=∑i=1ma(ki)a′(ki)/μki3\langle\Pi_{N}a,a'\rangle_{H}=\sum_{i=1}^{m}a(k_{i})a'(k_{i})/\mu_{k_{i}}^{3}, which is symmetric in a,a′a,a'; so ⟨ΠNa,a′⟩H=⟨a,ΠNa′⟩H\langle\Pi_{N}a,a'\rangle_{H}=\langle a,\Pi_{N}a'\rangle_{H}. The family a−ΠNaa-\Pi_{N}a vanishes at every kik_{i}, so by (b) ⟨ΠNa,a−ΠNa⟩H=0\langle\Pi_{N}a,a-\Pi_{N}a\rangle_{H}=0, and Elementary Identities in a Real Inner Product Space §expansion gives ∣a∣H2=∣ΠNa∣H2+∣a−ΠNa∣H2≥∣a−ΠNa∣H2|a|_{H}^{2}=|\Pi_{N}a|_{H}^{2}+|a-\Pi_{N}a|_{H}^{2}\ge|a-\Pi_{N}a|_{H}^{2}, whence ∣a−ΠNa∣H≤∣a∣H|a-\Pi_{N}a|_{H}\le|a|_{H}, both sides being nonnegative. For y∈Vy\in V, (2) with ζ=prNy\zeta=\mathrm{pr}_{N}y and c=y−ΠNyc=y-\Pi_{N}y gives

h(y)=h(ΠNy)+h(y−ΠNy).(3)h(y)=h(\Pi_{N}y)+h(y-\Pi_{N}y).\qquad(3)

(d) For a∈Ha\in H and a mode kk, (1) and The Cauchy-Schwarz Inequality in a Real Inner Product Space give ∣a(k)∣=μk3∣⟨a,ek⟩H∣≤μk3∣a∣H∣ek∣H≤μk3∣a∣H|a(k)|=\mu_{k}^{3}|\langle a,e_{k}\rangle_{H}|\le\mu_{k}^{3}|a|_{H}|e_{k}|_{H}\le\mu_{k}^{3}|a|_{H}. For ζ∈Rm\zeta\in\mathbb{R}^{m} the family ιNζ\iota_{N}\zeta has value ζi\zeta_{i} at kik_{i} by (a), so Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations and (1) give, as in (b),

∣ιNζ∣H2=⟨ιNζ,ιNζ⟩H=∑i=1mζiζiμki3=∑i=1mζi2μki3≤∑i=1mζi2=∥ζ∥2,|\iota_{N}\zeta|_{H}^{2}=\langle\iota_{N}\zeta,\iota_{N}\zeta\rangle_{H}=\sum_{i=1}^{m}\zeta_{i}\frac{\zeta_{i}}{\mu_{k_{i}}^{3}}=\sum_{i=1}^{m}\frac{\zeta_{i}^{2}}{\mu_{k_{i}}^{3}}\le\sum_{i=1}^{m}\zeta_{i}^{2}=\lVert\zeta\rVert^{2},

the inequality termwise because μki−3≤1\mu_{k_{i}}^{-3}\le1, and the last equality by Euclidean Norm on Rn\mathbb{R}^n. Both ∣ιNζ∣H|\iota_{N}\zeta|_{H} and ∥ζ∥\lVert\zeta\rVert being nonnegative,

∣ιNζ∣H≤∥ζ∥.(4)|\iota_{N}\zeta|_{H}\le\lVert\zeta\rVert .\qquad(4)

Put CN=∑i=1mμki3C_{N}=\sum_{i=1}^{m}\mu_{k_{i}}^{3}; its terms are at least 11 since 1≤μk1\le\mu_{k}, so μki3≤CN\mu_{k_{i}}^{3}\le C_{N} for each ii, and 1≤μk13≤CN1\le\mu_{k_{1}}^{3}\le C_{N} because 1≤m1\le m. For a∈Ha\in H, ιN(prNa)=ΠNa\iota_{N}(\mathrm{pr}_{N}a)=\Pi_{N}a, and the identity ∣a∣H2=∣ΠNa∣H2+∣a−ΠNa∣H2|a|_{H}^{2}=|\Pi_{N}a|_{H}^{2}+|a-\Pi_{N}a|_{H}^{2} of (c) gives ∣ΠNa∣H2≤∣a∣H2|\Pi_{N}a|_{H}^{2}\le|a|_{H}^{2}; so, by Euclidean Norm on Rn\mathbb{R}^n and the computation above with ζ=prNa\zeta=\mathrm{pr}_{N}a,

∥prNa∥2=∑i=1mμki3 a(ki)2μki3≤CN∑i=1ma(ki)2μki3=CN∣ΠNa∣H2≤CN∣a∣H2≤CN2∣a∣H2,\lVert\mathrm{pr}_{N}a\rVert^{2}=\sum_{i=1}^{m}\mu_{k_{i}}^{3}\,\frac{a(k_{i})^{2}}{\mu_{k_{i}}^{3}}\le C_{N}\sum_{i=1}^{m}\frac{a(k_{i})^{2}}{\mu_{k_{i}}^{3}}=C_{N}|\Pi_{N}a|_{H}^{2}\le C_{N}|a|_{H}^{2}\le C_{N}^{2}|a|_{H}^{2},

whence, both sides being nonnegative,

∥prNa∥≤CN∣a∣H(a∈H).(4′)\lVert\mathrm{pr}_{N}a\rVert\le C_{N}|a|_{H}\qquad(a\in H).\qquad(4')

Step 2 (Clause 1, and semicontinuity of ww). For a,a′∈Ha,a'\in H, prNa−prNa′=prN(a−a′)\mathrm{pr}_{N}a-\mathrm{pr}_{N}a'=\mathrm{pr}_{N}(a-a') by linearity of prN\mathrm{pr}_{N} (Step 1(a)), so (4') gives ∥prNa−prNa′∥≤CN∣a−a′∣H\lVert\mathrm{pr}_{N}a-\mathrm{pr}_{N}a'\rVert\le C_{N}|a-a'|_{H}. Given a real ε>0\varepsilon>0, every a′∈Ha'\in H with ∣a−a′∣H<ε/CN|a-a'|_{H}<\varepsilon/C_{N} thus satisfies dE(prNa,prNa′)<εd_{E}(\mathrm{pr}_{N}a,\mathrm{pr}_{N}a')<\varepsilon, so prN:H→Rm\mathrm{pr}_{N}:H\to\mathbb{R}^{m} is continuous on HH relative to HH in the sense of Continuous Map Between Metric Spaces. By Semicontinuity and Continuity Under Composition with a Continuous Map, applied with the metric spaces (Rm,dE)(\mathbb{R}^{m},d_{E}) and (H,dH)(H,d_{H}), with its map taken to be prN\mathrm{pr}_{N} on HH and its function taken to be ZZ on Rm\mathbb{R}^{m}: if ZZ is continuous on (Rm,dE)(\mathbb{R}^{m},d_{E}) then w=Z∘prNw=Z\circ\mathrm{pr}_{N} is continuous on HH (claim 3), which is continuity on HH in the sense of Real Hilbert Spaces: Standing Notation and Background §topology; this proves clause 1. Likewise, by claims 1 and 2 there, ww is upper semicontinuous on HH when ZZ is upper semicontinuous on Rm\mathbb{R}^{m}, and lower semicontinuous on HH when ZZ is lower semicontinuous on Rm\mathbb{R}^{m}.

Step 3 (Composition with an affine contraction). For a∈Ha\in H put La=a−ΠNaLa=a-\Pi_{N}a. By Step 1(a),(c), L:H→HL:H\to H is linear, ∣La∣H≤∣a∣H|La|_{H}\le|a|_{H}, and ⟨La,a′⟩H=⟨a,a′⟩H−⟨ΠNa,a′⟩H=⟨a,La′⟩H\langle La,a'\rangle_{H}=\langle a,a'\rangle_{H}-\langle\Pi_{N}a,a'\rangle_{H}=\langle a,La'\rangle_{H} for a,a′∈Ha,a'\in H. Let c∈Hc\in H, let χ∈C2(H)\chi\in C^{2}(H), and define χT:H→R\chi^{T}:H\to\mathbb{R} by χT(y)=χ(c+Ly)\chi^{T}(y)=\chi(c+Ly). We claim that χT∈C2(H)\chi^{T}\in C^{2}(H), with

DχT(y)=L Dχ(c+Ly),D2χT(y)=By,By(a,a′)=D2χ(c+Ly)(La,La′).D\chi^{T}(y)=L\,D\chi(c+Ly),\qquad D^{2}\chi^{T}(y)=B_{y},\quad B_{y}(a,a')=D^{2}\chi(c+Ly)(La,La').

Fix y∈Hy\in H and write x=c+Lyx=c+Ly, so c+L(y+v)=x+Lvc+L(y+v)=x+Lv for v∈Hv\in H. Differentiability, in the sense of Fréchet Differentiability and the Gradient on an Open Subset of a Real Inner Product Space §differentiable: given a real ε>0\varepsilon>0, let η>0\eta>0 be a radius for χ\chi at xx and ε\varepsilon; for v∈Hv\in H with ∣v∣H<η|v|_{H}<\eta one has ∣Lv∣H≤∣v∣H<η|Lv|_{H}\le|v|_{H}<\eta, hence

∣χT(y+v)−χT(y)−⟨LDχ(x),v⟩H∣=∣χ(x+Lv)−χ(x)−⟨Dχ(x),Lv⟩H∣≤ε∣Lv∣H≤ε∣v∣H.\bigl|\chi^{T}(y+v)-\chi^{T}(y)-\langle LD\chi(x),v\rangle_{H}\bigr|=\bigl|\chi(x+Lv)-\chi(x)-\langle D\chi(x),Lv\rangle_{H}\bigr|\le\varepsilon|Lv|_{H}\le\varepsilon|v|_{H}.

So χT\chi^{T} is differentiable at yy with gradient LDχ(x)LD\chi(x), and by uniqueness of the gradient (Fréchet Differentiability and the Gradient on an Open Subset of a Real Inner Product Space §gradient) DχT(y)=LDχ(x)D\chi^{T}(y)=LD\chi(x); as yy was arbitrary, χT\chi^{T} is differentiable on HH and DχT(y′)=LDχ(c+Ly′)D\chi^{T}(y')=LD\chi(c+Ly') for every y′∈Hy'\in H. Continuity of the gradient map: ∣LDχ(c+Ly′)−LDχ(c+Ly)∣H≤∣Dχ(c+Ly′)−Dχ(c+Ly)∣H|LD\chi(c+Ly')-LD\chi(c+Ly)|_{H}\le|D\chi(c+Ly')-D\chi(c+Ly)|_{H} and ∣(c+Ly′)−(c+Ly)∣H≤∣y′−y∣H|(c+Ly')-(c+Ly)|_{H}\le|y'-y|_{H}, so continuity of DχD\chi at xx (part of The Classes C1C^1 and C2C^2 on an Open Subset of a Real Inner Product Space §c1) yields continuity of DχTD\chi^{T} at yy with the same radius. The form ByB_{y}: symmetry, additivity and homogeneity follow from those of D2χ(x)D^{2}\chi(x) and linearity of LL, and ∣By(a,a′)∣≤∥D2χ(x)∥∣La∣H∣La′∣H≤∥D2χ(x)∥∣a∣H∣a′∣H|B_{y}(a,a')|\le\lVert D^{2}\chi(x)\rVert|La|_{H}|La'|_{H}\le\lVert D^{2}\chi(x)\rVert|a|_{H}|a'|_{H} by Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §bound; so By∈Sym(H)B_{y}\in\mathrm{Sym}(H) by Bounded Symmetric Bilinear Forms on a Real Inner Product Space: Norm, Order, Identity Form and Restriction §form. Second derivative, in the sense of The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §expansion with ρ=1\rho=1 (every ball lies in HH and χT\chi^{T} is differentiable on HH): given ε>0\varepsilon>0, let δ′′≤1\delta''\le1 be a radius witnessing that D2χ(x)D^{2}\chi(x) is a second derivative of χ\chi at xx (with ρ=1\rho=1, admissible since the condition does not depend on ρ\rho). For v,a∈Hv,a\in H with ∣v∣H<δ′′|v|_{H}<\delta'', ∣Lv∣H<δ′′|Lv|_{H}<\delta'', and by self-adjointness of LL,

∣⟨DχT(y+v)−DχT(y),a⟩H−By(v,a)∣=∣⟨Dχ(x+Lv)−Dχ(x),La⟩H−D2χ(x)(Lv,La)∣≤ε∣Lv∣H∣La∣H≤ε∣v∣H∣a∣H.\bigl|\langle D\chi^{T}(y+v)-D\chi^{T}(y),a\rangle_{H}-B_{y}(v,a)\bigr|=\bigl|\langle D\chi(x+Lv)-D\chi(x),La\rangle_{H}-D^{2}\chi(x)(Lv,La)\bigr|\le\varepsilon|Lv|_{H}|La|_{H}\le\varepsilon|v|_{H}|a|_{H}.

So ByB_{y} is a second derivative of χT\chi^{T} at yy, and by uniqueness of the Hessian (The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §hessian) D2χT(y)=ByD^{2}\chi^{T}(y)=B_{y}. Continuity of the Hessian map: for y′∈Hy'\in H and a,a′∈Ha,a'\in H, ∣(By′−By)(a,a′)∣=∣(D2χ(c+Ly′)−D2χ(x))(La,La′)∣≤∥D2χ(c+Ly′)−D2χ(x)∥ ∣a∣H∣a′∣H|(B_{y'}-B_{y})(a,a')|=|(D^{2}\chi(c+Ly')-D^{2}\chi(x))(La,La')|\le\lVert D^{2}\chi(c+Ly')-D^{2}\chi(x)\rVert\,|a|_{H}|a'|_{H} by Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §bound, so by the second part of that claim dSym(By′,By)≤dSym(D2χ(c+Ly′),D2χ(x))d_{\mathrm{Sym}}(B_{y'},B_{y})\le d_{\mathrm{Sym}}(D^{2}\chi(c+Ly'),D^{2}\chi(x)), and continuity of D2χD^{2}\chi at xx (part of The Classes C1C^1 and C2C^2 on an Open Subset of a Real Inner Product Space §c2) gives continuity of y↦Byy\mapsto B_{y}. This proves the claim.

Step 4 (Restriction of a C2C^{2} function to a coordinate slice). Let y0∈Hy_{0}\in H and χ∈C2(H)\chi\in C^{2}(H), put S(z)=ιNz+y0S(z)=\iota_{N}z+y_{0} for z∈Rmz\in\mathbb{R}^{m} and χ0(z)=χ(S(z))\chi_{0}(z)=\chi(S(z)). By (4) and linearity of ιN\iota_{N}, ∣S(z′)−S(z)∣H≤∥z′−z∥|S(z')-S(z)|_{H}\le\lVert z'-z\rVert. We claim that χ0\chi_{0} is of class C2C^{2} on Rm\mathbb{R}^{m} (C^k Maps on a Euclidean Open Set), with

∂iχ0(z)=⟨Dχ(S(z)),Ei⟩H,∂l∂iχ0(z)=D2χ(S(z))(El,Ei)(i,l∈[m]).\partial_{i}\chi_{0}(z)=\langle D\chi(S(z)),E_{i}\rangle_{H},\qquad \partial_{l}\partial_{i}\chi_{0}(z)=D^{2}\chi(S(z))(E_{l},E_{i})\qquad(i,l\in[m]).

Continuity: if Θ:H→R\Theta:H\to\mathbb{R} is continuous at S(z)S(z) with radius η\eta for a given ε\varepsilon, then Θ∘S\Theta\circ S is continuous at zz with the same radius η\eta. This applies to χ\chi (continuous by Basic Properties of Differentiability on an Open Subset of a Real Inner Product Space §continuous), to x↦⟨Dχ(x),Ei⟩Hx\mapsto\langle D\chi(x),E_{i}\rangle_{H} (as ∣⟨Dχ(x′)−Dχ(x),Ei⟩H∣≤∣Dχ(x′)−Dχ(x)∣H|\langle D\chi(x')-D\chi(x),E_{i}\rangle_{H}|\le|D\chi(x')-D\chi(x)|_{H} by The Cauchy-Schwarz Inequality in a Real Inner Product Space and ∣Ei∣H≤1|E_{i}|_{H}\le1, and DχD\chi is continuous) and to x↦D2χ(x)(El,Ei)x\mapsto D^{2}\chi(x)(E_{l},E_{i}) (as ∣(D2χ(x′)−D2χ(x))(El,Ei)∣≤∥D2χ(x′)−D2χ(x)∥|(D^{2}\chi(x')-D^{2}\chi(x))(E_{l},E_{i})|\le\lVert D^{2}\chi(x')-D^{2}\chi(x)\rVert by Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §continuity, and D2χD^{2}\chi is continuous). First partials, in the sense of Partial Derivative on a Euclidean Open Set: moving the iith coordinate of zz by tt moves S(z)S(z) to S(z)+tEiS(z)+tE_{i}. Given ε>0\varepsilon>0, let η>0\eta>0 be a radius for differentiability of χ\chi at S(z)S(z) with gradient Dχ(S(z))D\chi(S(z)) and ε/2\varepsilon/2, the gradient Dχ(S(z))D\chi(S(z)) being by definition the unique vector with this property (Fréchet Differentiability and the Gradient on an Open Subset of a Real Inner Product Space §gradient); for 0<∣t∣<η0<|t|<\eta, ∣tEi∣H≤∣t∣<η|tE_{i}|_{H}\le|t|<\eta, so ∣χ(S(z)+tEi)−χ(S(z))−t⟨Dχ(S(z)),Ei⟩H∣≤ε2∣t∣|\chi(S(z)+tE_{i})-\chi(S(z))-t\langle D\chi(S(z)),E_{i}\rangle_{H}|\le\tfrac{\varepsilon}{2}|t|, and dividing by ∣t∣|t| the difference quotient of χ0\chi_{0} is within ε2<ε\tfrac{\varepsilon}{2}<\varepsilon of ⟨Dχ(S(z)),Ei⟩H\langle D\chi(S(z)),E_{i}\rangle_{H}. Second partials: with ωi(z)=⟨Dχ(S(z)),Ei⟩H\omega_{i}(z)=\langle D\chi(S(z)),E_{i}\rangle_{H} one has ωi(z′)−ωi(z)=⟨Dχ(S(z)+tEl)−Dχ(S(z)),Ei⟩H\omega_{i}(z')-\omega_{i}(z)=\langle D\chi(S(z)+tE_{l})-D\chi(S(z)),E_{i}\rangle_{H} when z′z' is zz with llth coordinate moved by tt; given ε>0\varepsilon>0, a radius δ′′\delta'' witnessing The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §expansion for χ\chi at S(z)S(z) with ε/2\varepsilon/2 and the form D2χ(S(z))D^{2}\chi(S(z)), which is by definition the unique second derivative of χ\chi at S(z)S(z) (The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §hessian), used with its increment tEltE_{l} and direction EiE_{i}, gives for 0<∣t∣<δ′′0<|t|<\delta'' the bound ∣ωi(z′)−ωi(z)−t D2χ(S(z))(El,Ei)∣≤ε2∣t∣|\omega_{i}(z')-\omega_{i}(z)-t\,D^{2}\chi(S(z))(E_{l},E_{i})|\le\tfrac{\varepsilon}{2}|t|, so the difference quotient of ωi\omega_{i} is within ε2<ε\tfrac{\varepsilon}{2}<\varepsilon of D2χ(S(z))(El,Ei)D^{2}\chi(S(z))(E_{l},E_{i}). Thus χ0\chi_{0} is continuous with continuous first partials, i.e. of class C1C^{1}, and each ∂iχ0=ωi\partial_{i}\chi_{0}=\omega_{i} is continuous with continuous first partials, i.e. of class C1C^{1}; by clause 2 of C^k Maps on a Euclidean Open Set (with k=1k=1) and its scalar convention, χ0\chi_{0} is of class C2C^{2} on Rm\mathbb{R}^{m}, and the formulas hold.

Step 5 (Set-up for clauses 2 and 3). Let σ∈{1,−1}\sigma\in\{1,-1\}: in case σ=1\sigma=1 assume that ZZ is a viscosity subsolution of FN♯\mathcal{F}^{\sharp}_{N} on Rm\mathbb{R}^{m}, and in case σ=−1\sigma=-1 that it is a viscosity supersolution, both in the sense of The Galerkin Hamilton-Jacobi-Bellman Equations of the Wick-Square Problem in the Coordinates of a Cube of Modes §equation, that is of Viscosity Subsolution and Supersolution of a Second-Order Equation for the open set Rm\mathbb{R}^{m}. Then ZZ is upper semicontinuous on Rm\mathbb{R}^{m} if σ=1\sigma=1 and lower semicontinuous if σ=−1\sigma=-1, so by Step 2 ww is upper, respectively lower, semicontinuous on HH. By Basic Properties of the δ\delta-Envelopes: Semicontinuity, Duality, Closed Superlevel Sets, Bounds and Monotonicity §semicontinuous-case, ww is bounded above (if σ=1\sigma=1), respectively below (if σ=−1\sigma=-1), near each point of HH; this one-sided local bound comes from the semicontinuity of ww alone, and no bound on ZZ is assumed or used anywhere in this proof. So the standing hypothesis of Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §subsolution, respectively Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §supersolution, holds; and for every real δ>0\delta>0 the envelope wδσw^{\sigma}_{\delta}, meaning wδ−w^{-}_{\delta} if σ=1\sigma=1 and wδ+w^{+}_{\delta} if σ=−1\sigma=-1 (The δ\delta-Envelopes uδ−u^-_\delta and uδ+u^+_\delta of a Function on an Open Subset of a Hilbert Triple §minus, The δ\delta-Envelopes uδ−u^-_\delta and uδ+u^+_\delta of a Function on an Open Subset of a Hilbert Triple §plus), satisfies

wδσ(x)=w(x)−σδh(x)(x∈V).w^{\sigma}_{\delta}(x)=w(x)-\sigma\delta h(x)\qquad(x\in V).

The operator FN♯F^{\sharp}_{N} is a second-order equation operator on HH relative to (H,V,A)(H,V,A) by The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §operator. Write FσF^{\sigma} for its shift (FN♯)δ−(F^{\sharp}_{N})^{-}_{\delta} if σ=1\sigma=1 and (FN♯)δ+(F^{\sharp}_{N})^{+}_{\delta} if σ=−1\sigma=-1; by Second-Order Equation Operator on an Open Subset of a Hilbert Triple and Its δ\delta-Shifts §shifted, in both cases

Fσ(x,r,p,Y)=FN♯(x, r+σδh(x), p+σδAx, Y∣V+σδIV)F^{\sigma}(x,r,p,Y)=F^{\sharp}_{N}\bigl(x,\ r+\sigma\delta h(x),\ p+\sigma\delta Ax,\ Y|_{V}+\sigma\delta I_{V}\bigr)

for x∈D(A)x\in D(A), r∈Rr\in\mathbb{R}, p∈Hp\in H and Y∈Sym(H)Y\in\mathrm{Sym}(H), where Y∣V+(−δ)IVY|_{V}+(-\delta)I_{V} is the form Y∣V−δIVY|_{V}-\delta I_{V}.

Fix a real δ>0\delta>0, φ∈C2(H)\varphi\in C^{2}(H), a point x^∈V\hat{x}\in V at which wδσ−φw^{\sigma}_{\delta}-\varphi has a local maximum (if σ=1\sigma=1), respectively a local minimum (if σ=−1\sigma=-1), relative to VV, and a real ε>0\varepsilon>0. By Local Maximum of a Function Relative to a Subset of a Metric Space and Local Minimum of a Function Relative to a Subset of a Metric Space, with the metric dHd_{H} (Hilbert Triples: Standing Notation and Background §open-sets), there is a real r>0r>0 such that every x∈Vx\in V with ∣x−x^∣H<r|x-\hat{x}|_{H}<r satisfies

σ[(w(x)−σδh(x)−φ(x))−(w(x^)−σδh(x^)−φ(x^))]≤0.(5)\sigma\Bigl[\bigl(w(x)-\sigma\delta h(x)-\varphi(x)\bigr)-\bigl(w(\hat{x})-\sigma\delta h(\hat{x})-\varphi(\hat{x})\bigr)\Bigr]\le0.\qquad(5)

Put z^=prNx^∈Rm\hat{z}=\mathrm{pr}_{N}\hat{x}\in\mathbb{R}^{m} and y0=x^−ΠNx^=Lx^∈Vy_{0}=\hat{x}-\Pi_{N}\hat{x}=L\hat{x}\in V; then x^=ιNz^+y0\hat{x}=\iota_{N}\hat{z}+y_{0}, y0y_{0} vanishes at every kik_{i}, and z^i=x^(ki)\hat{z}_{i}=\hat{x}(k_{i}).

Step 6 (The point x^\hat{x} lies in D(A)D(A), and the shifted gradient has no high modes). Define T:H→HT:H\to H by Ty=ΠNx^+LyTy=\Pi_{N}\hat{x}+Ly. For y∈Vy\in V: Ty∈VTy\in V; prNTy=prNx^+prNy−prNy=z^\mathrm{pr}_{N}Ty=\mathrm{pr}_{N}\hat{x}+\mathrm{pr}_{N}y-\mathrm{pr}_{N}y=\hat{z} by Step 1(a), so w(Ty)=Z(z^)=w(x^)w(Ty)=Z(\hat{z})=w(\hat{x}); Ty=ιNz^+LyTy=\iota_{N}\hat{z}+Ly with Ly∈VLy\in V vanishing at every kik_{i}, so (2) gives h(Ty)=h(ΠNx^)+h(Ly)h(Ty)=h(\Pi_{N}\hat{x})+h(Ly); Tx^=x^T\hat{x}=\hat{x}; and Ty−x^=L(y−x^)Ty-\hat{x}=L(y-\hat{x}), so ∣Ty−x^∣H≤∣y−x^∣H|Ty-\hat{x}|_{H}\le|y-\hat{x}|_{H}.

Let bNb_{N} be the map assigning to y,y′∈Hy,y'\in H the number bN(y,y′)=δ∑i=1my(ki)y′(ki)/μki2b_{N}(y,y')=\delta\sum_{i=1}^{m}y(k_{i})y'(k_{i})/\mu_{k_{i}}^{2}. It is symmetric, additive and homogeneous in each argument (the operations being pointwise), and by Step 1(d) ∣bN(y,y′)∣≤(δ∑i=1mμki4)∣y∣H∣y′∣H|b_{N}(y,y')|\le\bigl(\delta\sum_{i=1}^{m}\mu_{k_{i}}^{4}\bigr)|y|_{H}|y'|_{H}; so bN∈Sym(H)b_{N}\in\mathrm{Sym}(H). Let Q(y)=12bN(y,y)Q(y)=\tfrac12 b_{N}(y,y). By Affine and Quadratic Functions on a Real Hilbert Space are of Class C2C^2 §form, Q∈C2(H)Q\in C^{2}(H) with DQ(y)=TbNyDQ(y)=T_{b_{N}}y, where ⟨TbNy,a⟩H=bN(y,a)\langle T_{b_{N}}y,a\rangle_{H}=b_{N}(y,a) for a∈Ha\in H (The Bounded Symmetric Operator Represented by a Bounded Symmetric Bilinear Form on a Real Hilbert Space §representation). By (2) with ζ=prNy\zeta=\mathrm{pr}_{N}y, Q(y)=δh(ΠNy)Q(y)=\delta h(\Pi_{N}y) for y∈Hy\in H. Put

ψσ(y)=−φ(Ty)+σQ(y)(y∈H).\psi_{\sigma}(y)=-\varphi(Ty)+\sigma Q(y)\qquad(y\in H).

By Step 3 (with c=ΠNx^c=\Pi_{N}\hat{x}, χ=φ\chi=\varphi) and Constants, Sums, Scalar Multiples and Differences of Differentiable Functions on an Open Subset of a Real Inner Product Space §scalar, Constants, Sums, Scalar Multiples and Differences of Differentiable Functions on an Open Subset of a Real Inner Product Space §sum, ψσ∈C2(H)\psi_{\sigma}\in C^{2}(H) and, as Tx^=x^T\hat{x}=\hat{x},

Dψσ(x^)=−L Dφ(x^)+σTbNx^.(6)D\psi_{\sigma}(\hat{x})=-L\,D\varphi(\hat{x})+\sigma T_{b_{N}}\hat{x}.\qquad(6)

For y∈Vy\in V, Q(y)=δh(ΠNy)Q(y)=\delta h(\Pi_{N}y) and (3) give ψσ(y)−σδh(y)=−φ(Ty)−σδh(Ly)\psi_{\sigma}(y)-\sigma\delta h(y)=-\varphi(Ty)-\sigma\delta h(Ly), and by the properties of TT above this equals (w(Ty)−σδh(Ty)−φ(Ty))−w(x^)+σδh(ΠNx^)\bigl(w(Ty)-\sigma\delta h(Ty)-\varphi(Ty)\bigr)-w(\hat{x})+\sigma\delta h(\Pi_{N}\hat{x}). Subtracting the same identity at y=x^y=\hat{x} (where Tx^=x^T\hat{x}=\hat{x}),

(ψσ(y)−σδh(y))−(ψσ(x^)−σδh(x^))=(w−σδh−φ)(Ty)−(w−σδh−φ)(x^).\bigl(\psi_{\sigma}(y)-\sigma\delta h(y)\bigr)-\bigl(\psi_{\sigma}(\hat{x})-\sigma\delta h(\hat{x})\bigr)=\bigl(w-\sigma\delta h-\varphi\bigr)(Ty)-\bigl(w-\sigma\delta h-\varphi\bigr)(\hat{x}).

If y∈Vy\in V and ∣y−x^∣H<r|y-\hat{x}|_{H}<r, then Ty∈VTy\in V and ∣Ty−x^∣H<r|Ty-\hat{x}|_{H}<r, so (5) at x=Tyx=Ty shows that σ\sigma times the left side is at most 00. Thus for σ=1\sigma=1 the function ψ1−δh\psi_{1}-\delta h has a local maximum at x^\hat{x} relative to VV, and First- and Second-Order Conditions at a Local Extremum of a C2C^2 Function Penalised by hh on the Small Space §maximum (with λ=δ\lambda=\delta) gives x^∈D(A)\hat{x}\in D(A) and δAx^=Dψ1(x^)\delta A\hat{x}=D\psi_{1}(\hat{x}); for σ=−1\sigma=-1 the function ψ−1+δh\psi_{-1}+\delta h has a local minimum at x^\hat{x} relative to VV, and First- and Second-Order Conditions at a Local Extremum of a C2C^2 Function Penalised by hh on the Small Space §minimum gives x^∈D(A)\hat{x}\in D(A) and δAx^=−Dψ−1(x^)\delta A\hat{x}=-D\psi_{-1}(\hat{x}). In both cases x^∈D(A)\hat{x}\in D(A) and σδAx^=Dψσ(x^)\sigma\delta A\hat{x}=D\psi_{\sigma}(\hat{x}). Put

pσ=Dφ(x^)+σδAx^∈H.p_{\sigma}=D\varphi(\hat{x})+\sigma\delta A\hat{x}\in H .

Let k∉ΓNk\notin\Gamma_{N}. Then Lek=ekLe_{k}=e_{k} (Step 1(a)), so ⟨LDφ(x^),ek⟩H=⟨Dφ(x^),ek⟩H\langle LD\varphi(\hat{x}),e_{k}\rangle_{H}=\langle D\varphi(\hat{x}),e_{k}\rangle_{H} by self-adjointness of LL, and ⟨TbNx^,ek⟩H=bN(x^,ek)=0\langle T_{b_{N}}\hat{x},e_{k}\rangle_{H}=b_{N}(\hat{x},e_{k})=0 because eke_{k} vanishes at every kik_{i}. Pairing σδAx^=Dψσ(x^)\sigma\delta A\hat{x}=D\psi_{\sigma}(\hat{x}) with eke_{k} and using (6) gives σδ⟨Ax^,ek⟩H=−⟨Dφ(x^),ek⟩H\sigma\delta\langle A\hat{x},e_{k}\rangle_{H}=-\langle D\varphi(\hat{x}),e_{k}\rangle_{H}, i.e. ⟨pσ,ek⟩H=0\langle p_{\sigma},e_{k}\rangle_{H}=0, so by (1)

pσ(k)=0(k∉ΓN).(7)p_{\sigma}(k)=0\qquad(k\notin\Gamma_{N}).\qquad(7)

Moreover, by The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §domain and (1), ⟨Ax^,Ei⟩H=μkix^(ki)/μki3=z^i/μki2\langle A\hat{x},E_{i}\rangle_{H}=\mu_{k_{i}}\hat{x}(k_{i})/\mu_{k_{i}}^{3}=\hat{z}_{i}/\mu_{k_{i}}^{2} for i∈[m]i\in[m], hence

⟨Dφ(x^),Ei⟩H+σδz^iμki2=⟨pσ,Ei⟩H=pσ(ki)μki3.(8)\langle D\varphi(\hat{x}),E_{i}\rangle_{H}+\sigma\delta\frac{\hat{z}_{i}}{\mu_{k_{i}}^{2}}=\langle p_{\sigma},E_{i}\rangle_{H}=\frac{p_{\sigma}(k_{i})}{\mu_{k_{i}}^{3}}.\qquad(8)

Step 7 (The Euclidean test function). With S(z)=ιNz+y0∈VS(z)=\iota_{N}z+y_{0}\in V as in Step 4, define Ψ:Rm→R\Psi:\mathbb{R}^{m}\to\mathbb{R} by Ψ(z)=φ(S(z))+σδh(S(z))\Psi(z)=\varphi(S(z))+\sigma\delta h(S(z)). Since y0∈Vy_{0}\in V vanishes at every kik_{i}, (2) gives

Ψ(z)=φ0(z)+σδh(y0)+σδ2∑i=1mzi2μki2,φ0=φ∘S.\Psi(z)=\varphi_{0}(z)+\sigma\delta h(y_{0})+\frac{\sigma\delta}{2}\sum_{i=1}^{m}\frac{z_{i}^{2}}{\mu_{k_{i}}^{2}},\qquad \varphi_{0}=\varphi\circ S .

By Step 4, φ0\varphi_{0} is of class C2C^{2} on Rm\mathbb{R}^{m}; the remaining terms form a function built from constants and coordinate functions by sums, scalar multiples and products, which is smooth on Rm\mathbb{R}^{m} by claims 2 and 3 of Constants, Coordinate Functions, Sums and Products of CkC^k Functions on a Euclidean Open Set, and hence of class C2C^{2} by Smooth Map on a Euclidean Open Set; so Ψ\Psi is of class C2C^{2} on Rm\mathbb{R}^{m} by claim 3 of Constants, Coordinate Functions, Sums and Products of CkC^k Functions on a Euclidean Open Set. By claim 1 there (the iith partial of the llth coordinate function being 11 if l=il=i and 00 otherwise, directly from Partial Derivative on a Euclidean Open Set) and Step 4, for z∈Rmz\in\mathbb{R}^{m} and i,l∈[m]i,l\in[m],

∂iΨ(z)=⟨Dφ(S(z)),Ei⟩H+σδziμki2,∂l∂iΨ(z)=D2φ(S(z))(El,Ei)+{σδμki−2(l=i),0(l≠i).\partial_{i}\Psi(z)=\langle D\varphi(S(z)),E_{i}\rangle_{H}+\sigma\delta\frac{z_{i}}{\mu_{k_{i}}^{2}},\qquad \partial_{l}\partial_{i}\Psi(z)=D^{2}\varphi(S(z))(E_{l},E_{i})+\begin{cases}\sigma\delta\mu_{k_{i}}^{-2}&(l=i),\\0&(l\neq i).\end{cases}

For z∈Rmz\in\mathbb{R}^{m}, prNS(z)=z+prNy0=z\mathrm{pr}_{N}S(z)=z+\mathrm{pr}_{N}y_{0}=z, so w(S(z))=Z(z)w(S(z))=Z(z) and

Z(z)−Ψ(z)=(w−σδh−φ)(S(z)),S(z^)=x^,∣S(z)−x^∣H≤∥z−z^∥Z(z)-\Psi(z)=\bigl(w-\sigma\delta h-\varphi\bigr)(S(z)),\qquad S(\hat{z})=\hat{x},\qquad |S(z)-\hat{x}|_{H}\le\lVert z-\hat{z}\rVert

by (4). Hence every zz with dE(z,z^)<rd_{E}(z,\hat{z})<r satisfies σ[(Z−Ψ)(z)−(Z−Ψ)(z^)]≤0\sigma[(Z-\Psi)(z)-(Z-\Psi)(\hat{z})]\le0 by (5): Z−ΨZ-\Psi has a local maximum (if σ=1\sigma=1), respectively minimum (if σ=−1\sigma=-1), at z^\hat{z} relative to Rm\mathbb{R}^{m}. The defining property in Viscosity Subsolution and Supersolution of a Second-Order Equation, applied to the test function Ψ\Psi, gives

σ FN♯(z^,Z(z^),DΨ(z^),D2Ψ(z^))≤0.(9)\sigma\,\mathcal{F}^{\sharp}_{N}\bigl(\hat{z},Z(\hat{z}),D\Psi(\hat{z}),D^{2}\Psi(\hat{z})\bigr)\le0.\qquad(9)

Step 8 (Translation term by term). Let s=wδσ(x^)=w(x^)−σδh(x^)s=w^{\sigma}_{\delta}(\hat{x})=w(\hat{x})-\sigma\delta h(\hat{x}) and Xσ=D2φ(x^)∣V+σδIVX_{\sigma}=D^{2}\varphi(\hat{x})|_{V}+\sigma\delta I_{V}. Since x^∈D(A)\hat{x}\in D(A) (Step 6), Step 5 gives

Fσ(x^,s,Dφ(x^),D2φ(x^))=FN♯(x^,w(x^),pσ,Xσ),F^{\sigma}\bigl(\hat{x},s,D\varphi(\hat{x}),D^{2}\varphi(\hat{x})\bigr)=F^{\sharp}_{N}\bigl(\hat{x},w(\hat{x}),p_{\sigma},X_{\sigma}\bigr),

and w(x^)=Z(z^)w(\hat{x})=Z(\hat{z}). We compare the terms of The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §operator with those of The Galerkin Hamilton-Jacobi-Bellman Equations of the Wick-Square Problem in the Coordinates of a Cube of Modes §shifted, using S(z^)=x^S(\hat{z})=\hat{x}.

Gradient. By Step 7 and (8), ∂iΨ(z^)=⟨pσ,Ei⟩H\partial_{i}\Psi(\hat{z})=\langle p_{\sigma},E_{i}\rangle_{H} for every ii. The gradient DΨ(z^)D\Psi(\hat{z}) is the vector of the partials ∂iΨ(z^)\partial_{i}\Psi(\hat{z}) (Differential Calculus and Convexity on Euclidean Open Sets: Standing Notation §derivatives, Gradient of a Real-Valued Function on a Euclidean Open Set), and the squared Euclidean norm of a vector is the sum of the squares of its coordinates (Euclidean Norm on Rn\mathbb{R}^n). Hence 12∥DΨ(z^)∥2=12∑i=1m⟨pσ,Ei⟩H2=12∑k∈ΓN⟨pσ,ek⟩H2=12GN(pσ,pσ)\tfrac12\lVert D\Psi(\hat{z})\rVert^{2}=\tfrac12\sum_{i=1}^{m}\langle p_{\sigma},E_{i}\rangle_{H}^{2}=\tfrac12\sum_{k\in\Gamma_{N}}\langle p_{\sigma},e_{k}\rangle_{H}^{2}=\tfrac12 G_{N}(p_{\sigma},p_{\sigma}), with GNG_{N} as in The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §data.

Drift. By The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential §functions and claim 1 of Constants, Coordinate Functions, Sums and Products of CkC^k Functions on a Euclidean Open Set, ∂iΦN(z^)=(μki+2qki)z^i\partial_{i}\Phi_{N}(\hat{z})=(\mu_{k_{i}}+2q_{k_{i}})\hat{z}_{i}; since DΦN(z^)D\Phi_{N}(\hat{z}) and DΨ(z^)D\Psi(\hat{z}) are the vectors of partials (Gradient of a Real-Valued Function on a Euclidean Open Set) and the dot product is the sum of the coordinatewise products (Difference, Dot Product, and Orthogonality in Rn\mathbb{R}^n), and z^i=x^(ki)\hat{z}_{i}=\hat{x}(k_{i}),

DΦN(z^)⋅DΨ(z^)=∑i=1m(μki+2qki)x^(ki) ⟨pσ,Ei⟩H.D\Phi_{N}(\hat{z})\cdot D\Psi(\hat{z})=\sum_{i=1}^{m}\bigl(\mu_{k_{i}}+2q_{k_{i}}\bigr)\hat{x}(k_{i})\,\langle p_{\sigma},E_{i}\rangle_{H}.

On the other side, x^∈H−1\hat{x}\in H^{-1} and pσ∈Hp_{\sigma}\in H, so by The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §pairings ⟨Ax^,pσ⟩H\langle A\hat{x},p_{\sigma}\rangle_{H} is the lattice sum of k↦x^(k)pσ(k)/μk2k\mapsto\hat{x}(k)p_{\sigma}(k)/\mu_{k}^{2}. By (7) this family vanishes off ΓN\Gamma_{N}; for M≥NM\ge N one has ΓN⊆ΓM\Gamma_{N}\subseteq\Gamma_{M} (The Wick-Square Problem on the Torus: Standing Notation §cubes), so the MMth cube sum equals ∑k∈ΓNx^(k)pσ(k)/μk2\sum_{k\in\Gamma_{N}}\hat{x}(k)p_{\sigma}(k)/\mu_{k}^{2}, and by Cube Sums of Families on the Integer Lattice §lattice-sum so does the lattice sum. With (8),

⟨Ax^,pσ⟩H=∑i=1mx^(ki)pσ(ki)μki2=∑i=1mμkix^(ki)⟨pσ,Ei⟩H.\langle A\hat{x},p_{\sigma}\rangle_{H}=\sum_{i=1}^{m}\frac{\hat{x}(k_{i})p_{\sigma}(k_{i})}{\mu_{k_{i}}^{2}}=\sum_{i=1}^{m}\mu_{k_{i}}\hat{x}(k_{i})\langle p_{\sigma},E_{i}\rangle_{H}.

By The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §data and Step 1(a), BqNx^=ΠN(Bqx^)=∑i=1m2qkix^(ki)EiB^{N}_{q}\hat{x}=\Pi_{N}(B_{q}\hat{x})=\sum_{i=1}^{m}2q_{k_{i}}\hat{x}(k_{i})E_{i}, so by Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations ⟨BqNx^,pσ⟩H=∑i=1m2qkix^(ki)⟨Ei,pσ⟩H\langle B^{N}_{q}\hat{x},p_{\sigma}\rangle_{H}=\sum_{i=1}^{m}2q_{k_{i}}\hat{x}(k_{i})\langle E_{i},p_{\sigma}\rangle_{H}. Adding, by Elementary Identities in a Real Inner Product Space §bilinear,

⟨Ax^+BqNx^,pσ⟩H=DΦN(z^)⋅DΨ(z^).\langle A\hat{x}+B^{N}_{q}\hat{x},p_{\sigma}\rangle_{H}=D\Phi_{N}(\hat{z})\cdot D\Psi(\hat{z}).

Trace. The Hessian D2Ψ(z^)D^{2}\Psi(\hat{z}) is the matrix with entries ∂i∂lΨ(z^)\partial_{i}\partial_{l}\Psi(\hat{z}) (Differential Calculus and Convexity on Euclidean Open Sets: Standing Notation §derivatives, Hessian Matrix of a C^2 Function), so by Trace of a Real Square Matrix its trace is ∑i=1m∂i∂iΨ(z^)\sum_{i=1}^{m}\partial_{i}\partial_{i}\Psi(\hat{z}), and by Step 7, tr⁡D2Ψ(z^)=∑i=1m(D2φ(x^)(Ei,Ei)+σδμki−2)=∑i=1mXσ(Ei,Ei)\operatorname{tr}D^{2}\Psi(\hat{z})=\sum_{i=1}^{m}\bigl(D^{2}\varphi(\hat{x})(E_{i},E_{i})+\sigma\delta\mu_{k_{i}}^{-2}\bigr)=\sum_{i=1}^{m}X_{\sigma}(E_{i},E_{i}), since Xσ(Ei,Ei)=D2φ(x^)(Ei,Ei)+σδ⟨Ei,Ei⟩VX_{\sigma}(E_{i},E_{i})=D^{2}\varphi(\hat{x})(E_{i},E_{i})+\sigma\delta\langle E_{i},E_{i}\rangle_{V} by Bounded Symmetric Bilinear Forms on a Real Inner Product Space: Norm, Order, Identity Form and Restriction §identity and Bounded Symmetric Bilinear Forms on a Real Inner Product Space: Norm, Order, Identity Form and Restriction §restriction, and ⟨Ei,Ei⟩V=μki−2\langle E_{i},E_{i}\rangle_{V}=\mu_{k_{i}}^{-2} by Step 1. On the other side, with κ\kappa the enumeration of The Wick-Square Problem on the Sobolev Triple of Order Two: Standing Notation §noise, a bijection of N\mathbb{N} onto Zn\mathbb{Z}^{n} (an enumeration being by definition such a bijection, Properties of the Fourier Coefficients on the Torus, and the Realisation of Weighted Coefficient Families), the set J={j∈N:κ(j)∈ΓN}J=\{j\in\mathbb{N}:\kappa(j)\in\Gamma_{N}\} is finite and nonempty and κ\kappa maps it bijectively onto ΓN\Gamma_{N}; the terms Xσ(fjN,fjN)X_{\sigma}(f^{N}_{j},f^{N}_{j}) vanish for j∉Jj\notin J, since then fjN=0Hf^{N}_{j}=0_{H} and XσX_{\sigma} is bilinear. So every partial sum of index at least max⁡J\max J of the series defining TrfNXσ\mathrm{Tr}_{f^{N}}X_{\sigma} in Square-Summable Sequences in the Form Space of a Hilbert Triple and the Trace of a Form along Them §trace equals ∑j∈JXσ(eκ(j),eκ(j))=∑k∈ΓNXσ(ek,ek)=∑i=1mXσ(Ei,Ei)\sum_{j\in J}X_{\sigma}(e_{\kappa(j)},e_{\kappa(j)})=\sum_{k\in\Gamma_{N}}X_{\sigma}(e_{k},e_{k})=\sum_{i=1}^{m}X_{\sigma}(E_{i},E_{i}), and the series converges to this number. Hence tr⁡D2Ψ(z^)=TrfNXσ\operatorname{tr}D^{2}\Psi(\hat{z})=\mathrm{Tr}_{f^{N}}X_{\sigma}.

Running cost. ιNz^=ΠNx^∈H−1\iota_{N}\hat{z}=\Pi_{N}\hat{x}\in H^{-1}, so the running cost g∣H−1g|_{H^{-1}} of FN♯\mathcal{F}^{\sharp}_{N} takes the value g(ΠNx^)g(\Pi_{N}\hat{x}) at ιNz^\iota_{N}\hat{z}.

Combining these identities,

Fσ(x^,s,Dφ(x^),D2φ(x^))=γZ(z^)−ν2TrfNXσ+12GN(pσ,pσ)+⟨Ax^+BqNx^,pσ⟩H−g(ΠNx^)=FN♯(z^,Z(z^),DΨ(z^),D2Ψ(z^)),F^{\sigma}\bigl(\hat{x},s,D\varphi(\hat{x}),D^{2}\varphi(\hat{x})\bigr)=\gamma Z(\hat{z})-\tfrac{\nu}{2}\mathrm{Tr}_{f^{N}}X_{\sigma}+\tfrac12G_{N}(p_{\sigma},p_{\sigma})+\langle A\hat{x}+B^{N}_{q}\hat{x},p_{\sigma}\rangle_{H}-g(\Pi_{N}\hat{x})=\mathcal{F}^{\sharp}_{N}\bigl(\hat{z},Z(\hat{z}),D\Psi(\hat{z}),D^{2}\Psi(\hat{z})\bigr),

and therefore, by (9), σFσ(x^,s,Dφ(x^),D2φ(x^))≤0\sigma F^{\sigma}(\hat{x},s,D\varphi(\hat{x}),D^{2}\varphi(\hat{x}))\le0.

Step 9 (Conclusion). Take y=x^y=\hat{x}, which lies in D(A)=WD(A)=W by Step 6, the number s=wδσ(x^)s=w^{\sigma}_{\delta}(\hat{x}) of Step 8, the vector Dφ(x^)∈HD\varphi(\hat{x})\in H in the gradient slot of Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple and the form Y=D2φ(x^)∈Sym(H)Y=D^{2}\varphi(\hat{x})\in\mathrm{Sym}(H). Then ∣y−x^∣H=0|y-\hat{x}|_{H}=0, ∣wδσ(y)−wδσ(x^)∣=0|w^{\sigma}_{\delta}(y)-w^{\sigma}_{\delta}(\hat{x})|=0, ∣s−wδσ(x^)∣=0|s-w^{\sigma}_{\delta}(\hat{x})|=0, ∣Dφ(x^)−Dφ(x^)∣H=0|D\varphi(\hat{x})-D\varphi(\hat{x})|_{H}=0 and ∥Y−D2φ(x^)∥=∥0Sym∥=0\lVert Y-D^{2}\varphi(\hat{x})\rVert=\lVert0_{\mathrm{Sym}}\rVert=0 (Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §norm-axioms), all less than ε\varepsilon. If σ=1\sigma=1, Step 8 gives (FN♯)δ−(y,s,Dφ(x^),Y)≤0<ε(F^{\sharp}_{N})^{-}_{\delta}(y,s,D\varphi(\hat{x}),Y)\le0<\varepsilon; if σ=−1\sigma=-1, it gives −ε<0≤(FN♯)δ+(y,s,Dφ(x^),Y)-\varepsilon<0\le(F^{\sharp}_{N})^{+}_{\delta}(y,s,D\varphi(\hat{x}),Y). Since δ\delta, φ\varphi, x^\hat{x} and ε\varepsilon were arbitrary, ww is a viscosity subsolution of FN♯F^{\sharp}_{N} on HH when ZZ is a viscosity subsolution of FN♯\mathcal{F}^{\sharp}_{N} on Rm\mathbb{R}^{m}, by Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §subsolution, which is clause 2; and ww is a viscosity supersolution of FN♯F^{\sharp}_{N} on HH when ZZ is a viscosity supersolution of FN♯\mathcal{F}^{\sharp}_{N} on Rm\mathbb{R}^{m}, by Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §supersolution, which is clause 3.

Citations

Loading…

Dependencies

Uses0

Loading…

Comments

Log in to comment.

Loading…