76 integer,
intent(in) :: n, gdim
77 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
96 subroutine opcolv(a1, a2, a3, c, gdim, n)
97 integer,
intent(in) :: n, gdim
98 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
99 real(kind=
rp),
dimension(n),
intent(in) :: c
102 if (gdim .eq. 3)
then
118 subroutine opcolv3c(a1, a2, a3, b1, b2, b3, c, d, n, gdim)
119 integer,
intent(in) :: n, gdim
120 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
121 real(kind=
rp),
dimension(n),
intent(in) :: b1, b2, b3
122 real(kind=
rp),
intent(in) :: c(n), d
125 if (gdim .eq. 3)
then
141 subroutine opadd2cm(a1, a2, a3, b1, b2, b3, c, n, gdim)
142 integer,
intent(in) :: n, gdim
143 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
144 real(kind=
rp),
dimension(n),
intent(in) :: b1, b2, b3
145 real(kind=
rp),
intent(in) :: c
148 if (gdim .eq. 3)
then
150 a1(i) = a1(i) + b1(i)*c
151 a2(i) = a2(i) + b2(i)*c
152 a3(i) = a3(i) + b3(i)*c
156 a1(i) = a1(i) + b1(i)*c
157 a2(i) = a2(i) + b2(i)*c
164 subroutine opadd2col(a1, a2, a3, b1, b2, b3, c, n, gdim)
165 integer,
intent(in) :: n, gdim
166 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
167 real(kind=
rp),
dimension(n),
intent(in) :: b1, b2, b3
168 real(kind=
rp),
intent(in) :: c(n)
171 if (gdim .eq. 3)
then
173 a1(i) = a1(i) + b1(i)*c(i)
174 a2(i) = a2(i) + b2(i)*c(i)
175 a3(i) = a3(i) + b3(i)*c(i)
179 a1(i) = a1(i) + b1(i)*c(i)
180 a2(i) = a2(i) + b2(i)*c(i)
Collection of vector field operations operating on and . Note that in general the indices and ....
subroutine, public opcolv(a1, a2, a3, c, gdim, n)
subroutine, public opadd2col(a1, a2, a3, b1, b2, b3, c, n, gdim)
subroutine, public opchsign(a1, a2, a3, gdim, n)
for and .
subroutine, public opadd2cm(a1, a2, a3, b1, b2, b3, c, n, gdim)
subroutine, public opcolv3c(a1, a2, a3, b1, b2, b3, c, d, n, gdim)
integer, parameter, public rp
Global precision used in computations.