dot product of two vectors.
JW: rewrote this routine.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | n |
vector dimension. |
||
| real(kind=wp), | intent(in) | :: | x(*) |
x(n) input vector. |
||
| real(kind=wp), | intent(in) | :: | y(*) |
y(n) input vector. |
value of dot product mxvdot=trans(x)*y.
pure function mxvdot(n,x,y) integer,intent(in) :: n !!vector dimension. real(wp),intent(in) :: x(*) !! x(n) input vector. real(wp),intent(in) :: y(*) !! y(n) input vector. real(wp) :: mxvdot !! value of dot product `mxvdot=trans(x)*y`. mxvdot = dot_product(x(1:n),y(1:n)) end function mxvdot