oldtests/gauntgen.m

   1 gaunt[m_, n_, mu_, nu_,
   2    p_] := (-1)^(m + mu) (2 p + 1) Sqrt[
   3     Factorial[n + m] Factorial[
   4       nu + mu] Factorial[p - m - mu]/Factorial[n - m]/
   5       Factorial[nu - mu] / Factorial[p + m + mu]] ThreeJSymbol[{n,
   6      0}, {nu, 0}, {p, 0}] ThreeJSymbol[{n, m}, {nu,
   7      mu}, {p, -m - mu}]
   8
   9 lMax := 30
  10 For[n = 0, n <= lMax, n++,
  11  For[nu = 0, nu <= lMax, nu++,
  12   For[m = -n, m <= n, m++,
  13    For[mu = -nu, mu <= nu, mu++,
  14     For[q = 0, q <= Min[n, nu, (n + nu - Abs[m + mu])/2],q++,
  15      Print[StringForm["{`1`, `2`, `3`, `4`, `5`, `6`},",m,n,mu,nu,q,CForm[N[gaunt[m, n, mu, nu, n + nu - 2 q],32]]]]
  16      ]
  17     ]
  18    ]
  19   ]
  20  ]
  21