Raycast.jl provides a simple but fast 3D raycasting function. It is non-allocating and only uses integer-arithmetic.
using Raycast
p0 = (0,0,0) # any tuple of integers will work
p1 = (10, 5, -3) # any tuple of integers will work
prob = BresenhamProblem(p0, p1)
for p in prob
@show p
end] add https://github.com/dev10110/Raycast.jl
See the docs for details.