initialize1()

void WeylAlgebra::initialize1 ( )

private

Definition at line 175 of file weylalg.cpp.

{
  if (binomtable == nullptr)
    {
      int i, j;
 
      binomtable = newarray(int *, binomtop + 1);
      for (i = 0; i <= binomtop; i++)
        binomtable[i] = newarray_atomic(int, i + 1);
      binomtable[0][0] = 1;
      binomtable[1][0] = 1;
      binomtable[1][1] = 1;
      for (i = 2; i <= binomtop; i++)
        {
          binomtable[i][0] = 1;
          binomtable[i][i] = 1;
          for (j = 1; j < i; j++)
            binomtable[i][j] = binomtable[i - 1][j - 1] + binomtable[i - 1][j];
        }
 
      diffcoeffstable = newarray(int *, diffcoeffstop + 1);
      for (i = 0; i <= diffcoeffstop; i++)
        diffcoeffstable[i] = newarray_atomic(int, i + 1);
      diffcoeffstable[0][0] = 1;
      diffcoeffstable[1][0] = 1;
      diffcoeffstable[1][1] = 1;
      for (i = 2; i <= diffcoeffstop; i++)
        {
          diffcoeffstable[i][0] = 1;
          for (j = 1; j <= i; j++)
            diffcoeffstable[i][j] = i * diffcoeffstable[i - 1][j - 1];
        }
#if 0
//       // Display the binomial tables:
//       cout << "---binom table---" << endl;
//       for (i=0; i<=binomtop; i++)
//      {
//        for (j=0; j<=i; j++)
//          cout << "  " << binomtable[i][j];
//        cout << endl;
//      }
//       cout << "---diff table---" << endl;
//       for (i=0; i<=diffcoeffstop; i++)
//      {
//        for (j=0; j<=i; j++)
//          cout << "  " << diffcoeffstable[i][j];
//        cout << endl;
//      }
#endif
    }
}

References binomtable, binomtop, diffcoeffstable, diffcoeffstop, newarray, and newarray_atomic.

Referenced by initialize_weyl().