选择一段文字或公式,即可在这里留下建议。保存后的建议仅作者可见,审核后才会公开。
多重遍历平均 / 原文精读
第三单元(续):四个位置的循环逼近与最早两个增量
本篇展开 D Lemma 5.2,原文 pp.29–30,式 (5.7)–(5.9)。前篇已在 H 局部检测的前提下构造各精度的投影逼近。这里首先证明一个有限的四线性恒等式,再给出每对晚精度函数的相关估计,并检查终端余项确实趋于零。
每个原始区间 I 上,所有精度始终作用于同一个实际函数 Fe,s(I)。精度变大并不改变被逼近函数的深度下标。以下并不把对一个区间成立的比较自动求和到所有尺度;这里的全区间极限只用于原始有限树。
1. 每个精度已经具有的结论
固定 q∈(2,3)。精度级别 t=1,2,… 依次更新位置
ι(t)=(t−1)mod4∈{0,1,2,3}.
选择绝对常数 A>1、足够大的 k1,并令
kt+1=ktA.
在原求和区间 I 上写 s=s(I)。记第 t 级在位置 ι(t) 的逼近为
gt,I=PItFι(t),s.
不同精度的投影不必彼此嵌套;每个精度内部的空间深度投影族才具有前篇所证的继承性。若 t≤0,相应的 gt,I 约定为零。
前篇的构造给出
∥gt,I∥2,I≤1,dι(t),I(fι(t),s−(gt,I)phys)≤e−kt,
以及
∥(gt,I)phys∥∞,I≤Mt,Mt≤exp(Cqktddet).(Fit-bounds)
必要时把检测估计中的绝对指数 ddet 增大一次,使它同时控制次数上限与投影的物理上界。这个指数不依赖 q,而 Cq 可以依赖 q。
把不同精度、不同重启构造中的私有坐标作互不相交的标记,嵌入同一个物理标量坐标加私有坐标直和的 Hilbert 空间。每个旧投影在其余私有坐标上为零。这样 gt−gt−4 及下面的附加空间范数都有共同空间意义;形式仍只取物理分量。
以下在一个固定 I 中省略下标 I。所有形式都只使用物理分量。归一化实际函数满足 ∥fe,s∥∞,I≤1。
2. 两个指定位置的形式范数
对不同位置 i,j 的函数 ai,aj,定义
∥[ai,aj]∥H,I=∥zl∥2,I,∥zm∥2,I≤1sup∣HI(z0,z1,z2,z3)∣,
其中位置 i,j 固定为 ai,aj,{l,m} 为其余两个位置。比较两个这样的双线性测试时,范数取相应形式之差的上确界。
前篇的 dj,I 允许任意一个对手位置具有 Lq 归一化。因此若位置 j 的残差 r 满足 dj,I(r)≤ε,则任意另一位置的 h∈Lq(I) 都满足
∥[r,h]∥H,I≤ε∥h∥q,I.(Residual-test)
零函数情形直接成立;其余情形将 h 除以其归一化范数即可。
3. 两个不同精度的投影怎样换回实际函数
对应 D 式 (5.9),p.30。
设 v>w 且 ι(v)=ι(w),分别记这两个位置为 j,i。由双线性恒等式,
[gv,gw]−[fj,s,fi,s]=[gv−fj,s,gw]+[fj,s,gw−fi,s].
在第一项中使用晚精度残差,并把 gw 放在 Lq 对手位置;在第二项中使用早精度残差,并把有界实际函数放在 Lq 位置。式 (Residual-test) 给出
[gv,gw]−[fj,s,fi,s]H,I≤e−kvMw+e−kw≤Cq(e−kvexp(Cqkwddet)+e−kw).(Pair-replacement)
若其中一个函数已经是实际函数,只需一次残差测试,例如
∥[gv−fj,s,fi,s]∥H,I≤e−kv.
这里实际函数的 Lq 范数至多一,不产生新的复杂度因子。
取 A>ddet,并预留固定余量。由于 v≥w+1,
kv≥kwA,−kv+Cqkwddet≤−21kv
对足够大的起始 k1 成立,且以后所有 v>w 同时成立。于是式 (Pair-replacement) 的右侧至多
Cqexp(−cqkw)
,常数可以统一选取。这个比较解释了精度递推为何需要幂次增长,而不只要求 kt→∞。
4. 有限终端级别的精确展开
固定终端级别 n。对位置 e,令 ge,n 为不晚于 n 的最后一次该位置逼近;如果还未更新则为零。定义终端残差
re,n=Fe,s−ge,n.
同一位置的增量为 gt−gt−4,所以精确地有
Fe,s=1≤t≤nι(t)=e∑(gt−gt−4)+re,n.(Finite-telescope)
将四个式子代入四线性形式。对一个包含至少两个增量的展开项,记被选出的最早两个级别为 u<t。它们必在不同位置,因为每个位置的线性展开只选择其中一项。令
i=ι(u),j=ι(t),{a,b}={0,1,2,3}∖{i,j}.
定义四元组
Xiu,t=gu−gu−4,Xju,t=gt−gt−4,
Xau,t=Fa,s−ga,prev,Xbu,t=Fb,s−gb,prev,(Two-earliest)
其中 ge,prev 是位置 e 严格早于 t 的最后一次逼近,尚无则为零。
在其余位置 a,b,允许出现的增量只能晚于 t,也可以是终端残差。把这些可能性求和,式 (Finite-telescope) 正好给出 Fe,s−ge,prev。因此,每个原展开项按最早两级恰好归入一组,得到
HI(F0,s,F1,s,F2,s,F3,s)=1≤u<t≤nι(u)=ι(t)∑HI(Xu,t)+Remn,I.(Expansion)
没有增量的项为四个 re,n 的乘积形式。恰有一个位置出现增量时,将该位置所有级别的增量求和,得到 ge,n。所以余项准确地是五项之和:
Remn,I=HI(r0,n,r1,n,r2,n,r3,n)+e=0∑3HI(r0,n,…,ge,n,…,r3,n).
这里省略号只表示第 e 个位置换成 ge,n,其余三个位置保持各自终端残差。
5. 终端余项为什么趋于零
当终端 n≥4 时,四个位置都已有更新。取两个不同位置的终端残差,其最后更新级别 v,w 都在 {n−3,n−2,n−1,n} 中,且不同;设 v>w。晚残差通过精度 e−kv 的检测,早残差的物理 Lq 范数不超过 1+Mw。于是
∥[rι(v),n,rι(w),n]∥H,I≤e−kv(1+Mw)⟶0.
极限由上一节同样的 A>ddet 选择保证。
每个终端残差的归一化 L2 范数至多二,每个投影逼近的归一化 L2 范数至多一。余项的五项各含至少三个残差,因而总能选两个残差使用刚才的形式范数估计,再以另外两个函数的 L2 范数控制。故
Remn,I⟶0
对每个固定原区间成立。
若 U 是固定有限的原求和区间族,则
I∈U∑∣I∣∣Remn,I∣⟶0.
这里只交换极限与一个有限和。此时没有据此声称对任意无限区间族都可直接求和;后续还需独立控制各 (u,t) 项。
6. 晚精度的每对函数为何都具有小形式范数
D Lemma 5.2 还断言:可取绝对常数
0<θ0<A−7<1,
使对 t≥12,四元组 Xu,t 中任意两个不同位置的函数都满足
∥[Xlu,t,Xmu,t]∥H,I≤exp(−2ktθ0).(Small-pairs)
6.1 哪些因子在替换后抵消
把一个因子中出现的每个投影逼近换成对应的实际函数,原来已是实际函数的部分不变。
对于位置 j,
Xj=gt−gt−4
替换后为 Fj,s−Fj,s=0。由于 t≥12,这两次更新都存在。
对于位置 a,b,其严格早于 t 的最后更新在 t−1,t−2,t−3 中,所以
Xe=Fe,s−ge,prev
替换后也为零。只有早位置 i 在 u≤4 时可能没有前一次更新,因而替换后未必抵消。
任意两个位置总含至少一个属于 {j,a,b}。因此,只要逐项替换的误差受到控制,二者全部替换后的实际函数项必为零。
6.2 很早的 i 位置不能粗略使用最早精度误差
先设 u<t−4,并考虑一对包含 Xi=gu−gu−4 的函数。另一函数记为 V,其已有投影级别均至少为 t−4。这些级别都严格晚于 u 及存在的 u−4。
因为 V 替换成实际函数后系数之和为零,可以把它准确改写为
V=v∑εv(gv−Fι(v),s),
其中至多两个项、∣εv∣=1,并且这些 v 全属于 V 的同一个位置。于是
[Xi,V]=w∈{u,u−4}w≥1∑v∑ϵwεv[gw,gv−Fι(v),s].
每项直接用晚残差进行检测,不产生早残差的 e−kw 误差:
∥[gw,gv−Fι(v),s]∥H,I≤e−kvMw≤Cqe−cqkv.
所有 v≥t−4,所以总误差至多 Cqe−cqkt−4。
这一步避免了将一个非常早的 u 的固定小误差错误地当作随 t 衰减的量。
6.3 其余成对比较只含最近七级
仍在 u<t−4 时,若一对不含位置 i,它涉及的已有投影级别均至少为 t−4。展开为有限个“投影或实际函数”的二元组合,对投影—投影使用式 (Pair-replacement),对投影—实际函数使用一次残差测试。全部实际函数项因第 6.1 节的抵消而消失,剩下的误差至多 Cqe−cqkt−4。
若 u≥t−4,由于 ι(u)=ι(t),不能有 u=t−4,所以 u≥t−3。早位置的前次更新级别 u−4 至少为 t−7,其余位置已有投影的级别不更早。相同的成对替换因此给出
∥[Xl,Xm]∥H,I≤Cqexp(−cqkt−7)
对全部不同位置成立。每个因子至多两项,误差项数有一个固定的绝对上界。
6.4 将级别差写成当前精度的固定幂
递推式意味着
kt−7=ktA−7.
选一个绝对的 0<θ0<A−7。再把起始 k1 取得足够大,使
Cqexp(−cqktA−7)≤exp(−2ktθ0)
对以后所有 t≥12 成立,即得到式 (Small-pairs)。
θ0 不依赖区间、深度或输入分块数。它也可在 q 前由绝对次数和 A 确定;q 的依赖由起始精度吸收。
7. 各函数的大小与共同局部构造
每个 Xi,Xj 是两个范数至多一的投影之差,每个 Xa,Xb 是实际函数与一个投影之差。因此
∥Xeu,t∥2,I≤2(e=0,1,2,3).
这是附加 Hilbert 空间的范数界,当然也控制形式中所用物理分量。
一个固定四元组只涉及有界多个精度的空间分段构造。将这些分段取共同细分时,非空的二进根区间交等于其中最深的区间;所需共同根来自这些有限个既有根族。因而其总根长仍由原根长乘一个固定常数控制。这个处理不重置原始空间深度 s(I)。
原文范围与后续问题
- D p.29,循环更新规则:ι(t)、精度递推、各次投影始终作用于固定的 Fe,s(I)。
- D p.30,式 (5.7)–(5.8):最早两个增量的有限展开与五项终端余项。
- D p.30,式 (5.9) 与 Lemma 5.2 后半:两精度比较、极早位置的单独处理、t−7 与 θ0<A−7。
下一步是 D §§5.3–5.6:将两个残差位置换成延迟投影,保留这里得到的成对小相关。需要先由两个固定低位置结构函数确定有限预测函数族,再选择高位置投影。延迟的是投影深度,原操作数必须仍为同一个 Fe,s(I)。这部分尚未由本篇的逐区间比较自动推出。