sdf_sphere.Sphere Class Reference

Public Member Functions
	sDF (r, xc, yc, zc, x, y, z)
	gradSdf (xc, yc, zc, x, y, z)
	hessSdf (xc, yc, zc, x, y, z)

Detailed Description

Definition at line 20 of file sdf_sphere.py.

Member Function Documentation

gradSdf()

sdf_sphere.Sphere.gradSdf	(		xc,
			yc,
			zc,
			x,
			y,
			z
	)

Definition at line 29 of file sdf_sphere.py.

    def gradSdf(xc, yc, zc, x, y, z):
        a = (x-xc)**2 + (y-yc)**2 + (z-zc)**2
        c_val = np.sqrt(a)
        c_val_A = 1./c_val
        c_val_dx = c_val_A * (x-xc)
        c_val_dy = c_val_A * (y-yc)
        c_val_dz = c_val_A * (z-zc)
 
        c_val_dx = c_val_dx.reshape((-1,1))
        c_val_dy = c_val_dy.reshape((-1,1))
        c_val_dz = c_val_dz.reshape((-1,1))
        grad_array = np.hstack([c_val_dx,c_val_dy,c_val_dz])
        return grad_array
 

hessSdf()

sdf_sphere.Sphere.hessSdf	(		xc,
			yc,
			zc,
			x,
			y,
			z
	)

Definition at line 43 of file sdf_sphere.py.

    def hessSdf(xc, yc, zc, x, y, z):
        x = x-xc
        y = y-yc
        z = z-zc
        Hxx = -x**2/(x**2 + y**2 + z**2)**(3/2) + 1/np.sqrt(x**2 + y**2 + z**2)
        Hzx = -x*z/(x**2 + y**2 + z**2)**(3/2)
        Hxy = -x*y/(x**2 + y**2 + z**2)**(3/2)
        Hyy = -y**2/(x**2 + y**2 + z**2)**(3/2) + 1/np.sqrt(x**2 + y**2 + z**2)
        Hzy = -y*z/(x**2 + y**2 + z**2)**(3/2)
        Hzz = -z**2/(x**2 + y**2 + z**2)**(3/2) + 1/np.sqrt(x**2 + y**2 + z**2)
        # xx, yx, zx, yy, zy, zz
        Hxx = Hxx.reshape((-1,1))
        Hzx = Hzx.reshape((-1,1))
        Hxy = Hxy.reshape((-1,1))
        Hyy = Hyy.reshape((-1,1))
        Hzy = Hzy.reshape((-1,1))
        Hzz = Hzz.reshape((-1,1))
        hess_array = np.hstack([Hxx, Hxy, Hzx, Hyy, Hzy, Hzz])
 
        return hess_array
 
 
 
 
 

sDF()

sdf_sphere.Sphere.sDF	(		r,
			xc,
			yc,
			zc,
			x,
			y,
			z
	)

Definition at line 21 of file sdf_sphere.py.

    def sDF(r, xc, yc, zc, x, y, z):
        dx = np.subtract(x, xc)
        dy = np.subtract(y, yc)
        dz = np.subtract(z, zc)
        a = dx**2 + dy**2 + dz**2
        gap = np.sqrt(a) - r
        return gap
 

---

The documentation for this class was generated from the following file:

• users_modules/eshelbian_plasticity/examples/sdf_sphere.py