Data Types
type	mg_lvl_t
	Lists of blocks per refinement level. More...

type	mg_box_t
	Box data structure. More...

type	mg_buf_t
	Buffer type (one is used for each pair of communicating processes) More...

type	mg_comm_t

type	mg_bc_t

type	mg_timer_t

type	mg_t

interface	mg_subr_bc
	To fill ghost cells near physical boundaries. More...

Functions/Subroutines
elemental logical function, public	mg_has_children (box)
	Return .true. if a box has children. More...

integer function, public	mg_ix_to_ichild (ix)
	Compute the child index for a box with spatial index ix. With child index we mean the index in the children(:) array of its parent. More...

pure integer function, dimension(1), public	mg_get_child_offset (mg, id)
	Get the offset of a box with respect to its parent (e.g. in 2d, there can be a child at offset 0,0, one at n_cell/2,0, one at 0,n_cell/2 and one at n_cell/2, n_cell/2) More...

pure integer function, public	mg_highest_uniform_lvl (mg)

integer(i8) function, public	mg_number_of_unknowns (mg)
	Determine total number of unknowns (on leaves) More...

subroutine, public	mg_get_face_coords (box, nb, nc, x)
	Get coordinates at the face of a box. More...

integer function, public	mg_add_timer (mg, name)

subroutine, public	mg_timer_start (timer)

subroutine, public	mg_timer_end (timer)

subroutine, public	mg_timers_show (mg)

subroutine, public	mg_deallocate_storage (mg)
	Deallocate all allocatable arrays. More...

subroutine, public	mg_allocate_storage (mg)
	Allocate communication buffers and local boxes for a tree that has already been created. More...

subroutine, public	diffusion_solve (mg, dt, diffusion_coeff, order, max_res)
	Solve a diffusion equation implicitly, assuming a constant diffusion coefficient. The solution at time t should be stored in mg_iphi, which is on output replaced by the solution at time t+dt. More...

subroutine, public	diffusion_solve_vcoeff (mg, dt, order, max_res)
	Solve a diffusion equation implicitly, assuming a variable diffusion coefficient which has been stored in mg_iveps (also on coarse grids). The solution at time t should be stored in mg_iphi, which is on output replaced by the solution at time t+dt. More...

subroutine, public	diffusion_solve_acoeff (mg, dt, order, max_res)
	Solve a diffusion equation implicitly, assuming anisotropic diffusion coefficient which has been stored in mg_iveps1, mg_iveps2, mg_iveps3 (also on coarse grids). The solution at time t should be stored in mg_iphi, which is on output replaced by the solution at time t+dt. More...

subroutine, public	laplacian_set_methods (mg)

subroutine, public	vhelmholtz_set_methods (mg)

subroutine, public	vhelmholtz_set_lambda (lambda)

subroutine, public	mg_build_rectangle (mg, domain_size, box_size, dx, r_min, periodic, n_finer)

subroutine, public	mg_set_neighbors_lvl (mg, lvl)

subroutine, public	mg_set_next_level_ids (mg, lvl)

subroutine, public	mg_set_leaves_parents (boxes, level)
	Create a list of leaves and a list of parents for a level. More...

subroutine, public	mg_set_refinement_boundaries (boxes, level)
	Create a list of refinement boundaries (from the coarse side) More...

subroutine, public	mg_add_children (mg, id)

subroutine, public	mg_load_balance_simple (mg)
	Load balance all boxes in the multigrid tree, by simply distributing the load per grid level. This method will only work well for uniform grids. More...

subroutine, public	mg_load_balance (mg)
	Load balance all boxes in the multigrid tree. Compared to mg_load_balance_simple, this method does a better job of setting the ranks of parent boxes. More...

subroutine, public	mg_load_balance_parents (mg)
	Load balance the parents (non-leafs). Assign them to the rank that has most children. More...

subroutine, public	vlaplacian_set_methods (mg)

subroutine, public	mg_comm_init (mg, comm)
	Prolong from a parent to a child with index offset dix. This method could sometimes work better than the default prolongation, which does not take the variation in epsilon into account. More...

subroutine, public	sort_and_transfer_buffers (mg, dsize)

subroutine	sort_sendbuf (gc, dsize)
	Sort send buffers according to the idbuf array. More...

subroutine, public	mg_ghost_cell_buffer_size (mg, n_send, n_recv, dsize)
	Specify minimum buffer size (per process) for communication. More...

subroutine, public	mg_phi_bc_store (mg)
	Store boundary conditions for the solution variable, this can speed up calculations if the same boundary conditions are re-used. More...

subroutine, public	mg_fill_ghost_cells (mg, iv)
	Fill ghost cells at all grid levels. More...

subroutine, public	mg_fill_ghost_cells_lvl (mg, lvl, iv)
	Fill ghost cells at a grid level. More...

subroutine, public	mg_prolong_buffer_size (mg, n_send, n_recv, dsize)
	Specify minimum buffer size (per process) for communication. More...

subroutine, public	mg_prolong (mg, lvl, iv, iv_to, method, add)
	Prolong variable iv from lvl to variable iv_to at lvl+1. More...

subroutine, public	mg_prolong_sparse (mg, p_id, dix, nc, iv, fine)
	Prolong from a parent to a child with index offset dix. More...

subroutine, public	helmholtz_set_methods (mg)

subroutine, public	helmholtz_set_lambda (lambda)

subroutine, public	mg_set_methods (mg)

subroutine, public	mg_fas_fmg (mg, have_guess, max_res)
	Perform FAS-FMG cycle (full approximation scheme, full multigrid). More...

subroutine, public	mg_fas_vcycle (mg, highest_lvl, max_res, standalone)
	Perform FAS V-cycle (full approximation scheme). More...

subroutine	subtract_mean (mg, iv, include_ghostcells)

real(dp) function	get_sum (mg, iv)

subroutine, public	mg_apply_op (mg, i_out, op)
	Apply operator to the tree and store in variable i_out. More...

subroutine, public	mg_restrict_buffer_size (mg, n_send, n_recv, dsize)
	Specify minimum buffer size (per process) for communication. More...

subroutine, public	mg_restrict (mg, iv)
	Restrict all levels. More...

subroutine, public	mg_restrict_lvl (mg, iv, lvl)
	Restrict from lvl to lvl-1. More...

subroutine, public	ahelmholtz_set_methods (mg)

subroutine, public	ahelmholtz_set_lambda (lambda)

Variables
integer, parameter, public	dp = kind(0.0d0)
	Type of reals. More...

integer, parameter, public	i8 = selected_int_kind(18)
	Type for 64-bit integers. More...

integer, parameter, public	mg_laplacian = 1
	Indicates a standard Laplacian. More...

integer, parameter, public	mg_vlaplacian = 2
	Indicates a variable-coefficient Laplacian. More...

integer, parameter, public	mg_helmholtz = 3
	Indicates a constant-coefficient Helmholtz equation. More...

integer, parameter, public	mg_vhelmholtz = 4
	Indicates a variable-coefficient Helmholtz equation. More...

integer, parameter, public	mg_ahelmholtz = 5
	Indicates a anisotropic-coefficient Helmholtz equation. More...

integer, parameter, public	mg_cartesian = 1
	Cartesian coordinate system. More...

integer, parameter, public	mg_cylindrical = 2
	Cylindrical coordinate system. More...

integer, parameter, public	mg_spherical = 3
	Spherical coordinate system. More...

integer, parameter, public	mg_smoother_gs = 1

integer, parameter, public	mg_smoother_gsrb = 2

integer, parameter, public	mg_smoother_jacobi = 3

integer, parameter, public	mg_ndim = 1
	Problem dimension. More...

integer, parameter, public	mg_num_vars = 4
	Number of predefined multigrid variables. More...

integer, parameter, public	mg_max_num_vars = 10
	Maximum number of variables. More...

integer, parameter, public	mg_iphi = 1
	Index of solution. More...

integer, parameter, public	mg_irhs = 2
	Index of right-hand side. More...

integer, parameter, public	mg_iold = 3
	Index of previous solution (used for correction) More...

integer, parameter, public	mg_ires = 4
	Index of residual. More...

integer, parameter, public	mg_iveps = 5
	Index of the variable coefficient (at cell centers) More...

integer, parameter, public	mg_iveps1 = 5
	Indexes of anisotropic variable coefficients. More...

integer, parameter, public	mg_iveps2 = 6

integer, parameter, public	mg_iveps3 = 7

integer, parameter, public	mg_lvl_lo = -20
	Minimum allowed grid level. More...

integer, parameter, public	mg_lvl_hi = 20
	Maximum allowed grid level. More...

integer, parameter, public	mg_bc_dirichlet = -10
	Value to indicate a Dirichlet boundary condition. More...

integer, parameter, public	mg_bc_neumann = -11
	Value to indicate a Neumann boundary condition. More...

integer, parameter, public	mg_bc_continuous = -12
	Value to indicate a continuous boundary condition. More...

integer, parameter, public	mg_no_box = 0
	Special value that indicates there is no box. More...

integer, parameter, public	mg_physical_boundary = -1
	Special value that indicates there is a physical boundary. More...

integer, parameter, public	mg_max_timers = 20
	Maximum number of timers to use. More...

integer, parameter, public	mg_num_children = 2

integer, dimension(1, 2), parameter, public	mg_child_dix = reshape([0,1], [1,2])

integer, dimension(2, 1), parameter, public	mg_child_rev = reshape([2,1], [2,1])

integer, dimension(1, 2), parameter, public	mg_child_adj_nb = reshape([1,2], [1,2])

logical, dimension(1, 2), parameter, public	mg_child_low = reshape([.true., .false.], [1, 2])

integer, parameter, public	mg_num_neighbors = 2

integer, parameter, public	mg_neighb_lowx = 1

integer, parameter, public	mg_neighb_highx = 2

integer, dimension(1, 2), parameter, public	mg_neighb_dix = reshape([-1,1], [1,2])

logical, dimension(2), parameter, public	mg_neighb_low = [.true., .false.]

integer, dimension(2), parameter, public	mg_neighb_high_pm = [-1, 1]

integer, dimension(2), parameter, public	mg_neighb_rev = [2, 1]

integer, dimension(2), parameter, public	mg_neighb_dim = [1, 1]

real(dp), public, protected	vhelmholtz_lambda = 0.0_dp
	Module for implicitly solving diffusion equations. More...

real(dp), public, protected	helmholtz_lambda = 0.0_dp
	Module for load balancing a tree (that already has been constructed). The load balancing determines which ranks (MPI processes) allocated physical storage for boxes. The tree structure itself is present on all processes. More...

real(dp), public, protected	ahelmholtz_lambda = 0.0_dp
	Module which contains multigrid procedures for a Helmholtz operator of the form: div(D grad(phi)) - lambda*phi = f, where D has a smooth spatial variation and a component in each spatial direction. More...

Function/Subroutine Documentation

◆ ahelmholtz_set_lambda()

subroutine, public m_octree_mg_1d::ahelmholtz_set_lambda ( real(dp), intent(in) lambda )

Definition at line 3613 of file m_octree_mg_1d.t.

◆ ahelmholtz_set_methods()

subroutine, public m_octree_mg_1d::ahelmholtz_set_methods ( type(mg_t), intent(inout) mg )

Definition at line 3581 of file m_octree_mg_1d.t.

◆ diffusion_solve()

subroutine, public m_octree_mg_1d::diffusion_solve	(	type(mg_t), intent(inout)	mg,
		real(dp), intent(in)	dt,
		real(dp), intent(in)	diffusion_coeff,
		integer, intent(in)	order,
		real(dp), intent(in)	max_res
	)

Solve a diffusion equation implicitly, assuming a constant diffusion coefficient. The solution at time t should be stored in mg_iphi, which is on output replaced by the solution at time t+dt.

Definition at line 710 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ diffusion_solve_acoeff()

subroutine, public m_octree_mg_1d::diffusion_solve_acoeff	(	type(mg_t), intent(inout)	mg,
		real(dp), intent(in)	dt,
		integer, intent(in)	order,
		real(dp), intent(in)	max_res
	)

Solve a diffusion equation implicitly, assuming anisotropic diffusion coefficient which has been stored in mg_iveps1, mg_iveps2, mg_iveps3 (also on coarse grids). The solution at time t should be stored in mg_iphi, which is on output replaced by the solution at time t+dt.

Definition at line 806 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ diffusion_solve_vcoeff()

subroutine, public m_octree_mg_1d::diffusion_solve_vcoeff	(	type(mg_t), intent(inout)	mg,
		real(dp), intent(in)	dt,
		integer, intent(in)	order,
		real(dp), intent(in)	max_res
	)

Solve a diffusion equation implicitly, assuming a variable diffusion coefficient which has been stored in mg_iveps (also on coarse grids). The solution at time t should be stored in mg_iphi, which is on output replaced by the solution at time t+dt.

Definition at line 757 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ get_sum()

real(dp) function m_octree_mg_1d::get_sum	(	type(mg_t), intent(in)	mg,
		integer, intent(in)	iv
	)

Definition at line 3031 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ helmholtz_set_lambda()

subroutine, public m_octree_mg_1d::helmholtz_set_lambda ( real(dp), intent(in) lambda )

Definition at line 2719 of file m_octree_mg_1d.t.

◆ helmholtz_set_methods()

subroutine, public m_octree_mg_1d::helmholtz_set_methods ( type(mg_t), intent(inout) mg )

Definition at line 2698 of file m_octree_mg_1d.t.

◆ laplacian_set_methods()

subroutine, public m_octree_mg_1d::laplacian_set_methods ( type(mg_t), intent(inout) mg )

Definition at line 868 of file m_octree_mg_1d.t.

◆ mg_add_children()

subroutine, public m_octree_mg_1d::mg_add_children	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	id
	)

Parameters

[in] id Id of box that gets children

Definition at line 1411 of file m_octree_mg_1d.t.

◆ mg_add_timer()

integer function, public m_octree_mg_1d::mg_add_timer	(	type(mg_t), intent(inout)	mg,
		character(len=*), intent(in)	name
	)

Definition at line 575 of file m_octree_mg_1d.t.

◆ mg_allocate_storage()

subroutine, public m_octree_mg_1d::mg_allocate_storage ( type(mg_t), intent(inout) mg )

Allocate communication buffers and local boxes for a tree that has already been created.

Definition at line 657 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_apply_op()

subroutine, public m_octree_mg_1d::mg_apply_op	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	i_out,
		procedure(mg_box_op), optional	op
	)

Apply operator to the tree and store in variable i_out.

Definition at line 3192 of file m_octree_mg_1d.t.

◆ mg_build_rectangle()

subroutine, public m_octree_mg_1d::mg_build_rectangle	(	type(mg_t), intent(inout)	mg,
		integer, dimension(1), intent(in)	domain_size,
		integer, intent(in)	box_size,
		real(dp), dimension(1), intent(in)	dx,
		real(dp), dimension(1), intent(in)	r_min,
		logical, dimension(1), intent(in)	periodic,
		integer, intent(in)	n_finer
	)

Definition at line 1102 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_comm_init()

subroutine, public m_octree_mg_1d::mg_comm_init	(	type(mg_t), intent(inout)	mg,
		integer, intent(in), optional	comm
	)

Prolong from a parent to a child with index offset dix. This method could sometimes work better than the default prolongation, which does not take the variation in epsilon into account.

Initialize MPI if needed, and store MPI information

Parameters

[in] comm MPI communicator (default: MPI_COMM_WORLD)

Definition at line 1908 of file m_octree_mg_1d.t.

◆ mg_deallocate_storage()

subroutine, public m_octree_mg_1d::mg_deallocate_storage ( type(mg_t), intent(inout) mg )

Deallocate all allocatable arrays.

Definition at line 620 of file m_octree_mg_1d.t.

◆ mg_fas_fmg()

subroutine, public m_octree_mg_1d::mg_fas_fmg	(	type(mg_t), intent(inout)	mg,
		logical, intent(in)	have_guess,
		real(dp), intent(out), optional	max_res
	)

Perform FAS-FMG cycle (full approximation scheme, full multigrid).

Parameters

[in]	have_guess	If false, start from phi = 0
[out]	max_res	Store max(abs(residual))

Definition at line 2837 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_fas_vcycle()

subroutine, public m_octree_mg_1d::mg_fas_vcycle	(	type(mg_t), intent(inout)	mg,
		integer, intent(in), optional	highest_lvl,
		real(dp), intent(out), optional	max_res,
		logical, intent(in), optional	standalone
	)

Perform FAS V-cycle (full approximation scheme).

Parameters

[in]	highest_lvl	Maximum level for V-cycle
[out]	max_res	Store max(abs(residual))
[in]	standalone	Whether the V-cycle is called by itself (default: true)

Definition at line 2903 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_fill_ghost_cells()

subroutine, public m_octree_mg_1d::mg_fill_ghost_cells	(	type(mg_t)	mg,
		integer, intent(in)	iv
	)

Fill ghost cells at all grid levels.

Parameters

[in] iv Index of variable

Definition at line 2084 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_fill_ghost_cells_lvl()

subroutine, public m_octree_mg_1d::mg_fill_ghost_cells_lvl	(	type(mg_t)	mg,
		integer, intent(in)	lvl,
		integer, intent(in)	iv
	)

Fill ghost cells at a grid level.

Parameters

[in] iv Index of variable

Definition at line 2095 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_get_child_offset()

pure integer function, dimension(1), public m_octree_mg_1d::mg_get_child_offset	(	type(mg_t), intent(in)	mg,
		integer, intent(in)	id
	)

Get the offset of a box with respect to its parent (e.g. in 2d, there can be a child at offset 0,0, one at n_cell/2,0, one at 0,n_cell/2 and one at n_cell/2, n_cell/2)

Definition at line 516 of file m_octree_mg_1d.t.

◆ mg_get_face_coords()

subroutine, public m_octree_mg_1d::mg_get_face_coords	(	type(mg_box_t), intent(in)	box,
		integer, intent(in)	nb,
		integer, intent(in)	nc,
		real(dp), dimension(2), intent(out)	x
	)

Get coordinates at the face of a box.

Definition at line 555 of file m_octree_mg_1d.t.

◆ mg_ghost_cell_buffer_size()

subroutine, public m_octree_mg_1d::mg_ghost_cell_buffer_size	(	type(mg_t), intent(inout)	mg,
		integer, dimension(0:mg%n_cpu-1), intent(out)	n_send,
		integer, dimension(0:mg%n_cpu-1), intent(out)	n_recv,
		integer, intent(out)	dsize
	)

Specify minimum buffer size (per process) for communication.

Definition at line 1987 of file m_octree_mg_1d.t.

◆ mg_has_children()

elemental logical function, public m_octree_mg_1d::mg_has_children ( type(mg_box_t), intent(in) box )

Return .true. if a box has children.

Definition at line 497 of file m_octree_mg_1d.t.

◆ mg_highest_uniform_lvl()

pure integer function, public m_octree_mg_1d::mg_highest_uniform_lvl ( type(mg_t), intent(in) mg )

Definition at line 529 of file m_octree_mg_1d.t.

◆ mg_ix_to_ichild()

integer function, public m_octree_mg_1d::mg_ix_to_ichild ( integer, dimension(1), intent(in) ix )

Compute the child index for a box with spatial index ix. With child index we mean the index in the children(:) array of its parent.

Parameters

[in] ix Spatial index of the box

Definition at line 507 of file m_octree_mg_1d.t.

◆ mg_load_balance()

subroutine, public m_octree_mg_1d::mg_load_balance ( type(mg_t), intent(inout) mg )

Load balance all boxes in the multigrid tree. Compared to mg_load_balance_simple, this method does a better job of setting the ranks of parent boxes.

Note that in a typical application the load balancing of the leaves is already determined, then mg_load_balance_parents can be used.

Definition at line 1532 of file m_octree_mg_1d.t.

◆ mg_load_balance_parents()

subroutine, public m_octree_mg_1d::mg_load_balance_parents ( type(mg_t), intent(inout) mg )

Load balance the parents (non-leafs). Assign them to the rank that has most children.

Definition at line 1601 of file m_octree_mg_1d.t.

◆ mg_load_balance_simple()

subroutine, public m_octree_mg_1d::mg_load_balance_simple ( type(mg_t), intent(inout) mg )

Load balance all boxes in the multigrid tree, by simply distributing the load per grid level. This method will only work well for uniform grids.

Note that in a typical application the load balancing of the leaves is already determined, then mg_load_balance_parents can be used.

Definition at line 1483 of file m_octree_mg_1d.t.

◆ mg_number_of_unknowns()

integer(i8) function, public m_octree_mg_1d::mg_number_of_unknowns ( type(mg_t), intent(in) mg )

Determine total number of unknowns (on leaves)

Definition at line 542 of file m_octree_mg_1d.t.

◆ mg_phi_bc_store()

subroutine, public m_octree_mg_1d::mg_phi_bc_store ( type(mg_t), intent(inout) mg )

Store boundary conditions for the solution variable, this can speed up calculations if the same boundary conditions are re-used.

Definition at line 2036 of file m_octree_mg_1d.t.

◆ mg_prolong()

subroutine, public m_octree_mg_1d::mg_prolong	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	lvl,
		integer, intent(in)	iv,
		integer, intent(in)	iv_to,
		procedure(mg_box_prolong)	method,
		logical, intent(in)	add
	)

Prolong variable iv from lvl to variable iv_to at lvl+1.

Parameters

[in]	lvl	Level to prolong from
[in]	iv	Source variable
[in]	iv_to	Target variable
	method	Prolongation method
[in]	add	If true, add to current values

Definition at line 2561 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_prolong_buffer_size()

subroutine, public m_octree_mg_1d::mg_prolong_buffer_size	(	type(mg_t), intent(inout)	mg,
		integer, dimension(0:mg%n_cpu-1), intent(out)	n_send,
		integer, dimension(0:mg%n_cpu-1), intent(out)	n_recv,
		integer, intent(out)	dsize
	)

Specify minimum buffer size (per process) for communication.

Definition at line 2526 of file m_octree_mg_1d.t.

◆ mg_prolong_sparse()

subroutine, public m_octree_mg_1d::mg_prolong_sparse	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	p_id,
		integer, dimension(1), intent(in)	dix,
		integer, intent(in)	nc,
		integer, intent(in)	iv,
		real(dp), dimension(nc), intent(out)	fine
	)

Prolong from a parent to a child with index offset dix.

Parameters

[in]	p_id	Id of parent
[in]	dix	Offset of child in parent grid
[in]	nc	Child grid size
[in]	iv	Prolong from this variable
[out]	fine	Prolonged values

Definition at line 2668 of file m_octree_mg_1d.t.

◆ mg_restrict()

subroutine, public m_octree_mg_1d::mg_restrict	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	iv
	)

Restrict all levels.

Definition at line 3270 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_restrict_buffer_size()

subroutine, public m_octree_mg_1d::mg_restrict_buffer_size	(	type(mg_t), intent(inout)	mg,
		integer, dimension(0:mg%n_cpu-1), intent(out)	n_send,
		integer, dimension(0:mg%n_cpu-1), intent(out)	n_recv,
		integer, intent(out)	dsize
	)

Specify minimum buffer size (per process) for communication.

Definition at line 3214 of file m_octree_mg_1d.t.

◆ mg_restrict_lvl()

subroutine, public m_octree_mg_1d::mg_restrict_lvl	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	iv,
		integer, intent(in)	lvl
	)

Restrict from lvl to lvl-1.

Definition at line 3281 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_set_leaves_parents()

subroutine, public m_octree_mg_1d::mg_set_leaves_parents	(	type(mg_box_t), dimension(:), intent(in)	boxes,
		type(mg_lvl_t), intent(inout)	level
	)

Create a list of leaves and a list of parents for a level.

Parameters

[in]	boxes	List of boxes
[in,out]	level	Level type which contains the indices of boxes

Definition at line 1339 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_set_methods()

subroutine, public m_octree_mg_1d::mg_set_methods ( type(mg_t), intent(inout) mg )

Definition at line 2784 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_set_neighbors_lvl()

subroutine, public m_octree_mg_1d::mg_set_neighbors_lvl	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	lvl
	)

Definition at line 1257 of file m_octree_mg_1d.t.

◆ mg_set_next_level_ids()

subroutine, public m_octree_mg_1d::mg_set_next_level_ids	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	lvl
	)

Definition at line 1268 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_set_refinement_boundaries()

subroutine, public m_octree_mg_1d::mg_set_refinement_boundaries	(	type(mg_box_t), dimension(:), intent(in)	boxes,
		type(mg_lvl_t), intent(inout)	level
	)

Create a list of refinement boundaries (from the coarse side)

Definition at line 1377 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ mg_timer_end()

subroutine, public m_octree_mg_1d::mg_timer_end ( type(mg_timer_t), intent(inout) timer )

Definition at line 590 of file m_octree_mg_1d.t.

◆ mg_timer_start()

subroutine, public m_octree_mg_1d::mg_timer_start ( type(mg_timer_t), intent(inout) timer )

Definition at line 584 of file m_octree_mg_1d.t.

◆ mg_timers_show()

subroutine, public m_octree_mg_1d::mg_timers_show ( type(mg_t), intent(in) mg )

Definition at line 596 of file m_octree_mg_1d.t.

◆ sort_and_transfer_buffers()

subroutine, public m_octree_mg_1d::sort_and_transfer_buffers	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	dsize
	)

Definition at line 1931 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ sort_sendbuf()

subroutine m_octree_mg_1d::sort_sendbuf	(	type(mg_buf_t), intent(inout)	gc,
		integer, intent(in)	dsize
	)

Sort send buffers according to the idbuf array.

Parameters

[in] dsize Size of send buffer elements

Definition at line 1963 of file m_octree_mg_1d.t.

◆ subtract_mean()

subroutine m_octree_mg_1d::subtract_mean	(	type(mg_t), intent(inout)	mg,
		integer, intent(in)	iv,
		logical, intent(in)	include_ghostcells
	)

Definition at line 2998 of file m_octree_mg_1d.t.

Here is the call graph for this function:

◆ vhelmholtz_set_lambda()

subroutine, public m_octree_mg_1d::vhelmholtz_set_lambda ( real(dp), intent(in) lambda )

Definition at line 1021 of file m_octree_mg_1d.t.

◆ vhelmholtz_set_methods()

subroutine, public m_octree_mg_1d::vhelmholtz_set_methods ( type(mg_t), intent(inout) mg )

Definition at line 989 of file m_octree_mg_1d.t.

◆ vlaplacian_set_methods()

subroutine, public m_octree_mg_1d::vlaplacian_set_methods ( type(mg_t), intent(inout) mg )

Definition at line 1700 of file m_octree_mg_1d.t.

Variable Documentation

◆ ahelmholtz_lambda

real(dp), public, protected m_octree_mg_1d::ahelmholtz_lambda = 0.0_dp

Module which contains multigrid procedures for a Helmholtz operator of the form: div(D grad(phi)) - lambda*phi = f, where D has a smooth spatial variation and a component in each spatial direction.

The lambda used for the Helmholtz equation (should be >= 0)

Definition at line 488 of file m_octree_mg_1d.t.

◆ dp

integer, parameter, public m_octree_mg_1d::dp = kind(0.0d0)

Type of reals.

Definition at line 30 of file m_octree_mg_1d.t.

◆ helmholtz_lambda

real(dp), public, protected m_octree_mg_1d::helmholtz_lambda = 0.0_dp

Module for load balancing a tree (that already has been constructed). The load balancing determines which ranks (MPI processes) allocated physical storage for boxes. The tree structure itself is present on all processes.

Module which contains multigrid procedures for a variable-coefficient Laplacian operator, assuming the variation is smooth Module which contains multigrid procedures for a Helmholtz operator of the form: laplacian(phi) - lambda*phi = f The lambda used for the Helmholtz equation (should be >= 0)

Definition at line 443 of file m_octree_mg_1d.t.

◆ i8

integer, parameter, public m_octree_mg_1d::i8 = selected_int_kind(18)

Type for 64-bit integers.

Definition at line 33 of file m_octree_mg_1d.t.

◆ mg_ahelmholtz

integer, parameter, public m_octree_mg_1d::mg_ahelmholtz = 5

Indicates a anisotropic-coefficient Helmholtz equation.

Definition at line 48 of file m_octree_mg_1d.t.

◆ mg_bc_continuous

integer, parameter, public m_octree_mg_1d::mg_bc_continuous = -12

Value to indicate a continuous boundary condition.

Definition at line 97 of file m_octree_mg_1d.t.

◆ mg_bc_dirichlet

integer, parameter, public m_octree_mg_1d::mg_bc_dirichlet = -10

Value to indicate a Dirichlet boundary condition.

Definition at line 91 of file m_octree_mg_1d.t.

◆ mg_bc_neumann

integer, parameter, public m_octree_mg_1d::mg_bc_neumann = -11

Value to indicate a Neumann boundary condition.

Definition at line 94 of file m_octree_mg_1d.t.

◆ mg_cartesian

integer, parameter, public m_octree_mg_1d::mg_cartesian = 1

Cartesian coordinate system.

Definition at line 51 of file m_octree_mg_1d.t.

◆ mg_child_adj_nb

integer, dimension(1, 2), parameter, public m_octree_mg_1d::mg_child_adj_nb = reshape([1,2], [1,2])

Definition at line 115 of file m_octree_mg_1d.t.

◆ mg_child_dix

integer, dimension(1, 2), parameter, public m_octree_mg_1d::mg_child_dix = reshape([0,1], [1,2])

Definition at line 111 of file m_octree_mg_1d.t.

◆ mg_child_low

logical, dimension(1, 2), parameter, public m_octree_mg_1d::mg_child_low = reshape([.true., .false.], [1, 2])

Definition at line 117 of file m_octree_mg_1d.t.

◆ mg_child_rev

integer, dimension(2, 1), parameter, public m_octree_mg_1d::mg_child_rev = reshape([2,1], [2,1])

Definition at line 113 of file m_octree_mg_1d.t.

◆ mg_cylindrical

integer, parameter, public m_octree_mg_1d::mg_cylindrical = 2

Cylindrical coordinate system.

Definition at line 53 of file m_octree_mg_1d.t.

◆ mg_helmholtz

integer, parameter, public m_octree_mg_1d::mg_helmholtz = 3

Indicates a constant-coefficient Helmholtz equation.

Definition at line 42 of file m_octree_mg_1d.t.

◆ mg_iold

integer, parameter, public m_octree_mg_1d::mg_iold = 3

Index of previous solution (used for correction)

Definition at line 73 of file m_octree_mg_1d.t.

◆ mg_iphi

integer, parameter, public m_octree_mg_1d::mg_iphi = 1

Index of solution.

Definition at line 69 of file m_octree_mg_1d.t.

◆ mg_ires

integer, parameter, public m_octree_mg_1d::mg_ires = 4

Index of residual.

Definition at line 75 of file m_octree_mg_1d.t.

◆ mg_irhs

integer, parameter, public m_octree_mg_1d::mg_irhs = 2

Index of right-hand side.

Definition at line 71 of file m_octree_mg_1d.t.

◆ mg_iveps

integer, parameter, public m_octree_mg_1d::mg_iveps = 5

Index of the variable coefficient (at cell centers)

Definition at line 78 of file m_octree_mg_1d.t.

◆ mg_iveps1

integer, parameter, public m_octree_mg_1d::mg_iveps1 = 5

Indexes of anisotropic variable coefficients.

Definition at line 81 of file m_octree_mg_1d.t.

◆ mg_iveps2

integer, parameter, public m_octree_mg_1d::mg_iveps2 = 6

Definition at line 82 of file m_octree_mg_1d.t.

◆ mg_iveps3

integer, parameter, public m_octree_mg_1d::mg_iveps3 = 7

Definition at line 83 of file m_octree_mg_1d.t.

◆ mg_laplacian

integer, parameter, public m_octree_mg_1d::mg_laplacian = 1

Indicates a standard Laplacian.

Definition at line 36 of file m_octree_mg_1d.t.

◆ mg_lvl_hi

integer, parameter, public m_octree_mg_1d::mg_lvl_hi = 20

Maximum allowed grid level.

Definition at line 88 of file m_octree_mg_1d.t.

◆ mg_lvl_lo

integer, parameter, public m_octree_mg_1d::mg_lvl_lo = -20

Minimum allowed grid level.

Definition at line 86 of file m_octree_mg_1d.t.

◆ mg_max_num_vars

integer, parameter, public m_octree_mg_1d::mg_max_num_vars = 10

Maximum number of variables.

Definition at line 67 of file m_octree_mg_1d.t.

◆ mg_max_timers

integer, parameter, public m_octree_mg_1d::mg_max_timers = 20

Maximum number of timers to use.

Definition at line 105 of file m_octree_mg_1d.t.

◆ mg_ndim

integer, parameter, public m_octree_mg_1d::mg_ndim = 1

Problem dimension.

Definition at line 62 of file m_octree_mg_1d.t.

◆ mg_neighb_dim

integer, dimension(2), parameter, public m_octree_mg_1d::mg_neighb_dim = [1, 1]

Definition at line 134 of file m_octree_mg_1d.t.

◆ mg_neighb_dix

integer, dimension(1, 2), parameter, public m_octree_mg_1d::mg_neighb_dix = reshape([-1,1], [1,2])

Definition at line 125 of file m_octree_mg_1d.t.

◆ mg_neighb_high_pm

integer, dimension(2), parameter, public m_octree_mg_1d::mg_neighb_high_pm = [-1, 1]

Definition at line 129 of file m_octree_mg_1d.t.

◆ mg_neighb_highx

integer, parameter, public m_octree_mg_1d::mg_neighb_highx = 2

Definition at line 122 of file m_octree_mg_1d.t.

◆ mg_neighb_low

logical, dimension(2), parameter, public m_octree_mg_1d::mg_neighb_low = [.true., .false.]

Definition at line 127 of file m_octree_mg_1d.t.

◆ mg_neighb_lowx

integer, parameter, public m_octree_mg_1d::mg_neighb_lowx = 1

Definition at line 121 of file m_octree_mg_1d.t.

◆ mg_neighb_rev

integer, dimension(2), parameter, public m_octree_mg_1d::mg_neighb_rev = [2, 1]

Definition at line 132 of file m_octree_mg_1d.t.

◆ mg_no_box

integer, parameter, public m_octree_mg_1d::mg_no_box = 0

Special value that indicates there is no box.

Definition at line 100 of file m_octree_mg_1d.t.

◆ mg_num_children

integer, parameter, public m_octree_mg_1d::mg_num_children = 2

Definition at line 108 of file m_octree_mg_1d.t.

◆ mg_num_neighbors

integer, parameter, public m_octree_mg_1d::mg_num_neighbors = 2

Definition at line 120 of file m_octree_mg_1d.t.

◆ mg_num_vars

integer, parameter, public m_octree_mg_1d::mg_num_vars = 4

Number of predefined multigrid variables.

Definition at line 65 of file m_octree_mg_1d.t.

◆ mg_physical_boundary

integer, parameter, public m_octree_mg_1d::mg_physical_boundary = -1

Special value that indicates there is a physical boundary.

Definition at line 102 of file m_octree_mg_1d.t.

◆ mg_smoother_gs

integer, parameter, public m_octree_mg_1d::mg_smoother_gs = 1

Definition at line 57 of file m_octree_mg_1d.t.

◆ mg_smoother_gsrb

integer, parameter, public m_octree_mg_1d::mg_smoother_gsrb = 2

Definition at line 58 of file m_octree_mg_1d.t.

◆ mg_smoother_jacobi

integer, parameter, public m_octree_mg_1d::mg_smoother_jacobi = 3

Definition at line 59 of file m_octree_mg_1d.t.

◆ mg_spherical

integer, parameter, public m_octree_mg_1d::mg_spherical = 3

Spherical coordinate system.

Definition at line 55 of file m_octree_mg_1d.t.

◆ mg_vhelmholtz

integer, parameter, public m_octree_mg_1d::mg_vhelmholtz = 4

Indicates a variable-coefficient Helmholtz equation.

Definition at line 45 of file m_octree_mg_1d.t.

◆ mg_vlaplacian

integer, parameter, public m_octree_mg_1d::mg_vlaplacian = 2

Indicates a variable-coefficient Laplacian.

Definition at line 39 of file m_octree_mg_1d.t.

◆ vhelmholtz_lambda

real(dp), public, protected m_octree_mg_1d::vhelmholtz_lambda = 0.0_dp

Module for implicitly solving diffusion equations.

Module which contains multigrid procedures for a Laplacian operator Module which contains multigrid procedures for a Helmholtz operator of the form: div(D grad(phi)) - lambda*phi = f, where D has a smooth spatial variation The lambda used for the Helmholtz equation (should be >= 0)

Definition at line 384 of file m_octree_mg_1d.t.

Data Types

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ ahelmholtz_set_lambda()

◆ ahelmholtz_set_methods()

◆ diffusion_solve()

◆ diffusion_solve_acoeff()

◆ diffusion_solve_vcoeff()

◆ get_sum()

◆ helmholtz_set_lambda()

◆ helmholtz_set_methods()

◆ laplacian_set_methods()

◆ mg_add_children()

◆ mg_add_timer()

◆ mg_allocate_storage()

◆ mg_apply_op()

◆ mg_build_rectangle()

◆ mg_comm_init()

◆ mg_deallocate_storage()

◆ mg_fas_fmg()

◆ mg_fas_vcycle()

◆ mg_fill_ghost_cells()

◆ mg_fill_ghost_cells_lvl()

◆ mg_get_child_offset()

◆ mg_get_face_coords()

◆ mg_ghost_cell_buffer_size()

◆ mg_has_children()

◆ mg_highest_uniform_lvl()

◆ mg_ix_to_ichild()

◆ mg_load_balance()

◆ mg_load_balance_parents()

◆ mg_load_balance_simple()

◆ mg_number_of_unknowns()

◆ mg_phi_bc_store()

◆ mg_prolong()

◆ mg_prolong_buffer_size()

◆ mg_prolong_sparse()

◆ mg_restrict()

◆ mg_restrict_buffer_size()

◆ mg_restrict_lvl()

◆ mg_set_leaves_parents()

◆ mg_set_methods()

◆ mg_set_neighbors_lvl()

◆ mg_set_next_level_ids()

◆ mg_set_refinement_boundaries()

◆ mg_timer_end()

◆ mg_timer_start()

◆ mg_timers_show()

◆ sort_and_transfer_buffers()

◆ sort_sendbuf()

◆ subtract_mean()

◆ vhelmholtz_set_lambda()

◆ vhelmholtz_set_methods()

◆ vlaplacian_set_methods()

Variable Documentation

◆ ahelmholtz_lambda

◆ dp

◆ helmholtz_lambda

◆ i8

◆ mg_ahelmholtz

◆ mg_bc_continuous

◆ mg_bc_dirichlet

◆ mg_bc_neumann

◆ mg_cartesian

◆ mg_child_adj_nb

◆ mg_child_dix

◆ mg_child_low

◆ mg_child_rev

◆ mg_cylindrical

◆ mg_helmholtz

◆ mg_iold

◆ mg_iphi

◆ mg_ires

◆ mg_irhs

◆ mg_iveps

◆ mg_iveps1

◆ mg_iveps2

◆ mg_iveps3

◆ mg_laplacian