Hydro_Analysis

Frank Löffler

December 4, 2021

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_min*rho_max when computing centroid of rho*x 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   compact0
  descriptionproduct of cctk_delta_space
    descriptionto be computed in local mode and later used in global mode
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  tagscheckpoint=”no”
  timelevels1
 variable typeREAL




hydro_analysis_masses Hydro_Analysis_masses   compact0
  descriptionBaryonic masses at different radii
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  tagscheckpoint=”no”
  timelevels1
 vararray_sizerestmass_masses_nr
 variable typeREAL




hydro_analysis_masses_fractionsHydro_Analysis_masses_fractions  compact0
  descriptionFractional Baryonic masses at different radii
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  tagscheckpoint=”no”
  timelevels1
 vararray_sizerestmass_masses_nr
 variable typeREAL




hydro_analysis_total_rest_mass   compact0
total_rest_mass   descriptionTotal Baryonic mass
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  tagscheckpoint=”no”
  timelevels1
 variable typeREAL




hydro_analysis_masses_temps Hydro_Analysis_masses_temps   compact0
  descriptionTemporaries for the mass calculation
  dimensions3
  distributionDEFAULT
  group typeGF
  tagscheckpoint=”no”
  timelevels3
 vararray_sizerestmass_masses_nr+1
 variable typeREAL




3.0.2 PUBLIC GROUPS




  Group Names     Variable Names     Details   




hydro_analysis_rho_max Hydro_Analysis_rho_max   compact0
  descriptionvalue of the maximum of rho
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 variable typeREAL




hydro_analysis_rho_sum Hydro_Analysis_rho_sum   compact0
  descriptionvalue of the sum of rho
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 variable typeREAL




hydro_analysis_rho_max_loc Hydro_Analysis_rho_max_loc   compact0
  descriptioncoordinate location of the maximum of rho
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 vararray_size3
 variable typeREAL




hydro_analysis_rho_center_volume_weighted_gfHydro_Analysis_rho_center_volume_weighted_gf  compact0
  descriptiontemporary GF to obtain the coordinate location of the volume weighted center of mass via a reduction
  dimensions3
  distributionDEFAULT
  group typeGF
  tagscheckpoint=”no” tensortypealias=”u”
  timelevels3
 vararray_size3
 variable typeREAL




hydro_analysis_rho_center_volume_weighted Hydro_Analysis_rho_center_volume_weighted   compact0
  descriptioncoordinate location of the volume weightes center of mass
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 vararray_size3
 variable typeREAL




hydro_analysis_rho_max_origin_distance Hydro_Analysis_rho_max_origin_distance   compact0
  descriptionproper distance between the maximum of the density and the origin (along a straight coordinate line)
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 variable typeREAL








  Group Names     Variable Names     Details   




hydro_analysis_core_rho_sum Hydro_Analysis_core_rho_sum   compact0
  descriptionvalue of the sum of rho in the core region
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 variable typeREAL




hydro_analysis_core_rho_centroid_gfHydro_Analysis_core_rho_centroid_gf  compact0
  descriptiontemporary GF to obtain the coordinate location of the centroid location of rho*x in the core region via a reduction
  dimensions3
  distributionDEFAULT
  group typeGF
  tagscheckpoint=”no” tensortypealias=”4u”
  timelevels3
 vararray_size4
 variable typeREAL




hydro_analysis_core_rho_centroid Hydro_Analysis_core_rho_centroid   compact0
  descriptioncoordinate location of the centroid of rho*x in the core region
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 vararray_size3
 variable typeREAL




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