#include <iostream>
#include <iomanip>
#include <cmath>
#include "ChemPlugin.h"


int main(int argc, char** argv) {
   int nlevel = 4;
   std::cout << "Link ChemPlugin instances into a " << nlevel
             << " bifurcating tree" << std::endl << std::endl;

   // Create the ChemPlugin instances.
   int ninst = pow(2, nlevel) - 1;
   ChemPlugin *cp = new ChemPlugin[ninst];

   // Link the instances into a tree.
   int inst = 0, linked_inst = 1;
   for (int level=0; level<nlevel-1; level++) {
      for (int i=0; i<pow(2, level); i++) {
         cp[linked_inst++].Link(cp[inst]);
         cp[linked_inst++].Link(cp[inst]);
         inst++;
      }
   }

   // Report the number of links to each instance.
   for (int i=0; i<ninst; i++)
      std::cout << "Instance " << i << " is linked to " << cp[i].nLinks()
                << " instance(s)" << std::endl;
      
   // Any keystroke closes the console.
   std::cin.get();
   delete[] cp;
   return 0;
}