#2015: thorn Vectors does not provide a sum() reduction --------------------------+------------------------------------------------- Reporter: rhaas | Owner: Type: enhancement | Status: new Priority: minor | Milestone: Component: Cactus | Version: development version Resolution: | Keywords: Vectors --------------------------+-------------------------------------------------
Comment (by rhaas):
Thank you. I was looking at CarpetInterp2's lagrange interpolation. See https://bitbucket.org/eschnett/carpet/branch/rhaas/fasterp though this has not yet been speed tested compared to the previous implementation.
The sum() at the end is done via a loop over the elt function. I did a quick check to see what an omp simd reduction(+: var) would do and found that at least on my laptop it does a similar loop and no horizontal sum at all.