Task: Construct an explicit CSR matrix inside the Residual class and benchmark it against the current matrix-free approach.
This is primarily an exploratory issue. The matrix-based approach may be advantageous in contexts where the operator is applied many times (amortizing assembly cost) or where hardware prefers structured memory access. The matrix-free approach is expected to win on memory footprint, but the crossover point is worth understanding empirically.
Results should be documented and used to inform future architecture decisions.
Task: Construct an explicit CSR matrix inside the
Residualclass and benchmark it against the current matrix-free approach.This is primarily an exploratory issue. The matrix-based approach may be advantageous in contexts where the operator is applied many times (amortizing assembly cost) or where hardware prefers structured memory access. The matrix-free approach is expected to win on memory footprint, but the crossover point is worth understanding empirically.
Results should be documented and used to inform future architecture decisions.