Animated Graphs

xxxxxxxxxx

 

1

import warnings; warnings.simplefilter('ignore')

2

from sage.plot.animate import Animation

3

class MyAnimation(Animation):

4

    def make_image(self,frame,filename,**kwds):

5

        G=frame.plot()

6

        G.set_axes_range(floor(G.xmin()),ceil(G.xmax()),

7

                         floor(G.ymin()),ceil(G.ymax()))

8

        G.save_image(filename,**kwds)

9

def rotate_line(a,l): 

10

    return [[(cos(a)*x-sin(a)*y).n(),

11

             (sin(a)*x+cos(a)*y).n()] 

12

            for [x,y] in l]

13

def linegraph(g):

14

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

15

    return [[list(vg[eg[i][0]]),

16

             list(vg[eg[i][1]])] 

17

            for i in range(ne)]    

18

def rotate_bezier_path_cubegraph(m,k):

19

    g=graphs.CubeGraph(m) 

20

    lg=linegraph(g); LG=Graphics()

21

    rlg=[[rotate_line(j*pi/k,lg[i]) 

22

          for i in range(len(lg))] 

23

         for j in range(2*k)] 

24

    for i in range(2*k): 

25

        col=(.3+random()/3,0,.3+random()/3)

26

        LG+=bezier_path(rlg[i],color=col,

27

                        thickness=.3)

28

    return LG

29

MyAnimation([rotate_bezier_path_cubegraph(3,k) 

30

             for k in range(2,6)],

31

            figsize=(5,5),axes=False)

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

---

Messages