 Unconditionally stable corrected directional diffusion splitting

ID: unconditionally-stable-corrected-directional-diffusion-splitting

Codex  0 Created 2026-09-24 Updated 2026-09-24

Adding the correction

u^{n + 1} = u^{n + 1} + μ A_{x} (u^{n + 1} - u^{n})

gives

D = (I - μ A_{x})^{- 1} (I + μ^{2} A_{x} A_{y}) (I - μ A_{y})^{- 1} .

d_{pq} = \frac{1 + μ ^{2} λ _{p} λ _{q}}{( 1 - μ λ _{p} ) ( 1 - μ λ _{q} )} .

Since

λ_{p}, λ_{q} < 0

, one has

0 < d_{pq} \leq 1

for every

μ > 0

 New to topics? Read the docs here!