64 use math,
only:
rzero,
rone,
copy,
cmult,
cadd,
cfill,
invcol1,
vdot3,
add2, &
88 integer,
intent(in),
optional :: n
89 type(
field_t),
intent(inout) :: a
101 call rzero(a%x, size)
107 integer,
intent(in),
optional :: n
108 type(
field_t),
intent(inout) :: a
126 integer,
intent(in),
optional :: n
128 type(
field_t),
intent(inout) :: a
140 call copy(a%x, b%x, size)
146 integer,
intent(in),
optional :: n
147 type(
field_t),
intent(inout) :: a
148 real(kind=
rp),
intent(in) :: c
160 call cmult(a%x, c, size)
166 integer,
intent(in),
optional :: n
167 type(
field_t),
intent(inout) :: a
168 real(kind=
rp),
intent(in) :: s
180 call cadd(a%x, s, size)
186 integer,
intent(in),
optional :: n
187 type(
field_t),
intent(inout) :: a
188 real(kind=
rp),
intent(in) :: c
200 call cfill(a%x, c, size)
206 integer,
intent(in),
optional :: n
207 type(
field_t),
intent(inout) :: a
227 integer,
intent(in),
optional :: n
228 type(
field_t),
intent(in) :: u1, u2, u3
229 type(
field_t),
intent(in) :: v1, v2, v3
230 type(
field_t),
intent(out) :: dot
241 u1%x_d, u2%x_d, u3%x_d, &
242 v1%x_d, v2%x_d, v3%x_d, &
255 integer,
intent(in),
optional :: n
256 type(
field_t),
intent(inout) :: a
269 call add2(a%x, b%x, size)
276 integer,
intent(in),
optional :: n
277 type(
field_t),
intent(inout) :: a
278 type(
field_t),
intent(inout) :: b
290 call sub2(a%x, b%x, size)
297 integer,
intent(in),
optional :: n
298 type(
field_t),
intent(inout) :: c
299 type(
field_t),
intent(inout) :: b
300 type(
field_t),
intent(out) :: a
312 call sub3(a%x, b%x, c%x, size)
321 integer,
intent(in),
optional :: n
322 type(
field_t),
intent(inout) :: a
323 type(
field_t),
intent(inout) :: b
324 real(kind=
rp),
intent(in) :: c1
336 call add2s1(a%x, b%x, c1, size)
344 integer,
intent(in),
optional :: n
345 type(
field_t),
intent(inout) :: a
346 type(
field_t),
intent(inout) :: b
347 real(kind=
rp),
intent(in) :: c1
359 call add2s2(a%x, b%x, c1, size)
366 integer,
intent(in),
optional :: n
367 type(
field_t),
intent(inout) :: a
369 real(kind=
rp),
intent(in) :: c1
388 integer,
intent(in),
optional :: n
389 type(
field_t),
intent(inout) :: a
391 real(kind=
rp),
intent(in) :: c
403 call cmult2(a%x, b%x, c, size)
410 integer,
intent(in),
optional :: n
411 type(
field_t),
intent(inout) :: a
432 integer,
intent(in),
optional :: n
433 type(
field_t),
intent(inout) :: a
446 call col2(a%x, b%x, size)
453 integer,
intent(in),
optional :: n
454 type(
field_t),
intent(inout) :: a
468 call col3(a%x, b%x, c%x, size)
475 integer,
intent(in),
optional :: n
476 type(
field_t),
intent(inout) :: a
490 call subcol3(a%x, b%x, c%x, size)
497 integer,
intent(in),
optional :: n
498 type(
field_t),
intent(inout) :: a
501 real(kind=
rp),
intent(in) :: c1, c2
513 call add3s2(a%x, b%x, c%x, c1, c2, size)
520 integer,
intent(in),
optional :: n
521 type(
field_t),
intent(inout) :: a
535 call addcol3(a%x, b%x, c%x, size)
542 integer,
intent(in),
optional :: n
543 type(
field_t),
intent(inout) :: a
558 call addcol4(a%x, b%x, c%x, d%x, size)
564 integer,
intent(in),
optional :: n
578 sum =
glsum(a%x, size)
584 integer,
intent(in),
optional :: n
585 type(
field_t),
intent(in) :: a, b
586 real(kind=
rp) :: norm
598 norm =
glsc2(a%x, b%x, size)
604 integer,
intent(in),
optional :: n
605 type(
field_t),
intent(in) :: a, b, c
606 real(kind=
rp) :: norm
618 norm =
glsc3(a%x, b%x, c%x, size)
subroutine, public device_add2(a_d, b_d, n)
Vector addition .
subroutine, public device_addcol3(a_d, b_d, c_d, n)
Returns .
subroutine, public device_col2(a_d, b_d, n)
Vector multiplication .
subroutine, public device_add2s1(a_d, b_d, c1, n)
subroutine, public device_rzero(a_d, n)
Zero a real vector.
subroutine, public device_rone(a_d, n)
Set all elements to one.
subroutine, public device_add2s2(a_d, b_d, c1, n)
Vector addition with scalar multiplication (multiplication on first argument)
subroutine, public device_invcol1(a_d, n)
Invert a vector .
subroutine, public device_col3(a_d, b_d, c_d, n)
Vector multiplication with 3 vectors .
subroutine, public device_cadd(a_d, c, n)
Add a scalar to vector .
subroutine, public device_vdot3(dot_d, u1_d, u2_d, u3_d, v1_d, v2_d, v3_d, n)
Compute a dot product (3-d version) assuming vector components etc.
subroutine, public device_cmult2(a_d, b_d, c, n)
Multiplication by constant c .
subroutine, public device_cmult(a_d, c, n)
Multiplication by constant c .
real(kind=rp) function, public device_glsc2(a_d, b_d, n)
Weighted inner product .
subroutine, public device_sub3(a_d, b_d, c_d, n)
Vector subtraction .
real(kind=rp) function, public device_glsc3(a_d, b_d, c_d, n)
Weighted inner product .
real(kind=rp) function, public device_glsum(a_d, n)
Sum a vector of length n.
subroutine, public device_copy(a_d, b_d, n)
Copy a vector .
subroutine, public device_add3s2(a_d, b_d, c_d, c1, c2, n)
Returns .
subroutine, public device_subcol3(a_d, b_d, c_d, n)
Returns .
subroutine, public device_sub2(a_d, b_d, n)
Vector substraction .
subroutine, public device_cfill(a_d, c, n)
Set all elements to a constant c .
subroutine, public device_addcol4(a_d, b_d, c_d, d_d, n)
Returns .
subroutine, public device_invcol2(a_d, b_d, n)
Vector division .
subroutine, public device_addsqr2s2(a_d, b_d, c1, n)
Returns .
subroutine, public field_add2s1(a, b, c1, n)
Vector addition with scalar multiplication (multiplication on first argument)
subroutine, public field_invcol2(a, b, n)
Vector division .
subroutine, public field_cadd(a, s, n)
Add a scalar to vector .
subroutine, public field_invcol1(a, n)
Invert a vector .
subroutine, public field_col2(a, b, n)
Vector multiplication .
subroutine, public field_sub2(a, b, n)
Vector substraction .
subroutine, public field_sub3(a, b, c, n)
Vector subtraction .
subroutine, public field_cmult2(a, b, c, n)
Multiplication by constant c .
real(kind=rp) function, public field_glsc2(a, b, n)
subroutine, public field_cfill(a, c, n)
Set all elements to a constant c .
subroutine, public field_add2s2(a, b, c1, n)
Vector addition with scalar multiplication (multiplication on second argument)
real(kind=rp) function, public field_glsum(a, n)
subroutine, public field_rzero(a, n)
Zero a real vector.
subroutine, public field_addcol3(a, b, c, n)
Returns .
subroutine, public field_rone(a, n)
Set all elements to one.
subroutine, public field_add3s2(a, b, c, c1, c2, n)
Returns .
subroutine, public field_add2(a, b, n)
Vector addition .
subroutine, public field_col3(a, b, c, n)
Vector multiplication with 3 vectors .
subroutine, public field_copy(a, b, n)
Copy a vector .
subroutine, public field_subcol3(a, b, c, n)
Returns .
subroutine, public field_cmult(a, c, n)
Multiplication by constant c .
real(kind=rp) function, public field_glsc3(a, b, c, n)
subroutine, public field_addcol4(a, b, c, d, n)
Returns .
subroutine, public field_addsqr2s2(a, b, c1, n)
Returns .
subroutine, public field_vdot3(dot, u1, u2, u3, v1, v2, v3, n)
Compute a dot product (3-d version) assuming vector components etc.
subroutine, public cmult(a, c, n)
Multiplication by constant c .
subroutine, public cmult2(a, b, c, n)
Multiplication by constant c .
subroutine, public invcol2(a, b, n)
Vector division .
real(kind=rp) function, public glsc3(a, b, c, n)
Weighted inner product .
subroutine, public cadd(a, s, n)
Add a scalar to vector .
subroutine, public addsqr2s2(a, b, c1, n)
Returns .
subroutine, public add2s1(a, b, c1, n)
Vector addition with scalar multiplication (multiplication on first argument)
real(kind=rp) function, public glsc2(a, b, n)
Weighted inner product .
subroutine, public subcol3(a, b, c, n)
Returns .
subroutine, public rone(a, n)
Set all elements to one.
real(kind=rp) function, public glsum(a, n)
Sum a vector of length n.
subroutine, public sub3(a, b, c, n)
Vector subtraction .
subroutine, public addcol4(a, b, c, d, n)
Returns .
subroutine, public add2(a, b, n)
Vector addition .
subroutine, public cfill(a, c, n)
Set all elements to a constant c .
subroutine, public add3s2(a, b, c, c1, c2, n)
Returns .
subroutine, public addcol3(a, b, c, n)
Returns .
subroutine, public invcol1(a, n)
Invert a vector .
subroutine, public col2(a, b, n)
Vector multiplication .
subroutine, public copy(a, b, n)
Copy a vector .
subroutine, public col3(a, b, c, n)
Vector multiplication with 3 vectors .
subroutine, public vdot3(dot, u1, u2, u3, v1, v2, v3, n)
Compute a dot product (3-d version) assuming vector components etc.
subroutine, public rzero(a, n)
Zero a real vector.
subroutine, public sub2(a, b, n)
Vector substraction .
subroutine, public add2s2(a, b, c1, n)
Vector addition with scalar multiplication (multiplication on second argument)
integer, parameter neko_bcknd_device
integer, parameter, public rp
Global precision used in computations.