place_particles_lattice Interface

interface


Calls

interface~~place_particles_lattice~~CallsGraph interface~place_particles_lattice place_particles_lattice proc~place_particles_lattice place_particles_lattice interface~place_particles_lattice->proc~place_particles_lattice none~validate_position_final~3 validate_position_final proc~place_particles_lattice->none~validate_position_final~3 validate_position validate_position none~validate_position_final~3->validate_position

public module subroutine place_particles_lattice(this, central_density, xmin, xmax, ymin, ymax, zmin, zmax, npart_des, npart_out, stretch, thres, pos, pvol, nu, h, get_density, get_id, validate_position)

Arguments

Type IntentOptional Attributes Name
class(particles), intent(inout) :: this

particles object which this PROCEDURE is a member of
double precision, intent(in) :: central_density

Maximum baryon mass density of the system
double precision, intent(in) :: xmin

Left boundary of the lattice
double precision, intent(in) :: xmax

Right boundary of the lattice
double precision, intent(in) :: ymin

Left boundary of the lattice
double precision, intent(in) :: ymax

Right boundary of the lattice
double precision, intent(in) :: zmin

Left boundary of the lattice
double precision, intent(in) :: zmax

Right boundary of the lattice
integer, intent(in) :: npart_des

Desired particle number
integer, intent(out) :: npart_out

Real, output particle number
double precision, intent(in) :: stretch

Stretching factor fo the lattice. xmin to zmax are multiplied by it
double precision, intent(in) :: thres

(~rho_max)/thres is the minimum mass density considered when placing particles. Used only when redistribute_nu is .FALSE. . When redistribute_nu is .TRUE. thres= 100*nu_ratio
double precision, intent(inout), DIMENSION(:,:), ALLOCATABLE :: pos
double precision, intent(inout), DIMENSION(:), ALLOCATABLE :: pvol

Array storing the particle positions Array storing the particle volumes
double precision, intent(inout), DIMENSION(:), ALLOCATABLE :: nu

Array storing the particle baryon masses
double precision, intent(inout), DIMENSION(:), ALLOCATABLE :: h

Array storing the initial guess for the particle smoothing lengths
function get_density(x, y, z) result(density)

Returns the baryon mass density at the desired point

Arguments
Type IntentOptional Attributes Name
double precision, intent(in) :: x

coordinate of the desired point
double precision, intent(in) :: y

coordinate of the desired point
double precision, intent(in) :: z

coordinate of the desired point
Return Value double precision

Baryon mass density at

subroutine get_id(x, y, z, sqdetg, baryon_density, gamma_euler)

Returns the baryon mass density at the desired point

Arguments
Type IntentOptional Attributes Name
double precision, intent(in) :: x

coordinate of the desired point
double precision, intent(in) :: y

coordinate of the desired point
double precision, intent(in) :: z

coordinate of the desired point
double precision, intent(out) :: sqdetg
double precision, intent(out) :: baryon_density
double precision, intent(out) :: gamma_euler
procedure(validate_position_int), optional :: validate_position

Returns 1 if the position is not valid, 0 otherwise

Description

Places particles on a lattice containing a physical object