Procedures

Allocate memory for a new block and return a pointer to a new m_allocator object.

Constructor function for the cuda_tdsops_t class. See tdsops_t for details.

Curl of a vector field (u, v, w).

Wrapper for curl

Go through the block list from head to tail, deallocating each memory block in turn. Deallocation of a m_allocator object automatically deallocates its internal allocatable data array.

Divergence of a vector field (u, v, w).

Wrapper for divergence_v2p

Supports 1D, 2D, and 3D domain decomposition.

Return a pointer to the first available memory block, i.e. the current head of the block list. If the list is empty, allocate a new block with create_block first.

Utility function that returns a array made of the id of the block currently in the block list. Return the array [0] if block list is empty.

Get the physical location of a cell center with i,j,k local indices

Getter for local domain dimensions

Getter for domain dimensions

Extract data from field f optionally reordering into dir orientation. To output in same orientation as f, use call ...%get_field_data(data, f, f%dir)

Getter for the dimensions of an array directed along dir where data would be located on data_loc

Getter for the dimensions of field phi

Getter for the dimensions of field phi where data is located on data_loc

Getter for local domain dimensions

Get application storage directional index from cartesian index

Get cartesian index from application storage directional one

Converts a set of application storage directional index to an other direction. The two directions are defined by the reorder_dir variable, RDR_X2Y will go from storage in X to Y etc.

Getter for the main dimension a field oriented along dir with data on data_loc

Getter for the number of groups for fields in direction dir

Getter for the number of groups for fields phi

Getter for the main dimension of field phi

Getter for padded dimensions with structure in dir direction

Getter for padded dimensions for field phi Gets the field direction from the field itself

Returns RDR_?2? value based on two direction inputs

Getter for parameter SZ

Get the tds_n size based on the from_to value (and the mesh)

Gradient of a scalar field 'p'.

Wrapper for gradient_p2v

Returns whether or not the current rank is the root rank

Laplacian of a scalar field 'u'.

Completely initialise the mesh object. Upon initialisation the mesh object can be read-only and shouldn't be edited Takes as argument global information about the mesh like its length, number of cells and decomposition in each direction

Reports the div(u) at cell centres

Reports the enstrophy

Post-processes the divergence of velocity in spectral space, including scaling w.r.t. grid size.

Release memory block pointed to by HANDLE to the block list. It is pushed to the front of the block list, in other words it is made the head block.

Runs the solver forwards in time from t=t_0 to t=T, performing postprocessing/IO and reporting diagnostics.

Constructor function for the tdsops_t class.

Skew-symmetric form of convection-diffusion terms in the incompressible Navier-Stokes momemtum equations, excluding pressure terms. Inputs from velocity grid and outputs to velocity grid.

Thomas algorithm implementation. So much more easier than the distributed algorithm. It is intended to work only on a single rank so there is no MPI communication.

Spectral equivalence constants