m_omptgt_kernels_fieldops Module

Arguments

Return Value real(kind=dp)

Scalar product of two DIR_C fields, whose storage is the padded Cartesian box n_c with the physical points in 1:dims.

Arguments

Return Value real(kind=dp)

Arguments

Return Value real(kind=dp)

Arguments

Return Value real(kind=dp)

Arguments

Arguments

Arguments

Arguments

Arguments

Arguments

Arguments

Arguments

Arguments

Copies the first and last four points of every line into the start and end send buffers. n is the line length, n_u the padded extent.

Arguments

s(1:3, l) += sums of ga, gb and g*c over the points labelled l, as cuda/kernels/fieldops.f90: each thread accumulates runs of equal labels along its pencil and flushes a run when the label changes, so the atomics scale with the label changes rather than the points. Atomics instead of reduction(+:s), which flang cannot do on the GPU for an array whose shape is only known at run time.

Arguments

a, b, c += coeffs(1:3, l)*g at the points labelled l.

Arguments

Maximum and sum of |f|.

Arguments

Signed maximum and sum of the slice f(:, i_slice, :).

Arguments

Sum of the slice f(:, i_slice, :) with the weight w(j) on the points of stacked index j.

Arguments

Sets the x faces: a Dirichlet value c_start at the inlet (pencil index 1) and the convective outflow du/dt + Uc*du/dx = 0, Uc = c_end, at the outlet (pencil index nx). BC_NEUMANN (1) leaves a face untouched.

Arguments

As field_set_x_face, with the inlet taken from the inlet plane of f_start rather than a constant.

Arguments

Add the constant c to every point of the x plane i_plane of f.

Arguments

Writes two (ny, nz) planes into the x = 1 and x = nx planes of f. p_end is read only when apply_end is true.

Arguments

f += g on the x faces flagged by set_start/set_end.

Arguments

Sets the bottom (y = 1) and top (y = ny) faces where flagged, so the two y faces can differ (ABL: no-slip floor, free-slip lid).

Arguments

Copies the bottom and top y planes of f_start onto those of f.

Arguments

Writes two (nx, nz) planes into the y = 1 and y = ny planes of f. p_end is read only when apply_end is true.

Arguments

f += g on the y faces flagged by set_start/set_end.

Arguments

Sets the y plane plane to c.

Arguments

Drag law tau = drag_coeffu_c|u_h| from each column's own sampled velocity, written into one stress plane (component 1: x, 3: z). The velocity is interpolated between sample_plane and upper_plane.

Arguments

band(i, j, k) = f(i_start + i - 1, j, k) for a DIR_X vertex field

Arguments

f(i_start + i - 1, j, k) = lambda(i)band(i, j, k) + (1 - lambda(i))f(i_start + i - 1, j, k) for a DIR_X vertex field

Arguments

First pass of field_plane_sums: partial(j, z) is the sum over x of y index j at z, summed in a fixed order.

Arguments

Second pass of field_plane_sums: sums(j) is the sum over z of partial(j, z), in a fixed order.

Arguments

Reorders u_i in layout dir_from into u_o in layout dir_to.

Arguments

u += u_d, where u is DIR_X and u_d the same field in layout dir_from.

Arguments

vals(1:3, p) = this rank's part of the trilinear interpolation of the DIR_X VERT fields a, b and c at points(:, p), one thread per point; see field_sample_points. grid holds the local VERT dims, global VERT dims, rank offsets and periodic flags (1/0), three entries each, then the padded y-block count of the DIR_X layout.

Arguments

a, b, c += src(4:6, s)*exp(-dist^2/eps^2) within r_cut of the point src(1:3, s), over the global index box box(1:3, s)..box(4:6, s); see source_index_boxes. Box indices are wrapped across periodic boundaries and skipped where this rank does not own them. One thread per (source, z, y) of the box, each sweeping x; sources overlap, so the adds are atomic. grid holds the local VERT dims, global VERT dims, rank offsets and periodic flags (1/0), three entries each, then the padded y-block count.

Arguments