def grow(self,plans,length = None,side = None,force = False):
side = self.grow_side(plans,side)
if side is None:return False
gleng = self.grow_length(plans,length,side)
if gleng is None:return False
self.face_away(side)
v1,v2 = side[0]
c1 = v2.copy()
c2 = v1.copy()
c3,c4 = dpr.extrude_edge(c1,c2,gleng,side[1]['normal'])
newcorners = [c1,c2,c3,c4]
x,y,z = dpv.center_of_mass(newcorners)
xpj = dpv.project_coords(newcorners,dpv.x())
ypj = dpv.project_coords(newcorners,dpv.y())
l,w = xpj.y-xpj.x,ypj.y-ypj.x
margs = ((),{'x':x,'y':y,'l':l,'w':w,'shafted':False})
self.create_bbox(margs)
if not self.verify_growth(plans,margs):return False
splan = self.should_shaft(plans,margs)
cpairs = [(c2,c3),(c3,c4),(c4,c1)]
iwargs = []
ewargs = [[cp,
{'h':4.0,'w':0.5,'walltype':'exterior','room':margs}]
for cp in cpairs]
self.switch_walltype(plans,side)
if self.resolve_walls(plans,ewargs,iwargs,margs):
rps,eps,ips,sps = plans
rps.append(margs)
eps.extend(ewargs)
ips.extend(iwargs)
if splan:sps.append(splan)
return True
else:return False
def create_bbox(self, roomplan):
x, y = roomplan[1]['x'], roomplan[1]['y']
l, w = roomplan[1]['l'] - 0.01, roomplan[1]['w'] - 0.01
cns = dpr.square(l, w, dpv.vector(x, y, 0))
xpj = dpv.project_coords(cns, dpv.x())
ypj = dpv.project_coords(cns, dpv.y())
zpj = dpv.project_coords(cns, dpv.z())
bb = dbb.bbox(xpj, ypj, zpj)
roomplan[1]['bbox'] = bb
def create_bbox(self,roomplan):
x,y = roomplan[1]['x'],roomplan[1]['y']
l,w = roomplan[1]['l']-0.01,roomplan[1]['w']-0.01
cns = dpr.square(l,w,dpv.vector(x,y,0))
xpj = dpv.project_coords(cns,dpv.x())
ypj = dpv.project_coords(cns,dpv.y())
zpj = dpv.project_coords(cns,dpv.z())
bb = dbb.bbox(xpj,ypj,zpj)
roomplan[1]['bbox'] = bb
def _project_uv_flat(self,rng = None):
if rng is None:rng = range(len(self.faces))
for nf in rng:
face = self.faces[nf]
for fdx in face:
p = self.pcoords[fdx]
n = self.ncoords[fdx]
if dpv.near(n,dpv.nx()) or dpv.near(n,dpv.x()):
nu = p.copy().yz2d()
elif dpv.near(n,dpv.ny()) or dpv.near(n,dpv.y()):
nu = p.copy().xz2d()
elif dpv.near(n,dpv.nz()) or dpv.near(n,dpv.z()):
nu = p.copy().xy2d()
else:continue
self.ucoords[fdx] = nu
def grow(self, plans, length=None, side=None, force=False):
side = self.grow_side(plans, side)
if side is None: return False
gleng = self.grow_length(plans, length, side)
if gleng is None: return False
self.face_away(side)
v1, v2 = side[0]
c1 = v2.copy()
c2 = v1.copy()
c3, c4 = dpr.extrude_edge(c1, c2, gleng, side[1]['normal'])
newcorners = [c1, c2, c3, c4]
x, y, z = dpv.center_of_mass(newcorners)
xpj = dpv.project_coords(newcorners, dpv.x())
ypj = dpv.project_coords(newcorners, dpv.y())
l, w = xpj.y - xpj.x, ypj.y - ypj.x
margs = ((), {'x': x, 'y': y, 'l': l, 'w': w, 'shafted': False})
self.create_bbox(margs)
if not self.verify_growth(plans, margs): return False
splan = self.should_shaft(plans, margs)
cpairs = [(c2, c3), (c3, c4), (c4, c1)]
iwargs = []
ewargs = [[
cp, {
'h': 4.0,
'w': 0.5,
'walltype': 'exterior',
'room': margs
}
] for cp in cpairs]
self.switch_walltype(plans, side)
if self.resolve_walls(plans, ewargs, iwargs, margs):
rps, eps, ips, sps = plans
rps.append(margs)
eps.extend(ewargs)
ips.extend(iwargs)
if splan: sps.append(splan)
return True
else:
return False
def inpolyhedron(point,triangles):
pray = dry.ray(point,dpv.x())
isects = []
isects = dry.intersect_hits(pray,triangles)
#print('inininin',len(triangles),isects)
'''#
for x in range(len(triangles)):
t1,t2,t3 = triangles[x]
isct = pray.intersect_tri(t1,t2,t3)
if isct == 1:
t = pray.cast.x
print('t',t)
isects.append(t)
else:print('noinsct',x)
'''#
ins = len(isects) % 2 > 0
pray = dry.ray(point,dpv.y())
isects = []
isects = dry.intersect_hits(pray,triangles)
ins = ins or len(isects) % 2 > 0
#print('insideeee',ins,len(isects))
'''#
ax = dtl.plot_axes()
for x in range(len(triangles)):
ax = dtl.plot_polygon(list(triangles[x]),ax)
ax = dtl.plot_edges([
pray.origin,pray.origin.copy().translate(
pray.direction.copy().scale_u(100))],ax)
plt.show()
'''#
return ins
def roll_up(self,p,d):d.rotate(dpq.q_from_av(self.angle,dpv.y())) def roll_down(self,p,d):d.rotate(dpq.q_from_av(self.angle,dpv.ny()))
def polar_down(self,p,d):
if d.near(dpv.z()) or d.near(dpv.nz()):qv = dpv.y()
else:qv = d.cross(dpv.z())
d.rotate(dpq.q_from_av(-self.polar,qv))