mod_grid Module

Module containing all grid-related things. Contains subroutines to create the base grid, Gaussian grid and scale factors. An integral of $f(x)$ in $[a, b]$ can be approximated with $$f(x) \approx 0.5(b-a)\sum_{i=1}^n\bigl[w_i f\bigl(0.5(b-a)x_i + 0.5(a+b)\bigr)\bigr] !!$$ where $w_i$ and $x_i$ are the weights and nodes of the Gaussian quadrature. The Gaussian grid is hence set up in every interval $[a, b]$ across the nodes $j$ as $$x_i(j) = 0.5 * dx * w_i(j) + 0.5(a + b)$$

Contents

Interfaces

dx_func_i

Derived Types

grid_t

Functions

new_grid get_eps get_deps is_valid_custom_base_grid get_updated_number_of_gridpoints is_valid_dx

Subroutines

initialise set_custom_grid set_spacing_function set_base_grid set_gaussian_grid set_ef_grid generate_grid delete log_msg_by_gridpoint_change

Interfaces

interface

• private function dx_func_i(x)

  Arguments

  Type	Intent	Optional	Attributes		Name
  real(kind=dp),	intent(in)			::	x

  Return Value real(kind=dp)

Derived Types

type, public :: grid_t

Components

Type	Visibility	Attributes		Name		Initial
real(kind=dp),	public,	allocatable	::	base_grid(:)
real(kind=dp),	public,	allocatable	::	gaussian_grid(:)
real(kind=dp),	public,	allocatable	::	ef_grid(:)
procedure(dx_func_i),	private,	pointer, nopass	::	dx_func	=>	null()
type(settings_t),	private,	pointer	::	settings
logical,	private		::	is_initialised
logical,	private		::	uses_custom_base_grid
logical,	private		::	uses_custom_dx

Type-Bound Procedures

procedure, private :: set_base_grid
procedure, private :: set_gaussian_grid
procedure, private :: set_ef_grid
procedure, private :: generate_grid
procedure, public :: initialise
procedure, public :: set_custom_grid
procedure, public :: set_spacing_function
procedure, public :: get_eps
procedure, public :: get_deps
procedure, public :: delete

Functions

public function new_grid(settings) result(grid)

Arguments

Type	Intent	Optional	Attributes		Name
type(settings_t),	intent(in),		target	::	settings

Return Value type(grid_t)

private impure elemental function get_eps(this, x)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(in)			::	this
real(kind=dp),	intent(in)			::	x

Return Value real(kind=dp)

private impure elemental function get_deps(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(in)			::	this

Return Value real(kind=dp)

private function is_valid_custom_base_grid(settings, custom_grid)

Arguments

Type	Intent	Optional	Attributes		Name
type(settings_t),	intent(inout)			::	settings
real(kind=dp),	intent(in)			::	custom_grid(:)

Return Value logical

private function get_updated_number_of_gridpoints(settings, dx_func) result(updated_pts)

Arguments

Type	Intent	Optional	Attributes		Name
type(settings_t),	intent(in)			::	settings
procedure(dx_func_i)				::	dx_func

Return Value integer

private function is_valid_dx(dx)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	dx

Return Value logical

Subroutines

private subroutine initialise(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this

private subroutine set_custom_grid(this, custom)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this
real(kind=dp),	intent(in)			::	custom(:)

private subroutine set_spacing_function(this, dx_func)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this
procedure(dx_func_i)				::	dx_func

private subroutine set_base_grid(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this

private pure subroutine set_gaussian_grid(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this

private pure subroutine set_ef_grid(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this

private subroutine generate_grid(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this

private pure subroutine delete(this)

Arguments

Type	Intent	Optional	Attributes		Name
class(grid_t),	intent(inout)			::	this

private subroutine log_msg_by_gridpoint_change(old_pts, new_pts)

Arguments

Type	Intent	Optional	Attributes		Name
integer,	intent(in)			::	old_pts
integer,	intent(in)			::	new_pts