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