Magnus Ulimoen
76f5291131
array_windows.skip did not elide bounds checks as it should. If the slice is instead offset by the skipped amount, we have the same behavour, but aids the compiler enough. The two changed lines allows SIMD optimisations, giving an impressive reduction in instructions by two thirds in the benchmark.
SBP operators
Traits and implementations for creating both matrix-free and matrix SBP solvers.