""" build solid basket """
newSolid = g.addNode(3)
faces = list(set(faces + [cell]))
for face in faces: g.addArch(face,newSolid)
nodes = DOWNCELLS(g)(newSolid)
for i in range(2):
nodes = list(set(CAT(AA(DOWNCELLS(g))(nodes))))
for node in nodes:
g.addArch(newSolid,node)
'myprint("CELLSPERLEVEL(g)(0)",CELLSPERLEVEL(g)(0))'
'myprint("CELLSPERLEVEL(g)(1)",CELLSPERLEVEL(g)(1))'
'myprint("CELLSPERLEVEL(g)(2)",CELLSPERLEVEL(g)(2))'
'myprint("CELLSPERLEVEL(g)(3)",CELLSPERLEVEL(g)(3))'
'myprint("mapping",mapping)'
DRAW(g,[2.5,2.5,2.5])(CELLSPERLEVEL(g)(3))
return g
h = ngon(6)
h = pcsphere(tetra())
h = pcsphere(h,3)
g = hccmesh(h)
DRAW(g,[2.5,2.5,2.5])(CELLSPERLEVEL(g)(3))
h = pcsphere(h,4)
""" add 1-edges from tips to mapped 0-nodes of bottom faces"""
for node in list(set(CAT(AA(DOWNCELLS(g))(DOWNCELLS(g)(cell))))):
'myprint("node",(node))'
newArc = g.addNode(1)
g.addArch(mapping[node],newArc)
g.addArch(newCell,newArc)
myprint("CELLSPERLEVEL(g)(0)",CELLSPERLEVEL(g)(0))
myprint("CELLSPERLEVEL(g)(1)",CELLSPERLEVEL(g)(1))
myprint("CELLSPERLEVEL(g)(2)",CELLSPERLEVEL(g)(2))
myprint("mapping",mapping)
DRAW(g)()
# STEP 5: (d-1)-chain pivoting \\\\\\\\\\\\\\\\\\\\\\\\\\
# STEP 6: (d-1)-chain roofing \\\\\\\\\\\\\\\\\\\\\\\\\\\
# STEP 7: (d)-chain construction \\\\\\\\\\\\\\\\\\\\\\\\
return g
h = ngon(6)
#h = pcsphere(h)
h = pcsphere(tetra())
CAT(AA(UPCELLS(g))(hedges))))
wedgeFaces += sideFaces
for face in wedgeFaces: g.addArch(face,wedgeCell)
'myprint("CELLSPERLEVEL(g)(0)",CELLSPERLEVEL(g)(0))'
'myprint("CELLSPERLEVEL(g)(1)",CELLSPERLEVEL(g)(1))'
'myprint("CELLSPERLEVEL(g)(2)",CELLSPERLEVEL(g)(2))'
'myprint("CELLSPERLEVEL(g)(3)",CELLSPERLEVEL(g)(3))'
'myprint("mapping",mapping)'
#DRAW(g,[2.5,2.5,2.5])(CELLSPERLEVEL(g)(3))
return g
h = pcsphere(tetra(),trihedrons=True)
DRAW(h)()
##g = hccmesh(h)
##DRAW(g)()
##h = pcsphere(h,2,True)
##DRAW(h)()
##g = hccmesh(h)
##DRAW(g)()
##h = pcsphere(h,3,True)
##DRAW(h)()
##g = hccmesh(h)
##DRAW(g)()
##h = pyramid()
##DRAW(h)(CELLSPERLEVEL(h)(0)+CELLSPERLEVEL(h)(1))
wedgeFaces += sideFaces
for face in wedgeFaces:
g.addArch(face, wedgeCell)
'myprint("CELLSPERLEVEL(g)(0)",CELLSPERLEVEL(g)(0))'
'myprint("CELLSPERLEVEL(g)(1)",CELLSPERLEVEL(g)(1))'
'myprint("CELLSPERLEVEL(g)(2)",CELLSPERLEVEL(g)(2))'
'myprint("CELLSPERLEVEL(g)(3)",CELLSPERLEVEL(g)(3))'
'myprint("mapping",mapping)'
#DRAW(g,[2.5,2.5,2.5])(CELLSPERLEVEL(g)(3))
return g
h = pcsphere(tetra(), trihedrons=True)
DRAW(h)()
##g = hccmesh(h)
##DRAW(g)()
##h = pcsphere(h,2,True)
##DRAW(h)()
##g = hccmesh(h)
##DRAW(g)()
##h = pcsphere(h,3,True)
##DRAW(h)()
##g = hccmesh(h)
##DRAW(g)()
##h = pyramid()
##DRAW(h)(CELLSPERLEVEL(h)(0)+CELLSPERLEVEL(h)(1))
##h = hccmesh(h)
