Hello! Suppose I have the parameters: using StaticArrays using LinearAlgebra h = 0.05 eta2 = (0.01)h * (0.01)h points = rand(SVector{3,Float64},100000) v = points * 0.05 And I wish to evaluate this expression: maximum(abs.(h* dot.(v,points) ./ (dot.(points,points) .+ eta2))) W…

You can get another bit of speed like: julia> function maxvp_2(v,points,h,eta2) maxval = -Inf for i in eachindex(v,points) tmp = abs(h * dot(v[i],points[i]) / (dot(points[i],points[i]) + eta2)) tmp > maxval && (maxval = tmp) end …

How do I make this allocation free?

New to Julia

uniment January 21, 2023, 11:40pm 18

Seems related:

TLDR: performance of maximum has a lot of room for improvement, and there is an ongoing effort to fix it.

1 Like