This module defines non-local operators. More...

Detailed Description

This module defines non-local operators.

Data Types
type	nl_operator_index_t
	index type for non-local operators More...

type	nl_operator_t
	data type for non local operators More...

interface	op_is_available

Functions/Subroutines
subroutine, public	nl_operator_global_init (namespace)
	initialize global settings for non-local operators More...

subroutine, public	nl_operator_global_end ()

subroutine, public	nl_operator_init (op, label, symm)
	initialize an instance of a non-local operator by setting the label More...

subroutine, public	nl_operator_build (space, mesh, op, np, const_w, regenerate)
	Creates the nonlocal operators for the stencils used for finite differences. More...

subroutine, public	nl_operator_output_weights (this)

subroutine, public	nl_operator_end (op)

subroutine	nl_operator_clear_gpu_buffers (op)

integer pure function, public	nl_operator_get_index (op, is, ip)

subroutine, public	nl_operator_allocate_gpu_buffers (op)

subroutine, public	nl_operator_update_gpu_buffers (op)

integer pure function, public	nl_operator_np_zero_bc (op)

subroutine, public	nl_operator_remove_zero_weight_points (op, space, mesh)
	Removes the zero-weight points for constant weight stencils. More...

subroutine	group_by_pairs_sym (size, ldf, offsets, wre, ri, nri, npairs, wpair, pair_pos, pair_neg, wcenter)
	Take a list of weights and offsets and build pairs of symmetric points with common weights. More...

subroutine	group_by_pairs_antisym (size, ldf, offsets, wre, ri, nri, npairs, wpair, pair_pos, pair_neg)
	Take a list of weights and offsets and build pairs of symmetric points with common weights. More...

subroutine, public	nl_operator_build_symmetric_weights (op, max_size)
	Builds (or rebuild) the necessary arrays for symmetric and antisymmetric stencils. More...

subroutine	reallocate_array (ri, stencil_size, nri, old_size, new_size)
	Reallocate an ri array. More...

subroutine, public	dnl_operator_operate_batch (op, fi, fo, ghost_update, profile, points, factor, async)

subroutine, public	dnl_operator_operate (op, fi, fo, ghost_update, profile, points)

subroutine, public	dnl_operator_operate_diag (op, fo)

subroutine, public	znl_operator_operate_batch (op, fi, fo, ghost_update, profile, points, factor, async)

subroutine, public	znl_operator_operate (op, fi, fo, ghost_update, profile, points)

subroutine, public	znl_operator_operate_diag (op, fo)

Variables
integer, parameter, public	op_general = 1

integer, parameter, public	op_symmetric = 2

integer, parameter, public	op_antisymmetric = 3

integer, parameter	op_fortran = 0

integer, parameter	op_vec = 1

integer, parameter	op_min = OP_FORTRAN

integer, parameter	op_max = OP_VEC

integer, parameter	op_invmap = 1

integer, parameter	op_map = 2

integer, parameter, public	op_all = 3

integer, parameter, public	op_inner = 1

integer, parameter, public	op_outer = 2

integer	dfunction_global = -1

integer	zfunction_global = -1

integer	function_accel

logical	use_symmetries = .false.

Function/Subroutine Documentation

◆ nl_operator_global_init()

subroutine, public nl_operator_oct_m::nl_operator_global_init ( type(namespace_t), intent(in) namespace )

initialize global settings for non-local operators

Definition at line 276 of file nl_operator.F90.

◆ nl_operator_global_end()

subroutine, public nl_operator_oct_m::nl_operator_global_end

Definition at line 352 of file nl_operator.F90.

◆ nl_operator_init()

subroutine, public nl_operator_oct_m::nl_operator_init	(	type(nl_operator_t), intent(inout)	op,
		character(len=*), intent(in)	label,
		integer, intent(in), optional	symm
	)

initialize an instance of a non-local operator by setting the label

Definition at line 362 of file nl_operator.F90.

◆ nl_operator_build()

subroutine, public nl_operator_oct_m::nl_operator_build	(	class(space_t), intent(in)	space,
		type(mesh_t), intent(in), target	mesh,
		type(nl_operator_t), intent(inout)	op,
		integer, intent(in)	np,
		logical, intent(in), optional	const_w,
		logical, intent(in), optional	regenerate
	)

Creates the nonlocal operators for the stencils used for finite differences.

This code loops over grid points and finds the stencil to be used by looking at the coordinates of the points in the stencils. All the consecutive points that are using the sames shifts in relative coordinates are considered to be within one stencil application. The number of them is opri The map oprimap contains the map (mesh index ip) -> stencil (between 1 to opnri) The inverse map oprimap_inv contains the start and end indices for a given stencil, such that oprimap_inv(ii):oprimap_inv(ii+1) is the list of points sharing the same stencil index With this, all points within oprimap_inv(ii):oprimap_inv(ii+1) have the same value in oprimap

Parameters

[in]	np	Number of (local) points.
[in]	const_w	are the weights constant (independent of the point)

Definition at line 388 of file nl_operator.F90.

◆ nl_operator_output_weights()

subroutine, public nl_operator_oct_m::nl_operator_output_weights ( type(nl_operator_t), intent(inout) this )

Definition at line 678 of file nl_operator.F90.

◆ nl_operator_end()

subroutine, public nl_operator_oct_m::nl_operator_end ( type(nl_operator_t), intent(inout) op )

Definition at line 709 of file nl_operator.F90.

◆ nl_operator_clear_gpu_buffers()

subroutine nl_operator_oct_m::nl_operator_clear_gpu_buffers ( type(nl_operator_t), intent(inout) op )

private

Definition at line 746 of file nl_operator.F90.

◆ nl_operator_get_index()

integer pure function, public nl_operator_oct_m::nl_operator_get_index	(	type(nl_operator_t), intent(in)	op,
		integer, intent(in)	is,
		integer, intent(in)	ip
	)

Definition at line 786 of file nl_operator.F90.

◆ nl_operator_allocate_gpu_buffers()

subroutine, public nl_operator_oct_m::nl_operator_allocate_gpu_buffers ( type(nl_operator_t), intent(inout) op )

Definition at line 796 of file nl_operator.F90.

◆ nl_operator_update_gpu_buffers()

subroutine, public nl_operator_oct_m::nl_operator_update_gpu_buffers ( type(nl_operator_t), intent(inout) op )

Definition at line 817 of file nl_operator.F90.

◆ nl_operator_np_zero_bc()

integer pure function, public nl_operator_oct_m::nl_operator_np_zero_bc ( type(nl_operator_t), intent(in) op )

Definition at line 857 of file nl_operator.F90.

◆ nl_operator_remove_zero_weight_points()

subroutine, public nl_operator_oct_m::nl_operator_remove_zero_weight_points	(	type(nl_operator_t), intent(inout)	op,
		type(space_t), intent(in)	space,
		class(mesh_t), intent(in)	mesh
	)

Removes the zero-weight points for constant weight stencils.

Definition at line 873 of file nl_operator.F90.

◆ group_by_pairs_sym()

subroutine nl_operator_oct_m::group_by_pairs_sym	(	integer, intent(in)	size,
		integer, intent(in)	ldf,
		integer, dimension(:, :), intent(in)	offsets,
		real(real64), dimension(:), intent(in)	wre,
		integer, dimension(:, :), intent(in)	ri,
		integer, intent(in)	nri,
		integer, intent(out)	npairs,
		real(real64), dimension(:), intent(inout)	wpair,
		integer, dimension(:,:), intent(inout)	pair_pos,
		integer, dimension(:,:), intent(inout)	pair_neg,
		real(real64), intent(out)	wcenter
	)

private

Take a list of weights and offsets and build pairs of symmetric points with common weights.

Definition at line 914 of file nl_operator.F90.

◆ group_by_pairs_antisym()

subroutine nl_operator_oct_m::group_by_pairs_antisym	(	integer, intent(in)	size,
		integer, intent(in)	ldf,
		integer, dimension(:, :), intent(in)	offsets,
		real(real64), dimension(:), intent(in)	wre,
		integer, dimension(:, :), intent(in)	ri,
		integer, intent(in)	nri,
		integer, intent(out)	npairs,
		real(real64), dimension(:), intent(inout)	wpair,
		integer, dimension(:,:), intent(inout)	pair_pos,
		integer, dimension(:,:), intent(inout)	pair_neg
	)

private

Take a list of weights and offsets and build pairs of symmetric points with common weights.

Definition at line 987 of file nl_operator.F90.

◆ nl_operator_build_symmetric_weights()

subroutine, public nl_operator_oct_m::nl_operator_build_symmetric_weights	(	type(nl_operator_t), intent(inout)	op,
		integer, intent(in), optional	max_size
	)

Builds (or rebuild) the necessary arrays for symmetric and antisymmetric stencils.

Parameters

[in] max_size For reallocation

Definition at line 1053 of file nl_operator.F90.

◆ reallocate_array()

subroutine nl_operator_oct_m::reallocate_array	(	integer, dimension(:,:,:), intent(inout), allocatable	ri,
		integer, intent(in)	stencil_size,
		integer, intent(in)	nri,
		integer, intent(in)	old_size,
		integer, intent(in)	new_size
	)

private

Reallocate an ri array.

Definition at line 1139 of file nl_operator.F90.

◆ dnl_operator_operate_batch()

subroutine, public nl_operator_oct_m::dnl_operator_operate_batch	(	type(nl_operator_t), intent(inout), target	op,
		type(batch_t), intent(inout), target	fi,
		type(batch_t), intent(inout)	fo,
		logical, intent(in), optional	ghost_update,
		logical, intent(in), optional	profile,
		integer, intent(in), optional	points,
		real(real64), intent(in), optional	factor,
		logical, intent(in), optional	async
	)

Parameters

[in,out] fo this should be target, but old ifort 9.1 segfaults with it

Symmetry of the stencil can be used to group points by pairs

Definition at line 1222 of file nl_operator.F90.

◆ dnl_operator_operate()

subroutine, public nl_operator_oct_m::dnl_operator_operate	(	type(nl_operator_t), intent(inout)	op,
		real(real64), dimension(:), intent(inout), contiguous	fi,
		real(real64), dimension(:), intent(out), target, contiguous	fo,
		logical, intent(in), optional	ghost_update,
		logical, intent(in), optional	profile,
		integer, intent(in), optional	points
	)

Parameters

[in,out] fi fi(opnp_part)

Definition at line 1862 of file nl_operator.F90.

◆ dnl_operator_operate_diag()

subroutine, public nl_operator_oct_m::dnl_operator_operate_diag	(	type(nl_operator_t), intent(in)	op,
		real(real64), dimension(:), intent(out)	fo
	)

Definition at line 1887 of file nl_operator.F90.

◆ znl_operator_operate_batch()

subroutine, public nl_operator_oct_m::znl_operator_operate_batch	(	type(nl_operator_t), intent(inout), target	op,
		type(batch_t), intent(inout), target	fi,
		type(batch_t), intent(inout)	fo,
		logical, intent(in), optional	ghost_update,
		logical, intent(in), optional	profile,
		integer, intent(in), optional	points,
		real(real64), intent(in), optional	factor,
		logical, intent(in), optional	async
	)

Parameters

[in,out] fo this should be target, but old ifort 9.1 segfaults with it

Symmetry of the stencil can be used to group points by pairs

Definition at line 1979 of file nl_operator.F90.

◆ znl_operator_operate()

subroutine, public nl_operator_oct_m::znl_operator_operate	(	type(nl_operator_t), intent(inout)	op,
		complex(real64), dimension(:), intent(inout), contiguous	fi,
		complex(real64), dimension(:), intent(out), target, contiguous	fo,
		logical, intent(in), optional	ghost_update,
		logical, intent(in), optional	profile,
		integer, intent(in), optional	points
	)

Parameters

[in,out] fi fi(opnp_part)

Definition at line 2619 of file nl_operator.F90.

◆ znl_operator_operate_diag()

subroutine, public nl_operator_oct_m::znl_operator_operate_diag	(	type(nl_operator_t), intent(in)	op,
		complex(real64), dimension(:), intent(out)	fo
	)

Definition at line 2644 of file nl_operator.F90.

Variable Documentation

◆ op_general

integer, parameter, public nl_operator_oct_m::op_general = 1

Definition at line 181 of file nl_operator.F90.

◆ op_symmetric

integer, parameter, public nl_operator_oct_m::op_symmetric = 2

Definition at line 181 of file nl_operator.F90.

◆ op_antisymmetric

integer, parameter, public nl_operator_oct_m::op_antisymmetric = 3

Definition at line 181 of file nl_operator.F90.

◆ op_fortran

integer, parameter nl_operator_oct_m::op_fortran = 0

private

Definition at line 247 of file nl_operator.F90.

◆ op_vec

integer, parameter nl_operator_oct_m::op_vec = 1

private

Definition at line 247 of file nl_operator.F90.

◆ op_min

integer, parameter nl_operator_oct_m::op_min = OP_FORTRAN

private

Definition at line 247 of file nl_operator.F90.

◆ op_max

integer, parameter nl_operator_oct_m::op_max = OP_VEC

private

Definition at line 247 of file nl_operator.F90.

◆ op_invmap

integer, parameter nl_operator_oct_m::op_invmap = 1

private

Definition at line 253 of file nl_operator.F90.

◆ op_map

integer, parameter nl_operator_oct_m::op_map = 2

private

Definition at line 253 of file nl_operator.F90.

◆ op_all

integer, parameter, public nl_operator_oct_m::op_all = 3

Definition at line 257 of file nl_operator.F90.

◆ op_inner

integer, parameter, public nl_operator_oct_m::op_inner = 1

Definition at line 257 of file nl_operator.F90.

◆ op_outer

integer, parameter, public nl_operator_oct_m::op_outer = 2

Definition at line 257 of file nl_operator.F90.

◆ dfunction_global

integer nl_operator_oct_m::dfunction_global = -1

private

Definition at line 266 of file nl_operator.F90.

◆ zfunction_global

integer nl_operator_oct_m::zfunction_global = -1

private

Definition at line 267 of file nl_operator.F90.

◆ function_accel

integer nl_operator_oct_m::function_accel

private

Definition at line 268 of file nl_operator.F90.

◆ use_symmetries

logical nl_operator_oct_m::use_symmetries = .false.

private

Definition at line 269 of file nl_operator.F90.

Detailed Description

Data Types

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ nl_operator_global_init()

◆ nl_operator_global_end()

◆ nl_operator_init()

◆ nl_operator_build()

◆ nl_operator_output_weights()

◆ nl_operator_end()

◆ nl_operator_clear_gpu_buffers()

◆ nl_operator_get_index()

◆ nl_operator_allocate_gpu_buffers()

◆ nl_operator_update_gpu_buffers()

◆ nl_operator_np_zero_bc()

◆ nl_operator_remove_zero_weight_points()

◆ group_by_pairs_sym()

◆ group_by_pairs_antisym()

◆ nl_operator_build_symmetric_weights()

◆ reallocate_array()

◆ dnl_operator_operate_batch()

◆ dnl_operator_operate()

◆ dnl_operator_operate_diag()

◆ znl_operator_operate_batch()

◆ znl_operator_operate()

◆ znl_operator_operate_diag()

Variable Documentation

◆ op_general

◆ op_symmetric

◆ op_antisymmetric

◆ op_fortran

◆ op_vec

◆ op_min

◆ op_max

◆ op_invmap

◆ op_map

◆ op_all

◆ op_inner

◆ op_outer

◆ dfunction_global

◆ zfunction_global

◆ function_accel

◆ use_symmetries