📑 Circus of Programming Tricks 01. Animated Graphs.

tricks_01.html

xxxxxxxxxx

 

1

from sage.plot.animate import Animation

2

class MyAnimation(Animation):

3

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

4

        G=frame.plot()

5

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

6

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

7

        G.save_image(filename,**kwds)

8

def rotate_line(a,l): 

9

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

10

            for [x,y] in l]

11

def linegraph(g):

12

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

13

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

14

            for i in range(ne)]  

15

def col(): return (.3+random()/3,0,.3+random()/3)

16

def rotate_bezier_path_cubegraph(m,k,color='default',th=.5):

17

    g=graphs.CubeGraph(m) 

18

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

19

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

20

          for i in range(len(lg))] for j in range(2*k)] 

21

    for i in range(2*k): 

22

        if color=='default': color=col()

23

        LG+=bezier_path(rlg[i],color=color,thickness=th)

24

    return LG

25

@interact

26

def _(m1=[3,4,5,6],m2=[3,4,5,6],

27

      th1=[.5,.7,.9,1.1],th2=[.3,.7,1.1,1.5]):

28

    M=MyAnimation([sum(

29

        [rotate_bezier_path_cubegraph(m1,k,th=th1),

30

         rotate_bezier_path_cubegraph(m2,k,color='white',th=th2)]) 

31

                  for k in range(3,9)],figsize=(5,5),axes=False)

32

    M.show()

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

---

Messages