implement O(d^2) dawson F function, and de-clutter test_utpm.py by moving scalar function identities to a separate module