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