Jesse Chan

A related discussion on "paired" RK methods https://github.com/trixi-framework/Trixi.jl/issues/21 and associated implementation issues. These are somewhat simpler than multi-rate schemes but aim towards a similar goal.

Thanks for the explanation. I was using `BroadcastArray`, but that seemed to result in a lot of time being spent in `lazybroadasting.jl` (see https://github.com/trixi-framework/Trixi.jl/pull/757#issuecomment-892157256). I probably need to look a...

Sure - here's an example. I'm using LazyArrays for linear combinations of matrices. ```julia using LazyArrays using StructArrays n = 100 u = randn(n) du = similar(u) A, B, C,...

Gotcha. I don't think I can write `foo_lazy!` as a matrix multiplication easily, but it's helpful to know `getindex` is the bottleneck.

Was `reset_workers!` removed in the recent 0.5.1 update? I still need to reset the state manually from time to time but can't seem to access it anymore.

Thanks!

Just so its visible for others, this is what I'm using with PolyesterWeave 0.1 and ThreadingUtilities 0.4.6: ```julia PolyesterWeave.reset_workers!() for i = 1:Threads.nthreads()-1; ThreadingUtilities.initialize_task(i); end ```

Turns out this is just because the constructor used in the example calls `recursive_bottom_eltype` to determine `T`. Using the default constructor instead and specifying `T` manually fixes things, and broadcasting...

Oops, I copied and pasted the wrong example code. StructArrays were just in the MWE for similarity to the full application. Switching to the other constructor fixed an unrelated issue,...

Sorry - to clarify, the problem remains. I solved an unrelated issue recently, and it'd been too long since I looked at this and thought I'd solved it too.