cvit.tst 3.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128
  1. % Tests of Cvitanovic Package.
  2. % COPYRIGHT (C) 1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE UNIV.
  3. % CVITBUBLE TEST OF CVITANOVIC PACKAGE
  4. % AUTHOR A. KRYUKOV, ARODIONOV, A.TARANOV
  5. % VERSION 1.1
  6. % RELEASE 18-SEP-90
  7. index j1,j2,j3,j4,j5,j6,j7,j8,j9,j0;
  8. vecdim n$
  9. % Tests of the weels with buble
  10. % (Use notation from SIGSAM Bull, 1989, v.23, no.4, pp.15-24)
  11. g(l,j1,j2,j2,j1);
  12. g(l,j1,j2)*g(l1,j3,j1,j2,j3);
  13. g(l,j1,j2)*g(l1,j3,j1,j3,j2);
  14. g(l,j1,j2)*g(l1,j3,j3,j2,j1);
  15. g(l,j1,j2,j3,j4)*g(l1,j1,j2,j3,j4);
  16. g(l,j1,j2)*g(l1,j3,j4,j1,j2,j4,j3);
  17. g(l,j1,j2,j3,j4)*g(l1,j1,j4,j2,j3);
  18. g(l,j1,j2)*g(l1,j3,j4,j1,j4,j3,j2);
  19. g(l,j1,j2)*g(l1,j3,j4,j5,j1,j2,j3,j4,j5);
  20. g(l,j1,j2,j3,j4)*g(l1,j5,j1,j2,j3,j5,j4);
  21. g(l,j1,j2,j3,j4,j5,j1)*g(l1,j2,j5,j3,j4);
  22. g(l,j1,j2,j3,j4,j5,j1,j2,j5)*g(l1,j4,j3);
  23. g(l,j1,j2)*g(l1,j3,j4,j5,j6,j1,j2,j3,j4,j5,j6);
  24. g(l,j1,j2,j3,j4)*g(l1,j5,j6,j1,j2,j3,j4,j6,j5);
  25. g(l,j1,j2,j3,j4,j5,j6)*g(l1,j1,j2,j4,j3,j6,j5);
  26. g(l,j1,j2,j3,j4,j5,j6,j1,j2)*g(l1,j6,j3,j4,j5);
  27. g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5)*g(l1,j6,j7);
  28. g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3)*g(l1,j4,j5,j7,j6);
  29. g(l,j1,j2,j3,j4,j5,j6,j7,j2)*g(l1,j1,j3,j4,j5,j6,j7);
  30. % COPYRIGHT (C) 1988,1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE U.
  31. % CVITEST Test of CVITANOVIC PACKAGE
  32. % AUTHOR A. KRYUKOV, A.RODIONOV, A.TARANOV
  33. % VERSION 1.2
  34. % RELEASE 11-MAR-90
  35. %
  36. % Test for trace of Dirac matrices.
  37. %
  38. % All tests are the lattices with difference lines
  39. % (Use notation from SIGSAM Bull, 1989, v.4,no.23, pp.15-24)
  40. index m1,m2,m3,m4,m5,m6,m7,m8,m9,m0;
  41. index n1,n2,n3,n4,n5,n6,n7,n8,n9,n0;
  42. vecdim n;
  43. g(l,n1,n1);
  44. g(l,n1,m1,n1,m1);
  45. g(l,n1,n2,n2,n1);
  46. g(l,n1,n2,m1,n2,n1,m1);
  47. g(l,n1,n2,m1,m2,n2,n1,m2,m1);
  48. g(l,n1,n2,n3,n3,n2,n1);
  49. g(l,n1,n2,n3,m1,n3,n2,n1,m1);
  50. g(l,n1,n2,n3,m1,m2,n3,n2,n1,m2,m1);
  51. g(l,n1,n2,n3,m1,m2,m3,n3,n2,n1,m3,m2,m1);
  52. g(l,n1,n2,n3,m1,n3,n1,n2,m1);
  53. g(l,n1,n2,n3,m1,m2,n3,n1,n2,m1,m2);
  54. g(l,n1,n2,n3,m1,m2,m3,n2,n3,n1,m3,m1,m2);
  55. % COPYRIGHT (C) 1988,1990, INSTITUTE OF NUCLEAR PHYSICS, MOSCOW STATE U.
  56. % CVITWEEL TEST OF CVITANOVIC PACKAGE
  57. % AUTHOR A. KRYUKOV, ARODIONOV, A.TARANOV
  58. % VERSION 1.2
  59. % RELEASE 11-MAR-90
  60. index j1,j2,j3,j4,j5,j6,j7,j8,j9,j0;
  61. vecdim n$
  62. % Test of CVITANOVIC PACKAGE
  63. %
  64. % All tests are the weels with defferent spoke
  65. % (Use notation from SIGSAM Bull, 1989, v.23, no.4, pp.15-24)
  66. g(l,j1,j2,j2,j1);
  67. g(l,j1,j2,j3,j1,j2,j3);
  68. g(l,j1,j2,j3,j1,j3,j2);
  69. g(l,j1,j2,j3,j3,j2,j1);
  70. g(l,j1,j2,j3,j4,j1,j2,j3,j4);
  71. g(l,j1,j2,j3,j4,j1,j2,j4,j3);
  72. g(l,j1,j2,j3,j4,j1,j4,j2,j3);
  73. g(l,j1,j2,j3,j4,j1,j4,j3,j2);
  74. g(l,j1,j2,j3,j4,j5,j1,j2,j3,j4,j5);
  75. g(l,j1,j2,j3,j4,j5,j1,j2,j3,j5,j4);
  76. g(l,j1,j2,j3,j4,j5,j1,j2,j5,j3,j4);
  77. g(l,j1,j2,j3,j4,j5,j1,j2,j5,j4,j3);
  78. g(l,j1,j2,j3,j4,j5,j6,j1,j2,j3,j4,j5,j6);
  79. g(l,j1,j2,j3,j4,j5,j6,j1,j2,j3,j4,j6,j5);
  80. g(l,j1,j2,j3,j4,j5,j6,j1,j2,j4,j3,j6,j5);
  81. g(l,j1,j2,j3,j4,j5,j6,j1,j2,j6,j3,j4,j5);
  82. g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5,j6,j7);
  83. g(l,j1,j2,j3,j4,j5,j6,j7,j1,j2,j3,j4,j5,j7,j6);
  84. g(l,j1,j2,j3,j4,j5,j6,j7,j2,j1,j3,j4,j5,j6,j7);
  85. % Test of example that calculated incorrectly in earlier package.
  86. index ix,iy,iz;
  87. mass p1=mm, p2=mm, p3=mm, p4=mm, k1=0;
  88. mshell p1,p2,p3,p4,k1;
  89. vector q1,q2;
  90. operator ga,gb;
  91. for all p let ga(p)=g(la,p) + mm,
  92. gb(p)=g(lb,p) + mm;
  93. xx := g(la,ix)*g(la,iy)*(g(lb,ix)*gb(p1)*g(lb,iy)*gb(q2) +
  94. gb(p3)*g(lb,ix)*g(lb,iy));
  95. let q1=p1-k1, q2=p3+k1;
  96. xx;
  97. end;