逐步主鏈
S01 · definition
操作前
\text{declared BVP on } \Omega=(0,2)\times(0,1)操作後
(I):\ \Delta u-4\equiv 0\ \text{in}\ \Omega;\quad (B):\ u|_{\Gamma_L}=y^2,\ u|_{\Gamma_R}=y^2+2y+4,\ u|_{\Gamma_B}=x^2,\ u|_{\Gamma_T}=x^2+x+1條件:無
規則:mko_bvp · 依賴:G1, G2, mko_bvp
葉節點依據:mko_bvp + def_unfold
依 MKO,「候選滿足宣告的 BVP」由兩組不同的義務組成:(I) 內部 Poisson 方程 Δu=4 在 Ω 的每一點成立;(B) u 在四條邊的每一點等於 Dirichlet 資料。先把兩組義務分開寫下,之後各自取得證據,最後才合併。
S02 · implies
操作前
u=x^2+xy+y^2
操作後
u_x=2x+y,\ u_y=x+2y
條件:無
規則:poly_partial · 依賴:G3, poly_partial
葉節點依據:poly_partial
逐項求偏導:∂x(x²)=2x、∂x(xy)=y、∂x(y²)=0,得 u_x=2x+y;∂y(x²)=0、∂y(xy)=x、∂y(y²)=2y,得 u_y=x+2y。
S03 · implies
操作前
u_x=2x+y,\ u_y=x+2y
操作後
u_{xx}=2,\ u_{yy}=2條件:無
規則:poly_partial · 依賴:S02, poly_partial
葉節點依據:poly_partial
u_xx=∂x(2x+y)=2,u_yy=∂y(x+2y)=2。Laplacian 只需要這兩個純二階偏導。
S04 · implies
操作前
\Delta u-4=u_{xx}+u_{yy}-4操作後
\Delta u-4=2+2-4=0\ \text{(zero polynomial)};\ (I)\ \text{holds}條件:無
規則:poly_coeff_zero · 依賴:S03, S01, poly_coeff_zero
葉節點依據:poly_coeff_zero + numeric_evaluation + forall_intro
Δu−4=2+2−4=0。內部殘差是零多項式(常數 0),所以在 Ω 的每一點、事實上在 ℝ² 的每一點都等於 0。這是符號恆等,不是在內部取幾個點核對。記錄:u 滿足義務 (I)。
S05 · equivalent
操作前
u=g\ \text{on}\ \partial\Omega操作後
u=g_L\ \text{on}\ \Gamma_L\ \land\ u=g_R\ \text{on}\ \Gamma_R\ \land\ u=g_B\ \text{on}\ \Gamma_B\ \land\ u=g_T\ \text{on}\ \Gamma_T條件:無
規則:set_operations_def · 依賴:G1, S01, set_operations_def
葉節點依據:set_operations_def + and_intro_elim
∂Ω=Γ_L∪Γ_R∪Γ_B∪Γ_T(G1)。一個點在聯集中,當且僅當它在某一條邊上;所以「在 ∂Ω 的每一點 u 等於資料」等價於「在四條邊的每一條的每一點都相等」。四個合取項缺一不可。
S06 · implies
操作前
u|_{\Gamma_L}=u(0,y)操作後
u(0,y)=y^2=g_L(y)\ \ (\text{coefficients } (1,0,0)=(1,0,0))條件:0\le y\le 1
規則:edge_restrict · 依賴:S05, G2, G3, edge_restrict, poly_coeff_zero
葉節點依據:edge_restrict + poly_coeff_zero + numeric_evaluation
左邊 x=0:u(0,y)=0²+0·y+y²=y²。按 (y², y, 1) 的係數,限制結果為 (1,0,0),資料 g_L 也是 (1,0,0),逐係數相同,所以在整條左邊(每個 y∈[0,1])u=g_L。這是全邊恆等,不是抽查幾個 y。
S07 · implies
操作前
u|_{\Gamma_R}=u(2,y)操作後
u(2,y)=y^2+2y+4=g_R(y)\ \ (\text{coefficients } (1,2,4)=(1,2,4))條件:0\le y\le 1
規則:edge_restrict · 依賴:S05, G2, G3, edge_restrict, poly_coeff_zero
葉節點依據:edge_restrict + poly_coeff_zero + numeric_evaluation
右邊 x=2:u(2,y)=2²+2y+y²=y²+2y+4。係數 (1,2,4) 與 g_R 的 (1,2,4) 相同,所以在整條右邊 u=g_R。
S08 · implies
操作前
u|_{\Gamma_B}=u(x,0)操作後
u(x,0)=x^2=g_B(x)\ \ (\text{coefficients } (1,0,0)=(1,0,0))條件:0\le x\le 2
規則:edge_restrict · 依賴:S05, G2, G3, edge_restrict, poly_coeff_zero
葉節點依據:edge_restrict + poly_coeff_zero + numeric_evaluation
下邊 y=0:u(x,0)=x²+x·0+0²=x²。按 (x², x, 1) 的係數 (1,0,0) 與 g_B 相同,所以在整條下邊 u=g_B。
S09 · implies
操作前
u|_{\Gamma_T}=u(x,1)操作後
u(x,1)=x^2+x+1=g_T(x)\ \ (\text{coefficients } (1,1,1)=(1,1,1))條件:0\le x\le 2
規則:edge_restrict · 依賴:S05, G2, G3, edge_restrict, poly_coeff_zero
葉節點依據:edge_restrict + poly_coeff_zero + numeric_evaluation
上邊 y=1:u(x,1)=x²+x·1+1²=x²+x+1。係數 (1,1,1) 與 g_T 相同,所以在整條上邊 u=g_T。
S10 · implies
操作前
u=g_L,\ u=g_R,\ u=g_B,\ u=g_T\ \text{on their edges}操作後
(B)\ \text{holds for } u條件:無
規則:and_intro_elim · 依賴:S05, S06, S07, S08, S09
葉節點依據:and_intro_elim + set_operations_def
四條邊的恆等都已證明(S06 至 S09);依合取引入與 S05 的等價,u 在整個 ∂Ω 上等於 Dirichlet 資料。記錄:u 滿足義務 (B)。
S11 · implies
操作前
(I)\ \text{for } u,\ (B)\ \text{for } u操作後
u\ \text{satisfies the declared BVP}條件:無
規則:and_intro_elim · 依賴:S04, S10, S01
葉節點依據:and_intro_elim + mko_bvp
兩組義務都已分開記錄(S04 與 S10),依 MKO 此時才可說 u 滿足宣告的 BVP。附帶觀察:每個角點的資料來自兩條邊,而且一致,例如 (2,1) 處 g_R(1)=1+2+4=7=g_T(2)=4+2+1,與它們都是同一個函數的限制相符。
S12 · implies
操作前
v=2x^2+xy-2x+y
操作後
\Delta v-4=4+0-4=0;\ (I)\ \text{holds for } v條件:無
規則:poly_partial · 依賴:G3, poly_partial
葉節點依據:poly_partial + poly_coeff_zero
依 S12.1、S12.2,v_xx=4、v_yy=0,所以 Δv−4=0 是零多項式,v 也滿足內部義務 (I)。注意 v=u+h,其中 h=x²−y²−2x+y 是調和的(h_xx+h_yy=2−2=0),所以內部義務不能區分 u 與 v。
S13 · implies
操作前
v(0,0),\ v(2,0),\ v(0,1),\ v(2,1)
操作後
0,\ 4,\ 1,\ 7\ =\ g\ \text{at the four corners}條件:無
規則:numeric_evaluation · 依賴:G2, G3
葉節點依據:numeric_evaluation
四個角點:v(0,0)=0,資料 g_B(0)=0;v(2,0)=8+0−4+0=4,資料 g_B(2)=4;v(0,1)=0+0−0+1=1,資料 g_L(1)=1;v(2,1)=8+2−4+1=7,資料 g_R(1)=7。四個角點全部吻合,因為 h 在四個角點都為 0。只抽查角點會誤以為 v 符合邊界資料。
S14 · implies
操作前
v|_{\Gamma_B}=v(x,0)操作後
v(x,0)=2x^2-2x\neq x^2;\ v(1,0)=0\neq 1=g_B(1)
條件:0\le x\le 2
規則:edge_restrict · 依賴:S05, G2, G3, edge_restrict
葉節點依據:edge_restrict + numeric_evaluation
下邊 y=0:v(x,0)=2x²+0−2x+0=2x²−2x。係數 (2,−2,0) 與 g_B 的 (1,0,0) 不同。取下邊上的點 x=1∈[0,2]:v(1,0)=2−2=0,而 g_B(1)=1,所以 v 在下邊不等於資料。兩者之差 v(x,0)−x²=x(x−2) 只在 x=0 與 x=2(兩個角點)為零。
S15 · implies
操作前
v|_{\Gamma_L}=v(0,y)操作後
v(0,y)=y\neq y^2;\ v(0,\tfrac12)=\tfrac12\neq\tfrac14=g_L(\tfrac12)
條件:0\le y\le 1
規則:edge_restrict · 依賴:S05, G2, G3, edge_restrict
葉節點依據:edge_restrict + numeric_evaluation
左邊 x=0:v(0,y)=0+0−0+y=y,係數 (0,1,0) 與 g_L 的 (1,0,0) 不同;在 y=1/2∈[0,1],v=1/2 而 g_L=1/4。左邊同樣失敗。(右邊 v(2,y)=3y+4 與上邊 v(x,1)=2x²−x+1 也與資料不同,本稿不再逐一展開。)
S16 · implies
操作前
(I)\ \text{for } v,\ v(1,0)\neq g_B(1)操作後
(B)\ \text{fails for } v;\ v\ \text{does not satisfy the declared BVP}條件:無
規則:quantifier_negation · 依賴:S05, S12, S13, S14, quantifier_negation
葉節點依據:quantifier_negation + and_intro_elim + exists_intro
(B) 是「邊界上每一點都相等」的全稱命題;S14 給出一個不相等的邊界點,依量詞否定 (B) 對 v 不成立。BVP 要求 (I) 與 (B) 同時成立,所以 v 不滿足宣告的 BVP,儘管 (I) 成立(S12)而且四個角點都吻合(S13)。