CfbRingRebridgeMove.h

#ifndef MCMD_CFB_RING_REBRIDGE_MOVE_H
#define MCMD_CFB_RING_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/CfbRebridgeBase.h>  // base class
#include <util/containers/DArray.h>             // member template
#include <util/space/Vector.h>                   // member template parameter

namespace McMd
{

   using namespace Util;

   class McSystem;
   class Molecule;

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

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

      virtual void save(Serializable::OArchive& ar);
  
      virtual bool move();
   
   protected:
   
      DArray<Vector> oldPos_;
    
      int speciesId_;
   
      int nRegrow_;

   private:

      void deleteSequence(int sign, Molecule *molPtr, int endId, int *bonds, 
                         double &rosenbluth, double &energy);

      void addSequence(int sign, Molecule *molPtr, int beginId, int *bonds, 
                         double &rosenbluth, double &energy);

      int modId(int id, int n);

   };

   // Private inline methods.

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

}      
#endif