CircularPeriodicBoundary.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef CircularPeriodicBoundary_H
#define CircularPeriodicBoundary_H
 
#include "BaseBoundary.h"
 
class ParticleHandler;
 
class BaseParticle;
 
class CircularPeriodicBoundary : public BaseBoundary
{
public:
    CircularPeriodicBoundary();
    
    explicit CircularPeriodicBoundary(double innerRadius);
    
    ~CircularPeriodicBoundary() override;
    
    CircularPeriodicBoundary* copy() const override;
    
    void rotateParticle(BaseParticle* P, double angle);
    
    void createPeriodicParticle(BaseParticle* p, ParticleHandler& pH) override;
    
    void createPeriodicParticles(ParticleHandler& pH) override;
    
    bool checkBoundaryAfterParticleMoved(BaseParticle* P, ParticleHandler& pH);
    
    void checkBoundaryAfterParticlesMove(ParticleHandler& pH) override;
    
    void read(std::istream& is) override;
    
    void oldRead(std::istream& is);
    
    void write(std::ostream& os) const override;
    
    std::string getName() const override;
 
private:
    
    double innerRadius_;
    
    //A particle is between to Radii in plane (i) and has two (straight) walls (plane 0 defines centre)
    //If it is close to its inner Radius it should be copied once with a positive rotation of 2*pi/2^i
    //If it is close to one of its straight walls is should be rotative with +/- 2*pi/2^i
    
    //Particles can only apply forces to real particles
    
    //If a particle crosses a straight wall it should simply be shifted
    //If a particle crosses its inner Radius it should be copied
    //If a particle crosses its outer Radius it may need to be deleted
    
};
 
#endif