Try to project gas particles (SPH) to a mesh with octree structure

[zsw6666]

Hi, everyone! I am trying to project gas particles (SPH) to a mesh in the octree structure.
The particle data contains several np.array:
particle_position: (n, 3) array where n=5303345 is the number of particles
particle_density: (n, 1) array
I use the following code to do this job:
particle_ds=yt.load_particles(particle_data)
gas_ds.add_sph_fields() # to generate the fields of (‘io’, ‘smoothing_length’) and (‘io’, ‘density’). add_sph_fields only iterates 5303342 steps

right_edge=[xxx,xxx]
left_edge=[xxx,xxx]
particle_oct=particle_ds.octree(right_edge=right_edge,left_edge=left_edge,n_ref=10)
########################################
Then, I would like to access some fields of particle_oct:
cell_density=gas_oct[‘io’, ‘density’]
########################################
This raise the following error:
IndexError: boolean index did not match indexed array along dimension 0; dimension is 5303342 but corresponding boolean dimension is 5303345
########################################
How do I fix this error and get the correct gas_oct[‘io’,‘density’]?

[welcome]

Hi, and welcome to yt! Thanks for opening your first issue. We have an issue template that helps us to gather relevant information to help diagnosing and fixing the issue.

[zsw6666]

I just figured out that this error will be solved if I change the left_edge` and right_edge to a broader value.