60    def hessSdf(xc, yc, zc, x, y, z):
   61        x = x-xc
   62        y = y-yc
   63        z = z-zc
   64        Hxx = -x**2/(x**2 + y**2 + z**2)**(3/2) + 1/np.sqrt(x**2 + y**2 + z**2)
   65        Hzx = -x*z/(x**2 + y**2 + z**2)**(3/2)
   66        Hxy = -x*y/(x**2 + y**2 + z**2)**(3/2)
   67        Hyy = -y**2/(x**2 + y**2 + z**2)**(3/2) + 1/np.sqrt(x**2 + y**2 + z**2)
   68        Hzy = -y*z/(x**2 + y**2 + z**2)**(3/2)
   69        Hzz = -z**2/(x**2 + y**2 + z**2)**(3/2) + 1/np.sqrt(x**2 + y**2 + z**2)
   70        
   71        Hxx = Hxx.reshape((-1,1))
   72        Hzx = Hzx.reshape((-1,1))
   73        Hxy = Hxy.reshape((-1,1))
   74        Hyy = Hyy.reshape((-1,1))
   75        Hzy = Hzy.reshape((-1,1))
   76        Hzz = Hzz.reshape((-1,1))
   77        hess_array = np.hstack([Hxx, Hxy, Hzx, Hyy, Hzy, Hzz])
   78 
   79        return hess_array
   80 
   81