diffstar_lorene Module

public module subroutine allocate_diffstar_memory(this, d) Implementation →

Allocates allocatable arrays member of a diffstarlorene object

Arguments

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

Constructs a diffstarlorene object 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 construct_drs(this, id_file, eos_filename) Implementation →

Interface of the subroutine that constructs the \(\texttt{LORENE}\) Etdiffrot object

Arguments

private function construct_etdiffrot(c_id_file, c_eos_file) result(optr) bind(C, NAME= "construct_et_diffrot")

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that constructs the \(\texttt{LORENE}\) \(\texttt{Et_diffrot}\) object

FT 24.10.2021

Arguments

Return Value type(C_PTR)

C pointer pointing to the constructed \(\texttt{LORENE}\) \(\texttt{Et_diffrot}\) object

public module subroutine deallocate_diffstar_memory(this) Implementation →

Deallocates allocatable arrays member of a diffstarlorene object

Arguments

public module subroutine destruct_diffstarlorene(this) Implementation →

Destruct a diffstarlorene object

Arguments

public module subroutine destruct_drs(this) Implementation →

Destructs a \(\texttt{LORENE}\) Etdiffrot object

Arguments

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

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that destructs the \(\texttt{LORENE}\) \(\texttt{Et_diffrot}\) object

FT 24.10.2021

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

private subroutine get_diffstar_full(optr, x, y, z, lapse, shift_x, shift_y, shift_z, g_rr, g_tt, g_pp, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, baryon_density, energy_density, specific_energy, pressure, v_euler_x, v_euler_y, v_euler_z) bind(C, NAME= "0")

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that reads the full \(\texttt{LORENE}\) ID at the specified point. That is, imports 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 24.10.2021

Arguments

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

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that reads the hydro fields from \(\texttt{LORENE}\), 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 24.10.2021

Arguments

public module function get_diffstar_identifier(this) Implementation →

Arguments

Return Value double precision

private subroutine get_diffstar_mass_b(optr, x, y, z, g_rr, g_tt, g_pp, baryon_density, gamma_euler) bind(C, NAME= "0")

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that reads the hydro fields and the metric fields from \(\texttt{LORENE}\), 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 24.10.2021

Arguments

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

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that returns the baryon mass density $[\mathrm{kg}\, \mathrm{m}^{-3}]$ from \(\texttt{LORENE}\), at the specified point

FT 24.10.2021

Arguments

Return Value real(kind=C_DOUBLE)

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

private subroutine get_diffstar_params(optr, omega_c, mass, mass_grav, angular_momentum, tsw, grv2, grv3, r_circ, surface_area, r_mean, r_eq, r_eq_pi2, r_eq_pi, r_eq_3pi2, r_eq_pole, r_ratio, r_isco, f_isco, specific_energy_isco, specific_angular_momentum_isco, area_radius, ent_center, nbar_center, rho_center, energy_density_center, specific_energy_center, pressure_center, redshift_eqf, redshift_eqb, redshift_pole, eos, eos_id, gamma, kappa, npeos, gamma0, gamma1, gamma2, gamma3, kappa0, kappa1, kappa2, kappa3, logP1, logRho0, logRho1, logRho2, eos_table) bind(C, NAME= "0")

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that stores the physical parameters of the binary system from \(\texttt{LORENE}\) in the desired variables

FT 24.10.2021

Arguments

private subroutine get_diffstar_particles(optr, x, y, z, lapse, shift_x, shift_y, shift_z, g_rr, g_tt, g_pp, baryon_density, energy_density, specific_energy, pressure, v_euler_x, v_euler_y, v_euler_z) bind(C, NAME= "0")

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that reads the hydro fields and the metric fields * from \(\texttt{LORENE}\), 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 24.10.2021

Arguments

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

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that returns the baryon mass density $[\mathrm{kg}\, c^2 \mathrm{m}^{-3}]$ from \(\texttt{LORENE}\), at the specified point

FT 06.12.2022

Arguments

Return Value real(kind=C_DOUBLE)

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

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

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

FT 24.10.2021

Arguments

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

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that returns the diagonal components of the metric, all equal to the \(\texttt{LORENE}\) conformal factor to the 4th power.

FT 24.10.2021

Arguments

Return Value real(kind=C_DOUBLE)

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

public module function get_eos_loreneid(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 diffstarlorene member arrays named field

Arguments

Return Value double precision, DIMENSION(:), ALLOCATABLE

Desired diffstarlorene member array

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

Returns the component n of the diffstarlorene member arrays named field

Arguments

Return Value double precision

Component n of the desired diffstarlorene member array

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

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

Arguments

Return Value logical

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

public module subroutine nothing(this, flag, switch) Implementation →

Procedure that does nothing. It is used to instantiate a deferred idbase procedure which is not needed in TYPE diffstarlorene. It also serves as a placeholder in case the idbase procedure will be needed in the future.

Arguments

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

Interface to the \(\texttt{LORENE}\) method of class \(\texttt{Et_diffrot}\) with the same name, that returns 1 if the energy density is positive, and if the specific energy is positive, and if the pressure is nonpositive, at the specified point; it returns 0 otherwise

FT 24.10.2021

Arguments

Return Value integer(kind=C_INT)

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

public module subroutine print_diffstar_properties(this) Implementation →

Prints the DRS parameters to the standard output

Arguments

public module subroutine read_diffstar_properties(this) Implementation →

Imports the DRS parameters from \(\texttt{LORENE}\)

Arguments

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

Returns the \(\texttt{LORENE}\) 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_drslorene_pressure(this, x, y, z) result(res) Implementation →

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

Arguments

Return Value double precision

Pressure at $(x,y,z)$

public module subroutine read_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 ID in non diffstarlorene-member arrays with the same shape as the diffstarlorene member arrays

Arguments

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

Stores the hydro ID in the arrays needed to compute the constraints on the refined mesh

Arguments

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

Stores the ID in the diffstarlorene member arrays

Arguments

public module subroutine read_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_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler) Implementation →

Stores the hydro ID in the arrays needed to compute the baryon mass

Arguments

public module subroutine read_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 ID in the arrays needed to compute the SPH ID

Arguments

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

Stores the spacetime ID in multi-dimensional arrays needed to compute the BSSN variables and constraints

Arguments

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

Returns the \(\texttt{LORENE}\) 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)$