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


int main(int argc, char** argv) {
   int nx = 5, ny = 3;
   std::cout << "Link ChemPlugin instances into a " << nx << " by " << ny
             << " grid" << std::endl << std::endl;

   // Create the ChemPlugin instances.
   ChemPlugin cp_left, cp_right, cp_bottom, cp_top;
   ChemPlugin *cp = new ChemPlugin[nx*ny];

   // Link the instances into a tree.
   for (int j=0; j<ny; j++) {
      for (int i=0; i<nx; i++) {
         int ij = i + j*nx;
         cp[ij].Link(i == 0? cp_left : cp[ij-1]);
         cp[ij].Link(j == 0? cp_bottom : cp[ij-nx]);
         if (j == ny-1) cp_top.Link(cp[ij]);
      }
      cp_right.Link(cp[(j+1)*nx-1]);
   }

   // Report the number of links to each instance.
   std::cout << "Left boundary is linked to " << cp_left.nLinks()
             << " instances" << std::endl;
   std::cout << "Bottom boundary is linked to " << cp_bottom.nLinks()
             << " instances" << std::endl;
   for (int i=0; i<nx*ny; i++)
      std::cout << "Instance " << i << " is linked to " << cp[i].nLinks()
                << " instances" << std::endl;
   std::cout << "Top boundary is linked to " << cp_top.nLinks()
             << " instances" << std::endl;
   std::cout << "Right boundary is linked to " << cp_right.nLinks()
             << " instances" << std::endl;
      
   // Any keystroke closes the console.
   std::cin.get();
   delete[] cp;
   return 0;
}