3D Layer Graphs (Rotating) 2

xxxxxxxxxx

 

1

import warnings

2

warnings.filterwarnings("ignore",category=DeprecationWarning)

3

def col(): 

4

    c1,c2,c3=random(),random(),.8

5

    return sage.plot.colors.float_to_html(c1,c2,c3)

6

def graph_pe(g,n):

7

    PG=Graphics(); k=randint(70,120)/100

8

    vg=g.get_pos(); vgp=list(vg.values())

9

    eg=g.edges(); nv=len(vg); ne=len(eg)

10

    x0=mean([list(vg.values())[i][0] for i in range(nv)])

11

    y0=mean([list(vg.values())[i][1] for i in range(nv)])

12

    cp=[[vgp[i][0]-x0,vgp[i][1]-y0] for i in range(nv)]

13

    l=[[list(vg[eg[i][0]]),list(vg[eg[i][1]])] 

14

       for i in range(ne)]

15

    cl=[[[x-x0,y-y0] for [x,y] in l[i]] 

16

        for i in range(ne)]

17

    cp3d=[[[cos(j*pi/k)*x-sin(j*pi/k)*y,

18

            sin(j*pi/k)*x+cos(j*pi/k)*y,(j+1)/k] 

19

           for j in range(n)] for [x,y] in cp]

20

    cl3d=[[[cos(j*pi/k)*x-sin(j*pi/k)*y,

21

            sin(j*pi/k)*x+cos(j*pi/k)*y,(j+1)/k] 

22

           for [x,y] in cl[i]] 

23

          for j in range(n) for i in range(ne)]

24

    for i in range(nv): 

25

        PG+=line3d(cp3d[i],color=col())

26

    for i in range(n*ne): 

27

        PG+=line3d(cl3d[i],color=col())

28

    PG.show(frame=False)

29

g=graphs.CubeGraph(randint(3,7))

30

graph_pe(g,randint(3,6))

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages