self_force – SelfForce1D

type, public :: self_force

A global self-force type.

Variables

f a da_dt d2a_dt2 udota ioo_id

Type-Bound Procedures

set_force get_force set_accel set_daccel_dt set_d2accel_dt2 get_accel get_daccel_dt get_d2accel_dt2 output

Source Code

self_force

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.

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}$ .

procedure, public :: get_force

Routine to get the self-force values.

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})$ .

procedure, public :: set_accel

Routine to set the values of the 4-acceleration.

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}$ .

procedure, public :: set_daccel_dt

Routine to set the values of the time derivative of the 4-acceleration.

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}$ .

procedure, public :: set_d2accel_dt2

Routine to set the values of the second time derivative of the 4-acceleration.

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}$ .

procedure, public :: get_accel

Routine to get the values of the 4-acceleration.

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})$ .

procedure, public :: get_daccel_dt

Routine to get the values of the time derivative of the 4-acceleration.

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})$ .

procedure, public :: get_d2accel_dt2

Routine to get the values of the second time derivative of the 4-acceleration.

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})$ .

procedure, public :: output

Routine to perform output.

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.

Source Code

  type :: self_force
  !! A global self-force type.
    real(wp), dimension(4), private :: f = 0.0_wp
    !! A 4-vector containing the components of the self-force, \(f_{\mu}\).
    real(wp), dimension(4), private :: a = 0.0_wp
    !! A 4-vector containing the components of the 4-acceleration, \(a^{\mu}\).
    real(wp), dimension(4), private :: da_dt = 0.0_wp
    !! A 4-vector containing the components of the coordinate time derivatives
    !! of the 4-acceleration, \(\dot{a}^{\mu}\).
    real(wp), dimension(4), private :: 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(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(ip), private :: ioo_id = -1
    !! The file unit number used for output.
  contains
    procedure :: set_force
    !! Routine to set the self-force values.
    procedure :: get_force
    !! Routine to get the self-force values.
    procedure :: set_accel
    !! Routine to set the values of the 4-acceleration.
    procedure :: set_daccel_dt
    !! Routine to set the values of the time derivative of the 4-acceleration.
    procedure :: set_d2accel_dt2
    !! Routine to set the values of the second time derivative of the
    !! 4-acceleration.
    procedure :: get_accel
    !! Routine to get the values of the 4-acceleration.
    procedure :: get_daccel_dt
    !! Routine to get the values of the time derivative of the 4-acceleration.
    procedure :: get_d2accel_dt2
    !! Routine to get the values of the second time derivative of the
    !! 4-acceleration.
    procedure :: output
    !! Routine to perform output.
  end type self_force

self_force Derived Type

Contents

Variables

Type-Bound Procedures

Source Code

type, public :: self_force

Contents

Variables

Type-Bound Procedures

Source Code

Components

Type-Bound Procedures

procedure, public :: set_force

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

Arguments

procedure, public :: get_force

public subroutine get_force(this, force)

Arguments

procedure, public :: set_accel

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

Arguments

procedure, public :: set_daccel_dt

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

Arguments

procedure, public :: set_d2accel_dt2

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

Arguments

procedure, public :: get_accel

public subroutine get_accel(this, accel)

Arguments

procedure, public :: get_daccel_dt

public subroutine get_daccel_dt(this, daccel_dt)

Arguments

procedure, public :: get_d2accel_dt2

public subroutine get_d2accel_dt2(this, d2accel_dt2)

Arguments

procedure, public :: output

public subroutine output(this)

Arguments

Source Code