doc/example1.py

   1 #!/home/mjo/bin/sage
   2
   3 # example1.py
   4 #
   5 # In the paper, we display a lot of sage code to be entered at the
   6 # prompt. It's not much fun re-entering all of these formulas every
   7 # time you'd like to work on the paper, so a quick load('example1.py')
   8 # from the sage prompt will pull all of this stuff in for you.
   9 #
  10
  11 x,y,z = var('x,y,z')
  12
  13 f = 1 + x + x*y + x*y*z
  14
  15 g = [[[ f(x=a, y=b, z=c) for a in range(0,3) ] for b in range(0,3) ] for c in range(0,3) ]
  16
  17 def grid(x,y,z):
  18     return g[z][y][x]
  19
  20 def I(x,y,z):
  21     return grid(x,y,z)
  22
  23 def F(x,y,z):
  24     return grid(x-1, y, z)
  25
  26 def B(x,y,z):
  27     return grid(x+1, y, z)
  28
  29 def L(x,y,z):
  30     return grid(x, y-1, z)
  31
  32 def R(x,y,z):
  33     return grid(x, y+1, z)
  34
  35 def T(x,y,z):
  36     return grid(x, y, z+1)
  37
  38 def D(x,y,z):
  39     return grid(x, y, z-1)
  40
  41 def FL(x,y,z):
  42     return grid(x-1, y-1, z)
  43
  44 def FR(x,y,z):
  45     return grid(x-1, y+1, z)
  46
  47 def FD(x,y,z):
  48     return grid(x-1, y, z-1)
  49
  50 def FT(x,y,z):
  51     return grid(x-1, y, z+1)
  52
  53 def BL(x,y,z):
  54     return grid(x+1, y-1, z)
  55
  56 def BR(x,y,z):
  57     return grid(x+1, y+1, z)
  58
  59 def BD(x,y,z):
  60     return grid(x+1, y, z-1)
  61
  62 def BT(x,y,z):
  63     return grid(x+1, y, z+1)
  64
  65 def LD(x,y,z):
  66     return grid(x, y-1, z-1)
  67
  68 def LT(x,y,z):
  69     return grid(x, y-1, z+1)
  70
  71 def RD(x,y,z):
  72     return grid(x, y+1, z-1)
  73
  74 def RT(x,y,z):
  75     return grid(x, y+1, z+1)
  76
  77 def FLD(x,y,z):
  78     return grid(x-1, y-1, z-1)
  79
  80 def FLT(x,y,z):
  81     return grid(x-1, y-1, z+1)
  82
  83 def FRD(x,y,z):
  84     return grid(x-1, y+1, z-1)
  85
  86 def FRT(x,y,z):
  87     return grid(x-1, y+1, z+1)
  88
  89 def BLD(x,y,z):
  90     return grid(x+1, y-1, z-1)
  91
  92 def BLT(x,y,z):
  93     return grid(x+1, y-1, z+1)
  94
  95 def BRD(x,y,z):
  96     return grid(x+1, y+1, z-1)
  97
  98 def BRT(x,y,z):
  99     return grid(x+1, y+1, z+1)
 100