cdef extern from "math.h": cpdef double pxd_sqrt "sqrt"(double)