64 type(
space_t),
intent(inout) :: Xh
65 type(
mesh_t),
intent(inout) :: msh
66 type(
coef_t),
intent(inout) :: coef
67 real(kind=
rp),
intent(inout) :: au(xh%lx, xh%ly, xh%lz, msh%nelv)
68 real(kind=
rp),
intent(inout) :: av(xh%lx, xh%ly, xh%lz, msh%nelv)
69 real(kind=
rp),
intent(inout) :: aw(xh%lx, xh%ly, xh%lz, msh%nelv)
70 real(kind=
rp),
intent(inout) :: u(xh%lx, xh%ly, xh%lz, msh%nelv)
71 real(kind=
rp),
intent(inout) :: v(xh%lx, xh%ly, xh%lz, msh%nelv)
72 real(kind=
rp),
intent(inout) :: w(xh%lx, xh%ly, xh%lz, msh%nelv)
74 call this%compute(au, u, coef, msh, xh)
75 call this%compute(av, v, coef, msh, xh)
76 call this%compute(aw, w, coef, msh, xh)
subroutine ax_helm_compute_vector(this, au, av, aw, u, v, w, coef, msh, Xh)
Defines a Matrix-vector product.
subroutine, public addcol4(a, b, c, d, n)
Returns .
integer, parameter, public rp
Global precision used in computations.
Defines a function space.
Matrix-vector product for a Helmholtz problem.
Base type for a matrix-vector product providing .
Coefficients defined on a given (mesh, ) tuple. Arrays use indices (i,j,k,e): element e,...
The function space for the SEM solution fields.