Description: Inner (=scalar=dot) product of two pv-vectors.
The source code for the Fortran routine PVDPV is available: