mod_matrix_manager Module

Uses

Used by

- Descendants:
- smod_conduction_matrix
- smod_flow_matrix
- smod_hall_matrix
- smod_heatloss_matrix
- smod_regular_matrix
- smod_resistive_matrix
- smod_viscosity_matrix

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=dp),	public		::	h_quad(4)			quadratic basis functions
real(kind=dp),	public		::	dh_quad(4)			derivative of quadratic basis functions
real(kind=dp),	public		::	h_cubic(4)			cubic basis functions
real(kind=dp),	public		::	dh_cubic(4)			derivative of cubic basis functions

Interfaces

interface

public module subroutine add_bmatrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_regular_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_flow_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_resistive_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_heatloss_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_conduction_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_viscosity_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_hall_matrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

interface

public module subroutine add_hall_bmatrix_terms(x_gauss, weight, quadblock, settings, grid, background, physics)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=dp),	intent(in)			::	x_gauss
real(kind=dp),	intent(in)			::	weight
complex(kind=dp),	intent(inout)			::	quadblock(:,:)
type(settings_t),	intent(in)			::	settings
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics

Subroutines

public subroutine build_matrices(matrix_B, matrix_A, settings, grid, background, physics)

The quadblock is shifted along the main (tri)diagonal. We add dim_subblock instead of dim_quadblock to the indices, due to overlapping of the bottom-right part of the quadblock with the top-left part of the next grid interval.

Arguments

Type	Intent	Optional	Attributes		Name
type(matrix_t),	intent(inout)			::	matrix_B	the B-matrix
type(matrix_t),	intent(inout)			::	matrix_A	the A-matrix
type(settings_t),	intent(in)			::	settings	the settings object
type(grid_t),	intent(in)			::	grid
type(background_t),	intent(in)			::	background
type(physics_t),	intent(in)			::	physics