|
We can write the right hand side of equation (3.26) as the matrix vector product
|
The vector Uvec has Nvec = 5 Nx Nh Nz components and RHS is a Nvec by Nvec matrix. The forcing term Hvec is also a vector with Nvec components
-
The components of vectors and rows and columns of matrices are labelled by i,j,k and CFD index index running from 1 to 5
|
|
Formally we can write full equations:
|
If we had used the explicit Euler method, then equations would have read the simple:
-
which only needs sparse matrix vector multiplication as its computational kernel.
|
