Hydro_nalysis

Hydro_Analysis

Frank Löffler

December 12, 2024

Abstract

Provides basic hydro analysis routines and quantities.

1 Purpose

2 Parameters


hydro_analysis_average_multiple_maxima_locations	Scope: private	BOOLEAN

Description: when finding more than one global maximum location, average position and use result

		Default: false


hydro_analysis_comp_core_rho_centroid	Scope: private	BOOLEAN

Description: compute location of the centroid of rho*x in region r_core around densest point

		Default: no


hydro_analysis_comp_rho_max	Scope: private	BOOLEAN

Description: Look for the value and location of the maximum of rho

		Default: false


hydro_analysis_comp_rho_max_every	Scope: private	INT

Description: How often to look for the value and location of the maximum of rho

Range		Default: 1
0:0	Never
1:*	every so often


hydro_analysis_comp_rho_max_origin_distance	Scope: private	BOOLEAN

Description: Look for the proper distance between the maximum of the density and the origin (along a straight coordinate line)

		Default: false


hydro_analysis_comp_vol_weighted_center_of_mass	Scope: private	BOOLEAN

Description: Look for the location of the volume-weighted center of mass

		Default: false


hydro_analysis_core_rho_rel_min	Scope: private	REAL

Description: only include points where rho>rho_rel_minrho_max when computing centroid of rhox around densest point

Range		Default: 0.01
0:*	any positive value. For best results should be such that the region selected is smaller than r_core. 1/10 of rho_max should be fine.


hydro_analysis_interpolator_coordinates	Scope: private	STRING

Description: Coordinate system

Range		Default: cart3d
.*	must be a registered coordinate system


hydro_analysis_interpolator_name	Scope: private	STRING

Description: Name of the interpolator

Range		Default: uniform cartesian
see [1] below	from AEILocalInterp
see [1] below	from AEILocalInterp
see [1] below	from AEILocalInterp
uniform cartesian	from LocalInterp
.*	must be a registered interpolator

[1]

Lagrange polynomial interpolation (tensor product)

[1]

Lagrange polynomial interpolation (maximum degree)

[1]

Hermite polynomial interpolation


hydro_analysis_interpolator_options	Scope: private	STRING

Description: Options for the interpolator

Range		Default: order=2
.*	must be a valid option specification


hydro_analysis_r_core	Scope: private	REAL

Description: size of region around densest point in which to compute the centroid of rho*x

Range		Default: 4.0
0:*	any positive radius. suggested is 1/2 of stellar radius


hydro_analysis_rho_max_loc_only_positive_x	Scope: private	BOOLEAN

Description: Restrict location search for density maximum to positive values of x

		Default: false


hydro_analysis_rho_max_loc_use_rotatingsymmetry180	Scope: private	BOOLEAN

Description: Map found maxima into positive x half-plane assuming pi-symmetry

		Default: false


hydro_analysis_rho_max_origin_distance_npoints	Scope: private	INT

Description: Number of points along the straight line for measuring proper distance

Range		Default: 100
1:*	Any positive number


verbosity_level	Scope: private	INT

Description: how much information to ouptut to the logs

Range		Default: 1
	Output nothing
1	warn when finding multiple maxima
2	also output location of maxima


restmass_compute_masses	Scope: restricted	BOOLEAN

Description: Should we compute the masses?

		Default: no


restmass_masses_nr	Scope: restricted	INT

Description: number of radii within which to compute the rest mass

Range		Default: (none)
0:100	Positive and hard-limited


restmass_ref_radius_mass	Scope: restricted	REAL

Description: Radii within which mass will be computed

Range		Default: 10.0
0.0:	Positive


restmass_rho_min	Scope: restricted	REAL

Description: Points with rest-mass density below this value are excluded (it could be set equal to the atmosphere value)

Range		Default: 1.e-9
0:


timelevels	Scope: shared from HYDROBASE	INT

3 Interfaces

General

Implements:

hydro_analysis

Inherits:

grid

hydrobase

admbase

Grid Variables

3.0.1 PRIVATE GROUPS


Group Names	Variable Names	Details

grid_spacing_product	grid_spacing_product	compact	0
		description	product of cctk_delta_space
		description	to be computed in local mode and later used in global mode
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		tags	checkpoint=”no”
		timelevels	1
		variable type	REAL

hydro_analysis_masses	Hydro_Analysis_masses	compact	0
		description	Baryonic masses at different radii
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		tags	checkpoint=”no”
		timelevels	1
		vararray_size	restmass_masses_nr
		variable type	REAL

hydro_analysis_masses_fractions	Hydro_Analysis_masses_fractions	compact	0
		description	Fractional Baryonic masses at different radii
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		tags	checkpoint=”no”
		timelevels	1
		vararray_size	restmass_masses_nr
		variable type	REAL

hydro_analysis_total_rest_mass		compact	0
	total_rest_mass	description	Total Baryonic mass
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		tags	checkpoint=”no”
		timelevels	1
		variable type	REAL

hydro_analysis_masses_temps	Hydro_Analysis_masses_temps	compact	0
		description	Temporaries for the mass calculation
		dimensions	3
		distribution	DEFAULT
		group type	GF
		tags	checkpoint=”no”
		timelevels	3
		vararray_size	restmass_masses_nr+1
		variable type	REAL

3.0.2 PUBLIC GROUPS


Group Names	Variable Names	Details

hydro_analysis_rho_max	Hydro_Analysis_rho_max	compact	0
		description	value of the maximum of rho
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		variable type	REAL

hydro_analysis_rho_sum	Hydro_Analysis_rho_sum	compact	0
		description	value of the sum of rho
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		variable type	REAL

hydro_analysis_rho_max_loc	Hydro_Analysis_rho_max_loc	compact	0
		description	coordinate location of the maximum of rho
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		vararray_size	3
		variable type	REAL

hydro_analysis_rho_center_volume_weighted_gf	Hydro_Analysis_rho_center_volume_weighted_gf	compact	0
		description	temporary GF to obtain the coordinate location of the volume weighted center of mass via a reduction
		dimensions	3
		distribution	DEFAULT
		group type	GF
		tags	checkpoint=”no” tensortypealias=”u”
		timelevels	3
		vararray_size	3
		variable type	REAL

hydro_analysis_rho_center_volume_weighted	Hydro_Analysis_rho_center_volume_weighted	compact	0
		description	coordinate location of the volume weightes center of mass
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		vararray_size	3
		variable type	REAL

hydro_analysis_rho_max_origin_distance	Hydro_Analysis_rho_max_origin_distance	compact	0
		description	proper distance between the maximum of the density and the origin (along a straight coordinate line)
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		variable type	REAL


Group Names	Variable Names	Details

hydro_analysis_core_rho_sum	Hydro_Analysis_core_rho_sum	compact	0
		description	value of the sum of rho in the core region
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		variable type	REAL

hydro_analysis_core_rho_centroid_gf	Hydro_Analysis_core_rho_centroid_gf	compact	0
		description	temporary GF to obtain the coordinate location of the centroid location of rho*x in the core region via a reduction
		dimensions	3
		distribution	DEFAULT
		group type	GF
		tags	checkpoint=”no” tensortypealias=”4u”
		timelevels	3
		vararray_size	4
		variable type	REAL

hydro_analysis_core_rho_centroid	Hydro_Analysis_core_rho_centroid	compact	0
		description	coordinate location of the centroid of rho*x in the core region
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		vararray_size	3
		variable type	REAL

4 Schedule

This section lists all the variables which are assigned storage by thorn EinsteinAnalysis/Hydro_Analysis. Storage can either last for the duration of the run (Always means that if this thorn is activated storage will be assigned, Conditional means that if this thorn is activated storage will be assigned for the duration of the run if some condition is met), or can be turned on for the duration of a schedule function.

Storage

	Conditional:
	Hydro_Analysis_rho_max
	Hydro_Analysis_rho_max_loc
	Hydro_Analysis_total_rest_mass
	Hydro_Analysis_masses_fractions
	grid_spacing_product
	Hydro_Analysis_rho_sum
	Hydro_Analysis_rho_center_volume_weighted
	Hydro_Analysis_rho_center_volume_weighted_gf[3]
	Hydro_Analysis_rho_max_origin_distance
	Hydro_Analysis_core_rho_sum
	Hydro_Analysis_core_rho_centroid
	Hydro_Analysis_masses_temps[3]
	Hydro_Analysis_masses

Scheduled Functions

CCTK_BASEGRID (conditional)

hydro_analysis_init

initialize variables

	Language:	c
	Options:	global
	Type:	function

CCTK_POSTSTEP (conditional)

hydro_analysis

group for hydro_analysis routines

	Type:	group

Hydro_Analysis_CompCoreRhoCentroid (conditional)

hydro_analysis_compcorerhocentroid_getstorage

get temporary storage for duration of reduction

	Before:	hydro_analysis_compcorerhocentroid_preparereduction
	Language:	c
	Options:	global
	Type:	function

Hydro_Analysis_CompCoreRhoCentroid (conditional)

hydro_analysis_compcorerhocentroid_preparereduction

prepare core center of mass densities for reduction

	Language:	c
	Options:	global
		loop-local
	Type:	function

Hydro_Analysis_CompCoreRhoCentroid (conditional)

hydro_analysis_compcorerhocentroid_reduction

compute the global core center of mass reduction results

	After:	hydro_analysis_compcorerhocentroid_preparereduction
	Language:	c
	Options:	global
	Type:	function

Hydro_Analysis_CompCoreRhoCentroid (conditional)

hydro_analysis_compcorerhocentroid_freestorage

free temporary storage

	After:	hydro_analysis_compcorerhocentroid_reduction
	Language:	c
	Options:	global
	Type:	function

CCTK_POSTSTEP (conditional)

hydro_analysis_masses

computation of rest masses contained within given radii

	Type:	group

Hydro_Analysis_Masses (conditional)

hydro_analysis_masses_local

rest masses within given radii, local calculations

	Language:	fortran
	Type:	function

Hydro_Analysis_Masses (conditional)

hydro_analysis_masses_global

rest masses within given radii, reduction operation

	After:	hydro_analysis_masses_local
	Language:	fortran
	Options:	global
	Type:	function

Hydro_Analysis (conditional)

hydro_analysis_preparereduction

compute the local reduction results

	Language:	c
	Options:	global-early
		loop-local
	Type:	function

Hydro_Analysis (conditional)

hydro_analysis_reduction

compute the global reduction results

	After:	hydro_analysis_preparereduction
	Language:	c
	Options:	global
	Type:	function

Hydro_Analysis (conditional)

hydro_analysis_locationsearch

look for the location of the maximum density

	After:	hydro_analysis_reduction
	Type:	group

Hydro_Analysis_LocationSearch (conditional)

hydro_analysis_locationsearch_setup

prepare data structures for search

	Language:	c
	Options:	global
	Type:	function

Hydro_Analysis_LocationSearch (conditional)

hydro_analysis_locationsearch_search

search for the location of the maximum density

	After:	hydro_analysis_locationsearch_setup
	Language:	c
	Options:	global
		loop-local
	Type:	function

Hydro_Analysis_LocationSearch (conditional)

hydro_analysis_locationsearch_combine

communicate and verify the location of the maximum density

	After:	hydro_analysis_locationsearch_search
	Language:	c
	Options:	global
	Type:	function

Hydro_Analysis (conditional)

hydro_analysis_findseparation

compute the proper distance between the maximum of the density and the origin (along a straight coordinate line)

	After:	hydro_analysis_locationsearch
	Language:	fortran
	Options:	global
	Type:	function

Hydro_Analysis (conditional)

hydro_analysis_compcorerhocentroid

compute center of mass of core region

	After:	hydro_analysis_locationsearch
	Type:	group