scalar_schw_eff Module

  1. module_scalar_schw_eff_source.f90
  2. scalar_schw_eff

Module that defines an effective source class for a generic effective source for a scalar charge in orbit around a Schwarzschild black hole.

The implementation in submodule_scalar_schw_eff_source_implementation.f90 is an interface to a C++ class provided by Barry Wardell.


Uses

  • module~~scalar_schw_eff~~UsesGraph module~scalar_schw_eff scalar_schw_eff module~kinds kinds module~scalar_schw_eff->module~kinds module~effective_source effective_source module~scalar_schw_eff->module~effective_source iso_c_binding iso_c_binding module~scalar_schw_eff->iso_c_binding module~effective_source->module~kinds module~effective_source->iso_c_binding module~grid_function grid_function module~effective_source->module~grid_function module~world_tube world_tube module~effective_source->module~world_tube module~grid_function->module~kinds module~element element module~grid_function->module~element module~world_tube->module~kinds module~world_tube->module~grid_function module~element->module~kinds

Used by

  • module~~scalar_schw_eff~~UsedByGraph module~scalar_schw_eff scalar_schw_eff module~scalar_schw scalar_schw module~scalar_schw->module~scalar_schw_eff module~scalar_schw_eff_implementation scalar_schw_eff_implementation module~scalar_schw_eff_implementation->module~scalar_schw_eff program~test test program~test->module~scalar_schw_eff program~test->module~scalar_schw module~self_force_observer self_force_observer program~test->module~self_force_observer proc~scal_schw_eff_init scal_schw_eff_init proc~scal_schw_eff_init->module~scalar_schw module~scalar_schw_implementation scalar_schw_implementation module~scalar_schw_implementation->module~scalar_schw module~self_force_observer->module~scalar_schw proc~scal_schw_save_globals_2 scal_schw_save_globals_2 proc~scal_schw_save_globals_2->module~self_force_observer module~self_force_observer_implementation self_force_observer_implementation module~self_force_observer_implementation->module~self_force_observer

Contents

Interfaces

Derived Types


Interfaces

interface

  • public module subroutine scal_schw_eff_init(this, nmodes, nvars, l, m, mass)

    Interface of the constructor compatible with init.

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    On return, the constructed object.
    integer(kind=ip), intent(in) :: nmodes

    The number of modes.
    integer(kind=ip), intent(in) :: nvars

    The number of variables.

    integer(kind=c_int), intent(in), dimension(nmodes):: l

    A 1d array of size nmodes containing the -values of the modes.
    integer(kind=c_int), intent(in), dimension(nmodes):: m

    A 1d array of size nmodes containing the m-values of the modes.
    real(kind=wp), intent(in) :: mass

    The mass of the black hole.

interface

  • public module subroutine scal_schw_eff_set_time_window(this, tfac, dtfac_dt, d2tfac_dt2, do_smooth_after_lmax)

    Interface to the set_time_window routine compatible with set_time_window.

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    The routine is called on this object.
    real(kind=wp), intent(in) :: tfac

    The current value for the time window.
    real(kind=wp), intent(in) :: dtfac_dt

    The current value for the time derivative of the time window.
    real(kind=wp), intent(in) :: d2tfac_dt2

    The current value for the second time derivative of the time window.
    integer(kind=ip), intent(in) :: do_smooth_after_lmax

    If do_smooth_after_lmax use tfac=1, dtfac_dt=0, d2tfac_dt2=0. This allows for smooth turn on of the effective source for the modes with do_smooth_after_lmax, while the rest gets turned on instantaneously (e.g. when external initial data is available).

interface

  • public module subroutine scal_schw_eff_set_particle_pos(this, r, phi, ur, En, Lz, ar, aphi, dardt, daphidt, d2ardt2, d2aphidt2)

    Interface of the routine that sets the state of the particle.

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    The routine is called on this object.
    real(kind=wp), intent(in) :: r

    The radial coordinate of the particle (in Schwarzschild coordinates), .
    real(kind=wp), intent(in) :: phi

    The azimuthal angle of the particle, .
    real(kind=wp), intent(in) :: ur

    The radial component of the 4-velocity, .
    real(kind=wp), intent(in) :: En

    The energy per unit mass of the particle, .
    real(kind=wp), intent(in) :: Lz

    The angular momentum per unit mass of the particle .
    real(kind=wp), intent(in) :: ar

    The radial component of the 4-acceleration, .
    real(kind=wp), intent(in) :: aphi

    The -component of the 4-acceleration, .
    real(kind=wp), intent(in) :: dardt

    The derivative of with respect to Schwarzschild coordinate time.
    real(kind=wp), intent(in) :: daphidt

    The derivative of with respect to Schwarzschild coordinate time.
    real(kind=wp), intent(in) :: d2ardt2

    The second derivative of with respect to Schwarzschild coordinate time.
    real(kind=wp), intent(in) :: d2aphidt2

    The second derivative of with respect to Schwarzschild coordinate time.

interface

  • public module subroutine scal_schw_eff_evaluate_source(this, r, wt)

    Interface of the routine that evaluates the effective source.

    Read more…

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    The routine is called on this object.
    type(rgf), intent(in) :: r

    A real values grid function that contain the radial coordinate (in Schwarzschild coordinates).
    type(wtube), intent(in) :: wt

    A world-tube object that ensures that the effective source is only non-zero inside the world-tube.

interface

  • public module subroutine scal_schw_eff_get_singular(this, r, mode, psi)

    Interface of the routine that evaluates the singular field for a given mode and at a given radial coordinate.

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    The routine is called on this object.
    real(kind=wp), intent(in) :: r

    The radial coordinate (Schwarzschild coordinates).
    integer(kind=ip), intent(in) :: mode

    The mode.
    complex(kind=wp), intent(out), dimension(:):: psi

    A 1d-array of size nvars of complex values that on return contains the singular field for all variables.

interface

  • public module subroutine scal_schw_eff_get_dsingular_dt(this, r, mode, dpsidt)

    Interface of the routine that evaluates the time derivative of the singular field for a given mode and at a given radial coordinate.

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    The routine is called on this object.
    real(kind=wp), intent(in) :: r

    The radial coordinate (Schwarzschild coordinates).
    integer(kind=ip), intent(in) :: mode

    The mode.
    complex(kind=wp), intent(out), dimension(:):: dpsidt

    A 1d-array of size nvars of complex values that on return contains the time derivative of the singular field for all variables.

interface

  • public module subroutine scal_schw_eff_get_dsingular_dr(this, r, mode, dpsidr)

    Interface of the routine that evaluates the radial derivative of the singular field for a given mode and at a given radial coordinate.

    Arguments

    Type IntentOptional AttributesName
    class(scal_schw_eff), intent(inout) :: this

    The routine is called on this object.
    real(kind=wp), intent(in) :: r

    The radial coordinate (Schwarzschild coordinates).
    integer(kind=ip), intent(in) :: mode

    The mode.
    complex(kind=wp), intent(out), dimension(:):: dpsidr

    A 1d-array of size nvars of complex values that on return contains the radial derivative of the singular field for all variables.

Derived Types

type, public, extends(eff_source) :: scal_schw_eff

A class that interfaces with a C++ class that provides an effective source for a scalar charge on a generic accelerated orbit around a Schwarzschild black hole.

Components

TypeVisibility AttributesNameInitial
integer(kind=ip), public :: nmodes

The number of modes an effective source is provided for.
integer(kind=ip), public :: nvars

The number of variables an effective source is provided for.
type(cgf), public, dimension(:,:), allocatable:: source

A 2d-array of complex grid functions. When allocated the size is (nvars:nmodes).
real(kind=wp), public, dimension(:,:), allocatable:: sre

A 2d-array of reals to hold the real part of the effective source for all nodes in a DG element and for all modes. When allocated the size is (order+1:nmodes).
real(kind=wp), public, dimension(:,:), allocatable:: sim

A 2d-array of reals to hold the complex part of the effective source for all nodes in a DG element and for all modes. When allocated the size is (order+1:nmodes).

Type-Bound Procedures

procedure, public :: init => scal_schw_eff_init

The init routine is provided by scal_schw_eff_init.

procedure, public :: set_time_window => scal_schw_eff_set_time_window

The set_particle_pos routine is provided by scal_schw_eff_set_particle_pos.

procedure, public :: set_particle_pos => scal_schw_eff_set_particle_pos

The set_particle_pos routine is provided by scal_schw_eff_set_particle_pos.

procedure, public :: evaluate_source => scal_schw_eff_evaluate_source

The evaluate_source routine is provided by scal_schw_eff_evaluate_source.

procedure, public :: get_singular => scal_schw_eff_get_singular

The get_singular routine is provided by scal_schw_eff_get_singular.

procedure, public :: get_dsingular_dt => scal_schw_eff_get_dsingular_dt

The get_dsingular_dt routine is provided by scal_schw_eff_get_dsingular_dt.

procedure, public :: get_dsingular_dr => scal_schw_eff_get_dsingular_dr

The get_dsingular_dr routine is provided by scal_schw_eff_get_dsingular_dr.