initial.h

Go to the documentation of this file.

#ifndef INITIAL_H
#define INITIAL_H
 
#include <vector>
#include <cmath>
 
class Particle {
public:
    double x, y;  
 
    Particle(double x_init = 0.0, double y_init = 0.0);
 
    void updatePosition(double dx, double dy);
};
class SimulationBox {
private:
    double Lx;  
    double Ly;  
    double invLx;  
    double invLy;  
    double volume; 
 
public:
    SimulationBox(double Lx, double Ly);
 
    void applyPBC(Particle &p) const;
 
    double minimumImageDistance(const Particle &p1, const Particle &p2) const;
 
    double minimumImageDistanceSquared(const Particle &p1, const Particle &p2) const;
 
    double getLx() const;
 
    double getLy() const;
 
    double getV() const;
};
 
void initializeParticles(std::vector<Particle> &particles, const SimulationBox &box, int N, bool random = true, unsigned int seed = 0);
 
 
#endif // INITIAL_H