domain.py 文件源码

def get_d2rdX2(self,X):
        "Second derivative of radial coordinate with respect to compactified"
        d2Xdr2 = self.get_d2Xdr2(X)
        with np.errstate(invalid='ignore'):
            d2rdX2 = 1./d2Xdr2
        return d2rdX2