modified PVfunctions.cpp, considering massless limits, etc. carefully.