The first stage is
Because , this consists of independent tridiagonal systems of dimension . The Thomas algorithm solves each in operations, for a total of .
The second stage,
is multiplication by a matrix with only nonzero entries, so it also costs . Hence one complete step uses at most
arithmetic operations.

Articles by others on the same topic (0)

There are currently no matching articles.