Hydro_Analysis

Frank Löffler

April 20, 2017

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






timelevels
Scope: shared from HYDROBASE INT



3 Interfaces

General

Implements:

hydro_analysis

Inherits:

grid

hydrobase

Grid Variables

3.0.1 PUBLIC GROUPS




  Group Names     Variable Names     Details   




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




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




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




hydro_analysis_rho_center_volume_weighted_gf   compact0
Hydro_Analysis_rho_center_volume_weighted_gf  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   compact0
Hydro_Analysis_rho_center_volume_weighted   descriptioncoordinate location of the volume weightes center of mass
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 vararray_size3
 variable typeREAL




hydro_analysis_rho_max_origin_distance   compact0
Hydro_Analysis_rho_max_origin_distance   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   compact0
Hydro_Analysis_core_rho_sum   descriptionvalue of the sum of rho in the core region
  dimensions0
  distributionCONSTANT
  group typeSCALAR
  timelevels1
 variable typeREAL




hydro_analysis_core_rho_centroid_gf   compact0
Hydro_Analysis_core_rho_centroid_gf  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   compact0
Hydro_Analysis_core_rho_centroid   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_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
   

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

Hydro_Analysis (conditional)

  hydro_analysis_preparereduction

  compute the local reduction results

 

 Language:c
 Options: global-early
 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