bns_fuka Module

public module subroutine allocate_bnsfuka_memory(this, d) Implementation →

Allocates allocatable arrays member of a bnsfuka object

Arguments

public module subroutine construct_binary(this, fukafile) Implementation →

Interface of the subroutine that constructs the $\texttt{FUKA}$ $\texttt{Bin_NS}$ object

Arguments

private function construct_bns_fuka(fuka_file) result(optr) bind(C, NAME= "construct_bns_fuka")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that constructs the $\texttt{FUKA}$ $\texttt{Bin_NS}$ object

FT

Arguments

Return Value type(C_PTR)

C pointer pointing to the constructed $\texttt{FUKA}$ BNS_export object

public module subroutine construct_bnsfuka(derived_type, filename, eos_filenames) Implementation →

Prints a summary of the physical properties the system to the standard output and, optionally, to a formatted file whose name is given as the optional argument filename

Arguments

public module subroutine deallocate_bnsfuka_memory(this) Implementation →

Deallocates allocatable arrays member of a bnsfuka object

Arguments

public module subroutine destruct_binary(this) Implementation →

Destructs a $\texttt{FUKA}$ $\texttt{Bin_NS}$ object

Arguments

private subroutine destruct_bns_fuka(optr) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that destructs the $\texttt{FUKA}$ $\texttt{Bin_NS}$ object

FT

Arguments

public module subroutine destruct_bnsfuka(this) Implementation →

Destruct a bnsfuka object

Arguments

public module subroutine finalize(npart, pos, nlrf, u, pr, vel_u, theta, nstar, nu) Implementation →

Post-process the $\mathrm{SPH}$ $\mathrm{ID}$; for example, correct for the residual ADM linear momentum.

Arguments

public module function get_bns_identifier(this) Implementation →

Arguments

Return Value double precision

public module function get_eos1_fukaid(this) Implementation →

Arguments

Return Value integer

public module function get_eos2_fukaid(this) Implementation →

Arguments

Return Value integer

public module subroutine get_eos_parameters(this, i_matter, eos_params) Implementation →

Arguments

public module function get_field_array(this, field) result(field_array) Implementation →

Returns the bnsfuka member arrays named field

Arguments

Return Value double precision, DIMENSION(:), ALLOCATABLE

Desired bnsfuka member array

public module function get_field_value(this, field, n) result(field_value) Implementation →

Returns the component n of the bnsfuka member arrays named field

Arguments

Return Value double precision

Component n of the desired bnsfuka member array

private subroutine get_fuka_id(optr, x, y, z, lapse, shift_x, shift_y, shift_z, psi4, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, mass_density, specific_energy, pressure, v_euler_x, v_euler_y, v_euler_z) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that reads the full $\texttt{FUKA}$ $\mathrm{ID}$ at the specified point. That is, read_fukas the metric fields, the components of the extrinsic curvature [c/km], and the hydro fields.

• shift vector [c]
• baryon mass density [kg m^{-3}]
• energy density [kg c^2 m^{-3}]
• pressure [kg c^2 m^{-3}]
• specific internal energy [c^2]
• fluid 3-velocity with respect to the Eulerian observer [c]

FT

Arguments

private subroutine get_fuka_id_hydro(optr, x, y, z, mass_density, specific_energy, pressure, v_euler_x, v_euler_y, v_euler_z) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that reads the hydro fields from $\texttt{FUKA}$, at the specified point

• baryon mass density [kg m^{-3}]
• energy density [kg c^2 m^{-3}]
• pressure [kg c^2 m^{-3}]
• specific internal energy [c^2]
• fluid 3-velocity with respect to the Eulerian observer [c]

FT

Arguments

private subroutine get_fuka_id_k(optr, x, y, z, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that reads the components of the extrinsic curvature [c/km] from $\texttt{FUKA}$, at the specified point

FT

Arguments

private subroutine get_fuka_id_mass_b(optr, x, y, z, psi4, mass_density, gamma_euler) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that reads the hydro fields and the metric fields from $\texttt{FUKA}$, at the specified point, needed to compute the baryon mass.

• shift vector [c]
• baryon mass density [kg m^{-3}]
• fluid 3-velocity with respect to the Eulerian observer [c]

FT

Arguments

private subroutine get_fuka_id_params(optr, angular_vel, distance, mass1, mass2, massg1, massg2, radius1_min, radius1_max, radius2_min, radius2_max, adm_mass, komar_mass, adm_linear_momentum_x, adm_linear_momentum_y, adm_linear_momentum_z, adm_angular_momentum_z, COMx, COMy, COMz, area_radius1, center1_x, area_radius2, center2_x, ent_center1, rho_center1, energy_density_center1, ent_center2, rho_center2, energy_density_center2, eos_type_1, eos_file_1, eos_type_2, eos_file_2, gamma, kappa, n_poly, gamma0, gamma1, gamma2, gamma3, kappa0, kappa1, kappa2, kappa3, logP1, logRho0, logRho1, logRho2) bind(C, NAME= "0")

FT

Arguments

private subroutine get_fuka_id_particles(optr, x, y, z, lapse, shift_x, shift_y, shift_z, psi4, mass_density, specific_energy, pressure, v_euler_x, v_euler_y, v_euler_z) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that reads the hydro fields and the metric fields * from $\texttt{FUKA}$, at the specified point

• shift vector [c]
• baryon mass density [kg m^{-3}]
• energy density [kg c^2 m^{-3}]
• pressure [kg c^2 m^{-3}]
• specific internal energy [c^2]
• fluid 3-velocity with respect to the Eulerian observer [c]

FT

Arguments

private subroutine get_fuka_id_spacetime(optr, x, y, z, lapse, shift_x, shift_y, shift_z, psi4, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz) bind(C, NAME= "0")

Interface to the $\texttt{FUKA}$ method of class $\texttt{Bin_NS}$ with the same name, that reads the metric fields and the components of the extrinsic curvature [c/km] from $\texttt{FUKA}$, at the specified point

FT

Arguments

private function get_fuka_mass_density(optr, x, y, z) result(res) bind(C, NAME= "get_fuka_mass_density")

FT

Arguments

Return Value real(kind=C_DOUBLE)

Baryon mass density $[\mathrm{kg}\, \mathrm{m}^{-3}]$ at the desired point $(x,y,z)$

private function get_fuka_pressure(optr, x, y, z) result(res) bind(C, NAME= "get_fuka_pressure")

Interface to the $\texttt{FUKA}$ method of class bns_export with the same name, that returns the pressure $[\mathrm{kg}\, c^2 \mathrm{m}^{-3}]$ from $\texttt{LORENE}$, at the specified point

FT 27.05.2022

Arguments

Return Value real(kind=C_DOUBLE)

Pressure $[\mathrm{kg}\,c^2\, \mathrm{m}^{-3}]$ at the desired point $(x,y,z)$

private function get_fuka_spatial_metric(optr, x, y, z) result(res) bind(C, NAME= "get_fuka_g")

FT

Arguments

Return Value real(kind=C_DOUBLE)

Spatial metric component $g_{xx}=g_{yy}=g_{zz}$ at the point $(x,y,z)$

public recursive module subroutine initialize_id_bnsfuka(this, flag, switch) Implementation →

Initialize the $\texttt{FUKA}$ $\mathrm{BNS}$ $\mathrm{ID}$.

• If flag= [[utility:flag$sph]], set up the lattices around the stars for the $\mathrm{BNS}$ produced with $\texttt{FUKA}$.
• If flag= [[utility:flag$tpo]], allocate memory for the hydro grid functions.
• If flag > 0, assign its value to l_curr.
• If [[utility:flag$tpo]] < flag < 0, assign its value to tpo_curr.

Arguments

public module subroutine interpolate_fuka_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler) Implementation →

Stores the hydro $\mathrm{ID}$ in the arrays needed to compute the baryon mass

Arguments

public module subroutine interpolate_fuka_id_particles(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z) Implementation →

Stores the hydro $\mathrm{ID}$ in the arrays needed to compute the $\mathrm{SPH}$ $\mathrm{ID}$

Arguments

public module function interpolate_fuka_mass_density(this, x, y, z) result(res) Implementation →

Returns the $\texttt{FUKA}$ baryon mass density at a point $(x,y,z)$

Arguments

Return Value double precision

Baryon mass density at $(x,y,z)$

public module function interpolate_fuka_pressure(this, x, y, z) result(res) Implementation →

Returns the $\texttt{FUKA}$ pressure at a point $(x,y,z)$

Arguments

Return Value double precision

Pressure at $(x,y,z)$

public module function interpolate_fuka_spatial_metric(this, x, y, z) result(res) Implementation →

Returns the $\texttt{FUKA}$ conformally flat spatial metric component $g_{xx}=g_{yy}=g_{zz}$ at a point $(x,y,z)$

Arguments

Return Value real(kind=C_DOUBLE)

$g_{xx}=g_{yy}=g_{zz}$ at $(x,y,z)$

public module function is_hydro_positive(this, x, y, z) result(res) Implementation →

Returns 1 if the energy density or the specific energy or the pressure are negative, 0 otherwise

Arguments

Return Value logical

.TRUE. if the energy density or the specific energy or the pressure are negative, .FALSE. otherwise

public module function is_hydro_positive_interpolation(this, x, y, z) result(res) Implementation →

Returns 1 if the energy density or the specific energy or the pressure are negative, 0 otherwise

Arguments

Return Value logical

.TRUE. if the energy density or the specific energy or the pressure are negative, .FALSE. otherwise

private function positive_hydro(optr, x, y, z) result(res) bind(C, NAME= "is_fuka_hydro_positive")

FT

Arguments

Return Value integer(kind=C_INT)

1 if the energy density or the specific energy or the pressure are positve, 0 otherwise

public module subroutine print_bns_properties(this) Implementation →

Prints the $\mathrm{BNS}$ parameters to the standard output

Arguments

public module subroutine print_summary_bnsfuka(this, filename) Implementation →

Prints a summary of the physical properties of the $\mathrm{BNS}$ produced by $\texttt{FUKA}$ to the standard output and, optionally, to a formatted file whose name is given as the optional argument filename

Arguments

public module subroutine read_bns_properties(this) Implementation →

Imports the $\mathrm{BNS}$ parameters from $\texttt{FUKA}$

Arguments

public module subroutine read_fuka_id_full(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z) Implementation →

Stores the $\mathrm{ID}$ in non bnsfuka-member arrays with the same shape as the bnsfuka member arrays

Arguments

public module subroutine read_fuka_id_hydro(this, nx, ny, nz, pos, baryon_density, energy_density, specific_energy, pressure, u_euler) Implementation →

Stores the hydro $\mathrm{ID}$ in the arrays needed to compute the constraints on the refined mesh

Arguments

public module subroutine read_fuka_id_k(this, n, x, y, z, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz) Implementation →

Stores the components of the extrinsic curvature in arrays

Arguments

public module subroutine read_fuka_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler) Implementation →

Stores the hydro $\mathrm{ID}$ in the arrays needed to compute the baryon mass

Arguments

public module subroutine read_fuka_id_member(this, n, x, y, z) Implementation →

Stores the $\mathrm{ID}$ in the bnsfuka member arrays

Arguments

public module subroutine read_fuka_id_particles(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z) Implementation →

Stores the hydro $\mathrm{ID}$ in the arrays needed to compute the $\mathrm{SPH}$ $\mathrm{ID}$

Arguments

public module subroutine read_fuka_id_spacetime(this, nx, ny, nz, pos, lapse, shift, g, ek) Implementation →

Stores the spacetime $\mathrm{ID}$ in multi-dimensional arrays needed to compute the BSSN variables and constraints

Arguments

public module function read_fuka_mass_density(this, x, y, z) result(res) Implementation →

Returns the $\texttt{FUKA}$ baryon mass density at a point $(x,y,z)$

Arguments

Return Value double precision

Baryon mass density at $(x,y,z)$

public module function read_fuka_pressure(this, x, y, z) result(res) Implementation →

Returns the $\texttt{FUKA}$ pressure at a point $(x,y,z)$

Arguments

Return Value double precision

Pressure at $(x,y,z)$

public module function read_fuka_spatial_metric(this, x, y, z) result(res) Implementation →

Returns the $\texttt{FUKA}$ conformally flat spatial metric component $g_{xx}=g_{yy}=g_{zz}$ at a point $(x,y,z)$

Arguments

Return Value real(kind=C_DOUBLE)

$g_{xx}=g_{yy}=g_{zz}$ at $(x,y,z)$

public module subroutine run_kadath_reader(this, mpi_ranks, nx, ny, nz, xmin, xmax, ymin, ymax, zmin, zmax, coords, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, mass_density, specific_energy, pressure, v_eul_x, v_eul_y, v_eul_z, filename) Implementation →

Calls the MPI-parallelized vsion of the function KadathExportBNS from Kadath

Arguments

public module subroutine set_up_lattices_around_stars(this) Implementation →

Sets up two fine lattice, one around each star, to be able to interpolate the $\mathrm{ID}$ at the particle positions. It calls run_kadath_reader

Arguments