NullEvolve

Y.Z. <yosef@raven.phyast.pitt.edu>

January 09 2003

Abstract

1 Introduction

2 Physical System

3 Numerical Implementation

4 Using This Thorn

4.1 Obtaining This Thorn

4.2 Basic Usage

4.3 Special Behaviour

4.4 Interaction With Other Thorns

4.5 Examples

4.6 Support and Feedback

5 History

5.1 Thorn Source Code

5.2 Thorn Documentation

5.3 Acknowledgements

References

6 Parameters


dissip_fudge	Scope: private	INT

Description: do we want to fudge around the inner boundary

Range		Default: (none)
0:1	0 no fudging, 1 fudging


dissip_fudge_maxx	Scope: private	REAL

Description: where in x to stop fudging

Range		Default: .53
.5:.56


dissip_j	Scope: private	REAL

Description: dissipation factor for the J equation (time-integration)

Range		Default: 0.0
0.0:*	positive


dissip_jx	Scope: private	REAL

Description: dissipation factor for the J equation (radial marching)

Range		Default: 0.0
0.0:*	positive


dissip_mask_type	Scope: private	KEYWORD

Description: type of dissipation mask

Range		Default: one
one	set dissipation mask to ’one’ everywhere
see [1] below	”set 2B=1-(p+q)+—1-(p +q)—”
see [1] below	”set 2B=rD0-(p+q)+—rD 0-(p+q)—”
see [1] below	”set 2B=min(1, (rD0-(p+q)+—rD0- (p+q)—) / (rD0-1) )”
see [1] below	”set 2B=min(1, (rD0-(p+q)+—rD0- (p+q)—) / (rD0-rD1) )”

[1]

zero at eq, one at pole

[1]

zero at rD0, one at pole

[1]

zero at rD0, one at eq

[1]

zero at rD0, one at rD1


dissip_q	Scope: private	REAL

Description: dissipation factor for the Q equation (the auxiliary of U)

Range		Default: 0.0
0.0:*	positive


dissip_w	Scope: private	REAL

Description: dissipation factor for the W equation

Range		Default: 0.0
0.0:*	positive


id_amp	Scope: private	REAL

Description: Amplitude of compact J pulse

Range		Default: .00001
0.:*	An y positive number


id_l	Scope: private	INT

Description: l mode

Range		Default: 2
2:*	2 or larger


id_m	Scope: private	INT

Description: m mode

Range		Default: (none)
-100:100	-l to l


id_power	Scope: private	INT

Description: related to order of polynomial in compact J pulse

Range		Default: 3
3:6	best to choose between 3 and 6


jcoeff_r0	Scope: private	REAL

Description: real part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r0i	Scope: private	REAL

Description: imag part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r0r	Scope: private	REAL

Description: real part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r1	Scope: private	REAL

Description: imag part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r1i	Scope: private	REAL

Description: imag part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r1r	Scope: private	REAL

Description: real part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r3	Scope: private	REAL

Description: real part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r3i	Scope: private	REAL

Description: imag part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


jcoeff_r3r	Scope: private	REAL

Description: real part of constant for fitted linearized J initial data

Range		Default: (none)
:	any value obtained from Nigel’s fitting script


n_dissip_one_outside_eq	Scope: private	REAL

Description: determines dissipation radius rD1 = 1 + dd*N_dissip_one_outside_eq

Range		Default: (none)
:	any value


n_dissip_zero_outside_eq	Scope: private	REAL

Description: determines dissipation radius rD0 = 1 + dd*N_dissip_zero_outside_eq

Range		Default: 1
0:*	non-negative


null_dissip	Scope: private	REAL

Description: dissipation factor for the radial evolution stencil

Range		Default: 0.05
0.0:1.0	postive less than 1


psmax	Scope: private	REAL

Description: outer p bdry of CS pulse

Range		Default: 1.0
:	be reasonable though


psmin	Scope: private	REAL

Description: inner p bdry of CS pulse

Range		Default: -1.0
:	be reasonable though


qsmax	Scope: private	REAL

Description: outer q bdry of CS pulse

Range		Default: 1.0
:	be reasonable though


qsmin	Scope: private	REAL

Description: inner q bdry of CS pulse

Range		Default: -1.0
:	be reasonable though


use_rsylm	Scope: private	BOOLEAN

Description: do we use real spherical harmonics for compact pulse initial data

		Default: yes


wrot	Scope: private	REAL

Description: frequency of CS pule

Range		Default: 7.5
:	be reasonable though


xmax	Scope: private	REAL

Description: right boundary of CS pulse

Range		Default: .8
.5:1	choose .5 < xmin < xmax < 1


xmin	Scope: private	REAL

Description: left boundary of CS pulse

Range		Default: .56
.5:1	choose .5 < xmin < xmax < 1


boundary_data	Scope: restricted	KEYWORD

Description: Choose boundary data type

Range		Default: flat
whitehole	White hole boundary data
flat	flat boundary data
randomJ	random number boundary data


debug_skip_b_update	Scope: restricted	BOOLEAN

Description: Should the update of B be turned off?

		Default: no


debug_skip_cb_update	Scope: restricted	BOOLEAN

Description: Should the update of CB be turned off?

		Default: no


debug_skip_ck_update	Scope: restricted	BOOLEAN

Description: Should the update of CK be turned off?

		Default: no


debug_skip_evolution	Scope: restricted	BOOLEAN

Description: Should the evolution be turned off?

		Default: no


debug_skip_j_update	Scope: restricted	BOOLEAN

Description: Should the update of J be turned off?

		Default: no


debug_skip_nu_update	Scope: restricted	BOOLEAN

Description: Should the update of NU be turned off?

		Default: no


debug_skip_q_update	Scope: restricted	BOOLEAN

Description: Should the update of Q be turned off?

		Default: no


debug_skip_u_update	Scope: restricted	BOOLEAN

Description: Should the update of U be turned off?

		Default: no


debug_skip_w_update	Scope: restricted	BOOLEAN

Description: Should the update of W be turned off?

		Default: no


debug_verbose	Scope: restricted	BOOLEAN

Description: should debugging messages be printed?

		Default: no


diagnostics_coord_x	Scope: restricted	REAL

Description: the coordinate x at which diagnostics are to be measured

Range		Default: 0.0
0:1	will test only for x >= NullGrid::null_xin.


first_order_scheme	Scope: restricted	BOOLEAN

Description: Should the first order (angular) scheme be used?

		Default: yes


initial_j_data	Scope: restricted	KEYWORD

Description: What kind of initial data for J shall we pick?

Range		Default: vanishing_J_scri
vanishing_J	vanishing J
vanishing_J_scri	vanishing J at scri
rotating_pulse	rotating pulse
compact_J	compact suport J
smooth_J	smoothly vanishing J
extracted_J	J as given by extraction
polynomial_J	J as third power polynomial to vanish J_xx at scri
see [1] below	linearized fitted J based on some previous run

[1]

fitted\_linearized\_J


old_j_xderiv	Scope: restricted	BOOLEAN

Description: should we compute the x derivative of J with the old values?

		Default: no


timestep_method	Scope: shared from TIME	KEYWORD

7 Interfaces

General

Implements:

nullevolve

Inherits:

nullinterp

nullgrid

nullvars

time

Grid Variables

7.0.1 PRIVATE GROUPS


Group Names	Variable Names	Details

distmp		compact	0
	distmp_F	description	temporaries used for dissipation
	distmp_d2F	dimensions	2
	distmp_d4F	distribution	DEFAULT
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		group type	ARRAY
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		timelevels	1
		vararray_size	2
		variable type	COMPLEX

dissip_mask	dissip_mask	compact	0
		description	dissipation mask
		dimensions	2
		distribution	DEFAULT
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		group type	ARRAY
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		timelevels	1
		variable type	REAL

diagtmp	diagtmp	compact	0
		description	temporaries used for diagnostics
		dimensions	2
		distribution	DEFAULT
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		group type	ARRAY
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		timelevels	1
		vararray_size	2
		variable type	COMPLEX

aux_mask		compact	0
	auxiliary_maskn	description	auxiliary evolution masks
	auxiliary_masks	dimensions	2
		distribution	DEFAULT
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		group type	ARRAY
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		timelevels	1
		variable type	INT

eth4_mask		compact	0
	eth4_maskn	description	dissipation masks
		description	to protect the startup from taking eth4
	eth4_masks	dimensions	2
		distribution	DEFAULT
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		ghostsize	NULLGRID::N_ANG_GHOST_PTS
		group type	ARRAY
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		size	(NULLGRID::N_ANG_PTS_INSIDE_EQ+2*(NULLGRID::N_ANG_EV_OUTSIDE_EQ+NULLGRID::N_ANG_STENCIL_SIZE))
		timelevels	1
		variable type	INT

jrad		compact	0
	jcn_rad	dimensions	1
	jcs_rad	distribution	CONSTANT
		group type	ARRAY
		size	NULLGRID::N_RADIAL_PTS
		timelevels	1
		variable type	COMPLEX


Group Names	Variable Names	Details

dxjrad		compact	0
	dxjcn_rad	dimensions	1
	dxjcs_rad	distribution	CONSTANT
		group type	ARRAY
		size	NULLGRID::N_RADIAL_PTS
		timelevels	1
		variable type	COMPLEX

8 Schedule

This section lists all the variables which are assigned storage by thorn PITTNullCode/NullEvolve. 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

Always:	Conditional:
NullVars::realcharfuncs[2]	NullVars::cmplxcharfuncs_aux[2]
NullVars::cmplxcharfuncs_basic[2]
NullVars::null_mask
NullGrid::EG_mask NullGrid::EQ_mask NullGrid::EV_mask
eth4_mask dissip_mask
diagtmp aux_mask
Jrad dxJrad

Scheduled Functions

CCTK_POSTINITIAL (conditional)

nullevol_initial

null init data

	After:	harmidata_init_to_adm
		harmidata_pinit_to_adm
		mol_fillalllevels
		adm_bssn_init
	Options:	global
	Type:	group

CCTK_INITIAL (conditional)

nullevol_initializearrays

initialize all arrays to large values

	Before:	nullevol_initial
	Language:	fortran
	Type:	function

NullEvol_BoundaryInit (conditional)

nullevol_bdry_flat

give flat boundary conditions for the null metric

	Language:	fortran
	Options:	global
	Type:	function

NullEvol_BoundaryInit (conditional)

nullevol_bdry_whitehole

give white hole boundary conditions for the null metric

	Language:	fortran
	Options:	global
	Type:	function

NullEvol_BoundaryInit (conditional)

nullevol_bdry_randomj

give random j boundary conditions for the null metric

	Language:	fortran
	Options:	global
	Type:	function

NullEvol_Boundary (conditional)

nullevol_bdry_flat

give flat boundary conditions for the null metric

	Language:	fortran
	Options:	global
	Type:	function

NullEvol_Boundary (conditional)

nullevol_bdry_whitehole

give white hole boundary conditions for the null metric

	Language:	fortran
	Options:	global
	Type:	function

NullEvol_Boundary (conditional)

nullevol_bdry_randomj

give random j boundary conditions for the null metric

	Language:	fortran
	Options:	global
	Type:	function

CCTK_EVOL (conditional)

nullevol_resettop

reset top values

	Before:	nullevol_boundary
	Language:	fortran
	Options:	global
	Type:	function

CCTK_EVOL (conditional)

nullevol_boundary

boundary data

	After:	harmevol_to_adm
		mol_evolution
	Options:	global
	Type:	group

CCTK_EVOL (conditional)

nullevol_step

evolution

	After:	nullevol_boundary
	Language:	fortran
	Options:	global
	Storage:	distmp
	Type:	function

NullEvol_Initial (conditional)

nullevol_initialdata

give j on the initial null hypersurface

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

NullEvol_Initial (conditional)

nullevol_boundaryinit

boundary data for the characteristic data

	Before:	nullevol_initialdata
		nullevol_initialslice
	Type:	group

NullEvol_Initial (conditional)

nullevol_initialslice

construct null metric on the initial null hypersurface

	After:	nullevol_boundaryinit
		nullevol_initialdata
	Language:	fortran
	Options:	global
	Type:	function

NullEvol_Initial (conditional)

nullevol_diag

diagnostics of the characteristic code

	After:	nullevol_initialslice
	Language:	fortran
	Options:	global
	Storage:	diagtmp
	Type:	function

CCTK_EVOL (conditional)

nullevol_diag

diagnostics of the characteristic code

	After:	nullevol_step
	Language:	fortran
	Options:	global
	Storage:	diagtmp
	Type:	function