|
MvR的主要对角线元素也是正的,叫做dj,j(MvR)=√τhvβ+W(vj)> 0,j=1。。。,N.剩下的证据可以根据Itkin(2014b)中的定义和引理A.2再次证明。因此,基于这两个步骤,必须选择系数γ,以提供x,v[W(v)+ρs,vW(s)]√τhS(γ+1)>1。由于公式(34)中的两个步骤都在谱范数中收敛,并且是无条件稳定的,因此整个格式的无条件稳定性和收敛性如下。也就是说,整个方案保持了解的非负性。公式(34)中的空间近似值lhs用hS中的二阶近似值表示,而rhs部分中的第一行用一阶导数的一阶近似值表示。当OS=ABS+O(hS)时,在等式(34)的rhs的第一行,我们有一个产品√τOS,被忽略项的顺序是O(√τhS)。严格地说,整个计划。(34)提供了这种近似顺序。B命题2.3的证明与命题2.1相比,该方案有唯一的修改。也就是说,不是等式(34)中的第一步QIv+√τW(v)AB2,v五、*= α+V(τ)- Vk,α+=(pq+1)I- Q√τρs,vW(s)ABS+P√τW(v)AFv。我们现在使用QIv+√τW(v)AB2,v五、*= α+V(τ)- Vk,(39)α+=(pq+1)I- Q√τρs,vW(s)AB2,s+P√τW(v)AF2,v.下面我们想证明v=α+v(τ)- Vkis是一个正向量。假设我们以V=V(τ)开始迭代。然后V=VV(τ)+VV(τ),(40)V=pq- Q√τρs,vW(s)AB2,s,V=pq+P√τW(v)AF2,v.向量U=VV(τ)可以在vi上的循环中计算,i=1。。。,N.换句话说,如果V(τ)是一个N×N矩阵,其中行表示S坐标,列表示V坐标,则ui的第j列是矩阵和V(τ)的第j列的乘积。
|