Guide

julia> using PhysicalParticles, PhysicalTrees, UnitfulAstro
julia> pos = [PVector(1.0, 1.0, 1.0, u"kpc"), PVector(-1.0, -1.0, -1.0, u"kpc"), PVector(1.0, 0.0, -1.0, u"kpc"), PVector(-1.0, 0.0, 1.0, u"kpc"), PVector(0.0, 0.0, -1.0, u"kpc"), PVector(-1.0, 0.0, 0.0, u"kpc")];
julia> tree = octree(pos) Tree defined on worker 1: Distributed on workers: [1] units: Unitful.FreeUnits{N, D, nothing} where {N, D}[kpc, Gyr, A, K, cd, M⊙, mol] Number of topnodes: 953 Number of topleaves: 834 Domain factor: 1.0475284715284717e6 kpc^-1 Domain start list: [1] Domain end list: [834] Load list: [6] Work list: [0.0] Number of tree nodes: [953] Extent: xMin = -1.0 kpc, xMax = 1.0 kpc, yMin = -1.0 kpc, yMax = 1.0 kpc, zMin = -1.0 kpc, zMax = 1.0 kpc, SideLength = 2.002 kpc, Center = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc) -------------------------- Data info -------------------------- total: 6 cuts: [6] last communicate: Dict(1 => 0)
julia> particles = StructArray(Star(uAstro) for i in 1:6);
julia> assign_particles(particles, :Pos, pos)
julia> tree = octree(particles) Tree defined on worker 1: Distributed on workers: [1] units: Unitful.FreeUnits{N, D, nothing} where {N, D}[kpc, Gyr, A, K, cd, M⊙, mol] Number of topnodes: 953 Number of topleaves: 834 Domain factor: 1.0475284715284717e6 kpc^-1 Domain start list: [1] Domain end list: [834] Load list: [6] Work list: [0.0] Number of tree nodes: [953] Extent: xMin = -1.0 kpc, xMax = 1.0 kpc, yMin = -1.0 kpc, yMax = 1.0 kpc, zMin = -1.0 kpc, zMax = 1.0 kpc, SideLength = 2.002 kpc, Center = PVector(0.0 kpc, 0.0 kpc, 0.0 kpc) -------------------------- Data info -------------------------- total: 6 cuts: [6] last communicate: Dict(1 => 0)