tests/ehrhart/e14.in

   1 # Philippe's simple example
   2 #--------------------------
   3 # 1 <= i <= N
   4 # 1 <= j <= N
   5 # 1 <= k <= N
   6 # i+j+k-2 = M
   7 #--------------------------
   8 7 7
   9 #  i  j  k  N  M  1
  10 1  1  0  0  0  0 -1
  11 1 -1  0  0  1  0  0
  12 1  0  1  0  0  0 -1
  13 1  0 -1  0  1  0  0
  14 1  0  0  1  0  0 -1
  15 1  0  0 -1  1  0  0
  16 0 -1 -1 -1  0  1  2
  17 #--------------------------
  18 0 4
  19
  20 N M
  21 Philippe's simple example
  22
  23 #--------------------------
  24 # I have done the work by hand, and here is the right
  25 # answer :
  26 #
  27 # If 1 <= M <= N then EP(N,M) = 1/2 M^2 + 1/2 M
  28 # If N <= M <= 2N-1 then EP(N,M) = -3/2 N^2 - M^2 + 3NM - M
  29 # If 2N-1 <= M <= 3N-2 then EP(N,M) = 9/2 N^2 + 1/2 M^2 - 3NM - 3/2 N + 1/2 M