PointGroup.h

#ifndef UTIL_POINT_GROUP_H
#define UTIL_POINT_GROUP_H

#include <util/crystal/PointSymmetry.h>
#include <util/crystal/SymmetryGroup.h>
#include <util/containers/FSArray.h>
#include <util/space/IntVector.h>

#include <iostream>

namespace Util {

   class PointGroup : public SymmetryGroup<PointSymmetry>
   {

   public:

      template <int N>
      void makeStar(const IntVector& root, FSArray<IntVector, N>& star);

   // friends:

      friend 
      std::ostream& operator << (std::ostream& out, const PointGroup& A);

   };

   std::ostream& operator << (std::ostream& out, const PointGroup& A);

   // Template function definition

   template <int N>
   void PointGroup::makeStar(const IntVector& root, FSArray<IntVector, N>& star)
   {
      // Precondition
      if (star.capacity() < size()) {
         UTIL_THROW("Star array capacity is less than order of group");
      }

      IntVector vector;
      int       i, j;
      bool      found;
      star.clear();
      for (i = 0; i < size(); ++i) {
         vector = root*(*this)[i];
         found = false;
         for (j = 0; j < star.size(); ++j) {
            if (vector == star[j]) {
               found = true;
               break;
            }
         }
         if (!found) {
            star.append(vector);
         }
      }
   
   }

}
#endif