Oliver Schulz
Oliver Schulz
Haven't forgotten about this.
Still not forgotten, just overloaded ...
Sorry for the long silence on this - I've been thinking a bit: Now that view-like wrapper objects don't have to be stack-allocated anymore, we can just add a method...
CC @devmotion
Just ran into this in another context: #534 isn't quite the same, because it's about values that are vectors of real numbers, but the cause is the (recursive) behavior of...
Sounds good. Regarding on where do draw the line - BLAS would only accelerate `Real`, right (resp. do the wrong thing for `Complex` for this use case)? So we'd do...
I've been thinking - maybe we should drop use of `dot` altogether, though, even if it means a performance loss, since we'd gain precision. AFAIK, Julia's `sum()` uses pairwise summation,...
> `mapreduce(*, +, A, values(w))` doesn't use pairwise summation either [...] Yes, it does (AFAIK the default `Base.sum` is actually defined via `mapreduce, which operates pairwise). Here's a more extreme...
> Actually mapreduce(*, +, A, w) allocates a copy to store the result, so that's not an option. Oops, indeed - I naively thought it would be allocation-free. In that...
`sum(Base.Broadcast.Broadcasted(*, (A, values(w))))` seems to be allocation-free, but it's fairly slow.