Is there any way to get hold of a long double pi to the maximum possible precision? The following only gets one at the precision of double:

p = np.longdouble(np.pi)

One option would be to change the definitions of np.pi to longdouble, but then that contaminates all code that uses them into also using longdouble.