RingTetraRebridgeMove.h

#ifndef MCMD_RING_TETRA_REBRIDGE_MOVE_H
#define MCMD_RING_TETRA_REBRIDGE_MOVE_H

/*
* Simpatico - Simulation Package for Polymeric and Molecular Liquids
*
* Copyright 2010 - 2017, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/

#include <mcMd/mcMoves/base/GroupRebridgeBase.h>  // base class
#include <util/space/Vector.h>                     // template parameter
#include <util/containers/DArray.h>               // member template

namespace McMd
{

   using namespace Util;

   class McSystem;
   class Molecule;

   class RingTetraRebridgeMove : public GroupRebridgeBase
   {
   
   public:
   
      RingTetraRebridgeMove(McSystem& system);
   
      virtual void readParameters(std::istream& in);

      virtual bool move();
   
      virtual void loadParameters(Serializable::IArchive& ar);

      virtual void save(Serializable::OArchive& ar);
  
   protected:
   
      int speciesId_;

      int nAtom_;

      double upperBridge_;

      double lowerBridge_;

      DArray<Vector> R_;

      bool scanBridge(Molecule* molPtr, int aId, int &cId);

      int modId(int id, int n);

   };

   // Private inline methods.

   inline int RingTetraRebridgeMove::modId(int id, int n)
   {
      if (id >= n) id = id%n;
      if (id < 0) id = n - (-id)%n;
      return id;
   }

}      
#endif