73 integer,
intent(in) :: n, gdim
74 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
93 subroutine opcolv(a1, a2, a3, c, gdim, n)
94 integer,
intent(in) :: n, gdim
95 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
96 real(kind=
rp),
dimension(n),
intent(in) :: c
115 subroutine opcolv3c(a1, a2, a3, b1, b2, b3, c, d, n, gdim)
116 integer,
intent(in) :: n, gdim
117 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
118 real(kind=
rp),
dimension(n),
intent(in) :: b1, b2, b3
119 real(kind=
rp),
intent(in) :: c(n), d
122 if (gdim .eq. 3)
then
138 subroutine opadd2cm(a1, a2, a3, b1, b2, b3, c, n, gdim)
139 integer,
intent(in) :: n, gdim
140 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
141 real(kind=
rp),
dimension(n),
intent(in) :: b1, b2, b3
142 real(kind=
rp),
intent(in) :: c
145 if (gdim .eq. 3)
then
147 a1(i) = a1(i) + b1(i)*c
148 a2(i) = a2(i) + b2(i)*c
149 a3(i) = a3(i) + b3(i)*c
153 a1(i) = a1(i) + b1(i)*c
154 a2(i) = a2(i) + b2(i)*c
161 subroutine opadd2col(a1, a2, a3, b1, b2, b3, c, n, gdim)
162 integer,
intent(in) :: n, gdim
163 real(kind=
rp),
dimension(n),
intent(inout) :: a1, a2, a3
164 real(kind=
rp),
dimension(n),
intent(in) :: b1, b2, b3
165 real(kind=
rp),
intent(in) :: c(n)
168 if (gdim .eq. 3)
then
170 a1(i) = a1(i) + b1(i)*c(i)
171 a2(i) = a2(i) + b2(i)*c(i)
172 a3(i) = a3(i) + b3(i)*c(i)
176 a1(i) = a1(i) + b1(i)*c(i)
177 a2(i) = a2(i) + b2(i)*c(i)
Collection of vector field operations operating on and . Note that in general the indices and ....
subroutine opcolv3c(a1, a2, a3, b1, b2, b3, c, d, n, gdim)
subroutine opchsign(a1, a2, a3, gdim, n)
for and .
subroutine opcolv(a1, a2, a3, c, gdim, n)
subroutine opadd2col(a1, a2, a3, b1, b2, b3, c, n, gdim)
subroutine opadd2cm(a1, a2, a3, b1, b2, b3, c, n, gdim)
integer, parameter, public rp
Global precision used in computations.