Dot product of two 3-vectors (typically two VSH values evaluated at
the same ): . Note this is a plain component-wise
product, not a Hermitian inner product -- callers wanting
must conjugate VSH2
themselves first (as done throughout src/tests.f90 and the
examples/ drivers).
Returns: , complex scalar.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| complex(kind=dp), | intent(in), | DIMENSION(3) | :: | VSH1 |
First complex 3-vector . |
|
| complex(kind=dp), | intent(in), | DIMENSION(3) | :: | VSH2 |
Second complex 3-vector . |
FUNCTION DOT_DP(VSH1, VSH2) RESULT(DOT) IMPLICIT NONE COMPLEX(KIND=dp), DIMENSION(3), INTENT(IN) :: VSH1 COMPLEX(KIND=dp), DIMENSION(3), INTENT(IN) :: VSH2 COMPLEX(KIND=dp) :: DOT DOT = VSH1(1)*VSH2(1) + VSH1(2)*VSH2(2) + VSH1(3)*VSH2(3) RETURN END FUNCTION DOT_DP