TheVeggydude
commented
4 years ago Cannot be done through a single transformation matrix. The values for delta_S and delta_E depend on a multiplication between two elements of the original vector.
- Either the transformation matrix needs to be updated every iteration to include the param value in the correct cell,
- or the original vector needs to be multiplied by an artificial vector containing [I, I, 0, 0]^T
Moreover, the matrix is also sparse. This means that most of the cells are zero. A lot of spaces are, therefore, not used. So it may not be the most efficient solution.
The update loop is best done in a naive manner, as the matrix methods do not add any speed up, readability or ease of use.
The math for updating a single node could potentially be done by matrix/vector computations. Verify this.