1 | This factorized Operator can now be solved in three steps: |
2 | Each step involves solution of a 5 by 5 block tridiagonal matrix |
3 | The matrix is tridiagonal -- not pentadiagonal -- as numerical dissipation with Ñ4U is only added to explicit RHS term in this approximation |