self_force_base – SelfForce1D

Module that defines a global self-force object that can be used to pass information between different equations using the save_globals/load_globals mechanism.

Uses

- kinds
Help
Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Used by

Help
Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Variables

Derived Types

self_force

Subroutines

set_force get_force set_accel set_daccel_dt set_d2accel_dt2 get_accel get_daccel_dt get_d2accel_dt2 output

Variables

Type	Visibility	Attributes		Name		Initial
type(self_force),	public		::	sf			The global self_force object that is available by use association.

Derived Types

type, public :: self_force

A global self-force type.

Components

Type	Visibility	Attributes		Name		Initial
real(kind=wp),	private,	dimension(4)	::	f	=	0.0_wp	A 4-vector containing the components of the self-force, $f_{\mu}$ .
real(kind=wp),	private,	dimension(4)	::	a	=	0.0_wp	A 4-vector containing the components of the 4-acceleration, $a^{\mu}$ .
real(kind=wp),	private,	dimension(4)	::	da_dt	=	0.0_wp	A 4-vector containing the components of the coordinate time derivatives of the 4-acceleration, $\dot{a}^{\mu}$ .
real(kind=wp),	private,	dimension(4)	::	d2a_dt2	=	0.0_wp	A 4-vector containing the components of the second coordinate time derivatives of the 4-acceleration, $\ddot{a}^{\mu}$ .
real(kind=wp),	private		::	udota	=	0.0_wp	A scalar containing the inner product of the 4-velocity and the self-force $u^{\mu}f_{\mu}$ .
integer(kind=ip),	private		::	ioo_id	=	-1	The file unit number used for output.

Type-Bound Procedures

procedure, public :: set_force	Routine to set the self-force values.
procedure, public :: get_force	Routine to get the self-force values.
procedure, public :: set_accel	Routine to set the values of the 4-acceleration.
procedure, public :: set_daccel_dt	Routine to set the values of the time derivative of the 4-acceleration.
procedure, public :: set_d2accel_dt2	Routine to set the values of the second time derivative of the 4-acceleration.
procedure, public :: get_accel	Routine to get the values of the 4-acceleration.
procedure, public :: get_daccel_dt	Routine to get the values of the time derivative of the 4-acceleration.
procedure, public :: get_d2accel_dt2	Routine to get the values of the second time derivative of the 4-acceleration.
procedure, public :: output	Routine to perform output.

Subroutines

public subroutine set_force(this, ft, fr, ftheta, fphi)

Routine that sets the self-force variable.

Arguments

Type	Intent		Name
class(self_force),	intent(inout)	::	this	The routine is called on this self_force object.
real(kind=wp),	intent(in)	::	ft	The $t$ -component of the self-force, $f_t$ .
real(kind=wp),	intent(in)	::	fr	The $r$ -component of the self-force, $f_r$ .
real(kind=wp),	intent(in)	::	ftheta	The $\theta$ -component of the self-force, $f_{\theta}$ .
real(kind=wp),	intent(in)	::	fphi	The $\phi$ -component of the self-force, $f_{\phi}$ .

public subroutine get_force(this, force)

Routine that gets the self-force variable.

Arguments

Type	Intent	Optional	Attributes		Name
class(self_force),	intent(inout)			::	this	The routine is called on this self_force object.
real(kind=wp),	intent(out),		dimension(4)	::	force	On return this 1d array contains the 4 components of the force, $(f_t,f_r,f_{\theta},f_{\phi})$ .

public subroutine set_accel(this, at, ar, atheta, aphi, udota)

Routine that sets the 4-acceleration and the inner product of the 4-velovity and the self-force.

Arguments

Type	Intent		Name
class(self_force),	intent(inout)	::	this	The routine is called on this self_force object.
real(kind=wp),	intent(in)	::	at	The $t$ -component of the 4-acceleration, $a^t$ .
real(kind=wp),	intent(in)	::	ar	The $r$ -component of the 4-acceleration, $a^r$ .
real(kind=wp),	intent(in)	::	atheta	The $\theta$ -component of the 4-acceleration, $a^{\theta}$ .
real(kind=wp),	intent(in)	::	aphi	The $\phi$ -component of the 4-acceleration, $a^{\phi}$ .
real(kind=wp),	intent(in)	::	udota	The inner product of the 4-velocity and the self-force, $u^{\mu}f_{\mu}$ .

public subroutine set_daccel_dt(this, dat_dt, dar_dt, datheta_dt, daphi_dt)

Routine that sets the coordinate time derivative of the 4-acceleration.

Arguments

Type	Intent		Name
class(self_force),	intent(inout)	::	this	The routine is called on this self_force object.
real(kind=wp),	intent(in)	::	dat_dt	The time derivative of the $t$ -component of the 4-acceleration, $\dot{a}^t$ .
real(kind=wp),	intent(in)	::	dar_dt	The time derivative of the $r$ -component of the 4-acceleration, $\dot{a}^r$ .
real(kind=wp),	intent(in)	::	datheta_dt	The time derivative of the $\theta$ -component of the 4-acceleration, $\dot{a}^{\theta}$ .
real(kind=wp),	intent(in)	::	daphi_dt	The time derivative of the $\phi$ -component of the 4-acceleration, $\dot{a}^{\phi}$ .

public subroutine set_d2accel_dt2(this, d2at_dt2, d2ar_dt2, d2atheta_dt2, d2aphi_dt2)

Routine that sets the second coordinate time derivative of the 4-acceleration.

Arguments

Type	Intent		Name
class(self_force),	intent(inout)	::	this	The routine is called on this self_force object.
real(kind=wp),	intent(in)	::	d2at_dt2	The second time derivative of the $t$ -component of the 4-acceleration, $\ddot{a}^t$ .
real(kind=wp),	intent(in)	::	d2ar_dt2	The second time derivative of the $r$ -component of the 4-acceleration, $\ddot{a}^r$ .
real(kind=wp),	intent(in)	::	d2atheta_dt2	The second time derivative of the $\theta$ -component of the 4-acceleration, $\ddot{a}^{\theta}$ .
real(kind=wp),	intent(in)	::	d2aphi_dt2	The second time derivative of the $\phi$ -component of the 4-acceleration, $\ddot{a}^{\phi}$ .

public subroutine get_accel(this, accel)

Routine that gets the 4-acceleration.

Arguments

Type	Intent	Optional	Attributes		Name
class(self_force),	intent(inout)			::	this	The routine is called on this self_force object.
real(kind=wp),	intent(out),		dimension(4)	::	accel	On return this 1d array contains the 4 components of the 4-acceleration, $(a^t,a^r,a^{\theta},a^{\phi})$ .

public subroutine get_daccel_dt(this, daccel_dt)

Routine that gets the coordinate time derivative of the 4-acceleration.

Arguments

Type	Intent	Optional	Attributes		Name
class(self_force),	intent(inout)			::	this	The routine is called on this self_force object.
real(kind=wp),	intent(out),		dimension(4)	::	daccel_dt	On return this 1d array contains the 4 components of the time derivative of the 4-acceleration, $(\dot{a}^t,\dot{a}^r,\dot{a}^{\theta},\dot{a}^{\phi})$ .

public subroutine get_d2accel_dt2(this, d2accel_dt2)

Routine that gets the second coordinate time derivative of the 4-acceleration.

Arguments

Type	Intent	Optional	Attributes		Name
class(self_force),	intent(inout)			::	this	The routine is called on this self_force object.
real(kind=wp),	intent(out),		dimension(4)	::	d2accel_dt2	On return this 1d array contains the 4 components of the second time derivative of the 4-acceleration, $(\ddot{a}^t,\ddot{a}^r,\ddot{a}^{\theta},\ddot{a}^{\phi})$ .

public subroutine output(this)

Routine that performs output of the self-force.

Arguments

Type	Intent	Optional	Attributes		Name
class(self_force),	intent(inout)			::	this	The routine is called on this self_force object.

self_force_base Module

Contents

Variables

Derived Types

Subroutines

Uses

Used by

Contents

Variables

Derived Types

Subroutines

Variables

Derived Types

type, public :: self_force

Components

Type-Bound Procedures

Subroutines

public subroutine set_force(this, ft, fr, ftheta, fphi)

Arguments

public subroutine get_force(this, force)

Arguments

public subroutine set_accel(this, at, ar, atheta, aphi, udota)

Arguments

public subroutine set_daccel_dt(this, dat_dt, dar_dt, datheta_dt, daphi_dt)

Arguments

public subroutine set_d2accel_dt2(this, d2at_dt2, d2ar_dt2, d2atheta_dt2, d2aphi_dt2)

Arguments

public subroutine get_accel(this, accel)

Arguments

public subroutine get_daccel_dt(this, daccel_dt)

Arguments

public subroutine get_d2accel_dt2(this, d2accel_dt2)

Arguments

public subroutine output(this)

Arguments