Ejemplo n.º 1
0
            
        """ 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)
Ejemplo n.º 2
0
        """ 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())


Ejemplo n.º 3
0
                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))
Ejemplo n.º 4
0
            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)