CarpetProlongateTest

Erik Schnetter <schnetter@cct.lsu.edu>

July 10 2010

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


interp_nx	Scope: private	INT

Description: number of points for interpolated region

Range		Default: 4
0:*


interp_ny	Scope: private	INT

Description: number of points for interpolated region

Range		Default: 4
0:*


interp_nz	Scope: private	INT

Description: number of points for interpolated region

Range		Default: 4
0:*


interp_xmax	Scope: private	REAL

Description: xmax for interpolated region

Range		Default: +0.5
:


interp_xmin	Scope: private	REAL

Description: xmin for interpolated region

Range		Default: -0.5
:


interp_ymax	Scope: private	REAL

Description: ymax for interpolated region

Range		Default: +0.5
:


interp_ymin	Scope: private	REAL

Description: ymin for interpolated region

Range		Default: -0.5
:


interp_zmax	Scope: private	REAL

Description: zmax for interpolated region

Range		Default: +0.5
:


interp_zmin	Scope: private	REAL

Description: zmin for interpolated region

Range		Default: -0.5
:


interpolator	Scope: private	STRING

Description: The interpolator to use

Range		Default: Lagrange polynomial interpolation
	must be a registered interpolator


interpolator_options	Scope: private	STRING

Description: Options for the interpolator

Range		Default: order=2
	must be a valid options specification


power_t	Scope: private	INT

Description: Polynomial power of t coordinate

Range		Default: (none)
:


power_x	Scope: private	INT

Description: Polynomial power of x coordinate

Range		Default: (none)
:


power_y	Scope: private	INT

Description: Polynomial power of y coordinate

Range		Default: (none)
:


power_z	Scope: private	INT

Description: Polynomial power of z coordinate

Range		Default: (none)
:


prolongation	Scope: private	KEYWORD

Description: The prolongation operator to use

Range		Default: Lagrange
Lagrange
ENO
WENO
TVD
Lagrange_monotone

7 Interfaces

General

Implements:

carpetprolongatetest

Inherits:

grid

Grid Variables

7.0.1 PRIVATE GROUPS


Group Names	Variable Names	Details

scalar		compact	0
	u	description	Grid function
		dimensions	3
		distribution	DEFAULT
		group type	GF
		tags	ProlongationParameter=”CarpetProlongateTest::prolongation”
		timelevels	3
		variable type	REAL

scaled		compact	0
	uscaled	description	Scaled grid function
		dimensions	3
		distribution	DEFAULT
		group type	GF
		timelevels	3
		variable type	REAL

difference		compact	0
	u0	description	Error in grid function
	du	dimensions	3
		distribution	DEFAULT
		group type	GF
		timelevels	3
		variable type	REAL

interp_difference		compact	0
	interp_x	description	Error in interpolated grid array
	interp_y	dimensions	3
	interp_z	distribution	DEFAULT
	interp_u	group type	ARRAY
	interp_u0	size	INTERP_NX
		size	INTERP_NY
	interp_u0	size	INTERP_NZ
	interp_du	timelevels	1
		variable type	REAL

errornorm	errornorm	compact	0
		description	Norm of error in grid function
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		variable type	REAL

interp_errornorm	interp_errornorm	compact	0
		description	Norm of error in interpolated grid array
		dimensions	0
		distribution	CONSTANT
		group type	SCALAR
		timelevels	1
		variable type	REAL

8 Schedule

This section lists all the variables which are assigned storage by thorn Carpet/CarpetProlongateTest. 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:
scalar[3] scaled[3] difference[3]
interp_difference
errornorm interp_errornorm

Scheduled Functions

CCTK_INITIAL

carpetprolongatetest_init

set up initial data

	Language:	fortran
	Sync:	scalar
		scaled
	Type:	function
	Writes:	u(interior)
		uscaled(interior)

CCTK_EVOL

carpetprolongatetest_init

set up initial data

	Language:	fortran
	Reads:	grid::x
		grid::y
		grid::z
	Sync:	scalar
		scaled
	Type:	function
	Writes:	u(interior)
		uscaled(interior)

MoL_PostStep

carpetprolongatetest_diff

test data

	Language:	fortran
	Reads:	grid::x
		grid::y
		grid::z
		u
	Type:	function
	Writes:	u0(everywhere)
		du(everywhere)

CCTK_INITIAL

carpetprolongatetest_interpinit

set up interpolation

	Language:	fortran
	Options:	global-late
	Type:	function
	Writes:	interp_x(everywhere)
		interp_y(everywhere)
		interp_z(everywhere)

MoL_PostStep

carpetprolongatetest_interp

interpolate

	Language:	fortran
	Options:	global-late
	Reads:	interp_x
		interp_y
		interp_z
	Type:	function
	Writes:	interp_u(everywhere)

MoL_PostStep

carpetprolongatetest_interpdiff

test interpolated data

	After:	carpetprolongatetest_interp
	Language:	fortran
	Options:	global-late
	Reads:	interp_x
		interp_y
		interp_z
		interp_u
	Type:	function
	Writes:	interp_u0(everywhere)
		interp_du(everywhere)

MoL_PostStep

carpetprolongatetest_norminit

calculate error norm

	After:	carpetprolongatetest_diff
		carpetprolongatetest_interpdiff
	Language:	fortran
	Options:	global-late
	Type:	function
	Writes:	errornorm(everywhere)
		interp_errornorm(everywhere)

MoL_PostStep

carpetprolongatetest_normcalc

calculate error norm

	After:	carpetprolongatetest_norminit
	Language:	fortran
	Options:	global-late
		loop-local
	Reads:	errornorm
		interp_errornorm
		du
		interp_du
	Type:	function
	Writes:	errornorm(everywhere)
		interp_errornorm(everywhere)

MoL_PostStep

carpetprolongatetest_normreduce

calculate error norm

	After:	carpetprolongatetest_normcalc
	Language:	fortran
	Options:	global-late
	Reads:	errornorm
		interp_errornorm
	Type:	function
	Writes:	errornorm(everywhere)
		interp_errornorm(everywhere)