test_flume_weno_interpolation_fnl_kernels

Unit test RV-1 (device half) of the WENO interpolation tables: the device WENO primitive on one packed stencil.

Module-level (not a contains-internal procedure): device routines must be module procedures.

Source: src/tests/flume/unit/test_flume_weno_interpolation_fnl.F90

Dependencies

Contents

• interpolate
• interpolate_width

Subroutines

interpolate

Interpolate one packed stencil with the device WENO primitive.

subroutine interpolate(S, a, p, d, zeps, sigma, v, vr)

Arguments

Name	Type	Intent	Attributes	Description
S	integer(kind=I4P)	in		Number of stencils used.
a	real(kind=R8P)	in		Interpolation optimal weights.
p	real(kind=R8P)	in		Interpolation polynomials coefficients.
d	real(kind=R8P)	in		Smoothness indicators coefficients.
zeps	real(kind=R8P)	in		Parameter avoiding division by zero.
sigma	real(kind=R8P)	in		Descaler switch: 0 Jiang-Shu, 1 scale-invariant.
v	real(kind=R8P)	in		Packed stencil.
vr	real(kind=R8P)	out		Interface values.

Call graph

interpolate_width

Interpolate one packed stencil with the generic and with the per-width device WENO routines, plain and with the weights ratio: vr(:,1) generic, vr(:,2) per-width, vr(:,3) generic with ratio, vr(:,4) per-width with ratio, and wr(:,1), wr(:,2) the ratio of the last two. A width without a routine (S > 5) leaves NaN in the per-width outputs.

subroutine interpolate_width(S, a, p, d, zeps, sigma, v, vr, wr)

Arguments

Name	Type	Intent	Attributes	Description
S	integer(kind=I4P)	in		Number of stencils used.
a	real(kind=R8P)	in		Interpolation optimal weights.
p	real(kind=R8P)	in		Interpolation polynomials coefficients.
d	real(kind=R8P)	in		Smoothness indicators coefficients.
zeps	real(kind=R8P)	in		Parameter avoiding division by zero.
sigma	real(kind=R8P)	in		Descaler switch: 0 Jiang-Shu, 1 scale-invariant.
v	real(kind=R8P)	in		Packed stencil.
vr	real(kind=R8P)	out		Interface values, see above.
wr	real(kind=R8P)	out		Weights ratios, generic and per-width.

Call graph