expression_dirichlet::expression_dirichlet_t Type Reference

Dirichlet condition \( x = g(x, y, z, t) \) on \(\partial \Omega\), where \( g \) is a mathematical expression given in the case file. More...

Inheritance diagram for expression_dirichlet::expression_dirichlet_t:

Collaboration diagram for expression_dirichlet::expression_dirichlet_t:

Public Member Functions
procedure, pass(this)	init (this, coef, json)
	Constructor from JSON.
procedure, pass(this)	init_from_components (this, coef, str)
	Constructor from components.
procedure, pass(this)	free (this)
	Destructor.
procedure, pass(this)	finalize (this)
	Finalize.
procedure, pass(this)	apply_scalar (this, x, n, time, strong)
	Apply the condition to a scalar field.
procedure, pass(this)	apply_vector (this, x, y, z, n, time, strong)
	(No-op) Apply vector.
procedure, pass(this)	apply_scalar_dev (this, x_d, time, strong, strm)
	Apply the condition to a scalar field (device).
procedure, pass(this)	apply_vector_dev (this, x_d, y_d, z_d, time, strong, strm)
	(No-op) Apply vector (device).
procedure, pass(this)	update (this, time, strm)
	Bring g up to date with the current time.
procedure, pass(this)	init_base (this, coef)
	Constructor.
procedure, pass(this)	free_base (this)
	Destructor.
procedure, pass(this)	mark_facet (this, facet, el)
	Mark a facet on an element as part of the boundary condition.
procedure, pass(this)	mark_facets (this, facet_list)
	Mark all facets from a (facet, element) tuple list.
procedure, pass(this)	mark_zone (this, bc_zone)
	Mark all facets from a zone.
procedure, pass(this)	mark_labeled_zone (this, zone_index)
	Mark all facets from a labeled zone.
procedure, pass(this)	mark_labeled_zones (this, zone_indices)
	Mark all facets from a list of labeled zones.
procedure, pass(this)	finalize_base (this)
	Finalize the construction of the bc by populating the msk and facet arrays.
procedure, pass(this)	apply_scalar_generic (this, x, time, strong, strm)
	Apply the boundary condition to a scalar field. Dispatches to the CPU or the device version.
procedure, pass(this)	apply_vector_generic (this, x, y, z, time, strong, strm)
	Apply the boundary condition to a vector field. Dispatches to the CPU or the device version.
procedure, pass(this)	debug_mask_ (this, file_name)
	Write a field showing the mask of the bcs.
procedure(bc_apply_scalar), deferred, pass	apply_scalar (this, x, n, time, strong)
	Apply the boundary condition to a scalar field on the CPU.
procedure(bc_apply_vector), deferred, pass	apply_vector (this, x, y, z, n, time, strong)
	Apply the boundary condition to a vector field on the CPU.
procedure(bc_apply_scalar_dev), deferred, pass	apply_scalar_dev (this, x_d, time, strong, strm)
	Device version of apply_scalar on the device.
procedure(bc_apply_vector_dev), deferred, pass	apply_vector_dev (this, x_d, y_d, z_d, time, strong, strm)
	Device version of apply_vector on the device.
procedure(bc_destructor), deferred, pass	free (this)
	Deferred destructor.
procedure(bc_constructor), deferred, pass	init (this, coef, json)
	Deferred constructor.
procedure(bc_finalize), deferred, pass	finalize (this)
	Deferred finalizer.

Public Attributes
type(expression_t)	expr
	The compiled expression prescribing the boundary value.
real(kind=rp), dimension(:), allocatable	g
	The value in each point of the mask.
real(kind=rp), dimension(:), allocatable	xm
	Coordinates of the points of the mask.
real(kind=rp), dimension(:), allocatable	ym
real(kind=rp), dimension(:), allocatable	zm
type(c_ptr)	g_d = C_NULL_PTR
	Device pointer for g.
integer, dimension(:), allocatable	msk
	The linear index of each constrained local degree of freedom.
integer, dimension(:), allocatable	facet_node_msk
	The linear index of each node on the marked boundary facets.
integer, dimension(:), allocatable	facet
	A list of facet ids (1 to 6), one for each entry in facet_node_msk.
type(dofmap_t), pointer	dof => null()
	Map of degrees of freedom.
type(coef_t), pointer	coef => null()
	SEM coefficients.
type(mesh_t), pointer	msh => null()
	The mesh.
type(space_t), pointer	xh => null()
	The function space.
type(stack_i4t2_t)	marked_facet
	Index tuples (facet, element) marked as part of the boundary condition.
type(c_ptr)	msk_d = C_NULL_PTR
	Device pointer for msk.
type(c_ptr)	facet_node_msk_d = C_NULL_PTR
	Device pointer for facet_node_msk.
type(c_ptr)	facet_d = C_NULL_PTR
	Device pointer for facet.
integer	bc_type = -1
	Type of the boundary condition, from the BC_* constants.
logical	updated = .false.
	Indicates wether the bc has been updated, for those BCs that need additional computations.
character(len=:), allocatable	name
integer, dimension(:), allocatable	zone_indices

Detailed Description

The expression is compiled once, at setup, and evaluated on the host in the points of the mask only. The resulting values are scattered into the solution with the same kernel that the Blasius condition uses, so no expression evaluation happens on the device.

An expression that does not depend on time is evaluated once, in finalize, and from then on the condition costs exactly what a constant Dirichlet condition costs. A time dependent one is re-evaluated once per timestep, gated on the updated flag of the base type.

Definition at line 61 of file expression_dirichlet.f90.

Member Function/Subroutine Documentation

apply_scalar() [1/2]

procedure(bc_apply_scalar), deferred, pass bc::bc_t::apply_scalar ( class(bc_t), intent(inout)  this, real(kind=rp), dimension(n), intent(inout)  x, integer, intent(in)  n, type(time_state_t), intent(in), optional  time, logical, intent(in), optional  strong  )

pure virtualinherited

Definition at line 132 of file bc.f90.

apply_scalar() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::apply_scalar	(	class(expression_dirichlet_t), intent(inout)	this,
		real(kind=rp), dimension(n), intent(inout)	x,
		integer, intent(in)	n,
		type(time_state_t), intent(in), optional	time,
		logical, intent(in), optional	strong
	)

Parameters

[in,out]	x	The field onto which to apply the values.
[in]	n	The size of x.
[in]	time	The current time state.
[in]	strong	Whether the condition is applied strongly.

Definition at line 83 of file expression_dirichlet.f90.

apply_scalar_dev() [1/2]

procedure(bc_apply_scalar_dev), deferred, pass bc::bc_t::apply_scalar_dev ( class(bc_t), intent(inout), target  this, type(c_ptr), intent(inout)  x_d, type(time_state_t), intent(in), optional  time, logical, intent(in), optional  strong, type(c_ptr), intent(inout)  strm  )

pure virtualinherited

Definition at line 136 of file bc.f90.

apply_scalar_dev() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::apply_scalar_dev	(	class(expression_dirichlet_t), intent(inout), target	this,
		type(c_ptr), intent(inout)	x_d,
		type(time_state_t), intent(in), optional	time,
		logical, intent(in), optional	strong,
		type(c_ptr), intent(inout)	strm
	)

Parameters

[in,out]	x_d	Device pointer to the field.
[in]	time	The current time state.
[in]	strong	Whether the condition is applied strongly.
[in,out]	strm	The device stream to issue the work on.

Definition at line 87 of file expression_dirichlet.f90.

apply_scalar_generic()

procedure, pass(this) bc::bc_t::apply_scalar_generic ( class(bc_t), intent(inout)  this, type(field_t), intent(inout)  x, type(time_state_t), intent(in), optional  time, logical, intent(in), optional  strong, type(c_ptr), intent(inout), optional  strm  )

inherited

Parameters

x	The x comp of the field for which to apply the bc.
time	Current time state.
strong	Whether we are setting a strong or a weak bc.
strm	Device stream

Definition at line 125 of file bc.f90.

apply_vector() [1/2]

procedure(bc_apply_vector), deferred, pass bc::bc_t::apply_vector ( class(bc_t), intent(inout)  this, real(kind=rp), dimension(n), intent(inout)  x, real(kind=rp), dimension(n), intent(inout)  y, real(kind=rp), dimension(n), intent(inout)  z, integer, intent(in)  n, type(time_state_t), intent(in), optional  time, logical, intent(in), optional  strong  )

pure virtualinherited

Definition at line 134 of file bc.f90.

apply_vector() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::apply_vector	(	class(expression_dirichlet_t), intent(inout)	this,
		real(kind=rp), dimension(n), intent(inout)	x,
		real(kind=rp), dimension(n), intent(inout)	y,
		real(kind=rp), dimension(n), intent(inout)	z,
		integer, intent(in)	n,
		type(time_state_t), intent(in), optional	time,
		logical, intent(in), optional	strong
	)

Parameters

[in,out]	x	The x-component of the field.
[in,out]	y	The y-component of the field.
[in,out]	z	The z-component of the field.
[in]	n	The size of x, y and z.
[in]	time	The current time state.
[in]	strong	Whether the condition is applied strongly.

Definition at line 85 of file expression_dirichlet.f90.

apply_vector_dev() [1/2]

procedure(bc_apply_vector_dev), deferred, pass bc::bc_t::apply_vector_dev ( class(bc_t), intent(inout), target  this, type(c_ptr), intent(inout)  x_d, type(c_ptr), intent(inout)  y_d, type(c_ptr), intent(inout)  z_d, type(time_state_t), intent(in), optional  time, logical, intent(in), optional  strong, type(c_ptr), intent(inout)  strm  )

pure virtualinherited

Definition at line 138 of file bc.f90.

apply_vector_dev() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::apply_vector_dev	(	class(expression_dirichlet_t), intent(inout), target	this,
		type(c_ptr), intent(inout)	x_d,
		type(c_ptr), intent(inout)	y_d,
		type(c_ptr), intent(inout)	z_d,
		type(time_state_t), intent(in), optional	time,
		logical, intent(in), optional	strong,
		type(c_ptr), intent(inout)	strm
	)

Parameters

[in,out]	x_d	Device pointer to the x-component of the field.
[in,out]	y_d	Device pointer to the y-component of the field.
[in,out]	z_d	Device pointer to the z-component of the field.
[in]	time	The current time state.
[in]	strong	Whether the condition is applied strongly.
[in,out]	strm	The device stream to issue the work on.

Definition at line 90 of file expression_dirichlet.f90.

apply_vector_generic()

procedure, pass(this) bc::bc_t::apply_vector_generic ( class(bc_t), intent(inout)  this, type(field_t), intent(inout)  x, type(field_t), intent(inout)  y, type(field_t), intent(inout)  z, type(time_state_t), intent(in), optional  time, logical, intent(in), optional  strong, type(c_ptr), intent(inout), optional  strm  )

inherited

Parameters

x	The x comp of the field for which to apply the bc.
y	The y comp of the field for which to apply the bc.
z	The z comp of the field for which to apply the bc.
time	Current time state.
strong	Whether we are setting a strong or a weak bc.
Device	stream

Definition at line 128 of file bc.f90.

debug_mask_()

procedure, pass(this) bc::bc_t::debug_mask_ ( class(bc_t), intent(inout)  this, character(len=*), intent(in)  file_name  )

inherited

The mask will be marked with 1.

Parameters

file_name

The name of the fld file.

Definition at line 130 of file bc.f90.

finalize() [1/2]

procedure(bc_finalize), deferred, pass bc::bc_t::finalize ( class(bc_t), intent(inout), target  this )

pure virtualinherited

Definition at line 144 of file bc.f90.

finalize() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::finalize

(

class(expression_dirichlet_t), intent(inout), target

this

)

Tabulates the coordinates of the points of the mask, which is only known once the mask has been built, and evaluates the expression right away if it does not depend on time.

Definition at line 81 of file expression_dirichlet.f90.

finalize_base()

procedure, pass(this) bc::bc_t::finalize_base ( class(bc_t), intent(inout), target  this )

inherited

This will linearize the marked facet's indicies in the msk array.

Definition at line 121 of file bc.f90.

free() [1/2]

procedure(bc_destructor), deferred, pass bc::bc_t::free ( class(bc_t), intent(inout), target  this )

pure virtualinherited

Definition at line 140 of file bc.f90.

free() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::free

(

class(expression_dirichlet_t), intent(inout), target

this

)

Definition at line 79 of file expression_dirichlet.f90.

free_base()

procedure, pass(this) bc::bc_t::free_base ( class(bc_t), intent(inout)  this )

inherited

Definition at line 108 of file bc.f90.

init() [1/2]

procedure(bc_constructor), deferred, pass bc::bc_t::init ( class(bc_t), intent(inout), target  this, type(coef_t), intent(in), target  coef, type(json_file), intent(inout)  json  )

pure virtualinherited

Definition at line 142 of file bc.f90.

init() [2/2]

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::init	(	class(expression_dirichlet_t), intent(inout), target	this,
		type(coef_t), intent(in), target	coef,
		type(json_file), intent(inout)	json
	)

Parameters

[in]	coef	The SEM coefficients.
[in,out]	json	The JSON object configuring the boundary condition.

Definition at line 74 of file expression_dirichlet.f90.

init_base()

procedure, pass(this) bc::bc_t::init_base ( class(bc_t), intent(inout)  this, type(coef_t), intent(in), target  coef  )

inherited

Parameters

dof	Map of degrees of freedom.

Definition at line 106 of file bc.f90.

init_from_components()

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::init_from_components	(	class(expression_dirichlet_t), intent(inout), target	this,
		type(coef_t), intent(in), target	coef,
		character(len=*), intent(in)	str
	)

Parameters

[in]	coef	The SEM coefficients.
[in]	str	The expression prescribing the boundary value.

Definition at line 76 of file expression_dirichlet.f90.

mark_facet()

procedure, pass(this) bc::bc_t::mark_facet ( class(bc_t), intent(inout)  this, integer, intent(in)  facet, integer, intent(in)  el  )

inherited

Parameters

facet	The index of the facet.
el	The index of the element.

Definition at line 110 of file bc.f90.

mark_facets()

procedure, pass(this) bc::bc_t::mark_facets ( class(bc_t), intent(inout)  this, type(stack_i4t2_t), intent(inout)  facet_list  )

inherited

Parameters

facet_list

The list of tuples.

Definition at line 112 of file bc.f90.

mark_labeled_zone()

procedure, pass(this) bc::bc_t::mark_labeled_zone ( class(bc_t), intent(inout)  this, integer, intent(in)  zone_index  )

inherited

Parameters

zone_index

The index of the labeled zone to be marked.

Definition at line 116 of file bc.f90.

mark_labeled_zones()

procedure, pass(this) bc::bc_t::mark_labeled_zones ( class(bc_t), intent(inout)  this, integer, dimension(:), intent(in)  zone_indices  )

inherited

Parameters

zone_indices

The indices of the labeled zones to be marked.

Definition at line 118 of file bc.f90.

mark_zone()

procedure, pass(this) bc::bc_t::mark_zone ( class(bc_t), intent(inout)  this, class(facet_zone_t), intent(in)  bc_zone  )

inherited

Parameters

bc_zone

Boundary zone to be marked.

Definition at line 114 of file bc.f90.

update()

procedure, pass(this) expression_dirichlet::expression_dirichlet_t::update	(	class(expression_dirichlet_t), intent(inout)	this,
		type(time_state_t), intent(in), optional	time,
		type(c_ptr), intent(inout), optional	strm
	)

A no-op unless the expression depends on time and the values have not already been updated in this timestep.

Parameters

[in]	time	The current time state.
[in,out]	strm	The device stream to issue the copy on.

Definition at line 93 of file expression_dirichlet.f90.

Here is the call graph for this function:

Member Data Documentation

bc_type

integer bc::bc_t::bc_type = -1

inherited

Definition at line 96 of file bc.f90.

coef

type(coef_t), pointer bc::bc_t::coef => null()

inherited

Definition at line 82 of file bc.f90.

dof

type(dofmap_t), pointer bc::bc_t::dof => null()

inherited

Definition at line 80 of file bc.f90.

expr

type(expression_t) expression_dirichlet::expression_dirichlet_t::expr

Definition at line 63 of file expression_dirichlet.f90.

facet

integer, dimension(:), allocatable bc::bc_t::facet

inherited

Definition at line 78 of file bc.f90.

facet_d

type(c_ptr) bc::bc_t::facet_d = C_NULL_PTR

inherited

Definition at line 94 of file bc.f90.

facet_node_msk

integer, dimension(:), allocatable bc::bc_t::facet_node_msk

inherited

Definition at line 76 of file bc.f90.

facet_node_msk_d

type(c_ptr) bc::bc_t::facet_node_msk_d = C_NULL_PTR

inherited

Definition at line 92 of file bc.f90.

g

real(kind=rp), dimension(:), allocatable expression_dirichlet::expression_dirichlet_t::g

Definition at line 65 of file expression_dirichlet.f90.

g_d

type(c_ptr) expression_dirichlet::expression_dirichlet_t::g_d = C_NULL_PTR

Definition at line 71 of file expression_dirichlet.f90.

marked_facet

type(stack_i4t2_t) bc::bc_t::marked_facet

inherited

Definition at line 88 of file bc.f90.

msh

type(mesh_t), pointer bc::bc_t::msh => null()

inherited

Definition at line 84 of file bc.f90.

msk

integer, dimension(:), allocatable bc::bc_t::msk

inherited

Definition at line 74 of file bc.f90.

msk_d

type(c_ptr) bc::bc_t::msk_d = C_NULL_PTR

inherited

Definition at line 90 of file bc.f90.

name

character(len=:), allocatable bc::bc_t::name

inherited

Definition at line 101 of file bc.f90.

updated

logical bc::bc_t::updated = .false.

inherited

Definition at line 99 of file bc.f90.

xh

type(space_t), pointer bc::bc_t::xh => null()

inherited

Definition at line 86 of file bc.f90.

xm

real(kind=rp), dimension(:), allocatable expression_dirichlet::expression_dirichlet_t::xm

Definition at line 67 of file expression_dirichlet.f90.

ym

real(kind=rp), dimension(:), allocatable expression_dirichlet::expression_dirichlet_t::ym

Definition at line 68 of file expression_dirichlet.f90.

zm

real(kind=rp), dimension(:), allocatable expression_dirichlet::expression_dirichlet_t::zm

Definition at line 69 of file expression_dirichlet.f90.

zone_indices

integer, dimension(:), allocatable bc::bc_t::zone_indices

inherited

Definition at line 103 of file bc.f90.

---

The documentation for this type was generated from the following file:

• /home/runner/work/neko/neko/src/bc/expression_dirichlet.f90