ax_helm_compr_cpu Module Reference

CPU Helmholtz matrix-vector product on compressed geometric factors. More...

Data Types
type	ax_helm_compr_cpu_t
	CPU matrix-vector product for a Helmholtz problem, reading the geometric factors \( G_{ij} \) from their compressed form in coef_t. More...

Functions/Subroutines
subroutine	ax_helm_compr_compute (this, w, u, coef, msh, xh)
	Compute \( Ax \) product, taking three components of a vector field in an uncoupled manner.
subroutine	ax_helm_compr_lx (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m, lx)
	Generic CPU kernel for the Helmholz matrix-vector product.
subroutine	ax_helm_compr_lx14 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx13 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx12 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx11 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx10 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx9 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx8 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx7 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx6 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx5 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx4 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx3 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)
subroutine	ax_helm_compr_lx2 (w, u, dx, dy, dz, dxt, dyt, dzt, h1, g11, g22, g33, g12, g13, g23, ind, n, m)

Detailed Description

Note

The kernels are those of ax_helm_cpu with the geometric factors read at G(:,:,:,ind(e)), so a change to one belongs in the other.

Function/Subroutine Documentation

ax_helm_compr_compute()

subroutine ax_helm_compr_cpu::ax_helm_compr_compute	(	class(ax_helm_compr_cpu_t), intent(in)	this,
		real(kind=rp), dimension(xh%lx, xh%ly, xh%lz, msh%nelv), intent(inout)	w,
		real(kind=rp), dimension(xh%lx, xh%ly, xh%lz, msh%nelv), intent(in)	u,
		type(coef_t), intent(in)	coef,
		type(mesh_t), intent(in)	msh,
		type(space_t), intent(in)	xh
	)

Parameters

au	Result for the first component of the vector.
av	Result for the first component of the vector.
aw	Result for the first component of the vector.
u	The first component of the vector.
v	The second component of the vector.
w	The third component of the vector.
coef	Coefficients.
msh	Mesh.
Xh	Function space \( X_h \). Compute the product.
w	Vector of size (lx,ly,lz,nelv).
u	Vector of size (lx,ly,lz,nelv).
coef	Coefficients.
msh	Mesh.
Xh	Function space \( X_h \).

Note

Since this is a performance-crtical routine, it is implemented in several kernels corresponding to different polynmial orders.

Definition at line 101 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m, integer, intent(in)  lx  )

private

Parameters

w	Result.
u	Velocity field.
Dx	Derivative operator in first dimension.
Dy	Derivative operator in second dimension.
Dz	Derivative operator in third dimension.
Dxt	Derivative operator transpose in first dimension.
Dyt	Derivative operator transpose in second dimension.
Dzt	Derivative operator transpose in third dimension.
G11	Compressed geometric factor.
ind	Index of the compressed geometric factors of each element.
n	Number of elements.
m	Number of compressed geometric factors.
lx	Polynomial order.

Definition at line 248 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx10()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx10 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1007 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx11()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx11 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 855 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx12()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx12 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 697 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx13()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx13 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 532 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx14()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx14 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 362 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx2()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx2 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 2007 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx3()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx3 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1903 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx4()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx4 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1793 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx5()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx5 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1677 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx6()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx6 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1555 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx7()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx7 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1427 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx8()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx8 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1293 of file ax_helm_compr_cpu.f90.

ax_helm_compr_lx9()

subroutine ax_helm_compr_cpu::ax_helm_compr_lx9 ( real(kind=rp), dimension(lx, lx, lx, n), intent(inout)  w, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  u, real(kind=rp), dimension(lx, lx), intent(in)  dx, real(kind=rp), dimension(lx, lx), intent(in)  dy, real(kind=rp), dimension(lx, lx), intent(in)  dz, real(kind=rp), dimension(lx, lx), intent(in)  dxt, real(kind=rp), dimension(lx, lx), intent(in)  dyt, real(kind=rp), dimension(lx, lx), intent(in)  dzt, real(kind=rp), dimension(lx, lx, lx, n), intent(in)  h1, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g11, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g22, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g33, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g12, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g13, real(kind=rp), dimension(lx, lx, lx, m), intent(in)  g23, integer, dimension(n), intent(in)  ind, integer, intent(in)  n, integer, intent(in)  m  )

private

Definition at line 1153 of file ax_helm_compr_cpu.f90.