""" 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)