def get_drdX(self,X): "Derivative of radial coordinate with respect to compactified" dXdr = self.get_dXdr(X) with np.errstate(invalid='ignore'): drdX = 1./dXdr return drdX