oldtests/gauntmaxcumsum.c

   1 #include "../qpms/gaunt.h"
   2 #include <stdio.h>
   3 #include <stdlib.h>
   4
   5 //const int lMax = 30;
   6
   7 int main(int argc, char **argv)
   8 {
   9         int lMax;
  10         if (argc < 2)
  11                 lMax = 30;
  12         else lMax = atoi(argv[1]);
  13
  14         printf("// assuming lMax == %d:\n", lMax);
  15         size_t qmaxcumsum = 0;
  16         printf("0x%zx,\n", qmaxcumsum);
  17         for (int n = 0; n <= lMax; n++)
  18                 for (int m = -n; m <= n; m++)
  19                         for (int nu = 0; nu <= lMax; nu++)
  20                                 for (int mu = -nu; mu <= nu; mu++) {
  21                                         qmaxcumsum += gaunt_q_max(m, n, mu, nu) + 1;
  22                                         printf("0x%zx,\n", qmaxcumsum);
  23                                 }
  24         return 0;
  25 }