compute_and_print_bssn_constraints_grid Module Procedure

module procedure compute_and_print_bssn_constraints_grid module subroutine compute_and_print_bssn_constraints_grid(this, id, namefile, name_logfile, points)

Uses

    • mesh_refinement
    • matrix
    • utility
    • constants
    • McLachlan_refine
    • tensor
  • proc~~compute_and_print_bssn_constraints_grid~~UsesGraph proc~compute_and_print_bssn_constraints_grid compute_and_print_bssn_constraints_grid McLachlan_refine McLachlan_refine proc~compute_and_print_bssn_constraints_grid->McLachlan_refine constants constants proc~compute_and_print_bssn_constraints_grid->constants matrix matrix proc~compute_and_print_bssn_constraints_grid->matrix mesh_refinement mesh_refinement proc~compute_and_print_bssn_constraints_grid->mesh_refinement module~utility utility proc~compute_and_print_bssn_constraints_grid->module~utility tensor tensor proc~compute_and_print_bssn_constraints_grid->tensor module~utility->constants module~utility->matrix

Compute, store, analyze and print the \(\mathrm{BSSNOK}\) constraints to a formatted file. The computation is done by importing the hydro \(\mathrm{ID}\) on the gravity grid, without any information on the particles.

FT 1.02.2021

DEALLOCATE( u_coord ) DEALLOCATE( u_coord_l )

Arguments

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

bssn object to which this PROCEDURE is bound
class(idbase), intent(inout) :: id

idbase object used to read the hydrodynamical \(\mathrm{ID}\) to the mesh
character(len=*), intent(inout) :: namefile
character(len=*), intent(inout) :: name_logfile
double precision, intent(in), optional, DIMENSION(:,:,:,:), TARGET :: points

Calls

proc~~compute_and_print_bssn_constraints_grid~~CallsGraph proc~compute_and_print_bssn_constraints_grid compute_and_print_bssn_constraints_grid a_bssn3_ll a_bssn3_ll proc~compute_and_print_bssn_constraints_grid->a_bssn3_ll allocate_grid_function allocate_grid_function proc~compute_and_print_bssn_constraints_grid->allocate_grid_function deallocate_grid_function deallocate_grid_function proc~compute_and_print_bssn_constraints_grid->deallocate_grid_function g_bssn3_ll g_bssn3_ll proc~compute_and_print_bssn_constraints_grid->g_bssn3_ll gamma_u gamma_u proc~compute_and_print_bssn_constraints_grid->gamma_u gc gc proc~compute_and_print_bssn_constraints_grid->gc hc hc proc~compute_and_print_bssn_constraints_grid->hc interface~bssn_constraint_terms_interior bssn_constraint_terms_interior proc~compute_and_print_bssn_constraints_grid->interface~bssn_constraint_terms_interior k_phys3_ll k_phys3_ll proc~compute_and_print_bssn_constraints_grid->k_phys3_ll lapse lapse proc~compute_and_print_bssn_constraints_grid->lapse levels levels proc~compute_and_print_bssn_constraints_grid->levels mc mc proc~compute_and_print_bssn_constraints_grid->mc none~compute_4velocity_eul compute_4velocity_eul proc~compute_and_print_bssn_constraints_grid->none~compute_4velocity_eul none~compute_stress_energy compute_stress_energy proc~compute_and_print_bssn_constraints_grid->none~compute_stress_energy phi phi proc~compute_and_print_bssn_constraints_grid->phi rho rho proc~compute_and_print_bssn_constraints_grid->rho s s proc~compute_and_print_bssn_constraints_grid->s shift_u shift_u proc~compute_and_print_bssn_constraints_grid->shift_u trk trk proc~compute_and_print_bssn_constraints_grid->trk none~compute_4velocity_eul->lapse none~compute_4velocity_eul->shift_u g_phys3_ll g_phys3_ll none~compute_4velocity_eul->g_phys3_ll invert_4x4_matrix invert_4x4_matrix none~compute_4velocity_eul->invert_4x4_matrix proc~compute_g4 compute_g4 none~compute_4velocity_eul->proc~compute_g4 proc~determinant_sym4x4 determinant_sym4x4 none~compute_4velocity_eul->proc~determinant_sym4x4

Called by

proc~~compute_and_print_bssn_constraints_grid~~CalledByGraph proc~compute_and_print_bssn_constraints_grid compute_and_print_bssn_constraints_grid interface~compute_and_print_bssn_constraints_grid compute_and_print_bssn_constraints_grid interface~compute_and_print_bssn_constraints_grid->proc~compute_and_print_bssn_constraints_grid

Contents

Variables

Subroutines


Variables

Type Visibility Attributes Name Initial
type(grid_function_scalar), private :: HC_A
type(grid_function_scalar), private :: HC_derphi
type(grid_function_scalar), private :: HC_hand
type(grid_function_scalar), private :: HC_rho
type(grid_function_scalar), private :: HC_trK
type(grid_function), private :: Tmunu_ll
type(grid_function_scalar), private :: baryon_density
logical, private, parameter :: debug = .FALSE.
double precision, private :: detg4
type(grid_function_scalar), private :: energy_density
logical, private :: exist
integer, private :: fd_lim
character(len=:), private, ALLOCATABLE :: finalname_logfile
type(grid_function), private :: g4
double precision, private, DIMENSION(4, 4) :: g4temp
integer, private :: i
double precision, private, DIMENSION(4, 4) :: ig4
integer, private, DIMENSION(3) :: imax
integer, private, DIMENSION(3) :: imin
integer, private :: j
integer, private :: k
integer, private :: l
type(grid_function_scalar), private :: lorentz_factor
double precision, private :: min_abs_y
double precision, private :: min_abs_z
character(len=2), private :: n_reflev
character(len=:), private, ALLOCATABLE :: name_analysis
character(len=:), private, ALLOCATABLE :: name_constraint
integer, private :: nx
integer, private :: ny
integer, private :: nz
type(grid_function_scalar), private :: pressure
double precision, private, DIMENSION(:,:,:), POINTER :: pts_x
double precision, private, DIMENSION(:,:,:), POINTER :: pts_y
double precision, private, DIMENSION(:,:,:), POINTER :: pts_z
type(grid_function_scalar), private :: specific_energy
character(len=2), private :: tpo_id
type(grid_function), private :: u_euler_l
double precision, private :: u_euler_norm = zero
integer, private :: unit_logfile
type(grid_function), private :: v_euler
type(grid_function), private :: v_euler_l

Subroutines

subroutine compute_4velocity_eul()

Compute the components of the fluid -velocity wrt the Eulerian observer See Sec.7.3 in Alcubierre, "Introduction to 3+1 Numerical Relativity"

FT 25.04.2022

Arguments

None

subroutine compute_stress_energy()

Compute the components of the stress-energy tensor

FT 25.04.2022

Arguments

None