Each coordinate of the transformed vectors evolves independently because and are diagonal matrices. Reorder the unknowns by collecting the th components across all grid columns:
Then system (3) splits into the uncoupled systems
where
Thus the original problem reduces to one discrete sine transform followed by independent tridiagonal solves and an inverse sine transform.

Articles by others on the same topic (0)

There are currently no matching articles.